Next Article in Journal
Incorporation of Pilot Factors into Risk Analysis of Civil Aviation Accidents from 2008 to 2020: A Data-Driven Bayesian Network Approach
Next Article in Special Issue
Rapid Blade Shape Optimization for Contra-Rotating Propellers for eVTOL Aircraft Considering the Aerodynamic Interference
Previous Article in Journal
Service Life Modelling of Single Lap Joint Subjected to Cyclic Bending Load
Previous Article in Special Issue
Investigation of Improvement Design on Aileron Surface Flow State of High Lift Configuration in BWB
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Fast Aerodynamic Model for Aircraft Multidisciplinary Design and Optimization Process

by
Frédéric Moëns
ONERA, Université Paris-Saclay, F92190 Meudon, France
Aerospace 2023, 10(1), 7; https://doi.org/10.3390/aerospace10010007
Submission received: 21 November 2022 / Revised: 13 December 2022 / Accepted: 16 December 2022 / Published: 22 December 2022
(This article belongs to the Special Issue Aerodynamics Design)

Abstract

:
A multidisciplinary design analysis and optimization process is developed at ONERA for the design of tube and wing and blended wing–body aircraft configurations. This process is composed of different disciplinary modules (geometry, propulsion, aerodynamics, structure, handling qualities and flight mission), and the overall process considers different fidelity levels for these modules at each step of the design process. This article describes the low-fidelity aerodynamic module used during the preliminary design optimization process. Analytical formulations retained for lift and drag components are presented in the first part. Then, the performances estimated by the aerodynamic module on some reference configurations are compared with both numerical and experimental data, showing a quite good agreement for both tube and wing and blended wing–body configurations not only for global performance but also for individual drag components.

1. Introduction

Among the different options for more efficient aircraft configurations, the blended wing–body (BWB) seems to offer a very promising reduction in emissions (CO2, NOx). However, such a configuration needs a complete, integrated design process as all the different disciplines (cabin arrangement, aerodynamics, structure, propulsion, flight dynamics, etc.) strongly interact altogether [1]. In order to explore, in depth, the potential benefits of such configurations, a multidisciplinary design analysis and optimization (MDAO) process dedicated to blended wing–body aircraft configurations has been developed since 2015 at ONERA [2,3,4,5]. The design workflow covers the range of tools from level 0 (L0), considering (semi-) empirical methods, over level 1 (L1), taking into account low-level physics-based methods, to level 2 and 3 (L2, L3), considering high-fidelity methods (Figure 1).
The modules implemented are easily exchangeable to adapt the workflow for specific configurations or a higher level of fidelity in specific domains if different methods are needed. This design tool evolves continually due to the enrichment of the disciplinary modules integrated or the addition of new disciplinary modules. The last version is composed of six main disciplinary modules: aircraft geometry, propulsion, aerodynamics, structure and weight, mission and performance and, finally, handling qualities. Additionally, some acoustic analysis can be performed offline of the multidisciplinary process for the selected configurations (Figure 2). Those modules are integrated in the NASA OpenMDAO [6] framework to constitute the multidisciplinary design analysis and optimization (MDAO) process.
When the main characteristics of the aircraft are not fixed, a large number of configurations are evaluated during the exploratory phase of the overall aircraft design (OAD) process with a large variety of architectures and planforms. As an example, Figure 3 presents some BWB planforms evaluated for building up a multi-fidelity surrogate model for the first pre-design optimization loop, showing a large variety of architectures to be considered (from [4]). A first constraint for the aerodynamic module is to be able to manage 100% of these different architectures.
Then, for each configuration, the mission module generates an aerodynamic database. As an order of magnitude, these databases consider about 5000 aerodynamic flow conditions to be computed for each planform, which makes the use of an aerodynamic disciplinary module with a very fast computing time and a reliable performance level prediction mandatory.
Therefore, among the different numerical methods available for aerodynamic modeling, the ones belonging to the “low-fidelity” family (L0 or L1) are often used for OAD phases. For the present application, a dedicated L0 aerodynamic module based on analytical formulations derived from the theory or from the data analysis of past or present aircraft is developed for a fast evaluation of the aerodynamic performance for standard “Tube and Wing” (T&W), “Flying Wing” (FW) or “Blended Wing Body” (BWB) configurations for typical subsonic cruise flight conditions. This module considers only wing planform data, and no surface grid is needed.
The different formulations retained are described in the first part of the article. In the second part, the results issued from the module are compared to the numerical or experimental results of different aircraft configurations.

2. Standard Atmosphere Model

To characterize the evolution of the ambient static flow conditions throughout the atmosphere, the aerodynamic module integrates an analytic formulation of the International Standard Atmosphere (ISA) [7]. ISA is a mathematical model that divides the atmosphere into several layers (troposphere, stratosphere, mesosphere, etc.) that gives the evolution of the ambient pressure, temperature and density with the altitude. The flight altitudes considered in the module belong to the troposphere and the stratosphere. Within each layer or sublayer, the temperature is assumed to have a linear evolution with respect to the altitude, the pressure is calculated by the hydrostatic balance from Equation (1), and the air density is calculated assuming the air as a perfect gas from Equation (2):
d P d z = ρ g
ρ = P R T
For the air considered as a perfect gas, there is R = 287.04 J kg−1 K−1, and the standard gravity constant is g = 9.80665 m/s2.
Knowing the evolution of the temperature with the altitude (Figure 4) and considering some thermodynamic assumptions for the different layers, the integration of Equation (1) between the two altitudes that bound the given layer leads to some analytic formulations, as presented in Table 1.
Finally, the air dynamic viscosity is computed using Sutherland’s relation:
μ = μ 0 ( T T 0 ) 3 2   ( T 0 + 110.4 T + 110.4 )
where μ0 = 1.711 10−5 kg/(ms) and T0 = 273.15 K.

3. Aerodynamic Module

The aerodynamic module is dedicated to T&W, FW or BWB configurations in subsonic, attached flow conditions. It does not consider flow conditions with separation, buffet, high-lift configurations, supersonic flight and unconventional aircraft configurations, such as multiplane or boxed wing.
In the first stage, a reference wing is built based on geometrical characteristics such as the airfoil relative thickness, the leading-edge sweep angles, the local chord, etc., and reference geometrical parameters are then derived (reference surface, mean aerodynamic chord). For the performance evaluation, this wing is considered “aerodynamically optimized” with an elliptical span loading in order to have an estimation of the local lift coefficients.
In the second stage, the influence of other elements, such as the fuselage, winglets or nacelles, are considered as drag increments. Note that the span loading of the wing considering these elements is not computed by the module.
The module is suited to compare different wing planforms within a preliminary design loop. By construction, identical results will be obtained for two configurations with the same wing planform and the same airfoil thickness spanwise evolution. The optimization of geometrical details, such as the airfoil shape, camber or twist, is not possible with the use of this module and has, therefore, to be considered in the next step of the design process using more advanced L1 or L2 numerical methods.
The following chapters present the geometrical inputs used and the analytical formulations considered for the estimation of the aerodynamic performance (lift and drag).

3.1. Geometrical Inputs

3.1.1. Wing

The wing planform is defined by the use of NAIRF airfoil sections of NWSEG segments (NAIRF = NWSEG+1). By construction, the first airfoil is located at the symmetry plane Y(1) = 0, and the last section Y(NAIRF) corresponds to the wing tip. For each wing segment, the leading-edge sweep angle, an index (XTra) for laminar computation and the airfoil technology factor used for wave drag calculations (“Korn factor”, see chapter 3.5) are provided. Additionally, for each airfoil, the relative thickness (t/c) and an estimation of the maximum lift (Cl max) are provided. Figure 5 gives an example of data for a wing defined by five airfoil sections (or four wing segments).
It should be noted that no detailed geometrical inputs for airfoils are used (shape, twist, camber).
Then, the wing span is divided into NDY subsections. For a given subsection, the different geometrical characteristics are obtained by a linear interpolation of the given parameter between two given wing sections of the segment. This allows the evolution of the leading-edge lines XLE(Y) and YLE(Y), the chord law C(Y) and the different sweep evolutions with wing span (φ0(Y), φ25(Y), φ50(Y)) to be obtained.
The reference area considered for the aerodynamic coefficient calculations is the geometrical wing area obtained by the sum of the different trapezoid segments:
S REF = S Wing = i = 0 i = N W S E G   [ Y ( i + 1 ) Y ( i ) ] × [ C ( i + 1 ) + C ( i ) ] 2
The wing geometrical aspect ratio λ corresponds to:
λ = ( 2   Y ( N A I R F ) ) 2 S W i n g = b 2 S W i n g
The wing taper ratio ε is:
ε = C ( N A I R F ) C   ( 1 )
Finally, the value of the aerodynamic mean chord and its location on the half-wing can be obtained:
A M C = 2 S R e f 0 b / 2 C 2 ( y ) d y
X A M C = 2 S R e f 0 b / 2 X L E ( y ) × C ( y ) d y
Y A M C = 2 S R e f 0 b 2 Y L E ( y ) × C ( y ) d y

3.1.2. Fuselage

A fuselage is modeled as a slender circular cylinder of length LFUS and a diameter DFUS. The aerodynamic module considers the fuselage to be set on the symmetry plane but with no consideration about its relative position to the wing. Therefore, similar performance will be evaluated by the module for the different configurations presented in Figure 6. The wing dihedral angle is not taken into account.

3.1.3. Nacelles

The effect of airframe installation on performance is considered for the drag contribution only (there are no propulsion effects considered). The evaluation is therefore similar to a “through-flow nacelle” consideration. Nacelle elements are considered as cylinders. Single flux (NAC1 elements only) or double flux nacelles (NAC1 and NAC2 elements) can be considered (Figure 7), but it is supposed that all the nacelles are identical in the present version of the module.

3.2. Formulation for Aircraft Lift

3.2.1. Estimation of Wing Maximum Lift

By construction, the aerodynamic module considers an ideal elliptic span loading, KZ(Y), for the reference wing, which is fully defined by the load at the symmetry plane KZ0:
K Z ( Y ) = K Z 0 1 ( Y b 2 ) 2
where:
K Z 0 = 2 S W I N G π ( b 2 ) × C L
For a given Y section, the local airfoil lift coefficient corresponds to:
C l   ( Y ) = K Z ( Y ) C ( Y )
Wing maximum lift is considered to be reached when Cl (Y) = Cl max (Y) at one Y section.

3.2.2. Lift Slope

The aircraft C L α slope (α in rad) considered is the Polhamus formulation [8], with the effects of fuselage taken into account [9]:
C L α = π λ [ 1.07   [ 1 + D F U S b ] 2 ] 1 + 1 + λ 2 ( 1 + tan φ 50 2 M 2 4 ) ( 1 D F U S b )
where φ50 is the mean wing sweep angle at mid-chord. For FW or BWB configurations, DFUS is set to zero.
Some modules of the MDAO process require the CL(α) curve as the input, which means that it is necessary to know the zero-lift incidence of the aircraft α0. The aerodynamic module is not able to determine this parameter, but it can be estimated based on preliminary validation exercises on similar configurations. Otherwise, a value of α0 = 0° can be used.

3.3. Formulations for Aircraft Drag Evaluation

Among the different existing drag formulations, the one retained for the aerodynamic module is derived from the one presented by Gur [10]:
C D T o t a l = C D I n d u c e d + C D F r i c t i o n + C D A d d + C D W a v e + C D P a r a s i t i c
The following sections describe the formulations used for these different terms.

3.3.1. Lift-Induced Drag

Only the wing contribution to lift-induced drag is considered with the influence of the fuselage or winglet. The contribution of the tail surfaces to the total lift-induced drag is not considered in the module. The standard formulation is considered for the lift-induced drag coefficient of the reference wing:
C D I n d u c e d = C L 2 π   λ   O s w
where λ is the wing geometric aspect ratio defined in Equation (5), and Osw is the Oswald factor that characterizes the deviation from the ideal elliptic span loading. For a given wing planform, several methods are available for the determination of the Oswald factor (see [11] for instance). Most of them consider an Oswald factor such as:
  O s w = 1 1 + δ
where δ is a parameter based on wing geometrical characteristics. The one proposed in the aerodynamic module is based on the formulation used by Anderson [11,12]:
δ A n d e r s o n = [ 0.0015 + 0.016 ( ε 0.4 ) 2 ] × [ ( λ 1 M 2 ) 4.5 ]
leading to:
O s w A n d e r s o n = 1 1 + δ A n d e r s o n
However, this formulation is valid for unswept wings only, which is obviously not suitable for BWB configurations. To take the wing sweep angle into account, Hörner [12] proposes a simple correction:
O s w φ = O s w ( φ = 0 ) × cos φ
Some validations carried out on several CFD results analyzed by a far-field drag decomposition tool have shown that the lift-induced drag component computed for swept wings was surrounded by the evaluations made considering OswAnderson or Oswφ and that a mean value between these two formulations leads to an excellent agreement (Figure 8). Therefore, the formulation retained in the code is simply:
O s w W i n g = 1 + cos φ 2   O s w A n d e r s o n
When a fuselage or winglets are considered, the Oswald factor is modified, as described by Nita [11]:
O s w W i n g + F u s e l a g e + W L T = O s w W i n g × K F u s × K W L T
The correction factor for the fuselage is:
K F u s = ( 1 2 ( D F U S b ) 2 )
DFus is the fuselage diameter, and b is the total wing span.
The correction factor for the winglets is:
K W L T = 1 C o e f W L T ( 1 + 2   H W L T b ) 2
where CoefWLT is a correction factor derived from works from Bourdin [13] or Delavenne [14] to take into account the effect of the winglet cant angle on lift-induced drag:
C o e f W L T = 1.0 + 4.10 4 δ W L T + 1.10 5 δ W L T 2 3.10 8 δ W L T 3 5.10 10 δ W L T 4
HWLT is the winglet height, b is the total reference wing span without winglets, and δWLT is the winglet cant angle (in °) (see definitions in Figure 9).

3.3.2. Friction/Form Drag

According to the methodology described in [9,10], the skin friction and pressure drag of the different components are calculated using the following relation:
C D F r i c t i o n = C F · F F · S W E T S R E F
where CF is a flat plate skin friction coefficient, FF is the form factor of the component, and SWET is the wetted area of this component.
  • Friction coefficient CF
For a flat plate of length L, the friction coefficient is obtained considering the definition of the boundary layer momentum thickness. At the end of the plate, the integrated friction coefficient is:
C F = 2 Θ L C ( Y )
If the flow is laminar, the use of the Blasius relationships leads to:
C F L a m = 1.328 R e L
For a turbulent flow, the friction coefficient is obtained by the compressible Schlichting relation [9]:
C F T u r b = 0.455 ( l o g 10 R e L ) 2.58 ( 1 + 0.144 M 2 ) 0.65
For wings considering natural laminar flow (NLF) or hybrid laminar flow (HLF) technologies, it is possible to have an estimation of the maximum laminar extent X Lam max for the wing sections as a function of the local Reynolds number and the φ25 sweep angle (Figure 10, from [15]). This maximum laminar extent occurs for a design lift coefficient ( C L Adapt ), and the change in transition location due to the modification of the pressure gradients over the airfoil when the angle of attack changes has to be taken into account.
Though the module does not consider any airfoil shape, it is therefore necessary to consider a “suction” and a “pressure” side for the wing surfaces in order to take these changes in transition location into account. The module uses the simplified model presented in Figure 11.
Once the portion of laminar flow (XLam) is estimated for the considered wing section, the global friction coefficient is obtained via the computation of the boundary layer momentum thickness at the trailing-edge through the use of a “fictitious turbulent length” (LFict in Figure 12).
Firstly, the momentum thickness at the transition point (ΘTr) is computed using the Blasius relationships.
Θ T r = ( X L a m × C ( Y ) ) 0.664 R e ( L = X L a m × C ( Y ) )
Then, using the Michel relationships [16], the length LFict of a fully turbulent boundary layer with the same ΘTr momentum thickness is estimated as:
L F i c t = [ 1 0.02208   Θ T r ( R e ( L = X L a m × C ( Y ) X L a m × C ( Y ) ) 1 / 6 ] 1.2
Finally, ΘTE is computed as the momentum thickness of a turbulent flow that develops over a length of LTur:
L T u r = L F i c t + ( 1 X L a m ) × C ( Y )
Θ T E = 0.02208 L T u r ( R e L T u r ) 1 / 6
The friction coefficient for transitional flow on one surface is finally obtained by introducing ΘTE from Equation (32) into Equation (26).
  • Form Factors FF
There are several models in the literature available for form factors (see for instance [10] for a comparison of different formulations for FFWING or for body of revolutions). Table 2 presents the ones retained in the aerodynamic module and the corresponding SWET.

3.3.3. Interaction Effects for Nacelles

In addition to the simple geometrical effect on the friction drag, some interaction effects of the nacelles with a surface (wing or fuselage) are considered through the distance parameter ZNAC (Figure 13). The positive values of ZNAC are for standard nacelle arrangements. The negative values of ZNAC are for buried nacelles. On the basis of these geometrical characteristics, the interaction coefficient QN is defined according to the following statement given in [9]: “For a nacelle or external store mounted directly on the fuselage or wing, the interference factor QN is about 1.5. If the nacelle or store is mounted less than about one diameter away, the QN factor is about 1.3. If it is mounted much beyond one diameter, the QN factor approaches 1.0.”
For the positive values of ZNAC, a linear evolution of QN is considered in the module:
If   Z NAC > 0 ;   Q N = m a x ( 1 , 1.5 0.25   Z N A C D F a n )
For the negative values of ZNAC, the interference factor is considered at its maximum value (i.e., 1.5), but we take into account the modification of the wetted surface of the nacelle by imposing a minimum value of 1.0 to the factor. It leads to the following relationship:
If   Z NAC <   0 ;   Q N = m a x ( 1.0 ,   1.5 × ( 1 a c o s ( 1 + 2 Z N A C D F a n ) π ) )
Once this coefficient is evaluated, the global friction drag coefficient for the NENG nacelles is obtained by the following relation:
C D E N G = N E N G × ( Q N C D N A C 1 + C D N A C 2 )
where the individual friction drag coefficients for NAC1 and NAC2 are computed according to Equation (25) with the form factor and wetted areas given in Table 2.

3.4. Additional Drag Due to Lift

An additional profile of drag due to lift is considered for the wing. For each wing segment, the additive drag is estimated according to the relation from [18]:
C D A D D = 0.75   C D A D D R e f ( C L C L 0 C L M a x C L 0 ) 2   1 ( M c o s φ 25 ) 2   S S e g m e n t S R E F
where:
C D A D D R e f = [ 0.010   C L M a x 0.0046 ( 1 + 2.75 ( t c ) + 100 ( t c ) 4 ) ] c o s 3 φ 25
C L Max is the maximum lift coefficient of the wing estimated from Equation (12), and CL0 is the minimum drag lift coefficient provided in the data file.

3.5. Wave Drag

For aircraft missions at transonic flight conditions, the drag increase due to compressibility effects has to be considered. In the module, the estimation of this drag component is based on the Korn equation [10]. In the original work, Korn gives an estimation of the divergence Mach number for an airfoil as:
M D D = K A 1 10 C l ( t c )
The KA coefficient, referred to as the “Korn factor”, is an airfoil technology coefficient, depending on the nature of the airfoil. It is proposed to use KA = 0.95 for “modern” supercritical airfoils and KA = 0.87 for “conventional” airfoils. For wing segments corresponding to a transition area with the fuselage on a BWB configuration, KA = 0.90 can be used. However, Equation (38) is valid for 2D airfoils only. In order to have an estimation of MDD for a wing section, the simple swept wing theory, considering a normal to the leading-edge flow assumption, can be used:
M 2 D = M 3 D cos φ
C l = C L 3 D   1 c o s 2 φ
( t c ) 2 D = ( t c ) 3 D 1 cos φ
Introducing Equations (39)–(41) into Equation (38), the following relation for the 3D swept wing section is obtained:
M D D = K A cos φ 1 10   C l c o s 3 φ 1 c o s 2 φ ( t c )
Then, the local critical Mach number MCr is estimated according to:
M C r = M D D ( 0.1 80 ) 1 3 = M D D 0.108
Finally, for a given wing subsection, the wave drag contribution is obtained:
I f   M > M C r   ;    C D   W A V E = 20   ( M M C r ) 4   S e l e m e n t S R E F
The total wave drag of the wing can be found by adding the contributions of all the NDY wing subsections.

3.6. Parasitic Drag

Finally, an additional parasitic drag ( C D PARASITIC ) (due to protuberances, antenna, probes, paint, etc.) is considered. A general formulation used in the pre-design phase is to consider this contribution as a percentage of the total friction drag:
C D P A R A S I T I C = X P A R A × C D F R I C T I O N
For an aircraft performance estimation, XPARA = 0.025 is generally used. However, for a validation purpose with comparisons with the CFD results or experimental data in a wind tunnel, XPARA = 0.0 has to be used.

4. Module Validation

The following chapters present the validation of the aerodynamic model on standard tube and wing or blended wing–body configurations. The results obtained by the module are compared to the CFD results or experimental data in terms of overall performance and drag components when possible.
Note that all the drag coefficients presented are expressed in drag counts.

4.1. Standard Tube and Wing Configurations

4.1.1. ONERA NOVA Configuration

The first configuration considered is the Nextgen ONERA Versatile Aircraft (NOVA) [19] (Figure 14), which is designed with a downward winglet element and different engine configurations.
The aerodynamic flow conditions for this test case are a Mach number of M = 0.82 at an altitude of Z = 11,300 m. For this configuration, drag decompositions using the ffd72 tools [20] of a different configuration are available, which makes it possible to compare the lift-induced drag for the clean wing and the wing equipped with the winglet with a negative cant angle δWLT = −18°.
The comparison with the prediction from the aerodynamic module (Figure 15) shows a quite good prediction of this drag component, even though only four wing segments are used for the wing definition.
NOVA is a platform used to investigate the effect of the integration of ultra-high bypass ratio (UHBR) engines. Different engine arrangements are available (Figure 16) that make possible a first validation of the different nacelle interference ratios on drag.
Table 3 compares the drag increase for different engine installations for the NOVA aircraft between the CFD results and the estimation from the module. A quite good agreement can be found for the different cases proving that the formulation of the interference factor for the nacelle is realistic.

4.1.2. NASA Common Research Model (CRM)

The second configuration considered for the validation exercise is the NASA CRM model (Figure 17), for which both experimental and numerical data are available [21,22,23].
For the evaluation of the performance, the module considered only six wing sections (Figure 18).
The flow conditions considered are the ones used for the different AIAA Drag Prediction Workshop (DPW) exercises: model scale, 1/10; Mach number = 0.850; and Reynolds number = 5.36 × 106.
Note that for the CRM configuration, the reference area used for the force coefficients is not the geometrical wing one, as considered by the module (SWing = 411.806 m2), but a modified trapezoidal one (“Wimpress” area = 383.69 m2). For the present comparisons, the force coefficients considered for the CRM database (from CFD or experiments) are therefore corrected in order to deal with this difference in SREF with a factor of 0.93172.
The first validation exercise considered the effect of the fuselage on the CL(α) curve. Figure 19 compares the CL(α) curves from the aerodynamic module, with or without a fuselage taken into account, with the CFD results from ONERA obtained within the DPW framework. The module considers a zero-lift incidence of α0 = −1.55° to improve the comparison. A quite good agreement can be observed for the estimation of the lift slope of the wing–fuselage arrangement by the module with the CFD data.
Regarding the overall drag estimation, Figure 20 compares the CL(CD) curves computed by the module with the CFD results from ONERA [23]. The drag coefficients from the CFD results were obtained by the ffd72 far-field analysis tool [20] that allows the elimination of the artificial spurious drag from the numerical solution and gives the drag breakdown of the physical components (lift-induced, viscous and wave). It can be seen that the agreement is quite good (ΔCD = 1 d.c. maximum) between the module evaluation and the CFD results up to lift coefficients of around 0.50 (using SREF = SWING) as well as for the wing–fuselage configuration and for the wing–fuselage–horizontal tail planes case.
When considering the different drag components, namely, lift-induced, viscous and wave drag, a very good agreement can also be found (Figure 21). Some discrepancies can be observed for CL values higher than 0.50 for viscous and wave drag, but in these conditions, a separation can be found on the wing by CFD, which is not considered in the aerodynamic module formulations.
Finally, some experimental data are provided on the CRM website for a wing equipped with or without through-flow nacelles. Due to their relative position with the wing surface, nacelles are considered as out of interaction (QN = 1). It can be seen that the drag increment found experimentally is correctly predicted by the aerodynamic module (Figure 22).

4.1.3. CRM-NLF Configuration

The next validation exercise considers the CRM model with a modified wing with a natural laminar flow on the upper surface. This CRM-NLF wing has been designed for tests in the NASA NTF wind tunnel [24,25,26], and some experimental results are provided on the NASA CRM website. The data used for the CRM reference plane are modified for the inboard wing section in order to take the wing planform change (Figure 23) into account. Similar to the CRM reference case, the experimental data are corrected in order to take the change in the reference wing area into account.
The flow conditions used for the validation are similar to the NTF ones: a Mach number of M = 0.8565 and a Reynolds number of 15×106. The NLF capabilities of the aerodynamic module are considered for this validation exercise, but wing lower surfaces are considered turbulent, as during the wind tunnel tests. The estimated maximum location of the transition point on different wing sections is presented in Figure 24 at CL = 0.50 for the wind tunnel flow conditions. It can be seen that an extended natural laminar flow is possible (around 50%) for a large portion of the wing, except for some portions of the inboard wing, which are more limited in terms of NLF capabilities.
Then, these values are considered as the X Lam max parameter, as shown in Figure 11, the location of the transition point XLam is computed for each CL considering C L Adapt = 0.50. The comparison of the transition line estimated by the module at CL = 0.425, with infrared pictures available on the CRM website, is presented in Figure 25. It can be seen that, though not exact, the laminar flow extent estimated by the module is quite realistic, especially on the outboard wing.
Moreover, of most importance is the evaluation of the effect of this laminar flow extent on the overall performance. Figure 26 compares the CL(α) and the CL(CD) curves from the aerodynamic module with the experimental data. A zero-lift angle of attack of α0 = −0.90° is considered to improve the comparison for the CL(α) curve. For comparison, the performances estimated by the module in fully turbulent mode are presented in dashed lines. It can be seen that the agreement on the drag evaluation is quite good for NLF flow conditions around the different CL values available (between 1.5 d.c. and 3 d.c.).

4.1.4. NLF Regional Aircraft Configuration

The CRM-NLF configuration considered laminar flow on the upper wing only. In order to evaluate the performance prediction capabilities of the aerodynamic model for a wing in which laminar flow extents on both surfaces, we considered the AG2-NLF wing design by ONERA in the framework of the Clean Sky 2 AIRGREEN2 program. The reference aircraft considered is a 90-pax turboprop configuration with wing airfoils redesigned by ONERA at cruise conditions for natural laminar flow capabilities [27,28,29]. This configuration is referred to as AG2-NLF. Figure 27 presents the complete configuration, which considers some winglets, an under-carriage fairing and a wing–body junction Karman that are not modeled in the aerodynamic module. In order to take into account the effect of the Karman on the viscous drag, the first wing section is considered a thick airfoil (t/c = 40%) with a leading-edge sweep of 30° at the symmetry plane.
The flow conditions are a Mach number of M = 0.52 at an altitude of Z = 6100 m (20,000 ft). The NLF capabilities of the module are evaluated by comparison with 3D RANS results obtained by the ONERA elsA software with integrated transition prediction capabilities [30]. The results are presented in Figure 28. The zero-lift angles of attack considered for the CL(α) curves are α0 = −1.7° for turbulent conditions and α0 = −1.96° for the NLF aircraft. It can be seen that the agreement on the drag evaluation is quite good for both the turbulent and NLF flow conditions at around CL_Adapt = 0.50, with a difference of less than 1 point in LoD.

4.2. BWB Configurations

As the module is used for an MDAO design process of BWB configurations, it is important to check if the formulations used are valid for such configurations with high values of leading-edge sweep angle.

4.2.1. AVECA

The first BWB configuration considered is one of the optimized planforms of the AVECA project [31,32]. AVECA is a long-range BWB configuration designed by Airbus. Several wing planforms, considering some changes in the geometrical characteristics or different volume constraints for the passenger cabin or cargo, have been optimized by ONERA using the adjoint method. One of the optimized configurations is selected for the validation purpose considering only four wing segments or five wing sections (Figure 29), which is very crude. The flow conditions are a Mach number of M = 0.850 at an altitude of Z = 11,000 m. The results of the aerodynamic module are compared to the CFD results post-processed with the ffd72 tool in order to obtain the drag breakdown between the different drag components. Note that the engines were not considered for both CFD and the module.
Figure 30 compares the CL(α), with a zero-lift angle of attack of α0 = +0.175° to improve the comparison, and the CL(CD) curves from the aerodynamic module (continuous lines) with the CFD results (symbols). It can be seen that the agreement on the lift slope and on the total drag evaluation is quite good. Similarly, a very good correlation can be found on the different drag components (lift-induced, viscous and wave) presented in Figure 31, with a small underestimation of the viscous drag by the module.

4.2.2. NACOR-SMILE Configuration

The second BWB configuration considered for the validation exercise is the optimized SMILE planform from the CS2 NACOR program. This configuration is designed for a short–medium range (SMR) mission, similar to the A320 aircraft. Therefore, compared to the AVECA configuration, its span is much smaller (36 m for SMILE, 80 m for AVECA). It should be noted that the aerodynamic module is used within the OAD-MDAO definition phase, but the final optimization details are carried out using the CFD methods [33].
Figure 32 compares a CAD rendering of the final optimized shape with the data used by the aerodynamic module. Only six wing sections are considered to model the BWB planform. Winglet and nacelles are also considered for the performance evaluations by the module.
Figure 33 compares the CL(α), with a zero-lift angle of attack of α0 = −0.225° to improve the comparison, and the CL(CD) curves from the module with the CFD results. It can be seen that the agreement on the lift slope and on the drag evaluation is quite good. Similarly, a very good correlation can be found on the different drag components (lift-induced, viscous and wave) presented in Figure 34.

5. Computing Performances

The use of analytical formula allows quite fast computing times for this module, which makes it quite useful for its use within a multi-disciplinary optimization/design process where a large number of configurations are considered. As an example, for a complete aircraft configuration, as the CRM case presented, the data file considers about 90 lines only (including comments), and the computing time to obtain one complete CL(CD) curve (~120 values) is in the order of 1 sec on a working station (DELL Precision 3630).

6. Conclusions

In the framework of a multi-disciplinary design analysis and optimization process for BWB configurations, an analytic aerodynamic module was developed in order to evaluate the aerodynamic lift and drag components of different wing planforms and aircraft architectures with a fast restitution time.
The module considered analytic formulations derived from the theory, literature or from statistical data. Some validation exercises showed that it could be used to estimate the aerodynamic performances of T&W and BWB configurations within an MDAO pre-design process for subsonic cruise flight conditions with a quite satisfactory level of accuracy. The lift slope, the total and individual drag components (lift-induced, viscous and wave) estimated by the module were in quite good agreement with the reference data issued from the numerical and experimental data. In addition, the effects of different elements on the performance, such as fuselage, winglets, nacelles and tail surfaces, were well captured by the module. Finally, the performance of the wings using laminar flow technologies (NLF or HLF) could be estimated with a quite good level of accuracy for use in a pre-design phase.
However, it is important to note that this module cannot be used for detailed optimization (airfoil shape, wing twist or camber, nacelle positions). For instance, the module will estimate the same performances for two wings with similar planform and airfoil thickness evolution with span. These fine optimization steps should be considered in the next stage of the aircraft design process, using more advanced numerical methods once the general architecture is defined in the preliminary phase.

Funding

The work presented in this paper was funded by the internal research project CICAV funded by ONERA from 2015 to 2019 and by the EU NACOR Project, which has received funding from the Clean Sky 2 Joint Undertaking under the European Union’s Horizon 2020 research and innovation program under grant agreement N°. CS2-AIR-GAM-2018-2019-01. Part of the results used in this paper were carried out in the framework of AIRGREEN2 Project, which received funding from the Clean Sky 2 Joint Undertaking under the European Union’s Horizon 2020 research and innovation program, grant agreement N°. 807089—REG GAM 2018—H2020-IBA-CS2-GAMS-2017.

Institutional Review Board Statement

Not Applicable.

Informed Consent Statement

Not Applicable.

Data Availability Statement

Experimental and reference CFD data related to the CRM configurations are available at CRM (https://commonresearchmodel.larc.nasa.gov/ (accessed on 21 November 2022)) or at the Drag Prediction Workshop (https://aiaa-dpw.larc.nasa.gov/ (accessed on 21 November 2022)) websites.

Conflicts of Interest

The author declares no conflict of interest.

Nomenclature

bWing span (m)
cAirfoil chord (m)
d.c.Drag count (CD × 104)
tAirfoil thickness
zAltitude (m or km)
PPressure (Pa)
TTemperature (K)
VVelocity (m/s)
XCoordinate in the streamwise direction (m)
YCoordinate in the spanwise direction (m)
LoDLift over drag ratio
MMach number
MDDDivergence Mach number
MCrCritical Mach number
NNumber of components
ReReynolds number
FFForm factor
QNInterference factor for nacelles
OswOswald factor
SWINGWing area (m2)
SREFReference area (m2)
ClLocal (airfoil) lift coefficient
CLLift coefficient
CDDrag coefficient
KAAirfoil Korn factor
XTraIndex for turbulent or laminar flow computations
XLamLaminar flow extent (%)
NDYNumber of wing section subdivisions
AMCAerodynamic mean chord (m)
Greek Symbols
αAngle of attack (°)
δParameter for Oswald factor
δWLTWinglet cant angle (°)
λWing aspect ratio
ρAir density (kg/m3)
μDynamic viscosity (kg/(ms))
φSweep angle (°)
εWing taper ratio
Subscripts
LELeading-edge
TETrailing-edge
FUSFuselage
WLTWinglet
ENGEngines
NAC1Nacelle or Fan
NAC2Turbine
TAILTail surfaces

References

  1. Liebeck, R.H. Design of the Blended Wing Body Subsonic Transport. AIAA J. Aircr. 2004, 41, 10–25. [Google Scholar] [CrossRef] [Green Version]
  2. Gauvrit-Ledogar, J.; Defoort, S.; Tremolet, A.; Morel, F. Multidisciplinary Overall Aircraft Design Process Dedicated to Blended Wing Body Configurations. In Proceedings of the 2018 Aviation Technology, Integration, and Operations Conference, Atlanta, GA, USA, 25–29 June 2018; p. 3025. [Google Scholar]
  3. Tremolet, A.; Gauvrit-Ledogar, J.; Brevault, L.; Defoort, S.; Morel, F. Multidisciplinary Overall Aircraft Design and Optimization of Blended Wing Body Configurations. In Proceedings of the 2019 European Conference for Aeronautics and Space Sciences, Madrid, Spain, 1–4 July 2019. [Google Scholar]
  4. Iwanizki, M.; Wöhler, S.; Fröhler, B.; Zill, T.; Méheut, M.; Defoort, S.; Carini, M.; Gauvrit-Ledogar, J.; Liaboeuf, R.; Tremolet, A.; et al. Conceptual Design Studies of Unconventional Configurations. In Proceedings of the 3AF Aerospace Europe Conference 2020, Bordeaux, France, 25–28 February 2020. [Google Scholar]
  5. Gauvrit-Ledogar, J.; Tremolet, A.; Moens, F.; Meheut, M.; Defoort, S.; Liaboeuf, R.; Morel, F.; Paluch, B. Multidisciplinary Design Analysis and Optimization Process dedicated to Blended Wing Body Configurations. In Proceedings of the 33rd Congress of ICAS, Stockholm, Sweden, 4–9 September 2022. [Google Scholar]
  6. Gray, J.; Moore, K.; Naylor, B. OpenMDAO: An Open Source Framework for Multidisciplinary Analysis and Optimization. In Proceedings of the 13th AIAA/ISSMO Multidisciplinary Analysis Optimization Conference, Fort Worth, TX, USA, 13–15 September 2010; p. 9101. [Google Scholar]
  7. US Standard Atmosphere 1976, NASA-TM-X 74335, NOAA-NASA-USAF, Washington, DC, USA. 1976. Available online: https://ntrs.nasa.gov/citations/19770009539 (accessed on 21 November 2022).
  8. Lowry, J.; Polhamus, E. A Method for Predicting Lift Increments Due to Flap Deflection at Low Angles of Attack in Incompressible Flow; NASA TN-3911; NACA: Washington, DC, USA, 1957. [Google Scholar]
  9. Raymer, D.P. Aircraft Design—A Conceptual Approach. In AIAA Education Series, 2nd ed.; American Institute of Aeronautics and Astronautics Inc.: Washington, DC, USA, 1992. [Google Scholar]
  10. Gur, O.; Mason, W.H.; Schetz, J.A. Full Configuration Drag Estimation. In Proceedings of the AIAA Applied Aerodynamics Conference, San Antonio, TX, USA, 23–25 June 2009; AIAA: Reston, VA, USA, 2009; pp. 2009–4109. [Google Scholar]
  11. Nita, M.; Scholz, D. Estimating the Oswald Factor from Basic Aircraft Geometrical Parameters. In Proceedings of the Publikationen zum DLRK 2012 (Deutscher Luft-und Raumfahrtkongress, Berlin, Germany, 10–12 September 2012. [Google Scholar]
  12. Hörner, S.F. Fluid Dynamic Drag; Otterbein Press: Dayton, OH, USA, 1951. [Google Scholar]
  13. Bourdin, P. Numerical Predictions of Wing-Tip Effects on Lift Induced Drag. In Proceedings of the 23rd ICAS Congress, Toronto, ON, Canada, 8–13 September 2002. [Google Scholar]
  14. Delavenne, M.; Barriety, B.; Vetrano, F.; Ferrand, V.; Salaun, M. Parametric Analysis of an Active Winglet Concept for High Aspect Ratio Wing using CFD/CSM Computations. In Proceedings of the AIAA Aviation Forum, Virtual, 15–19 June 2020. [Google Scholar]
  15. Reneaux, J. Réduction de la traînée des avions de transport par laminarité étendue. In Proceedings of the AAAF 24e Colloque Aérodynamique Appliquée, Poitiers, France, 26–28 October 1987. [Google Scholar]
  16. Cousteix, J. Turbulence et Couche Limite; Cepadues Editions: Toulouse, France, 1989; ISBN 2854282108. [Google Scholar]
  17. Haftman, B.; Debbeler, F.J.; Gielen, H. Takeoff Drag Prediction for an Airbus A300–600 and A310 Compared with Flight Test Results. AIAA J. Aircr. 1988, 25, 1088–1096. [Google Scholar] [CrossRef]
  18. Torenbeek, E. Synthesis of Subsonic Airplane Design; Delft University Press: Delft, The Netherlands, 1986. [Google Scholar]
  19. Wiart, L.; Atinault, O.; Paluch, B.; Hue, D.; Grenon, R. Development of NOVA Aircraft Configurations for Large Engine Integration Studies. In Proceedings of the AIAA Paper, Dallas, TX, USA, 22–26 June 2015; pp. 2015–2254. [Google Scholar]
  20. Destarac, D. Drag Extraction from Numerical Solutions to the Equations of Fluid Dynamics: The Far-Field Philosophy. In Proceedings of the 43ème Colloque d’Aérodynamique Appliquée AAAF, Poitiers, France, 10–12 March 2008. [Google Scholar]
  21. Vassberg, J.C.; DeHaan, M.A.; Rivers, S.K.; Wahls, R.A. Development of a Common Research Model for Applied CFD Validations Studies. In Proceedings of the AIAA Paper, Honolulu, HI, USA, 18–21 August 2008; pp. 2008–6919. [Google Scholar]
  22. Rivers, M.B.; Dittberner, A. Experimental Investigations of the NASA Common Research Model in the NASA Langley NTF and the NASA Ames 11-ft Transonic Wind Tunnel, AIAA 2011-1126. In Proceedings of the 49th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, Orlando, FL, USA, 4–7 January 2011. [Google Scholar]
  23. Hue, D.; Esquieu, S. Computational Drag Prediction of the DPW4 Configuration using the Far-Field Approach. J. Aircr. 2011, 48, 1658–1670. [Google Scholar] [CrossRef] [Green Version]
  24. Lynde, M.; Campell, R. Computational Design and Analysis of a Transonic Natural Laminar Flow Wing for a Wind Tunnel Model. In Proceedings of the AIAA Aviation Forum 2017, 35th AIAA Applied Aerodynamics Conference, Denver, CO, USA, 5–7 June 2017. [Google Scholar] [CrossRef] [Green Version]
  25. Rivers, M.; Lynde, M.; Campbell, R.; Viken, S.; Chan, D.; Watkins, A.N.; Goodliff, S. Experimental Investigation of the NASA Common Research Model with a Natural Laminar Flow Wing in the NASA Langley National Transonic Facility. In Proceedings of the AIAA SciTech 2019 Forum, San Diego, CA, USA, 7–11 January 2019. [Google Scholar] [CrossRef]
  26. Lynde, M.; Campbell, R.; Viken, S. Additional Findings from the Common Research Model Natural Laminar Flow Wind Tunnel Test. In Proceedings of the AIAA AVIATION Forum 2019, Dallas, TX, USA, 17–21 June 2019. [Google Scholar] [CrossRef]
  27. de Gaspari, A.; Moens, F. Aerodynamic Shape Design and Validation of an Advanced High-Lift Device for a Regional Aircraft with Morphing Droop Nose. Hindawi Int. J. Aerosp. Eng. 2019, 7982168. [Google Scholar] [CrossRef] [Green Version]
  28. Moens, F. Augmented Aircraft Performance by the use of Morphing Technology for a Turboprop Regional Aircraft Wing. Biomimetics 2019, 4, 64. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Dimino, I.; Andreutti, G.; Moens, F.; Fonte, F.; Pecora, R.; Concilio, A. Integrated Design of a Morphing Winglet for Active Load Control and Alleviation of Turboprop Regional Aircraft. Appl. Sci. 2021, 11, 2439. [Google Scholar] [CrossRef]
  30. Moens, F.; Perraud, J.; Krumbein, A.; Toulorge, T.; Iannelli, P.; Hanifi, A. Transition Prediction and Impact on 3D High-Lift Wing Configuration. AIAA J. Aircr. 2008, 45, 1751–1766. [Google Scholar] [CrossRef]
  31. Meheut, M.; Grenon, R.; Carrier, G.; Defos, M.; Duffau, M. Aerodynamic Design of Transonic Flying Wing Configurations. In Proceedings of the CEAS/KAT net II Conference on Key Technologies, Bremen, Germany, 11–12 May 2009. [Google Scholar]
  32. Meheut, M.; Arntz, A.; Carrier, G. Aerodynamic Shape Optimizations of a Blended Wing Body Configuration for Several Wing Planforms. In Proceedings of the 30th Applied Aerodynamics Conference, New Orleans, LA, USA, 25–28 June 2012; AIAA: Reston, VA, USA, 2012; pp. 2012–3122. [Google Scholar]
  33. Bennehard, Q. Aerodynamic Shape Optimization of a Short Medium Range Blended Wing Body Aircraft. In Proceedings of the 56th 3AF International Conference on Applied Aerodynamics, Toulouse, France, 28–29 March 2022. [Google Scholar]
Figure 1. The different fidelity levels for flow physics modeling of the aerodynamic solvers.
Figure 1. The different fidelity levels for flow physics modeling of the aerodynamic solvers.
Aerospace 10 00007 g001
Figure 2. ONERA MDAO/OAD tool for BWB applications.
Figure 2. ONERA MDAO/OAD tool for BWB applications.
Aerospace 10 00007 g002
Figure 3. Extract of shapes considered for a DoE for a BWB design.
Figure 3. Extract of shapes considered for a DoE for a BWB design.
Aerospace 10 00007 g003
Figure 4. Evolution of the temperature with the altitude—ISA model.
Figure 4. Evolution of the temperature with the altitude—ISA model.
Aerospace 10 00007 g004
Figure 5. Geometrical data used for wing definition. (a) Geometrical data related to wing segments. (b) Geometrical data related to airfoil sections.
Figure 5. Geometrical data used for wing definition. (a) Geometrical data related to wing segments. (b) Geometrical data related to airfoil sections.
Aerospace 10 00007 g005
Figure 6. Fuselage arrangements with similar aerodynamic performance. (a) Front view. (b) Top view.
Figure 6. Fuselage arrangements with similar aerodynamic performance. (a) Front view. (b) Top view.
Aerospace 10 00007 g006
Figure 7. Nacelle types considered in the module.
Figure 7. Nacelle types considered in the module.
Aerospace 10 00007 g007
Figure 8. Lift-induced drag formulation for swept wings.
Figure 8. Lift-induced drag formulation for swept wings.
Aerospace 10 00007 g008
Figure 9. Definitions of the winglet parameters.
Figure 9. Definitions of the winglet parameters.
Aerospace 10 00007 g009
Figure 10. Maximum possible laminar flow extent as a function of Reynolds number and wing sweep angle. (a) NLF. (b) HLF.
Figure 10. Maximum possible laminar flow extent as a function of Reynolds number and wing sweep angle. (a) NLF. (b) HLF.
Aerospace 10 00007 g010
Figure 11. Simplified model used for transition location on wing surfaces. (a) Model for transition change on wing suction side. (b) Model for transition change on wing pressure side.
Figure 11. Simplified model used for transition location on wing surfaces. (a) Model for transition change on wing suction side. (b) Model for transition change on wing pressure side.
Aerospace 10 00007 g011
Figure 12. Estimation of Θ at airfoil trailing-edge for transitional flow.
Figure 12. Estimation of Θ at airfoil trailing-edge for transitional flow.
Aerospace 10 00007 g012
Figure 13. Definition of the ZNAC parameter. (a) “Standard”. (b) “Buried”.
Figure 13. Definition of the ZNAC parameter. (a) “Standard”. (b) “Buried”.
Aerospace 10 00007 g013
Figure 14. ONERA NOVA Aircraft.
Figure 14. ONERA NOVA Aircraft.
Aerospace 10 00007 g014
Figure 15. NOVA configuration—winglet effect (δWLT = −18°). Lines—aerodynamic module; symbols—CFD.
Figure 15. NOVA configuration—winglet effect (δWLT = −18°). Lines—aerodynamic module; symbols—CFD.
Aerospace 10 00007 g015
Figure 16. NOVA configurations used for UHBR engine integration.
Figure 16. NOVA configurations used for UHBR engine integration.
Aerospace 10 00007 g016
Figure 17. Common research model (CRM) reference configuration.
Figure 17. Common research model (CRM) reference configuration.
Aerospace 10 00007 g017
Figure 18. CRM wing—data used by the aerodynamic module.
Figure 18. CRM wing—data used by the aerodynamic module.
Aerospace 10 00007 g018
Figure 19. CRM configuration—fuselage effect on the CL(α) curve. Lines—aerodynamic module; symbols—CFD.
Figure 19. CRM configuration—fuselage effect on the CL(α) curve. Lines—aerodynamic module; symbols—CFD.
Aerospace 10 00007 g019
Figure 20. CRM configuration—CL(CD) curves (M = 0.85, Re = 5 × 106). Lines—aerodynamic module; symbols—CFD.
Figure 20. CRM configuration—CL(CD) curves (M = 0.85, Re = 5 × 106). Lines—aerodynamic module; symbols—CFD.
Aerospace 10 00007 g020
Figure 21. CRM configuration—drag components. Lines—aerodynamic module; symbols—CFD. (a) Lift-induced. (b) Wave and viscous.
Figure 21. CRM configuration—drag components. Lines—aerodynamic module; symbols—CFD. (a) Lift-induced. (b) Wave and viscous.
Aerospace 10 00007 g021
Figure 22. CRM configuration—Nacelle (TFN) effect. Lines—aerodynamic module; symbols—experiments.
Figure 22. CRM configuration—Nacelle (TFN) effect. Lines—aerodynamic module; symbols—experiments.
Aerospace 10 00007 g022
Figure 23. CRM-NLF geometry used by the aerodynamic module.
Figure 23. CRM-NLF geometry used by the aerodynamic module.
Aerospace 10 00007 g023
Figure 24. CRM-NLF configuration—estimated maximum transition location by the module (M = 0.8565, Re = 15 × 106, CL = 0.50).
Figure 24. CRM-NLF configuration—estimated maximum transition location by the module (M = 0.8565, Re = 15 × 106, CL = 0.50).
Aerospace 10 00007 g024
Figure 25. CRM-NLF configuration—transition line estimated by the module at CL = 0.425—comparison with experimental infrared image.
Figure 25. CRM-NLF configuration—transition line estimated by the module at CL = 0.425—comparison with experimental infrared image.
Aerospace 10 00007 g025
Figure 26. CRM-NLF configuration—global aerodynamic performances (M = 0.8565, Re = 15 × 106). Lines—aerodynamic module; symbols—experiments. (a) CL(α). (b) CL(CD).
Figure 26. CRM-NLF configuration—global aerodynamic performances (M = 0.8565, Re = 15 × 106). Lines—aerodynamic module; symbols—experiments. (a) CL(α). (b) CL(CD).
Aerospace 10 00007 g026
Figure 27. AG2-NLF configuration—data used by the aerodynamic module (4 wing sections).
Figure 27. AG2-NLF configuration—data used by the aerodynamic module (4 wing sections).
Aerospace 10 00007 g027
Figure 28. AG2-NLF configuration—global aerodynamic performances (M = 0.52, Re = 17.5 × 106). Lines—aerodynamic module; symbols—CFD. (a) CL(α). (b) CL(LoD).
Figure 28. AG2-NLF configuration—global aerodynamic performances (M = 0.52, Re = 17.5 × 106). Lines—aerodynamic module; symbols—CFD. (a) CL(α). (b) CL(LoD).
Aerospace 10 00007 g028
Figure 29. AVECA BWB configuration—data used by the aerodynamic module (5 wing sections).
Figure 29. AVECA BWB configuration—data used by the aerodynamic module (5 wing sections).
Aerospace 10 00007 g029
Figure 30. AVECA BWB configuration—global aerodynamic performances (M = 0.85, Re = 180 × 106). Lines—aerodynamic module; symbols—CFD. (a) CL(α). (b) CL(CD).
Figure 30. AVECA BWB configuration—global aerodynamic performances (M = 0.85, Re = 180 × 106). Lines—aerodynamic module; symbols—CFD. (a) CL(α). (b) CL(CD).
Aerospace 10 00007 g030
Figure 31. AVECA BWB configuration—drag components. Lines—aerodynamic module; symbols—CFD. (a) Lift-induced. (b) Wave and viscous.
Figure 31. AVECA BWB configuration—drag components. Lines—aerodynamic module; symbols—CFD. (a) Lift-induced. (b) Wave and viscous.
Aerospace 10 00007 g031
Figure 32. SMILE BWB configuration—data used by the aerodynamic module (6 wing sections).
Figure 32. SMILE BWB configuration—data used by the aerodynamic module (6 wing sections).
Aerospace 10 00007 g032
Figure 33. SMILE BWB configuration—global aerodynamic performances (M = 0.78, Alt. = 41,000 ft). Lines—aerodynamic module; symbols—CFD. (a) CL(α). (b) CL(CD).
Figure 33. SMILE BWB configuration—global aerodynamic performances (M = 0.78, Alt. = 41,000 ft). Lines—aerodynamic module; symbols—CFD. (a) CL(α). (b) CL(CD).
Aerospace 10 00007 g033
Figure 34. SMILE BWB configuration—drag components. Lines—aerodynamic module; symbols—CFD. (a) Lift-induced. (b) Wave and viscous.
Figure 34. SMILE BWB configuration—drag components. Lines—aerodynamic module; symbols—CFD. (a) Lift-induced. (b) Wave and viscous.
Aerospace 10 00007 g034
Table 1. Analytic standard atmosphere model used.
Table 1. Analytic standard atmosphere model used.
LayerAltitude (km)Temperature (K)Pressure (Pa)
TroposphereZ ≤ 11 T = 288.15 6.5 × Z P = 101325.0 × ( T 288.15 ) 5.2561
Stratosphere11 ≤ Z ≤ 20 T = 216.65   P = 22630.6 × 10 ( z 11 14.596 )
20 ≤ Z ≤ 32 T = 216.65 ( Z 20 )
Table 2. Form factors and wetted area considered in the aerodynamic module.
Table 2. Form factors and wetted area considered in the aerodynamic module.
ElementForm FactorSWet
Wing [17] F F W I N G = [ 3.4004 ( t c ) 0.4578 ( t c ) 2 + 13.0119 ( t c ) 3 ]   c o s 2 φ 50 + 1 2 SWING
Fuselage [9] F F F U S = 1 + 60 ( L F U S D F U S ) 3 + 0.0025 ( L F U S D F U S ) π LFUS DFUS
Winglets [17] F F W L T = 1 + 3.52 ( t c ) W L T c o s   φ W L T 2 SWLT
Tail surfaces [17] F F T A I L = 1 + 3.52 ( t c ) T A I L c o s   φ T A I L 2 STAIL
Nacelles [18] F F N A C = 1 + 0.35 ( D N A C L N A C ) 2π DNAC LNAC
Table 3. NOVA configuration—drag increase (drag counts) for different engine installations.
Table 3. NOVA configuration—drag increase (drag counts) for different engine installations.
ZNAC (m)Aerodynamic ModuleCFD
CL = 0.50 CDΔCDΔCD
Ref. (No Engine)--210----
Baseline+10.00232+22+28
Podded+0.75244+34+34
BLI−0.50226+16+18
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

Moëns, F. A Fast Aerodynamic Model for Aircraft Multidisciplinary Design and Optimization Process. Aerospace 2023, 10, 7. https://doi.org/10.3390/aerospace10010007

AMA Style

Moëns F. A Fast Aerodynamic Model for Aircraft Multidisciplinary Design and Optimization Process. Aerospace. 2023; 10(1):7. https://doi.org/10.3390/aerospace10010007

Chicago/Turabian Style

Moëns, Frédéric. 2023. "A Fast Aerodynamic Model for Aircraft Multidisciplinary Design and Optimization Process" Aerospace 10, no. 1: 7. https://doi.org/10.3390/aerospace10010007

APA Style

Moëns, F. (2023). A Fast Aerodynamic Model for Aircraft Multidisciplinary Design and Optimization Process. Aerospace, 10(1), 7. https://doi.org/10.3390/aerospace10010007

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