空气动力学数值方法：有限体积法(FVM)：高精度FVM算法.pdfVIP

空气动力学数值方法：有限体积法(FVM)：高精度FVM算

法

1空气动力学数值方法：有限体积法(FVM)：高精度FVM算

法

1.1绪论

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

有限体积法(FiniteVolumeMethod,FVM)起源于20世纪50年代，最初用于

解决流体力学中的偏微分方程。它是一种基于守恒定律的数值方法，通过将计

算域离散成一系列控制体积，然后在每个控制体积上应用守恒定律，从而将连

续的偏微分方程转化为离散的代数方程组。FVM在空气动力学、热力学、环境

工程、化学工程等多个领域有着广泛的应用，特别是在处理复杂的流体动力学

问题时，如湍流、边界层、激波等，FVM因其守恒性和稳定性而成为首选方法。

1.1.2高精度算法的重要性

在有限体积法中，高精度算法对于提高数值解的准确性和减少计算误差至

关重要。传统的低精度算法，如一阶迎风格式，虽然简单且计算效率高，但在

处理具有强梯度或激波的流场时，容易产生非物理的振荡和扩散。高精度算法，

如二阶迎风格式、WENO（WeightedEssentiallyNon-Oscillatory）格式、Roe平均

法等，能够在保持数值稳定性的同时，显著提高解的精度，减少非物理振荡，

更准确地捕捉流场中的细节特征，如激波、旋涡等。这对于空气动力学中的精

确模拟和预测，如飞机设计、风洞实验等，具有不可替代的作用。

1.2高精度有限体积法算法示例

1.2.1阶迎风格式

二阶迎风格式是一种常用的高精度算法，它通过引入上风向和下风向的差

分，提高了数值解的精度。下面是一个使用二阶迎风格式求解一维对流方程的

Python代码示例：

importnumpyasnp

#参数设置

nx=100#空间网格点数

1

nt=100#时间步数

dx=2/(nx-1)#空间步长

dt=0.025#时间步长

c=1#对流速度

#初始化网格和初始条件

u=np.ones(nx)

u[int(.5/dx):int(1/dx+1)]=2

#边界条件

u[0]=1.0

u[-1]=1.0

#二阶迎风格式求解

forninrange(nt):

un=u.copy()

foriinrange(1,nx):

ifun[i]un[i-1]:

u[i]=un[i]-c*dt/dx*(un[i]-un[i-1])

else:

u[i]=un[i]-c*dt/dx*(un[i+1]-un[i])

#输出结果

print(u)

1.2.2WENO格式

WENO格式是一种非线性高精度算法，特别适用于处理具有激波的流场。

它通过加权多个低精度格式的解，来获得一个高精度且非振荡的解。下面是一

个使用WENO格式求解一维Burgers方程的MATLAB代码示例：

%参数设置

nx=100;%空间网格点数

nt=100;%时间步数

dx=2/(nx-1);%空间步长

dt=0.025;%时间步长

c=1;%对流速度

%初始化网格和初始条件

u=ones(1,nx);

u(ceil(.5/dx):ceil(1/dx+1))=2;

%边界条件

u(1)=1.0;

u(end)=1.0;

2

%WENO格式求解

forn=1:nt

un=u;

fori=2:nx-1

%WENO权重计算

omega1=0.3;

omega2=0.6;

omega3=0.1;

%低精度格式解

uL=(un(i-2)+2*un(i-1)+un(i))/4;

uC=(un(i-1)+2*un(i)+un(i+1))/4;

uR=(un(i)+2*un(i+1)+un(i+2))/4;

%非线性权重计算

alpha1=(omega1/(eps+(uL-uC)^