Processing math: 0%
主管部门: 中国航天科技集团有限公司
主办单位: 中国航天空气动力技术研究院
中国宇航学会
中国宇航出版有限责任公司

钝锥再入飞行器黏性干扰效应建模

王沛, 李锦, 江定武, 毛枚良

王沛, 李锦, 江定武, 等. 钝锥再入飞行器黏性干扰效应建模[J]. 气体物理, 2022, 7(5): 29-37. DOI: 10.19527/j.cnki.2096-1642.0957
引用本文: 王沛, 李锦, 江定武, 等. 钝锥再入飞行器黏性干扰效应建模[J]. 气体物理, 2022, 7(5): 29-37. DOI: 10.19527/j.cnki.2096-1642.0957
WANG Pei, LI Jin, JIANG Ding-wu, et al. Viscous Interaction Model for Blunt Cone Reentry Vehicle[J]. PHYSICS OF GASES, 2022, 7(5): 29-37. DOI: 10.19527/j.cnki.2096-1642.0957
Citation: WANG Pei, LI Jin, JIANG Ding-wu, et al. Viscous Interaction Model for Blunt Cone Reentry Vehicle[J]. PHYSICS OF GASES, 2022, 7(5): 29-37. DOI: 10.19527/j.cnki.2096-1642.0957

钝锥再入飞行器黏性干扰效应建模

基金项目: 

国家自然科学基金 11972362

国家数值风洞项目 NNW

详细信息
    作者简介:

    王沛(1994-) 女, 硕士, 助工, 研究高超声速计算空气动力学、软件开发。E-mail: wangp_mail@foxmail.com

    通讯作者:

    江定武(1981-) 男, 副研究员, 研究高超声速计算空气动力学、软件开发。E-mail: dwjiang2010@foxmail.com

  • 中图分类号: V211.3

Viscous Interaction Model for Blunt Cone Reentry Vehicle

  • 摘要: 针对典型的钝锥外形, 采用统一气体动理学格式(UGKS)模拟了高度70~110 km下不同Mach数和攻角的流场, 进行了流场特性的分析, 并基于黏性干扰的理论成果, 将气动力特性与第3黏性干扰参数、攻角和Mach数等参数进行关联, 建立了气动力系数的黏性干扰模型, 给出了模型预测结果的相关性分析和准确性评估。经初步测试, 该模型预测结果与UGKS直接模拟结果具有良好的一致性, 对工程应用快速获取高空气动特性具有重要意义。
    Abstract: The flow around a blunt cone with different Mach numbers and angles of attack at 70~110 km was simulated by the unified gas kinetic scheme(UGKS). The aerodynamic characteristics were analyzed. Based on theoretical achievements of viscous interaction, the viscous interaction model of aerodynamic force coefficient was established, which associates aerodynamic characteristics with viscous interaction parameter, angle of attack and Mach number. The correction analysis and accuracy assessment of the prediction results were given. Preliminary tests show that the model prediction results are in good agreement with the UGKS results, which is significant for quickly obtaining the aerodynamic characteristics at high altitudes.
  • 对高超声速飞行器而言, 随着飞行高度增加, Mach数增加, Reynolds数减少, 表面边界层厚度增加, 形成位移效应, 导致飞行器有效几何外形发生变化, 影响边界层外部的无黏流场, 而外部流场的变化也会影响壁面附近的流场。边界层与外部压缩无黏流场之间这种强耦合作用形成黏性干扰效应。这种效应会引起当地压力、摩擦阻力以及热传导等参量的剧烈变化, 进而影响飞行器的整体气动性能, 是高超声速飞行时不可忽视的复杂物理效应。

    早期的黏性干扰研究主要采用理论分析的方法, 基于高超声速边界层相似律和高超声速无黏流的压力关联公式[1]: Bertram[2]结合摩阻与压力梯度的近似关系式, 研究了零厚度带攻角平板和三角翼的黏性干扰效应对压力、摩阻和气动力的影响; Whitfield等[3]完成了冷壁细长钝锥的黏性干扰效应理论及与实验的对比研究; Stollery[4]针对指数律外形给出了凹陷体、压缩拐角和膨胀拐角等若干实例的黏性干扰流场结果; Probstein等[5]评估了钝锥外形上横向曲率在弱黏性干扰区的影响; Yasuhara[6]研究了3/4幂次律轴对称外形上较弱至中等强度黏性干扰区中的相似流动。但是早期这些研究工作基本上只针对简单外形, 复杂外形面临更加复杂的流动机理, 这些黏性理论可能失效, 难以得到正确的结果。

    20世纪70年代, 美国对航天飞机OV102进行了一系列飞行实验, 获得了大量的气动力飞行数据[7], 给出了与气动特性参量相关的第1、第2、第3黏性干扰参数, 并得到第3黏性干扰参数数据关联结果更为准确的结论[8]。基于此, 国内学者在黏性干扰效应方面取得了一些研究成果。庄逢甘等[9]针对OV102外形进行了初步研究, 提出了较为成熟的超声速飞行器黏性干扰效应的工程计算方法; 龚安龙等[10]对轨道器再入阶段的黏性干扰效应进行了研究, 考察了黏性干扰相关参数的关联特性; 毛枚良等[11]对升力体简化外形的高超声速绕流流动进行了模拟, 认为黏性干扰导致的气动力系数的改变量与黏性干扰参数在特定范围内具有良好的线性关系; 李维东等[12]结合参考温度方法提出了一种能够考虑黏性干扰效应的高超声速乘波体气动性能的工程预测方法; 陈坚强等[13]针对OV102外形以及高升阻比复杂外形黏性干扰效应, 提出基于Mach数和第3黏性干扰参数的关联参数和公式, 具有实现高超声速气动力数据天地相关的能力; 张益荣等[14]针对典型高超声速翼身组合体外形, 建立了完全气体条件下纵向气动力系数的黏性干扰模型, 并进行了不确定度量化分析; 龚安龙等[15]针对高升阻比面对称飞行器, 建立了黏性干扰增量与黏性干扰参数之间的线性关系式, 并采用Mach数效应、真实气体效应等进行修正建立了天地换算公式; 罗长童等[16]利用专业化机器学习算法, 得到了不同风洞试验结果共同遵守的不变规律, 为建立黏性干扰等物理效应天地换算方法提供了新思路。

    国内外的研究现状表明, 传统黏性干扰的研究主要基于连续性假设成立的N-S方程, 但是真实飞行条件下随着速度和高度进一步变大, 气体分子远离平衡态, 稀薄效应变得不可忽略, 导致飞行器的气动特性相比于连续流有了较大差异, 而这一区域内目前仍然缺少合适的黏性干扰效应预测模型, 这种模型的建立对求解器提出了新要求。统一气体动理论(unified gas kinetic scheme, UGKS)是近些年来比较热门的一种跨流域多尺度求解方法, 在连续流区域可以得到精确的N-S解, 在稀薄区域也可以得到与DSMC一致的结果, 已经广泛应用于从低速到高速、从连续流到稀薄流的全速域全流域的流场模拟中[16]

    基于此, 本文采用UGKS求解器, 以典型外形钝锥为研究对象。预测其在高空(70~110 km)、高Mach数(≥10)条件下的气动特性, 并借助黏性干扰效应的理论成果, 将Mach数、攻角、黏性系数等参数与UGKS有黏解和Euler无黏解的差量进行关联, 建立了气动特性增量预测模型。这样后续进行流场预测时, 仅需要计算Euler无黏解, 就可通过该增量模型快速高效得到相应黏性条件下的气动力参数, 这极大简化了气动特性的求解过程。最后选取相关性曲线与Pearson乘积矩相关系数进行了数据相关性分析, 采用部分新状态验证了模型预测结果能够满足工程应用的精度要求。

    控制方程选用一般曲线坐标系下的三维可压缩Euler方程, 基于有限体积法进行离散, 采用隐式LUSGS方法进行迭代求解, 具体计算过程参见文献[17]。

    控制方程采用Shakhov模型方程, 其具体形式为

    f_t+\boldsymbol{u} \cdot \nabla f=\frac{f^{+}-f}{\tau}

    其中, f是气体在空间x, 速度u, 以及时间t的分布函数。τ为碰撞时间, 与黏性和压力有关, f+可以表示为

    f^{+}=g+g(1-{Pr}) \frac{8 \lambda^2}{5} \frac{\boldsymbol{c} \cdot \boldsymbol{q}}{\rho}\left(\lambda \boldsymbol{c}^2-\frac{5}{2}\right)

    其中, 第1项g是Maxwell分布函数, 第2项是基于原始BGK模型方程的修正项, 目的是为了得到合理的Prandtl数, cq分别为分子热运动速度矢量和热流矢量。

    方程有如下积分解

    \begin{array}{l} f(\boldsymbol{x}, \boldsymbol{u}, t)=\frac{1}{\tau} \int_0^t f^{+}\left[x-\boldsymbol{u}\left(t-t^{\prime}\right), t^{\prime}, \boldsymbol{u}\right] \mathrm{e}^{-\left(t-t^{\prime}\right) / \tau} \mathrm{d} t^{\prime}+\\ \;\;\;\;\;\;\;\;\;\;\;\;\mathrm{e}^{-t / \tau} f(x-\boldsymbol{u} t, 0, \boldsymbol{u}) \end{array}

    该积分解由初始分布函数和平衡态分布函数构成, 通过松弛时间τ将分子的碰撞和自由输运过程耦合起来, 从而实现分布函数在稀薄流和连续流区域的自动调节, 其具体表达式见文献[18]。

    宏观守恒量Q以及应力P、热流q等变量均可以通过在速度空间上对f求矩得到

    \begin{array}{c} \mathit{\boldsymbol{Q}} = \int f \mathit{\boldsymbol{\psi }}{\rm{d}}\mathit{\Xi }, \;\;\;\;\mathit{\boldsymbol{\psi }} = \left( {1, {\boldsymbol{u}}, \frac{1}{2}{\boldsymbol{u}} \cdot {\boldsymbol{u}}} \right)\\ {\boldsymbol{P}} = \int {\boldsymbol{c}} \cdot {\boldsymbol{c}}f{\rm{d}}\mathit{\Xi }, \;\;\;\;{\boldsymbol{q}} = \int {\frac{1}{2}} {\boldsymbol{c}} \cdot {\boldsymbol{c}} \cdot {\boldsymbol{c}}f{\rm{d}}\mathit{\Xi } \end{array}

    其中, dΞ=dudvdw为相空间的体积元。

    基于上述模型构造了一套隐式三维大规模并行软件, 并进行了多个算例的验证, 详见文献[19-20]。

    相较于返回舱, 钝锥适用于需要大横向机动能力、良好的机动性、低着陆误差和低减速负载的任务以及高进入速度和稀薄大气层的任务。本文采用钝锥外形进行高超声速黏性干扰模型的研究。钝锥球头半径取为600 mm, 半锥角10°, 全长L=3 600 mm, 力系数参考面积为1.0 m2, 力矩参考点位于质心, 坐标为(2 418 mm, 0 mm, 0 mm)。图 1给出了钝锥的计算网格示意图, 表 1给出了计算状态。

    图  1  计算网格示意图
    Figure  1.  Computational grid
    表  1  计算状态
    Table  1.  Computation states
    h/km Ma α/(°)
    70, 80, 85, 90, 100, 110 10, 15 10, 20, 30
    下载: 导出CSV 
    | 显示表格

    图 2给出了Mach数10、攻角10°时有黏(100 km)以及无黏两个状态的流场压力云图和速度矢量图, 图中变量均为无量纲结果。100 km时, 黏性边界层很厚, 黏性干扰效应明显改变了物面压力和空间流场结构。

    图  2  速度矢量及压力云图
    Figure  2.  Velocity vector and pressure contour

    图 3给出了两个高度下对称面及物面局部Knudsen数分布。图中局部Knudsen数变化范围为0.01~100, 而且随着高度进一步增加, 激波内部、壁面附近的局部Knudsen数进一步增加, 稀薄效应更加明显。这些局部区域内传统基于连续性假设的N-S方程已不再适用, 模拟的流场已不再准确。为了更好地捕捉流场细节, 预测真实流场, 需要采用跨流域求解器。

    图  3  对称面及物面局部Knudsen数分布
    Figure  3.  Local Kn number on symmetry plane and wall

    图 4给出了Mach数为10、攻角10°时不同高度对称面固壁压力分布比较, 在球头与锥的交接位置(x/L=0.167)存在明显的压力分布拐点。迎风面大面积区域(锥身)压力分布随着高度的增加呈增加趋势。背风面大面积区域压力分布随高度增加呈现先增后减的趋势, 但是量值较迎风面大约要小1个量级, 这与物理实际基本一致。

    图  4  对称面压力比较
    Figure  4.  Comparison of pressure on symmetry plane

    图 5, 6分别给出了Mach数为10、攻角为10°和30°时壁面压力和摩擦阻力对轴向力系数的贡献随黏性干扰参数的变化。随着黏性干扰参数的增加, 轴向力系数增加, 压力、黏性的贡献也在增加, 而且相比压力项, 黏性项的作用占主导, 且对比发现攻角变化对轴向力的增量影响较小。

    图  5  攻角10°时轴向力系数特性
    Figure  5.  Axial force coefficient characteristics at α=10°
    图  6  攻角30°时轴向力系数特性
    Figure  6.  Axial force coefficient characteristics at α=30°

    图 7给出了部分状态的俯仰力矩特性。高度100 km(V~0.33)及以下, 俯仰力矩变化缓慢, 符号为正, 压心在质心之前, 飞行器处于静不稳定状态。高度进一步增加到110 km, 攻角为10°,20°情况下, 俯仰力矩符号变为负数, 压心后移到质心之后, 飞行器处于静稳定状态。且Mach数对俯仰力矩的影响较小。

    图  7  俯仰力矩系数特性
    Figure  7.  Pitching moment coefficient characteristics

    以往的研究成果表明, 第3黏性参数能够较好地将来流条件与气动特性参量关联起来[13], 所以本文基于第3黏性干扰参数建立黏性干扰预测模型。该参数表达式为

    \bar{\nu}_{\infty}^{\prime}=\frac{M a_{\infty} \sqrt{C^{\prime}}}{\sqrt{R e_{L \infty}}}, C^{\prime}=\frac{\mu^{\prime} T_{\infty}}{\mu_{\infty} T^{\prime}}

    其中, T′表示边界层内的参考温度。从式中可以看出, Mach数一定的情况下, 高度增加, Reynolds数减小, 黏性干扰参数增加。为了更直观地表征黏性干扰参数与气动参量之间的关系, 所有的对比图均采用黏性干扰参数为横坐标。

    图 8给出了10°攻角对称面上几个典型流向位置处的压力改变量的变化情况。迎风面, 高度低于110 km时, 压力改变量与黏性干扰参数线性关系比较好。背风面, 基本看不到较好的线性关系。这说明, 即便是钝锥这种简单的外形, 压力改变量与黏性干扰参数之间也不能呈现线性关系。

    图  8  压力增量随黏性干扰参数的变化
    Figure  8.  Variation of pressure increment with the viscous interaction parameter

    以轴向力和俯仰力矩系数为例, 给出黏性干扰模型的分析建立过程。可以类似地得到法向力系数预测模型。

    图 9给出了轴向力和俯仰力矩系数变化量随第3黏性干扰参数的变化, 这里的系数变化量表示UGKS求解器得到的有黏解与Euler方程求解得到的无黏解的差量。

    图  9  力系数增量随黏性干扰参数的变化
    Figure  9.  Variation of force coefficient increment with the viscous interaction parameter

    从图中可以看到, 高度小于100 km时, 轴向力系数变化量与第3黏性干扰参数呈较好的线性关系, 但是随着高度的增加, 则逐渐偏离线性关系。因而, 本文采用二次多项式进行拟合, 并参照黏性干扰模型的推导过程[13], 初步给出了力系数的黏性干扰模型

    \begin{aligned} \Delta C\left(\bar{\nu}_{\infty}^{\prime}, M a_{\infty}, \alpha\right) \approx & a\left(M a_{\infty}, \alpha\right)+b\left(M a_{\infty}, \alpha\right) \bar{\nu}_{\infty}^{\prime}+\\ & c\left(M a_{\infty}, \alpha\right) \bar{\nu}_{\infty}^{\prime 2} \end{aligned}

    并进一步假设

    \begin{aligned} &a\left(M a_{\infty}, \alpha\right) \approx a_0+a_1\left(M a_{\infty}\right)+a_2(\alpha) \\ &b\left(M a_{\infty}, \alpha\right) \approx b_0+b_1\left(M a_{\infty}\right)+b_2(\alpha) \\ &c\left(M a_{\infty}, \alpha\right) \approx c_0+c_1\left(M a_{\infty}\right)+c_2(\alpha) \end{aligned}

    以往研究还发现[17], 气动力系数和黏性干扰效应引起的气动力系数增量均与Mach数平方的倒数项呈较好的线性关系, 因此模型选取1/Ma2项作为自由参数。根据钝锥外形的数值模拟结果, 通过数据拟合可以得到轴向力系数和俯仰力矩系数增量的黏性干扰预测模型

    \begin{gathered} \Delta C_{\mathrm{A}}\left(\bar{\nu}_{\infty}^{\prime}, M a_{\infty}, \alpha\right) \approx a^{\mathrm{A}}\left(M a_{\infty}, \alpha\right)+b^{\mathrm{A}}\left(M a_{\infty}, \alpha\right) \cdot \bar{\nu}_{\infty}^{\prime}+c^{\mathrm{A}}\left(M a_{\infty}, \alpha\right) \cdot \bar{\nu}_{\infty}^{\prime 2} \\ a^{\mathrm{A}}=-1.26 \times 10^{-2}-\frac{3.21 \times 10^{-3}}{M a_{\infty}^2}+6.14 \times 10^{-4} \cdot \alpha-5.38 \times 10^{-6} \cdot \alpha^2 \\ b^{\mathrm{A}}=1.31+\frac{2.72 \times 10^1}{M a_{\infty}^2}-2.23 \times 10^{-4} \cdot \alpha-7.08 \times 10^{-5} \cdot \alpha^2 \\ c^{\mathrm{A}}=-0.787-\frac{2.81 \times 10^1}{M a_{\infty}^2}+1.69 \times 10^{-2} \cdot \alpha-1.36 \times 10^{-4} \cdot \alpha^2 \\ \Delta C_{\mathrm{m}}\left(\bar{\nu}_{\infty}^{\prime}, M a_{\infty}, \alpha\right) \approx a^{\mathrm{m}}\left(M a_{\infty}, \alpha\right)+b^{\mathrm{m}}\left(M a_{\infty}, \alpha\right) \cdot \bar{\nu}_{\infty}^{\prime}+c^{\mathrm{m}}\left(M a_{\infty}, \alpha\right) \cdot \bar{\nu}_{\infty}^{\prime 2} \\ a^{\mathrm{m}}=5.22 \times 10^{-3}-\frac{2.22 \times 10^{-1}}{M a_{\infty}^2}-5.69 \times 10^{-4} \cdot \alpha-8.73 \times 10^{-6} \cdot \alpha^2 \\ b^{\mathrm{m}}=-2.66 \times 10^{-1}+\frac{8.49}{M a_{\infty}^2}+2.12 \times 10^{-2} \cdot \alpha-3.42 \times 10^{-4} \cdot \alpha^2 \\ c^{\mathrm{m}}=2.03 \times 10^{-1}-\frac{10.36}{M a_{\infty}^2}-3.76 \times 10^{-2} \cdot \alpha+7.54 \times 10^{-4} \cdot \alpha^2 \end{gathered}

    采用45°相关线对结果进行初步分析。图 10给出了黏性干扰模型预测数据(以model表示)与数值模拟计算数据(以UGKS表示)的相关性曲线, 不同攻角和Mach数条件下数据基本分布在相关线附近, 可见数据间相关性较好。

    图  10  黏性干扰模型与数值模拟预测结果相关性
    Figure  10.  Correlation between viscous interaction model and numerical simulation results

    根据上述相关性曲线, 采用相对正交距离考察黏性干扰模型数据拟合精准度

    \mathrm{d} r_i=\frac{d_i}{x_r}=\frac{\left|x_i-f\left(x_i\right)\right|}{\sqrt{2} x_r}

    其中, di为数据点到相关性曲线的正交距离, xr为数据点在相关性曲线上投影对应的横坐标, 示意图见图 11。图中横纵坐标分别为干扰模型与UGKS模拟得到的力/力矩系数增量。轴向力和俯仰力矩系数的相对正交距离计算结果如图 12所示。对于轴向力系数, 高度大于70 km时, 模型预测的力系数变化量的相对偏差都低于3%; 70 km时, 相对偏差最大可达14%, 但是此时变化量的量值(约为0.236)仅为基准量(Euler无黏解, 约为1.18)的1/5。基准量叠加变化量得到的轴向力系数预测值与UGKS计算值偏差大约仅为3%。俯仰力矩系数变化量的相对正交距离具有类似的情况, 在变化量较小的时候, 相对正交距离甚至可达28%。但是此时的变化量量值比基准量量值(Euler无黏解)低1~2个量级。

    图  11  正交距离示意图
    Figure  11.  Illustration of orthogonal distance
    图  12  力系数相对正交距离
    Figure  12.  Relative orthogonal distance of the force coefficient

    采用统计学中Pearson乘积矩相关系数进行进一步的相关性分析

    r = \frac{{\sum\limits_{i = 0}^n {\left( {{x_i} - \bar x} \right)} \left( {{y_i} - \bar y} \right)}}{{\sqrt {\sum\limits_{i = 0}^n {{{\left( {{x_i} - \bar x} \right)}^2}} \cdot \sum\limits_{i = 0}^n {{{\left( {{y_i} - \bar y} \right)}^2}} } }}

    其中, 与分别为变量xy的平均值, 系数r表征两个变量的相关程度, 其绝对值越接近1, 数据的相关性越大。计算可得轴向力和俯仰力矩模型预测结果与数值模拟结果的相关系数分别为0.999 930和0.999 841, 结果非常接近于1, 说明预测数据与计算数据的线性相关性很强。

    最后, 为了初步评估模型的准确性, 采用UGKS求解器和黏性干扰模型预测了Mach数20、攻角10°, 高度分别为100 km和110 km时的气动力增量, 具体结果及相对偏差见表 2, 本文建立的模型预测的轴向力增量与UGKS求解器模拟结果偏差小于1%, 俯仰力矩增量的预测误差在10%以内。初步说明模型预测的气动增量与UGKS黏性求解器的结果基本一致, 采用黏性干扰模型可以较好地实现不同Mach数下的数据相关, 可以为工程实践提供一定的依据。

    表  2  模型预测与UGKS模拟结果对比
    Table  2.  Comparison between model prediction and UGKS simulation
    case H/km UGKS simulation model prediction relative error
    dCA dCm dCA dCm dCA dCm
    1 100 3.116 -0.069 3.127 -0.063 0.36% -8.43%
    2 110 5.227 -0.241 5.261 -0.250 0.65% 3.52%
    下载: 导出CSV 
    | 显示表格

    本文基于统一气体动理学算法UGKS, 针对高超声速流场下的钝锥外形, 采用UGKS求解器得到黏性解, 采用Euler求解器得到无黏解, 两者相减得到气动特性增量, 然后通过第3黏性干扰参数将不同高度、Mach数和攻角下的气动特性增量进行关联, 建立了气动力系数的黏性干扰模型, 初步得到以下结论:

    (1) 在全流域范围内, 气动力系数增量与第3黏性干扰参数之间不能以简单的线性关系进行描述;

    (2) 基于UGKS求解器建立的黏性干扰模型, 在高度较高时考虑了稀薄气体效应的影响, 相比N-S求解器预测结果更加真实准确;

    (3) 利用第3黏性干扰参数和来流参数, 结合Euler无黏解以及全流域的黏性解, 可以建立具有良好精度的黏性干扰模型, 有助于高效获得飞行器全流域气动力特性。

    本文研究的钝锥飞行器外形比较简单, 与真实飞行器有一定差别。把本文的建模方法扩展至外形更加复杂的飞行器上是下一步努力的方向。

    致谢: 审稿人提出的宝贵意见对提高文章质量及完整性起到了很大的帮助和促进作用, 在此表示感谢。
  • 图  1   计算网格示意图

    Figure  1.   Computational grid

    图  2   速度矢量及压力云图

    Figure  2.   Velocity vector and pressure contour

    图  3   对称面及物面局部Knudsen数分布

    Figure  3.   Local Kn number on symmetry plane and wall

    图  4   对称面压力比较

    Figure  4.   Comparison of pressure on symmetry plane

    图  5   攻角10°时轴向力系数特性

    Figure  5.   Axial force coefficient characteristics at α=10°

    图  6   攻角30°时轴向力系数特性

    Figure  6.   Axial force coefficient characteristics at α=30°

    图  7   俯仰力矩系数特性

    Figure  7.   Pitching moment coefficient characteristics

    图  8   压力增量随黏性干扰参数的变化

    Figure  8.   Variation of pressure increment with the viscous interaction parameter

    图  9   力系数增量随黏性干扰参数的变化

    Figure  9.   Variation of force coefficient increment with the viscous interaction parameter

    图  10   黏性干扰模型与数值模拟预测结果相关性

    Figure  10.   Correlation between viscous interaction model and numerical simulation results

    图  11   正交距离示意图

    Figure  11.   Illustration of orthogonal distance

    图  12   力系数相对正交距离

    Figure  12.   Relative orthogonal distance of the force coefficient

    表  1   计算状态

    Table  1   Computation states

    h/km Ma α/(°)
    70, 80, 85, 90, 100, 110 10, 15 10, 20, 30
    下载: 导出CSV

    表  2   模型预测与UGKS模拟结果对比

    Table  2   Comparison between model prediction and UGKS simulation

    case H/km UGKS simulation model prediction relative error
    dCA dCm dCA dCm dCA dCm
    1 100 3.116 -0.069 3.127 -0.063 0.36% -8.43%
    2 110 5.227 -0.241 5.261 -0.250 0.65% 3.52%
    下载: 导出CSV
  • [1]

    Bertram M H. An approximate method for determining the displacement effects and viscous drag of laminar boundary layers in two-dimensional hypersonic flow[R]. NACA TN-2773, 1952.

    [2]

    Bertram M H. Hypersonic laminar viscous interaction effects on the aerodynamics of two-dimensional wedge and triangular planform wings[R]. NASA TND-3523, 1966.

    [3]

    Whitfield J, Griffith B. Hypersonic viscous drag effects on blunt slender cones[R]. AIAA 63-434, 1963.

    [4]

    Stollery J L. Hypersonic viscous interaction on curved surfaces[R]. AIAA 70-782, 1970.

    [5]

    Probstein R F, Elliott D. The transverse curvature effect in compressible axially symmetric laminar boundary-layer flow[J]. Journal of Aeronautical Science, 1956, 23(3): 208-224. DOI: 10.2514/8.3537

    [6]

    Yasuhara M. Axisymmetric viscous flow past very slender bodies of revolution[J]. Journal of Aeronautical Science, 1962, 29(6): 667-688.

    [7]

    Maus J R, Griffith B, Tolbert D, et al. Understanding space shuttle flight data by use of wind tunnel and CFD results[R]. AIAA 83-2745, 1983.

    [8]

    Wilhite A W, Arrington J P, McCandless R S. Performance aerodynamics of aero-assisted orbital transfer vehi-cles[R]. AIAA 84-0406, 1984.

    [9] 庄逢甘, 赵梦熊. 航天飞机的黏性干扰效应[J]. 气动实验与测量控制, 1988, 2(2): 1-11.
    [10] 龚安龙, 周伟江, 纪楚群, 等. 高超声速黏性干扰效应相关性研究[J]. 宇航学报, 2008, 29(6): 1706-1710.

    Gong A L, Zhou W J, Ji C Q, et al. Study on correlation of hypersonic viscous interaction[J]. Journal of Astronautics, 2008, 29(6): 1706-1710(in Chinese).

    [11] 毛枚良, 万钊, 陈亮中, 等. 高超声速流动黏性干扰效应研究[J]. 空气动力学学报, 2013, 31(2): 137-143.

    Mao M L, Wan Z, Chen L Z, et al. Studies of hypersonic viscous interaction[J]. Acta Aerodynamica Sinica, 2013, 31(2): 137-143(in Chinese).

    [12] 李维东, 韩汉桥, 陈文龙, 等. 考虑高空黏性干扰效应的乘波体气动性能工程预测方法研究[J]. 宇航学报, 2011, 32(6): 1217-1223. DOI: 10.3873/j.issn.1000-1328.2011.06.002

    Li W D, Han H Q, Chen W L, et al. An engineering prediction method for aerodynamic performance of waverider with hypersonic viscous interaction[J]. Journal of Astronautics, 2011, 32(6): 1217-1223(in Chinese). DOI: 10.3873/j.issn.1000-1328.2011.06.002

    [13] 陈坚强, 张益荣, 张毅锋, 等. 高超声速气动力数据天地相关性研究综述[J]. 空气动力学学报, 2014, 32(5): 587-599. https://www.cnki.com.cn/Article/CJFDTOTAL-KQDX201405005.htm

    Chen J Q, Zhang Y R, Zhang Y F, et al. Review of correlation analysis of aerodynamic data between flight and ground prediction for hypersonic vehicle[J]. Acta Aerodynamica Sinica, 2014, 32(5): 587-599(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-KQDX201405005.htm

    [14] 张益荣, 张毅锋, 解静, 等. 典型高超声速翼身组合体黏性干扰效应模型研究[J]. 空气动力学学报, 2017, 35(2): 186-191. https://www.cnki.com.cn/Article/CJFDTOTAL-KQDX201702005.htm

    Zhang Y R, Zhang Y F, Xie J, et al. Study of viscous interaction effect model for typical hypersonic wing-body figuration[J]. Acta Aerodynamica Sinica, 2017, 35(2): 186-191(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-KQDX201702005.htm

    [15] 龚安龙, 解静, 刘晓文, 等. 近空间高超声速气动力数据天地换算研究[J]. 工程力学, 2017, 34(10): 229-238. https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX201710026.htm

    Gong A L, Xie J, Liu X W, et al. Study on ground-to-flight extrapolation of near space hypersonic aerodynamic data[J]. Engineering Mechanics, 2017, 34(10): 229-238(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX201710026.htm

    [16] 罗长童, 胡宗民, 刘云峰, 等. 高超声速风洞气动力/热试验数据天地相关性研究进展[J]. 实验流体力学, 2020, 34(3): 78-89. https://www.cnki.com.cn/Article/CJFDTOTAL-LTLC202003005.htm

    Luo C T, Hu Z M, Liu Y F, et al. Research progress on ground-to-flight correlation of aerodynamic force and heating data from hypersonic wind tunnels[J]. Journal of Experiments in Fluid Mechanic, 2020, 34(3): 78-89(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-LTLC202003005.htm

    [17] 毛枚良. 高超声速复杂流动数值模拟实用算法研究[D]. 绵阳: 中国空气动力研究与发展中心, 2006.

    Mao M L. Study of practical algorithm for numerical simulation of complicated hypersonic flow[D]. Mianyang: China Aerodynamics Research and Development Center, 2006(in Chinese).

    [18]

    Xu K, Huang J C. A unified gas-kinetic scheme for continuum and rarefied flows[J]. Journal of Computational Physics, 2010, 229(20): 7747-7764.

    [19]

    Jiang D W, Mao M L, Li J, et al. An implicit parallel UGKS solver for flows covering various regimes[J]. Advances in Aerodynamics, 2019, 1: 8.

    [20]

    Li J, Jiang D W, Geng X R, et al. Kinetic comparative study on aerodynamic characteristics of hypersonic reentry vehicle from near-continuous flow to free molecular flow[J]. Advances in Aerodynamics, 2021, 3: 10.

图(12)  /  表(2)
计量
  • 文章访问数:  117
  • HTML全文浏览量:  37
  • PDF下载量:  8
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-10-08
  • 修回日期:  2021-10-31

目录

/

返回文章
返回