一种基于双界面函数的界面捕捉方法

基于代数重构思想,发展了一种新的双界面函数重构方法,并采用双正弦函数构造了双正弦界面重构方法(double sine interface capturing,DSINC).为验证不同界面函数对界面捕捉效果的影响,用数值方法求解了可压缩五方程模型,其中对流项的离散采用五阶WENO(weighted essentially non-oscillatory method)格式,时间积分采用三阶Runge--Kutta方法,通量计算分别考虑了HLL和HLLC方法,而状态方程采用Mie-Gr¨uneisen状态方程.在数值计算中,在界面附近,采用DSINC来获得体积分数的重构,而在远离界面的区域采用WENO格式来获得高阶插值状态.相比采用单界面函数的方法,如双曲正切界面重构方法(tangent of hyperbola for interface capturing,THINC),DSINC方法同样具有界面重构算法简单,在程序中添加方便等特点,两者区别在于,DSINC方法在重构过程中未知函数更易于求解,而无需求解复杂的非线性超越方程,这就使其具有易于向多维扩展的能力.一些典型的两相流动问题,如圆形水柱对流问题,两相三波点问题和激波-界面不稳定性问题等被用作不同界面函数对界面捕捉效果的影响对比.对比分析发现,DSINC与THINC在界面捕捉效果上大致保持一致,并在计算中表现出了较好的稳定性.双界面函数重构思想可以为多相流动界面的代数重构提供了一种新的思路.     

A NEW INTERFACE CAPTURING METHOD BASED ON DOUBLE INTERFACE FUNCTIONS

基于代数重构思想,发展了一种新的双界面函数重构方法,并采用双正弦函数构造了双正弦界面重构方法（double sine interface capturing,DSINC）.为验证不同界面函数对界面捕捉效果的影响,用数值方法求解了可压缩五方程模型,其中对流项的离散采用五阶WENO（weighted essentially non-oscillatory method）格式,时间积分采用三阶Runge——Kutta方法,通量计算分别考虑了HLL和HLLC方法,而状态方程采用Mie-Grüneisen状态方程.在数值计算中,在界面附近,采用DSINC来获得体积分数的重构,而在远离界面的区域采用WENO格式来获得高阶插值状态.相比采用单界面函数的方法,如双曲正切界面重构方法（tangent of hyperbola for interfacecapturing,THINC）,DSINC方法同样具有界面重构算法简单,在程序中添加方便等特点,两者区别在于,DSINC方法在重构过程中未知函数更易于求解,而无需求解复杂的非线性超越方程,这就使其具有易于向多维扩展的能力.一些典型的两相流动问题,如圆形水柱对流问题,两相三波点问题和激波——界面不稳定性问题等被用作不同界面函数对界面捕捉效果的影响对比.对比分析发现,DSINC与THINC在界面捕捉效果上大致保持一致,并在计算中表现出了较好的稳定性.双界面函数重构思想可以为多相流动界面的代数重构提供了一种新的思路.     

二维溃坝绕流的WENO格式数值模拟

WENO(Weighted essentially non-oscillatory scheme)格式是一类新的高精度无振荡差分格式.本文将WENO格式和Runge-Kutta时间离散的思想应用于二维浅水方程组的求解,数值模拟矩形河道中大坝瞬间局部溃倒,下游有障碍物的洪水演进过程,并对模拟结果进行了分析,表明采用WENO格式所建立的高分辨率模型能够有效地模拟溃坝波的演进过程.     

Mie-Grneisen状态方程可压缩多流体流动的PPM方法

采用流体体积分数的混合型多流体数值模型，将piecewise parabolic method (PPM)方法应用于可压缩多流体流动的数值模拟，拓展了以前提出的模型和数值方法，使它能够处理一般的Mie-Grneisen状态方程。采用双波近似和两层迭代算法求解一般状态方程的Riemann问题；并根据多流体接触界面无振荡原则设计高精度计算格式，对典型的纯界面平移问题可以从理论上证明本算法在接触间断附近压力和速度没有振荡，而且数值模拟结果表明界面数值耗散也被控制在2~3个网格之内。模拟了多种复杂的可压缩多流体流动，算例结果表明本文方法可以有效地处理接触间断、激波等物理问题，且具有耗散小精度高的特点。     

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

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

非线性流体-刚体结构相互作用问题的一种数值模拟方法

给出了一种模拟非线性流体-刚体结构相互作用问题的数值方法．文中假定结构承受大的刚体运动,流体流动受非线性有粘或无粘的场方程支配并满足自由表面和两相耦合界面上的非线性边界条件,利用任意拉氏-欧氏（ALE）网格系统构造了数值模型．采用所探讨的多块数值格式,允许可动重造网格间有独立的相对运动,从而克服了流体网格与固体大运动匹配的困难．通过数值离散化,导出了描述非线性流固耦合动力学的数值方程并应用耦合迭代过程对其作了求解．通过算例,说明了所提出数值模型的应用．     

基于VPM-THINC/QQ模型的波浪高保真模拟

为实现波浪传播的高保真数值模拟,采用包含单元均值和点值(volume-average/point-value method,VPM)的有限体积法求解纳维-斯托克斯方程和具有二次曲面性质和高斯积分的双曲正切函数(THINC method with quadratic surface representation and Gaussian quadrature,THINC/QQ)方法来重构自由面,建立以开源求解库OpenFOAM底层函数库为基础的VPM-THINC/QQ模型. 在本模型中添加推板造波法实现波浪的产生功能,采用松弛法实现消波功能,构建高精度黏性流数值波浪水槽. 分别采用VPM-THINC/QQ模型和InterFoam求解器(OpenFOAM软件包中广泛使用的多相流求解器)开展规则波的数值模拟,重点探究网格大小和时间步长等因素对波浪传播过程的影响,定量地分析波高衰减程度;为验证本模型的适应性,对长短波进行模拟. 结果表明,在相同网格大小或时间步长条件下,VPM-THINC/QQ模型的预测结果与参考值吻合较好,波高衰减较少,且无相位差,在波浪传播过程的模拟中呈现出良好的保真性. 本文工作 为波浪传播的模拟研究提供了一种高精度的黏性数值波浪水槽模型.      

无粘可压缩流动的改进高精度方法

为更准确捕捉复杂流场的流动细节,通过对WENO格式的光滑因子进行改进,发展了一种新的五阶WENO格式。对三阶ENO格式进行加权可以得到五阶WENO格式,但是不同的加权处理,WENO格式在极值处保持加权基本无振荡的效果不同,本文构造了二阶精度的局部光滑因子,及不含一阶二阶导数的高阶全局光滑因子,从而实现WENO格式在极值处有五阶精度。基于改进五阶WENO格式,对一维对流方程、一维和二维可压缩无粘问题进行算例验证,并与传统WENO-JS格式和WENO-Z格式进行比较。计算结果表明,改进五阶WENO格式有较高的精度和收敛速度,有较低的数值耗散,能有效捕捉间断、激波和涡等复杂流动。     

基于流场定常化的桥梁颤振分析简化数值方法

基于数值方法实现快速有效地分析评价大跨桥梁颤振稳定性.针对强迫振动法识别颤振导数试验中的大跨桥梁二维节段模型.利用RNG κ-ε湍流模型并采用有限体积法求解桥梁模型绕流二维不可压缩流体Navier-Stokes方程.通过计算桥梁断面模型在周期运动中少数离散时刻的气动力,利用最小二乘法计算颤振导数,采用SCANLAN方法求解颤振临界风速,最终进行颤振稳定性分析.通过该方法计算出了丹麦Great Belt East桥和我国虎门大桥的颠振临界风速.计算结果与已有试验结果十分接近,进而验证了该数值方法的有效性和可靠性.     

基于VPM-THINC/QQ模型的波浪高保真模拟

为实现波浪传播的高保真数值模拟,采用包含单元均值和点值(volume-average/point-value method,VPM)的有限体积法求解纳维-斯托克斯方程和具有二次曲面性质和高斯积分的双曲正切函数(THINC method with quadratic surface representation and Gaussian quadrature,THINC/QQ)方法来重构自由面,建立以开源求解库OpenFOAM底层函数库为基础的VPM-THINC/QQ模型.在本模型中添加推板造波法实现波浪的产生功能,采用松弛法实现消波功能,构建高精度黏性流数值波浪水槽.分别采用VPM-THINC/QQ模型和InterFoam求解器(OpenFOAM软件包中广泛使用的多相流求解器)开展规则波的数值模拟,重点探究网格大小和时间步长等因素对波浪传播过程的影响,定量地分析波高衰减程度;为验证本模型的适应性,对长短波进行模拟.结果表明,在相同网格大小或时间步长条件下,VPM-THINC/QQ模型的预测结果与参考值吻合较好,波高衰减较少,且无相位差,在波浪传播过程的模拟中呈现出良好的保真性.本文工作为波浪传播的模拟研究提供了一种高精度的黏性数值波浪水槽模型.     

Numerical modeling of wave interactions with coastal structures by a constrained interpolation profile/immersed boundary method

A high‐order difference method based multiphase model is proposed to simulate nonlinear interactions between water wave and submerged coastal structures. The model is based on the Navier–Stokes equations using a constrained interpolation profile (CIP) method for the flow solver, and employs an immersed boundary method (IBM) for the treatment of wave–structure interactions. A more accurate interface capturing scheme, the volume of fluid/weighed line interface calculation (VOF/WLIC) scheme, is adopted as the interface capturing method. A series of computations are performed to verify the application of the model for simulations of fluid interaction with various structures. These problems include flow over a fixed cylinder, water entry of a circular cylinder and solitary waves passing various submerged coastal structures. Computations are compared with the available analytical, experimental and other numerical results and good agreement is obtained. The results of this study demonstrate the accuracy and applications of the proposed model to simulate the nonlinear flow phenomena and capture the complex free surface flow. Copyright © 2015 John Wiley & Sons, Ltd.     

用隐式方法和WENO格式计算悬停旋翼跨声速无粘流场

寻找一种能够准确计算以涡为主要特征的复杂流场和克服尾迹耗散问题的数值方法,一直是旋翼空气动力学研究的热点和难点。本文发展了一种基于高阶迎风格式计算悬停旋翼无粘流场的隐式数值方法。无粘通量采用Roe通量差分分裂格式,为提高精度,使用五阶WENO格式进行左右状态插值,并与MUSCL插值进行比较。为提高收敛到定常解的效率,时间推进采用LU-SGS隐式方法。用该方法对一跨声速悬停旋翼无粘流场进行了数值计算,数值结果表明WENO-Roe的激波分辨率高于MUSCL-Roe,体现出了格式精度的提高对计算结果的改善,LU-SGS隐式方法的计算效率比5步Runge-Kutta显式方法的高。     

Improved THINC/SW scheme for computing incompressible two-phase flows

In this paper, we propose an improved tangent of hyperbola for interface capturing with slope weighting (THINC/SW) scheme for computing incompressible two-phase flows with surface tension on fixed Eulerian grids. The new scheme possesses the following major new properties in comparison with the original THINC/SW scheme: (i) providing a simple and accurate approach for generating a smooth level set (LS) field from the discontinuous volume-of-fluid field, (ii) determining the interface slope from the cogenerated LS field, and (iii) evaluating the surface tension force by the cogenerated LS field. We verified the proposed scheme with the widely used advection benchmark tests and multifluid simulations. Numerical results reveal that the new scheme can not only improve the solution quality of interface capturing but also increase the accuracy of curvature estimates and suppress the parasitic currents.     

Rotor wake capture improvement based on high-order spatially accurate schemes and chimera grids

A high-order upwind scheme has been developed to capture the vortex wake of a helicopter rotor in the hover based on chimera grids. In this paper, an improved fifth-order weighted essentially non-oscillatory (WENO) scheme is adopted to interpolate the higher-order left and right states across a cell interface with the Roe Riemann solver updating inviscid flux, and is compared with the monotone upwind scheme for scalar conservation laws (MUSCL). For profitably capturing the wake and enforcing the period boundary condition, the computation regions of flows are discretized by using the structured chimera grids composed of a fine rotor grid and a cylindrical background grid. In the background grid, the mesh cells located in the wake regions are refined after the solution reaches the approximate convergence. Considering the interpolation characteristic of the WENO scheme, three layers of the hole boundary and the interpolation boundary are searched. The performance of the schemes is investigated in a transonic flow and a subsonic flow around the hovering rotor. The results reveal that the present approach has great capabilities in capturing the vortex wake with high resolution, and the WENO scheme has much lower numerical dissipation in comparison with the MUSCL scheme.     

Numerical simulation of free surface flows with the level set method using an extremely high-order accuracy WENO advection scheme

Three-dimensional dynamic gas–liquid flow simulations that accurately track the phase interface are numerically challenging. This article presents a numerical study of the performance of the level-set phase–interface tracking method when combined with extremely high order (7th to 11th) weighted essentially non-oscillatory (WENO) advection schemes for gas–liquid free surface flows. Comparisons between simulation results and prior benchmark results suggest that such a combination of methods can be satisfactorily applied to the level-set and Navier-Stokes equations for free surface flow simulations when volume conservation is enforced at every time step, and minor numerical oscillations are suppressed through use of an artificial viscosity term. In particular, simulations of solid body rotation, the unsteady flow following an ideal dam break, tank sloshing, and the rise of a single bubble all agree with analytical or experimental results to within ± 3.12% when the level-set method is combined with an 11th order WENO scheme. Furthermore, use of an 11th order WENO advection scheme actually has a computational cost advantage because, for the same accuracy, it can be used on a coarser grid when compared with a more-common second-order advection scheme; computational savings of up to 87% are possible.     

Numerical experiments with several variant WENO schemes for the Euler equations

Numerical experiments with several variants of the original weighted essentially non‐oscillatory (WENO) schemes (J. Comput. Phys. 1996; 126 :202–228) including anti‐diffusive flux corrections, the mapped WENO scheme, and modified smoothness indicator are tested for the Euler equations. The TVD Runge–Kutta explicit time‐integrating scheme is adopted for unsteady flow computations and lower–upper symmetric‐Gauss–Seidel (LU‐SGS) implicit method is employed for the computation of steady‐state solutions. A numerical flux of the variant WENO scheme in flux limiter form is presented, which consists of first‐order and high‐order fluxes and allows for a more flexible choice of low‐order schemes. Computations of unsteady oblique shock wave diffraction over a wedge and steady transonic flows over NACA 0012 and RAE 2822 airfoils are presented to test and compare the methods. Various aspects of the variant WENO methods including contact discontinuity sharpening and steady‐state convergence rate are examined. By using the WENO scheme with anti‐diffusive flux corrections, the present solutions indicate that good convergence rate can be achieved and high‐order accuracy is maintained and contact discontinuities are sharpened markedly as compared with the original WENO schemes on the same meshes. Copyright © 2008 John Wiley & Sons, Ltd.     

A high‐resolution scheme for compressible multicomponent flows with shock waves

A simple methodology for a high‐resolution scheme to be applied to compressible multicomponent flows with shock waves is investigated. The method is intended for use with direct numerical simulation or large eddy simulation of compressible multicomponent flows. The method dynamically adds non‐linear artificial diffusivity locally in space to capture different types of discontinuities such as a shock wave, contact surface or material interface while a high‐order compact differencing scheme resolves a broad range of scales in flows. The method is successfully applied to several one‐dimensional and two‐dimensional compressible multicomponent flow problems with shock waves. The results are in good agreement with experiments and earlier computations qualitatively and quantitatively. The method captures unsteady shock and material discontinuities without significant spurious oscillations if initial start‐up errors are properly avoided. Comparisons between the present numerical scheme and high‐order weighted essentially non‐oscillatory (WENO) schemes illustrate the advantage of the present method for resolving a broad range of scales of turbulence while capturing shock waves and material interfaces. Also the present method is expected to require less computational cost than popular high‐order upwind‐biased schemes such as WENO schemes. The mass conservation for each species is satisfied due to the strong conservation form of governing equations employed in the method. Copyright © 2010 John Wiley & Sons, Ltd.     

加权型紧致格式与加权本质无波动格式的比较

线性紧致格式和加权本质无波动格式是两种典型的高阶精度数值格式,它们各有优缺点.线性紧致格式在具有高阶精度的同时,格式的分辨率也比较高,耗散低,是计算多尺度流场结构的较好格式,但是不能计算具有强激波的流场.加权本质无波动格式是一种高阶精度捕捉激波格式,鲁棒性好,但耗散比较高,分辨率也不理想.近年来,在莱勒的线性紧致格式基础上,采用加权本质无波动格式捕捉激波思想,发展了一系列加权型紧致格式.本文较全面地比较了加权型紧致格式和加权本质无波动格式,包括构造方法、鲁棒性、分辨率、耗散特性、收敛特性以及并行计算效率.结果表明,现有的加权型紧致格式基本保持了加权本质无波动格式的性质,对于气动力等宏观量的计算,比加权本质无波动格式没有明显的优势.      

An explicit finite volume element method for solving characteristic level set equation on triangular grids

Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow phenomena.This paper presents a characteristic level set equation which is derived from the two-dimensional level set equation by using the characteristic-based scheme.An explicit finite volume element method is developed to discretize the equation on triangular grids.Several examples are presented to demonstrate the performance of the proposed method for calculating interface evolutions in time.The proposed level set method is also coupled with the Navier-Stokes equations for two-phase immiscible incompressible flow analysis with surface tension.The Rayleigh-Taylor instability problem is used to test and evaluate the effectiveness of the proposed scheme.