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

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

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次幂律关系.     

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

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

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

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

一种新的光子多普勒速度频谱分析方法

光子多普勒速度计可给出飞层表面某一速度带内颗粒群速度随时间演化的频谱数据, 在冲击动力学实验尤其是微喷射及其混合研究中得到广泛应用. 本文提出一种新的光子多普勒频谱数据分析方法, 可推断出混合区厚度变化和前端等效颗粒尺度. 利用该方法, 对一些典型状态下喷射混合速度频谱开展分析, 获得了不同冲击压力、气体条件下颗粒度数据, 证实了气体环境下喷射颗粒的气动破碎现象, 以及破碎后尺度与初始条件的依赖性, 为喷射混合物理规律研究提供了重要依据.     

??????????????????У????

为避免用四阶紧致格式求解泊松方程所具有编程复杂和难以实现的困难，对 传统的五点二阶中心差分格式进行改进；通过增加对残差的校正计算，提出了一种新型具有 四阶精度的两步预估校正格式. 新格式虽需要增加一定的计算量，但它的格式精度高，边界 条件处理极简单，易于编程实现. 通过数值实验，结果证明上述格式的确具有易于编程和计 算精度高的优点. 预估校正格式很容易推广到其他复杂情形.     

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

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

高阶谱元区域分解算法及其在流动稳定性中的应用

以无时间分裂误差的区域分解Stokes谱元算法为基础构建整体稳定性分析方法.用Jacobian-free的Inexact-Newton-Krylov算法求解不可压缩Navier-Stokes方程的定常解,将Stokes算法的时间推进步作为Newton迭代的预处理,在此基础上采用Arnoldi方法计算大规模特征值问题,对复杂流动进行稳定性分析,该方法能统一处理定常和非定常计算,没有时间分裂误差,无需显式构造Jacobian矩阵,可以减少内存使用,降低计算量,并加速迭代收敛.对有分析解的Kovasznay流动的计算表明,高阶谱元法具有指数收敛的谱精度.对亚临界方腔对称驱动流的各种定常解的计算及其稳定性分析验证了方法的可行性.