风沙运动问题的SPH-FVM耦合方法数值模拟研究

针对离散颗粒模型和欧拉-欧拉双流体模型在求解气粒两相流动问题中存在的不足,提出了一种新方法——SPH-FVM耦合方法,并应用于风沙运动过程的数值模拟计算.新方法基于拟流体模型,采用光滑粒子流体动力学方法(SPH)对离散颗粒相进行求解,追踪单颗粒运动轨迹,采用有限体积法(FVM)求解连续气体相,捕捉流场特性,两相间通过曳力、压力梯度、体积分数等参量进行耦合,建立了两种不同坐标系下方法间的耦合框架.对SPH粒子所承载的物质属性进行了重定义,改造成了适用于离散颗粒相求解的光滑离散颗粒流体动力学模型(SDPH),阐述了SPH粒子与离散相中颗粒之间的关系,推导得到了拟流体SPH离散方程组.模拟了风沙运动中沙粒跃移过程和自由来流风速作用下沙粒的运动过程以及沙丘在风力作用下缓慢向前蠕移的过程,分析了颗粒的运动轨迹,流场中沙粒水平速度分布规律以及气体速度场在沙粒反作用下的变化情况等,与实验结果相符合,结果表明该方法不仅精度较高,而且计算量较小,适于求解风沙运动问题乃至其他气粒两相流动问题.     

不同界面组分分布对Richtmyer-Meshkov不稳定性的影响

基于二维非定常Euler方程,对平面激波与不同界面组分分布下氦气气柱作用过程所引起的Richtmyer-Meshkov不稳定性现象进行了数值模拟,探讨了激波冲击轻质气柱后气柱界面形态的演变及流场波系结构,定量分析了气柱特征尺度(气柱长度、高度和中轴宽度)和气柱体积压缩率随时间变化.此外,结合流场压强、速度、环量和气体混合率,多角度分析了激波驱动界面气体混合的流动机制,获得了不同界面组分分布对界面不稳定性的影响.结果表明,随着气柱界面从完全扩散界面向间断界面的过渡,界面两侧的声反射系数随之增大,使入射激波与气柱界面的作用由常规透射转变为非常规透射,反射激波逐渐加强,透射激波逐渐减弱,使得Richtmyer-Meshkov不稳定性随之增强;同时,界面两侧阿特伍德数的增大,加强了Rayleigh-Taylor不稳定性和Kelvin-Helmholtz不稳定性的发展.此外,界面不稳定性的加强使得流场环量增大,导致气体混合率的增长速率随之升高.     

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

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

微喷颗粒与气体混合过程的数值模拟研究

冲击波作用下金属与气体界面将发生微喷混合现象,即金属表面产生的微喷射物质在气体中的输运过程.提出采用散体颗粒分布代替微喷初始状态,基于气体-颗粒两相流模型对微喷混合现象进行了模拟研究.数值模拟给出了微喷混合的动力学演化过程,分析了初始气体压力和颗粒尺寸因素对混合层的影响规律;在数值模拟中发现了微喷颗粒的气动破碎现象,这可导致颗粒尺度明显减小,成为影响微喷混合演化性质的重要物理因素.本文模拟结果与相关实验结果取得一致,初步表明,气粒两相流模型是模拟微喷混合过程的一种有效方法.     

气粒两相流传热问题的光滑离散颗粒流体动力学方法数值模拟

在气粒两相流动问题中,颗粒间以及气体与颗粒间的传热问题不可忽略.光滑离散颗粒流体动力学(SDPH)模型作为一种新的求解气粒两相流动问题的方法,已经成功应用于模拟风沙运动等问题.在此基础上,提出了SDPH方法的热传导模型,模拟了气粒两相流动问题中的热传导过程以及颗粒蒸发过程.首先引入各相的能量方程,利用有限差分与光滑粒子流体动力学一阶导数相结合的方法,处理各相内部热传导项中的二阶导数问题,基于气粒两相间温度差及对流换热系数计算颗粒与气体间的热传导量,推导得到了含热传导模型的气粒两相流SDPH计算方程组,模拟计算了圆盘形颗粒团算例及鼓泡流化床内部热传导算例,并与双流体模型计算结果进行对比,结果基本符合;其次利用离散液滴模型中的颗粒蒸发传质传热定律计算颗粒的蒸发过程,数值模拟了颗粒射流蒸发过程,并与离散颗粒模型结果进行对比,两者符合得较好,验证了该方法的准确性及实用性.     

稠密可压缩气粒两相流动中的等熵声速计算建模及物理规律

气体相与颗粒相混合流场的声速研究, 由于具有重要的基础理论价值与广泛的工程应用背景, 逐渐受到人们重视. 针对稠密可压缩气粒两相流动, 综合考虑颗粒相所占空间体积以及颗粒间相互作用, 推导给出了新的等熵声速计算公式; 新公式包含了已有的纯气体、稀疏气粒两相流情形的计算公式作为其特例, 一方面验证了公式推导的正确性, 另一方面说明新公式更具有通用性; 分析了不同颗粒质量分数条件下的声速变化规律, 相应结果与普朗特的理论分析符合, 特别对于稠密气粒两相流动工况得到了一些新的物理认识; 开展了颗粒间相互作用建模参数的物理分析, 揭示了其对气粒两相流动声速的影响机理. 本文取得的成果为稠密可压缩气粒两相流动研究以及相关工程应用提供理论支撑.     

三维射流中颗粒与流体的双向耦合作用

采用基于欧拉一拉格朗日的双向耦合模型对三维气固两相平面射流中颗粒与流体的双向耦合作用进行了直接数值模拟.在考虑颗粒相的反作用后,气相运动采用直接耦合求解,计算颗粒场时,选取Stokes数为0.1的较小颗粒,采用Lagrangian方法跟踪其运动.重点考察了颗粒相与流体相之间的相互作用,分析了不同固相载率的颗粒对流场特性以及对自身扩散的影响.模拟结果表明由于颗粒的影响,在射流入口处流场最初生成的两个大涡沿横向被拉伸,而在射流下游区域,涡结构则沿流向被拉伸;在射流的下游,颗粒降低了流场中心区域的流向平均速度,削弱了流场中心区域的湍流强度.此外,跟单向耦合相比,双向耦合情况下的颗粒分布更加均匀,并且均匀程度随固相载率的增加而增大.     

水下超声速气体射流气液两相复杂流动研究

本文实验研究了水下超声速气体射流气液两相复杂流动。利用高速摄影仪和电子相机分别实时记录了过膨胀超声速工况水下气体射流的喷射状态和整体形貌,显示了不同工况水下高速气体垂直射流的演化过程和动态不稳定性形貌。研究结果表明:在射流的初始段存在与射流内部复杂波系相关的激波反馈特性,激波反馈特性发生之前存在能量积聚的高频低幅的胀鼓过程,二者均随机发生;在射流的主体段,在气水掺混和卷吸大规模能量交换作用下,射流呈现随机的偏摆效应,并且偏摆受环境流场影响明显。     

考虑双向耦合效应的格子Boltzmann气固两相湍流模型

在流体粒子概率密度函数输运方程中考虑颗粒对流体的反作用力,发展了考虑双向耦合效应的LB气固两相流模型,引入Smagorinsky亚格子模型模拟高雷诺数气相流场.对经典后台阶气固两相流动进行模拟,气相和颗粒相速度分布与实验结果进行比较,发现考虑双向耦合效应的LB气固两相流模型结果明显优于单向耦合结果.进一步研究不同惯性颗粒在流场中的弥散特性,小颗粒(St~O(0.1))对流体的跟随性较好,在流场中分布较为均匀;而St~O(1)的颗粒难被流场涡卷吸进入涡内,呈现倾向性弥散现象;大颗粒(St~O(10))由于自身惯性进入流场涡,在流场中分布较为均匀.     

DG方法求解可压缩气固两相流动

给出求解双向耦合可压缩气固两相流的间断有限元方法,对所得到的气固两相流方程组不需要采用分裂的方法离散,对气相、颗粒相方程及其对流部分和源项可以统一处理,两相都采用基于近似Riemann解的数值通量.数值模拟低压含尘激波管内的两相非平衡流动,并与平衡流、冻结流的结果进行比较.分析颗粒相的存在对气体运动的影响,及激波后松弛区域内两相间相互作用规律.发现颗粒质量比决定两相平衡后的最终状态,而颗粒直径决定两相流从非平衡到平衡的过渡过程,即不同尺寸颗粒对应的驰豫时间、松弛距离不同.结果表明：本文提出的计算方法对求解可压缩气固两相流是可行的,为研究复杂的气固两相流动问题奠定了基础.     

An Incompressible Three-Dimensional Multiphase Particle-in-Cell Model for Dense Particle Flows

A three-dimensional, incompressible, multiphase particle-in-cell method is presented for dense particle flows. The numerical technique solves the governing equations of the fluid phase using a continuum model and those of the particle phase using a Lagrangian model. Difficulties associated with calculating interparticle interactions for dense particle flows with volume fractions above 5% have been eliminated by mapping particle properties to an Eulerian grid and then mapping back computed stress tensors to particle positions. A subgrid particle, normal stress model for discrete particles which is robust and eliminates the need for an implicit calculation of the particle normal stress on the grid is presented. Interpolation operators and their properties are defined which provide compact support, are conservative, and provide fast solution for a large particle population. The solution scheme allows for distributions of types, sizes, and density of particles, with no numerical diffusion from the Lagrangian particle calculations. Particles are implicitly coupled to the fluid phase, and the fluid momentum and pressure equations are implicitly solved, which gives a robust solution.     

爆轰波对碰驱动平面锡飞层的动力学及动载行为特征研究

强冲击下的物质变形、破坏及诱发的轻重介质混合问题,是内爆压缩科学和工程应用领域的研究重点.本文针对爆轰波对碰条件下的复杂加载动力学过程及其动载破坏形态特征,开展数值模拟研究与极曲线理论分析.设计了爆轰波对碰驱动平面锡飞层的计算模型,获得了爆轰加载动力学过程及波系相互作用物理图像,分析了锡飞层对碰区自由表面速度历史的典型特征.给出了锡飞层中折射激波对碰发生马赫反射的临界条件,解读了三波结构的传播行为,阐明了对碰区内存在"一维正冲击"区域,一维区外存在单次斜冲击向两次斜冲击过渡的复杂加载动力学过程,提出了对碰区冲击动力学模型,揭示了影响对碰区动载行为特征的机理.数值模拟结果与极曲线理论分析结果相互印证,符合较好.本文的研究成果,将为深入理解和解读对碰区特殊的物质破坏及混合现象提供重要的理论支撑.     

Ejection of heavy elements from the stellar core to the periphery of the cloud of ejecta during a supernova explosion: A possible model of the processes

The possibility of simulating the processes during supernova explosions in laboratory conditions using powerful lasers (laboratory astrophysics) is investigated. The Chandra observations of ejecta in the Cassiopeia A supernova remnant are analyzed. Based on the DIANA and NUTCY numerical codes, we have performed 1D and 2D hydrodynamic simulations of the ejecta expansion dynamics for a supernova with a mass of ～5–15 solar masses within several hundred seconds after its explosion, including an initial asymmetry. We propose a model for the explosion and expansion of ejecta that illustrates strong inhomogeneities in the distribution of material to the extent that the Fe, Si, and S material from the stellar center turns out to be ejected to the periphery, the “star turns inside out,” in agreement with observations. Based on hydrodynamic similarity criteria, we consider possible supernova-simulating laser targets that will allow one to reproduce the physical processes that take place during the explosion of an astrophysical object, such as the shock propagation through the material, the growth of hydrodynamic instabilities at the boundaries of envelopes with different densities, etc.     

微喷加载装置中爆轰波阵面的测量

 根据金属材料微喷射激光全息诊断技术的要求,提出了一种测量微喷加载装置中爆轰波波前形状的实验方法,利用快响应光纤探针阵列来测量传爆药柱加载飞片的冲击波波前形状,并根据所得的实验数据用计算机模拟出爆轰波波前的形状,给出冲击波到达飞片上表面的时间分散性曲线。     

Experimental-calculation simulation of the ejection of particles from a shock-loaded surface

The possibilities of measuring complex Pylenie, which is intended for studying the parameters of the particles ejected from the free surface of a shock-loaded material and was created at the Institute of Experimental Gas Dynamics and Detonation Physics RFNC-VNIIEF, are demonstrated. The operation of the complex is based on the following three methods, which are based on different physical principles and supplement each other: laser-optical method (macro- and microfilming), pulsed X-ray method, and piezoelectric pressure sensors. This complex is used to study the ejection of particles from the free surface of lead samples when a shock wave with a pressure of 7, 16, and 23 GPa at its front reaches this surface. The effect of the surface roughness and the pressure amplitude at the shock wave front on the quantitative characteristics of the process is shown. The calculation-theoretical simulation is performed by two-dimensional numerical calculations of gasdynamic flows, and the calculation results are used to estimate the “ejecta” characteristics in terms of the developed phenomenological model of the process.     

金属表面几何缺陷微细结构对微喷射特性的影响

基于光滑粒子流体动力学方法, 分别采用实测样品几何缺陷模型和简化V形沟槽模型对铅的微喷射过程进行了模拟. 重点分析了金属表面几何缺陷微细结构对微喷射特性的影响, 并将数值计算结果与相应的实验测量值进行对比. 结果表明, 基于实测样品几何缺陷模型计算的最快喷射速度和累积喷射量与实验测量结果符合得很好. 进一步研究发现, 在实测样品几何缺陷诱导的微喷射过程中存在“二次汇聚喷射”现象, 与单次喷射相比, 该过程会诱导产生更高的喷射速度并显著影响微喷物的空间密度分布. 这说明除了受扰动波长、深度影响外,表面几何缺陷微细结构也是影响金属微喷射过程的重要因素.     

A high-order WENO-Z finite difference based particle-source-in-cell method for computation of particle-laden flows with shocks

A high-order particle-source-in-cell (PSIC) algorithm is presented for the computation of the interaction between shocks, small scale structures, and liquid and/or solid particles in high-speed engineering applications. The improved high-order finite difference weighted essentially non-oscillatory (WENO-Z) method for solution of the hyperbolic conservation laws that govern the shocked carrier gas flow, lies at the heart of the algorithm. Finite sized particles are modeled as points and are traced in the Lagrangian frame. The physical coupling of particles in the Lagrangian frame and the gas in the Eulerian frame through momentum and energy exchange, is numerically treated through high-order interpolation and weighing. The centered high-order interpolation of the fluid properties to the particle location is shown to lead to numerical instability in shocked flow. An essentially non-oscillatory interpolation (ENO) scheme is devised for the coupling that improves stability. The ENO based algorithm is shown to be numerically stable and to accurately capture shocks, small flow features and particle dispersion. Both the carrier gas and the particles are updated in time without splitting with a third-order Runge–Kutta TVD method. One and two-dimensional computations of a shock moving into a particle cloud demonstrates the characteristics of the WENO-Z based PSIC method (PSIC/WENO-Z). The PSIC/WENO-Z computations are not only in excellent agreement with the numerical simulations with a third-order Rusanov based PSIC and physical experiments in [V. Boiko, V.P. Kiselev, S.P. Kiselev, A. Papyrin, S. Poplavsky, V. Fomin, Shock wave interaction with a cloud of particles, Shock Waves, 7 (1997) 275–285], but also show a significant improvement in the resolution of small scale structures. In two-dimensional simulations of the Mach 3 shock moving into forty thousand bronze particles arranged in the shape of a rectangle, the long time accuracy of the high-order method is demonstrated. The fifth-order PSIC/WENO-Z method with the fifth-order ENO interpolation scheme improves the small scale structure resolution over the third-order PSIC/WENO-Z method with a second-order central interpolation scheme. Preliminary analysis of the particle interaction with the flow structures shows that sharp particle material arms form on the side of the rectangular shape. The arms initially shield the particles from the accelerated flow behind the shock. A reflected compression wave, however, reshocks the particle arm from the shielded area and mixes the particles.     

Coupled radiation and turbulent multiphase flow in an aluminised solid propellant rocket engine

We present in this paper numerical simulations of coupled radiative transfer and turbulent flows at high temperature and pressure, typical of multiphase flows encountered in aluminised solid propellant rocket engines. The radiating medium is constituted of gases and of liquid or solid particles of oxidised aluminum. The turbulent flow of the gaseous phase is treated by using a four equation, low Reynolds number, boundary-layer-type turbulence model. The distributions of concentrations, temperatures, and temperature fluctuation variances of particles are calculated from a Lagrangian approach and a turbulence dispersion model. Thermal and mechanical non-equilibrium between the gas and different classes of particles is allowed. A locally one dimensional, iteratively based, radiative transfer solver is developed to compute wall fluxes and radiative source terms. It is shown that the thermal boundary layer attenuates significantly the radiative fluxes coming from the outer regions. Particle radiation is found to be much more important than gas radiation. Turbulent dispersion of particles in the boundary layer induces a decrease of particle concentration in the region of maximum turbulent kinetic energy, and then, decreases the attenuation effect of wall fluxes due to the boundary layer. The effects of turbulent temperature fluctuations are found to be small in the problem under consideration.