Numerical simulation of 3D-microstructures in solidification processes based on the CAFE method

It was analyzed that the finite element-cellular automaton （CAFE） method was used to simulate 3D-microstructures in solidification processes. Based on this method, the 3D-microstructure of 9SMn28 free-cutting steel was simulated in solidification processes and the simulation results are consistent with the experimental ones. In addition, the effects of Gaussian distribution parameters were also studied. The simulation results show that the higher the mean undercooling, the larger the columnar dendrite zones, and the larger the maximum nucleation density, the smaller the size of grains. The larger the standard deviation, the less the number of minimum grains is. However, the uniformity degree decreases first, and then increases gradually.

Analysis of cranial aneurysm's 3D-computational hemodynamics numerical simulation and findings in operation

Objective With the cranial aneurysm’s 3D-computational hemodynamics numerical simulation,we can get the cranial aneurysms and its ’s mother arteries ’blood flow rate,path and wall shearing stress and so on. Taking use of the contrast analysis of the hemodynamics result and the findings in operation,we try to find hemodynamic character and mechanism of rupture

Numerical simulation and analysis of solid-liquid two-phase threedimensional unsteady flow in centrifugal slurry pump

Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.

Numerical simulations of full-wave fi elds and analysis of channel wave characteristics in 3-D coal mine roadway models

Currently, numerical simulations of seismic channel waves for the advance detection of geological structures in coal mine roadways focus mainly on modeling two- dimensional wave fields and therefore cannot accurately simulate three-dimensional （3-D） full-wave fields or seismic records in a full-space observation system. In this study, we use the first-order velocity-stress staggered-grid finite difference algorithm to simulate 3-D full-wave fields with P-wave sources in front of coal mine roadways. We determine the three components of velocity Vx, Vy, and Vz for the same node in 3-D staggered-grid finite difference models by calculating the average value of Vy, and Vz of the nodes around the same node. We ascertain the wave patterns and their propagation characteristics in both symmetrical and asymmetric coal mine roadway models. Our simulation results indicate that the Rayleigh channel wave is stronger than the Love channel wave in front of the roadway face. The reflected Rayleigh waves from the roadway face are concentrated in the coal seam, release less energy to the roof and floor, and propagate for a longer distance. There are surface waves and refraction head waves around the roadway. In the seismic records, the Rayleigh wave energy is stronger than that of the Love channel wave along coal walls of the roadway, and the interference of the head waves and surface waves with the Rayleigh channel wave is weaker than with the Love channel wave. It is thus difficult to identify the Love channel wave in the seismic records. Increasing the depth of the receivers in the coal walls can effectively weaken the interference of surface waves with the Rayleigh channel wave, but cannot weaken the interference of surface waves with the Love channel wave. Our research results also suggest that the Love channel wave, which is often used to detect geological structures in coal mine stopes, is not suitable for detecting geological structures in front of coal mine roadways. Instead, the Rayleigh channel wave can be used for the advance detection of geological structures in coal mine roadways.

Numerical simulation for influence of excavation and blasting vibration on stability of mined-out area

Dynamic analysis steps and general flow of fast lagrangian analysis of continua in 3 dimensions(FLAC3D) were discussed. Numerical simulation for influence of excavation and blasting vibration on stability of mined-out area was carried out with FLAC3D. The whole analytical process was divided into two steps, including the static analysis and the dynamic analysis which were used to simulate the influence of excavation process and blasting vibration respectively. The results show that the shape of right upper boundary is extremely irregular after excavation, and stress concentration occurs at many places and higher tensile stress appears. The maximum tensile stress is higher than the tensile strength of rock mass, and surrounding rock of right roof will be damaged with tension fracture. The maximum displacement of surrounding rock is 4.75 mm after excavation. However, the maximum displacement increases to 5.47 mm after the blasting dynamic load is applied. And the covering area of plastic zones expands obviously, especially at the foot of right upper slope. The analytical results are in basic accordance with the observed results on the whole. Damage and disturbance on surrounding rock to some degree are caused by excavation, while blasting dynamic load increases the possibility of occurrence of dynamic instability and destruction further. So the effective supporting and vibration reducing measures should be taken during mining.

Study of 3-D Numerical Simulation for Gas Transfer in the Goaf of the Coal Mining

In order to simulate field distribution rules,mathematical models for 3-D air flows and gas transfer in the goaf of the coal mining are established,based on theories of permeability and dynamic dispersion through porous media. A gas dispersion equation in a 3-D field is calculated by use of numerical method on a weighted upstream multi-element balance. Based on data of an example with a U type ventilation mode,surface charts of air pressure distribution and gas concentration are drawn by Graphtool software. Finally,a comparison between actually measured results in the model test and the numerical simulation results is made to proves the numerical implementation feasible.

Numerical Simulation of the Protective Effect of Complex Boundaries Toward Shock Waves in a 3D Explosive Field

A numerical method is presented that simulates 3D explosive field problems. A code MMIC3D using this method can be used to simulate the propagation and reflected effects of all kinds of rigid boundaries to shock waves produced by an explosive source. These numerical results indicate that the code MMIC3D has the ability in computing cases such as 3D shock waves produced by air explosion, vortex region of the shock wave, the Mach wave, and reflected waves behind rigid boundaries.

Numerical Simulation of 2D and 3D Sloshing Waves in a Regularly and Randomly Excited Container

In this paper,various aspects of the 2D and 3D nonlinear liquid sloshing problems in vertically excited containers have been studied numerically along with the help of a modified-transformation.Based on this new numerical algorithm,a numerical study on a regularly and randomly excited container in vertical direction was conducted utilizing four different cases： The first case was performed utilizing a 2D container with regular excitations.The next case examined a regularly excited 3D container with two different initial conditions for the liquid free surface,and finally,3D container with random excitation in the vertical direction.A grid independence study was performed along with a series of validation tests.An iteration error estimation method was used to stop the iterative solver（used for solving the discretized governing equations in the computational domain） upon reaching steady state of results at each time step.In the present case,this method was found to produce quite accurate results and to be more time efficient as compared to other conventional stopping procedures for iterative solvers.The results were validated with benchmark results.The wave elevation time history,phase plane diagram and surface plots represent the wave nonlinearity during its motion.

Quasi-3D Numerical Simulation of Tidal Hydrodynamic Field

Based on the 2D horizontal plane numerical model, a quasi-3D numerical model is established for coastal regions of shallow water. The characteristics of this model are that the velocity profiles;can be obtained at the same time when the equations of the value of difference between the horizontal current velocity and its depth-averaged velocity in the vertical direction are solved and the results obtained are consistent with the results of the 2D, model. The circulating flow in the rectangular area induced by wind is simulated and applied to the tidal flow field of the radial sandbanks in the South Yellow Sea. The computational results from this quasi-3D model are in good agreement with analytical results and observed data. The solution of the finite difference equations has been found to be stable, and the model is simple, effective and practical.

An explicit method for numerical simulation of wave equations: 3D wave motion

In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order （M=1） recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.

Numerical Simulation of Injection Molding Cooling Process Based on 3D Surface Model

The design of the cooling system of injection molds directly affects both productivity and the quality of the final part. Using the cooling process CAE system to instruct the mold design, the efficiency and quality of design can be improved greatly. At the same time, it is helpful to confirm the cooling system structure and optimize the process conditions. In this paper, the 3D surface model of mold cavity is used to replace the middle-plane model in the simulation by Boundary Element Method, which break the bottleneck of the application of the injection molding simulation softwares base on the middle-plane model. With the improvements of this paper, a practical and commercial simulation software of injection molding cooling process named as HsCAE3D6.0 is developed.

Numerical simulation of intermediate phase growth in Ti/Al alternate foils

To investigate the diffusion reaction between Ti/Al solid diffusion couple, Ti/Al alternate foils formed by hot pressing were annealed at 525, 550, 575 and 600 °C for time ranging from 1 to 40 h. The experimental results show that TiAl3 was the only observed phase at Ti/Al interface. The interface thermodynamics favored the preferential formation of TiAl3 in Ti/Al couple. The growth of TiAl3 layer occurred mainly towards Al foil side and exhibited a parabolic law. Using the interdiffusion coefficients calculated based on the contribution of grain boundary diffusion, the growth of TiAl3 was simulated numerically with the finite difference method, and the simulated results were in good agreement with the experimental ones.

NUMERICAL SIMULATION OF 3-D TURBULENT FLOWS OVER DREDGED TRENCHES

A 3- D free surface flow in open channels based on the Reynolds equations with the k-ε turbulence closure model is presented in this paper. Insted of the 'rigid lid' approximation, the solution of the free surface equation is implemented in the velocity-pressure iterative procedure on the basis of the conventional SIMPLE method. This model was used to compute the flow in rectangular channels with trenches dredged across the bottom. The velocity, eddy viscosity coefficient, turbulent shear stress, turbulent kinetic energy and elevation of the free surface can be obtained. The computed results are in good agreement with previous experimental data.

Mechanical Properties of Nano-CaCO_(3) and Basalt Fiber Reinforced Concrete:Experiments and Numerical Simulations

In this study,the compressive,split tensile,and flexural strengths of concrete with nano-CaCO_(3) only were compared with those of concrete with nano-CaCO_(3) and basalt fibers through field experiments,and the underlying mechanisms were analyzed by the Scanning Electron Microscope (SEM) techniques.On the mesoscale,a damage model of concrete was established based on the continuum progressive damage theory,which was used to investigate the influence of different lengths and contents of fibers on the mechanical properties of concrete.Then,the experimental and numerical simulation results were compared and analyzed to verify the feasibility of model.The results show that nano-CaCO_(3) can enhance the compressive strength of the concrete,with an optimal content of 2.0%.Adding basalt fibers into the nano-CaCO_(3) reinforced concrete may further enhance the compressive,split tensile,and flexural strengths of the concrete;however,the higher content of basalt fiber can not lead to higher performance of concrete.The optimal length and content of fiber are 6 mm and 0.20%,respectively.The SEM result shows that the aggregation of basalt fibers is detrimental to the mechanical properties of concrete.The numerical simulation results are in good agreement with the experimental results.

3-D Finite-Element Numerical Simulation of Centrifugal Induction Electrosalg Casting Solidification Process

A 3-D finite-element numerical simulation model of temperature field for CIESC casting solidification process was developed with the aid of ANSYS software and a series of corresponding experiments were made. The results showed that the good agreement was obtained between the numerical simulation and the experiments. Based on the numerical simulation results, the characteristics of temperature distribution in the castings during CIESC solidification process were analyzed and summarized. According to the G/R-1/2 method and numerical simulation results, there is no any shrinkage defect in the CIESC casting and structure or casting is fine and compact.

Numerical Simulation Study on the Detection of Weak Structural Plane of Rock Slope by Using 3D Electrical Resistivity Tomography

The weak structure plane is an important factor affecting the stability of rock slope, and detecting the spatial structure of the weak structural plane is beneficial to analyze the stability of the slope and estimate the quantity of the landslide. Based on 3D electrical resistivity tomography, a model of rock slope with weak structure plane is established, and the characteristics of three-dimensional resistivity imaging of weak structure plane under different ground water conditions are simulated. The results show that the weak structural plane has a better reflection in 3D electrical resistivity tomography;the distribution of weak structural plane of 3D resistivity imaging can be roughly determined under different ground water conditions;the three-dimensional electrical resistivity tomography is feasible in the detection of weak structural plane of rock slope.

Temperature Distributions in Single Cell of Polymer Electrolyte Fuel Cell Simulated by an 1D Multi-plate Heat-Transfer Model and a 3D Numerical Simulation Model

The purpose of this study is to verify an 1D multi-plate heat-transfer model estimating the temperature distribution on the interface between polymer electrolyte membrane and catalyst layer at cathode in single cell of polymer electrolyte fuel cell, which is named as reaction surface in this study, with a 3D numerical simulation model solving many governing equations on the coupling phenomena in the cell. The results from both models/simulations agreed well. The effects of initial operation temperature, flow rate, and relative humidity of supply gas on temperature distribution on the reaction surface were also investigated. It was found in both 1D and 3D simulations that, the temperature rise （i.e., Treact-Tini） of the reaction surface from initial operation temperature at 70℃ was higher than that at 80℃ irrespective of flow rate of supply gas. The effect of relative humidity of supply gas on Treact- Tini near the inlet of the cell was small. Compared to the previous studies conducted under the similar operation conditions, the Treact - Tini calculated by 1D multi-plate heat-transfer model in this study as well as numerical simulation using 3D model was reasonable.

Penetration resistance of Al_(2)O_(3) ceramic-ultra-high molecular weight polyethylene(UHMWPE)composite armor:Experimental and numerical investigations

To enhance the protective performance of ceramic composite armor,ballistic penetration experiments were conducted on Al_(2)O_(3) ceramic-ultra-high molecular weight polyethylene(UHMWPE)composite armor with different thickness configurations.The damage and failure modes of hard projectiles and ceramic-fiber composite targets were analyzed.The recovered projectiles and ceramic fragments were sieved and weighed at multiple stages,revealing a positive correlation between the degree of fragmentation of the projectiles and ceramics and the overall ballistic resistance of the composite targets.Numerical simulations were performed using the LS-DYNA finite element software,and the simulation results showed high consistency with the experimental results,confirming the validity of the material parameters.The results indicate that the projectile heads primarily exhibited crushing and abrasive fragmentation.Larger projectile fragments mainly resulted from tensile and shear stress-induced failure.The failure modes of the composite targets included the formation of ceramic cones and radial cracks under high-velocity impacts.The UHMWPE laminated plates exhibited interlayer separation caused by tensile waves,permanent plastic deformation of the rear surface bulging,and perforation failure primarily due to shear forces.Through extended numerical simulations,while maintaining the same areal density and configuration of9 mm Al_(2)O_(3) ceramic+12 mm UHMWPE laminated composite armor,the thickness configurations of the Al_(2)O_(3) ceramic and UHMWPE laminated backplates were varied,and various thicknesses of UHMWPE laminates were simulated as the cover layer for the ceramic panels.The simulation results indicated that the composite armor configuration of 10 mm Al_(2)O_(3) ceramic+8 mm UHMWPE composite armor increased energy absorption by13.48%.When altering the cover layer thickness,a 4 mm UHMWPE+9 mm Al_(2)O_(3)+8 mm UHMWPE composite armor demonstrated a 27.11%improvement in energy absorption,showing a relatively significant enhancement.

Dynamic response and numerical simulation of Al-Sc and Al-Ti alloys under high-speed impact

Al-Sc and Al-Ti semi-infinite targets were impacted by high-speed projectiles at velocities of 0.8, 1.0, 1.2 and 1.5 km/s, respectively. It is found that the Al-Sc targets demonstrate more excellent ability to resist high-speed impact. It is concluded that different microstructures of Al-Sc and Al-Ti alloys, including different grain sizes and secondary particles precipitated in the matrix, result in their greatly different capabilities of resisting impact. Furthermore, the effect of the size range ofnanoscale A13Sc precipitate in A1-Sc alloy on the resistance of high-speed impact was investigated. In addition, computer simulations and validation of these simulations were developed which fairly accurately represented residual crater shapes/geometries. Validated computer simulations allowed representative extrapolations of impact craters well beyond the laboratory where melt and solidification occurred at the crater wall, especially for hypervelocity impact （〉5 km/s）.