en
×

分享给微信好友或者朋友圈

使用微信“扫一扫”功能。
参考文献 1
DEJARNETTEF R,HAMILTONH H.Aerodynamic heating on 3D bodies including the effects of entropy layer swallowing[J]. Journal of Spacecraft and Rockets,1975,12(1):5-12.
参考文献 2
BRANDONH J,MASEKR V,DUNAVANTJ C.Aerodynamic heating to corrugation stiffened struc-tures in thick turbulent boundary layers[J].AIAA Journal,1975,13(11):1460-1466.
参考文献 3
INGERG R.Nonequilibrium boundary-layer effects on the aerodynamic heating of hypersonic waverider vehicles[J]. Journal of Thermophysics and Heat Transfer,1995,9(4):595-604.
参考文献 4
MILOSF S,SQUIRET H.Thermostructural analysis of X-34 wing leading-edge tile thermal protection system[J]. Journal of Spacecraft and Rockets,1999,36(2):189-198.
参考文献 5
程克明,吕英伟.飞行器持续气动加热的耦合性分析[J].南京航空航天大学学报,2000,32(2):150-155.
CHENGKeming,YingweiLÜ.An analysis of coupling feature in continuing gasdynamic heating over flight vehicles[J]. Journal of Nanjing University of Aeronautics & Astronautics,2000,32(2):150-155.
参考文献 6
任青梅,杨志斌,城主,等.气动加热与结构温度场耦合分析平台研发技术[J].强度与环境,2009,36(5):33-38.
RENQingmei,YANGZhibin,CHENGZhu,et al.Development of platform for analysis coupling aeroheating and structural temperature field[J]. Structure & Environment Engineering,2009,36(5):33-38.
参考文献 7
季卫栋,王江峰,樊孝峰,等.高超声速流场与结构温度场一体化计算方法[J].航空动力学报,2016,31(1):153-160.
JIWeidong,WANGJiangfeng,FANXiaofeng,et al.Algorithms for hypersonic fluid-structural-thermal integrated[J].Journal of Aerospace Power,2016,31(1):153-160.
参考文献 8
黄唐,毛国良,姜贵庆,等.二维流场、热、结构一体化数值模拟[J].空气动力学学报,2000,18(1):115-119.
HUANGTang,MAOGuoliang,JIANGGuiqing,et al.Two dimensional coupled flow-thermal-structural numerical simulation[J].Acta Aerodynamica Sinica,2000,18(1):115-119.
参考文献 9
WIETINGA R,DECHAUMPHAIP,BEY K S,et al.Application of integrated fluid-thermal-structural analysis methods[J].Thin-Walled Structures,1991,10(1):1-12.
参考文献 10
LIOUM S.A sequel to AUSM: AUSM+[J].Journal of Computational Physics, 1996, 129(2):364-382.
参考文献 11
LEERB V. Towards the ultimate conservative difference scheme V: A second-order sequel to Godunov's method[J]. Journal of Computational Physics,1979,32(1):101-136.
参考文献 12
MENTERF R. Two-equation eddy-viscosity turbul-ence models for engineering applications[J]. AIAA Journal,1994,32(8):1598-1605.
参考文献 13
DEESEJ E, AGARWALR K. Navier-Stokes calculations of transonic viscous flow about wing/body configurations[J].Journal of Aircraft,1988,25(12):1106-1112.
参考文献 14
ZHAOY.Computation of complex turbulent flow using matrix-free implicit dual time-stepping scheme and LRN turbulence model on unstructured grids[J].Computers & Fluids,2004,33(1):119-136.
参考文献 15
WIETINGA R.Experimental study of shock wave interference heating on a cylindrical leading edge[R].NASA-TM-100484, USA: NASA, 1987.
参考文献 16
MILLERC G.Experimental and predicted heating distributions for biconics at incidence in air at Mach 10[R].NASA Technical Paper 2234,USA: NASA,1984.
目录 contents

    摘要

    针对锥体热环境问题,提出了气动热与结构传热的分区迭代推进分析方法。其中流场采用有限体积法计算,空间离散采用AUSM+格式。时间推进采用显示多步Runge-Kutta格式,结构热传导采用有限元方法求解,而数据传递采用基于虚拟空间的插值方法。圆管验证算例分析显示,2 s时刻驻点处的热流密度和温度的计算值与试验值的相对误差分别为1.34%和4.95%。最后进行了直二次圆锥体的热环境分析,壁面初始热流密度值与试验值吻合得很好,其中驻点热流的计算值与试验值的相对误差为3.1%。耦合分析过程中驻点温度随时间的推移而升高,且上升趋势逐渐变缓,最终趋于稳态值。此外时间的变化对锥体表面压强的影响可忽略不计,而壁面热流却随时间的增加而降低。

    Abstract

    The division iterative marching method on aerodynamic heating and structural heat transfer for the thermal environment of the cone body is presented. The flow field is calculated by the finite volume method. Spatial discretization scheme uses AUSM+. Explicit multi-step Runge-Kutta is used to calculate time iteration scheme. However, the structural heat transfer is calculated by the finite element method. Besides, the data exchange on the coupled wall is conducted by the interpolation method based on the virtual space. The verification example on the circular tube is analyzed, and the relative errors between the calculated values and corresponding test values for the heat flux and temperature of the stagnation point are 1.34% and 4.95% respectively at 2 s. Finally, the analysis on thermal environment of the straight biconic body is conducted. The initial wall heat flux is well matched with the experimental result, and the relative error between calculated value and experimental value at stagnation point is 3.1%. The temperature at the stagnation point increases with the time, and the upward trend slows down gradually. Finally it tends to the steady-state value. The time almost has no influence on the wall pressure, but the wall heat flux decreases with the increase of the time.

    随着高超声速飞行器的发展,其飞行速度越来越快,气动加热问题日趋严[1,2,3],它成为了高超声速飞行器设计必须考虑的因素。若能准确计算飞行器外表面的热流密度及结构的温度场,将对飞行器的安全性及性能起到极为重要的作用。传统的高超声速飞行器热问题的分析方法是将气动热分析与结构传热分开考[4],即将多物理场问题人为地分割为单物理场问题处理。其在假定飞行器壁面温度的情况下计算壁面的热流密度,再将此热流密度作用在飞行器表面以分析结构温度场。这种未考虑气动热与结构传热之间耦合效应的分析方法必然会带来计算精度的问题,无法真正评价结构所处的热环境。由于气动热和结构传热本质上具有耦合特性,必须采用多物理场耦合分析方法求解。首先气动热会对结构传热产生影响,因为外部热流密度会造成结构温度升高;此外结构温度也会对气动加热产生影响,因为结构温度的升高导致边界层内的气体与壁面之间的温度梯度减小,热流密度也将减小。

    文献[5]进行了气动热与结构传热之间的耦合分析,获得了飞行器结构温度场随时间变化的情况;任青梅[6]搭建了气动热与结构传热的耦合分析平台,计算了飞行器头锥的热流密度和结构温度场。以上研究工作中气动热分析采用工程算法,这必然会影响物面热流密度和结构温度场的分析精度。季卫栋[7]利用有限体积法开展了锥体和翼面的气动热和结构温度场的一体化分析。黄唐[8]采用有限差分法进行了流场和结构温度场的离散,并分析了锥体驻点温度随时间变化的情况。文献[9]用有限元法研究气动加热和结构传热的一体化,其在间断附近采用网格加密的方法捕捉激波。以上学者采用数值方法求解气动热,并且流场和结构温度场采用相同的数值离散方法,在耦合面需要流场和结构网格节点一一对应,这就要求结构传热模型的网格尺度必须减小到与流体网格尺度相同,这大大增加了结构传热分析的计算量。本文针对锥体热环境问题,提出气动热与结构传热的分区迭代推进分析方法,并进行了方法的圆管验证算例分析,最后进行了直二次圆锥体的热环境分析。

  • 1 数值算法

    1
  • 1.1 气动热数值算法

    1.1

    不考虑体积力和内热源情况下,在直角坐标系下的流体动力学N-S方程的积分形式为

    tΩWdΩ+Ω(Fc-Fv)dS=0
    (1)

    式中:W为守恒变量;Fc为无粘通量;Fv为粘性通量。针对以上流动控制方程,采用有限体积法进行空间离散,对流通量采用AUSM+[10]空间离散格式,该格式的数值耗散较小,激波分辨率较高,且鲁棒性较强。AUSM+格式将无粘通量分为对流项Fe和压力项P,在单元界面处其表达式为

    Fi+1/2,j,k=ρVρuV+nxpρvV+nypρwV+nzpρHVi+1/2,j,k=Fe+P
    (2)
    Fe=Mai+1/2,j,kρcρcuρcvρcwρcHL/RP=0nxpnypnxp0i+1/2,j,k
    (3)

    AUSM+格式通过单元界面处的马赫数进行流场上、下游的判断:

    ()L/R=()LMai+1/2,j,k0()RMai+1/2,j,k<0
    (4)

    式中:Mai+1/2,j,k为单元界面处的马赫数。而单元界面处的压力项则可分解为

    Pi+1/2,j,k=ψ+PL+ψ-PR
    (5)

    为获得单调解,采用完全迎风的二阶MUSCL格[11]离散分裂后的无粘通量,并采用Minmod限制器使格式达到空间二阶精度。流体控制方程中的粘性项采用中心格式进行空间离散,湍流模型采用两方程Menter’s SST k-ω[12]。针对非定常问题,采用三阶精度的显示多步Runge-Kutta方[13]求解,其在n~n+1时间步内,有

    W0=WnW1=W0+ΔtR(W0)W2=34W0+14W1+14ΔtR(W1)W3=13W0+23W2+24ΔtR(W2)Wn+1=W3
    (6)

    采用隐式残差光顺技[14]加速收敛,并且针对高超声速气动热问题,流体导热系数和粘性系数通常采用Sutherland公式或分子动力学计算,其对物面热流密度的计算精度有重要影响。由于本文的算例不考虑化学反应、热力学状态(内能松弛)等,以上控制方程和数值算法不涉及高超声速条件下的真实气体效应。

  • 1.2 结构热分析数值算法

    1.2

    在无体积热源的假设下结构瞬态热传导的控制方程为

    ρ0cTt=xkxTx+ykyTy+zkzTz
    (7)

    式中:ρ0为结构材料密度;c为材料比热容;kx,kykz分别为材料3个方向的导热系数。其中比热容和导热系数一般为温度的函数。针对本文热防护系统的热分析问题,其外表面边界条件为壁面热流密度Qaero和壁面热辐射量Qrad,其表达式分别为

    -λTnwall=Qaero(x,y,z,t)
    (8)
    Qrad=εσ(Twall4-Tat4)
    (9)

    式中:T/n为壁面法向温度梯度;ε为壁面热辐射率;σ为玻尔兹曼常数;Twall为壁面温度;Tat为大气环境温度。对式(7)进行有限元离散可得到总体合成矩阵求解方程

    CT˙n+KTn=Pn
    (10)

    式中:C为热容矩阵;K为热传导矩阵;P为温度载荷列阵;TT˙分别为节点温度列阵和节点温度对时间的导数列阵。针对第n~n+1时间步,用Gaierkin格式离散得

    2K+3ΔtCTn=2Pn+Pn-1+3ΔtC-KTn-1
    (11)

    求解式(11)即可得到结构各个时刻节点的瞬时温度,且本文结构传热采用ABAQUS有限元软件进行分析。

  • 2 耦合分析模型及策略

    2

    锥体高速飞行时产生的巨大气动热会造成结构温度升高,而结构温度的升高导致激波层和边界层内的气体与壁面之间的温度梯度减小,这将导致壁面热流密度也将减小即气动热与结构传热之间存在着强烈的耦合效应,如图1所示。

    图1
                            气动热与结构传热之间的耦合关系

    图1 气动热与结构传热之间的耦合关系

    Fig.1 Coupled relationship between aerodynamic heating and structural heat transfer

    本文采用分区迭代推进的方法进行气动热与结构传热的分析,其推进求解策略如图2所示,流场与结构分为两个不同的区域,分开建模与求解,其中流场采用有限体积法(Finite volume method,FVM)求解,而结构热传导采用有限元法(Finite element method,FEM)求解,其特点为:

    图2
                            分区迭代推进策略

    图2 分区迭代推进策略

    Fig.2 Division iterative marching strategy

    (1) 在任意n~n+1时间步内FVM(或FEM)求解过程中壁面温度(或壁面热流密度)不变,即各子学科求解时边界条件冻结。

    (2) 各子学科在相同的时间节点进行数据交换,以保证耦合分析的协调性以及耦合时间精度。

    本文气动热与结构传热的分区迭代推进方法的主要分析流程为:

    (1) 首先建立数值分析模型,定义来流条件和壁面初始温度T0,进行定常流场的计算,计算结果作为瞬态耦合分析的流场初始条件。

    (2) 进行ti时刻的非定常流场分析,将计算获得的热流密度Qi通过插值方法传递给结构传热有限元模型外表面。

    (3) 进行ti时刻的结构有限元瞬态传热分析。

    (4) 判断ti是否达到瞬态分析的总时间ttotal,若ti达到ttotal,结束耦合计算;若ti未达到ttotal,将计算获得的壁面温度通过插值方法传递给流体模型。

    (5) 返回步骤(2),进行下一个时间步ti=tit的求解,直至ti=ttotal,结束耦合分析。

  • 3 数据传递方法

    3

    耦合分析中涉及到热流和壁面温度的数据传递,而实际分析中FEM网格尺度通常远大于FVM网格,即FVM和FEM耦合面节点不一致,如图3所示,故需要在耦合面进行插值以实现数据传递。本文采用如下的基于虚拟空间的插值方法进行热流密度的数据传递:

    图3
                            FVM节点与FEM节点对应关系

    图3 FVM节点与FEM节点对应关系

    Fig.3 Nodes relationship between FVM and FEM

    (1) 将各学科耦合面上节点从物理空间(x, y, z)通过坐标变换映射到二维虚拟空间(u, v),x=x(u, v),y=y(u, v)和z=z(u, v),即将三维曲面上的空间节点转换到二维虚拟空间平面上。

    (2) 在物理空间中搜索任意FEM网格节点ζi(x, y, z)附近的FVM网格节点ηi(x, y, z),并将其转换到虚拟空间得到ζi(u, v)和ηi(u, v)。

    (3) 将FVM网格节点坐标ηi(u, v)和相应的热流Qi(u, v)代入到二次插值函数Q(u, v),采用最小二乘法求解插值函数的系数,即

    Q(u,v)=au2+bv2+cuv+du+ev+f
    (12)

    (4) 将所有FEM网格节点ζi(u, v)代入到已知系数的插值函数Q(u, v)中,即可求得FEM网格节点插值热流密度qi(u, v)。

  • 4 验证算例

    4

    采用NASA进行的高超声速圆管风洞试验模[15]验证了本文的分区迭代推进分析方法的正确性。圆管内径R1=25.4 mm,外径R2=38.1 mm,圆管材料为不锈钢,其密度ρ=8 030 kg/m3,导热系数k=16.72 W/(m·K),比热容c=502.48 J/(kg·K),此外圆管初始壁面温度T0=294.4 K。来流马赫数Ma=6.47,来流温度T=241.5 K,来流压强p=648.1 Pa,迎角α=0°。二维分析模型如图4所示,其中流场网格量为60×100,壁面第1层网格高度Δh=1×10-5 m,并且在激波处进行了局部加密,而流场分析对流项采用有限体积法AUSM+空间离散格式,湍流模型选择Menter’s SST k-ω两方程模型。圆管结构有限元网格量为40×25,采用四节点平面单元模拟。采用基于虚拟空间的插值方法进行耦合数据的传递,分析类型为瞬态分析,耦合分析时间步长Δt=1×10-4 s,分析总时间ttotal=2 s,并且将定常流场作为耦合分析的初始条件。分析获得了圆管驻点温度随时间的变化曲线,如图5所示,从图5中可观察到曲线斜率在初始阶段很大,随时间推移逐渐减小,说明结构开始阶段升温很快,随后温度上升速度逐渐减缓,最后趋于一个稳定的温度值(稳态解)。图6为2 s时刻流场和圆管结构的温度云图,从图6可观察到流场数值解捕捉到了光滑的弓形激波,且圆管最高温度出现在驻点处,沿壁面下游温度逐渐降低。

    图4
                            流体和结构模型网格

    图4 流体和结构模型网格

    Fig.4 Meshes for fluid and structural models

    图5
                            驻点温度随时间变化情况

    图5 驻点温度随时间变化情况

    Fig.5 Temperature of stagnation point varying with time

    图6
                            2 s时刻流体和结构温度场

    图6 2 s时刻流体和结构温度场

    Fig.6 Temperature fields for fluid and structural models at 2 s

    7为2 s时刻圆管外壁面的相对热流密度分布和相对温度分析情况。从图7中可知壁面耦合分析结果分布曲线与试验结果分布曲线吻合良好,其中驻点处热流密度耦合分析结果Qstag为661 kW/m2,而试验结果为670 kW/m2,相对误差为1.34%;而驻点处温度耦合分析结果Tstag为442 K,而试验结果为465 K,相对误差为4.95%。

    图7
                            2 s时刻壁面热流密度和温度分布

    图7 2 s时刻壁面热流密度和温度分布

    Fig.7 Distribution for wall heat flux and temperature at 2 s

    由以上讨论可知采用本文耦合推进分析方法得到的圆管壁面热流密度及温度分布均与试验结果吻合良好,从而验证了本文提出的分区迭代推进分析方法的正确性与分析精度,可用于分析锥体高速飞行时的热环境问题。

  • 5 锥体热环境分析

    5
  • 5.1 初始气动热分析

    5.1

    分析模型来自于NASA技术报告 NASA TP 2334[16],其几何形状为直二次圆锥体,如图8所示。来流马赫数Ma=9.86,来流静压p=59.92 Pa,来流静温T=48.88 K,迎角α=0°,壁面温度Twall=300 K。划分了CFD结构网格,如图9所示,六面体网格总量约50万,壁面第1层网格高度小于1×10-5 m,以保证量纲一参数y+≤1。进行了初始定常流场求解,空间离散格式采用二阶精度的AUSM+格式,湍流模型采用Menter’s SST k-ω模型。

    图8
                            分析模型几何外形及尺寸

    图8 分析模型几何外形及尺寸

    Fig.8 Geometry and sizes of analysis model

    图9
                            CFD模型的网格及边界条件

    图9 CFD模型的网格及边界条件

    Fig.9 Mesh and boudary conditions of CFD model

    获得的流场压力云图和温度云图如图10所示。数值解捕捉到了光滑的离体激波,锥体驻点处的压强超过了7 000 Pa,而正激波后气体温度则超过了950 K,与300 K的初始壁面温度形成强烈的温度梯度,产生很高的热流密度值。锥体表面压力分布情况如图11(a)所示,而锥体表面热流密度计算值和试验值分布情况如图11(b)所示,驻点处热流密度试验值为443.2 kW/m2,而本文的数值计算值为429.4 kW/m2,与试验值的相对误差为3.1%,除了驻点以外的其他区域热流密度计算值与试验值均吻合得很好,这验证了本文分析模型的正确性。

    图10
                            头部压力云图和温度云图

    图10 头部压力云图和温度云图

    Fig.10 Pressure and temperature fields of head area

    图11
                            物面压强和热流分布

    图11 物面压强和热流分布

    Fig.11 Pressure and heat flux distributions of wall

  • 5.2 分区迭代分析

    5.2

    结构瞬态传热分析的有限元模型如图12所示。结构材料为不锈钢,密度ρ=8 030 kg/m3,导热系数k=16.27 W/(m·K),比热容c=502.48 J/(kg·K),结构初始温度T0=300 K。采用以上CFD分析获得的初始流场进行气动热和结构传热的瞬态耦合分析,耦合方法为串行耦合迭代方法,耦合面数据传递采用虚拟空间插值法,耦合时间步Δt=0.000 1 s,瞬态分析总时间ttotal=100 s。

    图12
                            结构传热分析有限元网格

    图12 结构传热分析有限元网格

    Fig.12 Finite element mesh for structural heat transfer analysis

    分析获得了锥体结构前缘点(流场驻点)处的温度随时间变化的情况,如图13所示。从图13中可以观察到随着时间的推移,结构前缘点温度逐渐上升,这是由于外部连续的气动热作用于结构物面,导致结构温度持续上升。此外,随着时间的推移前,缘点温度上升的趋势逐渐减小,即曲线斜率逐渐减小,时间趋于无穷大时斜率将趋于0,结构温度场达到稳定(稳态解)。

    图13
                            驻点处温度随时间的变化情况

    图13 驻点处温度随时间的变化情况

    Fig.13 History of heat flux at stagnation point

    14为25,50,75和100 s时刻结构剖面的温度云图,从图14中可观察到各时刻结构最高温度分别为498.0,546.0,576.1和597.9 K,结构最高温度均在锥体前缘驻点处,且随着时间的推移高温区域越来越大。分析获得了不同时刻锥体表面压强和热流密度的分布情况,如图15所示,从图15中可观察到时间的变化对锥体表面的压强几乎没有影响,而表面热流密度却随时间的增加而降低。这主要是因为随着时间的推移锥体表面温度逐渐升高,激波层及边界层与壁面之间的温度梯度将逐渐减小,这就造成了锥体表面热流密度的降低;而壁面压强主要受控于流体力学而非热力学,故壁面压强几乎不改变。

    图14
                            不同时刻结构剖面的温度云图

    图14 不同时刻结构剖面的温度云图

    Fig.14 Temperature field of structural profile at different anlysis time

    图15
                            不同时刻壁面压力和热流分布

    图15 不同时刻壁面压力和热流分布

    Fig.15 Pressure and heat flux distributions at different time

  • 6 结 论

    6

    (1) 本文提出了气动热与结构传热的分区迭代推进分析方法,其考虑了气动热和结构传热之间的耦合效应,适用于本文所研究的高超声速钝头体热环境的精确分析。而传统的非耦合分析方法未考虑结构温度升高对气动热的反馈效应。

    (2) 进行了圆管验证算例分析,2 s时刻圆管驻点处的热流密度和温度的计算值与试验值的相对误差分别为1.34%和4.95%,从而验证了本文分区迭代推进分析方法的合理性与分析精度。

    (3) 最后研究了直二次圆锥体的热环境,锥体驻点温度随时间的推移而升高,且上升趋势逐渐变缓,最终趋于稳态值。此外时间的变化对锥体表面的压强几乎没有影响,而表面热流密度却随时间的增加而降低。

  • 参考文献

    • 1

      DEJARNETTE F R,HAMILTON H H.Aerodynamic heating on 3D bodies including the effects of entropy layer swallowing[J]. Journal of Spacecraft and Rockets,1975,12(1):5-12.

    • 2

      BRANDON H J,MASEK R V,DUNAVANT J C.Aerodynamic heating to corrugation stiffened struc-tures in thick turbulent boundary layers[J].AIAA Journal,1975,13(11):1460-1466.

    • 3

      INGER G R.Nonequilibrium boundary-layer effects on the aerodynamic heating of hypersonic waverider vehicles[J]. Journal of Thermophysics and Heat Transfer,1995,9(4):595-604.

    • 4

      MILOS F S,SQUIRE T H.Thermostructural analysis of X-34 wing leading-edge tile thermal protection system[J]. Journal of Spacecraft and Rockets,1999,36(2):189-198.

    • 5

      程克明,吕英伟.飞行器持续气动加热的耦合性分析[J].南京航空航天大学学报,2000,32(2):150-155.

      CHENG Keming,LÜ Yingwei.An analysis of coupling feature in continuing gasdynamic heating over flight vehicles[J]. Journal of Nanjing University of Aeronautics & Astronautics,2000,32(2):150-155.

    • 6

      任青梅,杨志斌,城主,等.气动加热与结构温度场耦合分析平台研发技术[J].强度与环境,2009,36(5):33-38.

      REN Qingmei,YANG Zhibin,CHENG Zhu,et al.Development of platform for analysis coupling aeroheating and structural temperature field[J]. Structure & Environment Engineering,2009,36(5):33-38.

    • 7

      季卫栋,王江峰,樊孝峰,等.高超声速流场与结构温度场一体化计算方法[J].航空动力学报,2016,31(1):153-160.

      JI Weidong,WANG Jiangfeng,FAN Xiaofeng,et al.Algorithms for hypersonic fluid-structural-thermal integrated[J].Journal of Aerospace Power,2016,31(1):153-160.

    • 8

      黄唐,毛国良,姜贵庆,等.二维流场、热、结构一体化数值模拟[J].空气动力学学报,2000,18(1):115-119.

      HUANG Tang,MAO Guoliang,JIANG Guiqing,et al.Two dimensional coupled flow-thermal-structural numerical simulation[J].Acta Aerodynamica Sinica,2000,18(1):115-119.

    • 9

      WIETING A R,DECHAUMPHAI P,BEY K S,et al.Application of integrated fluid-thermal-structural analysis methods[J].Thin-Walled Structures,1991,10(1):1-12.

    • 10

      LIOU M S.A sequel to AUSM: AUSM+[J].Journal of Computational Physics, 1996, 129(2):364-382.

    • 11

      LEER B V. Towards the ultimate conservative difference scheme V: A second-order sequel to Godunov's method[J]. Journal of Computational Physics,1979,32(1):101-136.

    • 12

      MENTER F R. Two-equation eddy-viscosity turbul-ence models for engineering applications[J]. AIAA Journal,1994,32(8):1598-1605.

    • 13

      DEESE J E, AGARWAL R K. Navier-Stokes calculations of transonic viscous flow about wing/body configurations[J].Journal of Aircraft,1988,25(12):1106-1112.

    • 14

      ZHAO Y.Computation of complex turbulent flow using matrix-free implicit dual time-stepping scheme and LRN turbulence model on unstructured grids[J].Computers & Fluids,2004,33(1):119-136.

    • 15

      WIETING A R.Experimental study of shock wave interference heating on a cylindrical leading edge[R].NASA-TM-100484, USA: NASA, 1987.

    • 16

      MILLER C G.Experimental and predicted heating distributions for biconics at incidence in air at Mach 10[R].NASA Technical Paper 2234,USA: NASA,1984.

黄杰

机 构:南京航空航天大学飞行器先进设计技术国防重点学科实验室,南京,210016

Affiliation:Key Laboratory of Fundamental Science for National Defense-Advanced Design Technology of Flight Vehicle, Nanjing University of Aeronautics & Astronautics, Nanjing, 210016, China

姚卫星

机 构:南京航空航天大学机械结构力学及控制国家重点实验室,南京,210016

Affiliation:State Key Laboratory of Mechanics and Control of Mechanical Structures, Nanjing University of Aeronautics & Astronautics, Nanjing, 210016, China

角 色:通讯作者

Role:Corresponding author

邮 箱:wxyao@nuaa.edu.cn

作者简介:姚卫星,男,教授,博士生导师,E-mail:wxyao@nuaa.edu.cn。

陈珺

角 色:中文编辑

Role:Editor

html/njhkht/201901014/alternativeImage/1e96c6b2-9f32-47aa-b022-c5e42fc939c6-F001.jpg
html/njhkht/201901014/alternativeImage/1e96c6b2-9f32-47aa-b022-c5e42fc939c6-F002.jpg
html/njhkht/201901014/alternativeImage/1e96c6b2-9f32-47aa-b022-c5e42fc939c6-F003.jpg
html/njhkht/201901014/alternativeImage/1e96c6b2-9f32-47aa-b022-c5e42fc939c6-F004.jpg
html/njhkht/201901014/alternativeImage/1e96c6b2-9f32-47aa-b022-c5e42fc939c6-F005.jpg
html/njhkht/201901014/alternativeImage/1e96c6b2-9f32-47aa-b022-c5e42fc939c6-F006.jpg
html/njhkht/201901014/alternativeImage/1e96c6b2-9f32-47aa-b022-c5e42fc939c6-F007.jpg
html/njhkht/201901014/alternativeImage/1e96c6b2-9f32-47aa-b022-c5e42fc939c6-F008.jpg
html/njhkht/201901014/alternativeImage/1e96c6b2-9f32-47aa-b022-c5e42fc939c6-F009.jpg
html/njhkht/201901014/alternativeImage/1e96c6b2-9f32-47aa-b022-c5e42fc939c6-F010.jpg
html/njhkht/201901014/alternativeImage/1e96c6b2-9f32-47aa-b022-c5e42fc939c6-F011.jpg
html/njhkht/201901014/alternativeImage/1e96c6b2-9f32-47aa-b022-c5e42fc939c6-F012.jpg
html/njhkht/201901014/alternativeImage/1e96c6b2-9f32-47aa-b022-c5e42fc939c6-F013.jpg
html/njhkht/201901014/alternativeImage/1e96c6b2-9f32-47aa-b022-c5e42fc939c6-F014.jpg
html/njhkht/201901014/alternativeImage/1e96c6b2-9f32-47aa-b022-c5e42fc939c6-F015.jpg

图1 气动热与结构传热之间的耦合关系

Fig.1 Coupled relationship between aerodynamic heating and structural heat transfer

图2 分区迭代推进策略

Fig.2 Division iterative marching strategy

图3 FVM节点与FEM节点对应关系

Fig.3 Nodes relationship between FVM and FEM

图4 流体和结构模型网格

Fig.4 Meshes for fluid and structural models

图5 驻点温度随时间变化情况

Fig.5 Temperature of stagnation point varying with time

图6 2 s时刻流体和结构温度场

Fig.6 Temperature fields for fluid and structural models at 2 s

图7 2 s时刻壁面热流密度和温度分布

Fig.7 Distribution for wall heat flux and temperature at 2 s

图8 分析模型几何外形及尺寸

Fig.8 Geometry and sizes of analysis model

图9 CFD模型的网格及边界条件

Fig.9 Mesh and boudary conditions of CFD model

图10 头部压力云图和温度云图

Fig.10 Pressure and temperature fields of head area

图11 物面压强和热流分布

Fig.11 Pressure and heat flux distributions of wall

图12 结构传热分析有限元网格

Fig.12 Finite element mesh for structural heat transfer analysis

图13 驻点处温度随时间的变化情况

Fig.13 History of heat flux at stagnation point

图14 不同时刻结构剖面的温度云图

Fig.14 Temperature field of structural profile at different anlysis time

图15 不同时刻壁面压力和热流分布

Fig.15 Pressure and heat flux distributions at different time

image /

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

  • 参考文献

    • 1

      DEJARNETTE F R,HAMILTON H H.Aerodynamic heating on 3D bodies including the effects of entropy layer swallowing[J]. Journal of Spacecraft and Rockets,1975,12(1):5-12.

    • 2

      BRANDON H J,MASEK R V,DUNAVANT J C.Aerodynamic heating to corrugation stiffened struc-tures in thick turbulent boundary layers[J].AIAA Journal,1975,13(11):1460-1466.

    • 3

      INGER G R.Nonequilibrium boundary-layer effects on the aerodynamic heating of hypersonic waverider vehicles[J]. Journal of Thermophysics and Heat Transfer,1995,9(4):595-604.

    • 4

      MILOS F S,SQUIRE T H.Thermostructural analysis of X-34 wing leading-edge tile thermal protection system[J]. Journal of Spacecraft and Rockets,1999,36(2):189-198.

    • 5

      程克明,吕英伟.飞行器持续气动加热的耦合性分析[J].南京航空航天大学学报,2000,32(2):150-155.

      CHENG Keming,LÜ Yingwei.An analysis of coupling feature in continuing gasdynamic heating over flight vehicles[J]. Journal of Nanjing University of Aeronautics & Astronautics,2000,32(2):150-155.

    • 6

      任青梅,杨志斌,城主,等.气动加热与结构温度场耦合分析平台研发技术[J].强度与环境,2009,36(5):33-38.

      REN Qingmei,YANG Zhibin,CHENG Zhu,et al.Development of platform for analysis coupling aeroheating and structural temperature field[J]. Structure & Environment Engineering,2009,36(5):33-38.

    • 7

      季卫栋,王江峰,樊孝峰,等.高超声速流场与结构温度场一体化计算方法[J].航空动力学报,2016,31(1):153-160.

      JI Weidong,WANG Jiangfeng,FAN Xiaofeng,et al.Algorithms for hypersonic fluid-structural-thermal integrated[J].Journal of Aerospace Power,2016,31(1):153-160.

    • 8

      黄唐,毛国良,姜贵庆,等.二维流场、热、结构一体化数值模拟[J].空气动力学学报,2000,18(1):115-119.

      HUANG Tang,MAO Guoliang,JIANG Guiqing,et al.Two dimensional coupled flow-thermal-structural numerical simulation[J].Acta Aerodynamica Sinica,2000,18(1):115-119.

    • 9

      WIETING A R,DECHAUMPHAI P,BEY K S,et al.Application of integrated fluid-thermal-structural analysis methods[J].Thin-Walled Structures,1991,10(1):1-12.

    • 10

      LIOU M S.A sequel to AUSM: AUSM+[J].Journal of Computational Physics, 1996, 129(2):364-382.

    • 11

      LEER B V. Towards the ultimate conservative difference scheme V: A second-order sequel to Godunov's method[J]. Journal of Computational Physics,1979,32(1):101-136.

    • 12

      MENTER F R. Two-equation eddy-viscosity turbul-ence models for engineering applications[J]. AIAA Journal,1994,32(8):1598-1605.

    • 13

      DEESE J E, AGARWAL R K. Navier-Stokes calculations of transonic viscous flow about wing/body configurations[J].Journal of Aircraft,1988,25(12):1106-1112.

    • 14

      ZHAO Y.Computation of complex turbulent flow using matrix-free implicit dual time-stepping scheme and LRN turbulence model on unstructured grids[J].Computers & Fluids,2004,33(1):119-136.

    • 15

      WIETING A R.Experimental study of shock wave interference heating on a cylindrical leading edge[R].NASA-TM-100484, USA: NASA, 1987.

    • 16

      MILLER C G.Experimental and predicted heating distributions for biconics at incidence in air at Mach 10[R].NASA Technical Paper 2234,USA: NASA,1984.

  • 您是第位访问者
    网站版权 © 南京航空航天大学学报
    技术支持:北京勤云科技发展有限公司