精品基于当地笛卡儿架构的无网格方法课件.pptVIP

基于当地笛卡儿架构的无网格方法 刘高联 李孝伟 （上海市应用数学和力学研究所，710072） 1 引 言 传统的数值解法在处理诸如固体裂纹扩展和流体中激波、涡面等问题时，为了避免网格的过度变形导致计算失败，必须不断地重新构造计算区域的网格，从而大大降低计算效率并造成计算资源的浪费。所以，为了有效地解决这些问题，最近几年来已有不少学者将重点放在了发展所谓的无网格方法上。 无网格方法包括以下几个步骤： 1）在分析的区域中任意地布点； 2）为每一个样点配点，形成覆盖整个计算区域的一系列子域； 3）在每个子域中构造未知量的近似函数； 4）推导样点的离散型的控制方程； 5）求解控制方程。 其中，步骤2）和3）是非常重要的，它们通常决定着操作过程的复杂程度和数值解的精确性。现有的大部分方法是预设一个半径，以样点为圆心划出一个子域，然后通过该子域内的所有点上的函数值构造近似函数。我们知道，如果点的分布是杂乱无章的，划分子域时很难预设半径大小，并且在不同子域内点的数目不一样，不容易以通用的形式构造近似函数。 本文提出一种新的无网格方法，其主要特点是自动地在每一样点建立一个笛卡儿架构并选取相应的邻近点，然后运用全导数公式构造该样点的所有导数。该方法不需要任何网格单元，所以是彻底的无网格方法，并且整个过程非常容易通用化，弥补了前有方法的不足。 2 离散方法 我们选用下面的Poisson方程来说明本文方法的离散思想： （1）  在无网格方法中，需要运用某一点的一系列邻近点的函数值来近似该点的函数值及其导数。那么一个重要任务就是为某个样点选择适当的邻点。下面将说明本文方法的选点过程。首先，在所研究的区域中任意布点。然后 选取点0作为要在其上构造导数的样点。 引入一个笛卡儿架构0xy ，其原点与点0重 合, 如图1所示。如果某一点在笛卡儿架构 0xy的第四象限中，并且该点离点0和x轴最 近，那么我们将这一点定义为点1，可以看 出点1不会是点0本身。如果某一点在笛卡儿 架构0xy的第一象限中，并且该点离点0和 x轴最近，那么我们将这一点定义为点2。如 果某一点在笛卡儿架构0xy的第一象限中， 并且该点离点0和y轴最近，那么我们将这一 点定义为点3。依据同样的原则可以定义出 图1中的点4—点8。对于所研究的区域中的 所有样点均实行同样的操作，那么为每一个 图1 样点0上的笛卡儿架构及点0的邻点 样点配置8个邻点的工作就完成了。 下面以二阶导数项 为例来说明本文方法的离散过程。 引入辅助点 和点 ，它们分别为线段 、 与x 轴的交点。那么二阶 导数项 可以写为如下的离散形式： （2A） 其中， ， , , 。 既然点 和点 为辅助点，而不是真正的节点，那么就应当消去方程(2A)中的 和 。为此，应用如下的全导数公式： （注： 为 与x轴的夹角） 可以得到 （2B1）