计算力学无网格学习报告要素.doc

无网格数值求解学习报告 无网格法的概念 当前数值模拟欧拉方程和 N-S方程的空间离散方法主要有三种:有限差分法、有限体积法和有限元法。有限差分法的主要思想是将微分型控制方程应用到每个网格节点上，对某个节点的导数值用相关网格节点值进行差分求得。再利用各种时间推进格式进行迭代求得流场的稳定解。有限体积法的主要思想是把流体力学的积分型控制方程应用到每一个由网格剖分得到的控制体微元上，从而把积分方程的求解转化为代数方程或常微分方程的求解，再利用各种时间推进格式进行迭代，最后在每个控制体微元上得到流场的稳定解。有限元方法的数值基础是变分原理和加权余量法，其主要思想是在每个网格单元上选择若干合适的点作为解函数的插值节点，方程的近似解由各网格单元的近似函数逼近，而网格中的近似函数又由单元基函数的线性组合得到。以上三种空间离散方法都是在网格的基础上进行的。 无网格方法是一种只需要节点信息而不需将节点连接成网格单元的数值解法，这是无网格算法区别于网格算法之处。因此，无网格算法不能直接沿用网格算法中的有限差分和有限体积法等的离散格式。无网格算法的基本思想是在求解域内布置一系列的离散的点称之为“节点”，然后采用一种形函数(或称为核函数)近似拟合曲线，要求每个小域(称为“点云”)内节点的个数多于建立的方程的数目，在该点云上建立矛盾方程组，通过解矛盾方程组求得变量的梯度，进而可以求得问题的解。