MLSPH方法在内爆问题中的应用.pdfVIP

计算机与计算数学 155 6．6 水平集方程在四边形网格上的数值离散方法 程俊霞 在多介质大变形数值计算中，正交的笛卡儿网格很难满足计算要求，经常要在一般四边形网格上进行 数值计算，因此需要在一般四边形网格上实现波阵面的捕捉算法。水 平集方法不需显式地构造界面，容易处理界面拓扑结构的变化，因此 选用水平集方法来捕捉波阵面。无曲率的水平集方程、重新初始化方 程是齐次Hamilton．Jacobi方程，可以直接运用正格式进行离散。文中 只考虑曲率项是曲率的线性函数的情况，与曲率相关的水平集方程含 有水平集函数的二阶导数，不能采用正格式离散，采用伽辽金有限元 方法进行离散。与曲率相关的水平集方程是一个具有抛物性质的方 程，在求解时一般应采用隐式求解，但隐式求解的计算量非常大，为 了避免这个缺点，采用迭代格式求解 曲率相关的水平集方程。在正交 图 1 与曲率相关的水平集方程的 均匀的笛卡儿网格上，在厶 ，三2范数意义下，无曲率的水平集方程的 零等值线的演化过程 数值离散可以达到二阶精度，与曲率相关的水平集方程的数值离散可以达到强一阶精度。在西-o．5，0．5)， y(-o．5，0．5)的正方形区域，零等值线为圆心在原点，半径为0．4的圆， =一a a=1．0。80~80笛卡儿网 ． 格的与曲率相关的水平集方程的计算结果如图 1所示。从图l可以看出，零等值线随时间逐渐向圆心运动， 没有出现不稳定现象。此离散格式可用于一般的凸四边形网格的水平集方程的求解。 6．7 MLSPH方法在内爆问题中的应用 温万治 董师舜 杭旭登 在已有MLSPH(MovingLeastSquaresParticleHydrodynma ics)方法研究的基础上，进一步研究MLSPH 方法在流体力学计算方面的应用，探索它模拟计算内爆的适用性。 在二维坐标系中，移动最d-乘近似选取二次多项式基函数向量，权函数采用各向同性的3次B样条 函数，离散格式采用Dilts推导的一种局部守恒弱格式，求解方法与 SPH方法类同。为了减小激波间断处 产生的非物理振荡，在两个粒子之间引入人为黏性压力项。对时间的积分，采用预估矫正的二阶精度格式。 粒子间相互作用的交叉面积向量用三次样条积分公式计算。同时暂不考虑边界向量的计算，固壁边界采用 传统的虚粒子方法处理。 二维内爆压缩过程的计算模型为柱对称磁压缩内爆模型，此模型结构简单，只有一个圆筒形的铜壳。 初始时，铜壳具有不可压分布的向心速度，铜壳外界面压力为零，壳内是一个磁通量守恒的磁场，磁场压 力pH(D确 (0) (0) (r)]，作为铜壳内界面的压力边界条件，铜壳内部的磁场压力随铜壳内爆收缩而升高， 其中瑚(0)是初始的磁场压力，A(0是铜壳内部的面积。模型的初始几何、材料参数如下：初始铜壳内外半 径为4．8，5cm；外界面径向速度为-0．2cm／ps；铜壳用简化的金属状态方程： = (#--po)+(y--1)pe，其中， ,o0=8．93g／cm ，Co=0．394c ，2=2．685，PH(0)=0．168lxl0MPa。粒子初始呈轴对称三角形交叉分布：切 向9Oo分240个粒子，径向分布 9层粒子。 轴、Y轴取为固壁，固壁外切向设对称的 lO个虚粒子。由于 铜壳在向心运动过程中会发生变形，随粒子的半径变小，粒子的切向间距也变小，但是粒子的法向间距却 变大。而MLSPH2D方法要求粒子拟均匀分布，因此在内爆问题的计算过程中必须对粒子分布进行重分。 此模型采用粒子切向 “二合一”简单重分方法，即将切向方向相邻两粒子合二为一，由守恒关系计算新的 156 中国工程物理研究院科技年报 物理量。“二合一”重分时刻为T=4，9，14．2，17gs。计算到T=19．5 时，接近反弹时刻，此时铜壳内压 剧烈升高，内界面向外加速，计算程序出现异常。结果表明：直到 l8．5 ，计算图像仍保持轴对称性； 铜壳内界面半径和径向速度随时间的变化曲线与二维有限元方法(FEM)的计算结果符合较好。 MLSPH 方法要求计算过程中粒子呈现拟均匀分布，以保证合理的形函数的计算，应用于二维柱对称 内爆压缩过程的试算，