一种火星大气真实气体环境气动特性预测方法技术

本发明专利技术公开了一种火星大气真实气体环境气动特性预测方法，该方法采用化学非平衡模型计算得到火星探测器高超声速零攻角下的流场，从流场中提取出正激波位置处的温度及各气体组分的质量百分比，然后通过热力学关系式和质量百分比加权平均得到混合气体的等效比热比γeff值；将该等效比热比值γeff作为已知参数输入完全气体模型，采用CFD数值模拟方法模拟火星探测器在火星真实气体环境下的气动特性。该方法的准确性和可靠性通过美国的凤凰号火星探测器典型算例进行了验证，可准确高效地进行火星大气真实气体效应作用下的气动力性能快速预测。

A method for predicting the aerodynamic characteristics of the real gas environment in the Martian atmosphere

The invention discloses a Mars atmospheric environment of real gas aerodynamic performance prediction method, the method to calculate the Mars rover hypersonic flow under the zero angle of the chemical non-equilibrium model, extract the mass percentage of normal shock wave position and temperature of each gas component from the flow field, and then through the thermodynamic relations and mass percentage the weighted average equivalent specific heat ratio of mixed gases of eff gamma value; the ratio of equivalent specific heat gamma eff known as input parameters of perfect gas model, using CFD numerical simulation of Mars Mars probe in real gas environment simulation method for aerodynamic characteristics. The accuracy and reliability of the proposed method are verified by a typical example of Phoenix Mars probe in the United States, which can accurately and efficiently predict the aerodynamic performance of the Martian atmosphere under the effects of the real gas effect.

【技术实现步骤摘要】
一种火星大气真实气体环境气动特性预测方法
本专利技术涉及一种火星大气真实气体环境气动特性预测方法，可准确高效地进行火星大气真实气体效应作用下的气动力性能快速预测。
技术介绍
火星探测器以超高速飞行进入火星大气层，进入阶段的气动力精确预测是探测器气动布局设计的前提。火星探测器高超声速流动的主要特点是：以CO2气体为主要介质，呈现低雷诺数、高马赫数的特点，且流动常伴有真实气体效应。只有充分理解火星大气环境下的高超声速流动的特殊性，建立与之相适用的气动力预测方法，才能有效精确地预测火星探测器的气动力特性，进而确保火星探测着陆器的气动布局设计可靠。化学非平衡模型能够较好地模拟火星大气环境中的真实气体效应，但计算效率相对较低，所耗费的计算资源和计算时间较高。与化学非平衡模型相比，选择合适的比热比来等效伴有真实气体效应的非空气介质高超声速流动，也即等效比热比模型，可以作为火星探测器进入段气动力特性预测的一种简单有效的研究手段。等效比热比模型是对真实气体效应的近似，在满足工程设计精度要求的前提下，通过选取适当与准确的等效比热比γeff值，相较于采用化学非平衡气体模型的数值预测方法，能够非常有效地节省计算时间和成本。现有的等效比热比模型对于比热比值的选取，通常采用以下几种方法：(1)根据平衡流动和等效流动激波层内温度相匹配，将波后温度对应的比热比作为等效比热比；(2)根据平衡流动和等效流动激波前后密度比相匹配，通过正激波关系式得到等效比热比：(3)根据平衡流动和等效流动激波脱体距离相匹配而得到等效比热比。以上这些方法所得到的比热比等效值并不相同，且不能完全模拟出火星探测器在火星大气环境中的真实气体效应，所计算出的火星探测器气动力与化学非平衡模型的计算结果之间存在一定差异。
技术实现思路
本专利技术所要解决的技术问题是：克服现有技术的不足，提供一种既能保证计算精度，同时又能提高效率的火星大气真实气体环境气动特性预测方法。本专利技术的技术方案是：一种火星大气真实气体环境气动特性预测方法，该方法包括以下步骤：(1)、基于火星大气的化学非平衡模型，通过CFD数值模拟方法计算预设的高超声速、零攻角、零侧滑角状态下的火星探测器绕流流场参数，所述绕流流场参数包括绕流流场中各气体组分的质量百分比、流场压力P、流场温度T和流场密度ρ；(2)、根据步骤(1)所获得的预设的高超声速、零攻角、零侧滑角状态下的火星探测器绕流流场参数，找出火星探测器对称轴上气体温度T或当地密度与来流密度比ρ/ρ∞出现阶跃变化的位置，即绕流流场中正激波与探测器对称轴的交点；(3)、从绕流流场参数中提取正激波与探测器对称轴的交点处的气体温度T、各气体组分的质量百分比Qi，i＝1～N，并计算各气体组分的定压比容Cpi(T)，i＝1～N，N为火星气体组分数；(4)、通过各气体组分的质量百分比加权得到正激波与探测器对称轴交点处混合气体的定压比容Cp(T)：(5)、根据步骤(4)计算得到的正激波与探测器对称轴的交点处混合气体的定压比容Cp(T)，计算正激波与探测器对称轴的交点处混合气体的比热比γ：γ＝Cp(T)/(Cp(T)-R)，其中，R为火星大气的气体常数；(6)、以正激波与探测器对称轴的交点处的混合气体比热比γ作为该火星探测器绕流流场的等效比热比γeff，将该等效比热比γeff作为已知参数输入完全气体模型，采用CFD数值模拟方法模拟步骤(1)所述的高超声速、不同攻角、不同侧滑角的火星探测器在火星真实气体环境下的气动特性，所述气动特性包括静态的轴向力、法向力、质心俯仰力矩、以及压力分布。所述化学非平衡模型为8组分、9反应的的化学反应动力学模型，所述8组分是指：CO2，CO，O2，O，C，N2，N，NO，9反应见下表：表1火星大气化学非平衡模型化学反应类型步骤(3)采用七次多项式计算带有化学反应的火星大气各气体组分的定压比容Cpi(T)，所述七次多项式为Cpi(T)＝(a1T-2+a2T-1+a3+a4T+a5T2+a6T3+a7T4)×R式中，R为火星大气的气体常数，a1～a7为多项式系数，所述多项式系数根据气体组分和温度确定。本专利技术与现有技术相比的优点如下：(1)、本专利技术针对伴有真实气体效应的火星探测器高超声速流动，提出了一种新的等效比热比模型，建立了火星大气真实气体环境的气动特性预测方法，该方法准确高效，可为火星探测器的气动布局设计提供依据。(2)、本专利技术的等效比热比模型与现有的等效比热比模型相比，能够有效改善γeff的取值精度，从而提高整体气动力的计算精度。经美国凤凰号(Phoenix)火星探测器的典型算例证明，本专利技术的等效比热比模型与化学非平衡模型相比，轴向力系数CA最大偏差不超过1％，法向力系数CN最大偏差约为0.0025，质心俯仰力矩系数CMzg最大偏差为0.001。(3)、本专利技术的等效比热比模型与化学非平衡模型计算相比，能够显著提高求解效率，降低时间与经费成本。经比较证明本专利技术的等效比热比模型相对于化学非平衡模型，其求解效率能够提高约1个量级。附图说明图1为本专利技术火星大气真实气体环境气动特性预测方法流程图；图2(a)为采用化学非平衡模型计算得到的高超声速、零攻角、零侧滑角状态下的火星探测器绕流流场图；图2(b)为火星探测器绕流流场中正激波前后的密度变化图；图3(a)为火星探测器对称轴上正激波附近的温度分布图；图3(b)为火星探测器对称轴上正激波附近的密度比分布图。图4(a)为本专利技术与现有技术提供的轴向力系数CA比较结果；图4(b)为本专利技术与现有技术提供的法向力系数CN比较结果；图4(c)为本专利技术与现有技术提供的质心俯仰力矩系数CMzg比较结果。具体实施方式下面结合实例，说明本专利技术的具体实施方式。如图1所示，本专利技术提供了一种火星大气真实气体环境气动特性预测方法，该方法包括以下步骤：(1)、采用火星大气的化学非平衡模型，通过CFD数值模拟方法计算预设的高超声速、零攻角、零侧滑角状态下的火星探测器绕流流场参数，所述绕流流场参数包括绕流流场中各气体组分的质量百分比、流场压力P、流场温度T和流场密度ρ，所述预设的高超声速马赫数大于等于5；所述化学非平衡模型为8组分(CO2，CO，O2，O，C，N2，N，NO)、9反应的化学反应动力学模型。具体反应如下表，反应常数参见文献“Park,C.,ReviewofChemical-KineticsProblemsofFutureNASAMissions,Part2:MarsEntries,JournalofThermophysicsandHeatTransfer,8(1):9-23,1994”。表1火星大气化学非平衡模型化学反应类型(2)、根据步骤(1)所获得的预设的高超声速、零攻角、零侧滑角状态下的火星探测器绕流流场参数，找出火星探测器对称轴上气体温度T或当地流场密度与来流密度比ρ/ρ∞出现阶跃变化的位置，即为绕流流场中正激波与探测器对称轴的交点，具体为：首先，从步骤(1)获得的火星探测器绕流流场参数中提取火星探测器对称轴上的气体温度T或来流密度ρ；然后，绘制流场温度T或者当地密度与来流密度比ρ/ρ∞随火星探测器对称轴位置的变化曲线；最后，找出流场温度或者当地密度与来流密度比ρ/ρ∞出现阶跃变化的位置，即为正激波与探测器对称本文档来自技高网...

【技术保护点】
一种火星大气真实气体环境气动特性预测方法，其特征在于以下步骤：(1)、基于火星大气的化学非平衡模型，通过CFD数值模拟方法计算预设的高超声速、零攻角、零侧滑角状态下的火星探测器绕流流场参数，所述绕流流场参数包括绕流流场中各气体组分的质量百分比、流场压力P、流场温度T和流场密度ρ；(2)、根据步骤(1)所获得的预设的高超声速、零攻角、零侧滑角状态下的火星探测器绕流流场参数，找出火星探测器对称轴上气体温度T或当地密度与来流密度比ρ/ρ∞出现阶跃变化的位置，即绕流流场中正激波与探测器对称轴的交点；(3)、从绕流流场参数中提取正激波与探测器对称轴的交点处的气体温度T、各气体组分的质量百分比Qi，i＝1～N，并计算各气体组分的定压比容Cpi(T)，i＝1～N，N为火星气体组分数；(4)、通过各气体组分的质量百分比加权得到正激波与探测器对称轴交点处混合气体的定压比容Cp(T)：

【技术特征摘要】
1.一种火星大气真实气体环境气动特性预测方法，其特征在于以下步骤：(1)、基于火星大气的化学非平衡模型，通过CFD数值模拟方法计算预设的高超声速、零攻角、零侧滑角状态下的火星探测器绕流流场参数，所述绕流流场参数包括绕流流场中各气体组分的质量百分比、流场压力P、流场温度T和流场密度ρ；(2)、根据步骤(1)所获得的预设的高超声速、零攻角、零侧滑角状态下的火星探测器绕流流场参数，找出火星探测器对称轴上气体温度T或当地密度与来流密度比ρ/ρ∞出现阶跃变化的位置，即绕流流场中正激波与探测器对称轴的交点；(3)、从绕流流场参数中提取正激波与探测器对称轴的交点处的气体温度T、各气体组分的质量百分比Qi，i＝1～N，并计算各气体组分的定压比容Cpi(T)，i＝1～N，N为火星气体组分数；(4)、通过各气体组分的质量百分比加权得到正激波与探测器对称轴交点处混合气体的定压比容Cp(T)：(5)、根据步骤(4)计算得到的正激波与探测器对称轴的交点处混合气体的定压比容Cp(T)，计算正激波与探测器对称轴的交点处混合气体的比热比γ：γ＝Cp(T)/(Cp(T)-R)，其中...

【专利技术属性】
技术研发人员：詹慧玲，周伟江，刘周，龚安龙，纪楚群，
申请(专利权)人：中国航天空气动力技术研究院，
类型：发明
国别省市：北京,11