物理质团法一个普适的数值方法体系-计算物理.PDFVIP

第３４卷 第３期 计　 　 算　 　 物　 　 理 Ｖｏｌ．３４，Ｎｏ．３ ２０１７年５月 　 ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＣＯＭＰＵＴＡＴＩＯＮＡＬ ＰＨＹＳＩＣＳ　 Ｍａｙ，２０１７ 文章编号：１００１⁃２４６Ｘ（２０１７）０３⁃０２６１⁃１２ 物理质团法：一个普适的数值方法体系 马智博 （北京应用物理与计算数学研究所，北京　 １０００９４） 摘　 要：无网格方法根据分布于近邻空间各个方向的微元体物理信息构造离散方程，显著降低了空间导数计算对 微元体本身及微元体之间拓扑结构的条件限制，极大提高了拉氏方法的大变形计算能力． 由于不能利用微元体的 完备几何信息，不容易构造符合物理的无网格算法，对那些物理参数存在间断的模型对象，难以获得稳定和准确的 计算结果．本文基于对物理规律及数值模拟发展趋势的分析，提出符合物理且具有强普适性的无网格方法体系．基 于该方法的一维算例表明，即使物理参数存在强烈间断，数值结果也能很好地逼近问题的真解． 关键词：数值模拟；无网格方法；光滑粒子法；物理质团法 中图分类号：Ｏ２４２１ 文献标志码：Ａ ０　 引言 数值模拟包含建模（Ｍｏｄｅｌｉｎｇ）与模拟（Ｓｉｍｕｌａｔｉｏｎ）两大技术环节，模拟环节中数值计算的能力及计算结 ［１－３］ 果的准确性是数值模拟走向实用的决定因素 ．宏观尺度的数值模拟被广泛地用于工程系统的设计和可 靠性认证，而微观和介观尺度的数值模拟则更多地用于物质的微结构研究和宏观数值建模．用数值方法求解 非线性微分方程，是数值模拟的技术关键．对激波等物理量存在间断的宏观问题，导数不存在，一般按照Ｖｏｎ Ｎｅｕｍａｎｎ或Ｇｏｄｕｎｏｖ提出的思想构造数值方法．前者为间断捕获法，对整个计算域均按连续光滑假设实施导 ［４］ 数的近似估计，通过添加人为粘性项来捕捉物理量的突变及发生突变的位置 ．后者为间断分解法，对整个 ［５］ 计算域均按不连续假设给出Ｒｉｅｍａｎｎ解，也用统一算法来表现物理量的缓变或突变 ．Ｖｏｎ Ｎｅｕｍａｎｎ方法借 助添加的人为粘性项，通过宏观离散体和宏观时间步所对应的物理量来表现发生在微观时空下的现象．本文 将这类作法称为模拟环节中的模化（Ｍｏｄｕｌａｒｉｚｉｎｇ），将产生的添加项称为模化因子（Ｍｏｄｕｌａｒｉｚｉｎｇ Ｆａｃｔｏｒ）． 无网格方法不依赖微元体拓扑结构，能显著提高大变形计算的能力，但面对工程领域提出的高可信度要 求，仍有一些技术问题需要解决． 自１９７７ 年提出以来，基于 Ｖｏｎ Ｎｅｕｍａｎｎ 思想的纯拉氏无网格 ＳＰＨ （Ｓｍｏｏｔｈｅｄ Ｐａｒｔｉｃｌｅ Ｈｙｄｒｏｄｙｎａｍｉｃｓ）方法赢得了非常广泛的关注并经历了４０ 年的发展［６－７］ ，但该方法存在两 个不足：其一是动力学计算可能产生不稳定或非物理解，其二是在物质的自由边界附近，因微元体近邻信息 不完备而导致计算精度下降［８－１３］．这两个不足对大变形计算优势的发挥及数值模拟可信度的建立影响较 大．为了克服这些缺点，有学者对该方法的稳定性及核函数结构进行了研究与改进．Ｓｗｅｇｌｅ于１９９５年通过傅 ［１４］ 里叶分析法对稳定性进行研究，给出引起失稳的条件 ．Ｌｉｕ 于１９９５ 年提出的ＲＫＰＭ（Ｒｅｐｒｏｄｕｃｉｎｇ Ｋｅｒｎｅｌ ＰａｒｔｉｃｌｅＭｅｔｈｏｄ）方法，用以空间坐标作为自变量的多项式函数对ＳＰＨ核函数进行校正，通过对多项式基函数 核估计的相容性原则确定校正函数，该方法可有效提高计算域内或边界附近对函数的估计精度，但只是有限 提高了对导数的估计精度［１５－１６］．Ｃｈｅｎ于１９９９年提出的ＣＳＰＭ（Ｃｏｒｒｅｃｔｉｖｅ Ｓｍｏｏｔｈｅｄ Ｐａｒｔｉｃｌｅ Ｍｅｔｈｏｄ）方法，借 助ＳＰＨ核函数或其导数对任意函数及