Acessibilidade / Reportar erro

Wind Tunnel Testing on a Generic Model of a Hybrid Lifting Hull

ABSTRACT

In this research, an experimental investigation was carried out at the International Islamic University Malaysia - Low Speed Wind Tunnel facility on a generic model of a hybrid lifting hull. Based on the historical trends of non-rigid airships, the fineness ratio of the said hull has been selected equal to 4. Free stream velocity was kept at 20 m/s and, along with the estimation of aerodynamic parameters, longitudinal and lateral stability characteristics were determined over a range of angles of attack from −8° to +12° and angles of sideslip from −10° to +10°. Zero lift coefficient was obtained at −4.2°, and the corresponding value was found to be greater than that at zero angle of attack. The comparison of the experimental results with the existing analytical relationships of wing has revealed that such an airfoil shaped hull cannot be considered as a wing due to 37% less analytical value of lift coefficient than that obtained by CFD simulations of the said hull. Existing equation of form factor of hull for conventional airships was also revisited, and a correction factor equal to 1.16 in the fundamental drag equation of aircraft’s fuselage was also proposed for fineness ratio equal to 4. Trends of the experimental data and mabsmoraiscomparison of the same with the theoretical calculations and computational results posed some interesting findings. The longitudinal and directional stabilities of a hybrid lifting hull were found to be statically unstable.

KEYWORDS
Aerodynamics; Hybrid airship; Fineness ratio; Static stability; Wind tunnel testing

INTRODUCTION

The concept of lifting fuselage for aircraft and hybrid lifting hull for hybrid airship is derived from nature as a few marine animals do generate aerodynamic lift from the body (Ul Haque et al. 2016aUl Haque A, Asrar W, Omar AA, Sulaeman E, Mohamed Ali JS (2016a) Hydrodynamic contour of Steller sea lion — an inspiration for future hybrid buoyant aircraft. Anim Rev 3(1):1-9. doi: 10.18488/journal.ar/2016.3.1/101.1.1.9
https://doi.org/10.18488/journal.ar/2016...
). Vogel (1994)Vogel S (1994) Life in moving fluids: the physical biology of flow. Princeton: Princeton University Press. was among the first who noticed that, if one looks at a housefly or fruit fly from the side, the head, thorax and abdomen also seem to form a non-symmetrical airfoil-flatter on the bottom, more rounded on top. Regarding marine animals, there are some cases in which there is a requirement of negative lift (Vogel 2013Vogel S (2013) Comparative biomechanics: life’s physical world. Princeton: Princeton University Press.). For example, ducks have enough air in their plumage so they are awkwardly buoyant and may need negative lift (Prange and Schmidt-Nielsen 1970Prange HD, Schmidt-Nielsen K (1970) The metabolic cost of swimming in ducks. J Exp Biol 53:763-777.). Recently, Ul Haque et al. (2015a)Ul Haque A, Asrar W, Omar AA, Sulaeman E, Mohamed Ali JS (2015a) Cambered profile of a California sea lion’s body. J Exp Biol 218:1270-1271. doi: 10.1242/jeb.117556
https://doi.org/10.1242/jeb.117556...
argued that the lift generated by the body is perhaps free of cost lift and can be utilized if such a marine animal want to swim at constant level/height in the sea. Additional lift is also required for those flight segments like coming out from water, sharp turns etc.

In the field of aviation, it is not a new concept and, as per the review paper by Wood and Bauer (2011)Wood RM, Bauer SXS (2011) Flying wings/flying fuselages. Proceedings of the 39th Aerospace Sciences Meeting and Exhibit; Reno, USA., Vincent Burnelli has earlier designed a number aircraft based on the lifting fuselage concept. Perhaps, the lifting body designs usually have greater operational flexibility as the adjustment of dynamic lift can accommodate changes in vehicle’s weight due to the burning of fuel. In an earlier reference, Wood (2003)Wood RM (2003) The contributions of Vincent Justus Burnelli. Proceedings of the 41st AIAA Aerospace Sciences Meeting and Exhibit; Reno: USA. mentioned that, during the period of 1920 to 1955, about 57 aircraft were developed, but these were either single or 2-seat. Other hybrid concepts have aerodynamic lift coming from the wing as well as from the fuselage. This concept was later applied by Santos Dumont (1973)Santos Dumont A (1973) My airships; the story of my life. New York: Dover Publications. (Original work published in 1904). for the design of a winged hybrid airship in which partial aerodynamic lift is obtained from the aerodynamic contour of hull; he was the first to prove this concept by flight testing.

According to Becker (1958)Becker JV (1958) Preliminary studies of manned satellites. Winged configurations. Proceedings of the NACA Conference on High-Speed Aerodynamics; Hampton, USA., if we look back in the history, the first lifting-body concepts involved very blunt half cones. Later, the concepts evolved by Saltzman et al. (1999)Saltzman EJ, Wang KC, Iliff KW (1999) Flight-determined subsonic lift and drag characteristics of seven lifting-body and wing-body reentry vehicle configurations with truncated bases. Proceedings of the 37th Aerospace Sciences Meeting and Exhibit; Reno, USA. into higher fineness-ratio cones to achieve the capability of an unpowered horizontal landing. Numerous wind tunnel tests were performed on candidate versions of the half cones and shapes having flattened bottom surfaces. In 1962, Reed and Darlene (1977)Reed R, Darlene L (1977) Wingless flight: the lifting body story. NASA SP-4220; [accessed 2016 Aug 26]. http://history.nasa.gov/SP4220/contents.htm
http://history.nasa.gov/SP4220/contents....
have provided the details of an unpowered horizontal landing and controllable flight with a miniature lightweight-radio-controlled model of an M2 half-cone configuration. Unfortunately, none of them has explored the airfoil shaped hull from the aerodynamic and stability point of view.

Several studies can be found in the literature starting from the early 1920s, such as Munk (1924)Munk MM (1924) The aerodynamic forces on airship hulls. NACATR-184. and Rizzo (1924)Rizzo F (1924) A study of static stability of airships. NACA Technical Note 204., and the late 1970s and early 1980s, like DeLaurier and Schenck (1979)DeLaurier J, Schenck D (1979) Airship dynamic stability. Proceedings of the AIAA Lighter-Than-Air Systems Technology Conference; Palo Alto, USA. and Tischler et al. (1981)Tischler MB, Jex HR, Ringland RF (1981) Simulation of Heavy Lift Airship dynamics over large ranges of incidence and speed. Proceedings of the Lighter-Than-Air Systems Technology Conference; Annapolis, USA., as well as many others, about the aerodynamic and stability characteristics of airships. However, in the case of hybrid airships, there is limited experimental data available for aerodynamic lift generated by the aerodynamic lifting profile of hull (alone) and its stability characteristics. Most of the experimental data is related to the lifting bodies of hypersonic vehicles. For example, some experimental data can be found in an old NACA report by Gumse (1967)Gumse B (1967) Low-speed wind-tunnel tests of a full-scale M2-F2 lifting body model. NASA TM X-1347., which is based on W-F2 configuration. As per Ash (1972)Ash LG (1972) Flight test and wind tunnel performance characteristics of the X-24A lifting body; [accessed 2016 Aug 26]. http://oai.dtic.mil/oai/oai?verb=getRecord&metadataPrefix=html&identifier=AD0901465
http://oai.dtic.mil/oai/oai?verb=getReco...
, this was a modified design of M2-F2 with modifications carried out for the afterbody, the control surfaces, and the canopy location. Others include this on NASA HL-20 to develop a preliminary subsonic aerodynamic model for simulation studies of its lifting body Jackson and Christopher (1992)Jackson EB, Christopher IC (1992) Preliminary subsonic aerodynamic model for simulation studies of the HL-20 lifting body. NASA TM-4302.. Nevertheless, all of these bodies have fineness ratio (λ) greater than 6 and wings blended with the fuselage.

Hybrid airships are among the potential candidates for tourism industry and transportation of agricultural products. Such vehicles also have the potential to reduce the aviation transport gas emissions and many have come up with promising conceptual designs like Aeroscraft (2013)Aeroscraft (2013) Launching a new century of aviation innovation; [accessed 2013 Apr 2]. http://www.aerosml.com/imagegallery.html
http://www.aerosml.com/imagegallery.html...
, the airship Sky Freighter (Millennium Air Ship Inc. 2012Millennium Air Ship Inc. (2012) The airship “Sky Freighter”; [accessed 2012 May 8]. http://www.millenniumairship.com/products.html
http://www.millenniumairship.com/product...
) and Lockheed’s LEMV (Harrison 2010Harrison J (2010) Lockheed’s LEMV consolation; [accessed 2012 May 8]. http://edgefighter.com/2010/06/17/lockheedslemvconsolation
http://edgefighter.com/2010/06/17/lockhe...
). Unfortunately, the aerodynamics and stability characteristics of airfoil shaped hulls (alone) of airships have not been fully explored yet. One of the probable reasons is due to the non-availability of the experimental data which can reveal its aerodynamic behavior. In the present study, an effort was done to fill this gap by carrying out an experimental study on a generic model of a hybrid lifting hull (HLH). The selection of its λ value was done such that it falls within the known range of non-rigid airships, discussed in detail in the following section.

SELECTION OF FINENESS RATIO

HLH is a type of unconventional hull with voluminous volume for the buoyant lift. It is usually designed to be partially supported by buoyancy lift generated by the buoyant gas while the remaining weight is held up by the aerodynamic lift generated by the aerodynamic contour of the hull (Trenkle 2014Trenkle K (2014) Hybrid Certification Criteria (HCC) for transport category hybrid airships document number: 1008D0122. Revision: C; [accessed 2014 Jun 4]. https://www.federalregister.gov
https://www.federalregister.gov...
). However, the selection of λ is found to be quite trivial, especially in a scenario where there is no guideline/data bank from the historical trends. From the certification point of view, most of the certified airships are flexible in terms of structural anatomy, i.e. nonrigid airships. But from an aerodynamic point of view, high value of λ is always desired. However, a low value is desired for optimum structural weight Tanaka et al. (2005)Tanaka K, Hamaguchi Y, Maekawa S, Ishikawa T (2005) Structures and advanced composites for lighter-than-air aircraft. Proceedings of the 56th International Astronautical Congress; Fukuoka, Japan.. If we look at the historical trends of non-rigid airships (Table 1), the lowest value of λ was found in the first winged airship, which was designed and flown by Santos Dumont (1973)Santos Dumont A (1973) My airships; the story of my life. New York: Dover Publications. (Original work published in 1904).; DM-20, on the other hand, has the highest value.

Table 1
Fineness ratio of non-rigid airships.

In the present study, the aerodynamic contour of the generic model of HLH was designed by using RONCZ 1082 airfoil till 70% of the chord length, and the λ value was set at 4. The length and the diameter of the hull were 0.56 and 0.135 m, respectively. Against this λ value, the form factor of airship and aircraft are a bit close to each other. But as per the findings of Tanaka et al. (2005)Tanaka K, Hamaguchi Y, Maekawa S, Ishikawa T (2005) Structures and advanced composites for lighter-than-air aircraft. Proceedings of the 56th International Astronautical Congress; Fukuoka, Japan., the effects of pressure drag will be there for a hull having λ < 5.

EXPERIMENTAL SETUP

In order to find the aerodynamic and static stability behavior of a HLH model, tests were conducted at the International Islamic University Malaysia - Low Speed Wind Tunnel (IIUM-LSWT) at Reynolds number (Re) equal to 6.3 × 105 against free stream velocity equal to 20 m∙s–1. IIUM-LSWT is a closed-loop wind tunnel with a test section of dimensions 1.5 × 2.3 × 6 m and a maximum achievable speed equal to 50 m∙s–1 (Hasim et al. 2008Hasim F, Rusyadi R, Surya WI, Asrar W, Omar AA, Mohamed Ali JS, Kafafy R (2008) The IIUM low speed wind tunnel. Proceedings of the 2nd Engineering Conference on Sustainable Engineering Infrastructures Development & Management (EnCon2008); Kuching, Malaysia.). The dynamic pressure in the test section varies from −0.5 to 0.4% from the plane mean value, and the flow angularity holds within ±0.2° (Wiriadidjaja et al. 2012Wiriadidjaja S, Hasim F, Mansor S, Asrar W, Rafie ASM, Abdullah EJ (2012) Subsonic wind tunnels in Malaysia: a Review. Appl Mech Mater 225:566-571. doi: 10.4028/www.scientific.net/AMM.225.566
https://doi.org/10.4028/www.scientific.n...
).

It is well known that the scaled-down models for the wind tunnel testing are either simplified or have extra bolts and nuts for attachment purposes, so this drawback of small-scale models cannot reproduce a ditto copy of the full-scale model. Therefore, the scaled-down model of HLH of aspect ratio (AR) equal to 0.25 is manufactured in a single piece and polished to get a gloss finish. Such a wooden model cannot be made hollow from side due to issues related to its attachment with the strut. In order to fulfill the requirement of blockage ratio, defined by Pope and Rae (1984)Pope A, Rae WH (1984) Low-speed wind tunnel testing. 2nd edition. New York: Wiley., the size of HLH model is kept small, and its blockage ratio including that of the strut was just 2.5% of the cross-sectional area of the test section. The model is attached with the strut in the test section by manufacturing an adopter made of stainless steel. However, similarly to any other wind tunnel model, this adopter will create an unavoidable cavity in the model. Hence, it is placed inside the hull body to avoid any additional drag due to model-strut attachment. Moreover, due to the requirement of the estimation of yaw stability parameters, the option of half-model testing is not explored.

Sign conventions are used in the present study such that positive axial force is along the positive x axis (forward from tail to nose), and positive normal force is along the positive z axis (upwards). The lift and drag coefficients (CL and CD) are estimated in the wind-axis system, and the results so obtained were plotted without any curve fitting — wing span (b) and chord length (c) equal to 0.135 and 0.54 m, respectively. The area (S) equal to 0.042 m2 was the reference parameter used for the data reduction. All the moments were obtained about the moment reference center (MRC) of the balance, which is at 0.14 m from the nose of the HLH. In the case of airships, as well as in the hybrid ones, the gondola is usually located at the base of the hull. Therefore, the offset distance in the vertical z direction was defined accordingly.

Because the data measured by balance includes the weight components of the model in a special state, it is therefore necessary to deduct the weight components to obtain the real aerodynamic coefficients. The method for deducting the influence was collecting the measurement data without wind, i.e. wind-off and model-off (W0M0) condition, and then it is deducted in the test matrix. Tests were also conducted to estimate the drag and pitching moment of the strut and fairing, referred here as wind-on and model-off (W1M0) condition. As per the guidelines provided by Tucker (1990)Tucker VA (1990) Measuring aerodynamic interference drag between a bird body and the mounting strut of a drag balance. J Exp Biol 154:439-461. to obtain the real aerodynamic coefficients for wind-on and model-on (W1M1) condition, results so obtained from wind tunnel testing were subtracted from the W1M0 and W0M0 conditions. The pictorial views of tests conducted for W1M0 and W1M1 conditions are shown in Figs. 1a and 1b, respectively, and the major dimensions of the HLH in the side and top views are shown in Figs. 2a and 2b, respectively.

Figure 1
Pictorial view of W1M0 and W1M1 tests. (a) Strut alone test; (b) Model attached inside the test section.
Figure 2
Geometric details of HLH along with its surface grid and distribution of pressure coefficient at Re = 6.3 × 105. (a) Side view; (b) Top view; (c) Surface grid; (d) Distribution of pressure coefficient over the hull’s surface.

RESULTS AND DISCUSSION

AERODYNAMICS AND STATIC STABILITY DERIVATIVES

In order to get the 1st-order approximation of aerodynamics and static stability derivatives of the HLH, the steady-state simulations were run by using ANSYS® FLUENT software for which the SIMPLE scheme is employed for pressure velocity coupling along with the k-ω SST model. The details of the geometry used for the mentioned purpose are shown in Figs. 2a and 2b. Structured mesh is generated for the said purpose with additional focus on grid refinement at the nose (Fig. 2c). No slip condition is used for the wall surface. Since the temperature problem is not of interest, the adiabatic wall conditions with no slip boundary condition are employed. Velocity magnitude, along with its direction components, is specified for the velocity inlet. Outflow boundary is also used, for which all the flow variables are basically extrapolated for incompressible flow. For subsonic flows, it is required to apply the boundary conditions on the computational domain of radius equal to 20 times the length of the HLH. Therefore, the domain extents are kept at 20 times the length of the body to avoid the flow variations near the surface. Coupled explicit solver is used with 2nd-order accuracy. All the simulations were run on an Intel® hex-core dual processor, 3.33-GHz system, with 32 GB of RAM, and the results were overlapped with those obtained from the wind tunnel testing (Fig. 3). A stagnation point is observed at the nose of the HLH, as shown from the contours of the pressure coefficient (Fig. 2d).

Figure 3
Experimental results of aerodynamic and stability coefficients of HLH. (a) Variation in CL w.r.t α; (b) Variation inCD w.r.t α; (c) Variation in Cm w.r.t α; (d) Variation in Cn w.r.t β.

For longitudinal tests (α ≠ 0, β = 0, where α is the angle of attack and β is the side-slip angle), α values were changed from −8° to +12°, and, for β sweep test, the model strut assembly remained rested on tunnel turntable and rotated at a range from −10° to +l0° with 5.-increments. The graphs in Fig. 3 show the trends of aerodynamic and static stability characteristics of HLH in longitudinal as well as in the lateral direction. The lift curve slope (CLα) value is obtained by applying the 2nd-order polynomial fit on the curve shown in Fig. 3a. Its value comes out to be equal to 0.0067/deg. This figure also shows a continuous increase in CL for a defined range of α and a reduction in its value for a negative α. Moreover, there will certainly be a stall value of α beyond 12°; however, if a wing is also attached to it, then the critical value will be that of the wing. Thus, the stalling phenomenon which gives rise to reduction in CL and dramatic rise in CD is not observed in these graphs. The angle of attack corresponding to zero-lift condition (αOL) corresponding to CLOcondition is equal to −4.2° and its corresponding CD value, i.e. zero lift drag coefficient (CDO), is equal to 0.028 (Fig. 3b). It can also be observed from this figure that the CDat zero angle of attack (0.041) was lower than that at zero-lift condition. Moreover, the induced drag is perhaps responsible for high values of CD values at a higher angle of attack.

Figure 3c illustrates the pitching moment coefficient for defined values of the angle of attack. Pitching moment coefficient exhibits the same behavior as the drag, with a small difference at negative angles of attack. In addition, the slope of pitching moment (Cma) is negative for the entire angle of attack range, indicating that it is statically unstable throughout this range. Cmα is positive and its value is 0.016/deg, indicating that the HLH model under study is longitudinally unstable, and the coefficient of pitching moment at zero lift (Cmo) is 0.137 (Fig. 3c). It is an important design parameter as it is always added to the Cmo values of other components, such that the total Cmo can be made available for the fundamental stability equation. Similarly to the CLα value, the slope of yawing moment coefficient (C) was also obtained by curve fitting and its value is equal to −0.01/deg. An anomaly was seen for the β sweep curve, i.e. an unsymmetric pattern in the value of the yawing moment coefficient (Cn) for positive and negative β was observed. One of the probable reasons is due to the asymmetry in the model, made of wood. All these results are presented in Table 1 for quick reference.

Wind tunnel testing always reflects the real flow over defined shape with the help of reliable data. Such results can be employed to develop as well as to check the applicability of known analytical relationships. Therefore, the experimental results are then compared with the existing analytical relationships, specially the slopes.

COMPARISON WITH EXISTING ANALYTICAL RELATIONSHIPS

If we consider HLH as a wing, then CLα obtained by using Eq. 1 is equal to 0.00422/deg. It is important to note that Eq. 1, taken from Raymer (2012)Raymer DP (2012) Aircraft design: a conceptual approach. 5th edition. Reston: American Institute of Aeronautics and Astronautics., will provide CLα in per radian, and, for the purpose of comparison of results, this digit has been converted into per degree.

The CLα value predicted by using the above-mentioned relationship is 0.0054/deg. This value is lesser than that predicted by the experiment, which is equal to 0.006/deg. Nevertheless, the experimental value is consistent with that obtained from the computational fluid dynamics (CFD) results. The basic difference is due to the fact that this relationship is for estimation of lift produced by the wing (a lifting surface) and which considers the wing as a flat plate. In our case, the (lifting body) hull is not a flat plate, and a wing with AR equal to 0.25 is hardly seen in any aircraft. Moreover, the analytical relationships estimated 37% less value of lift coefficient than that obtained by using CFD. Therefore, the use of analytical relationship of the wing for estimation of aerodynamic lift produced by lifting hull will not be suitable. Experimental data of space shuttles are available, but the AR for such lifting body includes the area of wings attached to it. The CLα values of M2-F2 (AR = 0.712) and space shuttle prototype (AR = 2.26) are 0.02 and 0.0446/deg, respectively (Paulson et al. 1960Paulson JW, Shanks RE, Johnson JL (1960) Low-speed flight characteristics of reentry vehicles of the glide-landing type. NASA TMX-67563.). Hence, considering the HLH as a wing planform cannot be justified for the airfoil shaped hull being investigated for airship.

CDo is an important parameter required for aerodynamic analysis and design-related studies (Roskam and Lan 1997Roskam J, Lan CTE (1997) Airplane aerodynamics and performance. Lawrence: DAR Corporation.). Its exact estimation is very important for accurate prediction of performance parameters of any flying vehicle. The analytical expression for CDo is given in Eq. 2, taken from Nicolai and Carichner (2013)Nicolai LM, Carichner GE (2013) Fundamentals of aircraft and airship design. Vol. 2: airship design and case studies. Reston: American Institute of Aeronautics and Astronautics.. This relationship is basically the function of form factor (FF). For the case of airship’s hull, it was defined with the help of Eq. 3, and, for aircraft’s fuselage, Eq. 4 is usually used.

where: Swet is the wetted area and Sref is the reference one, both in m2; Cf is the coefficient of skin friction; Q is the dynamic pressure; FR is the fineness ratio. Against the calculated value, if Cf is equal to 0.00488, considering the hull of airship as fuselage (Eq. 3), then the value of FF comes out to be 2.35, which is quite high when compared with FF equal to 1.94, obtained from Eq. 4 (Raymer, 2012Raymer DP (2012) Aircraft design: a conceptual approach. 5th edition. Reston: American Institute of Aeronautics and Astronautics.). For airship’s hull of any λ value, this author suggested an adjustment factor of 0.8 in basic drag equation to include scaling effects for the airship’s hull. For the value of Swet/Sref equal to 2.62, regarding the geometry under study, the values of CDo obtained by using Eqs. (3) and (4) are 0.029 and 0.024, respectively. The predicted value of CDo obtained from the CFD was quite high and equal to 0.048. Although, Eq. 3, which is taken from Nicolai and Carichner (2013)Nicolai LM, Carichner GE (2013) Fundamentals of aircraft and airship design. Vol. 2: airship design and case studies. Reston: American Institute of Aeronautics and Astronautics., is not really a FF by the same definition as used by Hoerner (1965)Hoerner SF (1965) Fluid-dynamic drag: theoretical, experimental and statistical information. Bakerfield: published by the author., but it has provided quite close results to the experimental value. According to a recent review of FF by Ul Haque et al. (2016b)Ul Haque A, Asrar W, Omar AA, Sulaeman E, Mohamed Ali JS (2016b) Assessment of engine’s power budget for hydrogen powered hybrid buoyant aircraft. Propulsion and Power Research 5(1):34-44. doi: 10.1016/j.jppr.2016.01.008
https://doi.org/10.1016/j.jppr.2016.01.0...
, Eq. 3 is taken from Hoerner (1965)Hoerner SF (1965) Fluid-dynamic drag: theoretical, experimental and statistical information. Bakerfield: published by the author., in which the definition of FF is actually based on a constant frontal area as reference, and the term Swet/Sref was defined to be equal to 3FR (Eq. 5). In fact, it is a simple adjustment which allows the use of Sref for the calculation of CD by using Eq. 6:

Against the value of FR equal to 4, if we use Eq. 5, then its value is 0.137. This huge digit is due to the fact that Hoerner (1965)Hoerner SF (1965) Fluid-dynamic drag: theoretical, experimental and statistical information. Bakerfield: published by the author. had derived the frontal area as reference. Moreover, based on the experimental results and its comparison study, we suggest a correction factor of 1.16 for the FF value obtained by using Eq. 4 and the planform area as reference. It also shows that, in comparison with conventional hull, an airfoil shaped hull will have more drag count.

The comparison of analytical values of Cmo and Cmα and those obtained from the experimental data is perhaps not straight forward, because the existing relationships add the effect of the wing as well in Eqs. 7 and 8. If we closely inspect these 2 equations, shape effects were considered for the value of Cmo, but, surprisingly, it has not catered for estimating the value of the slope of pitching moment of the fuselage (Cmα fus) (Ul Haque et al. 2015bUl Haque A, Asrar W, Omar AA, Sulaeman E, Mohamed Ali JS (2015b) Preliminary aerodynamic and static stability analysis for hybrid buoyant aerial vehicles at low speeds using digital DATCOM. Can Aeronaut Space J 61(3):51-60. doi: 10.5589/q16-001
https://doi.org/10.5589/q16-001...
).

where: k2 – k1 is the correction factor to account for the fuselage slenderness ratio; wf is the average width of the fuselage section, in m; αZLW is the wing zero-lift angle relative to the fuselage reference line, in degrees; Cmα fus is the pitching moment coefficient of the fuselage at zero lift; CMGC is the pitching moment coefficient about the mean geometric chord.

The pitching moment contribution of the fuselage of general aviation aircraft can be approximated by Eq. 9, taken from Raymer (2012)Raymer DP (2012) Aircraft design: a conceptual approach. 5th edition. Reston: American Institute of Aeronautics and Astronautics.. In this relationship as well, the presence of wing can be observed from the factor C and the reference area of the wing (Sw, in m2). Wf is the maximum width of the fuselage, Lf is the length and Kf is the empirical pitching moment coefficient:

In this case, we use a low-fidelity tool like Aircraft Digital DATCOM®, then the results so obtained reveal that it has underpredicted Cmo and overpredicted Cmα, when compared with the experimental and computational values. One of the probable reasons is due to the interference of strut with HLH model. On the other side, CFD underpredicted the Cmo value when compared with the experiment, which is equal to −0.187. However, the Cmα value obtained from CFD (0.0204/deg) is higher than that of the experiment, 0.0165/deg. The opposite is true when CFD results are compared with the DATCOM® ones. Similarly to CLα, Cnβ is obtained and is equal to −0.01/deg. Aircraft Digital DATCOM® does not provide the directional stability contribution at individual angles but provides the value of Cnβ, which is equal to −0.0065/deg for this case of HLH. The comparison of CL and CD is done as Aircraft Digital DATCOM® does not compute lift for an aircraft configuration which has no wing attached to the fuselage. Even if we add a fictitious wing in the hull, based on our experience on similar hybrid lifting fuselage, the lift estimated by CFD is 433.3% higher than that predicted by using the Aircraft Digital DATCOM® , since that produced by the HLH is also responsible for generating the induced drag. Therefore, the comparison of CD is also out of question. However, the same is not true for the case of pitching moment as the analytical relation derived by Munk-Multhopp (see Eqs. 7 and 8) contains the fuselage camber incidence angle (if) to represent the curvature effects of the fuselage (Multhopp 1942Multhopp H (1942) Aerodynamics of the fuselage. NACA-TM-1036.).

The wind tunnel test data has provided the aerodynamic as well as static stability characteristics of an HLH. This data can be used in flight performance estimation and system design of hybrid buoyant aerial vehicles. Except for the value of CD, the CFD results underpredicted the lift and pitching moment coefficients. Moreover, additional combinations of pitch and yaw angles can also be analyzed in the future to be added to the aerodynamic and stability database of such unconventional shaped hull.

CONCLUSION

During α sweep tests, it was observed that the CL was increased with some initial value. If we consider a hull as a wing, then a poor result of lift curve slope is obtained. The existing formula for airship FF has provided quite reasonable comparison with the experimental value. Overall, the CFD results were in good agreement with the wind tunnel ones. Moreover, the HLH is statically unstable in longitudinal and directional modes in the range of specified angles of attack and side slip angles.

ACKNOWLEDGEMENTS

The support of the Ministry of Science, Technology and Innovation (MOSTI), Malaysia, under the grant 06-01-08-SF0189, is gratefully acknowledged. The authors are thankful to the team of IIUM-LSWT for providing assistance during the experimental phase.

AUTHOR’S CONTRIBUTION

  • Haque AU conceived the idea, conducted the experiments and completed the literature review. Asrar W post process the raw data of wind tunnel, drafting of the article along with the CFD work for the purpose of comparison of the experimental data of the test article. Omar AA has done the data analysis and interpretation. Suleiman E discussed the findings of the research work and critical revision of the article.

REFERENCES

  • Aeroscraft (2013) Launching a new century of aviation innovation; [accessed 2013 Apr 2]. http://www.aerosml.com/imagegallery.html
    » http://www.aerosml.com/imagegallery.html
  • Ash LG (1972) Flight test and wind tunnel performance characteristics of the X-24A lifting body; [accessed 2016 Aug 26]. http://oai.dtic.mil/oai/oai?verb=getRecord&metadataPrefix=html&identifier=AD0901465
    » http://oai.dtic.mil/oai/oai?verb=getRecord&metadataPrefix=html&identifier=AD0901465
  • Becker JV (1958) Preliminary studies of manned satellites. Winged configurations. Proceedings of the NACA Conference on High-Speed Aerodynamics; Hampton, USA.
  • DeLaurier J, Schenck D (1979) Airship dynamic stability. Proceedings of the AIAA Lighter-Than-Air Systems Technology Conference; Palo Alto, USA.
  • Gumse B (1967) Low-speed wind-tunnel tests of a full-scale M2-F2 lifting body model. NASA TM X-1347.
  • Harrison J (2010) Lockheed’s LEMV consolation; [accessed 2012 May 8]. http://edgefighter.com/2010/06/17/lockheedslemvconsolation
    » http://edgefighter.com/2010/06/17/lockheedslemvconsolation
  • Hasim F, Rusyadi R, Surya WI, Asrar W, Omar AA, Mohamed Ali JS, Kafafy R (2008) The IIUM low speed wind tunnel. Proceedings of the 2nd Engineering Conference on Sustainable Engineering Infrastructures Development & Management (EnCon2008); Kuching, Malaysia.
  • Hoerner SF (1965) Fluid-dynamic drag: theoretical, experimental and statistical information. Bakerfield: published by the author.
  • Jackson EB, Christopher IC (1992) Preliminary subsonic aerodynamic model for simulation studies of the HL-20 lifting body. NASA TM-4302.
  • Millennium Air Ship Inc. (2012) The airship “Sky Freighter”; [accessed 2012 May 8]. http://www.millenniumairship.com/products.html
    » http://www.millenniumairship.com/products.html
  • Multhopp H (1942) Aerodynamics of the fuselage. NACA-TM-1036.
  • Munk MM (1924) The aerodynamic forces on airship hulls. NACATR-184.
  • Nicolai LM, Carichner GE (2013) Fundamentals of aircraft and airship design. Vol. 2: airship design and case studies. Reston: American Institute of Aeronautics and Astronautics.
  • Paulson JW, Shanks RE, Johnson JL (1960) Low-speed flight characteristics of reentry vehicles of the glide-landing type. NASA TMX-67563.
  • Pope A, Rae WH (1984) Low-speed wind tunnel testing. 2nd edition. New York: Wiley.
  • Prange HD, Schmidt-Nielsen K (1970) The metabolic cost of swimming in ducks. J Exp Biol 53:763-777.
  • Raymer DP (2012) Aircraft design: a conceptual approach. 5th edition. Reston: American Institute of Aeronautics and Astronautics.
  • Reed R, Darlene L (1977) Wingless flight: the lifting body story. NASA SP-4220; [accessed 2016 Aug 26]. http://history.nasa.gov/SP4220/contents.htm
    » http://history.nasa.gov/SP4220/contents.htm
  • Rizzo F (1924) A study of static stability of airships. NACA Technical Note 204.
  • Roskam J, Lan CTE (1997) Airplane aerodynamics and performance. Lawrence: DAR Corporation.
  • Saltzman EJ, Wang KC, Iliff KW (1999) Flight-determined subsonic lift and drag characteristics of seven lifting-body and wing-body reentry vehicle configurations with truncated bases. Proceedings of the 37th Aerospace Sciences Meeting and Exhibit; Reno, USA.
  • Santos Dumont A (1973) My airships; the story of my life. New York: Dover Publications. (Original work published in 1904).
  • Tanaka K, Hamaguchi Y, Maekawa S, Ishikawa T (2005) Structures and advanced composites for lighter-than-air aircraft. Proceedings of the 56th International Astronautical Congress; Fukuoka, Japan.
  • Tischler MB, Jex HR, Ringland RF (1981) Simulation of Heavy Lift Airship dynamics over large ranges of incidence and speed. Proceedings of the Lighter-Than-Air Systems Technology Conference; Annapolis, USA.
  • Trenkle K (2014) Hybrid Certification Criteria (HCC) for transport category hybrid airships document number: 1008D0122. Revision: C; [accessed 2014 Jun 4]. https://www.federalregister.gov
    » https://www.federalregister.gov
  • Tucker VA (1990) Measuring aerodynamic interference drag between a bird body and the mounting strut of a drag balance. J Exp Biol 154:439-461.
  • Ul Haque A, Asrar W, Omar AA, Sulaeman E, Mohamed Ali JS (2015a) Cambered profile of a California sea lion’s body. J Exp Biol 218:1270-1271. doi: 10.1242/jeb.117556
    » https://doi.org/10.1242/jeb.117556
  • Ul Haque A, Asrar W, Omar AA, Sulaeman E, Mohamed Ali JS (2015b) Preliminary aerodynamic and static stability analysis for hybrid buoyant aerial vehicles at low speeds using digital DATCOM. Can Aeronaut Space J 61(3):51-60. doi: 10.5589/q16-001
    » https://doi.org/10.5589/q16-001
  • Ul Haque A, Asrar W, Omar AA, Sulaeman E, Mohamed Ali JS (2016a) Hydrodynamic contour of Steller sea lion — an inspiration for future hybrid buoyant aircraft. Anim Rev 3(1):1-9. doi: 10.18488/journal.ar/2016.3.1/101.1.1.9
    » https://doi.org/10.18488/journal.ar/2016.3.1/101.1.1.9
  • Ul Haque A, Asrar W, Omar AA, Sulaeman E, Mohamed Ali JS (2016b) Assessment of engine’s power budget for hydrogen powered hybrid buoyant aircraft. Propulsion and Power Research 5(1):34-44. doi: 10.1016/j.jppr.2016.01.008
    » https://doi.org/10.1016/j.jppr.2016.01.008
  • Vogel S (1994) Life in moving fluids: the physical biology of flow. Princeton: Princeton University Press.
  • Vogel S (2013) Comparative biomechanics: life’s physical world. Princeton: Princeton University Press.
  • Wiriadidjaja S, Hasim F, Mansor S, Asrar W, Rafie ASM, Abdullah EJ (2012) Subsonic wind tunnels in Malaysia: a Review. Appl Mech Mater 225:566-571. doi: 10.4028/www.scientific.net/AMM.225.566
    » https://doi.org/10.4028/www.scientific.net/AMM.225.566
  • Wood RM (2003) The contributions of Vincent Justus Burnelli. Proceedings of the 41st AIAA Aerospace Sciences Meeting and Exhibit; Reno: USA.
  • Wood RM, Bauer SXS (2011) Flying wings/flying fuselages. Proceedings of the 39th Aerospace Sciences Meeting and Exhibit; Reno, USA.

Publication Dates

  • Publication in this collection
    Dec 2016

History

  • Received
    20 Mar 2016
  • Accepted
    14 July 2016
Departamento de Ciência e Tecnologia Aeroespacial Instituto de Aeronáutica e Espaço. Praça Marechal do Ar Eduardo Gomes, 50. Vila das Acácias, CEP: 12 228-901, tel (55) 12 99162 5609 - São José dos Campos - SP - Brazil
E-mail: submission.jatm@gmail.com