结构力学数值方法：边界元法(BEM)：BEM在流体力学中的应用.pdf

结构力学数值方法：边界元法(BEM)：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能够提供高效且准确的解决方案。

1.3.1原理

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

分方程，从而避免了对无限域的直接数值模拟，减少了计算量。

1.3.2内容

外部流场分析：如飞机翼型周围的流场分析，海洋结构物周围的

流体动力学分析等。

无限域问题：如声波在无限空间中的传播，地震波在地球表面的

传播等。

半无限域问题：如地下水流在半无限域中的模拟，热传导在半无

限域中的分析等。

1.3.3示例

假设我们有一个二维的流体问题，需要求解拉普拉斯方程在无限域中的解。

我们可以使用BEM来处理这个问题。

importnumpyasnp

fromscipy.specialimporthankel1

#定义格林函数

defgreen_function(r,k):

return-0.5*np.log(r)+1j*hankel1(0,k*r)

#定义边界单元

classBoundaryElement:

def__init__(self,x1,y1,x2,y2):

self.x1,self.y1=x1,y1

2

self.x2,self.y2=x2,y2

self.length=np.sqrt((x2-x1)**2+(y2-y1)**2)

self.normal=np.array([y2-y1,x1-x2])/self.length