论文部分内容阅读
在多介质流体动力学数值模拟计算方法中,Euler和Lagrange方法都有各自的优越性和局限性,如何有效利用它们的优势是多介质流体力学数值模拟计算的重要研究方向之一。本文将以Euler方法为基础的MF PPM(Piecewise-Parabolic Method)程序和以Lagrange方法为基础的DEFEL(2-D Finite Elements Code,二维流体弹塑性动力有限元)程序,根据压力和法向速度连续准则进行耦合,发展了基于level set的GEL(Ghost-fluidEuler-Lagrange)方法。该方法在处理大变形流场和小变形结构以及复杂流动与多物体相互作用等问题时具有优越性。通过对二维算例的计算结果与文献结果比较,检验了GEL方法和耦合程序的正确性,并对球形和椭球封头的爆炸容器进行了数值模拟,通过与实验结果的比较分析,表明本文程序可以比较好的处理内爆引起的壳体流固耦合问题。 本文主要研究内容包括: 1.对Euler和Lagrange两种数值方法及其耦合计算方法的特点进行了调研,比较分析了它们的异同。对Euler流场计算的PPM计算方法,以及DEFEL二维动力有限元程序的计算流程方面进行了介绍;详细阐述了GEL耦合计算的基本理论和方法,完成了GEL二维耦合程序的编制,实现了高精度的PPM计算程序和二维弹塑性流有限元程序DEFEL的连接;应用GEL程序对一维黎曼问题(Riemann problem)和二维移动边界问题(cylinder lift-off problem)进行了数值模拟计算,计算结果与文献给出的结果一致,验证了耦合Euler-Lagrange计算方法和计算程序的正确。 2.应用本文程序对球形和椭球封头的圆柱形爆炸容器在中心起爆时爆炸冲击作用下流场与容器结构的耦合作用进行了数值分析,并与实验结果作比较,考核程序用于流固耦合问题分析的可靠性。并在球坐标系下重新计算了爆炸容器的两个算例。数值结果表明:本文程序较好的解决了炸药爆炸气体产物膨胀问题,给出了爆炸流场的冲击波传播和反射、壳体的结构振动计算结果,这些结果真实的反映了爆炸容器的工作过程,也与实验结果符合。 3.给出了GEL耦合方法在其它方面的一些应用,包括多流体与多物体相互作用的水下爆炸问题的数值模拟,气体作用于靶托模型问题数值模拟及运动钢球壳与气体相互作用问题的数值模拟。 4.给出了主要结论,并对开展下一步研究提出建议。