Three-Dimensional Direct Numerical Simulation of Flow past A Near-Wall Circular Cylinder:Combined Effects of Gap Ratio and Boundary Layer Thickness on Flow Profiles and Pressure Distribution

Three-dimensional direct numerical simulations of the wake flow downstream of a near-wall circular cylinder at different gap ratios and boundary layer thicknesses are carried out by using the iterative immersed boundary method.The non-dimensional gap between the cylinder and the wall,G/D=0.2,0.6 and 1.0,the non-dimensional boundary layer thickness,δ/D=0.0,0.7 and 1.6,the Reynolds number,Re=350,and the aspect ratio of the cylinder,L/D=25are adopted.High-resolution visualizations of the complex vortex structures at differentδ/D and G/D are presented.The transition of the streamwise vortex mode,the combined effects ofδ/D and G/D on the flow statistics,the pressure and shear stress distribution and the hydrodynamic forces are analyzed.Results show that with decreasing G/D and increasingδ/D,the gap flow and its vortex-shedding are significantly weakened,together with an elongated wake and an enlarged low-velocity area near the wall,leading to the wake mode transition from the two-sided to one-sided vortex-shedding.Different relative positions of the cylinder regarding the boundary layer alter the flow features of the shear layers.With an increase inδ/D,the front stagnation point shifts to the upper surface,and the distance between the flow divergence point and the maximum pressure position increases.The mean drag coefficient and r.m.s.values of drag and lift coefficients decrease with a decrease in G/D and an increase inδ/D,while the mean lift coefficient increases with decreasing G/D but decreases with increasingδ/D.

Effects of Fe solid solute on grain boundaries of bi-crystal Cu: A molecular dynamics simulation

Grain boundaries(GBs)play a crucial role on the structural stability and mechanical properties of Cu and its alloys.In this work,molecular dynamics(MD)simulations are employed to study the effects of Fe solutes on the formation energy,excess volume,dislocations and melting behaviors of GBs in CuFe alloys.It is illustrated that Fe solute affects the structural stability of Cu GBs substantially,the formation energy of GBs is reduced,but the thickness and melting point of GBs are increased,that is,the structural stability of Cu GBs is significantly improved owing to the Fe solutes.A strong scaling law exists between the formation energy,excess volume,thickness and melting point of GBs.Therefore,Fe solid solute plays an important role in the characteristics of GBs in bi-crystal Cu.

Numerical simulation of three-dimensional fracturing fracture propagation in radial wells

A fracture propagation model of radial well fracturing is established based on the finite element-meshless method.The model considers the coupling effect of fracturing fluid flow and rock matrix deformation.The fracture geometries of radial well fracturing are simulated,the induction effect of radial well on the fracture is quantitatively characterized,and the influences of azimuth,horizontal principle stress difference,and reservoir matrix permeability on the fracture geometries are revealed.The radial wells can induce the fractures to extend parallel to their axes when two radial wells in the same layer are fractured.When the radial wells are symmetrically distributed along the direction of the minimum horizontal principle stress with the azimuth greater than 15,the extrusion effect reduces the fracture length of radial wells.When the radial wells are symmetrically distributed along the direction of the maximum horizontal principal stress,the extrusion increases the fracture length of the radial wells.The fracture geometries are controlled by the rectification of radial borehole,the extrusion between radial wells in the same layer,and the deflection of the maximum horizontal principal stress.When the radial wells are distributed along the minimum horizontal principal stress symmetrically,the fracture length induced by the radial well decreases with the increase of azimuth;in contrast,when the radial wells are distributed along the maximum horizontal principal stress symmetrically,the fracture length induced by the radial well first decreases and then increases with the increase of azimuth.The fracture length induced by the radial well decreases with the increase of horizontal principal stress difference.The increase of rock matrix permeability and pore pressure of the matrix around radial wells makes the inducing effect of the radial well on fractures increase.

Three-Dimensional Simulation of Hydrodynamic Mechanism of Fluidized Bed Methanation

Organic solid waste(OSW)contains many renewable materials.The pyrolysis and gasification of OSW can realize resource utilization,and its products can be used for methanation reaction to produce synthetic natural gas in the specific reactor.In order to understand the dynamic characteristics of the reactor,a three-dimensional numerical model has been established by the method of Computational Fluid Dynamics(CFD).Along the height of the reactor,the particle distribution in the bed becomes thinner and the mean solid volume fraction decreases from 4.18%to 0.37%.Meanwhile,the pressure fluctuation range decreased from 398.76 Pa at the entrance to a much lower value of 74.47 Pa at the exit.In this simulation,three parameters of gas inlet velocity,operating temperature and solid particle diameter are changed to explore their influences on gas-solid multiphase flow.The results show that gas velocity has a great influence on particle distribution.When the gas inlet velocity decreases from 6.51 to 1.98 m/s,the minimum height that particles can reach decreases from 169 to 100 mm.Additionally,as the operating temperature increases,the particle holdup inside the reactor changes from 0.843%to 0.700%.This indicates that the particle residence time reduces,which is not conducive to the follow-up reaction.Moreover,with the increase of particle size,the fluctuation range of the pressure at the bottom of the reactor increases,and its standard deviation increases from 55.34 to 1266.37 Pa.

Temperature effect on nanotwinned Ni under nanoindentation using molecular dynamic simulation

Temperature effect on atomic deformation of nanotwinned Ni (nt-Ni) under localized nanoindentation is investigated in comparison with nanocrystalline Ni (nc-Ni) through molecular simulation.The nt-Ni exhibits enhanced critical load and hardness compared to nc-Ni,where perfect,stair-rod and Shockley dislocations are activated at (111),(111) and (111) slip planes in nt-Ni compared to only SSockley dislocation nucleation at (111) and (111) slip planes of nc-Ni.The nt-Ni exhibits a less significant indentation size effect in comparison with nc-Ni due to the dislocation slips hindrance of the twin boundary.The atomic deformation associated with the indentation size effect is investigated during dislocation transmission.Different from the decreasing partial slips parallel to the indenter surface in nc-Ni with increasing temperature,the temperaturedependent atomic deformation of nt-Ni is closely related to the twin boundary:from the partial slips parallel to the twin boundary (~10 K),to increased confined layer slips and decreased twin migration(300 K–600 K),to decreased confined layer slips and increased dislocation interaction of dislocation pinning and dissociation (900 K–1200 K).Dislocation density and atomic structure types through quantitative analysis are implemented to further reveal the above-mentioned dislocation motion and atomic structure alteration.Our study is helpful for understanding the temperature-dependent plasticity of twin boundary in nanotwinned materials.

Dynamic simulation insights into friction weakening effect on rapid long-runout landslides:A case study of the Yigong landslide in the Tibetan Plateau,China

This study proposed a novel friction law dependent on velocity,displacement and normal stress for kinematic analysis of runout process of rapid landslides.The well-known Yigong landslide occurring in the Tibetan Plateau of China was employed as the case,and the derived dynamic friction formula was included into the numerical simulation based on Particle Flow Code.Results showed that the friction decreased quickly from 0.64(the peak)to 0.1(the stead value)during the 5s-period after the sliding initiation,which explained the behavior of rapid movement of the landslide.The monitored balls set at different sections of the mass showed similar variation characteritics regarding the velocity,namely evident increase at the initial phase of the movement,followed by a fluctuation phase and then a stopping one.The peak velocity was more than 100 m/s and most particles had low velocities at 300s after the landslide initiation.The spreading distance of the landslide was calculated at the two-dimension(profile)and three-dimension scale,respectively.Compared with the simulation result without considering friction weakening effect,our results indicated a max distance of about 10 km from the initial unstable position,which fit better with the actual situation.

Parallelization strategies for resolved simulations of fluid-structure-particle interactions

Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.

Detailed analysis and simulation verification for reconstruction of plasma optical boundary with spectrometric technique on HL-2M tokamak

The plasma optical boundary reconstruction technique based on Hommen's theory is promising for future tokamaks with high parameters. In this work, we conduct detailed analysis and simulation verification to estimate the ‘logic loophole' of this technique. The finite-width effect and unpredictable errors reduce the technique's reliability, which leads to this loophole. Based on imaging theory, the photos of a virtual camera are simulated by integrating the assumed luminous intensity of plasma. Based on Hommen's theory, the plasma optical boundary is reconstructed from the photos. Comparing the reconstructed boundary with the one assumed, the logic loophole and its two effects are quantitatively estimated. The finite-width effect is related to the equivalent thickness of the luminous layer, which is generally about 2-4 cm but sometimes larger. The level of unpredictable errors is around 0.65 cm. The technique based on Hommen's theory is generally reliable, but finite-width effect and unpredictable errors have to be taken into consideration in some scenarios. The parameters of HL-2M are applied in this work.

Absorbing Boundary Conditions for Simulating Water Waves Near Solid Bodies

The objective of this paper is to present a new method for designing absorbing or non-reflective boundary conditions (ABC) or (NRBC), illustrated by the case study of the modelling of a solid body in water, specifically the capillary gravity waves generated by its motion at the surface. The study analyses the flow of an inviscid, barotropic, and compressible fluid around the stationary solid body. The dynamic behaviour of the fluid is analysed using a two-dimensional coupled Neumann-Kelvin model extended with capillarity and inertia terms. For computational purposes, it is necessary to truncate the unbounded spatial domain with artificial boundaries and then introduce appropriate absorbing boundary conditions. The propagation of short wavelength waves in a convective fluid medium with significant differences in properties between the interior and the surface of the fluid presents a number of difficulties in the design of these conditions. The results are illustrated numerically and commented upon.

Three-dimensional acoustic propagation model for shallow waters based on an indirect boundary element method

This work has a two-fold purpose.On the one hand,the theoretical formulation of a three-dimensional(3D)acoustic propagation model for shallow waters with a constant sound speed is presented,based on the boundary element method(BEM),which uses a half-space Green function instead of the more conventional free-space Green function.On the other hand,a numerical implementation is illustrated to explore the formulation in simple idealized cases,controlled by a few parameters,which provides necessary tests for the accuracy and performance of the model.The half-space Green's function,which has been previously used in scattering and diffraction,adds terms to the usual expressions of the integral operators without altering their continuity properties.Verifications against the wavenumber integration solution of the Pekeris waveguide suggest that the model allows an adequate prediction for the acoustic field.Likewise,numerical experiments in relation to the necessary mesh size for the description of the water-marine sediment interface lead to the conclusion that a transmission loss prediction with acceptable accuracy can be obtained with the use of a limited mesh around the desired evaluation region.

Numerical Simulation of Nonlinear Three-Dimensional Waves in Water of Arbitrary Varying Topography

The numerical simulation is based on the authors' high-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth. Corresponding finite-difference equations and general conditions for open and fixed natural boundaries with an arbitrary reflection coefficient and phase shift are also given in this paper. The systematical tests of numerical simulation show that the theoretical models, the finite-difference algorithms and the boundary conditions can give good calculation results for the wave propagating in shallow and deep water with an arbitrary slope varying from gentle to steep.

Three-Dimensional and Cross-sectional Characteristics of Normal Grain Growth Based on Monte Carlo Simulation

An appropriate Monte Carlo method was developed to simulate the three-dimensional normal grain growth more completely. Comparative investigation on the three-dimensional and the cross-sectional characteristics of normal grain growth was done. It was found that the time exponent of grain growth determined from cross-section exhibits the same rule of increasing slowly with time and approaching the theoretical value n = 0.5 of steadygrain growth as the three-dimensional (3-D) system. From change of the number of grains per unit area with timemeasured in cross-section, the state of 3-D normal grain growth may be predicted. The gtain size distribution incross-section is different from that in 3-D system and can not express the evolution characteristic of the 3-D distribution. Furthermore, there exists statistical connection between the topological parameters in cross-section and thosein three-dimensions.

Three-dimensional visualization and virtual reality simulation role in hepatic surgery:Further research warranted

Artificial intelligence(AI)is the study of algorithms that enable machines to analyze and execute cognitive activities including problem solving,object and word recognition,reduce the inevitable errors to improve the diagnostic accuracy,and decision-making.Hepatobiliary procedures are technically complex and the use of AI in perioperative management can improve patient outcomes as discussed below.Three-dimensional(3D)reconstruction of images obtained via ultrasound,computed tomography scan or magnetic resonance imaging,can help surgeons better visualize the surgical sites with added depth perception.Preoperative 3D planning is associated with lesser operative time and intraoperative complications.Also,a more accurate assessment is noted,which leads to fewer operative complications.Images can be converted into physical models with 3D printing technology,which can be of educational value to students and trainees.3D images can be combined to provide 3D visualization,which is used for preoperative navigation,allowing for more precise localization of tumors and vessels.Nevertheless,AI enables surgeons to provide better,personalized care for each patient.

Three-Dimensional Finite Element Numerical Simulation and Physical Experiment for Magnetism-Stress Detecting in Oil Casing

The casing damage has been a big problem in oilfield production. The current detection methods mostly are used after casing damage, which is not very effective. With the rapid development of China's offshore oil industry, the number of offshore oil wells is becoming larger and larger. Because the cost of offshore oil well is very high, the casing damage will cause huge economic losses. What's more, it can also bring serious pollution to marine environment. So the effective methods of detecting casing damage are required badly. The accumulation of stress is the main reason for the casing damage. Magnetic anisotropy technique based on counter magnetostriction effect can detect the stress of casing in real time and help us to find out the hidden dangers in time. It is essential for us to prevent the casing damage from occurring. However, such technique is still in the development stage. Previous studies mostly got the relationship between stress and magnetic signals by physical experiment, and the study of physical mechanism in relative magnetic permeability connecting the stress and magnetic signals is rarely reported. The present paper uses the ANSYS to do the three-dimensional finite element numerical simulation to study how the relative magnetic permeability works for the oil casing model. We find that the quantitative relationship between the stress' s variation and magnetic induction intensity's variation is: Δδ =K* ΔB, K = 8.04×109, which is proved correct by physical experiment.

Three-dimensional discrete element numerical simulation of Paleogene salt structures in the western Kuqa foreland thrust belt

Taking the Paleogene salt strata in the west of Kuqa foreland thrust belt as study object, the deformation features of salt structure in the compression direction and perpendicular to the compression direction were examined to find out the control factors and formation mechanisms of the salt structures. By using the three-dimensional discrete element numerical simulation method, the formation mechanisms of typical salt structures of western Kuqa foreland thrust belt in Keshen and Dabei work areas were comprehensively analyzed. The simulation results show that the salt deformation in Keshen and Dabei work areas is of forward spread type, with deformation concentrated in the piedmont zone;the salt deformation is affected by the early uplift near the compression end, pre-existing basement faults, synsedimentary process and the initial salt depocenter;in the direction perpendicular to the compression direction, salt rocks near the compression end have strong lateral mobility with the velocity component moving towards the middle part, and the closer to the middle, the larger the velocity will be, so that salt rocks will aggregate towards the middle and deform intensely, forming complex folds and separation of salt structures from salt source, and local outcrop with thrust faults. Compared with 2 D simulation, 3 D simulation can analyze salt structures in the principal stress direction and direction perpendicular to the principal stress, give us a full view of the formation mechanisms of salt structures, and guide the exploration of oil and gas reservoirs related to salt structures.

THE SYSTEM SIMULATION OF THREE-DIMENSIONAL RADAR

To provide a test platform for Electronic Warfare (EW) system, it is needed to simulate the radar received Intermediate Frequency (IF) signals and radar system functions.This letter gives a description of a radar system simulation software developed for frequencyphase scanning three-dimensional (3-D) radar. Experimental results prove that the software could be used for system evaluation and for training purposes as an attractive alternative to real EW system.

Three-Dimensional Simulations of RESET Operation in Phase-Change Random Access Memory with Blade-Type Like Phase Change Layer by Finite Element Modeling

An optimized device structure for reducing the RESET current of phase-change random access memory （PCRAM） with blade-type like （BTL） phase change layer is proposed. The electrical thermal analysis of the BTL cell and the blade heater contactor structure by three-dimensional finite element modeling are compared with each other during RESET operation. The simulation results show that the programming region of the phase change layer in the BTL cell is much smaller, and thermal electrical distributions of the BTL cell are more concentrated on the TiN/GST interface. The results indicate that the BTL cell has the superiorities of increasing the heating efficiency, decreasing the power consumption and reducing the RESET current from 0.67mA to 0.32mA. Therefore, the BTL cell will be appropriate for high performance PCRAM device with lower power consumption and lower RESET current.

Three-Dimensional Rigid-Plastic FEM Simulation of Bulk Forming Processes with New Contact and Remeshing Techniques

Some techniques such as die surface description, contact judgement algorithm and remeshing are proposed to improve the robustness of the numerical solution. Based on these techniques, a three-dimensional rigid-plastic FEM code has been developed. Isothermal forging process of a cylindrical housing has been simulated. The simulation results show that the given techniques and the FEM code are reasonable and feasible for three-dimensional bulk forming processes.

Numerical simulations and comparative analysis of two- and three-dimensional circulating fluidized bed reactors for CO2 capture

Carbon dioxide(CO2),the main gas emitted from fossil burning,is the primary contributor to global warming.Circulating fluidized bed reactor(CFBR)is proved as an energy-efficient method for post-combustion CO2 capture.The numerical simulation by computational fluid dynamics(CFD)is believed as a promising tool to study CO2 adsorption process in CFBR.Although three-dimensional(3D)simulations were proved to have better predicting performance with the experimental results,two-dimensional(2D)simulations have been widely reported for qualitative and quantitative studies on gas-solid behavior in CFBR for its higher computational efficiency recently.However,the discrepancies between 2D and 3D simulations have rarely been evaluated by detailed study.Considering that the differences between the 2D and 3D simulations will vary substantially with the changes of independent operating conditions,it is beneficial to lower computational costs to clarify the effects of dimensionality on the numerical CO2 adsorption runs under various operating conditions.In this work,the comparative analysis for CO2 adsorption in 2D and 3D simulations was conducted to enlighten the effects of dimensionality on the hydrodynamics and reaction behaviors,in which the separation rate,species distribution and hydrodynamic characteristics were comparatively studied for both model frames.With both accuracy and computational costs considered,the viable suggestions were provided in selecting appropriate model frame for the studies on optimization of operating conditions,which directly affect the capture and energy efficiencies of cyclic CO2 capture process in CFBR.

Simulation-Based Construction of Three-Dimensional Process Model for Punching Cartridge Cases

A method of constructing three-dimensional process model for the punching cartridge cases is presented based on DEFORM simulation analysis. Using DEFORM software,the finite element simulation models for the punching and forming process of cartridge cases are established,and the corresponding simulation result model of each intermediate procedure is obtained by continuously performing the forming process simulation. The simulation model cannot annotate size and process information due to poor interface between DEFORM software and CAD software. Thus,a 3D annotation module is developed with secondary development technology of UG NX software. Consequently,the final process model with dimension and process information is obtained. Then,with the current 3D process management system,the 3D punching and forming process design of cartridge cases can be completed further. An example is also provided to illustrate that the relative error between the simulation process model and the physical model is less than 2%,which proves the validity and reliability of the proposed method in this study.