鸭式布局火箭弹头部脱体激波与舵面干扰数值模拟研究(3)
时间:2022-08-09 22:14 来源:毕业论文 作者:毕业论文 点击:次
该方程是连续性方程的微分形式,适用于可压和不可压流动,代表流体密度,是时间,是速度矢量,代表进入连续相的质量。 2。1。2 动量守恒方程 文献综述 动量守恒是流体运动时应遵循的另一个普遍定律,描述为:在一给定的流体系统,其动量的时间变化率等于作用于其上的外力总和。可压缩粘性流体的动量守恒方程为: 上式也称为Navier-stokes 方程。Naver-stokes方程是流体力学中描述粘性牛顿流体的方程,是目前为止尚未被完全解决的方程,目前只有百分之一个特解。N-S方程式描述不可压缩流体动量守恒的运动方程。 2。1。3 能量守恒方程 将热力学第一定律应用于流体运动,即为能量方程: 2。2 模型的选择 本例中选取K-omega模型。K-omega模型是一种用于描述湍流运动的方程模型。K-omega 模型对近壁区域及尾流和绕流计算有比较高的精度,对于考虑粘性的流体优先考虑k-omega模型。方程为: 2。3空间离散格式和边界条件 本文采用AUSM 格式,AUSM格式是一种迎风数值计算格式,它的主要思想是在流场的传播中包含对流影响和声波影响。AUSM格式包括AUSM+格式和AUSMDV格式。其中AUSM+格式更适用于本文。 边界条件包括远场边界条件和壁面边界条件。本文远场边界条件为无反射边界条件,无反射边界条件是一种在计算中出现的特殊情况。在科学和工程计算中,大量出现无界或半无界区域上各种波动方程的数值求解问题,由于计算机容量有限,一般只能在有界的区域上计算,为此需引进人工边界,在这些边界上应加相应的边界条件。这些边界条件应满足两个基本要求:新的有界区域上的解应是原来的无界区域上解的近似。对波动方程来说这就要求在人工边界条件上不产生人工反射。因此这类边界条件称为无反射边界条件。本文壁面边界条件为无滑移边界条件,即壁面处法向和切向速度为零来:自[优E尔L论W文W网www.youerw.com +QQ752018766- 2。4网格 网格的合理设计和高质量的生成是计算流体力学计算的前提条件。计算网格按网格点之间的邻近关系可分为结构网格、非结构网格、混合网格。本文采用结构网格,结构网格的网格点之间的邻近关系是有序而规则的除了边界点外,内部网格点都有相同的邻近网格数,其单元是二维的四面体和三维的六面体。 并特别注意控制物面层的高度,以及控制第一层网格高度。网格要求全部右手系,且不能出现负体积,这样可以提高计算精度。 2。5算例分析 2。5。1 算例条件 算例中计算一个二维钝头体在高超声速下的流场模拟情况,网格数目为1万,边界条件分为远场和物面,远场为无反射边界条件,物面边界条件为无滑移边界条件。来流条件为马赫数为15。622、速度为2108。61m/s、温度为45。17K、 压强为23。622Pa、密度为1。814*10-3kg/m3 、粘度为2。851*10-5 Pa·s、 雷诺数为2。684*106 /m 。壁面温度为300。33K,圆柱体半径为0。0381m。 (责任编辑:qin) |