修正并行FPM算法求解两类非线性动力学方程及其模拟研究.docxVIP

修正并行FPM算法求解两类非线性动力学方程及其模拟研究

一、绪论

1.1研究背景与意义

在科学与工程领域中，众多复杂的物理现象和系统行为常需借助非线性动力学方程进行描述。这些方程广泛应用于力学、物理学、生物学、化学、天文学等多个学科，成为研究非线性动力系统中各种运动状态定量和定性规律，特别是运动模式演化行为的关键工具。例如，在力学领域，描述结构大变形、材料非线性本构关系时会用到非线性动力学方程；在物理学中，研究量子力学中的薛定谔方程、电磁学中的麦克斯韦方程组等，它们在特定情况下也呈现出非线性特征；在生物学里，可用于刻画生物种群的动态变化、神经传导等复杂过程；在化学领域，化学反应动力学中的一些模型也涉及非线性动力学方程；天文学中，天体的运动、星系的演化等问题同样离不开它的应用。然而，大多数非线性动力学方程难以获得精确的解析解，数值求解便成为获取其近似解的重要手段，对深入理解和研究相关领域的现象与问题起着至关重要的作用。

传统的数值方法，如有限元法、有限差分法等，在处理一些简单问题时表现出一定的优势，但在面对复杂几何形状、边界条件以及高维问题时，往往存在局限性。例如，有限元法需要对求解域进行网格划分，当求解域形状复杂时，网格生成和重划分过程繁琐且耗时，并且在处理大变形问题时，网格容易发生畸变，导致计算精度下降甚至计算失败；有限差分法在处理不规则边界时，边界条件的处理较为困难，可能会引入较大的误差。随着计算机技术的不断发展，并行计算技术为提高数值计算效率提供了新的途径。并行计算能够充分利用多处理器或多核处理器的计算资源，将计算任务分解为多个子任务，同时进行计算，从而显著缩短计算时间。在这种背景下，无网格有限点集方法（FPM）应运而生。FPM作为一种新兴的数值方法，无需对求解域进行网格划分，直接利用离散的点来近似求解域，避免了网格生成和重划分的难题，在处理复杂几何形状和边界条件时具有显著优势。此外，FPM还具有较高的计算精度和灵活性，能够更好地适应各种复杂问题的求解需求。

本文提出的修正并行FPM算法，是在传统FPM算法基础上的改进与创新。通过对算法的修正，进一步提高了其计算精度和稳定性。将并行计算技术引入FPM算法中，充分发挥并行计算的优势，显著提升了算法的计算效率，使其能够更快速地求解大规模的非线性动力学方程。对于麦克斯韦方程组和非线性薛定谔方程这两类典型的非线性动力学方程，利用修正并行FPM算法进行数值模拟研究。麦克斯韦方程组描述了电场、磁场的相互作用和传播规律，在电磁学、通信工程等领域有着广泛的应用；非线性薛定谔方程则在量子力学、非线性光学等领域占据重要地位，用于描述微观粒子的波函数演化以及光在非线性介质中的传播等现象。通过对这两类方程的数值模拟，能够更深入地了解修正并行FPM算法在不同类型非线性动力学方程求解中的性能和优势，为其在实际工程和科学研究中的应用提供有力的理论支持和实践经验。同时，也为相关领域的研究人员提供了一种高效、准确的数值计算方法，有助于推动相关领域的发展和进步。

1.2研究现状

1.2.1无网格有限点集方法

无网格方法作为计算力学领域近年来备受关注的数值方法，起源于20世纪90年代。与传统依赖网格的数值方法不同，它在处理复杂几何结构问题时，能在精度和计算效率间实现更优平衡。其发展历程中，1977年Lucy和Monaghan首次提出基于拉格朗日公式的光滑质点流体动力法（SPH），这是无网格方法的重要开端，最初用于解决无边界天体物理问题。此后，众多学者对其深入研究，Swegle等指出SPH方法不稳定的原因，并提出黏度系数保证运算稳定；Dyk则提出应力粒子法改善其稳定性。目前，SPH方法已广泛应用于水下爆炸数值模拟、弹丸侵彻混凝土数值模拟、高速碰撞等材料动态响应的数值模拟等领域。我国学者张锁春对SPH方法进行了综述，贝新源等将其用于高速碰撞问题，宋顺成等用于模拟弹丸侵彻混凝土。

Nayroles首先提出移动最小二乘法（MLS）并应用于边值问题的求解，进而提出模糊单元法（DEM），MLS的提出为无网格方法的发展奠定了基础。陈美娟、程玉民等提出了改进的移动最小二乘法。无网格方法目前在完全流体力学、有限弹性问题、非线性波以及高维度问题等方面都取得了良好应用。在完全流体力学领域，因其具有处理大变形、非连续场、自由表面和微观结构等优势，已发展成为一种有效的流体力学数值方法，被应用于波动传播、空气动力学、起伏型变流动和颗粒流等各种流体问题中。

有限点集方法（FPM）作为无网格方法的一种，是基于粒子离散化的数值方法，将连续域离散化为大量粒子，通过粒子间相互作用模拟连续域内物理场。与传统有限元方法相比，无需预先定义网格，在处理复杂几何形状和边界条件时优势