Study of a staggered fourth‐order compact scheme for unsteady incompressible viscous flows

The objective of this paper is the development and assessment of a fourth‐order compact scheme for unsteady incompressible viscous flows. A brief review of the main developments of compact and high‐order schemes for incompressible flows is given. A numerical method is then presented for the simulation of unsteady incompressible flows based on fourth‐order compact discretization with physical boundary conditions implemented directly into the scheme. The equations are discretized on a staggered Cartesian non‐uniform grid and preserve a form of kinetic energy in the inviscid limit when a skew‐symmetric form of the convective terms is used. The accuracy and efficiency of the method are demonstrated in several inviscid and viscous flow problems. Results obtained with different combinations of second‐ and fourth‐order spatial discretizations and together with either the skew‐symmetric or divergence form of the convective term are compared. The performance of these schemes is further demonstrated by two challenging flow problems, linear instability in plane channel flow and a two‐dimensional dipole–wall interaction. Results show that the compact scheme is efficient and that the divergence and skew‐symmetric forms of the convective terms produce very similar results. In some but not all cases, a gain in accuracy and computational time is obtained with a high‐order discretization of only the convective and diffusive terms. Finally, the benefits of compact schemes with respect to second‐order schemes is discussed in the case of the fully developed turbulent channel flow. Copyright © 2008 John Wiley & Sons, Ltd.     

A diffuse-interface immersed boundary method for simulation of compressible viscous flows with stationary and moving boundaries

In this paper, a diffuse-interface immersed boundary method (IBM) is proposed for simulation of compressible viscous flows with stationary and moving boundaries. In the method, the solution of flow field and the implementation of boundary conditions are decoupled into two steps by applying the fractional step technique, ie, the predictor step and the corrector step. Firstly, in the predictor step, the intermediate flow field is resolved by a recently developed gas kinetic flux solver (GKFS) without consideration of the solid boundary. The GKFS is a finite volume approach that solves the Navier-Stokes equations for the flow variables at cell centers. In GKFS, the inviscid and viscous fluxes are evaluated as a single entity by reconstructing the local solution of continuous Boltzmann equation. Secondly, in the corrector step, the intermediate flow field is corrected by the present diffuse-interface IBM. During this process, the velocity field is firstly corrected by the implicit boundary condition–enforced IBM so that the no-slip boundary condition can be accurately satisfied. After that, the density correction is made by an iterative approach with the help of the continuity equation. Finally, the correction of the temperature field is made in the same way as that of the velocity field. Good agreements between the present simulations and the reference data in literature demonstrate the reliability of the proposed method.     

Analytical solution of a double moving boundary problem for nonlinear flows in one-dimensional semi-infinite long porous media with low permeability简

Based on Huang＇s accurate tri-sectional nonlin- ear kinematic equation （1997）, a dimensionless simplified mathematical model for nonlinear flow in one-dimensional semi-infinite long porous media with low permeability is presented for the case of a constant flow rate on the inner boundary. This model contains double moving boundaries, including an internal moving boundary and an external mov- ing boundary, which are different from the classical Stefan problem in heat conduction： The velocity of the external moving boundary is proportional to the second derivative of the unknown pressure function with respect to the distance parameter on this boundary. Through a similarity transfor- mation, the nonlinear partial differential equation （PDE） sys- tem is transformed into a linear PDE system. Then an ana- lytical solution is obtained for the dimensionless simplified mathematical model. This solution can be used for strictly checking the validity of numerical methods in solving such nonlinear mathematical models for flows in low-permeable porous media for petroleum engineering applications. Finally, through plotted comparison curves from the exact an- alytical solution, the sensitive effects of three characteristic parameters are discussed. It is concluded that with a decrease in the dimensionless critical pressure gradient, the sensi- tive effects of the dimensionless variable on the dimension- less pressure distribution and dimensionless pressure gradi- ent distribution become more serious; with an increase in the dimensionless pseudo threshold pressure gradient, the sensi- tive effects of the dimensionless variable become more serious; the dimensionless threshold pressure gradient （TPG） has a great effect on the external moving boundary but has little effect on the internal moving boundary.     

A projection method for solving incompressible viscous flows on domains with moving boundaries

In this paper, a projection method is presented for solving the flow problems in domains with moving boundaries. In order to track the movement of the domain boundaries, arbitrary‐Lagrangian–Eulerian (ALE) co‐ordinates are used. The unsteady incompressible Navier–Stokes equations on the ALE co‐ordinates are solved by using a projection method developed in this paper. This projection method is based on the Bell's Godunov‐projection method. However, substantial changes are made so that this algorithm is capable of solving the ALE form of incompressible Navier–Stokes equations. Multi‐block structured grids are used to discretize the flow domains. The grid velocity is not explicitly computed; instead the volume change is used to account for the effect of grid movement. A new method is also proposed to compute the freestream capturing metrics so that the geometric conservation law (GCL) can be satisfied exactly in this algorithm. This projection method is also parallelized so that the state of the art high performance computers can be used to match the computation cost associated with the moving grid calculations. Several test cases are solved to verify the performance of this moving‐grid projection method. Copyright © 2004 John Wiley Sons, Ltd.