Poiseuille—Benard流的出口边界条件及其谱元法计算

研究二维矩形管道中底部加热的不可压缩Ｐｏｉｓｅｕｉｌｌｅ－Ｂｅｎａｒｄ流的谱元法数值计算问题。讨论各种不同的出口边界条件的处理及其对谱地数值模拟的影响,通过干扰区、干扰幅度和计算时间的比较,确定比较理想的出口边界条件。     

交错网格上应用有限谱QUICK法计算方腔流动

王健平教授提出的有限谱方法是一种局域化谱方法，具有精度高、无相位差、应用灵活等特点，在以往的实践中取得了很大成功。本文在交错网格上对二维驱动方腔流问题进行计算，求解了二维不可压缩流动的涡流流函数方程。其中微分部分采用有限谱法进行处理，对流项的处理则应用了QUICK格式。本文计算了雷诺数为1000、5000、10000、20000等多种情况，将所得的结果进行分析，并将中线上的速度分别同已有的文献数据进行对比，从而，验证有限谱微分的正确性和其在实际应用中的可行性。     

Poiseuille-Bénard流的出口边界条件及其谱元法计算

研究二维矩形管道中底部加热的不可压缩Ｐｏｉｓｅｕｉｌｌｅ－Ｂｅｎａｒｄ流的谱元法数值计算问题．讨论各种不同的出口边界条件的处理及其对谱元法数值模拟的影响．通过干扰区、干扰幅度和计算时间的比较,确定比较理想的出口边界条件．     

A THEORETICAL ANALYSIS OF THE FLOW DISTURBANCE INDUCED BY A SPHERE OSCILLATING IN INCOMPRESSIBLE FLUID

为了改进基于不可压缩流场的声类比法的气动声数值预测方法,首先要明确扰动在可压缩和不可压缩流体媒介中的传播特性. 推导了震荡小球在不可压缩流体中产生的小扰动的理论解,分析其速度场与压力场的特点,并与可压缩情况的解进行比较. 结果显示,速度场中包含传播速度为无穷大和有限值的分量;而压力场只有传播速度为无穷大的分量. 当流体黏性趋于零或小球震荡频率趋于无穷大时,其流场与经典声学中震荡小球声辐射问题的近场声一致,这表明震荡小球产生的近场扰动为不可压缩流场,即伪声.     

配置点谱方法-人工压缩法(SCM-ACM)求解不可压缩流体流动

开发了配置点谱方法SCM（spectral collocation method）与人工压缩法ACM（artificial compressibility method）相结合的方法SCM-ACM,用于求解不可压缩粘性流动问题。选取典型的方腔顶盖驱动流为研究测试对象,首先建立人工压缩格式的控制方程,其次采用SCM离散控制方程的空间偏微分项,推导出矩阵形式的代数方程,最后测试了SCM-ACM代码的有效性。结果显示,SCM-ACM能够有效求解不可压缩流动问题,并继承了谱方法的指数收敛特性,且具有ACM求解过程简单及易于实施的特点。     

配置点谱方法-人工压缩法(SCM-ACM)求解同心圆筒内流体流动

发展了配置点谱方法SCM(Spectral collocation method)和人工压缩法ACM(Artificial compressibility method)相结合的SCM-ACM数值方法,计算了柱坐标系下稳态不可压缩流动N-S方程组。选取典型的同心圆筒间旋转流动Taylor-Couette流作为测试对象,首先,采用人工压缩法获得人工压缩格式的非稳态可压缩流动控制方程;再将控制方程中的空间偏微分项用配置点谱方法进行离散,得到矩阵形式的代数方程;编写了SCM-ACM求解不可压缩流动问题的程序;最后,通过与公开发表的Taylor-Couette流的计算结果对比,验证了求解程序的有效性。结果证明,本文发展的SCM-ACM数值方法能够用于求解圆筒内不可压缩流体流动问题,该方法既保留了谱方法指数收敛的特性,也具有ACM形式简单和易于实施的特点。本文发展的SCM-ACM数值方法为求解柱坐标下不可压缩流体流动问题提供了一种新的选择。     

边界元法分析粘性液体小幅晃动

本文采用了边界元法对容器中粘性、不可压缩液体小幅晃动进行数值分析。在频域内考虑二维线性化Navier-Stokes方程,以问题的物理变量作为数值分析的未知函数,并推导了该问题分析的边界积分方程。自由面上的动力学条件为法向正应力和切向剪应力为零,这两个条件本身是线性的,避免了采用无粘势理论边界条件的非线性,固壁面上采用流体质点与固壁质点速度相等的条件,在数值计算过程中,结合有限差分法对边界条件进行了处理,由此建立了问题的一个边界元数值求解过程。     

粘性不可压缩流的壁涡数值边界条件

本文对具有固壁边界的二维粘性不可压缩定常流提出了一种四阶壁涡公式,并构造了一类非平凡的精确解。用两种具有代表性的内点差分格式及六种常用的壁涡公式(包括本文提出的一种),对本文所构造的雷诺数直到2000的一系列流动进行了数值求解,并以所构造的精确解为基准,分析和比较了这六种壁涡格式的精度、稳定性和收敛速度。计算结果表明,本文所提出的壁涡公式是这六种公式中最好的一种。本文数值实验中的一些结果对粘性不可压缩流固壁数值边界条件的处理可能具有一定的普遍意义。     

水/气多介质问题的界面处理方法

针对不可压缩可压缩水/气多介质问题, 提出一种新的界面处理方法。在可压缩水/气界面处构造Riemann问题, 在水中设音速趋于无穷大, 求解Riemann问题得到不可压缩可压缩水/气界面处流体的准确流动状态; 然后以此状态结合GFM(ghost fluid method)方法分别为2种流体定义界面边界条件, 将两相流问题转化为单相流问题计算, 通过求解level set方程来跟踪界面的位置。对各种不同的界面边界条件定义方法进行了比较, 数值模拟结果表明算法能准确地捕捉各类间断的位置, 证明了算法的有效性和稳健性。     

不可压缩粘性流动的CBS有限元解法

对于二维不可压缩粘性流动,首先通过坐标变换的方式得到了的不含对流项的NS方程,并给出了CBS有限元方法求解的一般过程。结合一类同时含有压力和速度的出口边界条件,对方腔顶盖驱动流、后向台阶绕流和圆柱绕流进行了计算。所得结果与基准解符合良好,验证了CBS算法对于定常、非定常粘性不可压缩流动问题的可行性和所用出口边界条件的无反射特性。特别的,对于圆柱绕流,Re=100时非定常升、阻力系数及漩涡脱落等非定常都得到了较好地模拟,为一进步研究自激振动等更加复杂的非定常流动问题奠定了基础。     

Consistent inlet and outlet boundary conditions for particle methods

In this paper, simple and consistent open boundary conditions are presented for the numerical simulation of viscous incompressible laminar flows. The present approach is based on an arbitrary Lagrangian-Eulerian particle method using upwind interpolation. Three kinds of inlet/outlet boundary conditions are proposed for particle methods, a pressure specified inlet/outlet condition, a velocity profile specified inlet/outlet condition, and a fully developed flow outlet condition. These inlet/outlet conditions are realized by using boundary particles and modification to the physical value such as velocity. Poiseuille flows, flows over a backward-facing step, and flows in a T-shape branch are calculated. The results are compared with those of mesh-based methods such as the finite volume method. The method presented herein exhibits accuracy and numerical stability.     

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.     

A promising boundary element formulation for three‐dimensional viscous flow

In this paper, a new set of boundary‐domain integral equations is derived from the continuity and momentum equations for three‐dimensional viscous flows. The primary variables involved in these integral equations are velocity, traction, and pressure. The final system of equations entering the iteration procedure only involves velocities and tractions as unknowns. In the use of the continuity equation, a complex‐variable technique is used to compute the divergence of velocity for internal points, while the traction‐recovery method is adopted for boundary points. Although the derived equations are valid for steady, unsteady, compressible, and incompressible problems, the numerical implementation is only focused on steady incompressible flows. Two commonly cited numerical examples and one practical pipe flow problem are presented to validate the derived equations. Copyright © 2004 John Wiley & Sons, Ltd.     

不可压缩Stokes流动的PSPG无网格法

将应用于有限元法的Pressure-Stabilizing/Petrov-Galerkin(PSPG)稳定化机制引入到无网格法中,有效消除了由于速度和压力的插值模式违反LBB条件而导致的压力场的虚假振荡。采用与有限元法耦合的连续掺混法(Continuous Blending Method)施加本质边界条件,使得边界条件不仅在边界节点上而且在整条边界上都得到严格满足。给出了两个典型算例的数值模拟结果,表明了所建议无网格法模拟不可压缩Stokes流动的有效性。     

A pressure correction method for fluid–particle interaction flow: Direct‐forcing method and sedimentation flow

A direct‐forcing pressure correction method is developed to simulate fluid–particle interaction problems. In this paper, the sedimentation flow is investigated. This method uses a pressure correction method to solve incompressible flow fields. A direct‐forcing method is introduced to capture the particle motions. It is found that the direct‐forcing method can also be served as a wall‐boundary condition. By applying Gauss's divergence theorem, the formulas for computing the hydrodynamic force and torque acting on the particle from flows are derived from the volume integral of the particle instead of the particle surface. The order of accuracy of the present method is demonstrated by the errors of velocity, pressure, and wall stress. To demonstrate the efficiency and capability of the present method, sedimentations of many spherical particles in an enclosure are simulated. Copyright © 2010 John Wiley & Sons, Ltd.     

AN APPROXIMATE PROJECTION SCHEME FOR INCOMPRESSIBLE FLOW USING SPECTRAL ELEMENTS

An approximate projection scheme based on the pressure correction method is proposed to solve the Navier–Stokes equations for incompressible flow. The algorithm is applied to the continuous equations; however, there are no problems concerning the choice of boundary conditions of the pressure step. The resulting velocity and pressure are consistent with the original system. For the spatial discretization a high-order spectral element method is chosen. The high-order accuracy allows the use of a diagonal mass matrix, resulting in a very efficient algorithm. The properties of the scheme are extensively tested by means of an analytical test example. The scheme is further validated by simulating the laminar flow over a backward-facing step.     

The effect of Reynolds number on boundary layer turbulence

A new facility for studying high Reynolds number incompressible turbulent boundary layer flows has been constructed. It consists of a moderately sized wind tunnel, completely enclosed by a pressure vessel, which can raise the ambient air pressure in and around the wind tunnel to 8 atmospheres. This results in a Reynolds number range of about 20:1, while maintaining incompressible flow. Results are presented for the zero pressure gradient flat plate boundary layer over a momentum thickness Reynolds number range 1500–15?000. Scaling issues for high Reynolds number non-equilibrium boundary layers are discussed, with data comparing the three-dimensional turbulent boundary layer flow over a swept bump at Reynolds numbers of 3800 and 8600. It is found that successful prediction of these types of flows must include length scales which do not scale on Reynolds number, but are inherent to the geometry of the flow.     

Computation of unsteady viscous incompressible flows in generalized non‐inertial co‐ordinate system using Godunov‐projection method and overlapping meshes

Time‐dependent incompressible Navier–Stokes equations are formulated in generalized non‐inertial co‐ordinate system and numerically solved by using a modified second‐order Godunov‐projection method on a system of overlapped body‐fitted structured grids. The projection method uses a second‐order fractional step scheme in which the momentum equation is solved to obtain the intermediate velocity field which is then projected on to the space of divergence‐free vector fields. The second‐order Godunov method is applied for numerically approximating the non‐linear convection terms in order to provide a robust discretization for simulating flows at high Reynolds number. In order to obtain the pressure field, the pressure Poisson equation is solved. Overlapping grids are used to discretize the flow domain so that the moving‐boundary problem can be solved economically. Numerical results are then presented to demonstrate the performance of this projection method for a variety of unsteady two‐ and three‐dimensional flow problems formulated in the non‐inertial co‐ordinate systems. Copyright © 2002 John Wiley & Sons, Ltd.     

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.     

Gauge finite element method for incompressible flows

A finite element method for computing viscous incompressible flows based on the gauge formulation introduced in [Weinan E, Liu J‐G. Gauge method for viscous incompressible flows. Journal of Computational Physics (submitted)] is presented. This formulation replaces the pressure by a gauge variable. This new gauge variable is a numerical tool and differs from the standard gauge variable that arises from decomposing a compressible velocity field. It has the advantage that an additional boundary condition can be assigned to the gauge variable, thus eliminating the issue of a pressure boundary condition associated with the original primitive variable formulation. The computational task is then reduced to solving standard heat and Poisson equations, which are approximated by straightforward, piecewise linear (or higher‐order) finite elements. This method can achieve high‐order accuracy at a cost comparable with that of solving standard heat and Poisson equations. It is naturally adapted to complex geometry and it is much simpler than traditional finite element methods for incompressible flows. Several numerical examples on both structured and unstructured grids are presented. Copyright © 2000 John Wiley & Sons, Ltd.