空气动力学方程：简化欧拉方程的物理背景.pdfVIP

空气动力学方程：简化欧拉方程的物理背景

1空气动力学基础

1.1流体动力学基本概念

流体动力学是研究流体（液体和气体）在静止和运动状态下的行为的学科。

在空气动力学中，我们主要关注气体，尤其是空气。流体动力学的基本概念包

括：

流体的连续性：流体在流动过程中，其质量是守恒的，即流体不

能被创造或销毁。

流体的压缩性：气体的密度可以随压力和温度的变化而变化，而

液体的密度变化相对较小。

流体的粘性：流体内部层与层之间存在摩擦力，这种性质称为粘

性。

流体的涡旋性：流体在流动时可以形成涡旋，这是流体动力学中

非常重要的现象。

1.2连续性方程解析

连续性方程描述了流体质量的守恒。在三维空间中，连续性方程可以表示

为：

∂

+∇⋅=0

∂

其中，是流体的密度，是流体的速度矢量，是时间。这个方程说明了

在任意体积内，流体的质量随时间的变化率等于流体通过该体积边界流出和流

入的质量差。

1.2.1示例

22

假设我们有一个简单的二维流体流动，其中流体的速度分布为=−

和=2，密度为常数=1。我们可以计算连续性方程的左侧，以验证质量

守恒。

importnumpyasnp

fromscipy.ndimageimportgaussian_filter

#定义网格

x=np.linspace(-10,10,100)

y=np.linspace(-10,10,100)

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

1

#定义速度场

vx=X**2-Y**2

vy=2*X*Y

#计算连续性方程的左侧

rho=1#假设密度为常数

dt=0.1#时间步长

dx=x[1]-x[0]#空间步长

dy=y[1]-y[0]

#使用中心差分法计算速度的偏导数

d_vx_dx=(np.roll(vx,-1,axis=0)-np.roll(vx,1,axis=0))/(2*dx)

d_vy_dy=(np.roll(vy,-1,axis=1)-np.roll(vy,1,axis=1))/(2*dy)

#计算连续性方程的左侧

continuity_eq_left=-d_vx_dx-d_vy_dy

#对结果进行高斯平滑，以减少数值噪声

continuity_eq_left_smoothed=gaussian_filter(continuity_eq_left,sigma=1)

#输出结果

print(连续性方程左侧（平滑后）的最大值：,np.max(continuity_eq_left_smoothed))

print(连续性方程左侧（平滑后）的最小值：,np.min(continuity_eq_left_smoothed))

这个例子中，我们计算了连续性方程的左侧，并对其结果进行了高斯平滑，

以减少由中心差分法引入的数值噪声。理想情况下，连续性方程的左侧应该为

零，表示质量守恒。

1.3动量守恒方程介绍

动量守恒方程描述了流体动量随时间的变化，它由三个方程组成，分别对

应于，和方向。在没有外力作用的情况下，动量守恒方程可以表示为：

∂∂

+∇⋅=−

∂