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个网格之内。模拟了多种复杂的可压缩多流体流动，算例结果表明本文方法可以有效地处理接触间断、激波等物理问题，且具有耗散小精度高的特点。     

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

采用数值模拟方法，研究了高度和宽度比为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个网格之内；一维和二维算例表明，该数值方法可以有效地处理接触间断、激波和多维滑移线等物理问题，并能够比其它多流体数值方法更精细地模拟多流体交界面．     

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.     

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

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

Local Heating Effect of Flow Past a Circular Cylinder

Heating effects of air flows past a two-dimensional circular cylinder at low Reynolds numbers and low Mach numbers are investigated by numerical simulation. The cylinder wall is heated partially rather than heated on the whole surface as with previous researches. The heating effects are completely different for various heating locations on the cylinder surface. Heating either windward or leeward side stabilizes the flow and reduces or completely suppresses vortex shedding from the cylinder at supercritical Reynolds numbers, which is consistent with previous results of heating on the whole surface of the cylinder. However, as the lateral sides of the cylinder （perpendicular to the stream-wise direction） are heated, an adverse effect is found for the first time in that the flow is destabilized and vortex shedding can be excited at subcritical Reynolds numbers. As the lateral sides of the cylinder are cooled, the flow is stabilized.     

柱状交界面Richtmyer-Meshkov不稳定性的数值研究

采用基于有限体积法的二阶Godunov格式模拟了柱形密度交界面在同轴激波加速下的演化过程。得到了以下初步结果:在参数相同的情形下,内聚激波比中心爆炸波对界面的扰动更危险;内聚激波从轻质流体进入重质流体与从重质流体进入轻质流体相比,界面有更快的增长。周向波数大小对界面增长率有很大的影响,在计算的参数下,n=8～12有最大增长率,大于和小于这个范围的周向波长,增长率均明显减弱。     

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

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