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

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

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

采用非结构动网格方法对含多介质的流场进行数值模拟.采用改进的弹簧方法来处理由于边界运动而产生的网格变形.采用基于格心的有限体积方法求解守恒型的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）性质的三阶龙格-库塔法进行时间积分.该格式具备间断捕捉、程序实现简单、数值稳定性强、海岸动边界以及波浪破碎处理方便和可调参数少等优点.利用典型算例对数值模型进行了验证,计算结果与实验数据吻合较好.     

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

建立了求解二维全非线性布氏（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迭代进行求解.通过两个二维模型注水驱油的数值模拟,验证了本文方法的有效性.     

基于WENO重构的真正多维黎曼求解器

具有良好守恒性与网格适应性的有限体积格式在流体力学的数值计算中占有重要地位。其中，求解数值流通量是实施有限体积法的关键步骤。一维情形下，通过求解局部黎曼问题来获得数值流通量的相关理论已经比较成熟。但是在计算多维问题时，传统的维度分裂方法仅考虑沿界面法向传播的信息，这不仅影响格式的精度，还可能会造成数值不稳定性从而诱发非物理现象。本文基于对流-压力通量分裂方法来构造真正多维的黎曼求解器，通过求解网格顶点处的多维黎曼问题来实现格式的多维特性。采用五阶WENO重构方法来获得空间的高阶精度，时间离散采用三阶TVD龙格-库塔格式。一系列数值实验的结果表明，真正多维的黎曼求解器不仅具有更高的分辨率还能有效克服多维强激波模拟中的数值不稳定性。     

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时间。     

A fast nested multi‐grid viscous flow solver for adaptive Cartesian/Quad grids

A nested multi‐grid solution algorithm has been developed for an adaptive Cartesian/Quad grid viscous flow solver. Body‐fitted adaptive Quad (quadrilateral) grids are generated around solid bodies through ‘surface extrusion’. The Quad grids are then overlapped with an adaptive Cartesian grid. Quadtree data structures are employed to record both the Quad and Cartesian grids. The Cartesian grid is generated through recursive sub‐division of a single root, whereas the Quad grids start from multiple roots—a forest of Quadtrees, representing the coarsest possible Quad grids. Cell‐cutting is performed at the Cartesian/Quad grid interface to merge the Cartesian and Quad grids into a single unstructured grid with arbitrary cell topologies (i.e., arbitrary polygons). Because of the hierarchical nature of the data structure, many levels of coarse grids have already been built in. The coarsening of the unstructured grid is based on the Quadtree data structure through reverse tree traversal. Issues arising from grid coarsening are discussed and solutions are developed. The flow solver is based on a cell‐centered finite volume discretization, Roe's flux splitting, a least‐squares linear reconstruction, and a differentiable limiter developed by Venkatakrishnan in a modified form. A local time stepping scheme is used to handle very small cut cells produced in cell‐cutting. Several cycling strategies, such as the saw‐tooth, W‐ and V‐cycles, have been studies. The V‐cycle has been found to be the most efficient. In general, the multi‐grid solution algorithm has been shown to greatly speed up convergence to steady state—by one to two orders. Copyright © 2000 John Wiley & Sons, Ltd.     

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.