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 numerical method for flows in porous and homogenous fluid domains coupled at the interface by stress jump

A numerical method was developed for flows involving an interface between a homogenous fluid and a porous medium. The numerical method is based on the finite volume method with body‐fitted and multi‐block grids. A generalized model, which includes Brinkman term, Forcheimmer term and non‐linear convective term, was used to govern the flow in the porous medium region. At its interface, a shear stress jump that includes the inertial effect was imposed, together with a continuity of normal stress. Furthermore, the effect of the jump condition on the diffusive flux was considered, additional to that on the convective part which has been usually considered. Numerical results of three flow configurations are presented. The method is suitable for coupled problems with regions of homogeneous fluid and porous medium, which have complex geometries. Copyright © 2006 John Wiley & Sons, Ltd.     

研究了雷诺数Re=200, 1000, 线速度比$\alpha =0.5$,

研究了雷诺数Re=200, 1000, 线速度比$\alpha =0.5$, 2.0, 4.0, 强迫振荡频率$f_{s}=0.1\sim 2.0$情况下的旋转振荡圆柱绕流问题. 通 过基于非结构同位网格有限体积法对Navier-Stokes方程进行数值求解. 对流项、扩 散项和非恒定项的离散格式均具有二阶精度，利用SIMPLE算法处理压力-速度耦合. 计算得到了作用力系数随不同控制参数的变化规律. 通过对升力系数的频谱分析得到 自然脱落频率和强迫振荡频率下的作用力振幅. 通过对不同频率作用力幅值的分析， 得到频率之间的竞争关系，进而定量地给出了不同尾迹涡脱落模式的分区图.     

Taylor-Galerkin-based spectral element methods for convection-diffusion problems

Several explicit Taylor-Galerkin-based time integration schemes are proposed for the solution of both linear and non-linear convection problems with divergence-free velocity. These schemes are based on second-order Taylor series of the time derivative. The spatial discretization is performed by a high-order Galerkin spectral element method. For convection-diffusion problems an operator-splitting technique is given that decouples the treatment of the convective and diffusive terms. Both problems are then solved using a suitable time scheme. The Taylor-Galerkin methods and the operator-splitting scheme are tested numerically for both convection and convection-diffusion problems.     

An approach for choosing discretization schemes and grid size based on the convection-diffusion equation

A new approach for selecting proper discretization schemes and grid size is presented. This method is based on the convection-diffusion equation and can provide insight for the Navier-Stokes equation. The approach mainly addresses two aspects, i.e., the practical accuracy of diffusion term discretization and the behavior of high wavenumber disturbances. Two criteria are included in this approach. First, numerical diffusion should not affect the theoretical diffusion accuracy near the length scales of interest. This is achieved by requiring numerical diffusion to be smaller than the diffusion discretization error. Second, high wavenumber modes that are much smaller than the length scales of interest should not be amplified. These two criteria provide a range of suitable scheme combinations for convective flux and diffusive flux and an ideal interval for grid spacing. The effects of time discretization on these criteria are briefly discussed.     

基于制造解的非结构二阶有限体积离散格式的精度测试与验证

随着计算机技术的飞速进步,计算流体力学得到迅猛发展,数值计算虽能够快速得到离散结果,但是数值结果的正确性与精度则需要通过严谨的方法来进行验证和确认.制造解方法和网格收敛性研究作为验证与确认的重要手段已经广泛应用于计算流体力学代码验证、精度分析、边界条件验证等方面.本文在实现标量制造解和分量制造解方法的基础上,通过将制造解方法精度测试结果与经典精确解(二维无黏等熵涡)精度测试结果进行对比,进一步证实了制造解精度测试方法的有效性,并将两种制造解方法应用于非结构网格二阶精度有限体积离散格式的精度测试与验证,对各种常用的梯度重构方法、对流通量格式、扩散通量格式进行了网格收敛性精度测试.结果显示,基于Green-Gauss公式的梯度重构方法在不规则网格上会出现精度降阶的情况,导致流动模拟精度严重下降,而基于最小二乘(least squares)的梯度重构方法对网格是否规则并不敏感.对流通量格式的精度测试显示,所测试的各种对流通量格式均能达到二阶精度,且各方法精度几乎相同;而扩散通量离散中界面梯度求解方法的选择对流动模拟精度有显著影响.     

Multidimensional FEM‐FCT schemes for arbitrary time stepping

The flux‐corrected‐transport paradigm is generalized to finite‐element schemes based on arbitrary time stepping. A conservative flux decomposition procedure is proposed for both convective and diffusive terms. Mathematical properties of positivity‐preserving schemes are reviewed. A nonoscillatory low‐order method is constructed by elimination of negative off‐diagonal entries of the discrete transport operator. The linearization of source terms and extension to hyperbolic systems are discussed. Zalesak's multidimensional limiter is employed to switch between linear discretizations of high and low order. A rigorous proof of positivity is provided. The treatment of non‐linearities and iterative solution of linear systems are addressed. The performance of the new algorithm is illustrated by numerical examples for the shock tube problem in one dimension and scalar transport equations in two dimensions. Copyright © 2003 John Wiley & Sons, Ltd.     

An ‘assumed deviatoric stress-pressure-velocity’ mixed finite element method for unsteady,convective, incompressible viscous flow: Part II: Computational studies

In Part I of this paper we presented a mixed finite element method, for solving unsteady, incompressible, convective flows, based on assumed ‘deviatoric stress–velocity–pressure’ fields in each element, which have the features: (i) the convective term is treated by the usual Galerkin technique; (ii) the unknowns in the global system of finite element equations are the nodal velocities, and the ‘constant term’ in the arbitrary pressure field over each element; and (iii) exact integrations are performed over each element. In this paper we present numerical studies, both for steady as well as unsteady cases, of the problems: (a) the driven cavity, (b) Jeffry–Hamel flow in a channel, (c) flow over a ‘backward’ or ‘downstream’ facing step, and (d) flow over a square step. All these problems are two-dimensional in nature, although certain 3-D solutions are to be presented in a separate paper. The present results are compared with those which are available in the literature and are based on alternative approaches to treat incompressibility and convective acceleration. The possible merits of the present method are thus pointed out.     

Ransom test results from various two-fluid schemes: Is enforcing hyperbolicity a thermodynamically consistent option?

The basic elliptic ill-posedness of physical models and numerical schemes for two-fluid flows is a recurring issue that has motivated the introduction of numerous possible correction strategies. In practical applications physical terms are generally present and regularize the models (viscosity, drag, surface tension, etc.). Yet, many numerical schemes were developed with the stringent and self-imposed constraint that the convective part of the models to be solved had to be hyperbolic, regardless of the type and magnitude of the particular physical regularizing terms. This leads to consider the simplest possible two-fluid “backbone” models corrected with the simplest “universal” terms to ensure hyperbolicity.Among the proposed corrections is the introduction of an interfacial pressure, either closed by algebraic relations or by supplementary evolution equations. Concurrently with the shift to hyperbolic behavior, these techniques also affect other features of systems: Kelvin–Helmholtz type instabilities are notably quenched at all scales, a highly undesirable effect in many practical situations. Less commonly recognized are also distortions in the transfers between kinetic, reversible, and irreversible energies, sometimes up to thermodynamic inconsistency.The present work aims at comparing on the standard Ransom-faucet test the results from various available hyperbolic and elliptic schemes and models against an explicit double Lagrange-plus-remap discretization of the basic elliptic, one-pressure, compressible, six-equations system (i.e. with energy equations). Four features are examined on this test: the entropy preservation, the stretched stream profile, the volume fraction discontinuity, and the unstable character of the analytical solution for the simplest backbone model.The paper highlights the fact that the convective part of two-fluid models might not be necessarily hyperbolic provided that it is physically consistent and numerically robust. Observation of published results for Ransom’s test shows that by enforcing hyperbolicity regardless of thermodynamical consistency, numerical models remove instabilities at the volume fraction discontinuity, but at the expense of distorted profiles of the stretched stream due to excessive numerical diffusion and to spurious forces in the momentum equation. The present approach provides a form of neutral starting point before including dissipative terms: robust but not excessively diffusive, with accurate capture of the stretched stream and volume fraction discontinuity for any practical mesh refinement. Moreover, and consistently with the chosen elliptic model, this numerical scheme eventually generates the elliptic instabilities for late times or fine meshes (but remains robust under the appropriate time step restrictions). It can be supplemented by any kind of small-scale regularization term in order to introduce a cut-off under which physical or numerical stability may be necessary.     

Upwind discretization of the steady Navier–Stokes equations

A discretization method is presented for the full, steady, compressible Navier–Stokes equations. The method makes use of quadrilateral finite volumes and consists of an upwind discretization of the convective part and a central discretization of the diffusive part. In the present paper the emphasis lies on the discretization of the convective part. The solution method applied solves the steady equations directly by means of a non-linear relaxation method accelerated by multigrid. The solution method requires the discretization to be continuously differentiable. For two upwind schemes which satisfy this requirement (Osher's and van Leer's scheme), results of a quantitative error analysis are presented. Osher's scheme appears to be increasingly more accurate than van Leer's scheme with increasing Reynolds number. A suitable higher-order accurate discretization of the convection terms is derived. On the basis of this higher-order scheme, to preserve monotonicity, a new limiter is constructed. Numerical results are presented for a subsonic flat plate flow and a supersonic flat plate flow with oblique shock wave–boundary layer interaction. The results obtained agree with the predictions made. Useful properties of the discretization method are that it allows an easy check of false diffusion and that it needs no tuning of parameters.     

Stability for Rayleigh–Benard Convective Solutions of the Boltzmann Equation

We consider the Boltzmann equation for a gas in a horizontal slab, subject to a gravitational force. The boundary conditions are of diffusive type, specifying the wall temperatures, so that the top temperature is lower than the bottom one (Benard setup). We consider a 2-dimensional convective stationary solution, which for small Knudsen numbers is close to the convective stationary solution of the Oberbeck–Boussinesq equations, near and above the bifurcation point, and prove its stability under 2-d small perturbations, for Rayleigh numbers above and close to the bifurcation point and for small Knudsen numbers.     

Simulation of multiple shock–shock interference using implicit anti‐diffusive WENO schemes

Accurate computations of two‐dimensional turbulent hypersonic shock–shock interactions that arise when single and dual shocks impinge on the bow shock in front of a cylinder are presented. The simulation methods used are a class of lower–upper symmetric‐Gauss–Seidel implicit anti‐diffusive weighted essentially non‐oscillatory (WENO) schemes for solving the compressible Navier–Stokes equations with Spalart–Allmaras one‐equation turbulence model. A numerical flux of WENO scheme with anti‐diffusive flux correction is adopted, which consists of first‐order and high‐order fluxes and allows for a more flexible choice of first‐order dissipative methods. Experimental flow fields of type IV shock–shock interactions with single and dual incident shocks by Wieting are computed. By using the WENO scheme with anti‐diffusive flux corrections, the present solution indicates that good accuracy is maintained and contact discontinuities are sharpened markedly as compared with the original WENO schemes on the same meshes. Computed surface pressure distribution and heat transfer rate are also compared with experimental data and other computational results and good agreement is found. Copyright © 2009 John Wiley & Sons, Ltd.     

Mathematical derivation of a finite volume formulation for laminar flow in complex geometries

This paper treats the mathematical derivation of a novel formulation of the Navier–Stokes equation for general non-orthogonal curvilinear co-ordinates. The covariant velocity components are solved in this FVM formulation, which leads to the pressure-velocity coupling becoming relatively easy to handle at the expense of a more complicated expression of the convective and diffusive fluxes. When a velocity component is solved at a point P, the neighbouring velocities are projected in the direction of the velocity component at the point P. Thus the base vectors are changed at the neighbouring points. This renders a simpler expression for the covariant derivatives. Neither the Cristoffel symbol nor its derivatives need be computed. This contributes to the accuracy of the formulation. The procedure of changing the base vectors affects only the convected velocity. The convecting term (dot product of velocity and area) is calculated without any change of the base vectors. The same is true for the operator on the covariant velocity in the diffusion term. It is shown that when using upwind differencing the use of projected velocities gives better results than when curvature effects are included in the source term. The discretized equations are written in a form which enables the use of the tridiagonal matrix algorithm (TDMA). The equations can be solved using either the SIMPLEC or the PISO procedure. Two examples of laminar flows are given.     

NUMERICAL SIMULATION FOR NAVIER-STOKES EQUATIONS BY PROJECTION/CHARACTERISTIC-BASED OPERATOR-SPLITTING FINITE ELEMENT METHOD

提出了一种求解非定常不可压缩纳维-斯托克斯方程（N-S方程）的新型有限元法：基于投影法的特征线算子分裂有限元法.在每一个时间层上将N-S方程分裂成扩散项、对流项、压力修正项.对流项采用多步显式格式,且在每一个对流子时间步内采用更加精确的显式特征线-伽辽金法进行时间离散,空间离散采用标准伽辽金法.应用此算法对平面泊肃叶流、方腔流和圆柱绕流进行数值模拟,所得结果与基准解符合良好.尤其对于Re=10000的方腔流,给出了方腔中分离涡发展和运动的计算结果,并发现在该雷诺数下存在周期解,表明该算法能较好地模拟流体流动中的小尺度物理量以及流场中分离涡的运动.     

N-S方程基于投影法的特征线算子分裂有限元求解

提出了一种求解非定常不可压缩纳维-斯托克斯方程（N-S方程）的新型有限元法:基于投影法的特征线算子分裂有限元法.在每一个时间层上将N-S方程分裂成扩散项、对流项、压力修正项.对流项采用多步显式格式,且在每一个对流子时间步内采用更加精确的显式特征线-伽辽金法进行时间离散,空间离散采用标准伽辽金法.应用此算法对平面泊肃叶流、方腔流和圆柱绕流进行数值模拟,所得结果与基准解符合良好.尤其对于Re=10000的方腔流,给出了方腔中分离涡发展和运动的计算结果,并发现在该雷诺数下存在周期解,表明该算法能较好地模拟流体流动中的小尺度物理量以及流场中分离涡的运动.      

非定常不可压N-S方程的最小二乘算子分裂有限元法数值求解

采用最小二乘算子分裂有限元法求解非定常不可压N-S(Navier-Stokes)方程,即在每个时间层上采用算子分裂法将N-S方程分裂成扩散项和对流项,这样既能考虑对流占优特点又能顾及方程的扩散性质。扩散项是一个抛物型方程,时间离散采用向后差分格式,空间离散采用标准Galerkin有限元法。对流项的时间项采用后向差分格式,非线性部分用牛顿法进行线性化处理,再用最小二乘有限元法进行空间离散,得到对称正定的代数方程组系数矩阵。采用Re=1000的方腔流对该算法的有效性进行检验,表明其具有较高的精度,能够很好地捕捉流场中的涡结构。同时,对圆柱层流绕流进行了数值研究,通过流线图、压力场、阻力系数、升力系数及斯特劳哈数等结果的分析与对比,表明本文算法对于模拟圆柱层流绕流是准确和可靠的。     

AN IMPROVED MAC METHOD (SIMAC) FOR UNSTEADY HIGH-REYNOLDS FREE SURFACE FLOWS

A modified MAC method (SIMAC; semi-implicit marker and cell) is proposed which accurately treats unsteady high-Reynolds free surface problems. SIMAC solves the Navier–Stokes equations in primitive variables on a non-uniform staggered Cartesian grid by means of a finite difference scheme. The convective term is treated explicitly by employing a second-order upwind scheme in space (HLPA) and the Adams–Bashforth technique in time. The diffusive part is solved by means of the implicit approximate factorization technique. A multigrid technique based on the additive correction strategy is employed to solve the Poisson equation for the pressure. Finally, the free surface treatment is carried out using massless particles which divide the domain of integration into full and empty cells as in a standard MAC method. The algorithm is used for the analysis of large-amplitude water sloshing in rectangular unbaffled and baffled containers. Experimental tests have been carried out in order to validate the algorithm. Numerical results satisfactorily agree with experimental data for the whole range of filing conditions analysed here. © 1997 by John Wiley & Sons, Ltd.     

Multidimensional upwind schemes for the shallow water equations

A multidimensional discretisation of the shallow water equations governing unsteady free-surface flow is proposed. The method, based on a residual distribution discretisation, relies on a characteristic eigenvector decomposition of each cell residual, and the use of appropriate distribution schemes. For uncoupled equations, multidimensional convection schemes on compact stencils are used, while for coupled equations, either system distribution schemes such as the Lax–Wendroff scheme or scalar schemes may be used. For steady subcritical flows, the equations can be partially diagonalised into a purely convective equation of hyperbolic nature, and a set of coupled equations of elliptic nature. The multidimensional discretisation, which is second-order-accurate at steady state, is shown to be superior to the standard Lax–Wendroff discretisation. For steady supercritical flows, the equations can be fully diagonalised into a set of convective equations corresponding to the steady state characteristics. Discontinuities such as hydraulic jumps, are captured in a sharp and non-oscillatory way. For unsteady flows, the characteristic equations remain coupled. An appropriate treatment of the coupling terms allows the discretisation of these equations at the scalar level. Although presently only first-order-accurate in space and time, the classical dam-break problem demonstrates the validity of the approach. © 1998 John Wiley & Sons, Ltd.     

同位网格摄动有限体积格式求解浮力驱动方腔流

利用对流扩散方程的摄动有限体积格式,在Rayleigh数从10$^{3}$ 到10$^{8}$的范围内对浮力驱动方腔流动问题作了数值模拟. 对流扩散方程的摄动 有限体积格式具有一阶迎风格式的简洁形式,使用相同的基点,重构近似精度高,特别是两 相邻控制体中心到公共界面的距离相等或不相等,PFV格式公式相同等优点. 在数值模拟中, 无论均匀网格还是非均匀网格均获得与DSC方法、自适应有限元法、多重网格法等Benchmark 解相符较好的数值结果,证明UPFV格式对高Rayleigh数对流传热问题的适用性和有效性.     

Surface tension driven flows in micro-gravity conditions

The importance of convective flows generated by surface tension gradients, in comparison with the ones generated by other driving forces, has been investigated in connection with space technological applications involving fluid processes. A theoretical model of the boundary conditions at the interface, considered free and diffusive, has been derived in general tensor form to allow for the use of non orthogonal curvilinear co-ordinates. For the study of flow fields contained in enclosures, these co-ordinates are more suitable to fit all teh boundaries, in particular near the contact angle between the interface and the solid walls, thus giving more accurate numerical solutions. A computational procedure to solve the complete set of bulk and surface equations is proposed and applied to a simplified two dimensional flow in a rectangular enclosure with a temperature gradient between the lateral walls. The numerical results show the importance of considering the interface to be deformable and diffusive for an accurate evaluation of the convective flow in the fluid bulk.