可压流数值模拟中当地DFD方法的改进和应用

阐述了求解守恒型Euler 方程的当地DFD (Domain-Free Discretization) 方法的改进和应用。DFD 离散策略的核心,是解域内点上控制方程的离散形式可与解域外的一些点相关。通过边界附近解的近似形式,外部相关点上的流动变量值得到确定并强加相应的边界条件。与最初的当地DFD方法不同,在解的近似形式构建中,采用了CCST技术 (Curvature-Corrected Symmetry Technique),因此外部相关点上的密度和切向速度分别由等熵和等总焓关系确定。空间离散采用Galerkin 有限体积格式。最后,给出了固定和运动物体可压缩绕流的数值模拟结果,以验证改进的当地DFD方法的可靠性和数值解精度的提高。     

Higher‐order bounded differencing schemes for compressible and incompressible flows

In recent years, three higher‐order (HO) bounded differencing schemes, namely AVLSMART, CUBISTA and HOAB that were derived by adopting the normalized variable formulation (NVF), have been proposed. In this paper, a comparative study is performed on these schemes to assess their numerical accuracy, computational cost as well as iterative convergence property. All the schemes are formulated on the basis of a new dual‐formulation in order to facilitate their implementations on unstructured meshes. Based on the proposed dual‐formulation, the net effective blending factor (NEBF) of a high‐resolution (HR) scheme can now be measured and its relevance on the accuracy and computational cost of a HR scheme is revealed on three test problems: (1) advection of a scalar step‐profile; (2) 2D transonic flow past a circular arc bump; and (3) 3D lid‐driven incompressible cavity flow. Both density‐based and pressure‐based methods are used for the computations of compressible and incompressible flow, respectively. Computed results show that all the schemes produce solutions which are nearly as accurate as the third‐order QUICK scheme; however, without the unphysical oscillations which are commonly inherited from the HO linear differencing scheme. Generally, it is shown that at higher value of NEBF, a HR scheme can attain better accuracy at the expense of computational cost. Copyright © 2006 John Wiley & Sons, Ltd.     

Characteristic‐based inlet and outlet boundary conditions for incompressible flows

Determining boundary conditions (BCs) for incompressible flows is such a delicate matter that affects the accuracy of the results. In this research, a new characteristic‐based BC for incompressible Navier‐Stokes equations is introduced. Discretization of equations has been done via finite volume. Additionally, artificial compressibility correction has been employed to deal with equations. Ordinary extrapolation from inner cells of a domain was used as a traditional way to estimate pressure and velocities on solid wall and inlet/outlet boundaries. Here, this method was substituted by the newly proposed BCs based on the characteristics of artificial compressibility equations. To follow this purpose, a computer code has been developed to carry out series of numerical tests for a flow over a backward‐facing step and was applied to a wide range of Reynolds numbers and grid combinations. Calculation of convective and viscous fluxes was done using Jameson's averaging scheme. Employing the characteristic‐based method for determining BCs has shown an improved convergence rate and reduced calculation time comparing with those of traditional ones. Furthermore, with the reduction of domain and computational cells, a similar accuracy was achieved for the results in comparison with the ones obtained from the traditional extrapolation method, and these results were in good agreement with the ones in the literature.     

Enhancement of the accuracy of the finite volume particle method for the simulation of incompressible flows

A finite volume particle (FVP) method for simulation of incompressible flows that provides enhanced accuracy is proposed. In this enhanced FVP method, a dummy neighbor particle is introduced for each particle in the calculation and used for the discretization of the gradient model and Laplacian model. The error‐compensating term produced by introducing the dummy neighbor particle enables higher order terms to be calculated. The proposed gradient model and Laplacian model are applied in both pressure and pressure gradient calculations. This enhanced FVP scheme provides more accurate simulations of incompressible flows. Several 2‐dimensional numerical simulations are given to confirm its enhanced performance.     

A COMPUTATIONAL METHOD FOR THE NAVIER-STOKES EQUATIONS AT ALL SPEEDS

ＩｎｔｒｏｄｕｃｔｉｏｎＦｏｒｃｏｍｐｕｔａｔｉｏｎｏｆｔｈｅｖｉｓｃｏｕｓｃｏｍｐｒｅｓｓｉｂｌｅｆｌｏｗｓｔｈｅＮａｖｉｅｒ＿Ｓｔｏｋｅｓｅｑｕａｔｉｏｎｓａｒｅｇｅｎｅｒａｌｌｙｒｅｐｒｅｓｅｎｔｅｄｉｎｔｈｅｃｏｎｓｅｒｖａｔｉｏｎｌａｗｆｏｒｍａｓａｈｙｐｅｒｂｏｌｉｃｓｙｓｔｅｍ .Ｌａｃｋｉｎｇｔｈｅｍａｔｈｅｍａｔｉｃａｌｔｏｏｌｔｏａｎａｌｙｚｅｔｈｉｓｎｏｎｌｉｎｅａｒｓｙｓｔｅｍ ,ｔｈｅｎｕｍｅｒｉｃａｌｍｅｔｈｏｄｓｕｓｅｄｉｎｓｏｌｖｉｎｇｔｈｅｎｏｎｌｉｎｅａｒｈｙｐｅ…     

A lattice Boltzmann model for simulation of compressible flows

This paper presents a new model of lattice Boltzmann method for full compressible flows. On the basis of multi‐speed model, an extra potential energy distribution function is introduced to recover the full compressible Navier–Stokes equations with a flexible specific‐heat ratio and Prandtl number. The Chapman–Enskog expansion of the kinetic equations is performed, and the two‐dimension‐seventeen‐velocity density equilibrium distribution functions are obtained. The governing equations are discretized using the third order monotone upwind scheme for scalar conservation laws finite volume scheme. The van Albada limiter is used to avoid spurious oscillations. In order to verify the accuracy of this double‐distribution‐function model, the Riemann problems, Couette flows, and flows around a NACA0012 airfoil are simulated. It is found that the proposed lattice Boltzmann model is suitable for compressible flows, even for strong shock wave problem, which has an extremely large pressure ratio, 100,000. Copyright © 2014 John Wiley & Sons, Ltd.     

An immersed boundary method for incompressible flows using volume of body function

A simple and effective immersed boundary method using volume of body (VOB) function is implemented on unstructured Cartesian meshes. The flow solver is a second‐order accurate implicit pressure‐correction method for the incompressible Navier–Stokes equations. The domain inside the immersed body is viewed as being occupied by the same fluid as outside with a prescribed divergence‐free velocity field. Under this view a fluid–body interface is similar to a fluid–fluid interface encountered in the volume of fluid (VOF) method for the two‐fluid flow problems. The body can thus be identified by the VOB function similar to the VOF function. In fluid–body interface cells the velocity is obtained by a volume‐averaged mixture of body and fluid velocities. The pressure inside the immersed body satisfies the same pressure Poisson equation as outside. To enhance stability and convergence, multigrid methods are developed to solve the difference equations for both pressure and velocity. Various steady and unsteady flows with stationary and moving bodies are computed to validate and to demonstrate the capability of the current method. Copyright © 2005 John Wiley & Sons, Ltd.     

Application of time preconditioning and high‐order compact discretization method for low Mach number flows

The paper describes a combination of a preconditioning method with a high‐order compact discretization scheme for the purpose of solving the compressible Navier–Stokes equations in moderate and low Mach number regimes. When combined with properly modified characteristic boundary conditions, the proposed approach is very efficient from the point of view of convergence acceleration and accuracy of the results. The computations were performed in typical benchmark cases including the Burggraf flow for which an analytical solution exists, the flow over a backward facing step, and also the flow in 2D and 3D shear‐driven cavities. Depending on the test case, the results were obtained for the Mach number in the range M = 0.001 ? 0.5 and the Reynolds number Re = 1 ? 1000. Copyright © 2012 John Wiley & Sons, Ltd.     

A Mach‐uniform preconditioner for incompressible and subsonic flows

In this study, a novel Mach‐uniform preconditioning method is developed for the solution of Euler equations at low subsonic and incompressible flow conditions. In contrast to the methods developed earlier in which the conservation of mass equation is preconditioned, in the present method, the conservation of energy equation is preconditioned, which enforces the divergence free constraint on the velocity field even at the limiting case of incompressible, zero Mach number flows. Despite most preconditioners, the proposed Mach‐uniform preconditioning method does not have a singularity point at zero Mach number. The preconditioned system of equations preserves the strong conservation form of Euler equations for compressible flows and recovers the artificial compressibility equations in the case of zero Mach number. A two‐dimensional Euler solver is developed for validation and performance evaluation of the present formulation for a wide range of Mach number flows. The validation cases studied show the convergence acceleration, stability, and accuracy of the present Mach‐uniform preconditioner in comparison to the non‐preconditioned compressible flow solutions. The convergence acceleration obtained with the present formulation is similar to those of the well‐known preconditioned system of equations for low subsonic flows and to those of the artificial compressibility method for incompressible flows. Copyright © 2013 John Wiley & Sons, Ltd.     

Local pressure preconditioning method for steady incompressible flows

The convergence and accuracy characteristics of the preconditioned incompressible Euler and Navier–Stokes equations are studied. An object-oriented C++ numerical code has been developed for solving the inviscid and viscous, steady, incompressible flows problems. The code is based on the cell-centred finite volume method. In this scheme, two-dimensional incompressible Euler and Navier–Stokes equations are modified by a robust artificial compressibility (AC) and a local preconditioning matrix of pressure-sensor type. The preconditioned equations are solved with the Jameson's numerical approach, i.e. artificial dissipation and artificial viscosity terms under the form of a fourth- and second-order derivative, respectively. An explicit four-stage Runge–Kutta integration algorithm is applied to obtain the steady-state condition. The computed results include the steady-state solution of flow past the NACA-hydrofoils and a circular cylinder in free stream, for which the numerical results are compared with numerical works of other researchers. Good agreement is observed. The effects of AC parameter, artificial viscosity and dissipation factor, and local preconditioning coefficient on convergence rate and solution accuracy are tested by computing flow over the NACA0012 hydrofoil. In addition, some important design criteria of a preconditioner, such as stiffness reduction, hyperbolicity, symmetrisability, accuracy preservation for M → 0, and M-property have been examined analytically.