## Abstract

The paper presents results obtained by the author for wing-body interference. The lift-curve slopes of the wing-body combinations are considered. A 2D potential model for cross-flow around the fuselage and a discrete vortex method (DVM) are used. Flat wings of various forms and the circular and elliptical cross sections of the fuselage are considered. It was found that the value of the lift-curve slopes of the wing-body combinations may exceed the same value for an isolated wing. An experimental and theoretical data obtained by other authors earlier confirm this result. Investigations to optimize the wing-body combination were carried within the framework of the proposed model. It was revealed that the maximums of the lift-curve slopes for the optimal midwing configuration with elliptical cross-section body had a sufficiently large relative width (more than 30% of the span wing). The advantage of the wing-fuselage combination with a circular cross section over an isolated wing for wing aspect ratio greater than 6 can reach 7.5% at the relative diameter of fuselage equal to approximately 0.2.

### Keywords

- wing-fuselage combination
- lift-curve slope
- discrete vortex method
- 2D potential cross-flow model
- optimization

## 1. Introduction

An analysis of a lift-curve slope for wing-fuselage combinations currently plays an important role in studies of aerodynamics and the preliminary design of a modern aircraft.

Since the aircraft occurrence aircraft designers have been interested in the problems of the wing-body interference in aviation and missile technology. Initially, research is focused on the experimental study of specific wing-body combinations [1, 2, 3, 4, 5, 6]. First mathematical models of the wing-fuselage interference were offered later. The solution of the linearized problem of the ideal incompressible flow around arbitrary shape wings in the presence of the fuselage is a difficult task since it is necessary to solve the three-dimensional Laplace equation for the velocity potential which satisfies the boundary conditions on the surface of the wing-body combination and the boundary conditions at infinity. One of the few exact solutions was obtained by Golubinsky in the article [7]. The first theoretical calculations were based on the inversion of discrete vortices inside the cross-section body [8], on the solution of integral equations [9], on the application the thin body theory [10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21] or the strip method [14, 15, 22], and on the application of the velocity potential [23, 24, 25, 26] or the stream function [27] written in the Trefftz’s plane. The application of the velocity potential in the Trefftz’s plane [22, 23, 24, 25, 26, 28, 29, 30] gives the opportunity to get the distribution of lift along the wingspan. Let us pay attention to one important result that was first theoretically obtained by Multhopp and presented in a review by Ferrari [22]. It is important to note that this fact was experimentally confirmed by Jacobs and Ward [1]. This result shows that the value of the lift-curve slope of the wing-body combination at a certain relative value diameter of the fuselage * b* is the wingspan;

This chapter by no means covers all papers on the interference of the wing and fuselage. Author’s book [41] and paper [42] contain more detailed bibliography on the problems of the lift of the wing-body combination.

The main purpose of this paper is to give results of solving optimization problems for the values of the lift wing-body configurations and to demonstrate the conformity of computational author’s results with the known experimental and theoretical results of other authors.

## 2. The calculation method of the interference for wing-fuselage combination

The calculation method of the interference for the wing-fuselage combination [42] includes two methods: (1) a discrete vortex method (DVM) for the surface of the wing and (2) 2D potential model of the flow for cross-flow around fuselage [41].

The original three-dimensional problem (Figure 1) is divided into two parts. First part is the two-dimensional problem of the flow around the cross section of the fuselage (Figure 2), and the second part is a three-dimensional problem for the isolated wing. In the 2D problem, the flow around the cross section of the fuselage adds a pair of discrete point vortices. The added vortices are the consequence of lift on the wing. According to Zhukovsky’s theory about the lift of the wing, any lifting surface can be replaced by an equivalent Π-shaped vortex; free vortices at low angles of attack lie in the plane of the wing and extend to infinity. In our model, it is proposed each console part of the wing replaces one Π-shaped vortex lying in the plane of the wing. The Π-shaped vortex in the left-wing console is shown in Figure 1. The coordinate of the free vortex and its intensity can be found from the bond equation; after the lift-isolated wing by DVM will be defined. The inversion method (Figure 2) can be used to satisfy the boundary conditions of impermeability on the surface of the body cross section for the canonical body, and for the arbitrary two-dimensional cross section can use the panel method. An example solution for the potential flow around the elliptical cross section of the fuselage in the present of the pair vortices is shown in Figure 3.

In this formulation, the problem is reduced to solving the following system of algebraic linear equations:

where * L* is the number of control points (collocation points) equal to the number of attached vortices on the right-wing console,

n

*is the unit normal vector to the*

_{j}

*th control point on the surface of the wing, A*j

*is the matrix of the aerodynamic influence or the matrix of the induced velocities at the control points of the wing surface from all system of horseshoe vortices (left and right consoles for the isolated wing), and*

_{ij}

F

*is a column vector of the velocity induced in the*

_{j}

*th the control point on the wing surface by incoming flow and the flow from the cross section of the fuselage that includes either inversion of the vortices or sources and sinks providing satisfying conditions impermeability on the surface body from the free vortices left and right wing.*j

For small angles of attack of the wing-body combination

Right parts of the system of algebraic linear equations Eq. (1) can be represented also as

In Eqs. (2) and (3), * α* and wing deflection angle

*, respectively.*δ

For calculating the right parts of the system Eq. (1) for the problem with the fuselage of an arbitrary cross section of the body, the panel method that leads to the solution system of algebraic linear equations (4) is proposed:

where

Let us give the final formula for the components of the induced velocity, for example, for the case of the circular cross section of the fuselage in * j*th control point of the wing panel:

where * j*th control point wing panel induced velocity component along the

*-axis of the cylinder in cross-flow (see Figure 1). Coordinate inversion vortices are defined by Milne-Thomson’s theorem about the circle [43]. The coordinate*OZ

So the task of the wing-body interference is reduced to the solution Eq. (1) with the right part Eq. (6) or right-hand parts, obtained by solving the system (4) that provides the solution of the problem for the potential flow around an arbitrary contour of the panel method. The method of the successive iterations provides an agreement of the velocity field on the surface wing and the surface fuselage. Each iteration is reduced to the solution of systems of linear algebraic equations (1) with corrected right part Eq. (6). The zero iteration can select the solution for the isolated wing. For small angles of attack and wing deflection angle, the proposed model or the linear formulation allows to get the solution of two problems at once, which can be called αα-problem (fuselage and wing have the same angle of attack, angle of the wing deflection angle equal to zero) and δ0-problem (the fuselage has a zero angle of attack, and the wing has deflection angle δ not equal to zero). For this linear case, formulas for the coefficients of the normal forces of the wing and the body are of the form Eq. (7)

where values

## 3. Calculation results

The comparison of the calculation results obtained from the above theoretical model with calculations by the DVM for case αα-problem [34, 35, 36] is shown in Figures 4–9. The rectangular, triangular, and swept wings were considered. It may be noted is enough good agreement of calculated data.

The changing of coordinates of the aerodynamic center * x*/

_{AC}

*measured from the beginning of the mean aerodynamic chord is also shown in Figures 4, 6, and 8. The coefficient of the interference*c

where value

The comparison of the calculation results obtained from the above theoretical model with calculations by the numerical method of singularities for case δ0-problem [31] is presented in Figure 10.

The comparison of calculated data for the mathematical model described above and the calculated and experimental data of other researchers [45, 46, 47, 48, 49] is shown in Figures 11–13.

Figure 14 shows an influence of compressibility on the values of theoretical lift-curve slopes for case midwing monoplane combination with the rectangular and delta-shaped wing. Figure 15 also shows an influence of compressibility on values of theoretical lift-curve slopes for case high-wing monoplane combination with the rectangular and delta-shaped wing [33].

Of particular interest is the comparison of calculated and experimental data to prove that the lift-curve slope for the wing-body combination exceeds this value for an isolated wing. Figure 16 shows this comparison.

The area shown in color in Figure 16 indicates the advantage of the lift-curve slopes of the wing-body combinations over an isolated wing. Calculations and experiments show that with an increasing aspect ratio of the wing, this advantage will increase. This circumstance is important since the modern development of the aircraft industry tends to increase the aspect ratio of the wing. Another conclusion is that the maximum of the lift-curve slopes with a wing aspect ratio of 6 is achieved at relative fuselage diameter of approximately 0.2. Such a relative diameter of the fuselage allows the design of modern aircraft with a wide fuselage. Numerical studies have shown that with increasing aspect ratio of the wing and the ratio of the width to the height of the fuselage elliptical cross sections, the advantage of lift-curve slopes of the wing-body combinations over isolated wings becomes larger. The noted facts allow us to formulate and solve an optimization problem.

## 4. The formulation of the optimization problem

Note that in some theoretical and experimental papers devoted to the wing-body interference revealed a maximum dependence

We will use the formulation of the optimization problem as a nonlinear programming problem as follows:

where

where

## 5. Results of the optimization problem for lift-curve slope for midwing-body monoplane configuration

Figures 17 and 18 show results of the optimization problem for lift-curve slope for midwing-body monoplane configuration with circular cross-section fuselage vs. the aspect ratio of the rectangle wing. In Figure 16, the notation is used:

where

Figure 18 shows the results of the solution of the optimization problem for lift-curve slopes for midwing-body monoplane configuration with elliptical cross-section fuselage. Maximum values of the lift-curve slopes depend on the aspect ratio of the rectangular wing and the ratio of the axes of the ellipse. Figure 18 shows that the advantage of the wing-fuselage combination over an isolated wing is enhanced with increasing the aspect ratio of the rectangular wing and with increasing the ratio of the axes of the cross-section fuselage. The optimal ratio of the width of the body to the span of the wing can reach 30% and more!

Figure 19a shows the effect of the compressibility and the statistics for modern aircraft also (Figure 19b). Red color point shows the project of fifth-generation aircraft (project M-60, Russia). The feature of the project M-60 is a wide fuselage. As can be seen from Figure 19b, with the aspect ratio wing equal to 15, the optimal ratio of the width of the circular cross section to the wingspan can reach 20%!

## 6. Conclusions

The paper presents results obtained by the author for wing-body interference. The lift-curve slopes of the wing-body combinations are considered. A 2D potential model for cross-flow around the fuselage and the discrete vortex method for the wing were used. Flat wings of various forms and the circular and elliptical cross sections of the fuselage are considered. It was found that the value of the lift-curve slopes of the wing-body combinations may exceed the same value for an isolated wing. An experimental and theoretical data obtained by other authors earlier also confirms this result. Investigations to optimize the wing-body combination were carried within the framework of the proposed model. The proposed mathematical model for the solution optimization problem for the wing-body combination allows selecting the optimal geometric parameters for configuration to maximize the values of the lift-curve slopes of the wing-body combination.

It was revealed that the maximums of the lift-curve slopes for the optimal midwing configuration with elliptical cross-section body reach their values at sufficiently large relative width of the body (more than 30% of the span wing!). The advantage of the wing-fuselage combination with a circular cross section over an isolated wing at the wing aspect ratio greater than 6 can reach 7.5% at the relative diameter of fuselage equal to approximately 0.2. The advantage of the wing-fuselage combination with the elliptical cross section with the ratio of axes of the body equal to 2.5 over an isolated wing with aspect ratio equal to 12 is that it can reach 29% at relative width of fuselage equal approximately to 0.35!

## Acronyms and abbreviations

NACA | National Advisory Committee for Aeronautics |

AGARD | Advisory Group for Aerospace Research and Development |

ZFF | Zeitschrift für Flugwissenschaften |

TsAGI | Central Aerohydrodynamic Institute |

J. Aircraft | Journal of Aircraft |

R&M | Reports and Memoranda |

J. Aeron. Sc. | Journal of the Aeronautical Sciences |

ARC CP | Aeronautical Research Council Current Papers |

## References

- 1.
Jacobs EN, Ward KE. Interference of wing and fuselage from tests of 209 combinations in the N.A.C.A. Variable Density Tunnel. NACA Report; 1935. 540, p. 37 - 2.
Pitts WC, Nielsen JN, Kaattari GE. Lift and center pressure of wing-body-tail combinations at subsonic, transonic and supersonic speeds. NACA Report; 1957. 1307, p. 70 - 3.
Schneider W. Experimental investigation of wing-body interferences in the mach number range from 0.5 to 2.0. Transonic Aerodynamics. AGARD CP No. 35; Göttingen, Germany: Aerodynamische Versuchsanstalt Göttingen; 1968. pp. 20-1-20-23 - 4.
Kirby DA, Hepworth AG. Low-speed wind-tunnel tests on some slender airbus configurations. R&M; 1971. No. 3747, p. 55 - 5.
Körner H. Berechnung der potentialtheoretischen Strömung um Flügel-Rumpf-Kombinationen und Vergleich mit Messungen. ZFF; Vol. 20. 1972. pp. 351-368 - 6.
Korner SF, Borst HV. Fluid-Dynamic Lift. 2d ed. Published by Mrs. L.A; 1985. Hoerner, 482 p - 7.
Golubinsky AI. Exact solution of interference problem of wing with fuselage at subsonic flow. In: Proceedings TsAGI; 1961. 810, pp. 23-36 - 8.
Lebedev BF. Approximate method of calculating load distribution on wing and fuselage at subsonic speeds. In: Proceedings TsAGI; 1958. p. 719 - 9.
Dorodnitsyn AA. Influence of fuselage on distribution of loads on wing span. In: Proceedings TsAGI; 1944. p. 546 - 10.
Keldysh VV. Interference of flat swept low-aspect ratio wing and body. In: Proceedings TsAGI; 1959. 759, pp. 1-23 - 11.
Keldysh VV. Lift small aspect ratio wing with body. In: Scientists Note TsAGI; VI, 5; 1975. pp. 15-28 - 12.
Keldysh VV. Lift and longitudinal moment of low-aspect ratio wing with body of rotation located near it. In: Scientists Note TsAGI; VIII, 3; 1977. pp. 19-31 - 13.
Flax AH. Comment on “correlation of wing-body combination lift data”. Journal of Aircraft. 1974; 11 (5):393-394 - 14.
Flax AH. Integral relations in linearized theory of wing-body interference. Journal of the Aeronautical Sciences. 1953; 20 (7):483-490 - 15.
Flax AH, Lawrence HR. The aerodynamics of low-aspect-ratio wings and wing-body combinations. In: Third Anglo-American Aeron. Conf. Brighton; 1951; pp. 363-398, 398A-398J - 16.
Lawrence HR, Flax AH. Wing-body interference at subsonic and supersonic speeds—Survey and new development. Journal of the Aeronautical Sciences. 1954; 21 (5):289-324 - 17.
Nielsen JN. Missile Aerodynamics. New York-Toronto-London: McGraw-Hill Book Company, Inc.; 1960. 450 p - 18.
Spreiter JR. Aerodynamic properties of slender wing-body combination at subsonic, transonic and supersonic speeds. NACA TN; 1948. 1662 - 19.
Spreiter JR. Aerodynamic properties of the cruciform wing and body combinations at subsonic, transonic, and supersonic speeds. NACA TN; 1949. 1897 - 20.
Spreiter JR. Aerodynamic forces on slender plane and cruciform wing and body combinations, NACA Report. 1950. 962, pp. 271-287 - 21.
Spreiter JR, Sacks AH. A theoretical study of the aerodynamics of slender cruciform-wing arrangements and their wakes. NACA Report; 1956. 1296, pp. 81-106 - 22.
Ferrari C. Aerodynamic components of aircraft at high speeds vol. VII. In: Donovan AF, Lawrence HR, editors. High-Speed Aerodynamics and Jet Propulsion. Princeton, New Jersey: Princeton University Press; 1957. pp. 228-435 - 23.
Karafoli E. Aerodynamics of an Aircraft Wing. Moscow: AN USSR Publishing; 1956. 479 p. (In Russian) - 24.
Ryzhkov UA. Approximate method of calculating forces acting on body in supersonic gas flow and problem of aerodynamic interference of wing-fuselage. In: Proceedings TsAGI; 1959. 759, pp. 24-37 - 25.
Lennertz J. On the mutual reaction of wings and body. NACA TM, 400; 1927 - 26.
Lennertz J. Beitrag zur theoretischen Behandlung des gegenseitigen Einflusses von Tragfläche und Rumpf. Abhandlungen aus dem Aerodyn. In: Durand WF, editor. Aerodynamics Vol. IV. Berlin: Springer; 1934 - 27.
Nikolsky AA. About lifting properties and induced drag of wing-fuselage system. In: Proceedings TsAGI; 1981. 2122, pp. 94-95 - 28.
Weber J. Interference problems on wing-fuselage combinations. Part I: Lifting unswept wing attached to a cylindrical fuselage at zero incidences in midwing position. RAE Technical Report 69130 (ARC 31532); 1969 - 29.
Weber J. Second-order small-perturbation theory for finite wings in incompressible flow. ARC R&M 3759; 1972 - 30.
Weber J, Joyce MG. Interference problems on wing-fuselage combinations Part III symmetrical swept wing at zero incidence attached to a cylindrical fuselage. ARC CP; 1975. 1333, p. 84 - 31.
Vernigora VN. Numerical study of interference of thin trapezoidal wing and cylindrical fuselage at subsonic speeds. In: Proceedings TsAGI; 1983. 2176, pp. 26-43 - 32.
Kholyavka VI. Interference of aircraft parts. Kharkiv: Aviation. Inst.; 1967. p. 107 - 33.
Ganiev FI, Urkova ME. Calculation results of aerodynamic characteristics of wing-body combination moving with high subsonic speed. In: Proceedings Air Force Academy named after Zhukovsky; 1971. 1302, pp. 220-241 - 34.
Kolesnikov GА. Calculation of aerodynamic characteristics of wing-body combination at unsteady motion. In: Proceedings TsAGI; 1961. 810, pp. 78-96 - 35.
Kolesnikov GА. Order and results of calculation of aerodynamic characteristics of wing-body combinations at unsteady motion. In: Proceedings TsAGI; 1965. 954, pp. 3-28 - 36.
Belotserkovsky SM, Skripach BK. Aerodynamic Derivatives Wing and the Aircraft at Subsonic Speeds. Moscow: Nauka; 1975. 424 p. (In Russian) - 37.
Woodward FA. Analysis and design of wing-body combinations at subsonic and supersonic speeds. Journal of Aircraft. 1968; 5 (5):528-534 - 38.
Woodward FA. An improved method for the aerodynamic analysis of wing-body-tail configuration in subsonic and supersonic flow. NASA CR-2228, Part I and Part II; 1973. pp. 126 and 315 - 39.
Cvitanović I, Virag Z, Krizmanić S. Analysis of potential flow around wing-body configurations. In: 4th International Congress of Croatian Society of Mechanics (4th ICCSM); 18-20 September 2003; Bizovac, Croatia; 2003. pp. 1-10 - 40.
Davari AR, Soltani MR, Askari F, Pajuhande HR. Effects of wing geometry on wing-body-tail interference in subsonic flow. Scientia Iranica. 2011; 18 (3):407-415 - 41.
Frolov VA. Methods for Calculating Lift for Wing-Body Combination. Saarbrucken: LAP Lambert Academic Publishing; 2011. 141 p. (in Russian) - 42.
Frolov VA. Review results on wing-body. In: Proceedings of the International Conference on Mechanical, System and Control Engineering (ICMSE 2016); 18-20 May 2016; Moscow, Russia: MATEC Web of Conferences ICMIE; 2016. 75, 09006. DOI: 10.1051/matecconf/2016750900 - 43.
Milne-Thomson LM. Theoretical Aerodynamics. London: Constable and Company; 1958 - 44.
Mendenhall MR, Nielsen JN. Effect of symmetrical vortex shedding on the longitudinal aerodynamic characteristics of wing-body-tail combinations. NASA CR, 2473; 1975. p. 115 - 45.
Giesing JP, McDonnell-Douglas Rep. Douglas Aircraft Co., Report DAC, Vol. I, 67212. 1968 - 46.
Giesing JP, Kalman WP, Rodden WP. Subsonic unsteady aerodynamics for general configurations; Part I—Direct application of the nonplanar doublet-lattice method; Part II—Application of the doublet-lattice method and the method of images to lifting surface/body interference. USAF FDL-TR-71-5; 1971 - 47.
Martina AP. The interference effect of a body on the spanwise load distribution of two 45° sweptback wings of aspect ratio 8.02 from low-speed test NACA TN. 1956; p. 3730 - 48.
Gregoriou G. Zur gegenseitigen Beeinflussung eines endlich Langen Rumpfes und eines Flügels in Mitteldeckeranordnung bei Unterschallströmung. BMVg-FBWT-73-33; 1973 - 49.
Labrujere TE, Loeve W, Slooff JW. An approximate method for the calculation of pressure distribution on wing-body combinations at subcritical speeds. AGARD-CP, No. 071; 1970