《两相流与热物理-刘婕》第二章第一节.pptVIP

边界元方法由于它直接在边界上进行离散计算，因而特别适合于高精度边界计算的自由面和运动界面问题，其次，对一类问题，由于它仅用边界的计算代替了整个区域的计算，又具有降维计算、节约计算量与节约计算机内存的计算效益。因而，边界元方法得到了迅猛的发展。 转化的方法：（1）格林公式法；（2）变分方法；（3）加权残差方法。 其余方法可参阅C.A.Brebbia教授写的专著[9]。 原则上可以应用格林公式和推广的格林公式[10]将偏微分方程（包括椭圆型、抛物型、双曲型方程）转变为边界积分方程，对于边值问题，能转变为纯边界积分方程，充分发挥降维计算的边界元计算效益。 [9] J. Dominquez , C.A. Brebbia, “BOUNDARY ELEMENTS - An Introduction Course”， Computational Mechanics Publications, Southampton Boston, 1989. [10] 梁昆淼,“数学物理方法”，p348-359, 人民教育出版社，北京，1962。 [11] P.W. Partridge, C.A. Brebbia, L.C. Wrobel, The dual reciprocity boundary element method, Elsevier Applied Science, London, 1992. 对于非边值问题由于数理方程方程的性质决定了不能转化为纯边界积分方程，需要对积分方程的非均匀项（余留的区域积分项）进行进一步的处理。 上个世纪90年代发展起来的双倒易边界元方法能将上述积分方程的非均匀项全部转变为边界积分方程[11]，为边界元方法用于非边值问题开拓了新的途径。 然而经双倒易原理转化的边界积分方程，形式上是边界积分，实质上进行了区域节点的计算，离散后的系数矩阵是满秩的非稀疏阵，增加了计算量，为了发挥边界元计算的优点，吸引了计算数学研究工作者在这方面做更深入的研究，己经获得了一批研究成果。 边界元方法优点 2 .仅用边界计算替代整个区域计算，具有降维计算、减少计算量与节约内程等好的计算效益。 转化方法: a 格林公式方法， b 变分方法， c 加权残数方法 格林公式方法 1.奥高公式 设S为包含体积V的逐片光滑的曲面，P P x,y,z , Q Q x,y,z , R R x,y,z 及其一阶偏导数在V+S上连续，则有奥高公式： 可将奥高公式写成下述矢量形式： 其中 2.第一格林公式 我们应用奥高公式于下述面积分 得到第一格林公式 1.4.1 或 2.第二格林公式 用同样的办法我们得到 将此式与 1.4.1 相减，我们得到： 亦即： 1.4.2 4.第三格林公式 考虑二维空间，取 r 为 p, q 点间距离，p 为固定观察点，q为积分区域任一点。 在 － 中，满足下述方程： x y . q p r 4.第三格林公式 应用格林第二公式： 4.第三格林公式 当p点不在A内，而在 上，这样挖掉的是半园。 格林第二公式在挖掉 后的 内是可以用的。于是我们同样有下述推导： x y . q p r 4.第三格林公式 4.第三格林公式 p 4.第三格林公式 如果 p点不在闭区域 ＝ + 则无需作上述处理。 于是我们得到： x y . q p r 第三格林公式小结 1.4.3 1.4.2描述泰勒泡自由面问题的积分方程组 对 应用第三格林公式，令 立即得到 1.4.4 第三格林公式 1.4.3 1.4.2描述泰勒泡自由面问题的积分方程组 对 应用第三格林公式，令 立即得到 第三格林公式 将 转化为边界积分 考虑 令 代入第二格林公式得： 第二格林公式 1.4.2 流函数的边界积分方程 将 代入右式，即得： 1.4.5 1.4.3 积分方程组的离散化 离散化方法：常单元法，线性单元法，二次单元法 常单元法：用N段包含节点 的光滑直线 来细分 ，在每段线段 上取 为常值。 则积分方程组 1.4.4 1.4.5 可离散为 1.4.6 1.4.7 若记 则方程组 1.4.6 1.4.7 能分别写成下列线性常微分方程组 线性代数方程组 及其矢量形式： 1.4.8 1.4.9 1.4.10 1.4.11 1.4.4 自由面迭代 自由面迭代方程 – 法向应力平衡条件 1.2.14 对 1.2.14 式积分得： 1.4.12 1.4.13 松弛迭代格式 1.4.14 松弛因子 1.4.4 数值结果和讨论 可供比较的结果 F.P. Bretherton摄动分析结果 D.A. Reinelt P.G. Saffman FDM 结果 E.I. Shen K.