CN103471723A - Novel method for forecasting daytime temperature and night temperature of stratosphere balloon - Google Patents

Novel method for forecasting daytime temperature and night temperature of stratosphere balloon Download PDF

Info

Publication number
CN103471723A
CN103471723A CN2013104066179A CN201310406617A CN103471723A CN 103471723 A CN103471723 A CN 103471723A CN 2013104066179 A CN2013104066179 A CN 2013104066179A CN 201310406617 A CN201310406617 A CN 201310406617A CN 103471723 A CN103471723 A CN 103471723A
Authority
CN
China
Prior art keywords
stratosphere balloon
balloon
stratosphere
formula
alpha
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
CN2013104066179A
Other languages
Chinese (zh)
Other versions
CN103471723B (en
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.)
Beihang University
Original Assignee
Beihang University
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 Beihang University filed Critical Beihang University
Priority to CN201310406617.9A priority Critical patent/CN103471723B/en
Publication of CN103471723A publication Critical patent/CN103471723A/en
Application granted granted Critical
Publication of CN103471723B publication Critical patent/CN103471723B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)

Abstract

A novel method for forecasting the daytime temperature and the night temperature of a stratosphere balloon comprises the steps of 1 providing assumption conditions of the novel method for forecasting the daytime temperature and the night temperature of the stratosphere balloon; 2 according to the use environment of the stratosphere balloon and basic knowledge of the engineering thermodynamics, determining the energy exchanging mechanism of the stratosphere balloon and all energy exchanging items; 3 establishing heat analysis models of the energy exchanging items respectively to obtain specific parameterized expression forms of all the energy exchanging items; 4 according to the principle of energy balance, an energy balance equation of the stratosphere balloon in daytime and an energy balance equation of the stratosphere balloon at night are established respectively; 5 substituting all the energy exchanging items obtained in the step 3 into the energy balance equations in the step 4, and then a simple iteration method in a numerical solution method is used for obtaining the daytime temperature of the stratosphere balloon and the night temperature of the stratosphere balloon. According to the novel method forecasting the daytime temperature and the night temperature of the stratosphere balloon, only bag cloth material parameters, the suspension height and the suspension time are needed, and then the daytime temperature and the night temperature of the stratosphere balloon can be obtained.

Description

A kind of new method of predicting the day and night temperature of stratosphere balloon
Technical field
The invention provides a kind of new method of predicting the day and night temperature of stratosphere balloon, belong to the aeronautical technology field.
Background technology
Stratosphere balloon is a kind of low cost, high efficiency aircraft, can carry out for a long time reconnaissance mission, carry out goods input etc., has boundless application potential and prospect, and therefore, correlation technique becomes the focus of studying this year.The temperature of stratosphere balloon when suspending is one of very crucial input parameter for the designer, and especially most important to choosing of its capsule cloth material, therefore, how to obtain effective temperature is all the problem that the engineering staff need to solve all the time.Can directly measure its temperature by the stratosphere balloon flight test, but experimentation cost is higher, and the cycle of wasting time and energy is long, also is difficult to all situations are all measured by test; Numerical simulation method need to be set up complicated finite element model, calculation of complex, and counting yield is low; More existing analytic methods, the factor that model is considered is comprehensive not, and theory solves also more complicated.Therefore, the present invention proposes a kind of analytic method of predicting the day and night temperature of stratosphere balloon, the method is very simple and practical, the capsule cloth material parameter, hoverheight and the time that only need a small amount of stratosphere balloon, just can be easy to obtain the day and night temperature of stratosphere balloon, visible the present invention has Important Academic meaning and engineering using value.
Summary of the invention
The invention provides a kind of new method of predicting the day and night temperature of stratosphere balloon, it is easy that the method has calculating, the precision advantages of higher, and its technical scheme is as follows:
A kind of new method of predicting the day and night temperature of stratosphere balloon of the present invention, the method concrete steps are as follows:
The assumed condition of the new method of the day and night temperature of step 1, proposition prediction stratosphere balloon.
Assumed condition comprises:
(1), when the stratosphere balloon stand under load, the actual similar water-drop-shaped of shape, for simplified model, be reduced to by stratosphere balloon the circle that radius is R.The characteristic dimension L of stratosphere balloon can mean with diameter;
(2) the capsule cloth due to stratosphere balloon is very thin, in micron dimension, therefore, supposes that the temperature of capsule cloth thickness direction is identical, does not consider the heat transmission between capsule cloth internal layer and skin;
(3) do not consider the loss of the inner helium of stratosphere balloon, the weight of supposing stratosphere balloon is constant;
(4) due to earth radius R efor 6371.23km, and the sun apart from the earth the distance be 1.5 * 10 8km, therefore, suppose that when sunshine arrives earth overhead be directional light;
(5) suppose that the earth is the thermal equilibrium body;
(6) do not consider the variation of intensity of solar radiation with solar radiation angle and atmospheric extinction coefficient;
(7) do not consider the impact of carbon dioxide and ozone infrared radiation.
Step 2, according to the environment for use of stratosphere balloon and the ABC of Engineering Thermodynamics, determine the energy exchange mechanism of stratosphere balloon and each energy exchange item.
Stratosphere balloon mainly carries out energy exchange by radiation and convection current, comprising: the direct solar radiation item Q that the stratosphere balloon episphere absorbs 1(t) sun scattered radiation item Q that, stratosphere balloon absorbs 2(t) earth reflected radiation item Q that, the stratosphere balloon lower semisphere absorbs 3(t) earth infrared radiation item Q that, the stratosphere balloon lower semisphere absorbs 4(t), the external infrared radiation item Q of stratosphere balloon 5(t) the internal reflection infrared radiation item Q that, stratosphere balloon absorbs 6(t), the convection heat transfer item Q between stratosphere balloon and outside atmosphere 7and the convection heat transfer item Q between stratosphere balloon and internal gas (t) 8(t).
Step 3, each energy exchange item proposed for step 2 are set up respectively its thermal model, obtain the concrete concrete Parameter Expression form of each energy exchange item.
The area of the infinitesimal of stratosphere balloon can be expressed as
In formula, dA is the unit dimension on stratosphere balloon, the latitude that θ is spherical coordinate system or solar radiation angle or sunshine and imfinitesimal method between angle,
Figure BDA00003792174900022
longitude for spherical coordinate system.
The direct solar radiation that infinitesimal on the stratosphere balloon episphere absorbs is
dQ 1(t)=α 1τI 0cosθdA (2)
In formula, α 1for the capsule cloth of the stratosphere balloon absorptivity to direct solar radiation, τ is the atmosphere coefficient of transparency, I 0for solar constant.
The sun scattered radiation that infinitesimal on stratosphere balloon absorbs is
dQ 2(t)=α 1I 1dA (3)
In formula, I 1for sun diffuse radiation.
The earth reflected radiation that infinitesimal on the stratosphere balloon lower semisphere absorbs is
d Q 3 ( t ) = α 1 ρ g τ I 0 ( R e R e + H ) 2 cos θ cos φdA - - - ( 4 )
In formula, ρ gfor ground surface reflectance mean value, R efor earth radius, the hoverheight of H stratosphere balloon, φ is the angle between sunshine direction and ground ball center peace fluid layer balloon center line connecting direction.
The earth infrared radiation that infinitesimal on the stratosphere balloon lower semisphere absorbs is
d Q 4 ( t ) = α 2 τ I 0 R e 2 1 - ρ g 4 ( R e + H ) 2 cos θdA - - - ( 5 )
In formula, α 2for the capsule cloth of the stratosphere balloon absorptivity to infrared radiation.
The external infrared radiation of infinitesimal on stratosphere balloon is
d Q 5 ( t ) = 2 ϵσ T 1 4 dA - - - ( 6 )
In formula, the emissivity of the capsule cloth that ε is stratosphere balloon to infrared radiation, σ is Si Difen-Boltzmann constant, T 1capsule Bu Wendu for stratosphere balloon.
The internal reflection infrared radiation that infinitesimal on stratosphere balloon absorbs is
d Q 6 ( t ) = α 2 ϵσ T 1 4 1 - r dA - - - ( 7 )
In formula, the reflectivity that r is stratosphere balloon.
Infinitesimal and the convection heat transfer between outside atmosphere on stratosphere balloon are
dQ 7(t)=h 1(T 1-T 2)dA (8)
In formula, h 1for the convection transfer rate between stratosphere balloon and outside atmosphere, T 2temperature for outside atmosphere.
Convection heat transfer between infinitesimal on stratosphere balloon and inner helium is
dQ 8(t)=h 2(T 1-T 3)dA (9)
In formula, h 2for the convection transfer rate between stratosphere balloon and inner helium, T 3temperature for inner helium.
Sun diffuse radiation I in formula (3) 1can be expressed as
I 1 = 1 2 I 0 cos θ cos 2 θ 2 - - - ( 10 )
Included angle between sunshine direction in formula (4) and ground ball center peace fluid layer balloon center line connecting direction can be expressed as
φ = π | 12 - t | 12 , ( 6 a . m . ≤ t ≤ 18 p . m . ) - - - ( 11 )
In formula, t means the time, and from the 6:00 AM to the afternoon, 18 is the cycle on a daytime.
H in formula (8) and formula (9) 1and h 2can be expressed as respectively
h 1 = N u 1 λ 1 L - - - ( 12 )
h 2 = N u 1 λ 1 L - - - ( 13 )
In formula, Nu 1for the Nusselt number of forced convection, λ 1for the heat-conduction coefficient of forced convection, Nu 2for the Nusselt number of free convection, λ 2heat-conduction coefficient for free convection.
In the heat of stratosphere balloon is analyzed, solar constant I 0for 1367W/m 2, atmosphere coefficient of transparency τ is 0.9, ground surface reflectance mean value ρ gbe that 0.3, Si Difen-Boltzmann constant σ is 5.6696 * 10 -8w(m 2k 4), the capsule cloth of stratosphere balloon equates with absorptivity the emissivity of infrared radiation, i.e. ε=α 2.
Respectively formula (2) to formula (9) is carried out to integration, and the substitution correlation parameter, can obtain the total amount of each energy exchange item:
Q 1 ( t ) = 1230.3 π R 2 α 1 , ( 6 a . m . ≤ t ≤ 18 p . m . ) 0 - - - ( 13 )
Q 2 ( t ) = 1139.2 π α 1 R 2 , ( 6 a . m . ≤ t ≤ 18 p . m . ) 0 - - - ( 14 )
Q 3 ( t ) = 366.8 π R 2 α 1 ( R e R e + H ) 2 cos π | 12 - t | 12 , ( 6 a . m . ≤ t ≤ 18 p . m . ) 0 - - - ( 15 )
Figure BDA00003792174900047
Q 5 ( t ) = ∫ A 2 ϵσ T 1 4 dA = 2 AϵAσ T 1 4 = 45.4 π R 2 α 2 × 10 - 8 T 1 4 - - - ( 17 )
Q 6 ( t ) = ∫ A α 2 ϵσ T 1 4 1 - r dA = A α 2 ϵσ T 1 4 1 - r = 22.7 π R 2 α 2 ϵ × 10 - 8 T 1 4 1 - r - - - ( 19 )
Q 7 ( t ) = ∫ A h 1 ( T 1 - T 2 ) dA = A h 1 ( T 1 - T 2 ) = 4 π R 2 h 1 ( T 1 - T 2 ) - - - ( 20 )
Q 8 ( t ) = ∫ A h 2 ( T 1 - T 3 ) dA = A h 2 ( T 1 - T 3 ) = 4 π R 2 h 2 ( T 1 - T 3 ) - - - ( 21 )
Wherein, the infinitesimal described in step 3 refers to the small unit intercepted on stratosphere balloon, and area is infinitely small.
Step 4, according to principle of energy balance, set up respectively stratosphere balloon by day with the energy-balance equation at night.
For daytime, the gross energy of the capsule cloth of stratosphere balloon is
△Q film(t)=Q 1(t)+Q 2(t)+Q 3(t)+Q 4(t)-Q 5(t)+Q 6(t)-Q 7(t)-Q 8(t) (22)
The gross energy of the inside helium of stratosphere balloon is
△Q gas(t)=Q 8(t) (23)
When the energy exchange of stratosphere balloon reaches equilibrium state, △ Q film(t)=0, △ Q gas(t)=0, substitution formula (22) and formula (23) can obtain capsule cloth and inner helium energy-balance equation by day about stratosphere balloon:
Q 1 ( t ) + Q 2 ( t ) + Q 3 ( t ) + Q 4 ( t ) - Q 5 ( t ) + Q 6 ( t ) - Q 7 ( t ) - Q 8 ( t ) = 0 Q 8 ( t ) = 0 - - - ( 24 )
For night, there is no direct solar radiation, sun scattered radiation and earth reflected radiation, i.e. Q 1(t)=0, Q 2(t)=0, Q 3(t)=0, therefore, the capsule cloth of stratosphere balloon and inner helium at the energy-balance equation at night are
Q 4 ( t ) - Q 5 ( t ) + Q 6 ( t ) - Q 7 ( t ) - Q 8 ( t ) = 0 Q 8 ( t ) = 0 - - - ( 25 )
Energy-balance equation (24), (25) of step 5, each energy exchange item substitution step 4 that step 3 is obtained, the simple process of iteration in the recycling method of value solving can try to achieve stratosphere balloon by day with the temperature at night.
Formula (13) to formula (21) is updated to respectively in formula (24) and formula (25), can obtains
( 22.7 π R 2 α 2 ϵ × 10 - 8 1 - r - 45.4 π R 2 α 2 × 10 - 8 ) T 1 4 - 4 π R 2 h 1 T 1 + [ 366.8 π R 2 α 1 ( R e R e + H ) 2 cos π | 12 - t | 12 + 2369.5 π R 2 α 1 + 214.0 π R 2 α 2 ( R e R e + H ) 2 + 4 π R 2 h 1 T 2 ] = 0 T 1 - T 3 = 0
(26)
( 22.7 π R 2 α 2 ϵ × 10 - 8 1 - r - 45.4 π R 2 α 2 × 10 - 8 ) T 1 4 - 4 π R 2 h 1 T 1 + [ 214.0 π R 2 α 2 ( R e R e + H ) 2 + 4 π R 2 h 1 T 2 ] = 0 T 1 - T 3 = 0
(27)
From (26) and formula (27), when the energy exchange of stratosphere balloon reaches equilibrium state, the capsule Bu Wendu of the stratosphere balloon of day and night equates with inner helium temperature.By the simple process of iteration in method of value solving, equation (26) and (27) are solved respectively, can obtain stratosphere balloon by day with capsule Bu Wendu and the inner helium temperature at night.
A kind of new method of predicting the day and night temperature of stratosphere balloon of the present invention, be characterized in very simple and practical, the capsule cloth material parameter, hoverheight and the time that only need a small amount of stratosphere balloon, just can be easy to obtain the day and night temperature of stratosphere balloon.
The accompanying drawing explanation
True form schematic diagram when Fig. 1 a is the stratosphere balloon stand under load.
Simplification shape schematic diagram when Fig. 1 b is the stratosphere balloon stand under load.
Fig. 2 is stratosphere balloon energy exchange schematic diagram.
Fig. 3 is the stratosphere balloon surface infinitesimal schematic diagram that is subject to solar light irradiation
The solar radiation angle schematic diagram that Fig. 4 is stratosphere balloon surface infinitesimal.
Fig. 5 is the relative position schematic diagram between stratosphere balloon and the earth and the sun.
Fig. 6 is to be the FB(flow block) of the method for the invention.
In figure, symbol description is as follows:
R in Fig. 1 b is the stratosphere balloon radius.
Q in Fig. 2 1(t) the direct solar radiation item absorbed for the stratosphere balloon episphere, Q 2(t) sun scattered radiation item absorbed for stratosphere balloon, Q 3(t) earth reflected radiation item absorbed for the stratosphere balloon lower semisphere, Q 4(t) earth infrared radiation item absorbed for the stratosphere balloon lower semisphere, Q 5(t) be the external infrared radiation item of stratosphere balloon, Q 6(t) the internal reflection infrared radiation item absorbed for stratosphere balloon, Q 7(t) be the convection heat transfer item between stratosphere balloon and outside atmosphere, Q 8(t) be the convection heat transfer item between stratosphere balloon and internal gas.
DA in Fig. 3 is the unit dimension on stratosphere balloon, the normal direction that n is the infinitesimal on stratosphere balloon, the latitude that θ is spherical coordinate system or solar radiation angle or sunshine and imfinitesimal method between angle,
Figure BDA00003792174900071
longitude for spherical coordinate system.
θ in Fig. 4 for for the latitude of spherical coordinate system or solar radiation angle or sunshine and imfinitesimal method between angle;
φ in Fig. 5 is the angle between sunshine direction and ground ball center peace fluid layer balloon center line connecting direction.
Embodiment
The FB(flow block) that Fig. 6 is the method for the invention, the present invention divides five steps to realize, is specially:
The assumed condition of the new method of the day and night temperature of step 1, proposition prediction stratosphere balloon.
Assumed condition comprises:
(1), during the stratosphere balloon stand under load, the actual similar water-drop-shaped of shape (as shown in Figure 1a), for simplified model, be reduced to by stratosphere balloon the circle (as shown in Figure 1 b) that radius is R.The characteristic dimension L of stratosphere balloon can mean with diameter;
(2) the capsule cloth due to stratosphere balloon is very thin, in micron dimension, therefore, supposes that the temperature of capsule cloth thickness direction is identical, does not consider the heat transmission between capsule cloth internal layer and skin;
(3) do not consider the loss of the inner helium of stratosphere balloon, the weight of supposing stratosphere balloon is constant;
(4) due to earth radius R efor 6371.23km, and the sun apart from the earth the distance be 1.5 * 10 8km, therefore, suppose that when sunshine arrives earth overhead be directional light;
(5) suppose that the earth is the thermal equilibrium body;
(6) do not consider the variation of intensity of solar radiation with solar radiation angle and atmospheric extinction coefficient;
(7) do not consider the impact of carbon dioxide and ozone infrared radiation.
Step 2, according to the environment for use of stratosphere balloon and the ABC of Engineering Thermodynamics, be defined as the energy exchange mechanism of stratosphere balloon and each energy exchange item.
Stratosphere balloon mainly carries out energy exchange by radiation and convection current, comprising: the direct solar radiation item Q that the stratosphere balloon episphere absorbs 1(t) sun scattered radiation item Q that, stratosphere balloon absorbs 2(t) earth reflected radiation item Q that, the stratosphere balloon lower semisphere absorbs 3(t) earth infrared radiation item Q that, the stratosphere balloon lower semisphere absorbs 4(t), the external infrared radiation item Q of stratosphere balloon 5(t) the internal reflection infrared radiation item Q that, stratosphere balloon absorbs 6(t), the convection heat transfer item Q between stratosphere balloon and outside atmosphere 7and the convection heat transfer item Q between stratosphere balloon and internal gas (t) 8(t).Provided stratosphere balloon energy exchange schematic diagram in Fig. 2.
Step 3, each energy exchange item proposed for step 1 are set up respectively its thermal model, obtain the concrete concrete Parameter Expression form of each energy exchange item.
In Fig. 2 and Fig. 3, provided the infinitesimal of stratosphere balloon, its area can be expressed as
Figure BDA00003792174900083
in formula, dA is the unit dimension on stratosphere balloon, the latitude that θ is spherical coordinate system or solar radiation angle or sunshine and imfinitesimal method between angle,
Figure BDA00003792174900084
longitude for spherical coordinate system.
The direct solar radiation that infinitesimal on the stratosphere balloon episphere absorbs is
DQ 1(t)=α 1τ I 0in cos θ dA (2) formula, α 1for the capsule cloth of the stratosphere balloon absorptivity to direct solar radiation, τ is the atmosphere coefficient of transparency, I 0for solar constant.
The sun scattered radiation that infinitesimal on stratosphere balloon absorbs is
DQ 2(t)=α 1i 1in dA (3) formula, I 1for sun diffuse radiation.
The earth reflected radiation that infinitesimal on the stratosphere balloon lower semisphere absorbs is
d Q 3 ( t ) = α 1 ρ g τ I 0 ( R e R e + H ) 2 cos θ cos φdA - - - ( 4 )
In formula, ρ gfor ground surface reflectance mean value, R efor earth radius, the hoverheight of H stratosphere balloon, φ is the angle between sunshine direction and ground ball center peace fluid layer balloon center line connecting direction, the position relationship between the earth, stratosphere balloon and sun three is as shown in Figure 5.The solar radiation angle schematic diagram that Fig. 4 is stratosphere balloon surface infinitesimal.
The earth infrared radiation that infinitesimal on the stratosphere balloon lower semisphere absorbs is
d Q 4 ( t ) = α 2 τ I 0 R e 2 1 - ρ g 4 ( R e + H ) 2 cos θdA - - - ( 5 )
In formula, α 2for the capsule cloth of the stratosphere balloon absorptivity to infrared radiation.
The external infrared radiation of infinitesimal on stratosphere balloon is
d Q 5 ( t ) = 2 ϵσ T 1 4 dA - - - ( 6 )
In formula, the emissivity of the capsule cloth that ε is stratosphere balloon to infrared radiation, σ is Si Difen-Boltzmann constant, T 1capsule Bu Wendu for stratosphere balloon.
The internal reflection infrared radiation that infinitesimal on stratosphere balloon absorbs is
d Q 6 ( t ) = α 2 ϵσ T 1 4 1 - r dA - - - ( 7 )
In formula, the reflectivity that r is stratosphere balloon.
Infinitesimal and the convection heat transfer between outside atmosphere on stratosphere balloon are
dQ 7(t)=h 1(T 1-T 2)dA (8)
In formula, h 1for the convection transfer rate between stratosphere balloon and outside atmosphere, T 2temperature for outside atmosphere.
Convection heat transfer between infinitesimal on stratosphere balloon and inner helium is
dQ 8(t)=h 2(T 1-T 3)dA (9)
In formula, h 2for the convection transfer rate between stratosphere balloon and inner helium, T 3temperature for inner helium.
Sun diffuse radiation I in formula (3) 1can be expressed as
I 1 = 1 2 I 0 cos θ cos 2 θ 2 - - - ( 10 )
Included angle between sunshine direction in formula (4) and ground ball center peace fluid layer balloon center line connecting direction can be expressed as
φ = π | 12 - t | 12 , ( 6 a . m . ≤ t ≤ 18 p . m . ) - - - ( 11 )
In formula, t means the time, and from the 6:00 AM to the afternoon, 18 is the cycle on a daytime.
H in formula (8) and formula (9) 1and h 2can be expressed as respectively
h 1 = N u 1 λ 1 L - - - ( 12 )
h 2 = N u 1 λ 1 L - - - ( 13 )
In formula, Nu 1for the Nusselt number of forced convection, λ 1for the heat-conduction coefficient of forced convection, Nu 2for the Nusselt number of free convection, λ 2heat-conduction coefficient for free convection.
The heat-conduction coefficient of forced convection can be expressed as
λ 1 = 0.0241 ( T 2 273.15 ) 0.9 - - - ( 14 )
The Nusselt number of forced convection can be expressed as
Nu 1=(0.037Re 4/5-817.0)Pr 1/3 (15)
In formula, Re is Reynolds number, the Prandtl number that Pr is atmosphere.
The Prandtl number of atmosphere can be expressed as
Pr=0.804-3.25×10 -4T 2 (16)
Reynolds number can be expressed as
Re = VL v - - - ( 17 )
In formula, the flowing velocity that V is the stratosphere balloon ambient atmosphere, the characteristic dimension that L is stratosphere balloon, the kinetic viscosity that ν is the stratosphere atmosphere.
The kinetic viscosity of stratosphere atmosphere can be expressed as
v = 1.46 × 10 - 6 T 2 1.5 ρ 1 ( T 2 + 110.4 ) - - - ( 18 )
In formula, ρ 1density for atmosphere.
The temperature of normal atmosphere and density T 2and ρ 1with the stratosphere balloon hoverheight, H can be expressed as respectively
T 2 = 288.15 - 6.5 H H ∈ [ 0,11 km ) T 2 = 216.65 H ∈ [ 11 km , 20 km ) T 2 = 216.65 + ( H - 20 ) H ∈ [ 20 km , 32 km ) - - - ( 19 )
ρ 1 = 1.225 × ( 288.15 - 6.5 H 288.15 ) 4.25588 H ∈ [ 0,11 km ) ρ 1 = 0.364 × e - H - 11 6.34162 H ∈ [ 11 km , 20 km ) ρ 1 = 0.088 × [ 216.65 + ( H - 20 ) 216.65 ] - 35.1632 H ∈ [ 20 km , 32 km ) - - - ( 20 )
While due to stratosphere balloon, reaching energy equilibrium, the free convection energy exchange of its inner helium and capsule cloth is 0, so do not need to provide the expression of the heat-conduction coefficient of the Nusselt number of free convection and free convection again.
In the heat of stratosphere balloon is analyzed, solar constant I 0for 1367W/m 2, atmosphere coefficient of transparency τ is 0.9, ground surface reflectance mean value ρ gbe that 0.3, Si Difen-Boltzmann constant σ is 5.6696 * 10 -8w(m 2k 4), the capsule cloth of stratosphere balloon equates with absorptivity the emissivity of infrared radiation, i.e. ε=α 2.
Respectively formula (2) to formula (9) is carried out to integration, and the substitution correlation parameter, can obtain the total amount of each energy exchange item:
Q 1 ( t ) = 1230.3 π R 2 α 1 , ( 6 a . m . ≤ t ≤ 18 p . m . ) 0 - - - ( 21 )
Q 2 ( t ) = 1139.2 π α 1 R 2 , ( 6 a . m . ≤ t ≤ 18 p . m . ) 0 - - - ( 22 )
Q 3 ( t ) = 366.8 π R 2 α 1 ( R e R e + H ) 2 cos π | 12 - t | 12 , ( 6 a . m . ≤ t ≤ 18 p . m . ) 0 - - - ( 23 )
Figure BDA00003792174900114
Q 5 ( t ) = ∫ A 2 ϵσ T 1 4 dA = 2 AϵAσ T 1 4 = 45.4 π R 2 α 2 × 10 - 8 T 1 4 - - - ( 25 )
Q 6 ( t ) = ∫ A α 2 ϵσ T 1 4 1 - r dA = A α 2 ϵσ T 1 4 1 - r = 22.7 π R 2 α 2 ϵ × 10 - 8 T 1 4 1 - r - - - ( 26 )
Q 7 ( t ) = ∫ A h 1 ( T 1 - T 2 ) dA = A h 1 ( T 1 - T 2 ) = 4 π R 2 h 1 ( T 1 - T 2 ) - - - ( 27 )
Q 8 ( t ) = ∫ A h 2 ( T 1 - T 3 ) dA = A h 2 ( T 1 - T 3 ) = 4 π R 2 h 2 ( T 1 - T 3 ) - - - ( 28 )
Wherein, the infinitesimal described in step 3 refers to the small unit intercepted on stratosphere balloon, and area is infinitely small.
Step 4, according to principle of energy balance, set up respectively stratosphere balloon by day with the energy-balance equation at night.
For daytime, the gross energy of the capsule cloth of stratosphere balloon is
△Q film(t)=Q 1(t)+Q 2(t)+Q 3(t)+Q 4(t)-Q 5(t)+Q 6(t)-Q 7(t)-Q 8(t) (29)
The gross energy of the inside helium of stratosphere balloon is
△Q gas(t)=Q 8(t) (30)
When the energy exchange of stratosphere balloon reaches equilibrium state, △ Q film(t)=0, △ Q gas(t)=0, substitution formula (29) and formula (30) can obtain capsule cloth and inner helium energy-balance equation by day about stratosphere balloon:
Q 1 ( t ) + Q 2 ( t ) + Q 3 ( t ) + Q 4 ( t ) - Q 5 ( t ) + Q 6 ( t ) - Q 7 ( t ) - Q 8 ( t ) = 0 Q 8 ( t ) = 0 - - - ( 31 )
For night, there is no direct solar radiation, sun scattered radiation and earth reflected radiation, i.e. Q 1(t)=0, Q 2(t)=0, Q 3(t)=0, therefore, the capsule cloth of stratosphere balloon and inner helium at the energy-balance equation at night are
Q 4 ( t ) - Q 5 ( t ) + Q 6 ( t ) - Q 7 ( t ) - Q 8 ( t ) = 0 Q 8 ( t ) = 0 - - - ( 32 )
The energy-balance equation of step 5, each energy exchange item substitution step 4 that step 3 is obtained, the simple process of iteration in the recycling method of value solving can try to achieve stratosphere balloon by day with the temperature at night.
Formula (21) to formula (28) is updated to respectively in formula (31) and formula (32), can obtains
( 22.7 π R 2 α 2 ϵ × 10 - 8 1 - r - 45.4 π R 2 α 2 × 10 - 8 ) T 1 4 - 4 π R 2 h 1 T 1 + [ 366.8 π R 2 α 1 ( R e R e + H ) 2 cos π | 12 - t | 12 + 2369.5 π R 2 α 1 + 214.0 π R 2 α 2 ( R e R e + H ) 2 + 4 π R 2 h 1 T 2 ] = 0 T 1 - T 3 = 0
(33)
( 22.7 π R 2 α 2 ϵ × 10 - 8 1 - r - 45.4 π R 2 α 2 × 10 - 8 ) T 1 4 - 4 π R 2 h 1 T 1 + [ 214.0 π R 2 α 2 ( R e R e + H ) 2 + 4 π R 2 h 1 T 2 ] = 0 T 1 - T 3 = 0
(34)
From (33) and formula (34), when the energy exchange of stratosphere balloon reaches equilibrium state, the capsule Bu Wendu of the stratosphere balloon of day and night equates with inner helium temperature.By the simple process of iteration in method of value solving, equation (33) and (34) are solved respectively, can obtain stratosphere balloon by day with capsule Bu Wendu and the inner helium temperature at night.

Claims (1)

1. a new method of predicting the day and night temperature of stratosphere balloon, it is characterized in that: the method concrete steps are as follows:
The assumed condition of the new method of the day and night temperature of step 1, proposition prediction stratosphere balloon;
Assumed condition comprises:
(1) when the stratosphere balloon stand under load, the actual similar water-drop-shaped of shape, for simplified model, be reduced to by stratosphere balloon the circle that radius is R, and the characteristic dimension L of stratosphere balloon means with diameter;
(2) the capsule cloth due to stratosphere balloon is very thin, and in micron dimension, therefore, the temperature of establishing capsule cloth thickness direction is identical, does not consider the heat transmission between capsule cloth internal layer and skin;
(3) do not consider the loss of the inner helium of stratosphere balloon, the weight of establishing stratosphere balloon is constant;
(4) due to earth radius R efor 6371.23km, and the sun apart from the earth the distance be 1.5 * 10 8therefore km, is directional light while establishing sunshine arrival earth overhead;
(5) establishing the earth is the thermal equilibrium body;
(6) do not consider the variation of intensity of solar radiation with solar radiation angle and atmospheric extinction coefficient;
(7) do not consider the impact of carbon dioxide and ozone infrared radiation;
Step 2, according to the environment for use of stratosphere balloon and the ABC of Engineering Thermodynamics, determine the energy exchange mechanism of stratosphere balloon and each energy exchange item;
Stratosphere balloon carries out energy exchange by radiation and convection current, comprising: the direct solar radiation item Q that the stratosphere balloon episphere absorbs 1(t) sun scattered radiation item Q that, stratosphere balloon absorbs 2(t) earth reflected radiation item Q that, the stratosphere balloon lower semisphere absorbs 3(t) earth infrared radiation item Q that, the stratosphere balloon lower semisphere absorbs 4(t), the external infrared radiation item Q of stratosphere balloon 5(t) the internal reflection infrared radiation item Q that, stratosphere balloon absorbs 6(t), the convection heat transfer item Q between stratosphere balloon and outside atmosphere 7and the convection heat transfer item Q between stratosphere balloon and internal gas (t) 8(t);
Step 3, each energy exchange item proposed for step 2 are set up respectively its thermal model, obtain the concrete concrete Parameter Expression form of each energy exchange item;
The cartographic represenation of area of the infinitesimal of stratosphere balloon is
In formula, dA is the unit dimension on stratosphere balloon, the latitude that θ is spherical coordinate system or solar radiation angle or sunshine and imfinitesimal method between angle,
Figure FDA00003792174800012
longitude for spherical coordinate system;
The direct solar radiation that infinitesimal on the stratosphere balloon episphere absorbs is
dQ 1(t)=α 1τI 0cosθdA (2)
In formula, α 1for the capsule cloth of the stratosphere balloon absorptivity to direct solar radiation, τ is the atmosphere coefficient of transparency, I 0for solar constant;
The sun scattered radiation that infinitesimal on stratosphere balloon absorbs is
dQ 2(t)=α 1I 1dA (3)
In formula, I 1for sun diffuse radiation;
The earth reflected radiation that infinitesimal on the stratosphere balloon lower semisphere absorbs is
d Q 3 ( t ) = α 1 ρ g τ I 0 ( R e R e + H ) 2 cos θ cos φdA - - - ( 4 )
In formula, ρ gfor ground surface reflectance mean value, R efor earth radius, the hoverheight of H stratosphere balloon, φ is the angle between sunshine direction and ground ball center peace fluid layer balloon center line connecting direction;
The earth infrared radiation that infinitesimal on the stratosphere balloon lower semisphere absorbs is
d Q 4 ( t ) = α 2 τ I 0 R e 2 1 - ρ g 4 ( R e + H ) 2 cos θdA - - - ( 5 )
In formula, α 2for the capsule cloth of the stratosphere balloon absorptivity to infrared radiation;
The external infrared radiation of infinitesimal on stratosphere balloon is
d Q 5 ( t ) = 2 ϵσ T 1 4 dA - - - ( 6 )
In formula, the emissivity of the capsule cloth that ε is stratosphere balloon to infrared radiation, σ is Si Difen-Boltzmann constant, T 1capsule Bu Wendu for stratosphere balloon;
The internal reflection infrared radiation that infinitesimal on stratosphere balloon absorbs is
d Q 6 ( t ) = α 2 ϵσ T 1 4 1 - r dA - - - ( 7 )
In formula, the reflectivity that r is stratosphere balloon;
Infinitesimal and the convection heat transfer between outside atmosphere on stratosphere balloon are
dQ 7(t)=h 1(T 1-T 2)dA (8)
In formula, h 1for the convection transfer rate between stratosphere balloon and outside atmosphere, T 2temperature for outside atmosphere;
Convection heat transfer between infinitesimal on stratosphere balloon and inner helium is
dQ 8(t)=h 2(T 1-T 3)dA (9)
In formula, h 2for the convection transfer rate between stratosphere balloon and inner helium, T 3temperature for inner helium;
Sun diffuse radiation I in formula (3) 1be expressed as
I 1 = 1 2 I 0 cos θ cos 2 θ 2 - - - ( 10 )
Included angle between sunshine direction in formula (4) and ground ball center peace fluid layer balloon center line connecting direction is expressed as
φ = π | 12 - t | 12 , ( 6 a . m . ≤ t ≤ 18 p . m . ) - - - ( 11 )
In formula, t means the time, and from the 6:00 AM to the afternoon, 18 is the cycle on a daytime;
H in formula (8) and formula (9) 1and h 2be expressed as respectively
h 1 = N u 1 λ 1 L - - - ( 12 )
h 2 = N u 1 λ 1 L - - - ( 13 )
In formula, Nu 1for the Nusselt number of forced convection, λ 1for the heat-conduction coefficient of forced convection, Nu 2for the Nusselt number of free convection, λ 2heat-conduction coefficient for free convection;
In the heat of stratosphere balloon is analyzed, solar constant I 0for 1367W/m 2, atmosphere coefficient of transparency τ is 0.9, ground surface reflectance mean value ρ gbe that 0.3, Si Difen-Boltzmann constant σ is 5.6696 * 10 -8w(m 2k 4), the capsule cloth of stratosphere balloon equates with absorptivity the emissivity of infrared radiation, i.e. ε=α 2;
Respectively formula (2) to formula (9) is carried out to integration, and the substitution correlation parameter, obtains the total amount of each energy exchange item:
Q 1 ( t ) = 1230.3 π R 2 α 1 , ( 6 a . m . ≤ t ≤ 18 p . m . ) 0 - - - ( 13 )
Q 2 ( t ) = 1139.2 π α 1 R 2 , ( 6 a . m . ≤ t ≤ 18 p . m . ) 0 - - - ( 14 )
Q 3 ( t ) = 366.8 π R 2 α 1 ( R e R e + H ) 2 cos π | 12 - t | 12 , ( 6 a . m . ≤ t ≤ 18 p . m . ) 0 - - - ( 15 )
Figure FDA00003792174800036
Q 5 ( t ) = ∫ A 2 ϵσ T 1 4 dA = 2 AϵAσ T 1 4 = 45.4 π R 2 α 2 × 10 - 8 T 1 4 - - - ( 17 )
Q 6 ( t ) = ∫ A α 2 ϵσ T 1 4 1 - r dA = A α 2 ϵσ T 1 4 1 - r = 22.7 π R 2 α 2 ϵ × 10 - 8 T 1 4 1 - r - - - ( 19 )
Q 7 ( t ) = ∫ A h 1 ( T 1 - T 2 ) dA = A h 1 ( T 1 - T 2 ) = 4 π R 2 h 1 ( T 1 - T 2 ) - - - ( 20 )
Q 8 ( t ) = ∫ A h 2 ( T 1 - T 3 ) dA = A h 2 ( T 1 - T 3 ) = 4 π R 2 h 2 ( T 1 - T 3 ) - - - ( 21 )
Wherein, the infinitesimal described in step 3 refers to the small unit intercepted on stratosphere balloon, and area is infinitely small;
Step 4, according to principle of energy balance, set up respectively stratosphere balloon by day with the energy-balance equation at night;
For daytime, the gross energy of the capsule cloth of stratosphere balloon is
△Q film(t)=Q 1(t)+Q 2(t)+Q 3(t)+Q 4(t)-Q 5(t)+Q 6(t)-Q 7(t)-Q 8(t) (22)
The gross energy of the inside helium of stratosphere balloon is
△Q gas(t)=Q 8(t) (23)
When the energy exchange of stratosphere balloon reaches equilibrium state, △ Q film(t)=0, △ Q gas(t)=0, substitution formula (22) and formula (23) obtain capsule cloth and inner helium energy-balance equation by day about stratosphere balloon:
Q 1 ( t ) + Q 2 ( t ) + Q 3 ( t ) + Q 4 ( t ) - Q 5 ( t ) + Q 6 ( t ) - Q 7 ( t ) - Q 8 ( t ) = 0 Q 8 ( t ) = 0 - - - ( 24 )
For night, there is no direct solar radiation, sun scattered radiation and earth reflected radiation, i.e. Q 1(t)=0, Q 2(t)=0, Q 3(t)=0, therefore, the capsule cloth of stratosphere balloon and inner helium at the energy-balance equation at night are
Q 4 ( t ) - Q 5 ( t ) + Q 6 ( t ) - Q 7 ( t ) - Q 8 ( t ) = 0 Q 8 ( t ) = 0 - - - ( 25 )
The energy-balance equation of step 5, each energy exchange item substitution step 4 that step 3 is obtained, the simple process of iteration in the recycling method of value solving try to achieve stratosphere balloon by day with the temperature at night;
Formula (13) to formula (21) is updated to respectively in formula (24) and formula (25),
( 22.7 π R 2 α 2 ϵ × 10 - 8 1 - r - 45.4 π R 2 α 2 × 10 - 8 ) T 1 4 - 4 π R 2 h 1 T 1 + [ 366.8 π R 2 α 1 ( R e R e + H ) 2 cos π | 12 - t | 12 + 2369.5 π R 2 α 1 + 214.0 π R 2 α 2 ( R e R e + H ) 2 + 4 π R 2 h 1 T 2 ] = 0 T 1 - T 3 = 0
(26)
( 22.7 π R 2 α 2 ϵ × 10 - 8 1 - r - 45.4 π R 2 α 2 × 10 - 8 ) T 1 4 - 4 π R 2 h 1 T 1 + [ 214.0 π R 2 α 2 ( R e R e + H ) 2 + 4 π R 2 h 1 T 2 ] = 0 T 1 - T 3 = 0
(27)
From (26) and formula (27), when the energy exchange of stratosphere balloon reaches equilibrium state, the capsule Bu Wendu of the stratosphere balloon of day and night equates with inner helium temperature; By the simple process of iteration in method of value solving, equation (26) and (27) are solved respectively, obtain stratosphere balloon by day with capsule Bu Wendu and the inner helium temperature at night.
CN201310406617.9A 2013-09-09 2013-09-09 A kind of new method predicting the day and night temperature of stratosphere balloon Expired - Fee Related CN103471723B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310406617.9A CN103471723B (en) 2013-09-09 2013-09-09 A kind of new method predicting the day and night temperature of stratosphere balloon

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310406617.9A CN103471723B (en) 2013-09-09 2013-09-09 A kind of new method predicting the day and night temperature of stratosphere balloon

Publications (2)

Publication Number Publication Date
CN103471723A true CN103471723A (en) 2013-12-25
CN103471723B CN103471723B (en) 2015-12-09

Family

ID=49796666

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310406617.9A Expired - Fee Related CN103471723B (en) 2013-09-09 2013-09-09 A kind of new method predicting the day and night temperature of stratosphere balloon

Country Status (1)

Country Link
CN (1) CN103471723B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105303053A (en) * 2015-11-12 2016-02-03 中国科学院光电研究院 Level flight process average temperature calculating method of high altitude balloon with solar cell
CN105404774A (en) * 2015-11-12 2016-03-16 中国科学院光电研究院 Method for computing distributed temperature in level flight of high altitude balloon with solar cell

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2429371A1 (en) * 2000-10-05 2002-04-11 W. Scott Andrus Temperature measuring balloon
CN201254288Y (en) * 2008-06-20 2009-06-10 金玉 Electric energy hot air balloon
CN101720422A (en) * 2007-06-20 2010-06-02 维萨拉公司 Method for compensating for temperature measurement error in a sond.

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2429371A1 (en) * 2000-10-05 2002-04-11 W. Scott Andrus Temperature measuring balloon
CN101720422A (en) * 2007-06-20 2010-06-02 维萨拉公司 Method for compensating for temperature measurement error in a sond.
CN201254288Y (en) * 2008-06-20 2009-06-10 金玉 Electric energy hot air balloon

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
QIUMIN DAI ET AL.: "Performance simulation of high altitude scientific balloons", 《ADVANCES IN SPACE RESEARCH》 *
秦玉梅 等: "高空气球加强带热合温度数值仿真", 《计算机仿真》 *
雷岩鹏 等: "火星超压气球温度分布研究", 《航空学报》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105303053A (en) * 2015-11-12 2016-02-03 中国科学院光电研究院 Level flight process average temperature calculating method of high altitude balloon with solar cell
CN105404774A (en) * 2015-11-12 2016-03-16 中国科学院光电研究院 Method for computing distributed temperature in level flight of high altitude balloon with solar cell
CN105303053B (en) * 2015-11-12 2018-10-12 中国科学院光电研究院 Aerostat with solar cell is flat to fly over journey mean temperature computational methods
CN105404774B (en) * 2015-11-12 2018-10-12 中国科学院光电研究院 Aerostat with solar cell is flat to fly over journey districution temperature computational methods

Also Published As

Publication number Publication date
CN103471723B (en) 2015-12-09

Similar Documents

Publication Publication Date Title
Zhang et al. Thermal performance analysis of a high-altitude solar-powered hybrid airship
Holstein‐Rathlou et al. Winds at the Phoenix landing site
Newman et al. Modeling the Martian dust cycle, 1. Representations of dust transport processes
Yao et al. A heat transient model for the thermal behavior prediction of stratospheric airships
Heavens et al. Structure and dynamics of the Martian lower and middle atmosphere as observed by the Mars Climate Sounder: 2. Implications of the thermal structure and aerosol distributions for the mean meridional circulation
Guerlet et al. Evolution of the equatorial oscillation in Saturn's stratosphere between 2005 and 2010 from Cassini/CIRS limb data analysis
Du et al. Thermal performance analysis and comparison of stratospheric airships with rotatable and fixed photovoltaic array
Lv et al. Effects of solar array on the thermal performance of stratospheric airship
Liu et al. Experimental and numerical study on the aero-heating characteristics of blunted waverider
Liu et al. Thermal performance of stratospheric airship with photovoltaic array
Hecht et al. Maui Mesosphere and Lower Thermosphere (Maui MALT) observations of the evolution of Kelvin‐Helmholtz billows formed near 86 km altitude
Hollingsworth et al. Extratropical cyclones, frontal waves, and Mars dust: Modeling and considerations
Li et al. Modeling and analysis of floating performances of stratospheric semi-rigid airships
Wang et al. A comprehensive numerical model examining the thermal performance of airships
Ai et al. Heat generation effects from shear friction along Xianshui river strike-slip fault in western Sichuan, China
Józsa et al. Solving Problems in Thermal Engineering
Xiong et al. Simplified analytical model for predicting the temperature of balloon on high-altitude
Zhu et al. Thermal performance of high-altitude solar powered scientific balloon
Zhang et al. Influences of initial launch conditions on flight performance of high altitude balloon ascending process
Soria-Salinas et al. Wind retrieval from temperature measurements from the rover environmental monitoring Station/Mars science laboratory
Li et al. Experimental investigation of transient thermal behavior of an airship under different solar radiation and airflow conditions
Yates et al. Magnetosphere‐ionosphere‐thermosphere coupling at Jupiter using a three‐dimensional atmospheric general circulation model
Samanta et al. Computational modeling and experiments of natural convection for a titan Montgolfiere
Alam et al. A multi-node model for transient heat transfer analysis of stratospheric airships
CN103471723B (en) A kind of new method predicting the day and night temperature of stratosphere balloon

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20151209

Termination date: 20200909