空气动力学方程:简化欧拉方程在飞机设计中的应用.pdf
空气动力学方程:简化欧拉方程在飞机设计中的应用
1空气动力学基础
1.1流体动力学概述
流体动力学是研究流体(液体和气体)在运动状态下的行为及其与固体边
界相互作用的学科。在飞机设计中,流体动力学尤为重要,因为它帮助工程师
理解飞机在不同飞行条件下的空气动力特性。流体动力学的核心是三大守恒定
律:质量守恒、动量守恒和能量守恒,这些定律构成了欧拉方程的基础。
1.1.1质量守恒定律
质量守恒定律,也称为连续性方程,表述为在任何封闭系统中,流体的质
量是恒定的。在流体动力学中,这意味着流体通过任意截面的流量是相等的。
对于不可压缩流体,连续性方程可以简化为:
∂
+∇⋅=0
∂
其中,是流体密度,是流体速度矢量,是时间。
1.1.2动量守恒定律
动量守恒定律描述了流体在运动过程中,其动量的变化率等于作用在流体
上的外力。在流体动力学中,这通常表示为纳维-斯托克斯方程,但在飞机设计
中,我们经常使用简化版的欧拉方程,忽略了粘性效应。欧拉方程可以表示为:
∂
⋅
+∇=−∇+
∂
其中,是流体压力,是重力加速度。
1.1.3能量守恒定律
能量守恒定律在流体动力学中表现为能量方程,它描述了流体内部能量的
变化率等于热能的产生率加上机械能的产生率。在飞机设计中,我们关注的是
流体的总能量,包括内能和动能。能量方程可以简化为:
∂
⋅
+∇=−∇⋅−⋅+⋅
∂
其中,是流体的总能量密度。
1
1.2连续性方程解析
连续性方程是流体动力学中最基本的方程之一,它确保了流体在任何点的
质量守恒。对于不可压缩流体,连续性方程简化为流体速度矢量的散度为零:
∇⋅=0
1.2.1示例
假设我们有一个二维不可压缩流体的流动,流体速度可以表示为=
,。我们可以使用Python的NumPy库来计算流体速度的散度,
并验证连续性方程。
importnumpyasnp
#定义流体速度函数
defvelocity_field(x,y):
u=x**2-y**2
v=2*x*y
returnu,v
#创建网格
x=np.linspace(0,1,100)
y=np.linspace(0,1,100)
X,Y=np.meshgrid(x,y)
#计算流体速度
u,v=velocity_field(X,Y)
#计算散度
divergence=np.gradient(u)[0]+np.gradient(v)[1]
#输出结果
print(流体速度的散度:\n,divergence)
在这个例子中,我们定义了一个流体速度场,然后计算了速度场的散度。
如果流体是不可压缩的,我们期望散度的值接近于零。
1.3动量守恒方程介绍
动量守恒方程描述了流体在运动过程中,其动量的变化率等于作用在流体
上的外力。在飞机设计中,我们通常关注的是流体对飞机表面的作用力,即升
力和阻力。简化欧拉方程忽略了粘性效应,适用于高速流动,如超音速飞行。
2
1.3.1示例
考虑一个简单的二维流体流动,我们可以使用简化欧拉方程来计算流体速
度的变化。假设流体在x方向上受到一个恒定的压力梯度,我们可以使用
Python