二维移动最小二乘粒子动力学方法MLSPH无网格方法研究.pdf

二维移动最小二乘粒子动力学方法(MLSPH) 无网格方法研究· 温万治1，董师舜1，杭旭登1，付汉清2 (1．北京应用物理与计算数学研究所计算物理实验室北京；2．中国工程物理研究院北京) 蕾要：随着对数值模拟要求的提高，更精细以及特殊的物理现象要求能够模拟，无网格方法足一条很 好的途径。移动最小二乘粒子动力学方法(MLSPH)是一种无网格的粒子计算方法。本文给出了移动最 小二乘粒子动力学方法的控制方程、MLSPH守恒格式等，重点研究了MLSPH守恒格式中MLS交叉面 积向量的数值求解方法。文中对激波管问题和平面Nob．问题进行了模拟，得到了合理的结果。 关键词=无网格方法；移动最小二乘；守恒格式；交叉面积向量 1引言 随着对数值模拟要求的提高，更精细以及特殊的物理现象(如损伤、断裂、微喷射、涡流等)要 方法(SPH)，从而引起SPH方法广泛发展，并且取得了很大进展。但直到上世纪九十年代，国际上 才愈来愈重视无网格数值模拟方法研究。首先，SPH方法在爆炸力学和高速碰撞方面得到广泛应用， 它能计算大变形，但是存在精度低及计算不稳定性问题，不少研究者在应用中对此进行了改进和发展。 任意不规则网格上的有限差分方法。1 式中，带动了弱形式无网格方法的迅速发展。1999年，Los SPH方法中，显著地提高了计算精度和稳定性，取得了很好的效果，取名为移动最小二乘粒子动力学 Least Particles 方法(MLSPH：MovingSquares 迅速发展，受到国际计算力学界的高度重视。最近，LawrenceLivemore实验室的Jin Yao发表了爆轰 模拟的MLS无网格方法研究文章眺10J。MLSPH方法是一种无网格的粒子计算方法，它在计算区域内 的拟均匀分布的离散点(粒子)上，引入具有紧支域的权函数，然后求得移动最小二乘插值的形函数， 再借助形函数来求解流体力学的MLSPH方程组。本文给出了二维MLSPH方法的控制方程、移动最 的数值求解方法。文中对激波管问题和平面Noh问题进行了模拟，得到了合理的结果。 2控制方程 在连续介质、局部热动平衡、各向同性等基本假定下，考虑可压、无粘性、无热传导的流体，并 且没有外源，忽略重力，则控制方程如下： 质量守恒方程： 112 ||一 一” 和～办 肌 动量守恒方程： I|一 p． (2) 兹一出 扣 能量守恒方程： !!!二圭2兰：兰!：一土v．(p访)、‘7 (3) dt P 状态方程： P=(y一1)pe 其中：f为时间；P为流体密度；磊为流体速度矢量；P为流体压力；e为单位质量的内能。 初始与边界条件：初始给定各物质区的Po、玩、eo、Po等物理量。考虑固壁(包括对称轴) 边界、物质界面、自由面边界等。 3MLSPH守恒格式 3．1移动最小二乘近似 对于定义在闭域D上的函数“，在点，．处的邻域内，己知给定点‘∈D，J=1，2，…，M上的函数 值“，，则在点，．处，其移动最小二乘近似函数”6p)可表示为【8】： “6(r)=∑谚(r)uj (4) 其中，形函数力(，)为 谚(r)=p71(r)4_1(r)·p(弓)形(，) N 4p)=Zp(rj)prq彤p) 』士1 Wj(r)是