材料力学数值方法：边界元法(BEM)在流体力学中的应用.pdfVIP

材料力学数值方法：边界元法(BEM)在流体力学中的应用

1绪论

1.1边界元法(BEM)简介

边界元法（BoundaryElementMethod,BEM）是一种数值分析方法，主要用

于解决偏微分方程问题。与有限元法（FEM）相比，BEM主要关注于问题的边

界条件，将问题域的内部信息转化为边界上的积分方程，从而大大减少了计算

的自由度，提高了计算效率。BEM在处理无限域、半无限域以及外部流场问题

时具有独特的优势。

1.1.1原理

BEM基于格林函数（Green’sfunction）和格林定理（Green’stheorem），

将偏微分方程转化为边界积分方程。对于流体力学中的拉普拉斯方程或泊松方

程，BEM可以有效地处理其边界条件，如速度边界条件、压力边界条件等。

1.1.2内容

格林函数的构建：格林函数是BEM的核心，它描述了在边界上施

加单位源或汇时，场量在边界上的响应。

边界积分方程的形成：通过格林定理，将场量的内部信息转化为

边界上的积分表达式。

边界离散化：将连续的边界转化为离散的边界单元，每个单元上

应用边界积分方程。

数值求解：通过数值方法求解离散后的边界积分方程，得到边界

上的未知量。

1.2流体力学中的数值方法概述

流体力学中的数值方法主要包括有限差分法（FDM）、有限元法（FEM）、

边界元法（BEM）和有限体积法（FVM）。每种方法都有其适用范围和特点，

BEM在处理流体边界问题时表现出色。

1.2.1原理

流体力学数值方法的核心是将连续的流体动力学方程转化为离散形式，以

便于计算机求解。这些方程包括连续性方程、动量方程和能量方程等。

1

1.2.2内容

有限差分法：将偏微分方程转化为差分方程，适用于规则网格。

有限元法：将问题域划分为有限个单元，每个单元内假设场量的

分布形式，适用于复杂几何。

边界元法：关注边界条件，将问题转化为边界上的积分方程，适

用于无限域或外部流场。

有限体积法：基于守恒原理，将问题域划分为体积单元，适用于

非结构化网格。

1.3BEM与流体力学的结合点

BEM在流体力学中的应用主要集中在外部流场的分析，特别是对于无限域

或半无限域的流体问题，BEM可以提供高效的解决方案。此外，BEM在处理流

体与固体的相互作用问题时也表现出色，如船舶水动力学、水下结构物的流体

动力学分析等。

1.3.1原理

在流体力学中，BEM通过构建格林函数，将流体动力学方程转化为边界积

分方程，从而避免了对无限域内部的直接求解，减少了计算量。

1.3.2内容

外部流场分析：适用于无限域或半无限域的流体问题，如飞机翼

的气动分析、船舶的水动力学分析等。

流固耦合问题：处理流体与固体的相互作用，如水下结构物的流

体动力学分析、风力涡轮机叶片的气动弹性分析等。

复杂边界条件：BEM可以有效地处理复杂的边界条件，如非线性

边界条件、多连通域的边界条件等。

1.3.3示例

假设我们有一个二维不可压缩流体的外部流场问题，边界上已知速度分布，

需要求解流体内部的压力分布。我们可以使用BEM来解决这个问题。

1.3.3.1数据样例

边界上的速度分布数据：

boundary_velocities={

left:0.0,#左边界速度

right:1.0,#右边界速度

top:0.0,#上边界速度

2

bottom:0.0#下边界速度

}

1.3.3.2代码示例

构建格林函数和边界积分方程的伪代码：

#定义格林函数

defgreen_function(r):

return-1/(2*pi*r)

#定义边界积分方程

defboundary_integral_equation(boundary,boundary