CN109840378A - 复杂发射条件火箭动态起飞燃气流场网格模型生成方法 - Google Patents

复杂发射条件火箭动态起飞燃气流场网格模型生成方法 Download PDF

Info

Publication number
CN109840378A
CN109840378A CN201910094844.XA CN201910094844A CN109840378A CN 109840378 A CN109840378 A CN 109840378A CN 201910094844 A CN201910094844 A CN 201910094844A CN 109840378 A CN109840378 A CN 109840378A
Authority
CN
China
Prior art keywords
rocket
subdomain
grid
dynamic
flow field
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910094844.XA
Other languages
English (en)
Other versions
CN109840378B (zh
Inventor
陈劲松
林禹
杜小坤
王明华
吴新跃
贾延奎
曾玲芳
余佳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Academy of Launch Vehicle Technology CALT
Beijing Institute of Space Launch Technology
Original Assignee
China Academy of Launch Vehicle Technology CALT
Beijing Institute of Space Launch Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Academy of Launch Vehicle Technology CALT, Beijing Institute of Space Launch Technology filed Critical China Academy of Launch Vehicle Technology CALT
Priority to CN201910094844.XA priority Critical patent/CN109840378B/zh
Publication of CN109840378A publication Critical patent/CN109840378A/zh
Application granted granted Critical
Publication of CN109840378B publication Critical patent/CN109840378B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Testing Of Engines (AREA)

Abstract

本发明公开了一种复杂发射条件火箭动态起飞燃气流场网格模型生成方法,首先基于火箭箭体(包含一级发动机喷管)、发射系统的简化实体模型生成包含火箭箭体、发射系统的包络计算域;然后基于火箭箭体、发射系统特征结构、特征尺度将整个包络计算域分区划分,生成火箭箭体、发射系统的包络计算子域,以及过渡子域及外围子域;接着以一级发动机喷口直径为参考尺度,生成子域网格;最后对动态网格模型质量、分辨率、可行性综合评估。本发明方法能够解决包括复杂发射系统、箭体结构条件下动态燃气流场网格模型生成难题,同时可有效保证网格质量和数量,节约了计算资源,提高了计算效率。

Description

复杂发射条件火箭动态起飞燃气流场网格模型生成方法
技术领域
本发明涉及复杂发射条件火箭动态起飞燃气流场网格模型开发方法,属于燃气流场仿真领域。
背景技术
火箭在起飞过程中,喷管喷出的高温高速燃气流会对底部发射系统产生剧烈的冲击和烧蚀作用,因此需要在设计阶段明确起飞过程发射系统所处的力、热环境,用于指导发射系统的结构强度和热防护设计。工程上通常采用燃气流场仿真来获取燃气流对发射系统的力、热影响。
为准确模拟火箭动态起飞过程的燃气流场特性,需要采取非定常燃气流场仿真方法,这就需要用到动网格生成技术。对复杂发射条件的动态燃气流场仿真而言,其动网格生成具有以下困难:1.发射系统结构复杂性造成网格质量难以保证;2.箭体与发射系统在高度方向有交错,底部生成动网格的包络面难以划分;3.箭体下方及侧向空间紧张,对包络箭体子域划分与子域内网格的生成难度较稳态网格生成更大。
为保证计算结果的准确性,通常保留流场内各结构件的真实构型,并控制结构件表面的网格尺度,这样会带来网格量偏大的后果,而且非定常计算的求解量远大于定常计算,所以动态起飞燃气流场非定常仿真计算往往需要大量的计算资源和漫长的时间。为了避免计算资源和时间的无谓浪费,必须合理控制燃气流场计算动网格模型的网格质量和数量。因此需要一种可行的复杂发射条件火箭动态起飞燃气流场网格模型生成方法。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供复杂发射条件火箭动态起飞燃气流场网格模型生成方法,解决了紧凑空间下,箭体底部包络面生成的难题,有效保证网格质量和数量,节约了计算资源,提高了计算效率。
本发明的技术解决方案是:
复杂发射条件火箭动态起飞燃气流场网格模型生成方法,包括如下步骤:
步骤1:对火箭箭体及发射系统的实体模型进行简化,所述火箭箭体包含一级发动机喷管;
步骤2:生成包含火箭箭体、发射系统的包络计算域,并划分子域;
步骤3:以一级发动机喷口直径为参考尺度,生成子域网格;
步骤4:对动态网格模型质量、分辨率、可行性综合评估,如果合理,则燃气流场网格模型生成结束,否则,对网格模型进行修正,直至网格模型合理。
所述步骤1中,简化包括以下几方面:
1)保留燃气流正冲范围的结构表面轮廓;
2)对结构表面小缝隙采取填充处理;
3)对背风区对流场几乎没有影响的凸起或凹陷外轮廓采取抹平处理;
4)对迎风区小尺度的结构凸起或凹陷采取抹平处理。
所述步骤2中,生成包络计算域的过程如下:
创建一个包含步骤1中简化实体模型的包络体,通过布尔运算将简化实体模型从体中差掉,得到的就是实体结构外的流场计算区域,即包含火箭箭体、发射系统的包络计算域。
所述步骤2中,划分子域的过程如下:
将包络计算域划分为运动子域和静止子域,运动子域是包络火箭箭体的子域,静止子域包括包络发射系统的子域和包络无特殊结构扰动的外围子域。
为保证动网格铺层质量,在火箭箭体与发射系统之间还设置特殊过渡子域。
所述步骤3中,子域网格生成包括运动子域和静止子域网格生成、特殊过渡子域网格动态生成。
运动子域和静止子域形成的区域中心为燃气流冲击和扰动核心区域;
运动子域和静止子域网格生成时,首先以喷管出口大小为参考,确定燃气流冲击和扰动核心区域网格尺度,核心区域网格尺度为喷管出口直径的1/25到1/100,随着网格位置与冲击扰动区域的距离增加,网格尺度逐渐增大。
基于火箭起飞弹道规划特殊过渡子域网格动态生成方法,当火箭起飞姿态接近垂直起飞状态时,采用在火箭底部特殊包络剖面铺层的方式生成动网格;当火箭起飞姿态严重偏离垂直起飞状态时,需要随计算进程重构生成特殊过渡子域网格。
在火箭底部特殊包络剖面铺层的方式生成动网格的方法如下:
在垂直火箭运动方向的过渡子域边界用呈“之字形”或“几字形”的特殊剖面包络,对位于特殊剖面包络之间的域采用铺层方式生成新网格,如果过渡子域包络于其它域之间时,平行火箭运动方向采用滑移交错网格差分求解的方式实现流场信息的传递。
动态网格模型质量的评估判据如下:网格内切圆与外接圆的比值大于0.1,则认为网格质量合格,否则不合格。
与现有技术相比,本发明具有如下有益效果:
(1)为保证网格质量,本发明通过简化模型来避免尖劈、缝隙等特征影响网格质量,同时可采用细化局部网格的方法来优化复杂结构区域的网格质量。
(2)本发明在定常燃气流场网格模型生成的基础上加入了对包络箭体的运动子域、包络发射系统的静止子域划分和箭体与发射系统间的特殊过渡子域划分。在生成过渡子域动网格铺层包络面时,为了解决箭体与发射系统在高度方向的交错问题,底部采取了组合“之字形”或组合“几”字形特殊剖面包络,对箭体和发射系统结构进行避让的同时保证了动网格能够正常在包络面铺层,解决了紧凑空间下,箭体底部包络面生成的难题。
(3)在子域网格生成时,以喷管出口大小为参考,确定燃气流冲击和扰动核心区域网格尺度,随着网格位置与冲击扰动区域的距离增加,网格尺度逐渐增大。从而使得本发明在保证足够的网格分辨率条件下,有效控制了网格数量,节约了计算资源,提高了计算效率。
附图说明
图1为本发明流程图;
图2为“之字形”包络面示意图;
图3为“几字形”包络面示意图;
图4为火箭起飞燃气流场网格模型示意图。
具体实施方式
本发明提出了一种复杂发射条件火箭动态起飞燃气流场网格模型生成方法,主要是在定常燃气流场网格模型生成的基础上加入了对包络箭体的运动子域划分和箭体与发射系统间的过渡子域划分。在生成过渡子域动网格铺层包络面时,为了解决箭体与发射系统在高度方向的交错问题,底部采取了组合“之字形”或组合“几”字形特殊剖面包络,对箭体和发射系统结构进行避让的同时保证了动网格能够正常在包络面铺层。
本发明具体步骤如图1所示,具体按如下步骤开展:
步骤1:对火箭箭体、一级发动机喷管及发射系统的实体模型进行简化
简化主要包括以下几方面:1.保留燃气流正冲范围的结构表面轮廓,以保证计算模型的准确性;2.对结构表面小缝隙采取填充处理;3.对背风区对流场几乎没有影响的凸起或凹陷外轮廓采取抹平处理;4.对迎风区小尺度的结构凸起或凹陷采取抹平处理。
步骤2:生成包含火箭箭体、发射系统的包络计算域,并基于火箭箭体、发射系统特征结构、特征尺度将整个包络计算域划分子域
在完成实体模型简化后,创建一个包含全部实体模型的包络体,通过布尔运算将实体模型从体中差掉,得到的就是实体结构外的流场计算区域,即包含火箭箭体、发射系统的包络计算域。
之后完成包络计算子域的划分,与定常计算的网格模型不同的是,动网格模型要将计算域划分为运动子域和静止子域,对火箭动态起飞流场计算而言,运动子域主要是包络火箭箭体的子域,静止子域包括包络发射系统的子域和无特殊结构扰动的外围包络子域。为保证动网格铺层质量,一般还要在火箭箭体与发射系统之间生成特殊过渡子域。
步骤3:子域网格生成
对上一步骤中划分的包络子域进行网格生成。对流场仿真计算而言,网格生成的重点是控制网格的分辨率和质量,在网格分辨率控制方面,以喷管出口大小为参考,确定燃气流冲击和扰动核心区域网格尺度,核心区域网格尺度可控制在喷管出口直径的1/25到1/100之间,而随着网格位置与冲击扰动区域的距离增加,网格尺度可逐渐增大。在网格质量控制方面,网格质量的评判标准有多种,本发明采用网格内切圆与外接圆的比值来判定,要求比值大于0.1。对于动态网格模型而言,由于要在箭体和发射装置间剖分包络面,留给网格调整的空间更小,网格质量更难控制。为保证网格质量,一方面通过简化模型来避免尖劈、缝隙等特征影响网格质量,另一方面可采用细化局部网格的方法来优化复杂结构区域的网格质量。
步骤4:基于火箭起飞弹道规划特殊过渡子域网格动态生成方法特殊过渡子域网格的生成方法根据火箭起飞弹道的不同可以分为两类,一种是火箭起飞姿态接近垂直起飞状态的情况,这种情况可以采用在火箭底部特殊包络剖面铺层的方式生成动网格,具体步骤是如图2和图3所示,先在垂直火箭运动方向的过渡子域边界用呈“之字形”或“几字形”的特殊剖面包络,其主要目的是对箭体和发射系统结构进行避让的同时保证动网格能够在火箭运动方向正常铺层,完成特殊剖面包络后可对位于特殊剖面包络之间的域采用铺层方式生成新网格,而如果过渡子域包络于其它域之间时,平行火箭运动方向采用滑移交错网格差分求解的方式实现流场信息的传递。对另一种火箭起飞姿态严重偏离垂直起飞状态的情况,则无法采用特殊包络面铺层的方式来生成动网格,需要随计算进程重构生成特殊过渡子域网格。
步骤5:对动态网格模型质量、分辨率、可行性综合评估
在完成全部网格生成后,对网格综合性能进行评估。对于网格质量和分辨率满足步骤3要求,且动网格能够顺利生成或重构的网格模型,则可开展动态起飞燃气流场仿真计算,若发现那一项没有满足要求,则继续对网格模型进行修正。如网格质量不满足要求,则采用细化局部网格的方法来优化复杂结构区域的网格质量。
生成的火箭起飞燃气流场网格模型示意图如图4所示。
本发明在定常燃气流场网格模型生成的基础上加入了对包络箭体的运动子域、包络发射系统的静止子域划分和箭体与发射系统间的特殊过渡子域划分。在生成过渡子域动网格铺层包络面时,为了解决箭体与发射系统在高度方向的交错问题,底部采取了组合“之字形”或组合“几”字形特殊剖面包络,对箭体和发射系统结构进行避让的同时保证了动网格能够正常在包络面铺层,解决了紧凑空间下,箭体底部包络面生成的难题。在子域网格生成时,以喷管出口大小为参考,确定燃气流冲击和扰动核心区域网格尺度,随着网格位置与冲击扰动区域的距离增加,网格尺度逐渐增大。从而使得本发明在保证足够的网格分辨率条件下,有效控制了网格数量,节约了计算资源,提高了计算效率。
本发明未详细说明部分属本领域技术人员公知常识。

Claims (10)

1.复杂发射条件火箭动态起飞燃气流场网格模型生成方法,其特征在于包括如下步骤:
步骤1:对火箭箭体及发射系统的实体模型进行简化,所述火箭箭体包含一级发动机喷管;
步骤2:生成包含火箭箭体、发射系统的包络计算域,并划分子域;
步骤3:以一级发动机喷口直径为参考尺度,生成子域网格;
步骤4:对动态网格模型质量、分辨率、可行性综合评估,如果合理,则燃气流场网格模型生成结束,否则,对网格模型进行修正,直至网格模型合理。
2.根据权利要求1所述的复杂发射条件火箭动态起飞燃气流场网格模型生成方法,其特征在于:所述步骤1中,简化包括以下几方面:
1)保留燃气流正冲范围的结构表面轮廓;
2)对结构表面小缝隙采取填充处理;
3)对背风区对流场几乎没有影响的凸起或凹陷外轮廓采取抹平处理;
4)对迎风区小尺度的结构凸起或凹陷采取抹平处理。
3.根据权利要求1所述的复杂发射条件火箭动态起飞燃气流场网格模型生成方法,其特征在于:所述步骤2中,生成包络计算域的过程如下:
创建一个包含步骤1中简化实体模型的包络体,通过布尔运算将简化实体模型从体中差掉,得到的就是实体结构外的流场计算区域,即包含火箭箭体、发射系统的包络计算域。
4.根据权利要求3所述的复杂发射条件火箭动态起飞燃气流场网格模型生成方法,其特征在于:所述步骤2中,划分子域的过程如下:
将包络计算域划分为运动子域和静止子域,运动子域是包络火箭箭体的子域,静止子域包括包络发射系统的子域和包络无特殊结构扰动的外围子域。
5.根据权利要求4所述的复杂发射条件火箭动态起飞燃气流场网格模型生成方法,其特征在于:为保证动网格铺层质量,在火箭箭体与发射系统之间还设置特殊过渡子域。
6.根据权利要求5所述的复杂发射条件火箭动态起飞燃气流场网格模型生成方法,其特征在于:所述步骤3中,子域网格生成包括运动子域和静止子域网格生成、特殊过渡子域网格动态生成。
7.根据权利要求6所述的复杂发射条件火箭动态起飞燃气流场网格模型生成方法,其特征在于:运动子域和静止子域形成的区域中心为燃气流冲击和扰动核心区域;
运动子域和静止子域网格生成时,首先以喷管出口大小为参考,确定燃气流冲击和扰动核心区域网格尺度,核心区域网格尺度为喷管出口直径的1/25到1/100,随着网格位置与冲击扰动区域的距离增加,网格尺度逐渐增大。
8.根据权利要求6所述的复杂发射条件火箭动态起飞燃气流场网格模型生成方法,其特征在于:基于火箭起飞弹道规划特殊过渡子域网格动态生成方法,当火箭起飞姿态接近垂直起飞状态时,采用在火箭底部特殊包络剖面铺层的方式生成动网格;当火箭起飞姿态严重偏离垂直起飞状态时,需要随计算进程重构生成特殊过渡子域网格。
9.根据权利要求8所述的复杂发射条件火箭动态起飞燃气流场网格模型生成方法,其特征在于:在火箭底部特殊包络剖面铺层的方式生成动网格的方法如下:
在垂直火箭运动方向的过渡子域边界用呈“之字形”或“几字形”的特殊剖面包络,对位于特殊剖面包络之间的域采用铺层方式生成新网格,如果过渡子域包络于其它域之间时,平行火箭运动方向采用滑移交错网格差分求解的方式实现流场信息的传递。
10.根据权利要求1所述的复杂发射条件火箭动态起飞燃气流场网格模型生成方法,其特征在于:动态网格模型质量的评估判据如下:网格内切圆与外接圆的比值大于0.1,则认为网格质量合格,否则不合格。
CN201910094844.XA 2019-01-31 2019-01-31 复杂发射条件火箭动态起飞燃气流场网格模型生成方法 Active CN109840378B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910094844.XA CN109840378B (zh) 2019-01-31 2019-01-31 复杂发射条件火箭动态起飞燃气流场网格模型生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910094844.XA CN109840378B (zh) 2019-01-31 2019-01-31 复杂发射条件火箭动态起飞燃气流场网格模型生成方法

Publications (2)

Publication Number Publication Date
CN109840378A true CN109840378A (zh) 2019-06-04
CN109840378B CN109840378B (zh) 2023-02-03

Family

ID=66884521

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910094844.XA Active CN109840378B (zh) 2019-01-31 2019-01-31 复杂发射条件火箭动态起飞燃气流场网格模型生成方法

Country Status (1)

Country Link
CN (1) CN109840378B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112307558A (zh) * 2020-10-23 2021-02-02 中国运载火箭技术研究院 一种三维曲面导流通道、导焰舱及导流通道的设计方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104376151A (zh) * 2014-10-30 2015-02-25 北京宇航系统工程研究所 一种火箭发动机真空干扰羽流场仿真方法
CN105354401A (zh) * 2015-12-24 2016-02-24 中国人民解放军装备学院 一种多喷管火箭或导弹尾焰流场计算方法
US9977848B1 (en) * 2014-07-10 2018-05-22 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Method and system for predicting rocket nozzle deformation during engine start-up and shut-down transients
CN108304684A (zh) * 2018-03-22 2018-07-20 北京航天动力研究所 一种火箭发动机尾喷射流仿真方法及系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9977848B1 (en) * 2014-07-10 2018-05-22 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Method and system for predicting rocket nozzle deformation during engine start-up and shut-down transients
CN104376151A (zh) * 2014-10-30 2015-02-25 北京宇航系统工程研究所 一种火箭发动机真空干扰羽流场仿真方法
CN105354401A (zh) * 2015-12-24 2016-02-24 中国人民解放军装备学院 一种多喷管火箭或导弹尾焰流场计算方法
CN108304684A (zh) * 2018-03-22 2018-07-20 北京航天动力研究所 一种火箭发动机尾喷射流仿真方法及系统

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112307558A (zh) * 2020-10-23 2021-02-02 中国运载火箭技术研究院 一种三维曲面导流通道、导焰舱及导流通道的设计方法
CN112307558B (zh) * 2020-10-23 2024-02-09 中国运载火箭技术研究院 一种三维曲面导流通道、导焰舱及导流通道的设计方法

Also Published As

Publication number Publication date
CN109840378B (zh) 2023-02-03

Similar Documents

Publication Publication Date Title
Rotunno On the evolution of thunderstorm rotation
CN110688791B (zh) 一种钝头体流场激波自适应结构网格生成方法
CN106406359A (zh) 一种基于虚拟目标的固定翼无人机跟踪地面目标制导方法
McCroskey Vortex wakes of rotorcraft
CN109885908A (zh) 一种新型羽翅仿生通风扑翼系统及多涡干扰机理分析方法
Preidikman et al. Time-domain simulations of linear and nonlinear aeroelastic behavior
CN109376403A (zh) 一种基于笛卡尔自适应重构技术的二维结冰数值模拟方法
CN116384290B (zh) 一种考虑真实气体效应的高超声速飞行器动导数预测方法
Oruc et al. Coupled flight dynamics and CFD simulations of helicopter/ship dynamic interface
CN105676672A (zh) 扑翼飞行器复合飞行策略仿真建模方法
CN109840378A (zh) 复杂发射条件火箭动态起飞燃气流场网格模型生成方法
CN112613246A (zh) 一种固体火箭发动机在飞行过载下的两相流仿真方法
CN107748494A (zh) 一种超空泡航行体运动状态转移控制方法
CN109871586A (zh) 适用于共轴旋翼地面效应的气动干扰数值确定方法和系统
CN110104164A (zh) 一种用于跨声速机翼的前加载-吸气组合流动控制方法
CN109871603A (zh) 复杂发射技术条件的燃气流场预示方法
CN115983075A (zh) 一种无人机回收双向流固耦合的数值模拟方法
Hariharan et al. Helicopter Rotor Aerodynamic Modeling in Hover Based on First-Principles: State-of-the-Art and Remaining Challenges
Hosseinjani et al. Numerical simulation of the wake structure and thrust/lift generation of a pitching airfoil at low reynolds number via an immersed boundary method
CN113935114A (zh) 一种用于水上飞机着水数值仿真的规则波造波方法
Syono On the motion of a typhoon (I)
CN109900154A (zh) 一种火箭发射燃气流场准稳态预示方法及装置
CN109858150A (zh) 复杂发射技术条件的燃气流场网格模型生成方法
CN109858153A (zh) 火箭发射燃气流冲击条件下流固耦合边界条件确定方法
Xu et al. Investigation of active flow control on aerodynamic performance of HALE UAV airfoil

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant