Three-dimensional transient numerical simulation for intake process in the engine intake port-valve-cylinder system

This paper presents a KIVA-3 code based numerical model for three-dimensional transient intake flow in the intake port-valve-cylinder system of internal combustion engine using body-fitted technique, which can be used in numerical study on internal combustion engine with vertical and inclined valves, and has higher calculation precision. A numerical simulation (on the intake process of a two-valve engine with a semi-sphere combustion chamber and a radial intake port) is provided for analysis of the velocity field and pressure field of different plane at different crank angles. The results revealed the formation of the tumble motion, the evolution of flow field parameters and the variation of tumble ratios as important information for the design of engine in-take system.

Three-dimensional transient numerical simulation for gas exchange process in a four-stroke motorcycle engine

Three-dimensional transient numerical simulation of gas exchange process in a four-stroke motorcycle engine with a semi-spherical combustion chamber with two tilt valves was studied. Combination of the grid re-meshing method and the snapper technique made the valves move smoothly. The flow structure and pattern in a complete engine cycle were described in detail. Tumble ratios around the x-axis and y-axis were analyzed. Comparison of computed pressure with experimental pressure under motored condition revealed that the simulation had high calculation precision; CFD simulation can be regarded as an im-portant tool for resolving the complex aerodynamic behavior in motorcycle engines.

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.

Numerical simulation of three-dimensional tidal current in the Bobal Sea

The model equations with tbree-dimensional, time-dependent, nonlinear Navier-Stokes equations are transformed by sigma-transformation.On the basis of the process splitting technique, the fluid flow problems are divided into two parts:the vertically-intopated equations (external mode) and the vertical structure equations(internal mace). The first set of equations being the propagation of the tidal weves and the ADI numerical scheme has ben chosen to solve them. Conerning the vertical structure equations, they are solved by means of leapfrog stepping procedure.The main features of the tide and associated tidal current in the Bohai Sea are examined with this 3-D model.To have a good reproduction of vertical structure, the column is divided into 10 layers and the M2 tidal current is computed in detail. The simulation reveal the spetial structure and some important characteristics of the tidal current of the Bohai Sea. The application of the 3-D madel to forecasting of the tidal current in the Bobal Sea has been Performed as an illustration.

Numerical Simulation of the Three-Dimensional Wave-Induced Currents on Unstructured Grid

By coupling the three-dimensional hydrodynamic model with the wave model, numerical simulations of the three- dimensional wave-induced current are carried out in this study. The wave model is based on the numerical solution of the modified wave action equation and eikonal equation, which can describe the wave refraction and diffraction. The hydrodynamic model is driven by the wave-induced radiation stresses and affected by the wave turbulence. The numerical implementation of the module has used the finite-volume schemes on unstructured grid, which provides great flexibility for modeling the waves and currents in the complex actual nearshore, and ensures the conservation of energy propagation. The applicability of the proposed model is evaluated in calculating the cases of wave set-up, longshore currents, undertow on a sloping beach, rip currents and meandering longshore currents on a tri-cuspate beach. The results indicate that it is necessary to introduce the depth-dependent radiation stresses into the numerical simulation of wave-induced currents, and comparisons show that the present model makes better prediction on the wave procedure as well as both horizontal and vertical structures in the wave-induced current field.

A three-dimensional Eulerian method for the numerical simulation of high-velocity impact problems

In the present paper, a three-dimensional （3D） Eulerian technique for the 3D numerical simulation of high-velocity impact problems is proposed. In the Eulerian framework, a complete 3D conservation element and solution element scheme for conservative hyperbolic governing equations with source terms is given. A modified ghost fluid method is proposed for the treatment of the boundary conditions. Numerical simulations of the Taylor bar problem and the ricochet phenomenon of a sphere impacting a plate target at an angle of 60~ are carried out. The numerical results are in good agreement with the corresponding experimental observations. It is proved that our computational technique is feasible for analyzing 3D high-velocity impact problems.

Numerical Simulation of Transient Thermal Stress Field for Laser Metal Deposition Shaping

To decrease thermal stress during laser metal deposition shaping(LMDS)process,it is of great importance to learn the transient thermal stress distribution regularities.Based on the“element life and death”technique of finite element analy- sis(FEA),a three-dimensional multi-track and multi-layer numerical simulation model for LMDS is developed with ANSYS parametric design language(APDL)for the first time,in which long-edge parallel reciprocating scanning paths is introduced. Through the model,detailed simulations of thermal stress during whole metal cladding process are conducted,the generation and distribution regularities of thermal stress are also discussed in detail.Using the same process parameters,the simulation results show good agreement with the features of samples which fabricated by LMDS.

Numerical simulation of transient flow performance during different periods in centrifugal pumpNumerical simulation of transient flow performance during different periods in centrifugal pump

The instantaneous variations of the hydraulic characteristics take place in centrifugal pumps during their start-up,shutdown and other variable speed operations.In this paper,the variable speed method was proposed to simulate the transient internal flow field and the external performance of the pump during starting and stopping periods.The terms of accelerations due to variable speeds in the flow governing equations were analyzed in a multiple reference of frame(MRF).A transient CFD simulation was performed for a typical centrifugal pump by using ANSYS-CFX with the standard k-εturbulence model.The entire simulation process was composed of four stages:start-up,normal run,shutdown and post-shutdown.The function of rotating speed with regard to time was set by CEL language directly into the impeller domain in the pre-processor of the software to conduct variable speed simulation.The variations of the flow field in the centrifugal pump were obtained from the transient simulation.The changing laws of flow rate,head and other performance parameters over time were also analyzed and summarized.

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.

Numerical Simulation of Multi-track and Multi-layer Temperature Field on Laser Direct Metal Shaping

To improve the mechanical properties of the parts fabricated by Laser Direct Metal Shaping(LDMS),it is of great significance to understand the distribution regularities of transient temperature field during LDMS process.Based on the“el- ement birth and death”technique of finite element method,a three-dimensional multi-track and multi-layer model for the transient temperature field analysis of LDMS is developed by ANSYS Parametric Design Language(APDL)for the first time.In the fab- ricated modal,X-direction parallel reciprocating scanning paths is introduced.Using the same process parameters,the simulation results show good agreement with the microstructure features of samples which fabricated by LDMS.

Numerical simulation of soft palate movement and airflow in human upper airway by fluid-structure interaction method

In this paper, the authors present airflow field characteristics of human upper airway and soft palate movement attitude during breathing. On the basis of the data taken from the spiral computerized tomography images of a healthy person and a patient with Obstructive Sleep Apnea-Hypopnea Syndrome （OSAHS）, three-dimensional models of upper airway cavity and soft palate are reconstructed by the method of surface rendering. Numerical simulation is performed for airflow in the upper airway and displacement of soft palate by fluid-structure interaction analysis. The reconstructed threedimensional models precisely preserve the original configuration of upper airways and soft palate. The results of the pressure and velocity distributions in the airflow field are quantitatively determined, and the displacement of soft palate is presented. Pressure gradients of airway are lower for the healthy person and the airflow distribution is quite uniform in the case of free breathing. However, the OSAHS patient remarkably escalates both the pressure and velocity in the upper airway, and causes higher displacement of the soft palate. The present study is useful in revealing pathogenesis and quantitative mutual relationship between configuration and function of the upper airway as well as in diagnosingdiseases related to anatomical structure and function of the upper airway.

Pre-reinforcement grout in fractured rock masses and numerical simulation for optimizing shrinkage stoping configuration

Proper room and pillar sizes are both critical factors for safe mining and high ore recovery rate in shrinkage stoping mining of underground metal mines. The rock masses of Tangdan copper mine of China are fractured, which needs much reinforcement and support prior to mining. Cement-sodium silicate grout technology was selected, then its related parameters such as grout pressure, diffusion radius and time were calculated and proposed. In order to test the effect of the pressured grout in the fractured No.4 ore block, field experiments were conducted. To optimize stoping configuration, three-dimensional numerical simulation with ANSYS and FLAC 3 D softwares was proposed. The results show that the drilling porosity and mechanical properties of the rock masses are increased obviously. After grout, ore recovery rate is increased by 10.2 % employing the newly designed stoping configuration compared with the previous. Last, analyzed from the surface movements, roof subsidence and the maximum principal stress of the pillars, the mining safety is probable of being ensured.

A numerical simulation of surface wave excitation in a rectangular planar-type plasma source

The principle of surface wave plasma discharge in a rectangular cavity is introduced simply based on surface plasmon polariton theory. The distribution of surface-wave electric field at the interface of the plasma-dielectric slab is investigated by using the three-dimensional finite-difference time-domain method （3D-FDTD） with different slotantenna structures. And the experimental image of discharge with a novel slot antenna array and the simulation of the electric field with this slot antenna array are both displayed. Combined with the distribution of surface wave excitation and experimental results, the numerical simulation performed by using 3D-FDTD is shown to be a useful tool in the computer-aided antenna design for large area planar-type surface-wave plasma sources.

Three-dimensional tidal current numerical model of the Oujiang Estuary

The characteristics of three-dimensional （3-D） tidal current in the Oujiang Estuary are investigated according to in situ observations. The Oujiang Estuary has features of irregular coastline, complex topography, many islands, moveable boundary, and submerged dyke, therefore, σ 3-D numerical model oil an unstructured triangular grid has been degeloped. The σ coordinate transforination, the moveable boundary and submerged dyke treatment techniques were employed in the model so it is suitable for the tidal simulations in the Oujing Estuary with submerged dyke and moveable boundary problems. The model is evaluated with the in situ data, and the results show that the calculated water elevations at 19 stations and currents at 19 profiler stations are in good agreement with measured data both in magnitude and phase. This numerical model is applied to the 3-D tidal circulation simulations of experiments in stopping flow transport through the South Branch of the Oujiang Estuary, and the feasibility to cutoff the flow in the South Branch of the Oujiang Estuary is demonstrated by numerical simulation experiments. The developed numerical model simulated the 3-D tidal current circulations in complicated coastal and estuarine waters very well.

3-D NUMERICAL SIMULATION OF FRACTURE PROCESSES IN HETEROGENEOUS BRITTLE MATERIALS

By using the lattice model combined with finite element methods andstatistical techniques, a numerical approach is developed to establish mechanical models ofthree-dimensional heterogeneous brittle materials. A special numerical code is introduced, in whicha lattice model and statistical approaches are used to simulate the initial heterogeneity ofmaterial properties. The size of displacement-load step is adap-tively determined so that only fewelements would fail in each load step. When the tensile principal strain in an element exceeds theultimate strain of this element, the element is considered broken and its Young's modulus is set tobe very low. Some important behaviors of heterogeneous brittle materials are indicated using thiscode. Load-displacement curves and figures of three-dimensional fracture patterns are alsonumerically obtained, which are similar to those observed in laboratory tests.

Numerical simulation of flow interference between four cylinders in in-line square arrangement

In the present paper,two-and three-dimensional numerical simulations of the flow interference between four cylinders in an in-line square arrangement at Re = 200 are performed.Assisted with the two-dimensional(2-D) numerical simulation,the mean and fluctuating forces,Strouhal number(St) and vortex shedding pattern in the wake for each cylinder were analyzed with the spacing ratio(L /D) ranging from 1.5 to 6.0.It was found that,four different vortex modes(viz.,flip-flopping,shielding anti-phase-synchronized,in-phasesynchronized and anti-phase-synchronized) gradually appear with the increase of the L/D ratio.The average drag coefficient of the upstream cylinders is larger than that of the downstream cylinders,while the downstream cylinders usually undergo serious fluctuating forces.When the L/D ratio ranges from 3.0 to 4.0,the dominant frequency of the drag coefficient is equal to the value of St of upstream cylinders.This indicates that a simultaneous resonance in the in-flow and cross-flow directions may occur for some single structures of a multi-body oscillating system.For the 3-D numerical simulation,the L/D and aspect ratios are kept constant as 5.0 and 10,respectively.It was found that some vortices are formed in the wake of the upstream cylinders.Besides,with the same spacing ratio,the calculated drag coefficient and lift coefficient fluctuation are slightly larger than the 2-D results,but with a phase difference.

Three-dimensional numerical manifold method for heat conduction problems with a simplex integral on the boundary

The three-dimensional numerical manifold method(3D-NMM),which is based on the derivation of Galerkin's variation,is a powerful calculation tool that uses two cover systems.The 3D-NMM can be used to handle continue-discontinue problems and extend to THM coupling.In this study,we extended the 3D-NMM to simulate both steady-state and transient heat conduction problems.The modelling was carried out using the raster methods(RSM).For the system equation,a variational method was employed to drive the discrete equations,and the crucial boundary conditions were solved using the penalty method.To solve the boundary integral problem,the face integral of scalar fields and two-dimensional simplex integration were used to accurately describe the integral on polygonal boundaries.Several numerical examples were used to verify the results of 3D steady-state and transient heat-conduction problems.The numerical results indicated that the 3D-NMM is effective for handling 3D both steadystate and transient heat conduction problems with high solution accuracy.

Simulations of the Transient Flow Generated from a Started Flat Plate

Transient operations are commonly founded in fluid machineries such as the starting, stopping, and variations of rotor speeds, etc. Flow generated from a started fiat plate is of fundamental importance. Experiments have been done to observe the flow evolution in current researches. And in order to explore the flow in more detailed scale, some vortex methods with high resolution and other numerical methods were developed to solve various related problems by some researchers. But the promotion of vortex method to engineering application is rare due to its complexity and difficulty in specifying the boundary conditions. In order to build up a method of numerical study for such problems, a simplified model is built up with a flat plate. The development of two-dimensional viscous incompressible flow generated from an impulsively started and uniformly accelerated infinitesimally thin flat plate is simulated numerically. A dynamic mesh（DM） method based on the spring analogue and local remeshing is applied to realize the mesh motion caused by the started plate. Researches show that the mesh quality will decline under large grid shear force during the updating process. To conquer this problem, a region near the plate is separated to guarantee the mesh quality at location of interest which is the innovation of the present paper. All computations at least cover a period during which the plate translates 6 times its length. The simulated instantaneous velocity profiles, flow structures and drag coefficients under several Reynolds numbers （20 ≤ Re ≤ 126） and accelerations （20 m/s2≤ a ≤ 152 m/s2） are presented and compared with existing results in literatures. Comparisons are found to be satisfactory, confirming the validity of the current proposed method（region separated DM）. The proposed DM method is firstly used to study the transient flow generated from a started flat plate and can be used in further study of transient characteristics during transient operations of turbo machineries.

On Transient Simulation of Field Equations

We investigate under what conditions transient simulation could be used to integrate backward in time so that the initial field could be recovered from later histories. In this paper we use realistic examples and find that, in long histories, traces of the initial field would be present only in the exact analytical solutions. We conclude that the recovery of initial field is possible only if the equations could be solved analytically or only short time periods are involved. In practice, it is not possible to detect those traces by measurements or observations. If numerical procedures are used, truncation and discretization errors are always present. Fine-tuning of system parameters used or transforming time into another pseudo time frame may allow numerical integration to be carried out backward in time. But numerical instability is still a problem. Large spurious increases found by numerical procedures are most likely due to numerical inaccuracy and instability.

Modeling and Simulation for Transient Thermal Analyses Using a Voltage-in-Current Latency Insertion Method

This article presents a modeling and simulation method for transient thermal analyses of integrated circuits(ICs)using the original and voltage-in-current(VinC)latency insertion method(LIM).LIM-based algorithms are a set of fast transient simulation methods that solve electrical circuits in a leapfrog updating manner without relying on large matrix operations used in conventional Simulation Program with Integrated Circuit Emphasis(SPICE)-based methods which can significantly slow down the solution process.The conversion from the thermal to electrical model is performed first by using the analogy between heat and electrical conduction.Since electrical inductance has no thermal equivalence,a modified VinC LIM formulation is presented which removes the requirement of the insertion of fictitious inductors.Numerical examples are presented,which show that the modified VinC LIM formulation outperforms the basic LIM formulation,in terms of both stability and accuracy in the transient thermal simulation of ICs.