多介质流体非守恒律欧拉方程组的数值计算方法

对多介质流体在界面处满足的Euler方程进行了探讨 ,方程组中增加了描述材料参数间断性质的对流形式非守恒律方程组。以波传播算法为基础 ,通过Roe方程近似求解Riemann问题 ,同时采用相同的数值差分格式求解流体动力学Euler方程组和界面方程组。该方法可以有效消除多介质流体在界面处压力、速度可能出现的非物理振荡。给出了部分典型一维和二维数值计算结果。     

A HIGH RESOLUTION CELL-CENTERED ALE METHOD FOR LARGE DEFORMATION MULTIMATERIAL INTERFACE

针对基于积分形式的Euler方程组耦合质量组份模型方程而发展的多介质整体ALE方法耗散大的问题,采用基于微分形式的Euler方程组所发展的高分辨率界面反耗散的思想来控制界面处的数值耗散,发展了一种二维平面中积分意义下的耦合质量组份方程和体积组份方程的界面反耗散的高分辨率多介质中心型ALE方法,从而高分辨率地模拟大变形物质界面。     

多介质可压缩流体动力学界面捕捉方法

研究多介质流界面捕捉方法的主要目的是消除多介质流体在界面处压力、速度可能出现的非物理振荡现象 ,并通过流体动力学方程和界面捕捉方程的耦合 ,将多介质流体动力学计算形式上转化为单介质流体计算 ,从而可以采用对计算单介质有效的高精度计算方法来处理多介质流动问题。推广了Shyue界面捕捉和其等效方程的推导方法 ,给出的结果可以适用于具有状态方程 p =( ,e,a1 , ,an) +( ,e,b1 , ,bn)e的介质 ,并通过了数值试验验证。     

空气自由场中强爆炸冲击波传播二维数值模拟

编写了适用于模拟具有高密度比、高压力比的强激波问题的二维柱对称多介质流体计算程序。利用有限体积方法求解流体的Euler方程组,采用level set方法捕捉爆炸产物与空气的运动界面,并通过求解物质界面两侧Riemann问题的精确解来计算爆炸产物与空气之间的数值通量。研制了三角形网格自适应技术来实现网格的自动加密和粗化,在保证捕捉激波峰值的前提下有效地提高了计算效率。利用计算程序对1 kt TNT当量的空气自由场强爆炸问题进行数值模拟,计算得到的峰值超压、冲击波到达时间等物理参数与点爆炸理论结果基本一致。

     

界面不稳定性的自适应欧拉数值计算

采用欧拉网格自适应算法数值模拟Richtmyer Meshkov和Rayleigh Taylor不稳定多介质流界面,获得了高精度界面特征。对不同流体引入不同位标函数跟踪界面运动,将位标函数方程与流体动力学方程耦合求解,在笛卡儿坐标系中运用二阶精度有限体积算法,保持流场守恒条件下,通过采用多层网格级对笛卡儿网格嵌套细化,从而实现多介质流体界面的高精度自适应跟踪。给出的方法逻辑简单,可以大大节省CPU时间。     

可压缩多介质粘性流体的数值计算

将考虑热传导和粘性情况下的Navier Stokes方程描述的物理过程分解成3个子过程进行数值计算,即把整个流量计算分解成无粘性流量、粘性流量和热流量3部分,采用多介质流体高精度parabolic piecewise method（PPM）方法、二阶空间中心差方法和两步Rung-Kutta时间推进方法相结合进行数值计算。给出了激波管中Riemann问题和二维、三维Richtmyer-Meshkov界面不稳定性的Navier Stokes方程和Euler方程对比计算结果,显示了粘性对界面不稳定性的影响。     

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

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

欧拉坐标系下具有锐利相界面的可压缩多介质流动数值方法研究

可压缩多介质流动问题的数值模拟在国防和工业领域内均具有重要的研究价值,诸如武器设计、爆炸安全防护等,通常具有大变形、高度非线性等特点,是一项极具挑战性的研究课题. 本文提出了一种基于 Euler 坐标系的非结构网格、具有锐利相界面的二维和三维守恒型多介质流动数值方法,可用于模拟可压缩流体和弹塑性固体在极端物理条件下的大变形动力学行为. 利用分片线性的水平集函数重构出单纯形网格内分段线性的相界面,并在混合网格内构建出具有多种介质的相界面几何结构,理论上可以处理全局任意种介质、局部 3 种介质的多介质流动问题. 利用传统的有限体积格式来计算单元边界上同种介质间的数值通量,并通过在相界面法向上求解局部一维多介质 Riemann 问题的精确解来计算不同介质间的数值通量,保证了相界面上的通量守恒. 提出了一种非结构网格上的单元聚合算法,消除了由于网格被相界面分割成较小碎片、违反 CFL 条件,进而可能带来数值不稳定的问题. 针对一维多介质 Riemann 问题、激波与气泡相互作用问题、浅埋爆炸问题、空中强爆炸冲击波和典型坑道内冲击波传播问题开展了数值模拟研究,将计算结果与相关的理论、实验结果进行比对,验证了数值方法的正确性和可靠性.      

NUMERICAL SCHEME OF MULTI-MATERIAL COMPRESSIBLE FLOW WITH SHARP INTERFACE ON EULERIAN GRIDS

可压缩多介质流动问题的数值模拟在国防和工业领域内均具有重要的研究价值,诸如武器设计、爆炸安全防护等,通常具有大变形、高度非线性等特点,是一项极具挑战性的研究课题. 本文提出了一种基于 Euler 坐标系的非结构网格、具有锐利相界面的二维和三维守恒型多介质流动数值方法,可用于模拟可压缩流体和弹塑性固体在极端物理条件下的大变形动力学行为. 利用分片线性的水平集函数重构出单纯形网格内分段线性的相界面,并在混合网格内构建出具有多种介质的相界面几何结构,理论上可以处理全局任意种介质、局部 3 种介质的多介质流动问题. 利用传统的有限体积格式来计算单元边界上同种介质间的数值通量,并通过在相界面法向上求解局部一维多介质 Riemann 问题的精确解来计算不同介质间的数值通量,保证了相界面上的通量守恒. 提出了一种非结构网格上的单元聚合算法,消除了由于网格被相界面分割成较小碎片、违反 CFL 条件,进而可能带来数值不稳定的问题. 针对一维多介质 Riemann 问题、激波与气泡相互作用问题、浅埋爆炸问题、空中强爆炸冲击波和典型坑道内冲击波传播问题开展了数值模拟研究,将计算结果与相关的理论、实验结果进行比对,验证了数值方法的正确性和可靠性.     

基于边界元法的非均质薄板弯曲问题的解

本文用边界元法研究非均质无限域弹性薄板弯曲问题.在数值实施过程中,对于夹杂和基体分别形成边界积分方程.通过离散边界积分方程,得到相应的方程组,然后结合界面条件,最终获得问题的求解方程组.在界面的相关量求得之后,可以根据需要来求解基体和夹杂中的有关位置的弯矩.数值结果与已有的解做了对比.     

Riemann solvers with evolved initial conditions

The scope of this paper is three fold. We first formulate upwind and symmetric schemes for hyperbolic equations with non‐conservative terms. Then we propose upwind numerical schemes for conservative and non‐conservative systems, based on a Riemann solver, the initial conditions of which are evolved non‐linearly in time, prior to a simple linearization that leads to closed‐form solutions. The Riemann solver is easily applied to complicated hyperbolic systems. Finally, as an example, we formulate conservative schemes for the three‐dimensional Euler equations for general compressible materials and give numerical results for a variety of test problems for ideal gases in one and two space dimensions. Copyright © 2006 John Wiley & Sons, Ltd.     

Stability of Contact Discontinuities for the 1-D Compressible Navier-Stokes Equations

In this paper, we study the large-time asymptotic behavior of solutions of the one-dimensional compressible Navier-Stokes system toward a contact discontinuity, which is one of the basic wave patterns for the compressible Euler equations. It is proved that such a weak contact discontinuity is a metastable wave pattern, in the sense introduced in [24], for the 1-D compressible Navier-Stokes system for polytropic fluid by showing that a viscous contact wave, which approximates the contact discontinuity on any finite-time interval for small heat conduction and then runs away from it for large time, is nonlinearly stable with a uniform convergence rate provided that the initial excess mass is zero. This result is proved by an elaborate combination of elementary energy estimates with a weighted characteristic energy estimate, which makes full use of the underlying structure of the viscous contact wave.     

一种在网格内部捕捉间断的Walsh函数有限体积方法

传统有限体积或有限元方法假定流动变量在单元内连续,间断仅限于控制体的交界面上,因此它们无法在控制体内部捕捉间断.本文摒弃控制体内流动变量连续的假设,将自身具有间断特点的Walsh基函数应用于有限体积方法,把控制体内的流场变量表示成间断基函数的组合形式.按照Walsh基函数在控制体内引入的间断数目和位置,将控制体单元虚分...     

Development of an adaptive discontinuity‐capturing hyperbolic finite element model

In this paper we present a five‐parameter Taylor–Galerkin finite element model to simulate Euler equations in a domain of two dimensions. The introduced free parameters are theoretically determined by employing M‐matrix theory to obtain a physically correct and non‐oscillatory solution in regions containing a sharp solution profile. To improve the computational efficiency and solution accuracy, grids are adaptively added to obtain solutions with fewer mesh points. The discontinuity‐capturing finite element model has been validated against test cases, reproducing analytical solutions to the gas dynamic problems under the current investigation. Copyright © 2004 John Wiley & Sons, Ltd.     

Two-dimensional self-similar flows generated by the interaction between a shock and low-density gas regions

A plane time-dependent flow generated by the interaction between a normal shock and a low-density gas region occupying a quarter of the plane is theoretically investigated. Numerical simulation is performed on the basis of the Euler equations. It is established that after the shock has come in contact with the low-density region two-dimensional self-similar flows of different type can develop. On regular interaction the original shock is refracted on the low-density region with the matching of the accelerated and original shock and the refracted contact discontinuity at a common point. On irregular interaction a complicated flow occurs; it includes curved and oblique shocks, a contact discontinuity with points of inflection, multiple matching points, a high-pressure jet, and a layered vortex. The jet and vortex structures are investigated in detail. The tendency of the gasdynamic structure development with variation in the control parameters of the problem is determined. A simplified, near-analytical technique for estimating the slopes of the main shocks and the gas parameters behind them is proposed.     

Finite element and implicit Runge–Kutta implementation of an acoustics–convection upstream resolution algorithm for the time‐dependent two‐dimensional Euler equations

The second of a two‐paper series, this paper details a solver for the characteristics‐bias system from the acoustics–convection upstream resolution algorithm for the Euler and Navier–Stokes equations. An integral formulation leads to several surface integrals that allow effective enforcement of boundary conditions. Also presented is a new multi‐dimensional procedure to enforce a pressure boundary condition at a subsonic outlet, a procedure that remains accurate and stable. A classical finite element Galerkin discretization of the integral formulation on any prescribed grid directly yields an optimal discretely conservative upstream approximation for the Euler and Navier–Stokes equations, an approximation that remains multi‐dimensional independently of the orientation of the reference axes and computational cells. The time‐dependent discrete equations are then integrated in time via an implicit Runge–Kutta procedure that in this paper is proven to remain absolutely non‐linearly stable for the spatially‐discrete Euler and Navier–Stokes equations and shown to converge rapidly to steady states, with maximum Courant number exceeding 100 for the linearized version. Even on relatively coarse grids, the acoustics–convection upstream resolution algorithm generates essentially non‐oscillatory solutions for subsonic, transonic and supersonic flows, encompassing oblique‐ and interacting‐shock fields that converge within 40 time steps and reflect reference exact solutions. Copyright © 2005 John Wiley & Sons, Ltd.     

Kinetic theory based multi-level adaptive finite difference WENO schemes for compressible Euler equations

In this paper we proposed the kinetic framework based fifth-order adaptive finite difference WENO schemes abbreviated as WENO-AO-K schemes to solve the compressible Euler equations, which are quasi-linear hyperbolic equations that can admit discontinuous solutions like shock and contact waves. The formulation of the proposed schemes is based on the kinetic theory where one can recover the Euler equations by applying a suitable moment method strategy to the Boltzmann equation. The kinetic flux vector splitting strategy is used in WENO-AO framework, which produces the computationally expensive error and exponential functions. Thus, to reduce the computational cost, a physically more relevant peculiar velocity based splitting strategy is used, which is more efficient than the kinetic flux vector splitting. High order of accuracy in time is achieved using the third-order total variation diminishing Runge–Kutta (TVD-RK) scheme. Several one- and two-dimensional test cases are solved for the compressible Euler equations using the proposed fifth-order WENO-AO-K schemes and the results are compared with conventional WENO-AO scheme. Proposed schemes capture the complex flow features in a smooth region accurately, and discontinuity is also well resolved. Error analysis of the proposed schemes shows optimal convergence rates in various norms.     

Upwinded finite difference schemes for dense snow avalanche modeling

Avalanche dynamics models are used by engineers and land‐use planners to predict the reach and destructive force of snow avalanches. These models compute the motion of the flowing granular core of dense snow avalanches from initiation to runout. The governing differential equations for the flow height and velocity can be approximated by a hyperbolic system of equations of first‐order with respect to time, formally equivalent to the Euler equations of a one‐dimensional isentropic gas. In avalanche practice these equations are presently solved analytically by making restrictive assumptions regarding mountain topography and avalanche flow behaviour. In this article the one‐dimensional dense snow avalanche equations are numerically solved using the conservative variables and stable upwinded and total variation diminishing finite difference schemes. The numerical model is applied to simulate avalanche motion in general terrain. The proposed discretization schemes do not use artificial damping, an important requirement for the application of numerical models in practice. In addition, non‐physical M‐wave solutions are not encountered as in previous attempts to solve this problem using Eulerian finite difference methods and non‐conservative variables. The simulation of both laboratory experiments and a field case study are presented to demonstrate the newly developed discretization schemes. Copyright © 2000 John Wiley & Sons, Ltd.     

Ideal gas flows with separation zones and time-dependent contact discontinuities of complicated shape

Some examples of flows with separation zones andmovable contact discontinuities obtained as a result of the numerical integration of the time-dependent equations for an ideal gas are presented. The examples concern a steady annular separation zone on the blunt nose of a body in a supersonic flow, periodic shedding of unsteady discontinuities from a cylinder in a steady uniform subsonic flow with a supercritical Mach number, and the complicated deformation of a contact (tangential) discontinuity, namely, the boundaries of a two-dimensional jet, either subsonic or supersonic, flowing into a cocurrent subsonic low-velocity flow. A multiple increase in the difference grid capacity in the numerical integration of the Euler equations indicates the absence of a noticeable scheme viscosity effect in the examples calculated. The inviscid nature of the separation flows obtained is also confirmed by their well-known counterparts constructed in the ideal incompressible fluid approximation. The time-average velocity fields of the two-dimensional jet and the intensity of its sound field are in reasonable agreement with the available data.     

Development of level set method with good area preservation to predict interface in two‐phase flows

A two‐step conservative level set method is proposed in this study to simulate the gas/water two‐phase flow. For the sake of accuracy, the spatial derivative terms in the equations of motion for an incompressible fluid flow are approximated by the coupled compact scheme. For accurately predicting the modified level set function, the dispersion‐relation‐preserving advection scheme is developed to preserve the theoretical dispersion relation for the first‐order derivative terms shown in the pure advection equation cast in conservative form. For the purpose of retaining its long‐time accurate Casimir functionals and Hamiltonian in the transport equation for the level set function, the time derivative term is discretized by the sixth‐order accurate symplectic Runge–Kutta scheme. To resolve contact discontinuity oscillations near interface, nonlinear compression flux term and artificial damping term are properly added to the second‐step equation of the modified level set method. For the verification of the proposed dispersion‐relation‐preserving scheme applied in non‐staggered grids for solving the incompressible flow equations, three benchmark problems have been chosen in this study. The conservative level set method with area‐preserving property proposed for capturing the interface in incompressible fluid flows is also verified by solving the dam‐break, Rayleigh–Taylor instability, bubble rising in water, and droplet falling in water problems. Good agreements with the referenced solutions are demonstrated in all the investigated problems. Copyright © 2010 John Wiley & Sons, Ltd.