Next Article in Journal
Systems Integration Framework for Hybrid-Electric Commuter and Regional Aircraft
Next Article in Special Issue
A Lunar-Orbiting Satellite Constellation for Wireless Energy Supply
Previous Article in Journal
Research on Enhanced Orbit Prediction Techniques Utilizing Multiple Sets of Two-Line Element
Previous Article in Special Issue
Short-Arc Horizon-Based Optical Navigation by Total Least-Squares Estimation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Conceptual Design of Hybrid Aerial Vehicle for Venus Exploration

1
Mechanical and Aerospace Engineering Department, New Mexico State University, Las Cruces, NM 88003, USA
2
Astronomy Department, New Mexico State University, Las Cruces, NM 88003, USA
*
Author to whom correspondence should be addressed.
Aerospace 2023, 10(6), 534; https://doi.org/10.3390/aerospace10060534
Submission received: 17 April 2023 / Revised: 27 May 2023 / Accepted: 1 June 2023 / Published: 3 June 2023
(This article belongs to the Special Issue Advanced Spacecraft/Satellite Technologies)

Abstract

:
The conceptual design of a hybrid aerial vehicle for the exploration of the upper Venus atmosphere is presented. The vehicle will float like a balloon and harvest solar energy which is stored in batteries. The neutral buoyancy reduces the energy consumption and makes the vehicle robust and durable. Energy stored in the batteries can be used for powered flight with good horizontal and vertical mobility to explore aspects of the atmosphere. The vehicle is intended to operate near 55.3 km altitude and to explore the cloud layer of the planet. The vehicle takes its inspiration from the Stingray inflatable wing by Prospective Concepts. Based on a trade study, the wing span was set to 25 m. Equations are developed for the altitude, gas and skin temperature, and skin stress during neutrally buoyant flight. To keep the equations in a simplified analytical form, the complex compartmentalized gas pockets of the vehicle are lumped into a single gas sphere. The equations take into account the volumetric expansion of the structure and the requirement that the differential pressure needs to be large enough to allow for brief periods of powered flight without significant structural deformation. An aerodynamic analysis provides the lift and drag coefficient curves and indicates that the vehicle is pitch-stable. A powered flight analysis shows that an airspeed of 30 m/s can be maintained for 31 min at 55 km and 69 min at 69 km altitude.

1. Introduction

Despite the fact that Venus is Earth’s closest planetary neighbour, the present knowledge of the Venus atmosphere composition and physics is limited. The Venus atmosphere consists mainly of carbon dioxide, CO 2 . Dense sulphuric acid clouds obscure the surface and reflect most of the incoming sunlight. The pressure and temperature at 55 km altitude are similar to sea level conditions on Earth [1] and may support life. An improved understanding of the Venus atmosphere could potentially provide new insight into the Earth’s greenhouse effect and long-term climate change on Earth-like planets.
In-situ measurements for the stable layer above the middle clouds at about 50–53.5 km altitude were obtained by the Vega probes in 1985. Horizontal winds in the Venus atmosphere can be as high as 100 m/s. The rotation rate of Venus is the lowest of all planets in the solar system (243 Earth days). As a result super-rotation occurs, which refers to a situation where the wind speed exceeds the equatorial tangential velocity of the planet. An unstable region with strong vertical shear develops at night at an altitude of roughly 63 km. At low latitudes, the unstable region disappears in the morning. At higher latitudes, the unstable region may persist for most of the day. As a result of the extreme surface conditions, the lifetime of Venus landers has been very short [2] and new studies have proposed minimum requirements for future landers [3]. Navigation is challenging as Venus has no magnetic field.
New concepts have been proposed for in-situ measurements of the Venus atmosphere [2,4,5]. Buoyant vehicles such as balloons are robust and resilient, making them well suited for Venus missions. For example, a variable altitude balloon (or aerobot) was proposed by Hall et al. [6] for exploring the atmosphere between 52 and 62 km altitude. To alter the buoyancy of the vehicle (and hence the altitude), helium gas is exchanged between a high-pressure balloon and a zero-pressure balloon. the balloon trajectory was computed with a model by Carlson and Horn [7]. Agrawal et al. [8] proposed a variable altitude balloon-type aerobot that releases ballast and lifting gas to adjust its altitude. The objective is to explore the habitability of the planet’s cloud layer at 48–60 km altitude. Other balloon-type concepts were considered by the 2020 Venus Flagship Mission Study [5]. The general intent of the missions is to investigate the composition, physical processes, and meteorology of the atmosphere. Although balloons offer many advantages (such as simplicity and zero-power requirement for remaining afloat), more manoeuvrable vehicle designs are attractive for an active exploration of the atmosphere. A larger number of researchers investigated solar-powered aircraft for flight in the Venus atmosphere [9,10,11,12,13,14]. For example, Herkenhoff et al. [13] estimated the available solar power and Acosta et al. [14] computed the required wing and thrust loading. Dynamic soaring refers to a flight trajectory that extracts energy from velocity shear and this concept was explored by David [15] and Bullock et al. [16]. The proposed fixed-wing aircraft are often intended for operation above the clouds where the atmospheric conditions are more Earth-like. The large uncertainty with respect to the atmospheric conditions at such altitudes does however raise questions with respect to the stability and control of aircraft. Concepts with aerodynamic lift (that are manoeuvrable such as aircraft) and sufficient buoyancy to maintain altitude when floating, appear to be a good compromise for Venus exploration. They are potentially more resistant to gusts. Vertical or sheer wind gusts could have an effect on hybrid airships; however, they are believed to be less likely on Venus. The main effect would be on the control system, trying to maintain heading. Young et al. [17] proposed a twin-hull Venus airship concept for flight at low altitudes. The Venus atmospheric manoeuvrable platform (VAMP) was designed for operation between 55 and 70 km [18].
This paper addresses the conceptual design of a powered aerodynamic buoyant vehicle with horizontal and vertical mobility. The vehicle geometry was modelled after the Stingray inflatable wing by Prospective Concepts [19,20]. Efficient aerodynamics (low induced drag) make the Stingray efficient (low power consumption) and provide good horizontal and vertical mobility. Of relevance, a stingray-shaped buoyant inflatable flying wing with flapping wings was developed by Festo [21]. To summarize, the buoyant aircraft concept was chosen for the present study, because it offers many advantages over balloons and aircraft such as robustness and resilience (vehicle will not crash if propulsion system fails or vehicle looses flight control), manoeuvrability (compared to balloons), energy efficiency (low power consumption in forward flight), and compactness when deflated (to fit into aeroshell).
The paper is organized as follows: The atmospheric environment at the proposed cruise altitude of 55.3 km is characterized and the envisioned operational strategy is discussed. The vehicle geometry and component mass are discussed and the battery and solar cells are sized. Differential equations for the vehicle altitude and gas temperature are proposed and utilized to investigate the vehicle flight path. The skin thickness is obtained from a stress analysis. The lift, drag, and moment coefficient are computed over a relevant angle of attack range.
Based on the battery charge and airspeed, the maximum powered flight time for two different altitudes are computed. Finally, a brief paragraph outlining steps that can be taken for a more detailed analysis, as well as a summary and conclusions are offered.

2. Mission Environment

The Venus atmosphere density, pressure and temperature for altitudes between z = 50 km and 100 km are plotted in Figure 1. Near 55 km, an altitude of interest to atmospheric scientists [5], the pressure and temperature are similar to Earth at sea level. Therefore, the average cruise altitude for the present conceptual vehicle design was set to 55.3 km. This information is based on the Venus International Reference Atmosphere (VIRA) and other data, as described by Zasova et al. [1]. The data in this reference are provided for different solar latitudes. For the present analysis, the data were averaged over a day–night cycle. This approach is justified for a conceptual design. The solar irradiation, I, is plotted in Figure 1b. The experimental data (symbols) were obtained from Tomasko et al. [22]. The solid line constitutes a fourth-order polynomial fit. Table 1 summarizes the atmospheric conditions and solar irradiation at 55.3 km.
The solar flux in Venus orbit can be as high as 2622 W/m 2 [23]. The dense Venus clouds between 40 km and 70 km altitude reflect more than 75% of the solar radiation back into space [24] and make Venus the planet with the highest albedo, r e , in the solar system. Here, albedo refers to the percentage of reflected solar radiation. For comparison, the albedo is 68.9% for Venus [25] and 43.4% for Earth [26]. The one-third of radiation that is not reflected is absorbed by Venus and remitted as infrared (IR) radiation (153 W/m 2 according to the Spacecraft Thermal Control Book [27]). The IR radiation is absorbed by the clouds and scattered in all directions, including back down towards the surface. This causes a strong greenhouse effect. The clouds of Venus are composed of sulphuric acid and, as a result, all exposed surfaces of atmospheric vehicles need to be corrosion-resistant.

3. Envisioned Operational Strategy

The vehicle will float with the predominant winds and circumnavigate the planet near the equator. The vehicle will pass over the sunlit side and then enter the dark side (see Figure 2). Compared to aircraft, the vehicle will require no power for remaining airborne which will make it robust and tolerant to upset conditions. For flight at an altitude of 55.3 km and a wind speed of ≈60 m/s [28], it will take the vehicle ≈178 h to circumnavigate the planet. This time period must be considered a rough estimate as the uncertainty of the wind velocity measurements is high and the wind velocity will vary depending on environmental factors. To the best of the author’s knowledge, presently no reliable data on wind gusts in the Venus atmosphere are available.
The vehicle will receive direct solar irradiation from above, reflected (albedo) radiation from below, and IR radiation from the surface of the planet, all of which will increase the vehicle temperature. Natural and forced convection (in forward flight) as well as IR radiation from the vehicle into the environment will affect the vehicle temperature. Regarding the convective heat transfer, the skin temperature will increase when the surrounding air is warmer than the vehicle and decrease otherwise. When the vehicle emerges from the dark side, thermal expansion of the buoyant gas, as a consequence of solar heating, will make the vehicle climb. This change in volume will deform the geometry and affect the aerodynamic characteristics of the vehicle during powered flight. To minimize the deformation, the vehicle volume can be split up into buoyant gas pockets such as seen for the Prospective Concepts [20] and Festo [21] vehicles. The detailed analysis can be part of a later more detailed design. Figure 3 illustrates the resulting mission profile.
The vehicle geometry will resemble that of a flying wing. Solar cells will be mounted to the upper wing surface to harvest solar power during day time. The solar power will be stored in a battery and used for propulsion. This will allow the vehicle to change position or climb. The stored power will also be used for scientific measurements, communication and heating of critical components during night time.

4. Vehicle Conceptualization

4.1. Vehicle Geometry

As already mentioned, neutrally buoyant inflatable aerodynamic vehicles, such as the (VAMP) by Northrop Grumman [29], are promising for exploring the Venus atmosphere. Their compactness makes them potentially easier to store in rocket fairings and aeroshells. Planetary entry, inflation, and insertion in the target altitude are not considered in this paper. Inflatable aircraft have a long history that goes back to the Goodyear Inflatoplane of the 1950s. We based our design on the Stingray aircraft of the 1990s by Prospective Concepts [19,20]. The Stingray has a maximum takeoff weight of 840 kg, a maximum airspeed of 36 m/s, a wingspan of 13 m, and a volume of 680 m 2 . Based on published photos, we approximated the platform shape and airfoil geometry (Figure 4) and replicated the shapes using the computer-aided design (CAD) software SolidWorks 2021. The airfoil has 23% relative thickness which makes the vehicle voluminous. Artistic impressions of our CAD model are provided in Figure 5. The solar panels are highlighted in black and the propeller is modelled by a disk. Similar to the Prospective Concepts [20] and Festo [21] vehicles the drawings show chord-wise gas pockets that are aligned in the span-wise direction. The necessary supply of buoyant gas could be stored in a carbon-fibre gas tank.
Regarding the internal structure, larger structural members will be based on the tensairity concept (inflatable beams) [30]. A preliminary outline of the main structural members and control surfaces is shown in Figure 6. The ribs will be made from pockets that will be filled with buoyant gas similar to the chambers on parafoils. This will help maintain the airfoil shape of the inflatable wing during flight. The design of the gas pockets can be part of a later more detailed development effort.

4.2. Skin Materials

The skin will be made from a 0.025 mm thick fluorinated ethylene propylene (FEP) layer that is bonded to a 0.1 mm thick Kapton layer. A skin stress analysis shows that this bi-laminate film can withstand the expected maximum stresses without bursting. A similar film was proposed by Hall et al. [6]. To reduce absorption and provide oxidation resistance, the FEP layer is metallized on the back with silver and Inconel. The film has a low absorption and high emission coefficient. With a surface density area of 196 g/m 2 [31,32], the film is relatively lightweight. The properties of the film are summarized in Table 2.

4.3. Vehicle Sizing

The vehicle volume, V, and surface area, S f , were obtained from the CAD model for different wingspans, b. A cubic regression for the volume,
V = 2.09 × 10 5 0.00025 b 5.73 × 10 6 b 2 + 0.053 b 3 ,
and a quadratic regression for the surface area,
S f = 0.00051 0.00052 b + 1.16 b 2 ,
were then used to approximate the data points (Figure 7).
The mass of the displaced atmospheric gas is,
m a = V ρ a ,
where ρ a = 0.8868 kg/m 3 is the density of the atmosphere at 55.3 km. The density of the buoyant gas, ρ g , is obtained using the ideal gas law,
ρ g = p g R s , g T g ,
where p g and T g are the pressure and temperature. We chose helium with a gas constant of R s , g = 2076.9 J/(kg K) as buoyant gas. The mass of the buoyant gas is
m g = V ρ g .
The mass of the vehicle’s skin or film is obtained by multiplying the skin surface area, S f , given by Equation (2), with the film surface density, m/A = 196 g/m 2 (Table 2). The mass of the film is then
m f = S f × 2 m A .
In Figure 8, m a (mass of displaced atmospheric gas), m f (skin mass), m g (buoyant gas), and m a m f m g (available mass for structure, instruments, battery, motor, etc.) are plotted versus the wingspan, b. For this analysis, p g and T g were taken as the atmospheric values, p a = 50.1 kPa and T a = 299.1 K at z = 55.3 km. Because of the different scaling of volume ( b 3 ) and wing area ( b 2 ), the available mass, m a m f m g , increases sharply with wingspan.
During powered flight with horizontal airspeed u, the stagnation pressure must not exceed the internal pressure to prevent deformation of the wing shape,
p g p a + 1 2 ρ a u 2 .
Keeping the volume constant, the internal gas pressure can be raised by increasing the gas temperature or total mass (see Equation (4)). For the latter case, the available mass will be reduced. These trends are illustrated in Figure 9 for a wingspan of b = 25 m.
From a mission perspective, a low skin mass is desirable as it allows for a larger payload mass or increased buoyancy and thus extended operating altitude range. Based on the previous analysis, we decided on a wingspan of b = 25 m. For this wingspan, the skin mass is m f = 141.7 kg and the available mass at cruise altitude is 530.1 kg for a total mass of 671.8 kg. This mass is similar to that of earlier missions, such as Venera 7 which had a launch mass of 1180 kg and a landing mass of 500 kg. The volume and surface area for this wingspan are V = 833.4 m 3 and S f = 722 m 2 . The wing aspect ratio is 1.73. The maximum cruise velocity was set to 30 m/s, which is close to the nominal value for the VAMP [29]. From Figure 8, the mass of the displaced gas is m a = 739 kg, and the mass of the buoyant gas is m g = 67.2 kg.

4.4. Navigation, Communication, and Science Payload

Data published on the internet and in relevant reports [5,35,36,37] was utilized to estimate the mass of the different components. Necessary devices will include an on-board computer, an inertial navigation system, aerodynamic sensors, and equipment for communicating with a mother-ship in orbit. The science payload will likely include a mass spectrometer, a visible imager, and other cameras. The total mass of the equipment is estimated as 20 kg and the maximum power consumption is estimated as 55 W.

4.5. Solar Panel and Battery Sizing

Assuming that during night time only the navigation and communication instruments are active, the power consumption (on-board computer, inertial navigation system, aero-sensors, radio-communication devices, etc.) is roughly ≈16.7 W [5,35,36,37]. For transit over the dark side of the planet of ≈89 h (see Section 3), this will require E = 5.36 MJ of energy. Assuming that all instruments as well as the science payload are active during daytime, the power will be 55 W [5,35,36,37] and the energy consumed over 89 h will be E = 17.6 MJ. Consequently, a minimum of 23 MJ of energy must be generated while transiting the sunlit side of Venus to carry out a basic balloon-like science mission without powered flight. Since the intent is to actively explore aspects of the atmosphere, provisions are made to harvest a total of 68.9 MJ to allow for powered flight. The power harvested by solar cells is,
P s = η s C s S I sin [ m i n ( θ , π ) ] ,
where S is the wing platform area (approximated here as S S f / 2 = 361 m 2 ), C S is the percentage of the wing surface covered by solar cells, and I = 312 W/m 2 is the average irradiation for the chosen cruise altitude (from Figure 1b). The angle θ is the sun angle with respect to the solar panel. When assuming that the bank and pitch angle are zero and that the solar panel surface is flat, the sun angle only depends on the longitude,
θ = ω t ,
where t is the elapsed time. The angular velocity is,
ω = u R p l a n e t + z ,
where R p l a n e t = 6051.8 km and z = 55.3 km are the radius of Venus and the cruise altitude, respectively. For an assumed wind speed of u w = 60 m/s, an angular velocity of ω = 9.82 × 10 6 rad/s is obtained. The axial tilt of Venus is only 3 degrees and therefore not considered in this analysis. In Figure 10, the harvested solar power is plotted versus time. Time was normalized by the orbit period, T orbit = 2 π / ω .
The energy harvested during the transit through the sunlit side is obtained by integration of Equation (8),
E = 0 T / 2 P s d t = 2 η s C s S I ω .
The conversion efficiency of flexible polymorphous solar cells can be as high as η s = 5 %. For an energy harvest of 68.9 MJ, Equation (11) suggests that the solar cells must cover at least 6% of the wing area, or C S = 0.06 . Some types of flexible solar cells have density areas as low as 50 g/m 2 [38]. Based on this number, the solar panel mass is estimated as 1.1 kg. To address design uncertainties as well as other issues, such as solar cell degradation and unfavourable sun incidence angles, a larger area of the wing may be covered with solar cells.
A commercial off-the-shelf 23 Ah lithium polymer (LiPo) battery with 22.2 V nominal voltage has a mass of 2.48 kg and can store approximately 510.6 Wh or 1.8 MJ of energy. A total of 38 batteries are required to collect all of the harvested energy. The total mass of the battery pack will then be 38 × 2.48 = 93.7 kg. LiPo batteries can operate over a temperature range of −20 C to 60 C. A thermal management system might be needed to ensure that the battery pack temperature remains within the operating range.

5. Thermal and Trajectory Analysis

5.1. Trajectory Model

Starting from the zero-pressure balloon equations by Carlson and Horn [7], differential equations for the altitude, film and gas temperature of pressurized balloons with flexible skin were developed. It is understood that the present analysis does at best provide guidance for the design of a buoyant flying wing. This approach was chosen for the conceptual analysis because the gas pockets of the vehicle have not yet been designed and an analysis in closed analytical form for a complex design with gas pockets is not possible.

5.1.1. Vertical Force Balance

Both buoyancy force and aerodynamic drag are acting in the vertical direction, z (positive up). The drag force is given by,
D = 1 2 ρ a C D d z d t d z d t S ,
where ρ a is the density of the outside gas and S is the platform area. As the vertical motion is normal to the platform area, the drag coefficient is C D 1 (i.e., bluff body drag). The forces acting on the vehicle result in a vertical acceleration,
( m t + C m ρ a V ) d 2 z d t 2 = g ( ρ a V m t ) 1 2 ρ a C D d z d t d z d t S ,
with total mass,
m t = m g + m f + m P + S ,
and gravitational acceleration, g = 8.7 m/s 2 . The combined payload and structure mass is m P + S = 470 kg and the virtual mass coefficient is C m = 1 [39]. The vehicle volume, V, is made up of the buoyant gas volume, V g , and the components and structure volume. To simplify the analysis, the latter is neglected, and the gas volume is computed as,
V = V g = m g R s , g T g p g .

5.1.2. Volume Change

The outside air pressure decreases with altitude. If the inside gas pressure remains constant, during a climb the pressure difference, Δ p = p g p a , will increase. To allow for an analysis in closed analytical form, the skin is assumed to be a sphere with an equivalent radius of
R = 3 V 4 π .
As part of a more detailed design, this admittedly rather crude model could be replaced by a finite-element analysis for the actual geometry with gas pockets. The Laplace law for a thin-walled sphere provides the wall stress,
σ w = Δ p R 2 t ,
based on the wall thickness, t. From Hooke’s law,
σ w = E ϵ ,
where E is Young’s modulus and
ϵ = Δ r R ,
is the strain. For the present analysis, the Young’s modulus was taken as E = 2.08 GPa (average for FEP and Kapton (see Section 6)). Equations (17)–(19) can be combined to
E Δ r R = Δ p R 2 t .
Then, the change in radius due to the pressure difference is
Δ r = Δ p R 2 2 E t ,
and the equivalent volume is
V = 4 3 π ( R + Δ r ) 3 .

5.1.3. Heat Balance for Film

The film temperature, T f , changes according to
m f c p , f d T f d t = d Q f d t ,
where m f and c p , f are the mass and heat capacity of the film, respectively, and
d Q f d t = S I α c ( 1 + r e ) incoming solar radiation and albedo + 2 ϵ σ ( T g 4 T f 4 ) + 2 h g ( T g T f ) radiative and convective heat transfer with buoyant gas + 2 h a ( T a T f ) 2 ϵ σ T f 4 convective heat transfer with atmosphere and radiation to outside + ϵ σ T B B 4 Venus radiation
is the heat flux to and from the film. Positive terms in Equation (24) signify incoming thermal fluxes and negative terms signify outgoing thermal fluxes. In Equation (24), I is the solar irradiation (Figure 1b), α c = 0.09 is the absorption coefficient, r e = 0.689 is the albedo, ϵ = 0.48 is the infrared radiation emissivity, σ = 5.67 × 10 8 W/(m 2 K 4 ) is the Stefan–Boltzmann constant, h g and h a are the convective heat transfer coefficients from the film to the buoyant gas and the atmosphere, respectively, and T B B = 226.6 K is the black body temperature of Venus. All terms for which both the top and bottom surface contribute, scale with two times the platform area. The effect of the solar panel on the absorption and emission coefficient is neglected.

5.1.4. Convective Heat Transfer

During unpowered buoyant flight, the forced convective heat transfer is assumed to be zero and only natural convection is considered. The heat transfer coefficient from the film to the outside atmosphere is given by [40],
h a = k a N u 2 R .
The CO 2 thermal conductivity is obtained from Sutherland’s law [41],
k a = k 0 T a T 0 3 / 2 T 0 + S k T a + S k ,
with k 0 = 0.0146 W/(mK), T 0 = 273 K, and S k = 1800 K [41]. In addition, R is the radius associated with the vehicle volume (see Equation (16)). Churchill et al. [42] provided a relationship for the Nusselt number,
N u = 2 + 0.589 R a 1 4 1 + 0.43 P r a 9 16 4 9 .
The Prandtl number for CO 2 is P r a = 0.726 and the Rayleigh number is,
R a = G r P r a ,
with Grashof number,
G r = g β ( T f T a ) ( 2 R ) 3 ν a 2 .
The thermal expansion coefficient is computed as β = 1 / ( T f T a ) . The dynamic viscosity of CO 2 is obtained from Sutherland’s law,
μ a = μ 0 T a T 0 3 / 2 T 0 + S μ T a + S μ ,
with μ 0 = 1.370 × 10 5 Ns/m 2 , S μ = 222 K, and T 0 = 273 K [41]. The kinematic viscosity is obtained from ν a = μ a / ρ a . The thermal conductivity and kinematic viscosity of helium are obtained from curve fits of published data over a relevant range of temperatures. The Prandtl number is taken as P r g = 0.71 and held constant.

5.1.5. Heat Balance for Lifting Gas

The temperature of the buoyant helium gas, T g , changes according to
d Q g d t = S 2 ϵ σ ( T g 4 T f 4 ) 2 h g ( T g T f ) radiative and convective heat transfer with film .
For the analysis it is assumed that the vehicle is not moving in the horizontal direction relative to the atmosphere. Therefore, in accordance with analyses for zero-pressure balloons (e.g., Carlson and Horn [7]) only natural convection is considered. The gas pressure is kept above the ambient pressure to allow for short segments of powered flight at a maximum airspeed of u,
p g = p a + 1 2 ρ a u 2 .
The change of the gas pressure is,
d p g = d p a + 1 2 u 2 d ρ a ,
where the density derivative can be substituted with the ideal gas law,
d p g = d p a + 1 2 ρ a u 2 d p a p a d T a T a .
Making use of the expression for the change in atmosphere hydrostatic pressure, p a = ρ a g z ,
d p g = 1 + 1 2 ρ a u 2 p a ( ρ a g d z ) 1 2 ρ a u 2 d T a T a
is obtained. The specific volume of the buoyant gas is,
1 ρ g = R s , g T g p g ,
and the change of the specific volume is,
d 1 ρ g = 1 ρ g d T g T g d p g p g .
The first law of thermodynamics for the buoyant gas can then be formulated,
d q g = c v , g d T g + p g ρ g d T g T g d p g p g ,
where c v , g is the specific heat at constant volume. Substitution of Equation (34) into Equation (37) gives,
d q g = c v , g d T g + R s , g T g d T g T g + 1 + 1 2 u 2 R s , a T a ρ a ρ g g d z R s , g T g + ρ a ρ g 1 2 u 2 R s , g T g d T a T a .
Differentiation with respect to time and rearranging terms, Equation (38) becomes,
d q g d t = c v , g d T g d t + R s , g d T g d t + 1 + 1 2 u 2 R s , a T a ρ a ρ g g d z d t + ρ a ρ g 1 2 u 2 T a d T a d t .
This equation can be further simplified by making use of c p , g = c v , g + R s , g and substitutions for the vertical velocity,
v = d z d t ,
and,
d T a d t = d T a d z d z d t = d T a d z v .
Equation (39) then becomes
d q g d t = c p , g d T g d t + 1 + 1 2 u 2 R s , a T a ρ a ρ g g v + ρ a ρ g 1 2 u 2 T a d T a d z v .
This equation is combined with Equation (30) to obtain
d q g d t = 1 m g d Q g d t = c p , g d T g d t + ρ a ρ g g + 1 2 u 2 T a g R s , a + d T a d z v .
For this approach, the amount of helium gas inside the vehicle has to be adjusted as will be shown later.
For u = 0 , Equation (43) simplifies to
1 m g d Q g d t = c p , g d T g d t + ρ a ρ g g v = c p , g d T g d t + p a M a T g p g M g T a g v .
with molar masses M a and M g . Assuming p g = p a , the expression for the zero-pressure balloon [7] is recovered,
1 m g d Q g d t = c p , g d T g d t + M a T g M g T a g v .
The present model is built on simplifying assumptions. By definition, the solar array will absorb most of the impinging sunlight. Only a small portion of the absorbed sunlight will be converted to electricity. The rest is rejected as infrared heat which will heat up the airship and gas. Inefficiencies of the components within the airship also cause heating. If no thermal system is used to reject the waste heat, it will raise the equilibrium temperature of the airship and lifting gas. These effects have to be included in a more complete analysis.

5.2. Results

To obtain the vehicle trajectory over a complete orbit, Equations (13), (23) and (43) were integrated numerically with the implicit Euler backward scheme. The computational time step was Δ t = 10 s. Two cases were investigated: For the first case (variable mass), Equation (43) was solved and for the second case (fixed mass), Equation (44) was solved. The initial conditions were z = 55 km and T g = T f = 250 K for the first case and 55.3 km and T g = T f = 235.4 K for the second case. The combined structural and payload mass was set to m P + S = 520 kg. The pressure differential to the outside was chosen large enough to allow for cruise at u = 30 m/s (see Equation (7)).

5.2.1. Case 1: Variable Mass

Results obtained from Equation (43) for one orbit are illustrated in Figure 11. The film and gas temperature are directly related to the solar irradiation (Figure 11a). Because the film is very thin, the film temperature is very close to the gas temperature. During daytime, the vehicle climbs by a small amount and, as a result, the outside atmospheric temperature, pressure and density do not change much (Figure 11b). As expected, the inside gas pressure is 0.5 ρ a u 2 above the outside atmospheric pressure. Buoyant gas has to be removed (about 3 kg per orbit) to regulate the inside gas pressure during the climb (Figure 11c). Gas also has to be removed during the descend. A practical approach for removing lifting gas is to pump it back into a gas tank. Alternatively, helium can be dumped into the atmosphere and then replenished from a gas tank which would, however, shorten the mission time. During the transit through the night, all variables remain constant.

5.2.2. Case 2: Fixed Mass

Releasing helium gas during each orbit reduces the total mission time. The cycling of helium between a reservoir tank and the vehicle requires a high-pressure pump which may not be practical or incur a mission penalty. Therefore, based on Equation (44) a scenario with constant amount of buoyant gas was analysed as well. Figure 12 shows results for one orbit.
Compared to case 1, the vehicle climbs up higher in the atmosphere (Figure 12a). As a result, the outside temperature as well as pressure and density drop more than for case 1 (Figure 12a,b). The difference between the inside and outside pressure remains large enough to satisfy Equation (7). As the vehicle climbs, the volume extends but the gas mass remains constant (Figure 12c).

6. Skin Stress Analysis

An annotated drawing of the bi-laminate skin material is provided in Figure 13. Assuming negligible surface curvature, when a force F is applied, the elongation, Δ l , of the two layers and the strain, ϵ = Δ l / l 0 , will be identical. The force F is represented as the sum of two independent stresses σ 1 and σ 2 multiplied by the corresponding cross sectional areas A 1 and A 2 ,
F = σ 1 A 1 + σ 2 A 2 .
The force F can also be represented as an average stress σ ¯ multiplied by the total cross sectional area of the bi-laminate film,
F = σ ¯ ( A 1 + A 2 ) .
By equating the two expressions for the force, the relation
σ 1 A 1 + σ 2 A 2 = σ ¯ ( A 1 + A 2 )
is obtained. With the definition of axial stress, σ = E ϵ , where E is Young’s modulus, an average Young’s modulus for the bi-laminate film can be computed,
E ¯ = E 1 A 1 + E 2 A 2 A 1 + A 2 .
To prevent bursting of the skin, the stress in each layer needs to be kept below the material yield strength. The yield stress is 12 MPa for FEP [31] and 69 MPa for Kapton [32]. Based on Equation (17), the stress in the vehicle skin, σ ¯ , was evaluated over one full orbit for both cases. The individual stresses were computed from ϵ = σ ¯ / E ¯ , σ 1 = E 1 ϵ , and σ 2 = E 2 ϵ .
We chose a 0.025 mm thick layer of FEP bonded to a 0.1 mm thick layer of Kapton. For case 1, the differential pressure, p g p a , is constant and the stress is 1.85 MPa in the FEP material and 11.1 MPa in the Kapton. These values are well below the yield stress. For case 2, the differential pressure changes and the stress reaches a maximum of 11.2 MPa for the FEP and 66.9 MPa for the Kapton (Figure 14). Both values are high but below the yield stress. Ultimately, the present analysis provides initial guidance for a more detailed design. An accurate skin stress analysis for the actual skin geometry (and internal chambers or pockets similar to parafoils) will have to rely on finite element analysis and incorporate safety factors to account for uncertainty with regard to the atmospheric conditions (e.g., wind gusts).

7. Aerodynamic Analysis

The vehicle will be equipped with a propeller for powered flight (Figure 5 and Figure 6). Aerodynamic analyses were performed with XFLR5 and VSPAERO (part of the NASA open vehicle sketch pad, OpenVSP). The atmospheric conditions were obtained from VIRA. At 55.3 km altitude, the atmospheric density is ρ a = 0.887 kg/m 3 and the dynamic viscosity is μ a = 1.49 × 10 5 Ns/m 2 . For an assumed airspeed of u = 30 m/s and mean aerodynamic chord of c = 16.11 m, the flight Reynolds number is R e = 29 × 10 6 .

7.1. Two-Dimensional Airfoil Section

The two-dimensional Stingray airfoil (see Figure 4) was analysed with XFLR5. The lift curve in Figure 15a has an almost constant slope between −10 and +5 degrees angle of attack. The zero-lift angle of attack is α −2 degrees. For α 2.5 degrees, the lift coefficient is C l 0.5 , and the drag coefficient is C d 0.006 (Figure 15b). Lift-drag polars for different Reynolds numbers are provided in Figure 15c. A low-Reynolds number effects result in a considerable drag increase for the smallest Reynolds number.

7.2. Three-Dimensional Vehicle

The geometries for the three-dimensional analysis in XFLR5 and VSPAERO are slightly different from each other (see Figure 16) but similar to that of the Stingray vehicle by Prospective Concepts [19,20]. Prospective Concepts carried out over 300 successful flight tests with the Stingray, giving us the confidence that the aerodynamics were sound.
VSPAERO provides the span-wise lift distributions, d C L / d b (Figure 17a). The lift distributions peak near the wing tips and generally increase with the angle of attack, α . A contour plot of the wall pressure coefficient, C p , for the top surface at α = 2 degrees is shown in Figure 17b. The pressure is increased at the leading edge and reduced over much of the top surface.
The lift, drag, and moment coefficient obtained from XFLR5 and VSPAERO are plotted in Figure 18. In the figure, “VLM” refers to vortex lattice method and “Panel” refers to panel method. For the inviscid method, the boundary layers are disregarded. Both XFLR5 and VSPAERO predict a zero-lift angle of attack around α 0 = −2 degrees. The lift coefficient slope and drag coefficient are smaller for VSPAERO compared to XFLR5. Both programs predict a negative moment coefficient slope, d C M / d α < 0 , which is required for the pitch stability of flying wings. Higher-fidelity computational fluid dynamics (CFD) approaches and/or wind tunnel experiments would allow us to evaluate the quality of the VSPAERO and XFLR5 results. A summary of some of the relevant aerodynamic parameters is provided in Table 3. The low lift coefficient slope can be explained by the low aspect ratio of only 1.9. The buoyancy requirement (large volume, thick airfoil, large mean aerodynamic chord, low aspect ratio) competes with the requirement for good aerodynamics (thin airfoil, large aspect ratio). The present design is a compromise between the two.

7.3. Powered Flight Analysis

The vehicle is intended to allow for dirigible powered flight. For the following analysis, the horizontal airspeed, u, is relative to the surrounding atmosphere which is moving with an unknown wind speed. Based on the results from the aerodynamic analysis (Section 7), the power required for straight and level neutrally buoyant flight at constant altitude (case 1) and the altitude and airspeed for flight with aerodynamic lift (case 2) were determined based on an available battery charge of E = 68.9 MJ. Higher performance can be achieved if the battery size is increased. The present analysis is preliminary in nature and only serves as a proof of concept. A deeper analysis that includes inefficiencies in power generation, charging and discharging, propulsion, etc., will be required for a more detailed vehicle design.

7.3.1. Case 1: Straight and Level Neutrally Buoyant Flight at Constant Altitude

For the following analysis, the vehicle is assumed to be neutrally buoyant and to cruise with a horizontal velocity, u, relative to the wind at 55.3 km. No lift is generated and as a result, the lift-induced drag is zero. To simplify the analysis, the effect of convective cooling on the gas temperature and thus buoyancy is neglected. The required power for straight and level flight at constant altitude is,
P r e q = 1 η p D u = 1 η p C D , C L = 0 1 2 ρ a u 3 S .
where η p = 0.8 is the assumed propulsive efficiency and ρ a = 0.8868 kg/m 3 is the atmospheric density at cruise altitude. The drag coefficient is normalized by the platform area, S. The viscous drag coefficient at zero lift is obtained from the VSPAERO analysis, C D , C L = 0 = 0.00684. As a “sanity check”, the drag coefficient was also computed based on an average skin-friction coefficient,
C D , C L = 0 = 2 C f r .
Based on a correlation from Corke [43],
C f r = 0.455 ( l o g 10 R e ) 2.58 ,
where C D , C L = 0 = 0.00509 is obtained for R e = 29 × 10 6 . The difference between this number and the drag coefficient from VSPAERO is ≈34%. Therefore, the present results were based on the VSPAERO drag coefficient. The vehicle can cruise until the battery is completely drained, resulting in a flight time of
t f l i g h t = E P r e q .
In Figure 19a, the required power and maximum flight time are plotted over the airspeed (i.e., the velocity relative to the wind). Flight at 30 m/s requires a power of 36.7 kW and can be sustained for around 31.2 min. The power consumption for lower airspeeds is less, thus allowing for longer flight times. For a more detailed analysis, the energy required for the acceleration to cruise speed can be included in the analysis.

7.3.2. Case 2: Altitude and Airspeed for Flight with Aerodynamic Lift

The lift force adds to the buoyancy force. The combined force is balanced by the weight of the skin, structure, and payload,
g V ( ρ a ρ g ) + 1 2 ρ a u m a x 2 C L , α = 10 degrees S = g ( m f + m P + S ) .
For this analysis, the buoyant gas density and volume of the vehicle during the case 2 night flight at 55.3 km are held constant (i.e., ρ g = 0.1033 kg/m 3 and V = 844.6 m 3 ). The maximum lift force is obtained at the maximum airspeed, u max , and maximum lift coefficient, C L , max . Since the maximum lift coefficient is unknown, the value predicted by VSPAERO for α = 10 degrees was chosen, C L , α = 10 degrees = 0.362 . For the following analysis, m P + S = 520 kg was assumed. The density of the atmosphere, ρ a , decreases with altitude. This reduces the buoyancy force and the lift force. Assuming unaccelerated flight (thrust equals drag), Equation (50) can be used to obtain the maximum airspeed, u max , if P r e q is replaced by the maximum power, P max , and the drag coefficient (base drag and lift-induced drag) for α = 10 degrees is used, C D , α = 10 degrees = 0.0277 . The power depends on the flight time via P max = E / t flight . The set of equations is solved for u max and ρ a and the corresponding altitude is taken from Figure 1a. In Figure 19b, the altitude and airspeed are plotted versus the flight time. If all stored energy is consumed over 69 min, the vehicle can sustain an airspeed of ≈30 m/s in straight and level flight at an altitude of ≈69.2 km (≈13.9 km above neutrally buoyant altitude). Lower battery discharge rates increase the flight time but reduce the ceiling and maximum airspeed. The energy required for the climb from the cruise altitude to the higher altitude was not considered and can be included in a more detailed analysis.

8. Future Analysis

The present analysis is based on many simplifying assumptions. To assess the feasibility of the proposed vehicle with greater certainty, a more detailed design and first-principles-based analysis will be required. The external skin geometry and internal chambers or pockets will need to be designed. The operating environment (wind gusts, solar radiation, etc.) and aerodynamic loads have to be characterized with greater certainty. For the thermal analysis, the solar cells have to be included. The internal structure of the vehicle has to be designed. Both the skin and internal structure under load will have to be analysed with finite element analysis. The margin for the onset of the fluid–structure interaction, which can lead to flutter of the skin during powered flight, has to be investigated. The propulsion system has to be designed and analysed. This list is by no means exhaustive. Many other areas will have to be addressed and offer exciting opportunities for future research.

9. Conclusions

A buoyant flying wing vehicle is proposed for the exploration of the upper Venus atmosphere. This vehicle shares positive attributes of both balloons (simplicity, robustness) and aircraft (manoeuvrability). The vehicle geometry is similar to that of the Stingray by Prospective Concepts. Based on a sizing analysis, the wingspan was set to 25 m. For this size, it is expected that the deflated vehicle fits into the aeroshell of a typical launcher. The available combined structure and payload mass is 520 kg. Solar cells mounted to the upper surface of the wings will produce sufficient power for operation during transit through the night and for powered flight. Differential equations for the altitude, as well as skin and gas temperature, were developed and solved numerically. If helium is pumped back into a storage tank or dumped overboard during each orbit, the altitude can be maintained better, but the mission time is shortened. Therefore, the decision was made to maintain a constant helium mass. The vehicle will then climb and descend by about 500 m during each orbit. A stress analysis shows that a skin made from a 0.025 mm thick fluorinated ethylene propylene film that is bonded to a 0.1 mm thick Kapton film will be strong enough. An aerodynamic analysis was performed to obtain the lift, drag, and moment coefficient distributions. The moment coefficient slope was found to be negative, which is a requirement for tailless flying wings. If all solar energy harvested over one day is spent, neutrally buoyant powered flight at 30 m/s can be sustained for about 31.2 min at 55.3 km. Assuming that the volume and amount of buoyant gas remain the same, the vehicle can fly at 30 m/s for 69 min at an altitude of 69.2 km. A more detailed analysis, that includes the energy required for acceleration and climb as well as the energy consumed by the instruments and thermal management, is required for sizing the battery pack and solar panels. Future work will focus on a more detailed analysis of the aerodynamics and propulsion system and wind tunnel testing of a scaled model of the inflatable wing.

Author Contributions

Conceptualization, J.R., A.G. and N.C.; methodology, J.R., A.M. and A.G.; software, J.R. and A.G; validation, J.R. and A.G.; formal analysis, J.R. and A.M.; investigation, J.R., A.M. and E.N.; resources, A.G.; data curation, N.C.; writing—original draft preparation, J.R., A.M. and A.G.; writing—review and editing, A.G. and N.C.; visualization, J.R., A.M. and E.N.; supervision, A.G. and N.C.; project administration, A.G. and N.C; funding acquisition, A.G. and N.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Aeronautics and Space Administration (NASA) through an Established Program to Stimulate Competitive Research (EPSCoR) Rapid Response Research (R3) cooperative agreement No. 80NSSC21M0172 with Paulo Oemig serving as principal investigator.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zasova, L.; Moroz, V.; Linkin, V.; Khatuntsev, I.; Maiorov, B. Structure of the Venusian atmosphere from surface up to 100 km. Cosm. Res. 2006, 44, 364–383. [Google Scholar] [CrossRef]
  2. Venus Exploration Analysis Group. Venus Technology Plan. 2014. Available online: https://www.lpi.usra.edu/vexag/documents/reports/Venus-Technology-Plan-140617.pdf (accessed on 16 April 2023).
  3. Izraelevitz, J.S.; Hall, J.L. Minimum-Mass Limits for Streamlined Venus Atmospheric Probes. J. Spacecr. Rocket. 2020, 57, 683–691. [Google Scholar] [CrossRef]
  4. Venus Exploration Analysis Group. Goals, Objectives, and Investigations for Venus Exploration. 2014. Available online: https://www.lpi.usra.edu/vexag/documents/reports/GOI-140625.pdf (accessed on 16 April 2023).
  5. Gilmore, M.; Beauchamp, P.; Lynch, R.; Amato, M. Venus Flagship Mission Planetary Decadal Study Final Report. 2020. Available online: https://science.nasa.gov/science-red/s3fs-public/atoms/files/Venus%20Flagship%20Mission.pdf (accessed on 16 April 2023).
  6. Hall, J.L.; Pauken, M.; Schutte, A.; Krishnamoorthy, S.; Aiazzi, C.; Izraelevitz, J.; Lachenmeier, T.; Turner, C. Prototype Development of a Variable Altitude Venus Aerobot. In Proceedings of the AIAA Aviation 2021 Forum, Virtual Event, 2–6 August 2021; AIAA Paper 2021–2696. AIAA: Reston, VA, USA, 2021. [Google Scholar]
  7. Carlson, L.A.; Horn, W.J. New thermal and trajectory model for high-altitude balloons. J. Aircr. 1983, 20, 500–507. [Google Scholar] [CrossRef]
  8. Agrawal, R.; Buchanan, W.P.; Arora, A.; Girija, A.P.; De Jong, M.; Seager, S.; Petkowski, J.J.; Saikia, S.J.; Carr, C.E.; Grinspoon, D.H.; et al. Mission architecture to characterize habitability of venus cloud layers via an aerial platform. Aerospace 2022, 9, 359. [Google Scholar] [CrossRef]
  9. Landis, G.; Colozza, A.; LaMarre, C. Atmospheric Flight on Venus. In Proceedings of the 40th AIAA Aerospace Sciences Meeting & Exhibit, Reno, NV, USA, 14–17 January 2002; AIAA Paper 2002–819. AIAA: Reston, VA, USA, 2002. [Google Scholar]
  10. Landis, G.A.; LaMarre, C.; Colozza, A. Venus atmospheric exploration by solar aircraft. Acta Astronaut. 2005, 56, 750–755. [Google Scholar] [CrossRef] [Green Version]
  11. Landis, G.A.; Colozza, A.J. Evaluation of Long Duration Flight on Venus; NASA/TM-2006-214452; NASA Glenn Research Center: Cleveland, OH, USA, 2006.
  12. Xiongfeng, Z.; Zheng, G.; Zhongxi, H. Sun-seeking eternal flight solar-powered airplane for Venus exploration. J. Aerosp. Eng. 2015, 28. [Google Scholar] [CrossRef]
  13. Herkenhoff, B.; Fisher, J.; Hassanalian, M. Performance Analysis of Solar Fixed-wing Drones on Venus. In Proceedings of the 18th Meeting of the Venus Exploration Analysis Group (VEXAG), Virtual Event, 16–17 November 2020. [Google Scholar]
  14. Acosta, G.A.; Grow, D.; Hassanalian, M. Design and Analysis of a High Pressure and High Temperature Fixed Wing Space Drone for Venus Exploration. In Proceedings of the AIAA Aviation 2019 Forum; Dallas, TX, USA, 17–21 June 2019, AIAA Paper 2019–3465; AIAA: Reston, VA, USA, 2019. [Google Scholar]
  15. David, L. A Venus Aircraft Could Be in NASA’s Plans. 2018. Available online: https://www.space.com/41033-venus-atmosphere-aircraft-nasa-contract.html (accessed on 29 June 2018).
  16. Bullock, M.; Elston, J.; Stachura, M.; Lebonnois, S. Long Duration In Situ Science in Venus’ Clouds Enabled by Dynamic Soaring. In Proceedings of the 18th Meeting of the Venus Exploration Analysis Group (VEXAG), Virtual Event, 16–17 November 2020. [Google Scholar]
  17. Young, L.A.; Aiken, E.W.; Derby, M.; Johnson, J.; Navarrete, J.; Klem, J.; Demblewski, R.; Andrews, J.; Torres, R. Engineering Studies into Vertical Lift Planetary Aerial Vehicles. In Proceedings of the American Helicopter Society International Meeting on Advanced Rotorcraft Technology and Life Saving Activities, Montreal, QC, Canada, 11–13 June 2002. [Google Scholar]
  18. Agis, F. Venus Inflatable Airplane Concept. 2015. Available online: https://wordlesstech.com/venus-inflatable-airplane-concept/ (accessed on 14 May 2015).
  19. Breuer, J.; Ockels, W.; Luchsinger, R. An Inflatable Wing Using the Principle of Tensairity. In Proceedings of the 48th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference, Honolulu, HI, USA, 23–26 April 2007; AIAA Paper 2007–2117. AIAA: Reston, VA, USA, 2019. [Google Scholar]
  20. Flying Stingray, Handouts 4/11/03 English, Prospective Concepts Company. Available online: https://www.prospective-concepts.ch/ (accessed on 9 March 2023).
  21. Festo Air Ray. A Remote-Controlled Hybrid Construction with Flapping-Wing Mechanism. Festo Company. Available online: https://www.festo.com/us/en/e/about-festo/research-and-development/bionic-learning-network/highlights-from-2006-to-2009/air-ray-id_33851/ (accessed on 9 March 2023).
  22. Tomasko, M.; Doose, L.; Smith, P.H.; Odell, A. Measurments of the Flux of sunlight in the atmosphere of Venus. J. Geophys. Res. 1980, 85, 8167–8186. [Google Scholar] [CrossRef]
  23. Moroz, V. Stellar magnitude and albedo data of Venus. In Venus; The University of Arizona Press: Tucson, AZ, USA, 1983; pp. 27–35. [Google Scholar]
  24. Titov, D.V.; Bullock, M.A.; Crisp, D.; Renno, N.O.; Taylor, F.W.; Zasova, L.V. Radiation in the Atmosphere of Venus. Explor. Venus Terr. Planet 2007, 176, 121. [Google Scholar] [CrossRef] [Green Version]
  25. NASA. Venus Fact Sheet. 2021. Available online: https://nssdc.gsfc.nasa.gov/planetary/factsheet/venusfact.html (accessed on 20 April 2022).
  26. NASA. Earth Fact Sheet. 2021. Available online: https://nssdc.gsfc.nasa.gov/planetary/factsheet/earthfact.html (accessed on 2 November 2022).
  27. Gilmore, D. (Ed.) Spacecraft Thermal Control Handbook; AIAA: Reston, VA, USA, 2002; Volume 1, Chapter 2. [Google Scholar] [CrossRef]
  28. Colozza, A.; Landis, G. Solar Powered Flight on Venus; NASA/CR-2004-213052; NASA Glenn Research Center: Cleveland, OH, USA, 2004.
  29. Lee, G.; Polidan, R.S.; Ross, F. Venus Atmospheric Maneuverable Platform (VAMP)—A Low Cost Venus Exploration Concept. In Proceedings of the American Geophysical Union, Fall Meeting 2015, San Francisco, CA, USA, 14–18 December 2015. [Google Scholar]
  30. Luchsinger, R.H.; Pedretti, M.; Reinhard, A. Pressure induced stability: From pneumatic structures to Tensairity. J. Bionic Eng. 2004, 1, 141–148. [Google Scholar] [CrossRef]
  31. Van Nuenen, P. FEP Datasheet and FEP Properties. 2019. Available online: https://fep-film.com/fep-datasheet-fep-properties/ (accessed on 2 November 2022).
  32. DuPont. DuPont™ Kapton® PST. Available online: https://www.dupont.com/content/dam/dupont/amer/us/en/ei-transformation/public/documents/en/EI-10150-Kapton-PST-Data-Sheet.pdf (accessed on 2 November 2022).
  33. Matweb. Chemours 100 FEP Fluorinated Ethylene Propylene. Available online: https://www.matweb.com/search/datasheet_print.aspx?matguid=3dbaaa8dbb114c57996acd6738a7efc1&n=1 (accessed on 2 November 2022).
  34. Matweb. DuPont™ Kapton® 500VN Polyimide Film, 125 Micron Thickness. Available online: https://www.matweb.com/search/datasheet_print.aspx?matguid=338573ad1bdf4586aa17fab95f3a57d7 (accessed on 2 November 2022).
  35. Parker Hannifin Corporation. MicroStrain Sensing Product Datasheet: 3DM-GX5-GNSS/INS. Available online: https://www.microstrain.com/inertial-sensors/3dm-gx5-45 (accessed on 2 November 2022).
  36. Innovative Sensor Technology IST AG. Application Note RTD Platinum Sensor. Available online: https://www.farnell.com/datasheets/2340312.pdf (accessed on 2 November 2022).
  37. Emerson Automation Solutions. Paine Miniature-Satellite Series Pressure Transmitter. Available online: https://www.emerson.com/documents/automation/product-data-sheet-212-miniature-space-satellite-pressure-transmitter-paine-en-80396.pdf (accessed on 2 November 2022).
  38. Pagliaro, M.; Ciriminna, R.; Palmisano, G. Flexible solar cells. ChemSusChem Chem. Sustain. Energy Mater. 2008, 1, 880–891. [Google Scholar] [CrossRef]
  39. Brennen, C.E. A Review of Added Mass and Fluid Inertial Forces; Department of the Navy: Port Hueneme, CA, USA, 1982. Available online: https://resolver.caltech.edu/CaltechAUTHORS:BREncel82 (accessed on 15 April 2023).
  40. Holman, J. Heat Transfer, 10th ed.; McGraw Hill: New York, NY, USA, 2010; pp. 346–347. [Google Scholar]
  41. White, F.M. Viscous Fluid Flow; McGraw-Hill: New York, NY, USA, 2006. [Google Scholar]
  42. Churchill, S.W.; Chu, H.H. Correlating equations for laminar and turbulent free convection from a vertical plate. Int. J. Heat Mass Transf. 1975, 18, 1323–1329. [Google Scholar] [CrossRef]
  43. Corke, T.C. Design of Aircraft; Prentice Hall: Hoboken, NJ, USA, 2003; pp. 75–76. [Google Scholar]
Figure 1. (a) Venus atmosphere density, pressure, and temperature. (b) Solar irradiation.
Figure 1. (a) Venus atmosphere density, pressure, and temperature. (b) Solar irradiation.
Aerospace 10 00534 g001
Figure 2. Vehicle circumnavigation of the planet.
Figure 2. Vehicle circumnavigation of the planet.
Aerospace 10 00534 g002
Figure 3. Anticipated mission profile.
Figure 3. Anticipated mission profile.
Aerospace 10 00534 g003
Figure 4. Stingray airfoil.
Figure 4. Stingray airfoil.
Aerospace 10 00534 g004
Figure 5. Venus vehicle front and side view.
Figure 5. Venus vehicle front and side view.
Aerospace 10 00534 g005
Figure 6. Conceptual design of internal structure and control surfaces.
Figure 6. Conceptual design of internal structure and control surfaces.
Aerospace 10 00534 g006
Figure 7. Vehicle volume and surface area as a function of wingspan.
Figure 7. Vehicle volume and surface area as a function of wingspan.
Aerospace 10 00534 g007
Figure 8. Mass analysis.
Figure 8. Mass analysis.
Aerospace 10 00534 g008
Figure 9. Available mass as a function of airspeed and gas temperature for a 25 m span.
Figure 9. Available mass as a function of airspeed and gas temperature for a 25 m span.
Aerospace 10 00534 g009
Figure 10. Harvested solar power over one orbit for C s = 0.2 .
Figure 10. Harvested solar power over one orbit for C s = 0.2 .
Aerospace 10 00534 g010
Figure 11. Vehicle trajectory and properties for case 1 (variable mass). (a) Altitude and temperatures. (b) Pressures and densities. (c) Buoyant gas volume and mass.
Figure 11. Vehicle trajectory and properties for case 1 (variable mass). (a) Altitude and temperatures. (b) Pressures and densities. (c) Buoyant gas volume and mass.
Aerospace 10 00534 g011
Figure 12. Vehicle trajectory and properties for case 2 (fixed mass). (a) Altitude and temperatures. (b) Pressures and densities. (c) Buoyant gas volume and mass.
Figure 12. Vehicle trajectory and properties for case 2 (fixed mass). (a) Altitude and temperatures. (b) Pressures and densities. (c) Buoyant gas volume and mass.
Aerospace 10 00534 g012
Figure 13. Bi-laminate film with tensile stress.
Figure 13. Bi-laminate film with tensile stress.
Aerospace 10 00534 g013
Figure 14. Stress in skin layers for second case.
Figure 14. Stress in skin layers for second case.
Aerospace 10 00534 g014
Figure 15. (a) Lift and (b) drag coefficient. (c) Lift-drag polar.
Figure 15. (a) Lift and (b) drag coefficient. (c) Lift-drag polar.
Aerospace 10 00534 g015
Figure 16. Geometries for XFLR5 and VSPAERO analysis.
Figure 16. Geometries for XFLR5 and VSPAERO analysis.
Aerospace 10 00534 g016
Figure 17. (a) Lift coefficient per span and (b) top surface wall pressure coefficient for α = 2 degrees.
Figure 17. (a) Lift coefficient per span and (b) top surface wall pressure coefficient for α = 2 degrees.
Aerospace 10 00534 g017
Figure 18. (a) Lift, (b) drag, and (c) moment coefficient.
Figure 18. (a) Lift, (b) drag, and (c) moment coefficient.
Aerospace 10 00534 g018
Figure 19. (a) Power required and maximum flight time for neutrally buoyant flight at 55.3 km. (b) Altitude and airspeed for powered flight with aerodynamic lift.
Figure 19. (a) Power required and maximum flight time for neutrally buoyant flight at 55.3 km. (b) Altitude and airspeed for powered flight with aerodynamic lift.
Aerospace 10 00534 g019
Table 1. Atmospheric conditions and solar irradiation for 55.3 km altitude.
Table 1. Atmospheric conditions and solar irradiation for 55.3 km altitude.
Temperature, T299 K
Pressure, p50.1 kPa
Density, ρ 0.887 kg/m 3
Solar irradiation, I318 W/m 2
Table 2. Properties of bi-laminate film.
Table 2. Properties of bi-laminate film.
Surface density, m / A [31,32]196 g/m 2
Absorption coefficient, α c [6]0.09
Emission coefficient, ϵ [6]0.48
Specific heat, c p [33,34]1047 J/(kg K)
Table 3. Aerodynamic parameters from VSPAERO.
Table 3. Aerodynamic parameters from VSPAERO.
Zero lift angle of attack, α 0 −1.906 d e g
Drag coefficient at zero lift, C D , C L = 0 0.00684
Lift coefficient at α = 10 degrees, C L , α = 10 degrees 0.362
Drag coefficient at α = 10 degrees, C D , α = 10 degrees 0.0277
Lift coefficient slope, d C L / d α 0.0399
Moment coefficient slope, d C M / d α −0.0129
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rosales, J.; Miller, A.; Nunez, E.; Gross, A.; Chanover, N. Conceptual Design of Hybrid Aerial Vehicle for Venus Exploration. Aerospace 2023, 10, 534. https://doi.org/10.3390/aerospace10060534

AMA Style

Rosales J, Miller A, Nunez E, Gross A, Chanover N. Conceptual Design of Hybrid Aerial Vehicle for Venus Exploration. Aerospace. 2023; 10(6):534. https://doi.org/10.3390/aerospace10060534

Chicago/Turabian Style

Rosales, Jesus, Addison Miller, Edgar Nunez, Andreas Gross, and Nancy Chanover. 2023. "Conceptual Design of Hybrid Aerial Vehicle for Venus Exploration" Aerospace 10, no. 6: 534. https://doi.org/10.3390/aerospace10060534

APA Style

Rosales, J., Miller, A., Nunez, E., Gross, A., & Chanover, N. (2023). Conceptual Design of Hybrid Aerial Vehicle for Venus Exploration. Aerospace, 10(6), 534. https://doi.org/10.3390/aerospace10060534

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop