空气动力学数值方法：有限差分法(FDM)：一维瞬态对流方程的有限差分解法.pdfVIP

空气动力学数值方法：有限差分法(FDM)：一维瞬态对流方

程的有限差分解法

1空气动力学数值方法：有限差分法(FDM)：一维瞬态对流

方程的有限差分解法

1.1绪论

1.1.1有限差分法的基本概念

有限差分法(FiniteDifferenceMethod,FDM)是一种广泛应用于偏微分方程数

值求解的方法，尤其在空气动力学领域中，用于模拟流体的运动和热传递过程。

FDM的基本思想是将连续的偏微分方程在空间和时间上离散化，将连续的域分

割成有限数量的网格点，然后在这些网格点上用差分近似代替导数，从而将偏

微分方程转换为一组代数方程，通过求解这些代数方程来得到原问题的近似解。

1.1.1.1示例：一维瞬态热传导方程的有限差分解

假设我们有一维瞬态热传导方程：

2

∂∂

=

2

∂∂

其中，是温度，是热扩散率。我们使用中心差分和向前差分来离散化这

个方程。

importnumpyasnp

#参数设置

alpha=0.1#热扩散率

L=1.0#域的长度

T=1.0#时间长度

dx=0.1#空间步长

dt=0.001#时间步长

nx=int(L/dx)+1#空间网格点数

nt=int(T/dt)#时间步数

#初始化温度分布

u=np.zeros(nx)

u[0]=100#左边界条件

#离散化方程

1

forninrange(nt):

u[1:nx]=u[1:nx]+alpha*dt/dx**2*(u[2:nx+1]-2*u[1:nx]+u[0:nx-1])

#输出最终温度分布

print(u)

1.1.2维瞬态对流方程的物理意义

一维瞬态对流方程描述了流体中物质或能量随时间在一维空间中的对流过

程。其数学形式为：

∂∂

+=0

∂∂

其中，是流体中某物理量的浓度或速度，是对流速度。这个方程表明，

物理量随时间的变化率与它在空间方向上的变化率成反比，且由对流速度决

定。在空气动力学中，这可以用来模拟空气中的污染物扩散、热量传递或流体

的速度分布。

1.1.2.1示例：一维瞬态对流方程的有限差分解

我们使用向前差分来离散化时间导数，中心差分来离散化空间导数，以求

解上述对流方程。

importnumpyasnp

#参数设置

c=1.0#对流速度

L=1.0#域的长度

T=1.0#时间长度

dx=0.1#空间步长

dt=0.001#时间步长

nx=int(L/dx)+1#空间网格点数

nt=int(T/dt)#时间步数

#初始化物理量分布

u=np.zeros(nx)

u[0]=1.0#左边界条件

#离散化方程

forninrange(nt):

u[1:nx]=u[1:nx]-c*dt/dx*(u[1:nx]-u[0:nx-1])

#输出最终物理量分布

print(u)

这个例子中，我们模拟了一维空间中物理量随时间的对流过程，通过调整

对流速度、空间步长和时间步长，可以观察到物理量在空间中的不同分布