A semi‐implicit scheme for large Eddy simulation of piston engine flow and combustion

A semi‐implicit scheme is presented for large eddy simulation of turbulent reactive flow and combustion in reciprocating piston engines. First, the governing equations in a deforming coordinate system are formulated to accommodate the moving piston. The numerical scheme is made up of a fourth‐order central difference for the diffusion terms in the transport equations and a fifth‐order weighted essentially nonoscillatory (WENO) scheme for the convective terms. A second‐ order Adams–Bashforth scheme is used for time integration. For higher density ratios, it is combined with a predictor–corrector scheme. The numerical scheme is explicit for time integration of the transport equations, except for the continuity equation which is used together with the momentum equation to determine the pressure field and velocity field by using a Poisson equation for the pressure correction field. The scheme is aimed at the simulation of low Mach number flows typically found in piston engines. An efficient multigrid method that can handle high grid aspect ratio is presented for solving the pressure correction equation. The numerical scheme is evaluated on two test engines, a laboratory four‐stroke engine with rectangular‐shaped engine geometry where detailed velocity measurements are available, and a modified truck engine with practical cylinder geometry where lean ethanol/air mixture is combusted under a homogeneous charge compression ignition (HCCI) condition. Copyright © 2012 John Wiley & Sons, Ltd.     

Efficient discretization of pressure‐correction equations on non‐orthogonal grids

An alternative discretization of pressure‐correction equations within pressure‐correction schemes for the solution of the incompressible Navier–Stokes equations is introduced, which improves the convergence and robustness properties of such schemes for non‐orthogonal grids. As against standard approaches, where the non‐orthogonal terms usually are just neglected, the approach allows for a simplification of the pressure‐correction equation to correspond to 5‐point or 7‐point computational molecules in two or three dimensions, respectively, but still incorporates the effects of non‐orthogonality. As a result a wide range (including rather high values) of underrelaxation factors can be used, resulting in an increased overall performance of the underlying pressure‐correction schemes. Within this context, a second issue of the paper is the investigation of the accuracy to which the pressure‐correction equation should be solved in each pressure‐correction iteration. The scheme is investigated for standard test cases and, in order to show its applicability to practical flow problems, for a more complex configuration of a micro heat exchanger. Copyright © 2003 John Wiley & Sons, Ltd.     

注塑成型三维流动的数值模拟

建立了非等温、粘性、不可压缩、非牛顿流体流动的控制方程。为了避免同时求解耦合的压力场、速度场,本文通过修改Galerkin方法的变分方程,导出了关于压力场的拟Poisson方程,用迭代法独立地求解连续性方程、动量方程,并进行速度一粘度迭代求出最终的压力场、速度场。由于直接使用Galerkin方法求解能量方程容易引起温度场的振荡,本文采用隐式格式及“上风”法离散能量方程,用超松驰迭代法求解温度场的代数方程组。比较了模拟结果与等温管道流动的解析解及法兰的实际注射结果,算例表明本文方法可以预测注射成型流动过程中的一些重要特征。与传统Galerkin方法相比,本文方法可以减少内存,提高数值方法的稳定性。     

Solution of the discretized incompressible Navier-Stokes equations with the GMRES method

We describe some experiences using interative solution methods of GMRES type to solve the discretized Navier-Stokes equations. The discretization combined with a pressure correction scheme leads to two different systems of equations: the momentum equations and the pressure equation. It appears that a fast solution method for the pressure equation is obtained by applying the recently proposed GMRESR method, or GMRES combined with a MILU preconditioner. The diagonally scaled momentum equations are solved by GMRES(m), a restarted version of GMRES.     

A simplified marker and cell method for unsteady flows on non-staggered grids

The SMAC (simplified marker and cell) time-advancing method for solving the unsteady incompressible Navier-Stokes equations on non-staggered grids is developed in generalized co-ordinate systems. The primitive variable formulation uses Cartesian velocities and pressure, all defined at the centre of the control volume, as the dependent variables. A special elliptic flux correction at the faces of the finite volume is utilized in discretizing the continuity equation to suppress pressure oscillations. The test flows considered are a polar cavity flow starting from rest and the flow around a circular cylinder. The numerical results are compared with experimental results and results obtained by the well-known SIMPLEC and PISO methods. The comparisons show that the elliptic flux correction technique works well in suppressing pressure oscillations and that the SMAC method is more efficient than the SIMPLEC and PISO methods for both steady and unsteady flows.     

Numerical analysis of pressure oscillations over axisymmetric spiked blunt bodies at Mach 6.80

The pressure oscillations over a forward facing spike attached to an axisymmetric blunt body are simulated by solving time-dependent compressible Navier–Stokes equations. The governing fluid flow equations are discretized in spatial coordinates employing a finite volume approach which reduces the equations to semidiscretized ordinary differential equations. Temporal integration is performed using the two-stage Runge–Kutta time stepping scheme. A global time step is used to obtain a time-accurate numerical solution. The numerical computation is carried out for a freestream Mach number of 6.80 and for spike length to hemispherical diameter ratios of 0.5, 1.0 and 2.0. The flow features around the spiked blunt body are characterized by a conical shock wave emanating from the spike tip, a region of separated flow in front of the hemispherical cap, and the resulting reattachment shock wave. Comparisons of the numerical results are made with the available experimental results, such as schlieren pictures and the surface pressure distribution along the spiked blunt body. They are found to be in good agreement. Spectral analysis of the computed pressure oscillations are performed employing fast Fourier transforms. The surface pressure oscillations over the spike and phase plots exhibit a behaviour analogous to that of the Van der Pol equation for a self-sustained oscillatory flow. Received 28 February 2001 / Accepted 17 January 2002     

A composite multigrid method for calculating unsteady incompressible flows in geometrically complex domains

A time-accurate, finite volume method for solving the three-dimensional, incompressible Navier-Stokes equations on a composite grid with arbitrary subgrid overlapping is presented. The governing equations are written in a non-orthogonal curvilinear co-ordinate system and are discretized on a non-staggered grid. A semi-implicit, fractional step method with approximate factorization is employed for time advancement. Multigrid combined with intergrid iteration is used to solve the pressure Poisson equation. Inter-grid communication is facilitated by an iterative boundary velocity scheme which ensures that the governing equations are well-posed on each subdomain. Mass conservation on each subdomain is preserved by using a mass imbalance correction scheme which is secondorder-accurate. Three test cases are used to demonstrate the method's consistency, accuracy and efficiency.     

A numerical method for 3D viscous incompressible flows using non-orthogonal grids

This paper presents a numerical method for fluid flow in complex three-dimensional geometries using a body-fitted co-ordinate system. A new second-order-accurate scheme for the cross-derivative terms is proposed to describe the non-orthogonal components, allowing parts of these terms to be treated implicitly without increasing the number of computational molecules. The physical tangential velocity components resulting from the velocity expansion in the unit tangent vector basis are used as dependent variables in the momentum equations. A coupled equation solver is used in place of the complicated pressure correction equation associated with grid non-orthogonality. The co-ordinate-invariant conservation equations and the physical geometric quantities of control cells are used directly to formulate the numerical scheme, without reference to the co-ordinate derivatives of transformation. Several two- and three-dimensional laminar flows are computed and compared with other numerical, experimental and analytical results to validate the solution method. Good agreement is obtained in all cases.     

Numerical simulation of oscillations of a monodisperse gas-particle mixture in a nonlinear wave field

This paper presents a numerical simulation procedure for the dynamics of a monodisperse gas-particle mixture in the nonlinear wave field of an acoustic resonator using a two-temperature two-velocity model ignoring phase transitions, particle collision, and possible coagulation. It is assumed that viscosity is present only in the carrier medium described by the Navier-Stokes equations for a compressible gas. The dispersed phase is described by the equation of conservation of mass, momentum, and energy. A monotonic solution is obtained by solving the equations of motion for the carrier medium and dispersed phase in generalized moving coordinates using the explicit McCormack method with splitting in the spatial directions and a conservative correction scheme. The method can be used to study nonlinear oscillations of two-phase mixtures in the vicinity of the first three eigenfrequencies in a flat channel.     

A finite element analysis of laminar unsteady flow of viscoelastic fluids through channels with non-uniform cross-sections

The effects of non-Newtonian behaviour of a fluid and unsteadiness on flow in a channel with non-uniform cross-section have been investigated. The rheological behaviour of the fluid is assumed to be described by the constitutive equation of a viscoelastic fluid obeying the Oldroyd-B model. The finite element method is used to analyse the flow. The novel features of the present method are the adoption of the velocity correction technique for the momentum equations and of the two-step explicit scheme for the extra stress equations. This approach makes the computational scheme simple in algorithmic structure, which therefore implies that the present technique is capable of handling large-scale problems. The scheme is completed by the introduction of balancing tensor diffusivity (wherever necessary) in the momentum equations. It is important to mention that the proper boundary condition for pressure (at the outlet) has been developed to solve the pressure Poisson equation, and then the results for velocity, pressure and extra stress fields have been computed for different values of the Weissenberg number, viscosity due to elasticity, etc. Finally, it is pertinent to point out that the present numerical scheme, along with the proper boundary condition for pressure developed here, demonstrates its versatility and suitability for analysing the unsteady flow of viscoelastic fluid through a channel with non-uniform cross-section.     

Development of a hybrid finite volume/element solver for incompressible flows

In this paper, we report our development of an implicit hybrid flow solver for the incompressible Navier–Stokes equations. The methodology is based on the pressure correction or projection method. A fractional step approach is used to obtain an intermediate velocity field by solving the original momentum equations with the matrix‐free implicit cell‐centred finite volume method. The Poisson equation derived from the fractional step approach is solved by the node‐based Galerkin finite element method for an auxiliary variable. The auxiliary variable is closely related to the real pressure and is used to update the velocity field and the pressure field. We store the velocity components at cell centres and the auxiliary variable at cell vertices, making the current solver a staggered‐mesh scheme. Numerical examples demonstrate the performance of the resulting hybrid scheme, such as the correct temporal convergence rates for both velocity and pressure, absence of unphysical pressure boundary layer, good convergence in steady‐state simulations and capability in predicting accurate drag, lift and Strouhal number in the flow around a circular cylinder. Copyright © 2007 John Wiley & Sons, Ltd.     

A PRESSURE-CORRECTION METHOD FOR THE SOLUTION OF INCOMPRESSIBLE VISCOUS FLOWS ON UNSTRUCTURED GRIDS

An unstructured grid, finite volume method is presented for the solution of two-dimensional viscous, incompressible flow. The method is based on the pressure-correction concept implemented on a semi-staggered grid. The computational procedure can handle cells of arbitrary shape, although solutions presented herein have been obtained only with meshes of triangular and quadrilateral cells. The discretization of the momentum equations is effected on dual cells surrounding the vertices of primary cells, while the pressure-correction equation applies to the primary-cell centroids and represents the conservation of mass across the primary cells. A special interpolation scheme s used to suppress pressure and velocity oscillations in cases where the semi-staggered arrangement does not ensure a sufficiently strong coupling between pressure and velocity to avoid such oscillations. Computational results presented for several viscous flows are shown to be in good agreement with analytical and experimental data reported in the open literature.     

Hybrid finite element/volume method for 3D incompressible flows with heat transfer

An implicit hybrid finite element (FE)/volume solver has been extended to incompressible flows coupled with the energy equation. The solver is based on the segregated pressure correction or projection method on staggered unstructured hybrid meshes. An intermediate velocity field is first obtained by solving the momentum equations with the matrix-free implicit cell-centred finite volume (FV) method. The pressure Poisson equation is solved by the node-based Galerkin FE method for an auxiliary variable. The auxiliary variable is used to update the velocity field and the pressure field. The pressure field is carefully updated by taking into account the velocity divergence field. Our current staggered-mesh scheme is distinct from other conventional ones in that we store the velocity components at cell centres and the auxiliary variable at vertices. The Generalized Minimal Residual (GMRES) matrix-free strategy is adapted to solve the governing equations in both FE and FV methods. The presented 2D and 3D numerical examples show the robustness and accuracy of the numerical method.     

Numerical simulation of viscous flow around unrestrained cylinders

A 2-D analysis is made for the dynamic interactions between viscous flow and one or more circular cylinders. The cylinder is free to respond to the fluid excitation and its motions are part of the solution. The numerical procedure is based on the finite volume discretization of the Navier–Stokes equations on adaptive tri-tree grids which are unstructured and nonorthogonal. Both a fully implicit scheme and a semi-implicit scheme in the time domain have been used for the momentum equations, while the pressure correction method based on the SIMPLE technique is adopted to satisfy the continuity equation. A new upwind method is developed for the triangular and unstructured mesh, which requires information only from two neighbouring cells but is of order of accuracy higher than linear. A new procedure is also introduced to deal with the nonorthogonal term. The pressure on the body surface required in solving the momentum equation is obtained through the Poisson equation in the local cell. Results including flow field, pressure distribution and force are provided for fixed single and multiple cylinders and for an unrestrained cylinder in steady incoming flow with Reynolds numbers at 200 and 500 and in unsteady flow with Keulegan–Carpenter numbers at 5 and 10.     

A PRESSURE-BASED ALGORITHM FOR HIGH-SPEED TURBOMACHINERY FLOWS

The steady state Navier–Stokes equations are solved in transonic flows using an elliptic formulation. A segregated solution algorithm is established in which the pressure correction equation is utilized to enforce the divergence-free mass flux constraint. The momentum equations are solved in terms of the primitive variables, while the pressure correction field is used to update both the convecting mass flux components and the pressure itself. The velocity components are deduced from the corrected mass fluxes on the basis of an upwind-biased density, which is a mechanism capable of overcoming the ellipticity of the system of equations, in the transonic flow regime. An incomplete LU decomposition is used for the solution of the transport-type equations and a globally minimized residual method resolves the pressure correction equation. Turbulence is resolved through the k–ε model. Dealing with turbomachinery applications, results are presented in two-dimensional compressor and turbine cascades under design and off-design conditions. © 1997 John Wiley & Sons, Ltd.     

On the stabilization of unphysical pressure oscillations in MPS method simulations

In this paper, pressure stability through the suppression of high‐frequency pressure oscillations in the moving particle semi‐implicit (MPS) method is presented. To obtain a stable pressure field, we improve the free‐surface particle search algorithm. Pressure stability follows from the suppression of high‐frequency pressure oscillations due to a correction in the Laplacian operator of the Poisson pressure equation and from the correction of the pressure gradient operator. The three proposed modifications are applied gradually and compared with the MPS method to show the improvements in the hydrostatic pressure and dam‐breaking problems. To validate the suppression of the high‐frequency numerical pressure oscillations, modified MPS methods with and without a removable wall are compared with published dam‐breaking experiment pressure measurements. Copyright © 2016 John Wiley & Sons, Ltd.     

A split operator scheme for ocean wave simulation

A coupled discrete spectral model was developed for the prediction of ocean waves by solving the energy conservation equation of the two-dimensional wave spectrum. The model includes the dispersion correction terms in the governing equation to account for the dispersive effect due to the frequency-dependent velocities of waves. A split operator scheme is used to deal with the numerical problems arising from different terms of the governing equation. The advection terms are solved by the proven accurate minimax characteristics method to avoid excessive numerical damping or oscillations. The dispersion correction terms are solved by central differencing. The source and sink terms are solved by a quasi-second-order explicit scheme with limitation on energy growth per time step to allow the use of a large time step. The model was verified by ideal test cases and wave-hindcasting studies under typhoon conditions in the South China Sea near Hong Kong.     

求解双曲守恒律方程的WENO型熵相容格式

通过在单元交界面处进行高阶WENO重构,得到了一种求解双曲型守恒律方程的WENO型熵相容格式。用该格式对一维Burgers方程和Euler方程进行数值模拟,结果表明,该格式具有高精度、基本无振荡性等特点。     

Some improvements on free surface simulation by the particle finite element method

The Lagrangian method has become increasingly popular in numerical simulation of free surface problems. In this paper, after a brief review of a recent Lagrangian method, namely the particle finite element method, some issues are discussed and some improvements are made. The least‐square finite element method is adopted to simplify the solving of the Navier–Stokes equations. An adaptive time method is derived to obtain suitable time steps. A mass correction procedure is imported to improve the mass conservation in long time calculations and time discretization scheme is adopted to decrease the pressure oscillations during the calculations. Finally, the method is used to simulate a series of examples and the results are compared with the commercial FLOW3D code. Copyright © 2008 John Wiley & Sons, Ltd.     

Forced oscillations of a gas bubble in a spherical volume of a compressible liquid

A spherically symmetric problem of oscillations of a single gas bubble at the center of a spherical flask filled with a compressible liquid under the action of pressure oscillations on the flask wall is considered. A system of differential-difference equations is obtained that extends the Rayleigh-Plesset equation to the case of a compressible liquid and takes into account the pressure-wave reflection from the bubble and the flask wall. A linear analysis of solutions of this system of equations is performed for the case of harmonic oscillations of the bubble. Nonlinear resonance oscillations and nearly resonance nonharmonic oscillations of the bubble caused by harmonic pressure oscillations on the flask wall are analyzed. Ufa Scientific Center, Russian Academy of Sciences, Ufa 450000. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 2, pp. 111–118, March–April, 1999.