Jtam.dvi JOURNAL OF THEORETICAL AND APPLIED MECHANICS 45, 2, pp. 425-436, Warsaw 2007 NUMERICAL INVESTIGATION OF AIRFOIL IN GROUND PROXIMITY Tomasz Abramowski Faculty of Maritime Technology, Technical University of Szczecin e-mail: tomasz.abramowski@ps.pl The paper presents numerical investigation of the ground effect which affects an air flow when a craft’s wing approaches the ground or sea surface. The effect has crucial significance for the wing-in-ground craft operation and decides about the purpose of building one at all. The two- dimensional viscous flowproblem is solved and the results of calculations are presented. Furthermore, some most important issues related to the basics of an ekranoplan application are very briefly discussed. On the basis of numerical calculations, an empirical formula has been proposed for quantitative assessment of the chord dominated ground effect. Key words: wing-in-ground effect, ekranoplan, numerical computations, fast transport 1. Introduction Wing-in-ground (WIG) boat is a promising mean of transport for the near future or at least one worth to be considered. Recently, people have been traveling at increasing speeds. Newmeans of transport are usually faster than the ones they replace. In marine transport, conventional displacement hulls reach their limits of performance and hardly can keep upwith the speed. It is very unlikely that any conventional boat would be able to operate at speeds higher than 100km/h with an acceptable fuel efficiency. Speeds of WIGs are much higher than ship speeds, and overall operational expenses are lower than those of planes.Higher speeds canbeachievedwithhydrofoil boats or so called surface effect ships. Definitely, a marine vehicle without any water contact would be the solution for a very low drag induced only by air. AWIG craft is a boat with wings that sails just above the water surface, floating ona cushionof high-pressureair regionbetween itswingand thewater surface. This cushion is created by the aerodynamic interaction between the 426 T. Abramowski wing and the surface, called the ground effect. A WIG boat is different from an aircraft because it cannot operate without the ground effect. There are many papers related toWIG’s construction or operation. Some of them are Afremov et al. (1996), Chawla (1988), who present wind tunnel results, also Kornev andMatveev (2003) who give details of their research on flight modelling problems. The problems with implementation of WIG crafts are presented by Taylor (2000). Experiments on aerodynamics of an airfoil subject to the ground effect have been recently carried out by Ahmeda and Sharmab (2004), who tested a symmetrical NACA 0015 in a low speed wind tunnel.The results of their investigationswere lift anddrag forces anddetailed characteristics of the flow. Another work is the paper of Rafiuddin (2005) who also presents experimental studies of the NACA 4415 airfoil in ground proximity. A numerical approach to the solution of the considered problem is presented by Park and Chun (1995). Fig. 1. AWIG craft recently introduced to operation byPacific Seaflight on Alaska 2. Overview of the ground effect When a wing approaches the ground two phenomena are actually involved in the increase of the lift force and reduction of the drag. The ground effect is a common name for both effects, which is sometimes confusing. These two phenomena are sometimes referred to as span dominated and chord domina- ted ground effect. The former results mainly in the reduction of the induced drag (D) and the latter in the increase of the lift (L). The designations span dominated and chord dominated are related to the the fact that the main parameter in the span dominated ground effect is the height-to-span ratio, whereas in the chord dominated ground effect it is the height-to-chord ratio. Numerical investigation of airfoil in ground proximity 427 As it was in thementioned above, the span dominated effect reduces drag. The drag of an aircraft can be split up into different contributions. The two main sources of drag are called the friction drag and induceddrag.The friction drag is caused by friction between the air and the surface of the craft, and is therefore dependent on its wetted area. The induced drag is sometimes called the lift induced drag because it is the drag due to generation of lift. When a wing generates positive lift, the static pressure on the lower side of thewing is higher than that on the upper side.At thewingtip, there is a complication: the high pressure area on the lower sidemeets the low pressure area on the upper side, therefore the air will flow from the lower side to the upper side around the wingtip. This is called the wingtip vortex. The energy that is stored in those vortices is lost and is experienced by the aircraft as drag. In free air, the vortices around the wing tips havemore space to develop than when they are bounded by the ground. The amount of induced drag is dependent on the spanwise lift distribution and the aspect ratio of thewing. A high aspect ratio wing has a lower induced drag than a low aspect ratio wing since its wingtip vortices areweaker. There is not enough space for the vortices to fully develop when awing is approaching the ground.The vortices are also pushed outward by theground, apparently the effective aspect ratio of thewingbecomeshigher than the geometric aspect ratio, Fig.2. Fig. 2. Span dominated ground effect Thechorddominatedground effect increases lift.Theair cushion is created byhighpressure thatbuildsupunder thewingwhen theground is approached. This is sometimes referred to as ram effect or ram pressure.When the ground distance becomes very small the air can even stagnate under the wing, giving the highest possible pressure, pressure coefficient unity. The chord dominated ground effect not always increases lift. It is possible, under certain conditions, 428 T. Abramowski that the lift reduces when an airfoil approaches the ground. This is the case when the bottom of the foil is convex and the angle of incidence is very low or negative, and in that case a venturi nozzle is created between the foil and the ground where high-speed low-pressure air sucks the airfoil down. This effect is used by race car designers to make it stick to the road at high speeds. To avoid this effect, the pressing side of the airfoil should be as flat as possible and the angle of attack should be positive. In the next part of the paper, only the chord dominated effect is studied in details. The Fluent Inc. numerical software was applied for calculations of viscous flow around a NACA airfoil. The results show significant influence of the flight height on the lift of the airfoil. 3. Governing equations and modelling assumptions For calculation of viscous flow around an airfoil, a numerical method based on solving equations describing the case under consideration, i.e. RANS equ- ations, is used.These equations have the following form for the incompressible, steady, two-dimensional flow ρ ( u ∂u ∂x +v ∂u ∂y ) ρ=F1− ∂P ∂x +µ (∂2u ∂x2 + ∂2u ∂y2 ) −ρ (∂u′u′ ∂x + ∂u′v′ ∂y ) (3.1) ρ ( u ∂v ∂x +v ∂v ∂y ) ρ=F2− ∂P ∂y +µ (∂2v ∂x2 + ∂2v ∂y2 ) −ρ (∂u′v′ ∂x + ∂v′v′ ∂y ) In the above equations, u, v are components of the mean velocity vector, P is the pressure, µ is the viscosity, u′, v′ are fluctuation parts of the velocity vector, F1,F2 are volumetric forces. Furthermore, themodelmust satisfy the continuity equation ∂u ∂x + ∂v ∂y =0 (3.2) TheSpalart-Allmaras turbulencemodel has been applied for themodelling of the Reynolds stresses. This is a single transport equation model solving directly for a modified turbulent viscosity, especially developed for aerospace applications. The overall of the applied technique is based on the finite volume method with segregated formulation and second-order upwind discretization schemes. The SIMPLE algorithm is employed for the coupling of velocity and pressure. The placement of the first grid pointwas established on the basis of a non- dimensional distance parameter y+, describing the local Reynolds number. Numerical investigation of airfoil in ground proximity 429 For the first estimation, y+ may be determined according to the theory of flat-plate flow, e.g. Schlichting (1968) y+ =0.172 (y L ) R0.9 n (3.3) where y is the distance from the wall, L – body length. The grid was con- structed in such a way that the parameter was maintained at y+ =100. The considered flow velocity (airfoil velocity) was assumed 30m/s, the Reynolds number of the flow was 2100000. This was close to the expected cruising velocity of a small WIG boat. The Mach number of flow was equal to 0.1, so then the compressibility effects could be neglected since they should be taken into account for a Mach number greater than 0.3. The flow was at first studied for the airfoil moving in unbounded air and then the results were compared against the calculations for the cases with the ground effect. The investigated heights of airfoil flights were introduced by relative ratios of the height-to-airfoil chord length (h/c). 4. Geometry, grid and computational domain The flow was studied for the NACA/Munk M15 airfoil with a flat shape of the pressing (bottom) side. This combined with the assumed angle of attack which was equal to 3◦, should result in the clearly visible chord dominated ground effect. The geometry of the airfoil is given in Fig.3. Fig. 3. Geometry of NACA/Munk M15 airfoil For discretization of the computational domain, an unstructured type of grid with quad elements was selected. The grid for the airfoil moving in free air is presented in Fig.4, while the grid used for calculations in the ground effect case is given in Fig.5. Inlet and outlet boundary conditions were specified on the outer sides of computational domain with necessary turbulence and flow parameters. The 430 T. Abramowski Fig. 4. Grid applied for airfoil moving in free air Fig. 5. Grid applied for airfoil moving close to surface boundary condition on the airfoil is a no-slip condition with zero relative spe- ed enforced. Just below the airfoil moving wall condition has been applied for better representation of ground modeling when the ground effect was in- vestigated, Fig.5. This should produce better representation of the reversing flow assumption. The velocity of the moving wall is the same as for the inlet condition, v=30m/s. Numerical investigation of airfoil in ground proximity 431 5. Results of calculations All applied physical models described above and grid topologies presented in Fig.5 and Fig.6 produced results which are presented in this Section. Com- puted pressure coefficients CP contours are given in Fig.6-Fig.9. The results of calculated distributions of CP values were compared to the potential flow solver, and the comparison is demonstrated in Fig.10. Due to restricted capa- bilities of available potential flow solver, the comparison was carried out only for the case of a flight in open air. The comparison of presented results with the potential flow solver in terms of the lift coefficient has given reasonable conformity. The lift coefficient for the unboundedair casewas found to be 0.6, while the one calculated with the help of potential solver 0.66. Fig. 6. Pressure coefficient contours CP . Airfoil moving in unbounded air The convergence history of the iteration process was investigated on the basis of scaled residuals of continuity, velocity and turbulence quantities. Usu- ally, it is alsoworth to study the quantity being the subject of interest, i.e. the lift force coefficient in the presented case. These values are shown in Fig.11, where the convergence history for the airfoil moving in unbounded air is pre- sented, while the history for the airfoil in the ground effect (c/h = 0.05) is given in Fig.12. The sufficient convergence was archived after about 900 iterations. On the basis of the function presented in Fig.12, an empirical formula has beenderived for practical andpreliminary calculations of the chord-dominated ground effect. It has the following form 432 T. Abramowski Fig. 7. Pressure coefficient contours CP . Airfoil moving in 0.25h/c Fig. 8. Pressure coefficient contours CP . Airfoil moving in 0.1h/c CLground =CL (h c )−0.11 (5.1) where CLground represents the wing lift coefficient in the ground effect mode, CL – lift coefficient in open air, h indicates height of the wing and c is the chord length. The above presented formula should be used with care, having in mind that it is based on numerical calculations. However, in terms of safety, we are on the safer side due to taking into account only one component of the overall Numerical investigation of airfoil in ground proximity 433 Fig. 9. Pressure coefficient contours CP . Airfoil moving in 0.05h/c Fig. 10. Calculated pressure coefficients CP – own calculations (right), compared with java-based, potential flow solver (left) available on-line at: http://www.aa.nps.navy.mil/∼jones/ ground effect, the chord dominated one. In a real WIG’s journey, the span dominated effect will occur aswell, which should produce some safetymargin. 6. Final remarks The increase of the numerically calculated lift coefficient in the groundmode was 40% when compared to the lift coefficient in unbounded air. The effect shows that wider application of WIG crafts has a great potential. Obviously, anyonewhoconsiders a small and fastmeanof transport in suchareas as rivers, 434 T. Abramowski Fig. 11. Convergence history of flow quantities for unbounded air case Fig. 12. Convergence history of flow quantities for ground effect case, h/c=0.05 lakes or gulfs should seriously have the WIG potential in mind. A company which is going to find a big, long range oceanic vehicle can also consider an ekranoplan. Cole (2002) presented a notable example of such a design idea. The results presented here can be applied at the preliminary design stage for the initial analysis of aWIG craft moving in the ground effect mode. The further work will deal with analysis of the second effect which affects total ground effect phenomena. This is the span dominated effect, and because it requires three dimensional calculations such analysis would take into account both effects and cumulative gains are expected to occur. For building a real WIG craft or at least amodel of it, one should also study the problemof flight stability and control. Numerical investigation of airfoil in ground proximity 435 Fig. 13. Lift coefficient for different values of h/c. Line of constant CL =0.6 is for airfoil in unbounded air References 1. AfremovA.Sh., ZhitmukA.P.,NikolaevE.P., SmolinaN.A., 1996,Hy- drodynamics of ekranoplanes, Proc. Int. Conf. 300th Anniversary of the Rus- sian Navy and Shipbuilding Nowadays, Saint-Petersburg, Russia 2. AgarwalR.K.,Deese J.E., 1984,Numerical solutions of theEuler equations for flow past an airfoil in ground effect,AIAA Paper, 1984-0051 3. AhmedaM.R., SharmabS.D., 2004,An investigationontheaerodynamicsof a symmetrical airfoil in groundeffect,Experimental Thermal andFluid Science, 29, 6, 633-647 4. ChawlaM.D., 1988,Wind tunnel investigationofwing-in-groundeffects,Col- lect. Techn. Pap. AIAA 6th Appl. Aerodyn. Conf., 147-153 5. Chen Y.S., Schweikhard W.G., 1985, Dynamic ground effects on a two- dimensional flat plate, Journal of Aircraft, 22, 638-640 6. ChubikovB.V.,PashinV.P.,TreshchevskyV.N.,MaskalikA.I., 1991, Ekranoplan - a High SpeedMarine Vehicle of NewType, Proc. of the Int. Conf. FAST’91, Trondheim 7. Cole W., 2002, The pelican: a big bird for the long haul, The Boening Fron- tiers, 01, available online at: http://www.boeing.com/news/frontiers/archive/2002/september/index.html 8. Grundy I.H., 1986,Airfoilsmoving close to a dynamicwater surface,J.Austr. Math. Soc., B27, 3, 327-345 9. Gryboś R., 1998,Podstawy mechaniki płynów, PWNWarszawa 10. IMO Resolution, 2002, Marine safety committee/circular 1054: interim stan- dards for wing-in-ground (WIG) craft 436 T. Abramowski 11. Kida T., Miyai Y., 1976, An alternative analytical method for ground effect airfoils,Aeronautical Quarterly, 292-303 12. Kornev N., Matveev K., 2003, Complex numerical modelling of dynamics and crashes of wing-in-ground vehicles, 41st Aerospace Sciences Meeting and Exhibit 13. MamadaH.,AndoS., 1973,Minimum induced drag of a hemicircular ground effect wing, Journal of Aircraft, 10, 660-663 14. Pacific Seaflight: http://www.pacificseaflight.com 15. Park I.R., Chun H.H., 1995, Numerical simulation of unsteady performance for 2D surface effect airfoils,Proc. of 95th Spring Meeting of Korea Committee for Ocean Resources and Eng., 10, 1 16. Plotkin A., Kennell C., 1981, Thickness induced lift on a thin airfoil in ground effect,AIAA Journal, 19, 1484-1486 17. Rafiuddin M.A., 2005, Aerodynamics of a cambered airfoil in ground effect, Int. Journal of Fluid Mechanics Research, 32, 2, 157-183 18. Schlichting H., 1968,Boundary Layer Theory, McGraw-Hill, New-York 19. Taylor G. K., 2000,Wise or otherwise. The dream or reality of commercial WIGvehicles,Proc. GEM2000 Int. Conf., St. PetersburgTechnicalUniversity, Russia 20. TheWIG page: http://www.se-technology.com 21. TuckE.O., 1980, Steady flowand static stability of airfoils in extreme ground effect, Journal of Eng. Mathematics, 15, 89-102 Numeryczne badanie efektu przypowierzchniowego profilu aerodynamicznego Streszczenie Artykuł prezentuje wyniki badań numerycznych jednego z dwóch zjawisk zmie- niających opływ skrzydła zbliżającego się do powierzchni gruntu lub morza. Efekt ten ma podstawowe znaczenie dla eksploatacji obiektu typu WIG, a przede wszyst- kim stanowi o sensie budowania takiej jednostki w ogóle. Uzyskano rozwiązanie dla przepływu dwuwymiarowego i zaprezentowanowyniki. Bardzo ogólnie zaprezentowa- no również podstawoweproblemy związane z zastosowaniemekranoplanu jako środka transportu. W oparciu o wyniki badań numerycznych zaproponowano wzór empi- ryczny do oceny ilościowej efektu przypowierzchniowego związanego ze zbliżaniem się profilu aerodynamicznego do powierzchni gruntu lubmorza. Zależność tąmożna wy- korzystać w praktyce, jeżeli tylko efekt związany z wirem wierzchołkowym skrzydła będzie miał pomijalny wpływ na cały efekt przypowierzchniowy. Manuscript received April 12, 2006; accepted for print January 3, 2007