多重网格法求解火焰稳定器钝头体管流问题

多重网格法是近年发展起来的一种软加速方法,用多重网格法可明显地加快数值计算的收敛速度,且网格越密其加速效果越明显。应用多重网格法求解有钝头体的火焰稳定器管内流动问题,计算结果与试验结果比较,吻合较好。     

Theoretical investigation of shock stand-off distance for non-equilibrium flows over spheres

We derived a theoretical solution of the shock stand-off distance for a non-equilibrium flow over spheres based on Wen and Hornung’s solution and Olivier’s solution. Compared with previous approaches, the main advantage of the present approach is allowing an analytic solution without involving any semi-empirical parameter for the whole non-equilibrium flow regimes. The effects of some important physical quantities therefore can be fully revealed via the analytic solution. By combining the current solution with Ideal Dissociating Gas (IDG) model, we investigate the effects of free stream kinetic energy and free stream dissociation level (which can be very different between different facilities) on the shock stand-off distance.     

乘波构型的钝化方法及其对性能影响研究

前缘钝化是解决乘波构型飞行器气动热问题的有效方法之一.按照乘波构型的设计特点,对已有的两种钝化方法分别进行了改进.采用CFD方法分析了前缘钝化及不同钝化半径对乘波构型性能的影响,得出了乘波构型气动力和气动热性能参数随钝化半径的变化规律.计算结果表明:在相同的钝化半径下,按改进的Tincher方法钝化后的乘波构型与按改进的Takashima方法钝化后的乘波构型相比:升阻比大、总的表面积小、最大热流密度基本一样,非驻点区乘波构型前缘的热流密度峰值较大.因而按改进的Tincher方法钝化后的乘波构型气动性能明显好于按改进的Takashima方法钝化后的乘波构型,而气动热性能则略差于后者.分析表明:钝化后的乘波构型性能不仅与钝化半径有关,而且受钝化方法的影响也很大.在对高超声速乘波飞行器进行布局设计时,应针对乘波构型的设计特点,采用合适的钝化方法,综合考虑钝化方法和钝化半径对气动力和气动热性能的影响效应,寻找最佳的钝化方案.研究结论可为高超声速乘波飞行器的外形设计提供一定的依据.     

钝体绕流流场分离涡的数值模拟

 采用近似因式分解和特殊Jacobin系数矩阵分裂法,求解可压缩Navier-Stokes方程,数值模拟矩形柱体绕流流场中分离涡的形成和发展。给出了Ma_∞=0.3,0.8;Re=10~4,10~6的计算结果。从结果中可清楚地看到柱体非对称涡形成的非定常过程。     

高超声速飞行器激波控制减阻技术

针对高超声速飞行器巨大的激波阻力，采用数值方法研究了由钝头体、气动杆和侧向喷流构成的组合模型的减阻性能。侧向喷流将弓形激波推离气动杆，组合模型的再附激波明显弱于传统气动杆模型，其阻力系数比气动杆模型低了33.52%，从而验证了本文组合模型优异的减阻效率。进行了组合模型的影响因素分析，随侧向喷流总压比和气动杆的长度的增加，再附激波强度减弱，减阻效率升高，但减阻效率的变化速率逐渐减小。随喷口位置向下游移动，再附激波逐渐增强，减阻效率降低，且减阻效率的变化速率逐渐增加。此外本文还研究了以上参数对流场结构及钝头体压力峰值位置的影响。     

飞行高度对高超声速钝锥边界层稳定性及转捩的影响

研究了飞行高度对高超声速钝锥边界层稳定性及转捩的影响。通过求解三维可压缩Navier-Stokes方程计算了来流Ma=6，半锥角为7°的钝锥在飞行高度20~40 km条件下的基本流场，利用线性稳定性理论（LST）研究了飞行高度对钝锥边界层流动稳定性的影响，最后采用e N方法进行了转捩预测。研究发现，随着飞行高度的增加，流向不稳定N s值和横流不稳定N cf 值均减小，由横流不稳定性引起的圆锥表面大部分区域转捩逐渐转变为流向扰动引起迎风面转捩横流扰动引起背风面转捩，继而横流扰动消失,流向不稳定波引起迎风面转捩。     

再入钝头体近尾流动计算

本文采用联立求解粘性剪切层和有旋元粘区的方法计算了再入钝头体层流近尾流场，有旋无粘区采用有旋特征线法计算，粘性剪切层采用改进的流管法计算，这两个区域间的基本流线上的压力，温度，速度等由两个区域迭代匹配确定，粘性剪切层的下边界区分流线采用的实验确定，回避了回流区复杂计算，计算结果与有关试验和文献一致，本文的目的是提供远尾流计算所需的颈部初剖面。     

高超音速化学非平衡层流粘性激波层钝头细长体绕流数值计算

本文研究了高超音速化学非平衡钝头细长体绕流的计算方法,给出了耦合求解连续方程和法向动量方程的具体做法,以及在高空情况下考虑滑移边界条件后的数学处理。计算结果表明,对钝头细长体,可以将计算推进到三十倍头部半径的后身区域。     

高超声速化学平衡流辐射流场数值解

本文用数值方法研究了高超声速大钝头体绕流辐射流场。采用粘性激波层控制方程和平衡化学反应气体模型。辐射部份包括连续辐射和线辐射。采用切平板辐射输运模型解决了空间推进的问题，计算结果表明：辐射传热与飞行器再入飞行速度，高度及头部半径密切相关，在计算的条件下，当再入速度小于７．６公里／秒时，辐射传热的影响可以忽略，当再入速度为１１公里／秒左右时，辐射传热可以与对流传热相比拟。当再入速度更大时，辐射传热远     

微波电热推力器钝头体方案的性能计算

进行了微波电热推力器（ＭＥＴ）钝头体方案的性能计算,对微波作了谐振模式ＴＭ０１１和波导模式ＴＭ０１两种选择。将计入电磁场作用的轴对称钝头体粘性流动方程,结合ＴＭ０１１或ＴＭ０１模式,与等离体区耦合求解,获得了ＴＥＴ比冲和推力等性能,结果与文献的实验及计算结果趋向一致。