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.

Effects of mesh style and grid convergence on numerical simulation accuracy of centrifugal pump

In order to evaluate the effects of mesh generation techniques and grid convergence on pump performance in centrifugal pump model, three widely used mesh styles including structured hexahedral, unstructured tetrahedral and hybrid prismatic/tetrahedral meshes were generated for a centrifugal pump model. And quantitative grid convergence was assessed based on a grid convergence index(GCI), which accounts for the degree of grid refinement. The structured, unstructured or hybrid meshes are found to have certain difference for velocity distributions in impeller with the change of grid cell number. And the simulation results have errors to different degrees compared with experimental data. The GCI-value for structured meshes calculated is lower than that for the unstructured and hybrid meshes. Meanwhile, the structured meshes are observed to get more vortexes in impeller passage.Nevertheless, the hybrid meshes are found to have larger low-velocity area at outlet and more secondary vortexes at a specified location than structured meshes and unstructured meshes.

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.

Overlapping Grid Technique for Numerical Simulation of a Fast-Cruising Catamaran Fitted with Active T-Foils

In marine engineering,appendages such as fin stabilizers and/or T-foils are made to rotate and to reduce the motion of ships.Research on the hydrodynamics of ships fitted with active appendages has significantly improved the design and control of such ships.However,most studies focus on fixed rather than rotating appendages,thereby ignoring the hydrodynamic unsteadiness of active appendages.To enhance the reliability and precision of the numerical simulations,we propose the use of overlapping grids for simulating advanced catamarans fitted with a pair of rotating T-foils under each bow.The fundamental purpose of the overlapping grid technique is to realize information exchange via regional overlap sharing in each subdomain of the computing domain,instead of using the method of boundary sharing,thus greatly alleviating the difficulty of generating the subdomain grid;moreover,the technique guarantees the quality of the subdomain grid.Within the main computational domain,a subdomain was allocated to accommodate the T-foil.Overlapping meshes near the interface between the two domains enable information flow during the simulation;the overlapping grids are updated at every iteration step because the subdomain rotates.The instantaneous trim and sinkage responses of the catamaran to the T-foil rotation were reproduced.From the moment the active T-foil stopped moving,there was no change in the ship’s sailing attitude,indicating that the response was in real time.By comparing with EFD data,the numerical results showed reasonable agreement,indicating the feasibility and effectiveness of the technique in simulating the hydrodynamics of ships fitted with active appendages.

The numerical simulation for a 3D two-phase anisotropic medium based on BISQ model

Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components （2D2C） or two-dimensions and three-components （2D3C）. There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.

The research analysis of aerodynamic numerical simulation of grid fin

This paper presents the results of an investigation to use arc-length mesh generation and finite volume TVD scheme to calculate Euler equations for predicting the effect of geometry parameters in reducing the drag force and improving the lift-drag ratio of grid fin in the supersonic flow regime. The effects of frame and web, whose cross section shape and thickness and spacing,on the aerodynamic character of the grid fin were studied. Calculations were made at Mach 2.5 and several angles of attack. The results were validated by comparing the computed aerodynamic coefficients against wind tunnel experimental data. Good agreement was found between computed and experimental results. The computed results suggest that parameters of the grid fin's frame have the greatest effect on the grid fin aerodynamic character, especially on its drag force. It was concluded proper choice of appropriate grid fin geometry parameters could reduce the drag force and improve the lift-drag ratios.

A New Method to Determine the Grid Directions in Reservoir Numerical Simulation

Grid direction selection and grid size design are two important elements that need to be considered in the grid direction design in reservoir numerical simulation. Reservoir engineers normally utilize geological data (such as the distribution of fractures, low permeability zones, faults and major stress) and simulation experiences to design the grid direction of simulation model qualitatively. The research of the paper indicates that the key to determine the grid direction is to determine the principal permeability direction. Under the circumstances of few static materials, a new grid direction determination method has been developed by using field data (well location map and inter-well permeability) on the bases of Darcy’s law and tensor analysis theory. The grid direction of WZ11-7 Oilfield simulation model has been determined using four production wells and two production zones (L1 and L3) in WZ11-7-2 well group, the results are in conformity with the geological studied major stress. Therefore, this method can give insights into the numerical simulation study.

2-D Current Field Numerical Simulation Integrating Yangtze Estuary with Hangzhou Bay

In this paper, integrating the Yangtze Estuary with the Hangzhou Bay, a 2-D velocity field model is established. In the model, fine self-adaptive grids are employed to adapt to the complicated coastal shape. The hydrodynamic equations satisfied by two contravariant components of velocity vector and surface elevation in non-orthogonal curvilinear coordinates are used. In each momentum equation the coefficients before the two partial derivatives of surface elevation with respect to variables of alternative direction coordinates have different orders of magnitude, i. e., the derivative with the larger coefficient may play a more important role than that with the smaller one. With this advantage, the ADI scheme can be easily employed. The hydrodynamic factors include tidal current, river runoff and wind-induced current. In terms of tidal current, seven main constituents in the area are considered in the open boundaries. The verifications of surface elevation process and current velocity process in the spring tide and in the neap tide show that the model can preferably reflect current fields in the area. Through the simulation of Lagrangian residual current fields in summer and in winter, the paths of the exchange of water and sediment between the Yangtze Estuary and the Hangzhou Bay are elementarily discussed.

A numerical simulation of seepage structure surface and its feasibility

According to Cubic law and incompressible fluid law of mass conservation, the seepage character of the fracture surface was simulated with the simulation method of fractal theory and random Brown function. Furthermore, the permeability coefficient of the single fracture was obtained. In order to test the stability of the method, 500 simulations were conducted on each different fractal dimension. The simulated permeability coefficient was analyzed in probability density distribution and probability cumulative distribution statistics. Statistics showed that the discrete degree of the permeability coefficient increases with the increase of the fractal dimension. And the calculation result has better stability when the fractal dimension value is relatively small. According to the Bayes theory, the characteristic index of the permeability coefficient on fractal dimension P（Dfi｜ Ri） is established. The index, P（Dfi｜ Ri）, shows that when the simulated permeability coefficient is relatively large, it can clearly represent the fractal dimension of the structure surface, the probability is 82%. The calculated results of the characteristic index verify the feasibility of the method.

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.

Numerical simulation of wind-driven circulation and pollutant transport in Taihu Lake based on a quadtree grid

In this study,a two-dimensional flow-pollutant coupled model was developed based on a quadtree grid.This model was established to allow the accurate simulation of wind-driven flow in a large-scale shallow lake with irregular natural boundaries when focusing on important smallscale localized flow features.The quadtree grid was created by domain decomposition.The governing equations were solved using the finite volume method,and the normal fluxes of mass,momentum,and pollutants across the interface between cells were computed by means of a Godunov-type Osher scheme.The model was employed to simulate wind-driven flow in a circular basin with non-uniform depth.The computed values were in agreement with analytical data.The results indicate that the quadtree grid has fine local resolution and high efficiency,and is convenient for local refinement.It is clear that the quadtree grid model is effective when applied to complex flow domains.Finally,the model was used to calculate the flow field and concentration field of Taihu Lake,demonstrating its ability to predict the flow and concentration fields in an actual water area with complex geometry.

Numerical simulation research of dimension variation in heavy plate during quenching and tempering process

The dimension variation of B610E steel during the quenching and tempering process was simulated by ABAQUS,the coefficient of heat transfer during quenching was verified by the buried thermocouple test and the trend of the dimension variation during the quenching process was also calculated by ABAQUS. It was shown by the comparison of simulated results and industrially measured results that only the simulation of thickness and width variation was accurate, while length variation needed further simulation. Besides, the dimension variation trend was identical with the measured results.

Three-dimensional Numerical Simulation of the Movement of the Flexible Body under Different Constraints

For the large deformation of the flexible body may cause the fluid grid distortion,which will make the numerical calculation tedious,even to end,the numerical simulation of the flexible body coupling with the fluid is always a tough problem.In this paper,the flexible body is under two kinds of constrained conditions and the ratio of length-diameter is 1:30.The Reynolds number of the airflow is 513,belonging to the area of low Reynolds number.The control equations of the coupling of flexible body with airflow are built and the adaptive grid control method is adopted to conduct the three-dimensional numerical simulation of the movement of the flexible body.The numerical results show that it is possible to simulate the characteristics of the flexible body's movement in the low Reynolds number airflow when the appropriate control equations are modeled and suitable equation-solving method is adopted.Unconstrained flexible body would turn over forward along the airflow's diffusion direction,while constrained flexible body in the flow field will make periodic rotation motion along the axis of the flexible body,and the bending deformation is more obvious than that of unconstrained flexible body.The preliminary three-dimensional numerical simulation can provide references for further research on the characteristics of the yam movement in high Reynolds number airflow.

Grid Convergence Property of Three-Dimensional Measurement-Integrated Simulation for Unsteady Flow behind a Square Cylinder with Karman Vortex Street

The purpose of this study was to clarify grid convergence property of three-dimensional measurement-integrated (3D-MI) simulation for a flow behind a square cylinder with Karman vortex street. Measurement-integrated (MI) simulation is a kind of the observer in the dynamical system theory by using CFD scheme as a mathematical model of the system. In a former study, two-dimensional MI (2D-MI) simulation with a coarse grid system showed a fairly good result in comparison with a 2D ordinary (2D-O) simulation, but the results were degraded with grid refinement. In this study, 3D-MI simulation and three-dimensional ordinary (3D-O) simulation were performed with three grid systems of different grid resolutions, and their grid convergence properties were compared. As a result, all 3D-MI simulations reproduced the vortex shedding frequency identical to that of the experiment, and the flow fields obtained were very close, within 5% difference between the results, while the results of the 3D-O simulations showed variation of the solution under convergence. It is shown that the grid convergence property of 3D-MI simulation is monotonic and better than that of 3D-O simulation, whereas those of 2D-O and 2D-MI simulations for streamwise velocity fluctuation are divergent. The solution of 3D-MI simulation with a relatively coarse grid system properly reproduces the basic three-dimensional structure of the wake flow as well as the drag and lift coefficients.

Productivity simulation of hydraulically fractured wells based on hybrid local grid refinement and embedded discrete fracture model

Using current Embedded Discrete Fracture Models(EDFM) to predict the productivity of fractured wells has some drawbacks, such as not supporting corner grid, low precision in the near wellbore zone, and disregarding the heterogeneity of conductivity brought by non-uniform sand concentration. An EDFM is developed based on the corner grid, which enables high efficient calculation of the transmissibility between the embedded fractures and matrix grids, and calculation of the permeability of each polygon in the embedded fractures by the lattice data of the artificial fracture aperture. On this basis, a coupling method of local grid refinement(LGR) and embedded discrete fracture model is designed, which is verified by comparing the calculation results with the Discrete Fracture Network(DFN) method and fitting the actual production data of the first hydraulically fractured well in Iraq. By using this method and orthogonal experimental design, the optimization of the parameters of the first multi-stage fractured horizontal well in the same block is completed. The results show the proposed method has theoretical and practical significance for improving the adaptability of EDFM and the accuracy of productivity prediction of fractured wells, and enables the coupling of fracture modeling and numerical productivity simulation at reservoir scale.

Numerical Analysis on the Effects of Submerged Depth of the Grid and Direction of Incident Wave on Gravity Cage

In this paper, the numerical model of the net cage with the grid mooring system in waves is set up by the lumped mass method and rigid kinematics theory, and then the motion equations of floating system, net system, mooring system, and floaters are solved by the Runge-Kutta fifth-order method. For the verification of the numerical model, a series of physical model tests have been carried out. According to the comparisons between the simulated and experimental results, it can be found that the simulated and experimental results agree well in each condition. Then, the effects of submerged depth of grid and direction of incident wave propagation on hydrodynamic behaviors of the net cage are analyzed. According to the simulated results, it can be found that with the increase of submerged depth of grid, the forces acting on mooring lines and bridle lines increase, while the forces on grid lines decrease; the horizontal motion amplitudes of floating collar decrease obviously, while the vertical motion amplitudes of floating collar change little. When the direction of incident wave propagation changes, forces on mooring lines and motion of net cage also change accordingly. When the propagation direction of incident wave changes from 0° to 45°, forces on the main ropes and bridle ropes increase, while the forces on the grid ropes decrease. With the increasing propagation direction of incident wave, the horizontal amplitude of the forces collar decreases, while the vertical amplitude of the floating collar has little variation.

Petrel2ANSYS: Accessible software for simulation of crustal stress fields using constraints provided by multiple 3D models employing different types of grids

Crustal stresses play an important role in both exploration and development in the oil and gas industry.However,it is difficult to simulate crustal stress distributions accurately,because of the incompatibilities that exist among different software.Here,a series of algorithms is developed and integrated in the Petrel2ANSYS to carry out two-way conversions between the 3D attribute models that employ corner-point grids used in Petrel and the 3D finite-element grids used in ANSYS.Furthermore,a modified method of simulating stress characteristics and analyzing stress fields using the finite-element method and multiple finely resolved 3D models is proposed.Compared to the traditional finite-element simulation-based approach,which involves describing the heterogeneous within a rock body or sedimentary facies in detail and simulating the stress distribution,the single grid cell-based approach focuses on a greater degree on combining the rock mechanics described by 3D corner-point grid models with the finely resolved material characteristics of 3D finite-element models.Different models that use structured and unstructured grids are verified in Petrel2ANSYS to assess the feasibility.In addition,with minor modifications,platforms based on the present algorithms can be extended to other models to convert corner-point grids to the finite-element grids constructed by other software.

Computational fluid dynamics simulations of respiratory airflow in human nasal cavity and its characteristic dimension study

To study the airflow distribution in human nasal cavity during respiration and the characteristic parameters of nasal structure, three-dimensional, anatomically accurate representations of 30 adult nasal cavity models were recons- tructed based on processed tomography images collected from normal people. The airflow fields in nasal cavities were simulated by fluid dynamics with finite element software ANSYS. The results showed that the difference of human nasal cavity structure led to different airflow distribution in the nasal cavities and variation of the main airstream passing through the common nasal meatus. The nasal resistance in the regions of nasal valve and nasal vestibule accounted for more than half of the overall resistance. The characteristic model of nasal cavity was extracted on the basis of characteristic points and dimensions deduced from the original models. It showed that either the geometric structure or the airflow field of the two kinds of models was similar. The characteristic dimensions were the characteristic parameters of nasal cavity that could properly represent the original model in model studies on nasal cavity.

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.

Slope reliability analysis based on Monte Carlo simulation and sparse grid method

In order to solve the problem of the reliability of slope engineering due to complex uncertainties, the Monte Carlo simulation method is adopted. Based on the characteristics of sparse grid, an interpolation algorithm, which can be applied to high dimensional problems, is introduced. A surrogate model of high dimensional implicit function is established, which makes Monte Carlo method more adaptable. Finally, a reliability analysis method is proposed to evaluate the reliability of the slope engineering, and is applied in the Sau Mau Ping slope project in Hong Kong. The reliability analysis method has great theoretical and practical significance for engineering quality evaluation and natural disaster assessment.