Euler-Lagrange耦合计算的GEL方法

研究了一种Euler-Lagrange耦合数值方法——GEL（Ghost-fluid Euler-Lagrange）方法，编写了GEL二维计算程序。     

水下爆炸特性的二维数值模拟研究

 对带有3 mm左右厚铝壳或钢壳的柱形TNT、RS-211、T/Γ（35/65）以及ROT901装药在水下采用一端点起爆后,水中冲击波超压分布进行了二维数值模拟研究。对该四种炸药的每一种,分别考虑了药量为1 kg、3 kg、5 kg和8 kg四种TNT当量的柱形装药情况。柱形装药的长径比为L/D=1.5左右。计算使用的程序为二维欧拉多流体网格法流体力学计算程序MFIC。给出了水中冲击波超压的空间和时间分布;对铝壳和钢壳两种情况的结果进行了分析、比较;同时给出了缩比模型与1∶1模型之间相应物理量之间的缩比关系。     

地下强爆炸辐射流体动力学过程的二维数值模拟

为评估强爆炸的毁伤效应,利用流固耦合界面处理方法,考虑辐射输运引起的能量沉积过程,建立同时求解辐射流体力学方程和流体弹塑性方程的二维耦合计算方法,给出地下强爆炸时,爆室内辐射波传播、辐射能量沉积到岩土表层、岩土表层膨胀、冲击波在岩土介质中传播、爆室内高压压缩爆室壁、爆室壁振荡增大的完整物理过程.     

突发事件用爆炸容器设计方法及寿命评估

 提出了应付突发事件用爆炸容器的设计方法和使用寿命的评估方法，对JBG系列抗爆球罐在爆炸载荷作用后的残余变形进行了计算，评估了多次重复载荷条件下的使用寿命，残余变形的验证实验结果与计算结果基本一致。     

圆柱形水下爆炸实验容器壁部强度设计研究

 基于薄壁壳理论和水下爆炸理论,对圆柱形水下爆炸实验容器在爆炸冲击波作用下弹性范围内的壁部应变进行了理论分析和实验研究。导出了圆柱形水下爆炸实验容器在爆炸冲击波作用下壁部弹性应变与容器直径、壁厚及内部爆炸药量之间的关系,并对计算结果进行了实验验证。实验表明公式求解结果与实验结果具有较好的一致性。     

柱形容器内部爆炸流场的数值模拟

采用NND差分格式,同时利用点爆炸解作为计算初值,对圆柱形容器内部的二维轴对称非定常流场进行了数值模拟,清楚地反映了内部流场的变化过程,得到了容器壁面上各点的压力-时间关系。对计算结果进行了讨论,与文献中结果作了比较。     

水中障碍物爆炸毁伤效应的二维数值模拟

 建立了水下障碍物的爆炸毁伤效应的物理模型和计算模型。建立了守恒方程组、混凝土的动态弹塑性本构关系和破坏准则,并确定了其在大应变、高应变率及高压条件下的的屈服准则。基于多物质流体的Euler算法,用面向对象的C++语言自行编制了NM-MMIC通用多物质二维流体弹塑性程序,对不同炸点位置、不同炸药量对轨条砦爆炸破坏过程进行数值模拟,并采用自行开发的后处理软件VISC2D实现了可视化,定量地反应了障碍物迎爆面所承受的压力变化过程,确定了轨条砦毁伤的最佳药量,并对毁伤程度进行了评价,得到了爆炸破障的毁伤规律。研究结果对设计和评价反水中障碍物的弹药具有一定的参考价值。     

抗爆容器内爆炸流场数值模拟

 采用计算流体动力学中的二阶精度TVD差分格式和特殊算子分裂法,按轴对称问题,对半球顶圆柱筒密闭式抗爆容器内部爆炸流场进行了数值模拟。计算得到的容器壁面载荷分布与实验结果基本一致。不同时刻爆炸流场压力分布图像清晰地描述了容器壁面的冲击波加载过程。     

多物质流体动力学方法与结构动力学方法结合的流固耦合计算技术

建立一种基于重叠网格空间覆盖判断的流固耦合计算技术,把流场计算的Euler多物质流体动力学计算过程和Lagrange结构动力学计算过程相结合以描述流固间的相互作用.介绍该计算技术的具体实现过程,讨论提高几何判断效率、避免耦合遗漏、流固计算时间步协调等处理技术.应用建立的流固耦合计算程序模拟爆炸容器中爆炸流场对容器内部结构、容器壳体的作用过程.计算结果表明其可以反映流场发展过程、流场结构相互作用和结构的运动变形和破坏过程.     

耦合方法在超高速碰撞数值模拟中的应用

 相比于传统的数值方法，无网格光滑粒子法（Smoothed Particle Hydrodynamics，SPH）更适合模拟超高速碰撞问题，它能够清晰地模拟材料界面，并且能够克服中网格大变形导致的数值不稳定。然而在SPH方法中，为了得到t时刻一个节点上的物理量，需要先找出与它相互作用的临近节点，导致计算速度很慢，计算所需的时间远远大于有限元方法所需的时间。采用SPH与Lagrange法相耦合的方法对3D超高速碰撞过程进行了模拟研究。计算结果表明：SPH与Lagrange法耦合方法的计算结果与SPH计算结果精度接近，并且能够大大地节省计算时间。     

The Modified Ghost Fluid Method Applied to Fluid-Elastic Structure Interaction

In this work, the modified ghost fluid method is developed to deal with 2D compressible fluid interacting with elastic solid in an Euler-Lagrange coupled system. In applying the modified Ghost Fluid Method to treat the fluid-elastic solid coupling, the Navier equations for elastic solid are cast into a system similar to the Euler equations but in Lagrangian coordinates. Furthermore, to take into account the influence of material deformation and nonlinear wave interaction at the interface, an Euler-Lagrange Riemann problem is constructed and solved approximately along the normal direction of the interface to predict the interfacial status and then define the ghost fluid and ghost solid states. Numerical tests are presented to verify the resultant method.     

强爆炸火球数值模拟中的算子分裂方法

采用算子分裂方法将辐射流体力学方程组分裂为对流项和刚性源项,设计一种高效求解刚性源项方程组的数值方法.数值实验表明:该方法对时间步长不敏感,计算精度能够满足工程计算要求.在不对方程组做任何近似的情况下,数值给出了较长时间内火球冲击波阵面压力峰值及阵面位置的变化,结果与实际强爆炸中的经验公式吻合较好.     

用SPH-MD耦合算法模拟微通道中的液氩流动

本文构造了SPH(光滑粒子动力学)-MD(分子动力学)耦合程序,在不同宽度微通道内,对不同驱动力下的液氩流动进行了模拟计算。结果表明,在通道宽度较窄时,流体速度分布和黏性系数分布均表现出很强的微尺度效应,随着通道宽度增大,这种微尺度效应逐渐减弱,流动特征逐渐接近于宏观Poiseuille流。     

Minimal Coupling Method and the Dissipative Scalar Field Theory

Quantum field theory of a damped vibrating string as the simplest dissipative scalar field is investigated by it's coupling to an infinite number of Klein–Gordon fields as the environment by introducing a minimal coupling method. Heisenberg equation containing a dissipative term proportional to the velocity is obtained for a special choice of coupling function and quantum dynamics for such a dissipative system is investigated. Some kinematical relations is calculated by tracing out the environment degrees of freedom. The rate of energy flowing between the system and it's environment is obtained.     

高速碰撞数值计算中的光滑粒子法

扼要讨论了改进的光滑粒子法的离散思想,给出了二维轴对称问题中连续介质力学守恒方程的离散过程及离散格式,提出了轴对称坐标下确定影响域内粒子数的方法.最后通过高速碰撞的系列算例说明,光滑粒子法不但适宜于计算大变形冲击力学问题,而且有着其它网格法所无法替代的优势.     

阻抗梯度飞片加载下的超高速发射二维数值模拟方法

阻抗梯度飞片准等熵加载和超高速发射的二维数值模拟,在计算方法上集中体现了多介质、多界面、大变形、高密度比等特点.采用多介质流体高精度PPM计算方法,以VOF为基础研制MFPPM2二维计算程序,数值模拟Sandia实验室的超高速发射实验模型,获得了与Sandia实验室数值计算一致的结果.     

内部爆炸载荷作用下混凝土动力响应的数值模拟

 目前,采用实验方法对爆炸载荷破坏效应进行全面、准确描述仍是一个难题。运用ANSYS/LS-DYNA非线性显式动力学有限元程序,采用流固耦合算法,对混凝土在内部爆炸载荷作用下的空腔形成和发展规律等动力响应问题进行了数值模拟。通过实际算例,比较研究了混凝土敷设钢板与否对成坑大小、损伤区、压力云图以及应力时程曲线等的影响,进一步研究了有钢板加固时混凝土在爆炸荷载作用下的破坏机理。研究表明：采用流固耦合算法可以避免由于物质的大变形引起的单元畸变和网格消蚀等问题,提高了计算效率;很好地模拟了混凝土内部爆炸时,中心爆炸空腔压缩面和自由面的拉伸层裂破坏过程,且与理论计算值吻合较好;通过表面敷设钢板,可以减小混凝土的拉应力峰值,抑制混凝土中裂纹的发展,从而提高其防震塌能力以及结构的完整性。     

基于MPS方法的汽泡凝结数值模拟研究

采用移动粒子半隐式(MPS)方法对静止过冷水中单个汽泡的凝结现象进行了数值模拟,研究了不同初始压力和初始直径时饱和蒸汽汽泡凝结过程,获得了凝结过程中汽泡形状、当量直径和压力的变化规律;汽泡初始压力为0.1～0.5MPa,初始直径为2mm、3mm和5mm;过冷水压力为0.1MPa,温度为70℃。结果表明两相界面不存在压差时,凝结过程中汽泡始终保持球形,汽泡当量直径逐渐减小,压力近似不变;相界面存在压差时,凝结过程中汽泡从球形逐渐变为心形、半月形,汽泡当量直径和压力会出现周期性振荡,且初始压力越大振荡幅度越大。     

液滴撞击液膜的流热耦合界面追踪方法数值模拟

采用界面追踪方法研究蒸馏过程中液滴撞击高温薄液膜的流动和传热特性,将数值结果与解析解和实验进行比较验证模型的正确性,研究气液界面和热流分布的演变过程.同时,分析液滴We数和无量纲液膜厚度对传热的影响.液滴撞击后的热流密度分布显示：液膜可分为撞击区、过渡区和静态区.由于液滴的撞击作用,强制对流是撞击区内主要的传热机制.增大液滴的韦伯数或减小无量纲液膜厚度会加强热量传递.随着液滴韦伯数的增加,冲击引起的扰动增强,在动量和能量共同作用下,平均热流密度明显增大,撞击区冠状水花越明显.无量纲液膜厚度越小,平均热流密度越大,且有更长的时间保持高热流密度换热.