三维超声速压力反问题的特征线求解技术

推进技术 ›› 2018, Vol. 39 ›› Issue (10): 2340-2350.

三维超声速压力反问题的特征线求解技术

赵玉新¹,蓝庆生¹,赵一龙¹,刘红阳²

国防科技大学高超声速冲压发动机技术重点实验室，湖南长沙 410073,国防科技大学高超声速冲压发动机技术重点实验室，湖南长沙 410073,国防科技大学高超声速冲压发动机技术重点实验室，湖南长沙 410073,中国空气动力研究与发展中心计算空气动力研究所，四川绵阳 621000

发布日期:2021-08-15
作者简介:赵玉新，男，博士，教授，研究领域为高/超声速气动设计、超声速边界层、超声速流动成像测量等。 E-mail:zyx_nudt@163.com 通讯作者:蓝庆生，男，硕士生，研究领域为高/超声速气动设计。
基金资助:
国家自然科学基金(11472304)。

A Characteristic Method for Solving Three-Dimensional Supersonic Pressure Inverse Problems

Science and Technology on Scramjet Laboratory，National University of Defense Technology，Changsha 410073，China,Science and Technology on Scramjet Laboratory，National University of Defense Technology，Changsha 410073，China,Science and Technology on Scramjet Laboratory，National University of Defense Technology，Changsha 410073，China and Computational Aerodynamics Institutes，China Aerodynamics Research and Development Center， Mianyang 621000，China

Published:2021-08-15

摘要/Abstract

摘要： 为了进一步探索三维超声速气动反问题的求解方法，基于双特征线理论及Butler解法，研究了三维超声速流线压力反问题的适定性。为了确保解的唯一性，提出在限定壁面膨胀/压缩方向时存在壁面压力与三维坐标的一一映射关系。基于该映射关系，提出了三维压力反问题的双特征线求解技术(iMOC-3D 求解器)。采用Prandtl-Meyer膨胀波、Busemann进气道的理论解，对iMOC-3D求解方法的膨胀、压缩过程进行了精度评估，误差均为1×10-4量级。为了进一步验证设计方法的可靠性和易控性，设计了进口为矩形和三角形的超声速喷管；通过设计壁面压力分布，完成了均匀膨胀的轴对称喷管设计，并将设计结果与数值模拟进行对比验证。研究表明:预设三维流线下游未知点的压力值，存在多个流动方向满足该压力条件，即该问题的解不唯一，因此三维超声速流线压力反问题是非适定的。对比验证表明:所设计的流场与CFD计算得到的等值线符合得较好，流场参数的最大误差为1%。因此，所提出的双特征线解法具有一定的可靠性，有望为三维超声速流道设计提供新思路。

关键词: 三维压力反问题；双特征线方法；三维超声速气动设计；数值方法

Abstract: In order to explore valid methods for three-dimensional supersonic aerodynamic inverse problem， well-posedness of three-dimensional supersonic streamline pressure inverse problem was explored based on the concept of Bicharacteristics and Butler algorithm. For ensuring the uniqueness of the solution of this inverse problem， a statement that the one-to-one mapping relationship between wall pressure and three-dimensional coordinates was existent with a local constrained expansion or compression direction was proposed. According to the numerical solution of the mapping relationship， an abbreviated iMOC-3D solver for the pressure inverse problem based on Bicharacteristics method was established. The accuracy order of expansion and compression progress of the scheme was tested by comparing numerical solutions with analytical results of Prandtl-Meyer expansion wave and Busemann flow. The reliability of iMOC-3D solution was proved to some extent by the order of relative error being 1×10-4. In order to further verify the reliability and controllability of the design method， supersonic nozzles with different inlets such as rectangle and triangle were designed， and an axisymmetric nozzle with prescribed pressure distribution was designed， which was verified by CFD numerical simulations. Research has shown that for three dimensional streamlines， there are multiple flow directions to meet the downstream pressure condition that has been preset. So the solution of this inverse problem is not unique. Thus， inverse problem for three-dimensional supersonic flow is not well-posed. Comparing simulations between CFD and the solver， it would be found that the solver’s numerical solutions agree well with CFD simulations and the maximum error of flow field is 1%. Thus， the proposed method would be reliable and might provide new ideas for three-dimensional supersonic aerodynamic design.

Key words: Three-dimensional pressure inverse problem；Bicharacteristics method；Three-dimensional supersonic aerodynamic design；Numerical method

赵玉新,蓝庆生,赵一龙,刘红阳. 三维超声速压力反问题的特征线求解技术[J]. 推进技术, 2018, 39(10): 2340-2350.

ZHAO Yu-xin¹,LAN Qing-sheng¹,ZHAO Yi-long¹,LIU Hong-yang². A Characteristic Method for Solving Three-Dimensional Supersonic Pressure Inverse Problems[J]. Journal of Propulsion Technology, 2018, 39(10): 2340-2350.

参考文献

[1] Beckel S, Garrett J, Gettinger C.Technologies for Robust and Affordable Scramjet Propulsion[R]. AIAA 2006-7980.
[2] Bulman M, Siebenhaar A. The Rebirth of Round Hypersonic Propulsion[R]. AIAA 2006-5035.
[3] 易仕和, 赵玉新, 何霖, 等. 超声速与高超声速喷管设计[M]. 北京:国防工业出版社, 2013.
[4] Murthy S N B, Thompson H D. Design of Optimized Three-Dimensional Nozzles[J]. Journal of Spacecraft & Rockets, 1966, 3(9): 1384-1393.
[5] Wang X, Damodaran M. Optimal Three-Dimensional Nozzle Shape Design Using CFD and Parallel Simulated Annealing[J]. Journal of Propulsion & Power, 2015, 18(1): 217-221.
[6] Xing X Q, Damodaran M. Design of Three-Dimensional Nozzle Shapes Using NURBS, CFD and Hybrid Optimization Strategies[R]. AIAA 2004-4368.
[7] Moelder S, Szpiro E J. Busemann Inlet for Hypersonic Speeds[J]. Journal of Spacecraft & Rockets, 2015, 3(8): 1303-1304.
[8] Jacobsen L S, Tam C J, Behdadnia R, et al. Starting and Operation of a Streamline-Traced Busemann Inlet at Mach 4[R]. AIAA 2006-4508.
[9] Smart M K. Design of Three-Dimensional Hypersonic Inlets with Rectangular-to-Elliptical Shape Transition [J]. Journal of Propulsion and Power, 1999, 15(3): 408-416.
[10] Gollan R J, Smart M K. Design of Modular Shape-Transition Inlets for a Conical Hypersonic Vehicle[J]. Journal of Propulsion and Power, 2013, 29(4): 832-838.
[11] 尤延铖, 梁德旺, 黄国平. 一种新型内乘波进气道初步研究[J]. 推进技术, 2006, 27(3): 252-256. (YOU Yan-cheng, LIANG De-wang, HUANG Guo-ping. Investigation of Internal Waverider-Derived Hypersonic Inlet[J]. Journal of Propulsion Technology, 2006, 27(3): 252-256.)
[12] Billig F S, Kothari A P. Streamline Tracing: Technique for Designing Hypersonic Vehicles[J]. Journal of Propulsion and Power, 2000, 16(3): 465-471.
[13] 卢鑫, 岳连捷, 肖雅彬, 等. 超燃冲压发动机尾喷管流线追踪设计[J]. 推进技术, 2011, 32(1). (LU Xin, YUE Lian-jie, XIAO Ya-bin, et al. Design of Scramjet Nozzle Based on Streamline Tracing Technique[J]. Journal of Propulsion Technology, 2011, 32(1).)
[14] Sobieczky H, Dougherty F C, Jones K. Hypersonic Waverider Design from Given Shock Waves[J]. International Hypersonic Waverider Symposium, University of Maryland Collage Park, 1990, 10: 17-19
[15] Sobieczky H, Zores B, Wang Z. et al. High Speed Flow Design Using Osculating Axisymmetric Flows[J]. Aerospace Science and Technology, 1997, 3: 182-187.
[16] 赵博. 乘波体飞行器气动外形设计及其数值验证[D]. 哈尔滨:哈尔滨工程大学, 2013.
[17] 尤延铖, 梁德旺. 基于内乘波概念的三维变截面高超声速进气道[J]. 中国科学(E 辑), 2009, (8): 1483-1494.
[18] 卢鑫, 岳连捷, 肖雅彬, 等. 超燃冲压发动机三维变截面尾喷管设计[C]. 无锡:第三届高超声速科技学术会议, 2010.
[19] 赵玉新. 超声速混合层时空结构的实验研究[D]. 长沙:国防科技大学, 2008.
[20] 全志斌, 徐惊雷, 莫建伟. 基于控制壁面压力分布的分级单边膨胀喷管设计及试验验证[J]. 推进技术, 2013, 34(3): 307-310. (QUAN Zhi-bin, XU Jing-lei, MO Jian-wei. Design and Experiment Validation of Single Expansion Ramp Nozzle Based on the Control of Wall Pressure Distribution[J]. Journal of Propulsion Technology, 2013, 34(3): 307-310.)
[21] 赵玉新, 刘红阳, 赵一龙. 三维超声速流动的反问题[C]. 哈尔滨:第八届全国高超声速科技学术会议, 2015.
[22] 刘红阳. 特征线追踪方法和应用 [D]. 长沙:国防科技大学, 2015.
[23] 李怡庆, 韩伟强, 尤延铖, 等. 压力分布可控的高超声速进气道/前体一体化乘波设计[J]. 航空学报, 2016, 37(9), 2711-2720.
[24] Cline M C, Hoffman J D. Comparison of Characteristic Schemes for Three-Dimensional, Steady, Isentropic Flow[J]. AIAA Journal, 1972, 10(11): 1452-1458.
[25] Ransom V H, Hoffman J D, Thompson H D. A Second-Order Numerical Method of Characteristics for Three-Dimensional Supersonic Flow[D]. USA: Purdue University, 1970.
[26] Ransom V H, Hoffman J D, Thompson HD. A Second- Order Bicharacteristics Method for Three-Dimensional, Steady, Supersonic Flow[J]. AIAA Journal, 1972, 10(12): 1573-1581.
[27] Butler D S. The Numerical Solution of Hyperbolic Systems of Partial Differential Equations in Three Independent Variables[J]. Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Science. The Royal Society, 1960, 255(1281): 232-252.
[28] 左克罗MJ, 霍夫曼JD. 气体动力学[M]. 北京:国防工业出版社, 1984.
[29] 童秉纲, 孔祥言, 邓国华. 气体动力学[M]. 北京:高等教育出版社, 1990.
[30] 蓝庆生, 赵玉新, 赵一龙, 等. 三维超声速流动的压力反问题[J]. 空气动力学学报, 2017, 35(3): 429-435.(编辑:张荣莉)　　　　　　　　　　　　　*　收稿日期:2017-10-09；修订日期:2017-11-13。基金项目:国家自然科学基金(11472304)。作者简介:赵玉新,男,博士,教授,研究领域为高/超声速气动设计、超声速边界层、超声速流动成像测量等。 E-mail:zyx_nudt@163.com通讯作者:蓝庆生,男,硕士生,研究领域为高/超声速气动设计。E-mail: lanqingsheng2015@163.com