CN113127795A - 一种近地空间飞行器反照地球热流的解析方法 - Google Patents

一种近地空间飞行器反照地球热流的解析方法 Download PDF

Info

Publication number
CN113127795A
CN113127795A CN202110393505.9A CN202110393505A CN113127795A CN 113127795 A CN113127795 A CN 113127795A CN 202110393505 A CN202110393505 A CN 202110393505A CN 113127795 A CN113127795 A CN 113127795A
Authority
CN
China
Prior art keywords
earth
heat flow
albedo
value
numerical
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
CN202110393505.9A
Other languages
English (en)
Other versions
CN113127795B (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.)
Suzhou Kachuang Information Technology Co ltd
Original Assignee
Suzhou Kachuang Information Technology Co ltd
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 Suzhou Kachuang Information Technology Co ltd filed Critical Suzhou Kachuang Information Technology Co ltd
Priority to CN202110393505.9A priority Critical patent/CN113127795B/zh
Publication of CN113127795A publication Critical patent/CN113127795A/zh
Application granted granted Critical
Publication of CN113127795B publication Critical patent/CN113127795B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

本发明涉及近地空间飞行器热分析技术领域，具体地说，涉及一种近地空间飞行器反照地球热流的解析方法，包括以下步骤：(1)构建地球反照热流双重积分一般表达式；(2)构建飞行器表面平行于地球球心截面、太阳辐射垂直飞行器表面的地球反照热流一般表达式；(3)构建地球反照热流解析表达式；(4)计算地球反照热流密度精确值和数值积分值；(5)误差分析与数值网格划分数选择。本发明将为近地空间飞行器热分析领域在地球反照热流密度数值计算精度评估中提供计算方法。

Description

一种近地空间飞行器反照地球热流的解析方法
技术领域
本发明涉及近地空间飞行器热分析技术领域，具体地说，涉及一种近地空间飞行器反照地球热流的解析方法。
背景技术
计算近地空间环境热流是飞行器热力分析的重要前提，而地球反照热流是空间环境热流的重要组成，目前均采用数值积分、Monte-Carlo等数值分析方法计算。为了更精确地进行热力分析，所采用数值计算方法的精度成为关键。解析法可得到各种数理分析的精确值，为数值分析方法提供精度验证依据。在相关文献中发现，现有技术方法并没有利用解析法数据对数值分析精确度进行验证。
发明内容
本发明的内容是提供一种近地空间飞行器反照地球热流的解析方法，其能够克服现有技术的某种或某些缺陷。
根据本发明的一种近地空间飞行器反照地球热流的解析方法，其包括以下步骤：
(1)构建地球反照热流双重积分一般表达式；
将地球视为漫反射规则灰球体，其表面按经纬线的方式近似划分为一系列的平面面元，用Γv表示对近地空间目标可视的地球球冠，ΓL表示被太阳照射的区域，则Γv∩ΓL是地球反照辐射的区域；dA1是地球表面上的一个微面元，U是dA1的几何中心，
Figure BDA0003017682660000011
是dA1的外向单位法矢量，则飞行器某表面平面dAt上的地球反照热流密度ESE可表示为：
Figure BDA0003017682660000021
式中，ES表示“太阳常数”，ρE表示地球反照率，O′是dAt的几何中心，φSE是
Figure BDA0003017682660000022
与太阳辐射向量
Figure BDA0003017682660000023
相反方向向量之间的夹角，βSE是
Figure BDA0003017682660000024
与
Figure BDA0003017682660000025
之间的夹角，γSE是dAt法向量
Figure BDA0003017682660000026
与
Figure BDA0003017682660000027
之间的夹角；
(2)构建飞行器表面平行于地球球心截面、太阳辐射垂直飞行器表面的地球反照热流一般表达式；
对于dAt平行于地球球心截面，太阳辐射向量
Figure BDA0003017682660000028
垂直于dAt的情形，给出Γv对dAt的照度表达式为：
Figure BDA0003017682660000029
式中忽略dAt在地球上投影的影响；
(3)构建地球反照热流解析表达式；
根据
Figure BDA00030176826600000210
与地球半径RE之间的关系，将照度表达式改写为：
Figure BDA00030176826600000211
因为：
Figure BDA00030176826600000212
式中
Figure BDA00030176826600000213
为dA1在地心直角坐标系投影的方位角；
作以下推导：
Figure BDA00030176826600000214
式中O表示地球球心，H表示O与O′之间的距离；
将
Figure BDA0003017682660000031
写为lUO′，则：
Figure BDA0003017682660000032
cosβSE可表示为：
Figure BDA0003017682660000033
又因为：
Figure BDA0003017682660000034
则ESE可作下列推导：
Figure BDA0003017682660000035
令ζ＝cosφSE，则：
Figure BDA0003017682660000036
令REH＝a，RE 2+H2＝b，则上式变为：
Figure BDA0003017682660000037
令
Figure BDA0003017682660000041
则：
Figure BDA0003017682660000042
Figure BDA0003017682660000043
Figure BDA0003017682660000044
(4)计算地球反照热流密度精确值和数值积分值；
(5)误差分析与数值网格划分数选择。
作为优选，步骤(4)中，地球反照热流密度精确值的计算方法为：
设定H值，将ES、ρE、RE代入地球反照热流解析表达式，即可计算H值高度下的地球反照热流密度精确值ESE-A。
作为优选，步骤(4)中，数值积分值的计算方法为：
将积分式中的积分区域Γv∩ΓL以类似地球经纬线的方式进行网格划分，其中经向m份、纬向n份，共得到m×n份，则数值积分值ESE-B为：
Figure BDA0003017682660000045
式中下标i表示第i个dA1以及其相关量；当m×n值越大，积分区域Γv∩ΓL划分越密，数值积分值ESE-B精度越高。
作为优选，步骤(5)中，地球反照热流密度精确值ESE-A与数值积分值ESE-B间的相对误差RESE为：
Figure BDA0003017682660000046
当RESE越接近于1，数值积分精度越高；根据相对误差要求，选择网格划分值m×n。
本发明基于地球反照热流一般计算表达式，考虑飞行器表面与地球、太阳辐射方位之间的特殊关系，对一般计算表达式进行简化，而后利用几何关系与积分原理，严格推导出地球反照热流解析表达式。本发明将为近地空间飞行器热分析领域在地球反照热流密度数值计算精度评估中提供计算方法。
本发明具体是将地球视为漫反射规则灰球体，飞行器表面平行于球心截面，太阳辐射垂直飞行器表面，建立地球反照热流密度解析表达式，进而精确分析飞行器表面地球反照热流密度，此方法可广泛应用于近地空间环境分析、飞行器热力分析等领域，分析地球反照热流密度数值方法的精度。
附图说明
图1为实施例1中一种近地空间飞行器反照地球热流的解析方法构建技术流程图；
图2为实施例1中地球视为漫反射规则灰球体的示意图；
图3为实施例1中构建地球反照热流双重积分一般表达式示意图；
图4为实施例1中构建地球反照热流解析表达式示意图。
具体实施方式
为进一步了解本发明的内容，结合附图和实施例对本发明作详细描述。应当理解的是，实施例仅仅是对本发明进行解释而并非限定。
实施例1
如图1所示，本实施例提供了一种近地空间飞行器反照地球热流的解析方法，其包括以下步骤：
(1)构建地球反照热流双重积分一般表达式；
如图2所示，将地球视为漫反射规则灰球体，其表面按经纬线的方式近似划分为一系列的平面面元，用Γv表示对近地空间目标可视的地球球冠，ΓL表示被太阳照射的区域，则Γv∩ΓL是地球反照辐射的区域；如图3所示，dA1是地球表面上的一个微面元，U是dA1的几何中心，
Figure BDA0003017682660000061
是dA1的外向单位法矢量，则飞行器某表面平面dAt上的地球反照热流密度ESE可表示为：
Figure BDA0003017682660000062
式中，ES表示“太阳常数”，ρE表示地球反照率，O′是dAt的几何中心，φSE是
Figure BDA0003017682660000063
与太阳辐射向量
Figure BDA0003017682660000064
相反方向向量之间的夹角，βSE是
Figure BDA0003017682660000065
与
Figure BDA0003017682660000066
之间的夹角，γSE是dAt法向量
Figure BDA0003017682660000067
与
Figure BDA0003017682660000068
之间的夹角；
(2)构建飞行器表面平行于地球球心截面、太阳辐射垂直飞行器表面的地球反照热流一般表达式；
如图4所示，对于dAt平行于地球球心截面，太阳辐射向量
Figure BDA0003017682660000069
垂直于dAt的情形，给出Γv(Γv∩ΓL＝Γv)对dAt的照度表达式为：
Figure BDA00030176826600000610
式中忽略dAt在地球上投影的影响；
(3)构建地球反照热流解析表达式；
根据
Figure BDA00030176826600000611
与地球半径RE之间的关系，将照度表达式改写为：
Figure BDA00030176826600000612
因为：
Figure BDA00030176826600000613
式中
Figure BDA00030176826600000614
为dA1在地心直角坐标系投影的方位角；
作以下推导：
Figure BDA0003017682660000071
式中O表示地球球心，H表示O与O′之间的距离；
将
Figure BDA0003017682660000072
写为lUO′，则：
Figure BDA0003017682660000073
cosβSE可表示为：
Figure BDA0003017682660000074
又因为：
Figure BDA0003017682660000075
则ESE可作下列推导：
Figure BDA0003017682660000076
令ζ＝cosφSE，则：
Figure BDA0003017682660000077
令REH＝a，RE 2+H2＝b，则上式变为：
Figure BDA0003017682660000081
令
Figure BDA0003017682660000082
则：
Figure BDA0003017682660000083
Figure BDA0003017682660000084
Figure BDA0003017682660000085
(4)计算地球反照热流密度精确值和数值积分值；
地球反照热流密度精确值的计算方法为：
设定H值，将ES、ρE、RE代入地球反照热流解析表达式，即可计算H值高度下的地球反照热流密度精确值ESE-A。
数值积分值的计算方法为：
将积分式中的积分区域Γv∩ΓL以类似地球经纬线的方式进行网格划分，其中经向m份、纬向n份，共得到m×n份，则数值积分值ESE-B为：
Figure BDA0003017682660000086
式中下标i表示第i个dA1以及其相关量；当飞行器表面平行于地球球心截面、太阳辐射垂直飞行器表面时，cosφSEi、cosβSEi、cosγSEi均大于0。
Figure BDA0003017682660000091
根据几何原理计算。dA1i实际为曲面，可按梯形近似计算。当m×n值越大，积分区域Γv∩ΓL划分越密，数值积分值ESE-B精度越高。
(5)误差分析与数值网格划分数选择。
地球反照热流密度精确值ESE-A与数值积分值ESE-B间的相对误差RESE为：
Figure BDA0003017682660000092
当RESE越接近于1，数值积分精度越高；根据相对误差要求，选择网格划分值m×n。
表1是本发明实施例中飞行器表面距地面高度200km，数值积分地球反照热流密度值与解析法精确值(399.2)的相对误差表。由表1可知随着地球表面网格划分数的增多，数值计算精确度逐步提高，但计算量增大，计算时间增多，在实际近地空间热分析中可依据此表格按照计算精度需要选择网格划分数。所以，利用图1技术路线图我们可以得到一种近地空间飞行器反照地球热流的解析方法。
表1 dAt上的地球反照热流密度的数值积分值(W/m2)
(距地面高度200km，解析法精确值为399.2)
Figure BDA0003017682660000093
以上示意性的对本发明及其实施方式进行了描述，该描述没有限制性，附图中所示的也只是本发明的实施方式之一，实际的结构并不局限于此。所以，如果本领域的普通技术人员受其启示，在不脱离本发明创造宗旨的情况下，不经创造性的设计出与该技术方案相似的结构方式及实施例，均应属于本发明的保护范围。

Claims (4)

1.一种近地空间飞行器反照地球热流的解析方法，其特征在于：包括以下步骤：
(1)构建地球反照热流双重积分一般表达式；
将地球视为漫反射规则灰球体，其表面按经纬线的方式近似划分为一系列的平面面元，用Γv表示对近地空间目标可视的地球球冠，ΓL表示被太阳照射的区域，则Γv∩ΓL是地球反照辐射的区域；dA1是地球表面上的一个微面元，U是dA1的几何中心，
Figure FDA0003017682650000011
是dA1的外向单位法矢量，则飞行器某表面平面dAt上的地球反照热流密度ESE可表示为：
Figure FDA0003017682650000012
式中，ES表示“太阳常数”，ρE表示地球反照率，O′是dAt的几何中心，φSE是
Figure FDA0003017682650000013
与太阳辐射向量
Figure FDA0003017682650000014
相反方向向量之间的夹角，βSE是
Figure FDA0003017682650000015
与
Figure FDA0003017682650000016
之间的夹角，γSE是dAt法向量
Figure FDA0003017682650000017
与
Figure FDA0003017682650000018
之间的夹角；
(2)构建飞行器表面平行于地球球心截面、太阳辐射垂直飞行器表面的地球反照热流一般表达式；
对于dAt平行于地球球心截面，太阳辐射向量
Figure FDA00030176826500000112
垂直于dAt的情形，给出Γv对dAt的照度表达式为：
Figure FDA0003017682650000019
式中忽略dAt在地球上投影的影响；
(3)构建地球反照热流解析表达式；
根据
Figure FDA00030176826500000110
与地球半径RE之间的关系，将照度表达式改写为：
Figure FDA00030176826500000111
因为：
Figure FDA0003017682650000021
式中
Figure FDA0003017682650000022
为dA1在地心直角坐标系投影的方位角；
作以下推导：
Figure FDA0003017682650000023
式中O表示地球球心，H表示O与O′之间的距离；
将
Figure FDA0003017682650000024
写为lUO′，则：
Figure FDA0003017682650000025
cosβSE可表示为：
Figure FDA0003017682650000026
又因为：
Figure FDA0003017682650000027
则ESE可作下列推导：
Figure FDA0003017682650000028
令ζ＝cosφSE，则：
Figure FDA0003017682650000031
令REH＝a，RE 2+H2＝b，则上式变为：
Figure FDA0003017682650000032
令
Figure FDA0003017682650000033
则：
Figure FDA0003017682650000034
Figure FDA0003017682650000035
Figure FDA0003017682650000036
(4)计算地球反照热流密度精确值和数值积分值；
(5)误差分析与数值网格划分数选择。
2.根据权利要求1所述的一种近地空间飞行器反照地球热流的解析方法，其特征在于：步骤(4)中，地球反照热流密度精确值的计算方法为：
设定H值，将ES、ρE、RE代入地球反照热流解析表达式，即可计算H值高度下的地球反照热流密度精确值ESE-A。
3.根据权利要求2所述的一种近地空间飞行器反照地球热流的解析方法，其特征在于：步骤(4)中，数值积分值的计算方法为：
将积分式中的积分区域Γv∩ΓL以类似地球经纬线的方式进行网格划分，其中经向m份、纬向n份，共得到m×n份，则数值积分值ESE-B为：
Figure FDA0003017682650000041
式中下标i表示第i个dA1以及其相关量；当m×n值越大，积分区域Γv∩ΓL划分越密，数值积分值ESE-B精度越高。
4.根据权利要求3所述的一种近地空间飞行器反照地球热流的解析方法，其特征在于：步骤(5)中，地球反照热流密度精确值ESE-A与数值积分值ESE-B间的相对误差RESE为：
Figure FDA0003017682650000042
当RESE越接近于1，数值积分精度越高；根据相对误差要求，选择网格划分值m×n。
CN202110393505.9A 2021-04-13 2021-04-13 一种近地空间飞行器反照地球热流的解析方法 Active CN113127795B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110393505.9A CN113127795B (zh) 2021-04-13 2021-04-13 一种近地空间飞行器反照地球热流的解析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110393505.9A CN113127795B (zh) 2021-04-13 2021-04-13 一种近地空间飞行器反照地球热流的解析方法

Publications (2)

Publication Number Publication Date
CN113127795A true CN113127795A (zh) 2021-07-16
CN113127795B CN113127795B (zh) 2023-05-23

Family

ID=76776088

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110393505.9A Active CN113127795B (zh) 2021-04-13 2021-04-13 一种近地空间飞行器反照地球热流的解析方法

Country Status (1)

Country Link
CN (1) CN113127795B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130226537A1 (en) * 2011-08-30 2013-08-29 Joon-Bum JEE Method and apparatus for providing one layer solar radiation model for calculation of insolation
CN104216864A (zh) * 2014-08-22 2014-12-17 航天东方红卫星有限公司 一种立方星的热设计方法
CN108827279A (zh) * 2018-04-11 2018-11-16 西北工业大学 一种基于地球反照的模拟式太阳敏感器抗干扰方法
CN109426654A (zh) * 2017-08-21 2019-03-05 哈尔滨工业大学 复杂结构变姿态卫星空间外热流获取方法
WO2020125839A1 (de) * 2018-12-18 2020-06-25 GRID INVENT gGmbH Elektronisches element und elektrisch angesteuertes anzeigeelement

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130226537A1 (en) * 2011-08-30 2013-08-29 Joon-Bum JEE Method and apparatus for providing one layer solar radiation model for calculation of insolation
CN104216864A (zh) * 2014-08-22 2014-12-17 航天东方红卫星有限公司 一种立方星的热设计方法
CN109426654A (zh) * 2017-08-21 2019-03-05 哈尔滨工业大学 复杂结构变姿态卫星空间外热流获取方法
CN108827279A (zh) * 2018-04-11 2018-11-16 西北工业大学 一种基于地球反照的模拟式太阳敏感器抗干扰方法
WO2020125839A1 (de) * 2018-12-18 2020-06-25 GRID INVENT gGmbH Elektronisches element und elektrisch angesteuertes anzeigeelement

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
邱义芬,陈拥华,袁修干,杨春信: "舱外航天服的轨道空间外热流计算方法", 航空学报 *

Also Published As

Publication number Publication date
CN113127795B (zh) 2023-05-23

Similar Documents

Publication Publication Date Title
CN104573251B (zh) 一种星载光学遥感器全视场表观光谱辐亮度确定方法
CN103616028B (zh) 一种基于单星敏感器的星光折射卫星自主导航方法
CN105160125B (zh) 一种星敏感器四元数的仿真分析方法
CN102261921B (zh) 一种修正大气折射对星敏感器精度影响的方法
CN101082497A (zh) 一种天体敏感器测量基准转换方法及装置
CN104406607A (zh) 一种多视场复合光学敏感器的标定装置及方法
CN104501835A (zh) 一种面向空间应用异构imu初始对准的地面试验系统及方法
CN105138005A (zh) 一种基于星间距离的相对轨道要素确定方法
CN105444778A (zh) 一种基于成像几何反演的星敏感器在轨定姿误差获取方法
CN105023281B (zh) 基于点扩散函数波前修正的星点像质心计算方法
Mitchell et al. A non-intrusive optical (NIO) approach to characterize heliostats in utility-scale power tower plants: Methodology and in-situ validation
CN109708667B (zh) 一种基于激光陀螺的双动态目标跟踪引导方法
CN115343744A (zh) 空中运动目标的光学单双星联合星上定位方法及系统
Linares et al. Photometric data from non-resolved objects for space object characterization and improved atmospheric modeling
CN112857306B (zh) 一种视频卫星任意视向点的连续太阳高度角确定方法
Dubyagin et al. Conditions of loss cone filling by scattering on the curved field lines for 30 keV protons during geomagnetic storm as inferred from numerical trajectory tracing
CN108072387A (zh) 一种低精度敏感器在轨偏差的地面校正方法及系统
CN111026139B (zh) 一种基于飞行轨迹的三维模型姿态调整控制方法
CN113127795A (zh) 一种近地空间飞行器反照地球热流的解析方法
US9182584B2 (en) Method and system for evaluating stare-time by a pointing system
Putin et al. Study of microaccelerations onboard the International Space Station with the DAKON-M convection sensor
Ziebart et al. GPS block IIR non-conservative force modeling: computation and implications
Hirth et al. Validation of models for the calculation of sun positions and mapped radiation on inclined surfaces
Harman et al. An Algorithm for Converting Static Earth Sensor Measurements Into Earth Observation Vectors
CN105486304B (zh) 太阳同步轨道星敏感器遇月亮的预报方法

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