无网格水波动力学数值模拟理论及其应用的研究.pdf

中文摘要 N．S方程是描述流体运动最基本的方程组。用数值方法，结合一定边界条件， 对其进行求解是目前解决水动力学问题的最为有效的方法之一。以往多数研究均 采用网格方法对N—s方程进行求解。本文采用无网格方法研究N。s方程的求解 问题，并结合水波动力学的具体算例，探讨无网格数值模型在实际工程中的应用 问题。这方面的课题不仅具有很高的理论意义，也有很强的工程实用价值和应用 前景。与网格数值方法相比，这类方法的优点就在于动量方程中不直接出现对流 项，能够有效地避免数值耗散现象。 本文首先在阅读大量相关文献的基础上，对水波动力学研究中N．s方程的求 解方法及其相关湍流封闭格式进行了总结，并对其中的重要方法做了简要介绍。 本文的一个主要内容是研究运用MPS方法建立水动力学数值水槽的理论问 题。在无网格条件下，粒子间的相互作用通过核子函数来表示，并通过这种相互 作用模型建立物理量光滑模型、梯度模型和扩散模型来离散控制方程。通过引入 粒子数密度的概念实现了流体的不可压缩性。 本文的另一项主要内容是通过将理论模型在计算机上编程，建立水动力学数 值水槽。求解压力Possion方程是其中的核心工作之一。本文通过深入研究，在 计算机上成功地编制了用MICCG法(预优不完全乔里斯基分解共轭梯度法)求 解大型稀疏对称正定矩阵的计算程序，使得Poisson方程的求解速度得以大幅度 提高。 结合具体的算例分析，探讨了模型模拟水流运动的具体过程。一是成功地模 拟了水柱的倒塌过程，得出了水柱倒塌后对固体建筑物的作用力。二是研究了无 网格数值水槽中的孤立波造波问题，将数值计算结果和解析解进行了对比分析， 证明了结果的合理性。 最后，为了使计算结果更为直观，还利用Matlab软件编制了相应的动态演 示程序，实现了计算结果的可视化。 ABSTRACT N-S thebasic are fluid of equations equations isonethemost describing flowing．It effective in ofwater tosolve them approachessolvingproblems dynamics by numericalmethodunder conditions．Mostofthe research boundary proper previous workused in N-S is methods this method grid solvingequations．IIl non·鲥d paper,a onseveral of adopted．Based computingexamples，the non—grid applications numericalmodelare isnot