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

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

基于MUSTA格式的全非线性Boussinesq波浪传播数值模型

建立了求解二维全非线性布氏（Boussinesq）水波方程的有限差分/有限体积混合数值格式. 针对守恒形式的控制方程,采用有限体积方法并结合 MUSTA格式计算数值通量, 剩余项则采用有限差分方法求解, 采用具有总变差减小（totalvariation diminishing, TVD）性质的三阶龙格-库塔法进行时间积分.该格式具备间断捕捉、程序实现简单、数值稳定性强、海岸动边界以及波浪破碎处理方便和可调参数少等优点.利用典型算例对数值模型进行了验证,计算结果与实验数据吻合较好.      

非结构动网格在多介质流体数值模拟中的应用

采用非结构动网格方法对含多介质的流场进行数值模拟.采用改进的弹簧方法来处理由于边界运动而产生的网格变形.采用基于格心的有限体积方法求解守恒型的ALE(Arbitrary Lagrangiall-Eulerian)方程,控制面通量的计算采用HLLC(Hartem,Lax,van Leer,Contact)方法,采用几何构造的方法使空间达到二阶精度,时间离散采用四阶Runge-Kutta方法.物质界面的处理采用虚拟流体方法.本文对含动边界的激波管、水下爆炸等流场进行数值模拟,取得较好的结果,不同时刻界面的位置和整个扩张过程被准确模拟.     

A FULLY NONLINEAR BOUSSINESQ MODELFOR WAVE PROPAGATION BASED ON MUSTA SCHEME

建立了求解二维全非线性布氏（Boussinesq）水波方程的有限差分/有限体积混合数值格式. 针对守恒形式的控制方程,采用有限体积方法并结合 MUSTA格式计算数值通量, 剩余项则采用有限差分方法求解, 采用具有总变差减小（totalvariation diminishing, TVD）性质的三阶龙格-库塔法进行时间积分.该格式具备间断捕捉、程序实现简单、数值稳定性强、海岸动边界以及波浪破碎处理方便和可调参数少等优点.利用典型算例对数值模型进行了验证,计算结果与实验数据吻合较好.     

全机绕流Euler方程多重网格分区计算方法

全机三维复杂形状绕流数值求解只能采用分区求解的方法,本文采用可压缩Euler方程有限体积方法以及多重网格分区方法对流场进行分区计算。数值方法采用改进的van Leer迎风型矢通量分裂格式和MUSCL方法,基于有限体积方法和迎风型矢通量分裂方法,建立一套处理子区域内分界面的耦合条件。各个子区域之间采用显式耦合条件,区域内部采用隐式格式和局部时间步长等,以加快收敛速度。计算结果飞机表面压力分布等气动力特性与实验值进行了比较,二者基本吻合。计算结果表明采用分析“V”型多重网格方法,能提高计算效率,加快收敛速度达到接近一个量级。根据全机数值计算结果和可视化结果讨论了流场背风区域旋涡的形成过程。     

顺序输送管道水力瞬变模拟的有限体积法

针对顺序输送水力瞬变中输送介质不连续的情况,提出采用有限体积法取代传统的特征线法进行求解,以减少插值误差。首先得到水击基本方程基于压力-流速的非守恒有限体积离散格式,通过Rieman求解器得到控制体界面的Rieman解。采用MUSCL-Hancock Primitive方法进行界面数值重构与时间推进,构建时空二阶精度的TVD格式。边界的处理采用Rieman不变量与构造虚拟边界节点相结合的方法,使其与整体精度保持一致。数值计算与实验对比,证明了本文算法能有效抑制非物理振荡,具有精度高、鲁棒性好的特点。     

基于聚合多重网格方法的跨音速颤振的数值模拟

采用流固耦合方法对跨音速颤振进行了数值模拟。流体方面在非结构网格上用有限体积方法求解了Euler方程;结构方面则求解了后掠机翼典型剖面的结构模态方程。时间推进采用双时间步长：对每一真实时间步,都通过基于聚合多重网格方法的伪时间步推进,对流体和结构方程交替迭代．得到一个稳态的流固耦合的解。文章最后给出了NACA64A010翼型剖面的跨音速颤振边界．与相关文献的计算结果符合良好。     

考虑磁头表面高度不连续性气膜润滑的数值模拟与有效算法

随着磁头滑块的飞行高度不断降低,给气体润滑方程的数值求解带来了诸如计算时间过长、甚至计算发散等方面的问题。为了获得1Tbit/in2的存储密度,磁头滑块尾部的最小飞行高度接近1.5nm。本文基于作者提出的修正气膜润滑方程的线性流率(LFR)模型,考虑磁头滑块表面高度的不连续性,建立了基于有限体积法的气膜润滑方程离散格式,并把网格自适应技术与多重网格法应用到离散方程的迭代算法中,发展了可模拟最小飞行高度为0.5nm时磁头滑块压力分布的数值模拟方法与有效算法。文中以一个具有复杂表面形状的磁头滑块为例,检验了计算方法与算法的有效性。数值结果表明:在磁头滑块最小飞行高度较低时,必须要考虑滑块表面高度的不连续性,否则就得不到收敛的数值计算结果;与FK-Boltzmann模型相比,LFR模型具有较高的计算效率,采用网格自适应技术与多重网格法能有效地提高求解气膜润滑方程的计算效率。     

缝洞型多孔介质中多相流的有限体积法数值模拟

本文研究的碳酸盐岩油藏储集体属于缝洞型多孔介质.这类缝洞型多孔介质由裂缝、溶蚀孔洞和低孔隙度低渗透率的基岩组成.裂缝是空隙流体流动的主要通道;溶蚀孔洞大小从几厘米到数米不等,渗透率和孔隙度都很高,是流体主要的储集空间.由于缝洞型多孔介质空隙空间的复杂性和强非均质性,数值计算中基本控制方程的空间离散应采用非结构化网格的计算模型.本文采用有限体积法模拟缝洞型多孔介质中多相流体的流动,并给出了相应的单元中心格式有限体积法的计算公式.裂缝介质和溶洞介质中单元间多相流体的流动考虑为高速非达西流,其质量通量采用Forchheimer定律计算.非线性方程的离散选取全隐式格式,并采用Newton-Raphson迭代进行求解.通过两个二维模型注水驱油的数值模拟,验证了本文方法的有效性.     

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

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

A numerical study on a Cartesian-based body-fitted adaptive grid method

ABSTRACT

A hybrid Cartesian-based body-fitted adaptive grid method for compressible Navier–Stokes equations is implemented and investigated. In this method, the body-fitted structured grids are generated around the geometries, and the left regions are filled with Cartesian grids. To transfer the data between the different grids, the donor cell searching technique is adopted. An unstructured data-based finite volume update procedure is used, and least squares method is suggested to retain the second order in the overlap region. The moving shock waves with different speeds and vortex passing through the interfaces of the hybrid Cartesian grid are used to explore the accuracy and conservation. A new technique is presented to deal with the non-physical stagnation of slowly moving shock wave around the interface of grid. Numerical examples are presented to demonstrate the results. The three-dimensional extension has also been shown by a benchmark problem.     

动态混合网格生成及隐式非定常计算方法

建立了一种基于动态混合网格的非定常数值计算方法. 混合网格由贴体的四边形网格、外场 的多层次矩形网格和中间的三角形网格构成. 当物体运动时，贴体四边形网格随物体运动而 运动，而外场的矩形网格保持静止，中间的三角形网格随之变形；当物体运动位移较大，导 致三角形网格的质量降低，甚至导致网格相交时，在局部重新生成网格. 新网格上的物理量 由旧网格上的物理量插值而得. 为了提高计算效率，采用了双时间步和子迭代相结合的隐式 有限体积格式计算非定常Navier-Stokes方程. 子迭代采用高效的块LU-SGS方法. 利用该 方法数值模拟了NACA0012振荡翼型的无黏和黏性绕流，得到了与实验和他人计算相当一致 的结果.     

An unstructured/structured multi‐layer hybrid grid method and its application

A multi‐layer hybrid grid method is constructed to simulate complex flow field around 2‐D and 3‐D configuration. The method combines Cartesian grids with structured grids and triangular meshes to provide great flexibility in discretizing a domain. We generate the body‐fitted structured grids near the wall surface and the Cartesian grids for the far field. In addition, we regard the triangular meshes as an adhesive to link each grid part. Coupled with a tree data structure, the Cartesian grid is generated automatically through a cell‐cutting algorithm. The grid merging methodology is discussed, which can smooth hybrid grids and improve the quality of the grids. A cell‐centred finite volume flow solver has been developed in combination with a dual‐time stepping scheme. The flow solver supports arbitrary control volume cells. Both inviscid and viscous flows are computed by solving the Euler and Navier–Stokes equations. The above methods and algorithms have been validated on some test cases. Computed results are presented and compared with experimental data. Copyright © 2006 John Wiley & Sons, Ltd.     

Simulation of sharp gas–liquid interface using VOF method and adaptive grid local refinement around the interface

The volume of fluid (VOF) method is used to perform two‐phase simulations (gas–liquid). The governing Navier–Stokes conservation equations of the flow field are numerically solved on two‐dimensional axisymmetric or three‐dimensional unstructured grids, using Cartesian velocity components, following the finite volume approximation and a pressure correction method. A new method of adaptive grid local refinement is developed in order to enhance the accuracy of the predictions, to capture the sharp gas–liquid interface and to speed up the calculations. Results are compared with experimental measurements in order to assess the efficiency of the method. Copyright © 2004 John Wiley & Sons, Ltd.     

Calculation of laminar recirculating flows using a local non-staggered grid refinement system

A grid-embedding technique for the solution of two-dimensional incompressible flows governed by the Navier-Stokes equations is presented. A finite volume method with collocated primitive variables is employed to ensure conservation at the interfaces of embedding grids as well as global conservation. The discretized equations are solved simultaneously for the whole domain, providing a strong coupling between regions of different refinement. The formulation presented herein is applicable to uniform or non-uniform Cartesian meshes. The method was applied to the solution of two scalar transport equations, to cavity flows driven by body and shear forces and to a sudden plane contraction flow. The numerical predictions are compared with the exact solutions when available and with experimental data. The results show that neither the convergence rate nor the stability of the method is affected by the presence of embedded grids. Embedded grids provide a better distribution of grid nodes over the computational domain and consequently the solution accuracy was improved. The grid-embedding technique proved also that significant savings in computing time could be achieved.     

A structured but non‐uniform Cartesian grid‐based model for the shallow water equations

In large‐scale shallow flow simulations, local high‐resolution predictions are often required in order to reduce the computational cost without losing the accuracy of the solution. This is normally achieved by solving the governing equations on grids refined only to those areas of interest. Grids with varying resolution can be generated by different approaches, e.g. nesting methods, patching algorithms and adaptive unstructured or quadtree gridding techniques. This work presents a new structured but non‐uniform Cartesian grid system as an alternative to the existing approaches to provide local high‐resolution mesh. On generating a structured but non‐uniform Cartesian grid, the whole computational domain is first discretized using a coarse background grid. Local refinement is then achieved by directly allocating a specific subdivision level to each background grid cell. The neighbour information is specified by simple mathematical relationships and no explicit storage is needed. Hence, the structured property of the uniform grid is maintained. After employing some simple interpolation formulae, the governing shallow water equations are solved using a second‐order finite volume Godunov‐type scheme in a similar way as that on a uniform grid. Copyright © 2010 John Wiley & Sons, Ltd.     

基于非结构化同位网格的SIMPLE算法

通过基于非结构化网格的有限体积法对二维稳态Navier—Stokes方程进行了数值求解。其中对流项采用延迟修正的二阶格式进行离散；扩散项的离散采用二阶中心差分格式；对于压力-速度耦合利用SIMPLE算法进行处理；计算节点的布置采用同位网格技术，界面流速通过动量插值确定。本文对方腔驱动流、倾斜腔驱动流和圆柱外部绕流问题进行了计算，讨论了非结构化同位网格有限体积法在实现SIMPLE算法时，迭代次数与欠松弛系数的关系、不同网格情况的收敛性、同结构化网格的对比以及流场尾迹结构。通过和以往结果比较可知，本文的方法是准确和可信的。     

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

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

Numerical simulation of free‐surface flow using the level‐set method with global mass correction

A new numerical method that couples the incompressible Navier–Stokes equations with the global mass correction level‐set method for simulating fluid problems with free surfaces and interfaces is presented in this paper. The finite volume method is used to discretize Navier–Stokes equations with the two‐step projection method on a staggered Cartesian grid. The free‐surface flow problem is solved on a fixed grid in which the free surface is captured by the zero level set. Mass conservation is improved significantly by applying a global mass correction scheme, in a novel combination with third‐order essentially non‐oscillatory schemes and a five stage Runge–Kutta method, to accomplish advection and re‐distancing of the level‐set function. The coupled solver is applied to simulate interface change and flow field in four benchmark test cases: (1) shear flow; (2) dam break; (3) travelling and reflection of solitary wave and (4) solitary wave over a submerged object. The computational results are in excellent agreement with theoretical predictions, experimental data and previous numerical simulations using a RANS‐VOF method. The simulations reveal some interesting free‐surface phenomena such as the free‐surface vortices, air entrapment and wave deformation over a submerged object. Copyright © 2009 John Wiley & Sons, Ltd.     

基于自适应直角网格的二维全速势方程有限体积解法A finite volume method for 2D Full-Potential equation on adaptive Cartesian grids

为满足亚声速和跨声速飞机概念设计中快速气动计算的需求,研究和发展一种基于自适应直角网格的非线性全速势方程有限体积解法。要点如下。（1）在几何自适应直角网格的基础上,使用结合单元融合的网格切割算法处理物面边界,提出一种修正非贴体切割网格的方法。（2）采用隐式格式结合GM RES算法求解该非线性位流方程,针对流场的自适应来捕捉激波。（3）采用镜像法处理物面边界处的无穿透条件,并提出解析的方法来修正镜像单元的值。（4）针对直角网格的特点,提出在库塔线上插入库塔单元的方法施加库塔条件。NACA0012翼型绕流的算例结果表明,该方法用于亚声速和跨声速气动计算能得到令人满意的结果,且自动化程度高、收敛速度快。