RICHTMYER-MESHKOV不稳定性的数值模拟

采用了自适应的非结构网格和基于有限体积法的二阶Godunov格式,数值模拟了在激波作用下两种不同密度流体的交界面的演化过程.着重讨论了Richtmyer Meshkov不稳定性以及斜压效应在交界面演化过程中的作用,并给出了交界面的扰动增长率.     

可压涡卷空间演化的迎风紧致差分数值模拟

从数值算法的耗散和色散特征的时空全离散Fourier分析出发,通过直接求解二维非定常可压Navier Stokes方程,将发展的5阶迎风紧致差分格式用于无约束可压平面受迫剪切层中基频涡卷空间演化过程的数值模拟.采用被动守恒标量等方法显示了基频涡卷的饱和、一次对并、二次对并等现象,据此探讨了入口来流亚谐扰动引起的初值效应问题,表明可压大尺度涡结构空间演化形态与受迫扰动方式之间存在关联.     

轴对称射流气动声场的数值模拟

构造了Kirchhoff积分和CFD计算相结合的算法,采用高精度差分格式对不同喷口马赫数的轴对称射流进行数值模拟,并以此作为近场声源,运用Kirchhoff积分方法研究远场气动噪声.数值结果表明,远场噪声具有方向性,噪声声压在离开对称轴20°处达到最大值.随着传播距离增大,噪声方向性逐渐减弱.     

求解双曲型守恒律方程的高精度迎风紧致群速度控制法

从迎风紧致逼进[1]出发,提出求解流体力学双曲型守恒律的一种高精度的数值方法,同时采用群速度控制方法捕捉激波。该方法在光滑区具有三阶精度。     

Rayleigh-Taylor和Richtmyer-Meshkov不稳定性的FCT方法数值模拟

介绍了二维流体力学不稳定性程序的FCT数值方法,数值模拟Rayleigh-Taylor和Richt-myer-Meshkov流体不稳定性,在线性阶段,与线性理论符合很好;在非线性阶段,与俄罗斯激波管实验的计算结果符合很好。计算结果表明:FCT方法有较高计算精度,给出了不稳定性发展的好的图象,适合于ICF烧蚀和内爆流体不稳定性问题的计算。     

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

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

超音速平板混合层被动标量的数值模拟

采用空间五阶精度紧致-WENO混合格式,数值模拟了Mc=0.2和1.0的二维空间发展超音速平板混合层被动标量混合过程.并与Clemens的实验结果进行比较,结果表明,被动标量时均值与实验结果符合很好,而被动标量脉动均方根值与实验值相比较小.还研究了可压缩性对标量混合过程的影响,随着对流马赫数的增加,被动标量时均值剖面变化不大,而脉动均方根峰值增大.     

圆柱界面不稳定性的多组分气体动理学数值计算

在假定单元内各组分同温同速的条件下,采用气体动理学格式(Gas-Kinetic Scheme,GKS)对空气/He和空气/R22圆柱界面不稳定性进行了数值计算,得到不同时刻的密度分布以及界面上特征位置的位移历史和平均速度。当激波穿过界面后,界面上特征位置的位移随时间逐渐增大,特征位置的平均速度与前人的实验结果和数值模拟结果吻合很好。对比结果表明,从微观气体运动角度出发的GKS方法对于界面不稳定性问题具有良好的模拟能力。     

有厚度平板尾缘可压缩剪切层中的涡结构数值模拟

采用一个新型Fu Ma高精度UCD5 SCD6紧致差分算法,通过直接求解二维Navier Stokes方程,成功实现了有厚度平板尾缘可压缩剪切层中涡结构的数值模拟,并考查了平板厚度对其的影响.计算对流马赫数M_c=0.3,平板厚度分别为1,2,3,4个参考长度.结果表明,增加平板厚度可促使平板尾缘可压缩剪切层中的涡提前卷起,有利于两股气流混合.     

激波作用下SF6气泡界面演化和射流发展的数值模拟

 数值研究了平面激波冲击氮气环境中SF₆气泡界面的Richtmyer-Meshkov不稳定性,重点关注其中的激波聚焦及射流的产生和发展过程。在入射激波马赫数为1.23的情况下,给出了压力、密度、数值纹影和涡量等物理量的演化图像,定量分析了流场中压力最大值、密度最大值、射流速度、环量和斜压力矩随时间的变化关系。计算结果表明,平面激波冲击SF₆气泡过程有很强的聚能效应,在气泡内部靠近下游极点处发生激波近似理想聚焦和点爆炸现象,直接导致出现二次波系以及向下游运动的细长射流结构。相比入射激波,二次波系产生斜压力矩和涡量的能力要弱得多。     

可压缩轴对称射流流场近区的数值模拟

采用高精度差分方法对二维轴对称可压缩射流流场进行了数值模拟,研究了大尺度涡结构的形成和发展过程及其在射流雪展过程中的作用,计算结果显示了射流失稳后首先出现Ｋｅｌｖｉｎ－Ｈｅｌｍｈｏｌｔｚ非稳定特征,之后向下游发展的过程中, 一笥效应的增长导致涡的配对和对并;指出射流场中压力波动的分布是大尺度涡结构分布的直接结果。     

三维均匀各向同性两相湍流的直接模拟

本文对三维气粒两相均匀各向同性湍流进行了直接数值模拟。气相控制方程组采用分布投影方法进行求解,微分方程采用六紧致阶差分格式和快速Fourier变换结合求解;计算颗粒场时,采用Lagrangian方法。由该方法得到的能谱和各统计量与由谱方法得到的对应值进行了比较,吻合十分理想,对不同Stocks数颗粒在流场内的瞬态分布也进行了初步模拟,并观察到局部富集现象,证明该方法是进行两相湍流直接数值模拟行之有效的方法。     

四阶紧致差分格式对静力模式中垂直网格计算频散性的影响

为了说明四阶紧致差分格式在大气和海洋数值模式中的潜在价值,提出一种通用方法,推导静力线性斜压适应方程组在微分和差分情况下的频散关系,水平尺度分100 km,10 km和1 km三种情况,从频率、水平群速和垂直群速方面,对采用二阶中央差和四阶紧致差分格式情况下,非跳点网格(N网格)、Lorenz网格(L网格)、Charney-Phillips网格(CP网格)、Lorenz时间跳点网格(LTS网格)和Charney-Phillips时间跳点网格(CPTS网格)的计算特性进行比较,发现采用高精度的四阶紧致差分格式总体上可以明显减少上述三种水平尺度波动在N网格、CP网格、L网格和CPTS网格上的频率、水平群速和垂直群速误差,但对LTS网格,采用四阶紧致差分格式,会使得计算水平群速和垂直群速误差变大.