参考学习资料 流体输配管网 ALE 译文.docx

订阅DeepLPro以编辑此文档。

访问www.DeepL.com/Pro，了解更多信息。

订阅DeepLPro以编辑此文档。

访问www.DeepL.com/Pro，了解更多信息。

高速撞击问题的ALE模拟。

M.Kucharˇ?k,J.Limpouch,R.Liska,P.Va′chal.

捷克技术大学，核科学和物理工程学院，Brˇehova7，捷克共和国布拉格。

2004年4月29日收到

我们在此介绍了我们新开发的二维四边形ALE代码。该代码被用于小型加速物体对散装目标的撞击过程中陨石坑形成的数值模拟。用拉格朗日流体力学代码对高速撞击问题进行模拟，在后期会导致拉格朗日网格严重变形，无法继续计算。然而，这种情况可以用任意拉格朗日-欧拉法（ALE）来处理。为了保持网格的质量，ALE方法在拉格朗日计算的几个时间步长后定期应用网格平滑。在改变网格后，必须将旧网格中的保守量子重新映射到新的更好的网格中。在重新映射后，可以继续进行拉格朗日计算。

PACS：52.65.Kj，47.11.+j。

关键字：等离子体流体力学、ALE方法、高速碰撞

一、简介

在自然界中，许多现象都可以理解为守恒定律。流体动力学和等离子体物理问题是由一组流体力学的质量、总能量和动量在每个方向上的守恒定律来模拟的。数学上，这些守恒定律用一组双曲偏方程来表示。一般来说，这些方程的求解必须用数值方法来完成。

人们可以使用两种不同的方法来求解流体力学方程，即欧拉法和拉格朗日法。在欧拉法的情况下，计算网格是静态的，流体在其中流动。在拉格朗日法中，计算网格是随流体移动的，它可以解决流体体积和形状变化较大的区域的解。不幸的是，它可能会发生移动网格变得纠结并失去其规律性。现在违反了数学假设，拉格朗日方法不能继续计算。

任意拉格朗日-欧拉法（arbitraryLagrangian-Eulerian，ALE）[1]是一种方法，如何处理这个概率。在经过几个时间步骤后，我们通过重新分区和重新把保守变量放到新的、更好的网格中来平滑计算网格。该算法为数值解增加了额外的di?Ll_FB00?usion，计算网格仍然合理，拉格朗日数值计算可以继续进行。

在本文中，我们介绍了逻辑正交网格上的二维ALE代码的描述，以及它在高速撞击问题上的应用。此前，高速度撞击问题已经模拟了激光PALS强脉冲与目标相互作用时的陨石坑形成过程[2]。强烈的激光束也可用于大粒子对固体目标的烧蚀加速。在此，我们用流体力学的ALE代码模拟了一个加速的固体铝盘对一个巨大目标的冲击。

捷克斯洛伐克物理学报，第54卷(2004年)，补编。CC391

M

M.Kucharˇ′?k,J.Limpouch,R.,Liska,P.Va′chal.

高速撞击问题的ALE模拟。

高速撞击问题的ALE模拟。

这里使用的模拟参数与PALS激光器实验中的条件几乎相似[3]。

描述质量、总能量和两个方向上的动量守恒的流体力学方程，可以写成这样的形式。

1dρdedv

ρdt=-?-v, ρdt=-p?-v, ρdt=-?p，（1）。

?其中t为时间，ρ为质量密度，v为速度矢量，p为压力，e为比内能。为了结束这个系统，我们必须添加状态方程（EOS），它涉及特定流体的密度、内能和压力。EOS最简单的可能性是理想气体的状态方程。它可以很容易地写成p=(γ1)sρ，其中s是内能的密度e=ρs。这种形式的状态方程可以用于许多测试计算，因为它的评估是非常快的，但它的用途是有限的等离子体模拟，因为等离子体di?Ll_FB00?ers从理想气体很多。为了近似真实的等离子体行为，我们采用了商态方程(QEOS)[4]，它在广泛的等离子体条件下有效。它通过应用Thomas-Fermi理论计算给定条件下的热力学特性，并应用Cowan模型计算离子。

?

2ALE代码

我们在逻辑正交四边形网格上的二维笛卡尔几何ALE代码在流体力学的各种问题上进行了测试，如Noh问题、Sedov问题、Rayleigh-Taylor不稳定问题。所有的数值结果都是合理的。

ALE代码的拉格朗日部分是基于在每个网格节点中计算出的区域力、次区域力和粘度力[5,6]，这些力定义了下一个时间步的新节点速度。力的分区部分取决于单元的几何形状和单元内的压力。它代表从网格单元到每个节点的力，由于内部的压力。分区压力取决于单元中的压力和单元角（单元中与特定节点相邻的四分之一）的压力之间的差值。它减少了人为的网格扭曲，防止了沙漏网格运动[6]，其值由输入参数之一的优点因子控制。因此，我们可以调节解的涡度。节点力的最后一部分，即粘度力，为溶液增加人工粘度。如果没有这部分，就没有足够的耗散力，溶液很快就会失效。我们在经典的Kurop