Van der Waals状态方程可压缩多流体流动的PPM方法

基于流体体积分数的混合型多流体数值模型,将Piecewise Parabolic Method（PPM）方法应用于可压缩多流体流动的数值模拟,采用双波近似求解多流体van der Waals状态方程的Riemann问题．模拟高密度比且含有激波的可压缩多流体流动,典型的纯界面平移问题模拟结果表明,在接触间断的界面附近,压力和速度没有任何的振荡且界面数值耗散都被控制在2—3个网格之内;一维和二维算例表明,该数值方法可以有效地处理接触间断、激波和多维滑移线等物理问题,并能够比其它多流体数值方法更精细地模拟多流体交界面．     

基于DNS的流动稳定性数值分析方法在双柱绕流中的应用

基于非结构谱元法对横向排列双圆柱绕流在固定间距情况下进行了整体线性稳定性分析,首先用Newton方法求得某个Re数的定常解,再以Arnoldi方法处理线化小扰动方程,进一步得到定常解的若干优势特征值及其特征向量.通过计算首次给出了无量纲圆心距T/D=1.2时系统的临界ReC数,并根据超临界情况下优势特征向量重构了周期流动.     

Mie-Grneisen状态方程可压缩多流体流动的PPM方法

采用流体体积分数的混合型多流体数值模型，将piecewise parabolic method (PPM)方法应用于可压缩多流体流动的数值模拟，拓展了以前提出的模型和数值方法，使它能够处理一般的Mie-Grneisen状态方程。采用双波近似和两层迭代算法求解一般状态方程的Riemann问题；并根据多流体接触界面无振荡原则设计高精度计算格式，对典型的纯界面平移问题可以从理论上证明本算法在接触间断附近压力和速度没有振荡，而且数值模拟结果表明界面数值耗散也被控制在2~3个网格之内。模拟了多种复杂的可压缩多流体流动，算例结果表明本文方法可以有效地处理接触间断、激波等物理问题，且具有耗散小精度高的特点。     

二维拉格朗日坐标系下气粒混合双向耦合对激波流场影响的计算

喷射颗粒与气体混合是内爆压缩领域的热点和难点. 针对喷射混合中的气粒双向耦合问题, 开展了理论建模、离散算法以及颗粒反馈对激波流场的影响研究. 建立了拉格朗日计算框架下的数学模型; 给出了耦合源项的离散算法; 开展了平面及汇聚构型条件下, 气粒双向耦合的数值模拟研究; 发现了颗粒反馈导致气体激波提速现象以及气区流场物理量分布形态的改变, 初步获得了量化分析结果. 本文建立的数学模型、计算方法和获得的新的物理认识, 为深入理解喷射混合现象、解决相关工程应用问题提供了重要理论支撑.     

Ejecta from periodic grooved Sn surface under unsupported shocks

Dynamic failure and ejection characteristics of a periodic grooved Sn surface under unsupported shock loading are studied using a smoothed particle hydrodynamics method. An "Eiffel Tower" spatial structure is observed, which is composed of high-speed jet tip, high-density jet slug, longitudinal tensile sparse zone, and complex broken zone between grooves. It is very different from the spike-bubble structure under supported shocks, and has been validated by detonation loading experiments. In comparison with that under supported shocks at the same peak pressure, the high-speed ejecta decreases obviously, whereas the truncated location of ejecta moves towards the interior of the sample and the total mass of ejecta increases due to the vast existence of low-speed broken materials. The shock wave profile determines mainly the total ejection amount, while the variation of V-groove angle will significantly alter the distribution of middle-and high-speed ejecta, and the maximum ejecta velocity has a linear correlation with the groove angle.     

高瑞利数下水平自然热对流的幂律关系

采用数值模拟方法，研究了高度和宽度比为1∶10的狭长矩形腔内的水平自然热对流. 根据对瑞利数(Rayleigh数)Ra在10⁴ <Ra<10¹¹内情形的计算结果，将流动分为三个不同的区间：线性区、连续过渡区、1/5次幂律区. 虽然流量和努塞尔数(Nusselt数)Nu随瑞利数的变化都包括了三个参数演化区间，但从一个区间到另外一个参数区间的转变时并不是同步的，其中努塞尔数的转变总是超前流量的转变. 对比前人的研究发现，流量1/3次幂律的结果是由于瑞利数不够高所致. 此外，模拟结果也表明Siggers等的理论分析过高估计了热通量强度，实际的温度边界层内努塞尔数和瑞利数为1/5次幂律关系.     

两种多介质流体可压缩流动界面捕捉方法的数值研究

本文比较研究两种多介质流体可压缩流动界面捕捉方法。其中混合模型界面捕捉方法将高精度PPM格式应用到多介质体积分数形式Euler方程组，用双波近似求解一般刚性气体状态方程Riemann问题；Level-Set界面捕捉方法应用无数值耗散的Ghost-Fluid-Method，并用WENO格式求解Euler运动方程和Level-Set演化方程。本文给出了一维气液接触间断，气液Riemann问题，激波接触间断相互作用，以及二维Richtmyer-Meshkov不稳定性的计算结果。     

Reynolds数对大攻角细长旋成体绕流滚转角效应影响的数值研究

使用低耗散的Roe格式,数值模拟了Reynolds数(Re)对大攻角细长旋成体绕流滚转角效应的影响.模型头部加了几何小扰动块以引发流场的不对称.在较大的Re数(Re=10 5)下,本文的计算结果与实验是相符的,此时细长体的滚转会导致双稳态、双周期现象,即侧向力随滚转角呈现类似方波形式的双周期变化,方波中侧向力基本保持不变的状态对应于流场的正则态,且两个正则态的侧向力方向相反,方波中侧向力基本保持不变的状态对应于流场的正则态,且两个正则态的侧向力方向相反;而在较小的Re数(Re=4 000)下,如果扰动足够大,细长体的滚转将导致不同的双稳态现象,此时两个正则态的侧向力方向相同,而在较小扰动下双稳态现象不再出现;Re数更小时(Re=1 000),即使在较大的扰动下,双稳态现象也不再出现,侧向力随滚动角仍是连续变化的.本文的计算结果表明,Re数越小,流场对头部扰动的感受性越弱.     

高阶谱元区域分解算法求解定常方腔驱动流

主要利用Jacobian-free的Newton-Krylov方法求解定常不可压缩Navier-Stokes方程,将基于高阶谱元法的区域分解Stokes算法的非定常时间推进步作为Newton迭代的预处理,回避了传统Newton方法Jacobian矩阵的显式装配,节省了程序内存,同时降低了Newton迭代线性系统的条件数,且没有非线性对流项的隐式求解,大大加快了收敛速度。对有分析解的Kovasznay流动的计算结果表明,本高阶谱元法在空间上有指数收敛的谱精度,且对定常解的Newton迭代是二次收敛的。本文模拟了二维方腔顶盖一致速度驱动流,同基准解符合得很好,表明本文方法是准确可靠的。本文还考虑了Re=800时方腔顶盖正弦速度驱动流,除得到已知的一个稳定对称解和一对稳定非对称解外,还获得了一对新的不稳定的非对称解。     

单模大扰动的Richtmyer-Meshkov不稳定性

采用高精度的多介质Ghost-Fluid方法,对马赫数为1.15的激波分别作用于单模大扰动Air-CO2、Air-SF6、Air-N2和Air-He界面后的Richtmyer-Meshkov不稳定现象进行了数值研究,得到了不同时刻扰动界面的演化图像,给出了流场的密度等值线和密度纹影图,同实验结果吻合较好。给出了界面的扰动增长随时间变化的情况,并同理论模型进行了对比。对激波从轻气体进入重气体的情况,扰动增长可采用Sadot模型描述线性阶段和早期非线性阶段;对于弱激波同密度接近的气体界面的相互作用,线性阶段时间较长,可用线性模型描述。