水下爆炸过程的高精度数值计算

采用三阶精度的Parabolic Piecewise Method(PPM)方法和Volume Of Fluid(VOF)方法相结合，运用Lagrange-Remapping算法，数值计算爆轰产物、水、空气多种介质之间的相互作用过程，该方法可以用来处理界面两边高密度比可压缩流动以及强剪切滑移运动等问题。     

界面不稳定性实验的数值研究

采用多介质流体的三阶精度Piecewise Parabolic Method (PPM) 计算方法对界面不稳定性实验模型进行数值模拟，通过对Lawrence Livermore National Laboratary (LLNL)实验室的果冻环实验模型的数值计算，获得了与其计算和实验图像基本一致的结果， 从而验证和确认了计算方法和计算程序. 在此基础上，对于冲击波物理与爆轰物理实验 室设计的果冻内外界面10模峰对峰、峰对谷振幅为1\,mm扰动的界面不稳定性实验模型，给出 了果冻内外界面位置、速度和加速度历史曲线，详细分析了果冻内外界面不稳定性的发展、 演化过程，并给出了两种实验模型实验结果和对应的数值模拟结果.     

MOF界面重构方法在超高速运动中的应用

应用MFPPM(Multi-Fluid Parabolic Piecewise Method)方法进行超高速运动数值模拟时,由于其使用的PPM(Parabolic Piecewise Method)方法通常采用几个过渡网格来描述间断面,使得界面分辨率不高。为了更精确地对界面进行描述,将MOF(Moment-of-Fluid)界面重构方法与具有自主知识产权的MFPPM程序相结合,并首次将MOF方法运用到超高速运动中,对瞬时起爆的TNT炸药内聚压缩不同形状(三角形、正方形、正五边形、正六边形)的气腔进行了数值模拟;分别对四种构型在t=3?s和t=6.5?s时的VOF和MOF结果进行了对比。结果表明:MOF方法不仅能有效地捕捉到超高速运动中的构型几何特征,且在界面处不需要多个网格进行过渡,提高了界面重构的精度和分辨率。本文研究为界面不稳定性等问题中的复杂界面重构提供了一种新的方法。     

柱形汇聚几何中内爆驱动金属界面不稳定性

采用自研的高保真度爆轰与冲击动力学程序，对柱形汇聚几何中内爆驱动金属材料界面不稳定性的动力学行为，进行了数值模拟研究。结果表明，首次冲击后至约12 μs，界面发展以RM（Richtmyer-Meshkov）不稳定性为主；12 μs后至冲击波聚心反弹加载前，界面聚心运动处于加速减速状态，界面发展由RT (Rayleigh-Taylor)不稳定性主导；冲击波聚心反弹加载后，界面发展又由RM不稳定性主导。另外，还研究了初始条件（初始振幅、初始波长、钢壳初始厚度和几何构型）对柱形内爆驱动金属材料界面不稳定性的影响。结果显示:初始振幅较大时振幅增长也较大；初始波长较小（模数较大）时振幅增长较小，而且存在一个截止波长；钢壳厚度会抑制扰动增长，也存在一个截止厚度；几何汇聚效应会使扰动增长速度更快。     

爆炸容器内冲击波系演化及壳体响应的数值研究

对中心装药爆炸后冲击波的产生、传播和壳体动态响应全过程进行了数值研究。认为RDX瞬时爆炸 ,爆炸近场采用自相似解 ;冲击波传播和波系演化采用PPM (the Piecewise Parabolic Method)格式求解Eu ler方程 ;壳体响应采用有限元方法求解拉氏坐标系下由虚功原理得到的动力学方程。壳体内壁面边界条件分别采用强耦合和弱耦合方法处理。结果表明 :(1)当装药量相同时 ,薄壁壳体振型比厚壁壳体复杂得多 ,振幅也大 ;(2 )当装药量不同 ,壳体厚度相同时 ,爆炸场冲击波的演化过程不同 ;(3)对少量装药 ,产生的冲击波强度低 ,壳体变形小 ,是否考虑内边界运动 ,对计算结果的影响不大 ;(4 )在本文条件下 ,爆炸容器封头顶点所受的载荷最大 ,是最易发生破坏的地方 ,侧壁与爆点所在横截面的交线 ,也易破坏。     

界面不稳定性引起混合过程的二维数值计算

在可压缩多介质流体动力学高精度欧拉计算方法多介质流体分段抛物方法(multi-fluid piecewise parabolic method, MFPPM)基础上,运用算子分裂技术,增加二阶空间中心差方法和两步Rung-Kutta时间推进方法计算动力学黏性以及热流部分对流场的影响,发展适用于NS(Navier-Stokes)方程的可压缩多介质黏性流体计算方法多介质黏性流体分段抛物方法(multi-viscousity-fluid piecewise parabolic method, MVPPM). 文中采用MVPPM对英国AWE(atomic weapons establishment)激波管实验进行二维计算,给出了与实验图像基本一致的计算结果;应用MFPPM和MVPPM分别对二维柱对称内爆动力学界面不稳定性及其后期混合过程进行数值模拟,给出内外界面演化、速度历史以及后期中心气穴不同半径内因RT(Rayleigh-Taylor)界面不稳定性引起的混合量分布情况,从计算结果比较可见黏性对物质界面处混合量的分布影响明显.      

可压缩多介质粘性流体和湍流的大涡模拟

在可压缩多介质粘性流体动力学高精度计算方法MVPPM(multi-viscous-fluid piecewise parabolicmethod)基础上,引入Smagorinsky和Vreman亚格子湍流模型,采用大涡数值模拟方法求解可压缩粘性流体NS(Navier-Stokes)方程,给出适用于可压缩多介质流体界面不稳定性发展演化至湍流阶段的计算方法和二维计算程序MVFT(multi-viscosity-fluid and turbulence)。在2种亚格子湍流模型下计算了LANL(Los Ala-mos National Laboratory)激波管单气柱RM不稳定性实验,分析了气柱的形状、流场速度以及涡的特征,通过与LANL实验和计算结果的比较可知,Vreman模型略优于Smagorinsky模型,MVFT方法和计算程序可用于对界面不稳定性发展演化至湍流阶段的数值模拟。     

再冲击载荷作用下流动混合的数值模拟

基于多介质的流体体积分数VOF（volume of fluid）方法和PPM（piecewise parabolic method）方法,给出和发展了可用于多介质粘性流体动力学的数值计算方法和计算代码MVPPM（multi-viscosity-fluid piecewise parabolic method）。为了检验和验证此计算代码,对某激波管实验再冲击载荷作用下的流体动力学不稳定性及其导致的混合过程进行了数值模拟,计算结果与实验结果一致。同时还研究了激波反射冲击作用下流体混合区的演化情况,在反射激波和混合区相互作用的瞬间,混合区的宽度明显减小,之后又迅速增大;另外,混合增长率与初始扰动的频谱有很大关联。通过对有粘性（分子动力学粘性）和无粘性结果的对比,发现粘性对混合区的影响很小。     

流场非均匀性对非平面激波诱导的Richtmyer-Meshkov不稳定性影响的数值研究

基于Navier-Stokes方程组，采用可压缩多介质黏性流动和湍流大涡模拟程序MVFT (multi-viscous-flow and turbulence)，模拟了均匀流场与初始密度呈现高斯函数分布的非均匀流场中马赫数为1.25的非平面激波加载初始扰动air/SF₆界面的Richtmyer-Meshkov (RM)不稳定性现象。数值模拟结果表明，初始流场非均匀性将会影响非平面激波诱导的RM不稳定性演化过程。反射激波加载前，非平面激波导致的界面扰动振幅随着流场非均匀性增强而增大；反射激波加载后，非均匀流场与均匀流场条件下的界面扰动振幅差异有所减小。进一步，定量分析流场中环量分布及脉动速度统计量揭示了前述规律的原因。此外，还与平面激波诱导的RM不稳定性进行了简单对比，发现由于非平面激波波阵面区域的涡量与激波冲击界面时产生的涡量的共同作用，使得非平面激波与平面激波诱导的界面失稳过程存在差异。     

微细三角柱流固耦合动态特性的三维数值研究

在三维双向流固耦合计算基础上使用动网格方法,通过求解不可压缩粘性流体的N-S方程,数值研究了低速流速下(2.5m/s~15m/s)气体绕流微细三角钝体的耦合动态特性及涡激振动特性,得到了运动微细三角柱的升力系数、涡脱频率和振动振幅。观测到了微细三角柱的"锁振"和"拍"现象。接着计算了相同流场条件下的相同尺寸二维三角柱双向流固耦合涡激振动,并将结果与文献及本文中三维模型进行对比发现,三维模型的升力系数的平均幅值、涡脱频率和斯特劳哈尔数较之本文及文献中的二维数值模型较小。     

Experimental study of the Richtmyer-Meshkov instability induced by a Mach 3 shock wave

An experimental investigation of a shock-induced interfacial instability (Richtmyer-Meshkov instability) is undertaken in an effort to study temporal evolution of interfacial perturbations in the late stages of development. The experiments are performed in a vertical shock tube with a square cross-section. A membraneless interface is prepared by retracting a sinusoidally shaped metal plate initially separating carbon dioxide from air, with both gases initially at atmospheric pressure. With carbon dioxide above the plate, the Rayleigh-Taylor instability commences as the plate is retracted and the amplitude of the initial sinusoidal perturbation imposed on the interface begins to grow. The interface is accelerated by a strong shock wave (M = 3.08) while its shape is still sinusoidal and before the Kelvin-Helmholtz instability distorts it into the well known mushroom-like structures; its initial amplitude to wavelength ratio is large enough that the interface evolution enters its nonlinear stage very shortly after shock acceleration. The pre-shock evolution of the interface due to the Rayleigh-Taylor instability and the post-shock evolution of the interface due to the Richtmyer-Meshkov instability are visualized using planar Mie scattering. The pre-shock evolution of the interface is carried out in an independent set of experiments. The initial conditions for the Richtmyer-Meshkov experiment are determined from the pre-shock Rayleigh-Taylor growth. One image of the post-shock interface is obtained per experiment and image sequences, showing the post-shock evolution of the interface, are constructed from several experiments. The growth rate of the perturbation amplitude is measured and compared with two recent analytical models of the Richtmyer-Meshkov instability.PACS: 52.35.Py, 52.35.Tc     

Richtmyer–Meshkov instability on a low Atwood number interface after reshock

The Richtmyer–Meshkov instability after reshock is investigated in shock tube experiments at the Wisconsin Shock Tube Laboratory using planar laser imaging and a new high-speed interface-tracking technique. The interface is a mixture of helium and argon (50% each by volume) stratified over pure argon. This interface has an Atwood number of 0.29 and a near single-mode, two-dimensional, standing wave perturbation with an average amplitude of 0.35?cm and a wavelength of 19.4?cm. The incident shock wave of Mach number 1.92 accelerates the interface before reflecting from the shock tube end wall with M =?1.70 and accelerating the interface in the opposite direction. The amplitude growth after reshock is reported for variations in this initial amplitude, and several amplitude growth rate models are compared to the experimental growth rate after reshock. A new growth model is introduced, based on a model of circulation deposition calculated from one-dimensional gas dynamics parameters. This model compares well with the amplitude growth rate after reshock and the circulation over one-half wavelength of the interface after the first shock wave and after reshock.     

可压缩多介质粘性流体的数值计算

将考虑热传导和粘性情况下的Navier Stokes方程描述的物理过程分解成3个子过程进行数值计算,即把整个流量计算分解成无粘性流量、粘性流量和热流量3部分,采用多介质流体高精度parabolic piecewise method（PPM）方法、二阶空间中心差方法和两步Rung-Kutta时间推进方法相结合进行数值计算。给出了激波管中Riemann问题和二维、三维Richtmyer-Meshkov界面不稳定性的Navier Stokes方程和Euler方程对比计算结果,显示了粘性对界面不稳定性的影响。     

Numerical modeling of slug flow initiation in a horizontal channels using a two-fluid model

This paper presents a methodology for modeling slug initiation and growth in horizontal ducts. Transient two-fluid equations are solved numerically using a class of high-resolution shock capturing methods. The advantage of this method is that slug formation and growth in a stratified regime can be calculated directly from the solutions to the flow field differential equations. In addition, by using high-resolution shock capturing methods that do not contain numerical diffusion, the discontinuity generated by slugging in the flow field can be modeled with good accuracy. The two-fluid model is shown to be well-posed mathematically only under certain conditions. Under these circumstances, the two-fluid model is capable of correctly predicting and modeling the flow physics. When ill-posed, an unbounded instability occurs in the flow field solution, and the instability amplitude increases exponentially with decreasing mesh sizes. This work shows that there are three zones associated with slug formation. In addition, long wavelength slugs are shown to initiate from short wavelength waves. These short waves are generated at the interface of the two phases by the Kelvin-Helmholtz hydrodynamic instability. The results obtained through numerical modeling show good agreement with experimental results.     

Numerical study of three-dimensional developments of premixed flame induced by multiple shock waves

The three-dimensional interactions of a perturbed premixed flame interface with a planar incident shock wave and its reflected shock waves are numerically simulated by solving the compressible, reactive Navier–Stokes equations with the high-resolution scheme and a single-step chemical reaction. The effects of the initial incident shock wave strength (Mach number) and the initial perturbation pattern of interface on the interactions are investigated. The distinct properties of perturbation growth on the flame interface during the interactions are presented. Our results show that perturbation growth is mainly attributed to the flame stretching and propagation. The flame stretching is associated with the larger-scale vortical flow due to Richtmyer–Meshkov instability while the flame propagation is due to the chemical reaction. The mixing properties of unburned/burned gases on both sides of the flame are quantitatively analyzed by using integral and statistical diagnostics. The results show that the large-scale flow due to the vortical motion always plays a dominating role during the reactive interaction process; however, the effect of chemistry becomes more important at the later stage of the interactions, especially for higher Mach number cases. The scalar dissipation due to the molecular diffusion is always small in the present study and can be negligible.     

High order multiplication perturbation method for singular perturbation problems

This paper presents a high order multiplication perturbation method for sin- gularly perturbed two-point boundary value problems with the boundary layer at one end. By the theory of singular perturbations, the singularly perturbed two-point boundary value problems are first transformed into the singularly perturbed initial value problems. With the variable coefficient dimensional expanding, the non-homogeneous ordinary dif- ferential equations （ODEs） are transformed into the homogeneous ODEs, which are then solved by the high order multiplication perturbation method. Some linear and nonlinear numerical examples show that the proposed method has high precision.