二维GEL耦合方法及其对爆炸容器的数值模拟

来源 :中国工程物理研究院 | 被引量 : 3次 | 上传用户:choww
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
在多介质流体动力学数值模拟计算方法中,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.给出了主要结论,并对开展下一步研究提出建议。
其他文献
本文主要研究理想的准素分解与Gr(o)bner基在多项式复合下的性质与计算.  设k[x1,…,xn]是域k上关于变量x1,…,xn的多项式环,I是k[x1,…,xn]上的零维理想,I的准素分解是将I分解成有
线性互补问题广泛应用于许多科学与工程领域,如经济平衡问题、非线性方程组问题、非线性规划问题等。近几十年来,人们提出了求解线性方程组的预处理方法。其基本思想是,先构