结构与非结构网格局部优化方法

在流体力学数值模拟中，最基本的有Lagrange方法和Euler方法。Lagrange方法可用来计算多介质系统，能够刻划多介质界面，但网格的扭曲，翻转，长宽比失调等网格大变形是一个突出问题。在Euler方法中，计算网格是固定的，但是，当系统中包含多种介质时，一定会出现在一个Euler网格中包含多种介质的情形，网格中的物理量的处理比较困难。为提高精度．一般将Lagrange方法和Euler方法结合。这时网格最优问题是一个重要的内容。     

Rayleigh-Taylor不稳定性的Runge-Kutta间断有限元模拟

 采用发展后的间断有限元方法,对Rayleigh-Taylor不稳定性进行了数值模拟。在计算中采用Level-Set方法进行界面追踪,用虚拟流体方法（Ghost Fluid Method,GFM）对界面附近物理量进行等压装配。对两个典型的Rayleigh-Taylor不稳定性算例的数值研究结果表明,采用该方法计算含有接触间断的多介质流体力学问题是有效的,在交界面附近不出现伪振荡,具有较高的分辨率。     

多介质流体界面不稳定性程序

研制了二维多介质流体程序，主要包括单介质内高精度流体力学计算，多介质混合网格内各种介质输运过程和压力驰豫平衡过程计算、实际状态方程的黎曼解计算。流体计算分别采用高分辨两步PPM（Parabolic Piecewise Method）算法、TVD（Total Variation Diminishing）算法和FCT（Flux Corrected—Transport）算法，流体界面追踪采用VOF（Volume-of-Fluid）。数值求解可压缩多流体方程组和可压缩VOF方程。二维界面追踪分别采用一阶精度Youngs方法和二阶精度Elivira方法，三维界面追踪采用一阶精度Youngs方法，     

高密度比多介质可压缩流动的PPM方法

介绍一种基于流体体积分数模型的界面捕捉方法,数值模拟高密度比及含强激波的多介质可压缩流动.将PPM方法应用到多介质体积分数形式的Euler方程组,用双波近似求解一般刚性气体状态方程Riemann问题,可方便处理界面强剪切的滑移线问题,并给出了水下激波和气泡相互作用以及多相Richtmyer Meshkov不稳定性的计算结果.     

可压缩多介质流体数值模拟中的Level-Set间断跟踪方法

针对可压缩多介质流体的数值模拟,发展了一种Level-Set间断追踪技术,用LS(Level-Set)函数追踪激波和捕捉界面,用Riemann问题解构造带状区域内的虚拟流体状态,对物理量的外推方法、间断附近虚拟流体的构造、间断推进速度的计算等问题进行了研究.最后对可压缩多介质流体一维和二维守恒律方程组进行数值模拟,数值计算采用通量重构的高精度WENO格式,计算结果令人满意.     

激波和单列水柱、气水固多界面相互作用的数值模拟

研究可压缩多介质流场的激波和多介质界面相互作用问题.在Descartes固定网格采用level-set方法追踪界面,气/气界面边界条件处理采用OGFM方法,采用修正的rGFM方法提高气/水和气/固界面处构造Riemann问题精度,将Riemann近似解得到的界面参数外推到两侧真实和虚拟流体,采用五阶WENO方法求解流场Euler方程和界面level-set方程,给出不同时刻流场数值纹影图像.结果表明：在可压缩流场嵌入固体和水、气体等目标,本文方法可较精确地分辨平面运动激波和单列水柱及包含气/气、气/水和气/固等界面作用后产生的复杂激波结构.和传统的分区与贴体变换方法不同,为Descartes网格包含多介质界面复杂流场计算提供新途径.     

一维多介质可压缩Euler方程的高精度RKDG有限元方法

采用RKDG有限元目的、Level Set目的和改进的带"Isentropic"修正的Ghost Fluid目的模拟了一维多介质可压缩Euler方程,其中Euler方程、Level Set方程和重新初始化方程都采用了三阶精度的RKDG有限元目的进行离散,并对一维两种介质可压缩流体进行了数值实验,得到了较高分辨率的计算结果.     

自适应坐标变换方法

提出求解流体力学Euler方程的自适应坐标变换方法,它们保持网格夹角近似不变、保持物质界面为Lagrange描述,且使得网格速度与流体速度的差在最小二乘意义下达到极小.所提出的新坐标系能自适应于流体流场的若干重要特性.     

Richtmyer-Meshkov不稳定性湍流混合数值模拟

采用PPM方法数值求解Euler方程；采用Shyue提出的考虑压力平衡的混合网格状态方程的处理方法，完成R-M不稳定性问题后期混合的数值模拟。界面不稳定性后期混合具有明显的三维特征，二维计算不能分辨后期混合流体团之间三维扭曲拉伸现象，因此要求三维数值模拟。另一方面，界面不稳定性后期，通过非线性作用，小尺度运动被充分激发，必须模拟从大尺度到小尺度的级串现象，因此数值模拟要求很高的空间分辨率，要求大规模数值计算。由此我们采用MPI、应用区域分解方法完成程序并行化，并行程序具有较好的可扩展性。     

基于CE/SE方法的二维Euler型多物质流体弹塑性问题计算

将CE/SE方法推广到二维固体流体弹塑性问题的数值计算,同时结合杂交粒子水平集方法追踪物质界面和合适的边界条件,提出一套完整的二维Euler型流体弹塑性计算方案.通过长钨杆侵彻装甲钢实验的数值模拟,对方法的精度和有效性进行验证.     

多介质流体力学计算的守恒型高分辨率格式

应用Lagrange坐标系下的守恒型格式计算多介质流体力学问题,在物质交界面附近采用一阶格式的数值通量,而在其余部分采用高分辨率格式的数值通量,不仅保持了高分辨率的良好性质,而且消除了一般的守恒型格式在界面附近所产生的震荡.     

多介质流模拟的Runge-Kutta控制体积间断有限元方法（英文）

构造可用于多介质流数值模拟的Runge-Kutta控制体积(RKCV)间断有限元方法.对于多介质流模拟,使用线性和非线性的Riemann问题解法器计算界面处的数值流通量.该方法是一种高精度的数值方法且可以保证流体的局部守恒.数值结果表明,即使是利用线性Riemann问题解法器的计算格式也可获得较好的数值结果.与Runge-kutta间断Galerkin方法的比较展示了本文构造算法的优势.     

界面不稳定性数值模拟中的虚拟流动方法

研究了流体界面不稳定性的一类数值模拟方法——虚拟流动方法(Ghost Fluid Method).在算法中直接针对多维问题设定虚拟区域的流动参数,在流体力学方程的计算中采用了非分裂型的高分辨SCB格式,最后利用该方法完成了R-M和R-T不稳定性问题的数值计算,得到了满意的计算结果.     

辐射扩散计算方法若干研究进展

辐射流体力学研究辐射的传输对流体运动的影响,并在此条件下研究流体的运动规律.实际应用问题中辐射流体力学所描述的是非常复杂的物理过程,数值模拟是主要的研究手段之一.模拟通常采用流体计算和辐射计算分裂求解的方法.讨论求解辐射扩散方程时迫切需要解决的一些计算方法问题,包括大变形网格上扩散计算格式与非线性迭代方法,并简要介绍部分研究进展.     

随机扰动下三维流体界面不稳定性的并行计算

对三维流体界面不稳定性的数值模拟引进了新的数值计算方法,并在MPI并行计算环境下进行了数值模拟.利用LevelSet方法确定界面位置,零水平集对应界面位置.对应离散LevelSet方程和界面两侧的两套Euler方程,借助于Ghost网格方法来完成离散.对最后网格点上的两套状态量的辨认依赖于该点的LevelSet值的符号.并进行了数值计算.     

多介质流动数值计算中的界面处理方法

 求解Riemann问题得到界面接触间断的流动状态,并以此构造带状区域的虚拟流体状态,对于多维问题设计了一种方便有效的算法。同时求解耦合的守恒形式欧拉方程组和非守恒界面捕捉方程,并用Level-Set函数捕捉界面,数值计算采用高分辨率MWENO格式。最后对可压缩多介质流动问题进行了数值模拟。     

基于Riemann解的二维流体力学Lagrange有限点无网格方法

在高维流体力学计算中,对于多介质大变形等一类问题,采用有网格方法常遇到较大的困难.针对二维问题,研究了一种无网格方法——Lagrange有限点方法:在求解区域上设置适当的离散点集,视其中每一点为流体力学Lagrange点;对于点集的任一点,确定邻点集合,并基于该点同邻点集合的联系,应用Godunov方法将流体力学Lagrange方程进行离散;考虑到算法的稳健性,方法中可设置较多邻点并采用最小二乘法.将该方法应用于典型的数值算例,取得了良好效果.     

运动激波和气泡串相互作用的初步数值模拟

通过对激波和流体界面相互作用诱导的大变形界面演化的数值模拟,验证Level set方法精确模拟多个流体界面的有效性.采用2阶迎风TVD求解欧拉方程得到流场解,采用5阶WENO求解Level set方程追踪多流体界面,采用GFM方法处理流体内界面.利用文[1]的计算结果校核本文程序.在此基础上,对运动激波和气泡串相互作用过程进行了初步数值模拟,得到了不同时刻运动激波和圆管内的两个气泡作用后的演化图象,包括压力和密度等值线分布.计算结果表明:针对推广后的多界面Level set方程,该方法仍可高质量地捕捉多个流体界面.     

应用多介质PPM方法计算斜激波与物质交界面的相互作用

利用多介质PPM方法研究斜激波与物质交界面的相互作用.采用与体积分数耦合的Euler方程组作为计算模型,用双波近似来求解一般刚性气体状态方程Riemann问题.通过体积分数的计算来获得界面的位置,在整个流场采用统一的高阶PPM格式进行计算.文中对斜激波与不同物质界面相互作用进行了数值模拟,并给出了交界面上由于斜压效应产生的涡列的演化过程,特别是强斜激波与不同物质界面的相互作用的情况.     

一维多介质可压缩流动数值方法

应用高精度界面追踪方法计算一般状态方程的多介质可压缩流动问题;应用LevelSet技术捕捉界面位置,在界面附近采用守恒数值离散,用双波近似求解一般状态方程Riemann问题,并采用统一高阶PPM格式进行内点和交界面点的计算.一维算例表明,该方法对于光滑区域以及多介质交界面具有二阶精度,能准确地模拟交界面的位置,交界面计算无数值振荡和数值耗散,并能处理一般状态方程的多介质可压缩流动问题.