CN109446559A - 一种浮空器辐射热环境特性估算方法 - Google Patents
一种浮空器辐射热环境特性估算方法 Download PDFInfo
- Publication number
- CN109446559A CN109446559A CN201811097912.XA CN201811097912A CN109446559A CN 109446559 A CN109446559 A CN 109446559A CN 201811097912 A CN201811097912 A CN 201811097912A CN 109446559 A CN109446559 A CN 109446559A
- Authority
- CN
- China
- Prior art keywords
- atm
- ground
- atmospheric transmittance
- height
- aerostatics
- 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.)
- Pending
Links
- 230000005855 radiation Effects 0.000 title claims abstract description 44
- 230000007613 environmental effect Effects 0.000 title claims abstract description 16
- 238000011156 evaluation Methods 0.000 title claims abstract description 11
- 238000002834 transmittance Methods 0.000 claims abstract description 36
- 238000000034 method Methods 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 claims description 5
- 238000009738 saturating Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 5
- 238000012360 testing method Methods 0.000 abstract description 2
- 238000013461 design Methods 0.000 abstract 1
- 230000035699 permeability Effects 0.000 description 3
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 1
- CBENFWSGALASAD-UHFFFAOYSA-N Ozone Chemical compound [O-][O+]=O CBENFWSGALASAD-UHFFFAOYSA-N 0.000 description 1
- 239000000443 aerosol Substances 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明属于浮空器热环境分析领域，尤其涉及一种浮空器辐射热环境特性估算方法。包括步骤计算大气层外边界辐射强度Isun；计算地面处太阳直射辐射强度Id，0；计算地面处大气透过率τatm，0，τatm，0＝Id，0/Isun；根据地面处大气透过率τatm，0计算高度H处大气透过率τatm；根据H处大气透过率τatm、大气层外边界辐射强度Isun确定高度H处太阳直射辐射强度Id，Id＝Isun·τatm；计算高度H处地面长波辐射大气透过率τIR,atm；确定地表温度Tground，并根据地面长波辐射透过率和发射率及玻尔兹曼定律确定地面长波辐射强度。本发明根据较少的测量数据估算目标场地放飞的浮空器在不同高度H处所受辐射热环境特性参数，减小了高度和地理环境给浮空器热特性计算带来的误差，在浮空器热设计、蒙皮选材、飞行试验安全等方面具有指导意义。
Description
技术领域
本发明属于浮空器热环境分析领域，尤其涉及一种浮空器辐射热环境特性估算方法。
背景技术
浮空器是一种依靠浮升气体升空、体积巨大的低动态飞行器，在通信中继、对地观测、空间探测等领域具有广阔应用前景。浮空器升降及驻空过程中，外界热环境变化导致内部浮升气体温度变化，会影响到其内部压力、飞行控制及安全。浮空器的环境热效应主要受对流和热辐射综合影响，其中热辐射主要受太阳直射辐射和地面长波辐射的影响。不同地区热环境迥异，尤其是辐射热环境存在较大差异，以往研究中多采用通用经验模型进行估算，难以反映地区热辐射环境差异给浮空器热特性带来的影响。戴秋敏(浮空器热环境与热特性研究，南京航空航天大学博士学位论文，2014年)根据气溶胶、水汽、臭氧、二氧化碳等参数的垂直分布情况拟合了新的太阳直射辐射强度和地面长波辐射强度模型，计算精度得到了提高，但该方法测量数据种类较多，对于一些测量条件有限的区域难以实现。
发明内容
本发明目的在于，提供一种浮空器辐射热环境特性估算方法，可根据较少的测量数据估算目标场地放飞的浮空器在不同高度处所受辐射热环境参数。
为达到上述发明目的，本发明所述浮空器辐射热环境特性估算方法，包括如下步骤：
S1，计算大气层外边界辐射强度Isun；
S2，计算地面处太阳直射辐射强度Id，0；
S3，计算地面处大气透过率τatm，0，τatm，0＝Id，0/Isun；
S4，根据地面处大气透过率τatm，0计算高度H处大气透过率τatm；
S41，高度H处大气透过率τatm与大气质量m呈指数关系，τatm＝0.5[exp(-am)+exp(-bm)]；其中，a，b均为待拟合系数；
S42，将高度0处大气透过率，即τatm地面处大气透过率τatm，0带入公式τatm＝0.5[exp(-am)+exp(-bm)]，利用最小二乘法对地面处大气透过率τatm，0与大气质量m进行拟合，确定系数a，b的值；
S34，将确定的系数a，b带入公式τatm＝0.5[exp(-am)+exp(-bm)]，得到高度H处大气透过率τatm＝0.5[exp(-am)+exp(-bm)]；
S5，根据H处大气透过率τatm、大气层外边界辐射强度Isun确定高度H处太阳直射辐射强度Id，Id＝Isun·τatm；
S6，计算高度H处地面长波辐射大气透过率τIR,atm；
式中，P0为地面处大气压强，Patm为高度H处大气压强；
S7，确定地表温度Tground，并根据地面长波辐射透过率和发射率及玻尔兹曼定律确定地面长波辐射强度。
S71，地表温度根据气温Tatm确定，Tground＝cTatm+d；
其中，c，d均为拟合系数；
S72，根据地表温度Tground、地面长波发射率εground和高度H处地面长波辐射大气透过率τIR,atm计算高度H处地面长波辐射强度IIR,ground，
式中，σ为斯蒂芬-玻尔兹曼常数，其值为σ＝5.67×10-8。
进一步，步骤S41中根据高度H处大气压强Patm以及所对应的太阳高度角h计算大气质量m，大气质量m计算公式为：
进一步，所述不同高度H的范围为地面至海拔32km处。
本发明根据较少的测量数据估算目标场地放飞的浮空器在不同高度H处所受辐射热环境特性参数，减小了高度和地理环境给浮空器热特性计算带来的误差，在浮空器热设计、蒙皮选材、飞行试验安全等方面具有指导意义。
附图说明
图1为本发明实施例中一种浮空器辐射热环境特性估算方法的流程示意图。
具体实施方式
下面结合附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例、基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例中所采用的太阳直射辐射表等设备均为市售。
实施例1
一种浮空器辐射热环境特性估算方法，如图1所示，包括如下步骤：
S1，计算大气层外边界辐射强度Isun；
根据不同日期的天序号n确定大气层外边界辐射强度Isun，大气层外边界辐射强度Isun计算公式为
其中I0为太阳常数，其值为1367W/m2。
S2，利用太阳直射辐射表计算地面处太阳直射辐射强度Id，0；
S3，根据大气层外边界辐射强度Isun及地面处太阳直射辐射强度Id，0计算地面处大气透过率τatm，0，τatm，0＝Id，0/Isun；
S4，根据地面处大气透过率τatm，0计算高度H处大气透过率τatm；
S41，不同高度处大气压不同，从而大气质量不同，先根据测量时刻的高度H处大气压强Patm以及所对应的太阳高度角h计算大气质量m，大气质量m计算公式为：
高度H处大气透过率τatm与大气质量m呈指数关系，τatm＝0.5[exp(-am)+exp(-bm)]；其中，a，b均为待拟合系数；
S42，将高度0处大气透过率，即τatm地面处大气透过率τatm，0带入公式τatm＝0.5[exp(-am)+exp(-bm)]，利用最小二乘法对地面处大气透过率τatm，0与大气质量m进行拟合，确定系数a，b的值；
S34，将确定的系数a，b带入公式τatm＝0.5[exp(-am)+exp(-bm)]，得到高度H处大气透过率τatm＝0.5[exp(-am)+exp(-bm)]；
所述不同高度H的范围为地面至海拔32km处。
S5，根据H处大气透过率τatm、大气层外边界辐射强度Isun确定高度H处太阳直射辐射强度Id，Id＝Isun·τatm；
S6，计算高度H处地面长波辐射大气透过率τIR,atm；
地面长波辐射透过率在地面处为1，透过率随海拔上升呈指数递减，与太阳直射辐射的大气透过率形式类似，故高度H处地面长波辐射大气透过率τIR,atm采用下式计算：
式中，系数a，b的值与步骤S4中确定的值相同，P0为地面处大气压强，Patm为高度H处大气压强；
优选地，所在高度H处的大气压强Patm根据下式确定：
S7，确定地表温度Tground，并根据H处地面长波辐射大气透过率τIR,atm和地面长波发射率εground及玻尔兹曼定律确定地面长波辐射强度；
S71，地表温度根据气温Tatm确定，Tground＝cTatm+d；
其中，c，d均为拟合系数；
S72，根据地表温度Tground、地面长波发射率εground和高度H处地面长波辐射大气透过率τIR,atm计算高度H处地面长波辐射强度IIR,ground，
式中，σ为斯蒂芬-玻尔兹曼常数，其值为σ＝5.67×10-8。
优选地，地面长波发射率εground根据不同地表类型确定：草地0.98，土壤0.93，城市0.92，冰雪0.99，水体0.95，沙漠0.88。
综上可知，本发明中的浮空器辐射热环境特性估算方法实现了利用较少的测量参数对目标地区放飞的浮空器在不同高度处受到的太阳直射辐射强度和地面长波辐射强度进行估算的效果。
Claims (3)
1.一种浮空器辐射热环境特性估算方法，其特征在于，包括如下步骤：
S1，计算大气层外边界辐射强度Isun；
S2，计算地面处太阳直射辐射强度Id，0；
S3，计算地面处大气透过率τatm，0，τatm，0＝Id，0/Isun；
S4，根据地面处大气透过率τatm，0计算高度H处大气透过率τatm；
S41，高度H处大气透过率τatm与大气质量m呈指数关系，τatm＝0.5[exp(-am)+exp(-bm)]；其中，a，b均为待拟合系数；
S42，将高度0处大气透过率，即τatm地面处大气透过率τatm，0带入公式τatm＝0.5[exp(-am)+ex(p-bm)]，利用最小二乘法对地面处大气透过率τatm，0与大气质量m进行拟合，确定系数a，b的值；
S34，将确定的系数a，b带入公式τatm＝0.5[exp(-am)+exp(-bm)]，得到高度H处大气透过率τatm＝0.5[exp(-am)+exp(-bm)]；
S5，根据H处大气透过率τatm、大气层外边界辐射强度Isun确定高度H处太阳直射辐射强度Id，Id＝Isun·τatm；
S6，计算高度H处地面长波辐射大气透过率τIR,atm；
式中，P0为地面处大气压强，Patm为高度H处大气压强；
S7，确定地表温度Tground，并根据H处地面长波辐射大气透过率τIR,atm和地面长波发射率εground及玻尔兹曼定律确定地面长波辐射强度；
S71，地表温度根据气温Tatm确定，Tground＝cTatm+d；
其中，c，d均为拟合系数；
S72，根据地表温度Tground、地面长波发射率εground和高度H处地面长波辐射大气透过率τIR,atm计算高度H处地面长波辐射强度IIR,ground，
式中，σ为斯蒂芬-玻尔兹曼常数，其值为σ＝5.67×10-8。
2.根据权利要求1所述一种浮空器辐射热环境特性估算方法，其特征在于，步骤S41中根据高度H处大气压强Patm以及所对应的太阳高度角h计算大气质量m，大气质量m计算公式为：
3.根据权利要求1所述一种浮空器辐射热环境特性估算方法，其特征在于，不同高度H的范围为地面至海拔32km处。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811097912.XA CN109446559A (zh) | 2018-09-20 | 2018-09-20 | 一种浮空器辐射热环境特性估算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811097912.XA CN109446559A (zh) | 2018-09-20 | 2018-09-20 | 一种浮空器辐射热环境特性估算方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109446559A true CN109446559A (zh) | 2019-03-08 |
Family
ID=65533111
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811097912.XA Pending CN109446559A (zh) | 2018-09-20 | 2018-09-20 | 一种浮空器辐射热环境特性估算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109446559A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110816881A (zh) * | 2019-11-01 | 2020-02-21 | 中国科学院光电研究院 | 一种浮空器热特性试验方法 |
| CN110816880A (zh) * | 2019-11-01 | 2020-02-21 | 中国科学院光电研究院 | 一种浮空器热特性试验系统 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4614429A (en) * | 1984-05-14 | 1986-09-30 | The United States Of America As Represented By The Secretary Of The Air Force | Atmospheric contrast transmittance monitor |
| CN101976275A (zh) * | 2010-09-21 | 2011-02-16 | 北京航空航天大学 | 飞机红外辐射与大气透过率建模方法 |
| US20120310448A1 (en) * | 2011-06-03 | 2012-12-06 | Sikorsky Aircraft Corporation | Aircraft environmental threat detection and mitigation |
| CN105243289A (zh) * | 2015-11-17 | 2016-01-13 | 上海无线电设备研究所 | 复杂目标红外辐射特性一体化建模方法 |
| CN106706133A (zh) * | 2016-12-31 | 2017-05-24 | 华中科技大学 | 一种点斑状目标姿态估计方法及系统 |
| CN106772682A (zh) * | 2016-12-31 | 2017-05-31 | 华中科技大学 | 一种动目标的红外辐射光谱特性仿真分析方法 |
| CN106770062A (zh) * | 2016-12-16 | 2017-05-31 | 中国科学院合肥物质科学研究院 | 一种大气透过率测量及标定方法 |
-
2018
- 2018-09-20 CN CN201811097912.XA patent/CN109446559A/zh active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4614429A (en) * | 1984-05-14 | 1986-09-30 | The United States Of America As Represented By The Secretary Of The Air Force | Atmospheric contrast transmittance monitor |
| CN101976275A (zh) * | 2010-09-21 | 2011-02-16 | 北京航空航天大学 | 飞机红外辐射与大气透过率建模方法 |
| US20120310448A1 (en) * | 2011-06-03 | 2012-12-06 | Sikorsky Aircraft Corporation | Aircraft environmental threat detection and mitigation |
| CN105243289A (zh) * | 2015-11-17 | 2016-01-13 | 上海无线电设备研究所 | 复杂目标红外辐射特性一体化建模方法 |
| CN106770062A (zh) * | 2016-12-16 | 2017-05-31 | 中国科学院合肥物质科学研究院 | 一种大气透过率测量及标定方法 |
| CN106706133A (zh) * | 2016-12-31 | 2017-05-24 | 华中科技大学 | 一种点斑状目标姿态估计方法及系统 |
| CN106772682A (zh) * | 2016-12-31 | 2017-05-31 | 华中科技大学 | 一种动目标的红外辐射光谱特性仿真分析方法 |
Non-Patent Citations (5)
| Title |
|---|
| QIUMIN DAI: "A simple model to predict solar radiation under clear sky condition", 《ADVANCES IN SPACE REASEARCH》 * |
| RODGER E.FARLEY: "BalloonAscent: 3-D Simulation Tool for the Ascent and Float of High-Altitude Balloons", 《HTTPS://DOI.ORG/10.2514/6.2005-7412》 * |
| 夏新林: "平流层浮空器的热特性与研究现状", 《航空学报》 * |
| 姜会飞: "地面温度与气温的统计分析", 《中国农业气象》 * |
| 字贵才: "超压气球驻空过程热力学特性分析", 《航天返回与遥感》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110816881A (zh) * | 2019-11-01 | 2020-02-21 | 中国科学院光电研究院 | 一种浮空器热特性试验方法 |
| CN110816880A (zh) * | 2019-11-01 | 2020-02-21 | 中国科学院光电研究院 | 一种浮空器热特性试验系统 |
| CN110816881B (zh) * | 2019-11-01 | 2022-11-15 | 中国科学院空天信息创新研究院 | 一种浮空器热特性试验方法 |
| CN110816880B (zh) * | 2019-11-01 | 2022-11-15 | 中国科学院空天信息创新研究院 | 一种浮空器热特性试验系统 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106169014B (zh) | 基于遥感数据的区域地表感热/潜热通量反演方法及系统 | |
| Vining et al. | Estimation of sensible heat flux from remotely sensed canopy temperatures | |
| Wang et al. | Assessment of uncertainties in eddy covariance flux measurement based on intensive flux matrix of HiWATER-MUSOEXE | |
| CN106932101B (zh) | Hj-1b卫星中红外和热红外通道联合的地表温度反演方法 | |
| Liu et al. | Full-scale measurements of convective coefficient on external surface of a low-rise building in sheltered conditions | |
| CN109580003A (zh) | 一种静止气象卫星热红外数据估算近地面大气温度方法 | |
| CN103438900B (zh) | 三线阵相机影像协同绝对辐射定标和校正方法 | |
| CN103810387A (zh) | 基于modis数据的地表蒸散发全遥感反演方法及系统 | |
| CN102176002A (zh) | 基于地表水热通量遥感反演的干旱监测方法及系统 | |
| Liu et al. | An improved physical split-window algorithm for precipitable water vapor retrieval exploiting the water vapor channel observations | |
| CN108020322A (zh) | 煤田火区的机载热红外高光谱遥感定量探测方法 | |
| CN109446559A (zh) | 一种浮空器辐射热环境特性估算方法 | |
| CN109959970B (zh) | 一种天空半球热红外大气下行辐射地面测量方法 | |
| Mahrt et al. | Bulk formulation of the surface heat flux | |
| CN108827963A (zh) | 基于多尺度热惯量模型反演的土壤表层含水量测算方法 | |
| CN110968955A (zh) | 一种蒸发比观测的时空拓展技术 | |
| Mitnik et al. | Microwave radiometer MTVZA-GY on new Russian satellite Meteor-M No. 2-2 and sudden stratospheric warming over Antarctica | |
| Smith et al. | Remote sensing of land surface temperature: the directional viewing effect | |
| Iwasaki | Relating lightning features and topography over the Tibetan Plateau using the World Wide Lightning Location Network data | |
| CN113343435A (zh) | 一种适应fy4a卫星上agri仪器的射出长波辐射计算方法 | |
| Zhou et al. | Land surface albedo estimation with Chinese GF-1 WFV data in Northwest China | |
| Er | A Study of Importance of Solar Calculations for" Two Colored Rainbow System" in Istanbul | |
| Seto et al. | Radiative characteristics at 89 and 36 GHz for satellite-based cloud water estimation over land | |
| Adeniyi et al. | Estimation of surface energy fluxes from bare ground in a tropical station using Priestley Taylor method | |
| Shan et al. | Calculation method of infrared temperature on the natural ground surface |
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 | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20211122 Address after: 471000 sub box F11, PO Box 086, Luoyang City, Henan Province Applicant after: Chinese People's Liberation Army 63660 Address before: 841700 box 190, box 21, Urumqi, Xinjiang Uygur Autonomous Region Applicant before: 63653 FORCES, PLA |
|
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190308 |