空气动力学数值方法：有限体积法(FVM)：二维气体动力学方程的FVM求解.pdfVIP

空气动力学数值方法：有限体积法(FVM)：二维气体动力学

方程的FVM求解

1空气动力学数值方法：有限体积法(FVM)：二维气体动力

学方程的FVM求解

1.1绪论

1.1.1有限体积法的起源与应用

有限体积法(FiniteVolumeMethod,FVM)是一种广泛应用于流体力学、热力

学和空气动力学领域的数值方法。它的起源可以追溯到20世纪50年代，最初

是为了求解流体动力学中的偏微分方程而发展起来的。FVM的核心思想是基于

控制体的概念，将计算域划分为一系列非重叠的体积，然后在每个体积上应用

守恒定律，从而将连续的偏微分方程转化为离散的代数方程组。这种方法不仅

能够保证守恒性，还具有较好的稳定性和准确性，因此在处理复杂的流体动力

学问题时表现出色。

在空气动力学中，FVM被用于求解二维和三维的气体动力学方程，包括欧

拉方程和纳维-斯托克斯方程。这些方程描述了气体在不同条件下的运动，如速

度、压力、温度和密度的变化。通过FVM，工程师和科学家能够模拟飞机、火

箭等飞行器在不同飞行条件下的气动特性，为设计和优化提供关键数据。

1.1.2维气体动力学方程简介

二维气体动力学方程通常包括连续性方程、动量方程和能量方程，它们是

描述气体在二维空间中运动的基本方程。这些方程可以写成守恒形式，即：

∂∂∂

++=0

∂∂∂

其中，是状态向量，包含了密度、动量、和总能量；和分别是

和方向上的通量向量，包含了与速度、压力和能量相关的项。

在有限体积法中，这些方程被应用于每个控制体，通过积分形式转化为离

散方程。例如，对于一个控制体，连续性方程可以转化为：

++=0​​​

其中，是控制体的体积，是控制体的边界，和是边界上的法向量分

量。

1

1.2有限体积法求解二维气体动力学方程

1.2.1网格划分

在应用FVM求解二维气体动力学方程之前，首先需要对计算域进行网格划

分。网格可以是结构化的（如矩形网格）或非结构化的（如三角形网格）。下面

是一个使用Python和matplotlib库进行简单矩形网格划分的例子：

importnumpyasnp

importmatplotlib.pyplotasplt

#定义网格参数

nx=50#x方向网格数

ny=50#y方向网格数

Lx=1.0#x方向计算域长度

Ly=1.0#y方向计算域长度

#创建网格

x=np.linspace(0,Lx,nx+1)

y=np.linspace(0,Ly,ny+1)

X,Y=np.meshgrid(x,y)

#绘制网格

plt.figure(figsize=(8,6))

plt.plot(X,Y,k-,lw=0.5)

plt.plot(X.T,Y.T,k-,lw=0.5)

plt.xlabel(x)

plt.ylabel(y)

plt.title(2DStructuredMesh)

plt.show()

1.2.2离散化过程

离散化过程是将连续的气体动力学方程转化为离散形式的关键步骤。在

FVM中，这通常通过在每个控制体上应用积分守恒定律来实现。下面是一个离

散化连续性方程的例子：

假设我们有一个矩形控制体，其边界为−/2到+/2和−/2

到+