Balloon Figure 1 Image from a Balloon from Actarus site

Figure 1 – A picture from the balloon Actarus project

A convenient way to gather atmospheric data at high altitudes is to use weather balloons (current record altitude for unmanned gas filled balloon seems to be 53 km of BU60-1 Japanese balloon). As this kind of devices is used for weather forecast and research the specific literature on the topic is widely available [1] and consolidated. Air data systems are a common payload across that kind of probe, from here our interest on the topic. Here we introduce the basic hardware, provide a basic calculation script for balloon sizing and describe the ground filling operations. The topics involved are different and heterogeneous so we will split topics along a mini-series of articles on balloons.

We will focus on the most common DIY setup, a sealed balloon filled with helium and with a rope suspended payload. Helium is lighter than air so it makes the balloon float.

System Layout

Any component you chose to use should be reliable, safe and disposable. Redundancy can be the winning choice. The payload can be put within a box or exposed to the atmosphere, in both cases consider the low ambient temperatures; temperatures can affect electronics, and batteries, nominal operation. In the tropopause the temperature can fall below -50°C.
Refer to Figure 2, the probe is composed of five main parts. The balloon, the rigging, the recovery parachute, the payload and the balloon filling gas.

The typical setup is shown in Figure 3. We will focus on latex sealed balloons, this kind of balloons will expand as the probe ascends and will burst when their maximum diameter is reached. Availability is good and they are affordable. Medium high altitudes, let’s say around 30 Km, can be reached without too much hassle. The balloon is secured to the hanging payload by means of strings, any string or accessory used for this purpose is part of the rigging. That is not an issue for small payloads but it is of relevance with heavy payloads; balloon wall is thin, it’s not able to support heavy loads, a concentrated load will pierce the envelope. The parachute will slow down the payload descent when the balloon explodes, during the ascent, the parachute is already exposed and hangs from the rigging.

Balloon Figure 2

Figure 2 – General Layout of a Sounding Balloon

The payload is composed of all the hardware that should fly with the probe, we can divide the equipment into two big families. Common equipment and mission specific hardware. The first family includes all the hardware that is absolutely necessary to fly, typical hardware listed in Table 1.

  • Power source, battery pack. Can be one single power source for every component of the probe or the power can be distributed among different devices. It is to be considered that power source is mission critical. For most important probes it’s possible to use dual failsafe power sources.
  • Position tracker. Any device needed to locate the payload after touchdown. With low weight probes, a GPS-based locator can be used. See for example. The probe will be located after the touchdown with a cell phone. Most complex probes can have a full duplex telemetry radio link active during the flight.

Table 1 – Common equipment

Besides that, you need mission specific hardware. You find a typical list of components in Table 2.

  • Camera. You can have also multiple cameras. For example, one looking downwards and another looking to the horizon. A useful feature is to have the possibility to shoot pictures at determined times. Pictures are usually stored on a memory card or solid state device, the storage limit has an impact on the possibility to have a full photo coverage of the flight.
  • Inertial Movement Unit. It provides acceleration, rotation and position information on the three axes. They are compact and can be equipped with three-axis magnetometers. (See for example)
  • Air Data Computer. Equipped with sensors for external temperature, pressure, and humidity.
  • Vertical speed probe. Apart inertial data it is possible to have a dedicated probe to measure vertical airspeed.

Table 2 – Mission-specific equipment

Balloon Figure 3 120623125936-cnn-weather-balloon-launch-story-top

Figure 3 – A Latex Weather Balloon

Balloon Math Modeling

To produce a really simple, and conservative, model we will work under strong simplifications; generated model can be used during the hardware selection phase but it is not suitable for accurate performance prediction nor for flight simulation.

Ascending Balloon [1] [2]

Refer to figure 5, let’s model the whole probe as a mass point. Our probe floats thanks to the buoyancy. The upthrust UT available using a sealed balloon is proportional to the weight of displaced air UT=\rho_{air}BVg where BV is the balloon volume and g is the gravitational acceleration.
The balloon, of diameter D will expand as the air density decreases, as a direct consequence the forces exerted on the balloon material will increase. In the specific case [3] the stress on the material section will be equal to \sigma_b=\frac{(pi-p)D}{4t}. Apart the punctual value the stress is to be noted that if the stress is greater than zero then we have that internal pressure p_i is greater than air pressure p. The pressure inside the balloon is not equal to the external pressure. To avoid the mechanical treatment of the mechanical expansion we will work under the hypothesis of “rigid balloon” or the so called not stressed balloon, the internal pressure is exactly the air pressure. During the ascension, the balloon will expand because of the temperature change to the air and to the helium. During the ascent the air temperature will change, the helium temperature will change too. The temperature distribution across the balloon volume will not be uniform, as a simplification in our model the temperature will be a determined parameter. The temperature across the atmosphere will be assumed to be distributed as per the U.S. 1976 Standard Atmosphere [4]. Solar radiation contributes to balloon energy balance [2] so will have an impact on expansion, we will neglect also that contribution, for the peace of the mind think that you are estimating the performance of a night launched balloon; bibliography reports [2] lower ascent ratios during the night.

Balloon Figure 4 Parachute

Figure 4 – Rigged Recovery Parachute

The weight force acting on the probe is W=m_t g, where m_t is the total probe mass. Another force acting on the probe is the drag exerted by the air the probe is moving in. The direction of drag is dependent on the lateral movement of the balloon. Laterals dynamic of the balloon is dominated by horizontal wind component. For example purposes we will consider a small angle of attack \alpha. Drag force DF expression is: [2]

(1)   \begin{equation*} DF={1 \over 2} C_{dr}\rho_{air}A_r W^2_r \end{equation*}

Where W^2_r is the still air rise rate in m/s, A_r is the cross area of the balloon and C_{dr} is the coefficient of drag for the balloon. The coefficient of drag is dependent on \alpha and Re. This drag coefficient should consider the contribute of the payload and rigging, but this former contributes are often neglected, most accurate way to determine coefficient value is by flight test [2] or wind tunnel testing.

A less accurate approximation is to consider the balloon as a perfect sphere, in such a way we can use the coefficient of drag for a sphere. We will operate with that assumption. The coefficient of drag of a sphere depends on the Reynolds number [2], so on the density; for this reason, the C_{dr} should be considered a function of Re and practically of altitude. In our model we will set the drag coefficient at 0.47; that will be a value too high for ground conditions and far more appropriate for low Re high altitude conditions. Setting a high will overestimate the drag force near the ground, that will give to us an optimistic or conservative estimation of balloon performances as the ascent rate.

Figure 5 – Probe Model

Given the mass of the payload and rigging m_p, the mass of the balloon m_b and the mass of helium contained within the balloon m_h we get that the net buoyancy force BF is:

(2)   \begin{equation*} BF=BV\rho g-(m_p+m_b+m_h)g \end{equation*}

Combing Equation 1 and 2 we get the probe rise rate in still air:

(3)   \begin{equation*} W_r=\sqrt{{2BF} \over {C_dA_r\rho}} \end{equation*}

 

Descending Balloon [1] [2]

At the balloon bursts, the probe will initiate to fall back to Earth. To avoid mushing of the payload the recovery parachute should be correctly sized. Reusing ascending equations and replacing the parachute coefficient of drag C_{dd} instead of the balloon C_d we get the expression of vertical descending speed.

(4)   \begin{equation*} W_d=\sqrt{{2m_t} \over {C_{dd}A_d\rho}} \end{equation*}

Where A_d is the area of the parachute. Equation 4 neglects the contribute of the payload drag, likely we will overestimate the descending speed.

Check your local regulations

Prior to any launch please check your local regulations about unmanned gas balloons launch. To date (15/04/2016) in Italy in necessary to notify the launch to ENAC, with this modules. The Italian Civil Aviation Authority names the ICAO [5] Annex 2 Appendix 4 as the source for determining what data should be communicated prior to launch to ENAC.
ICAO define different classes of unmanned balloons, classification is based on payload weight. For example, a balloon with a single package payload(anything different from the balloon) that weights less than 3kg is considered as a light balloon. To be considered light the payload should be suspended with a rope that requires an impact force less than 230N to be separated from the balloon. Another requirement is to have an area density of less than 13g per square centimeter. To calculate the area density you should divide the mass in grams of the payload package by the area if square centimeters of the package smallest surface. A light balloon has less operating limitations respect to heavy balloons.

Pick your equipment

Balloon Figure 6 Balloon classification

Figure 6 – Classification of unmanned free balloon from ICAO

The design is an iterative process, sometimes you should redo the same procedure over and over up to the solution convergence. At first, you need to decide you mission objective, the more constraints you have the faster you will found a solution.

For example purposes I selected the following requirements.
TBA, Target burst altitude, 32000m
m_{\cap{p}}, Payload weight, 0.6 kg
m_{\cap{b}}, Ballon weight, 0.6 kg
BBD, Balloon burst diameter by producer information, m
f_t, Fly time, 180 minutes
c_{dd}, Parachute coefficient of drag, first guess value 1.75
W_{\cap{r}}, Vertical ascending speed, 3 m/s
W_{\cap{d}}, Vertical descending speed at landing, 10 m/s

Note that in the list it is included the burst diameter, that data should be supplied by the balloon producer.
As a result of the preliminary design, we will determine the required balloon size, the mass of helium to load into the balloon at ground and the buoyancy at ground. You can see the desired design output parameters in Table 3.

FL, Net buoyancy force or free lift, N
m_h, Mass of helium, kg
m_p, Mass of balloon, kg
W_r, Vertical speed at ground at m_{\cap{p}}, m/s
r_r, Rise time m_{\cap{p}}, hour
f_t, Fall time m_{\cap{p}}, hour
max(m_{p}), Maximum allowable payload weight to reach TBA, kg
BG, Buoyancy at ground, N

Table 3 – Design output parameters

We have a model of our balloon and we have defined design output parameters, within the next post we will use a Scilab/Matlab script for sizing.

Example sizing

Here we walk through a preliminary design example. Recalling the model of the ascending balloon and the rigid balloon hypothesis we observe that net buoyancy FL during the ascent is constant. The calculated data will be used to launch, so inflate, properly the balloon. Measuring the balloon volume it is not an easy way to go, usually, the easier way is to measure free lift force at ground. You can use this [download id=”2565″] to calculate lift at ground, rise time, fall time, average speed during rising, and the required mass of helium. Input data will be burst altitude [m], payload mass [kg], balloon mass, balloon burst diameter [m], fly time in minutes, parachute coefficient of drag, vertical descending desired speed at ground[m/s]. To launch a preliminary design, with our example data, we invoke the SizeProbe function in the following way:

[FL, RiseTime, FallTime, wr,mHe]=SizeProbe(32000,0.6, 0.6,5.791,180,10,1.75)

The results will be:

mHe = 0.1983532
wr = 2.1052632
FallTime = 3200
RiseTime = 7600
FL = 0.3628609

Burst diameter should be read from producer data or found by experience, for example, see this table. You see that our preliminary design gives to us an indication that our balloon sizing can provide the required performances.

Another alternative and quick way to proceed is to use online design tools that are prepared by the balloon producers, for example.

I report here some other helpful online design calculator:

https://www.highaltitudescience.com/pages/balloon-performance-calculator

https://www.stratoflights.com

https://www.scientificsales.com/Meteorological-Weather-Sounding-Balloon-s/25.htm

Typically those tools are designed to be on the safer side, they will be pessimistic about performances; results are good for design but do not provide an accurate balloon behavior description. We have other aspects to consider before the launch, as the trajectory prediction. Balloon trajectory is highly dependent on the wind. Once again we’ve to consider that the available data on wind is far from to be complete or exact; often we desire to know a trajectory in the next future and so the simulations will be based on a wind forecast. Beside that we can find some trajectory predictor, for example, this in this link from CUSF. From the predictions, you will notice that balloons will cover wide distances, it is necessary to have one or more well-prepared recovery teams that can easily travel across the country.

We’ve exposed some aspects of balloon operation. Key factors are preparation, training, and performance estimation.

Happy landings!

References

[1] A. Gallice , F. G. Wienhold , C. R. Hoyle, F. Immler, and T. Peter(2011), Modeling the ascent of sounding balloons:derivation of the vertical air motion, Atmospheric Measurement Techniques Retrieved Online 01/07/2016
[2] Wang et al (2009), Vertical Air Motion from T-REX Radiosonde and Dropsonde Data, AMS  Retrieved Online 01/07/2016
[3] CAS/courses.d/Structures.d/IAST.Lect03.d/ Thin Walled Pressure Vessels, Retrieved Online 15/04/2016
and this doc
[4] U.S. 1976 Standard Atmosphere , Retrieved Online 15/04/2016
[5] ICAO (2006), Annex 2 Appendix 4, Retrieved Online 15/04/2016

0 Comments

Leave a Reply

Avatar placeholder