Next Article in Journal
Optimizing the Location of the Piezoelectric Actuator and Analyzing Its Effect on the Dynamics of Asymmetric Flexible Spacecraft
Next Article in Special Issue
Robust Design Optimization of Supersonic Biplane Airfoil Using Efficient Uncertainty Analysis Method for Discontinuous Problem
Previous Article in Journal
Nonlinear Time Series Analysis and Prediction of General Aviation Accidents Based on Multi-Timescales
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Low-Boom Design for Supersonic Transport with Canard and Forward-Swept Wings Using Equivalent Area Design Method

Department of Aeronautics and Astronautics, Graduate School of System Design, Tokyo Metropolitan University, Tokyo 191-0065, Japan
*
Author to whom correspondence should be addressed.
Current address: Aviation System Research Unit, Aviation Technology Directorate, Japan Aerospace Exploration Agency (JAXA), Tokyo 181-0015, Japan.
Current address: 2nd Space Systems Group, Aerospace Solutions Department, Aerospace Business Division, NEC Corporation, Tokyo 162-0055, Japan.
Aerospace 2023, 10(8), 717; https://doi.org/10.3390/aerospace10080717
Submission received: 13 June 2023 / Revised: 14 August 2023 / Accepted: 15 August 2023 / Published: 16 August 2023
(This article belongs to the Special Issue Research and Development of Supersonic Aircraft)

Abstract

:
Forward-swept wings can be expected to be lower-boom planforms with similar amount of drag as backward-swept wings because of their good lift distributions. In this study, the equivalent area distribution of a ten-seater supersonic forward-swept wing aircraft with a canard was designed to obtain design knowledge for leading boom reduction. The equivalent area distribution of the aircraft was calculated by solving the compressible Euler equation. A feasible target equivalent area distribution was generated based on Darden’s method and compared with the equivalent area distribution. To achieve a closer match in terms of lift and geometry with the target, the main wing planform and the position of the main wing along the body and vertical axes were modified. The low-boom performances were evaluated using the extended Burgers equation. The design results indicated that the forward-swept wing configuration with a canard could divide the single peak of the leading boom into two peaks. Thus, the sonic boom strength of the canard configuration was 2.5 PLdB lower than that of the configuration without the canard wing.

1. Introduction

Typical supersonic transport (SST) configurations have wings with large backward-swept angles such as a delta, arrow, and backward-swept wings. This design aims to reduce wave drag at transonic and supersonic speeds. Numerous studies regarding the aerodynamic design of SSTs were performed considering the utilization of wings with large backward-swept angle [1,2]. In the United States, the National Aeronautics and Space Administration (NASA) launched the X-59 QuessT (Quiet SuperSonic Technology) project in 2018 to demonstrate advanced low-boom flight technologies. X-59 is equipped with wings with large backward-swept angles. [1] discusses the comparison between computational fluid dynamics (CFD) and wind tunnel experiment concerning the X-59. In [2], a robust low-boom design method using reversed equivalent area based on off-body pressure distribution was proposed to consider off-track sonic boom. This method was applied to an SST configuration with an arrow wing and without a tail wing and engine. However, these wings with large sweep angles have difficulty significantly reducing sonic boom because of the lift caused by the horizontal tail and the tip of the main wing, located near the tail wing. Thus, the backward-swept wing cannot be realized with ideal equivalent cross-sectional area distribution along the longitudinal direction, which is called “Darden’s rule”, for sonic boom reduction. Therefore, innovative and unconventional configurations are expected to reduce drag and sonic boom simultaneously.
Forward-swept wings remain a potential concept for low-drag, low-boom SST configurations. In a supersonic cruise, the aerodynamic center of the forward-swept wing shifts forward. Therefore, the horizontal tail wing requires a positive lift. This indicates that a forward-swept wing can potentially solve the trimming issue [3]. However, there is little research concerning an SST with forward-swept wings. Originally, forward-swept wings had been studied for the improvement in the maneuverability of fighter jets and for the reduction in the drag under the transonic flight as backward-swept wings. In the 1980s, Grumman Aerospace Corporation built an advanced technology demonstrator called “X-29” with a forward-swept wing and it was flight tested in 1991 by NASA and the United States Air Force (USAF) for tomorrow’s fighter jet [4]. From an aerodynamic perspective, the integration of aeroelastically tailored composites, a close coupled, variable incidence canard, relaxed static stability, variable wing camber in conjunction with a thin supercritical airfoil, Reynolds number effects at high angles of attack and the implementation of three surface trimming were demonstrated [5]. In the 2010s, the German Aerospace Center (DLR) proposed a concept of a transonic transport “LamAiR” with forward-swept wings which had low drag due to a natural laminar flow [6].
We performed a parametric study to investigate the relationship between the sonic boom performance and backward-swept/forward-swept angle [7]. The results showed that the forward-swept wing could reduce the sonic boom strength owing to the stepwise-pattern pressure signature of the trailing boom while maintaining an equal level of aerodynamic drag as the backward-swept wing. In addition, optimal airfoil distributions for the forward-swept and backward-swept wings were designed using multi-objective optimization based on the evolutionary algorithm (EA). The forward-swept wing obtained a lower-boom solution than the conventional backward-swept wing [8].
In general, the sonic boom derives from the volume and lift. The sonic boom owing to the volume can be reduced by considering the equivalent area distribution owing to the volume along the body axis. However, there is limited room for the improvement in the equivalent area distribution owing to the volume, for the constraints to provide sufficient cabin and cargo space for passengers. On the other hands, the sonic boom owing to the lift can be reduced by considering the equivalent area distribution owing to the lift along the body axis which was directly affected by the lift distribution. In this study, we focused on the redesign of the lift distribution to reduce the sonic boom.
There are several methodologies for low-boom design. In [8], the sonic boom strength was defined as an objective function as well as drag, wing geometry was directly optimized using the EA. While this approach facilitates the discovery of low boom solutions, improving the local pressure signature artificially without a thorough review of the design space remains challenging. Another technique for low boom design involves a theoretical approach. One of the established theoretical methods is based on the Seebass George Darden (SGD) theory [9,10]. The SGD theory transforms the ideal pressure waveform associated with a low boom into an equivalent area distribution using Whitham’s F function [11]. In [12], the wing planform was optimized to minimize the drag and sonic boom strength. This was accomplished by employing an improved F function that considers trim capability in conjunction with the EA. These theoretical approaches enable the enhancement of the local pressure signature by comparing it to the target equivalent distribution.
From an aerodynamic perspective, there is scope for a forward-swept wing to reduce the sonic boom. This is a stronger pressure-increase that affects the leading part of the boom signature than that of the backward-swept wing. Thus, the optimal wing design could not significantly reduce the drag and sonic boom strength simultaneously [8]. Therefore, for a forward-swept wing, a reduction in the pressure-increase in the boom signature at the leading part or further improvement of the boom signature at the trailing part is required.
One of the leading boom reduction methods is the application of a canard wing configuration. The canard wing generates lift at the front of the body. This allows for a wider range of lift distribution designs. This indicates the feasibility of improving the boom signature at the leading part and reducing the trim drag. Several studies have been conducted on canard wings [13,14,15]. [13] discussed an integrated design method for wing-fuselage canard wing configurations using the adaptive range multi-objective genetic algorithm (ARMOGA) [16]. In [14], the lift-to-drag ratio and sonic boom strength in a supersonic cruise were improved by applying a canard wing to a supersonic business jet (SSBJ) with an arrow wing. In [15], only canard wings were designed considering their structural, thermal, and aeroelastic performance. Although this type of knowledge is effective as a fundamental understanding, there is negligible research focused on the aerodynamics of SSTs with canard and forward-swept wings.
Thus, this study aims to obtain design knowledge for an SST with a forward-swept wing to improve the pressure-increase affecting the leading part of the boom signature. We improved the equivalent area distribution of the SST-applied canard wing using a forward-swept wing. To improve the equivalent area distribution, we compared the target equivalent area distribution based on the SGD theory with the total equivalent area calculated from the CFD-based lift distribution and the cross-sectional distribution along the body axis.
Excluding this introduction, this paper consists of four chapters. In Chapter 2, we describe the details of the SSBJ as the design target. Chapter 3 introduces the methodology of this study, such as the evaluation methods of the aerodynamics and sonic boom and low-boom design method based on the equivalent area distribution. Chapter 4 presents the results and discussion. In this chapter, the equivalent area distributions and pressure signatures for several configurations were compared. Finally, Chapter 5 discusses the conclusions of this study.

2. Design Target

2.1. Design Target

In this study, a modified SSBJ with a forward-swept wing designed in [8] was defined as a “Baseline configuration”. In this study, the airframe configuration without the engine nacelle was considered and the wing planform and mounting position were modified from the SSBJ designed in [8]. The differences between the Baseline configuration and the SSBJ designed in [8] were as follows:
  • To shift the design range of the lift distribution forward for the sonic boom reduction, the horizontal tail wing was deleted, and a canard wing was added. The geometry and the mounting position of the canard wing were obtained from [14]. The airfoil of the canard wing is a biconvex airfoil. This is similar to the horizontal tail wing in [8]. The mounting angle of the canard wing is fixed at 2.0° for simplicity in this study.
  • To reduce the sonic boom behind the wing and maintain the radius of the fuselage, the overall lift should be increased to two times. Thus, the wing area was enlarged while the aspect ratio was maintained. The wing area was determined based on the specifications of an SSBJ with a canard wing [14] at a similar scale.
  • To reduce the difference between the target equivalent area distribution behind the trailing edge of the wing and the total equivalent area distribution, the mounting position of the wing was shifted 1.54 m backward (5% of the body length) along the body axis.
  • To prevent unforeseen variations in the sonic boom signature, the dihedral angle was set to zero.
An overview of the Baseline configuration and specifications of the wing and canard wing are shown in Figure 1, and Table 1 and Table 2, respectively.

2.2. Geometric Definition Guideline of Aircraft

The wing planform for both the main wing and the canard wing was designed under the condition that the equivalent area distribution owing to the lift at the aft of the aircraft was unaltered. This implies that the design was based on a constant wing area to maintain the total lift. There is a risk of inaccurate sonic boom estimation if the geometry varies and the distance between the wing tip and body axis increases after design. Thus, the length of the wing span was fixed by setting the wing area and aspect ratio constant.
The main wing was composed of an inboard wing connected to the fuselage as a strake with a backward-swept angle and an outboard wing with a forward-swept angle. The wing planform was defined using nine design variables, as shown in Figure 2.

3. Low-Boom Design Method Based on Equivalent Area Distribution

3.1. Overview

After calculating the difference d A e between the total equivalent area distribution A e and target equivalent area distribution A t of the design solution, we redesigned the design solution to minimize this difference. A flow chart of the low-boom design method is shown in Figure 3.

3.2. Aerodynamic Evaluation Method

3.2.1. Computational Fluid Dynamics (CFD)

The governing equation of the high-fidelity evaluation was the compressible Euler equation:
t V   q d V + S   H · n d S = 0
where q is the conservative vector, H is the inviscid flux vector, n is the outward normal vector on the boundary surface of the control volume, V is an element of the volume, and S is an element of the surface. A Simple Low-dissipation Advection Upstream Splitting Method (SLAU) scheme [17] was employed to evaluate the numerical function. In addition, the lower-upper symmetric Gaussian–Seidel (LU-SGS) implicit method [18] was used for time integration. The number of calculation grids was approximately 68 million. The FAST Aerodynamic Routines (FaSTAR) [19] was used as the aerodynamic solver. An unstructured hexahedral mesh was generated around the aircraft using HexaGrid [20,21]. It is an automatic mesh-generation software package developed by the Japan Aerospace Exploration Agency (JAXA). These calculations were performed using the JAXA Supercomputer System (JSS). The calculation grid dependency on the sonic boom strength was discussed in Appendix A.

3.2.2. Sonic Boom Evaluation Method

The strength of a sonic boom with a finite-thickness shock wave and pressure signature on the ground was obtained by solving the augmented Burgers equation after applying multipole analysis to the pressure distribution near the aircraft (based on the Euler simulation). In the multipole analysis [22], near-field pressure signatures were modified to attenuate only in the circumferential direction. This was achieved by introducing multipole distributions. Thus, the far-field pressure signature was independent of the location at which the initial condition was obtained. MPnoise [23] was used for this process. It is a multipole analysis tool developed by JAXA.
The augmented Burgers equation [24] is expressed in Equation (2). Here, p is the pressure fluctuation in the atmosphere, s is the coordinate along the ray, γ is the heat capacity ratio, ρ is the atmospheric density, c 0 is the speed of sound, and τ is the retardation time. In the Burgers equation, the gradient of the pressure difference along the ray [left-hand side of Equation (2)] is the sum of the non-linearity of the sound wave’s finite amplitude (the first term on the left side) and attenuation effect of the atmospheric heat viscosity (the second term on the left side). In this equation, B is expressed by Equation (3). V s denotes the Blokhintsev invariant. It is known to be conserved along a ray tube in atmospheric wind in the linear case. In the augmented Burgers equation, the stratification effect of the atmosphere [Equation (4)], the geometric expansion effect of the Mach cone [(Equation (5)] and the attenuation effect of the vibrational relaxation of atmospheric molecules [Equation (6)] are considered in addition to the Burgers equation. Here, ν is the chemical species, c ν is the increment in the speed of sound-related relaxation in ν , and τ ν is the relaxation time in ν .
p s = γ + 1 2 1 2 ρ c 0 3 p 2 τ 1 2 B B s p + a A B e q + b A B e q + c A B e q
B = V s p 2
a A B e q = 1 2 ρ c 0 ρ c 0 s
b A B e q = δ 2 c 0 3 2 p τ 2
c A B e q = ν c ν τ ν c 0 2 1 + τ ν τ 1 2 p τ 2
Xnoise [23] was used for this process. It is a far-field signature prediction tool developed by JAXA. The perceived level (PL) [25] was used as a measure of the sonic boom strength. Figure 4 shows some examples of PLdB noise levels. The PL was estimated using BoomMetre [23]. It is frequency analysis software developed by JAXA for sonic booms. This set of sonic boom evaluation procedures was compared with the results of free-flight experiments, and it showed good agreement [26,27].

3.3. Calculation Method of Equivalent Area Distribution

The total equivalent area distribution A e is the sum of the equivalent area distribution owing to the volume A V and that is owing to the lift A L . This is shown in Equation (7):
A e ( x ) = A V x + A L ( x )
Equation (8) shows the equivalent area distribution owing to the volume. Here, S M is the cross-sectional area of the aircraft cut by the Mach plane, and μ is the angle of the Mach plane.
A V ( x ) = S M ( x ) sin μ
The equivalent area distribution owing to the lift is obtained by integrating the lift distribution, as shown in Equation (9). Here, M is the Mach number under cruise flight, U is the cruising speed, and L is the overall lift of the aircraft.
A L ( x ) = M 2 1 ρ U 2 0 x L x d x

3.4. Calculation Method of Target Equivalent Area Distribution

Seebass, George, and Darden calculated the target equivalent area distribution based on the sonic boom minimization theory, considering the non-uniformity of the atmosphere and blunt body by [9,10]. In this theory, Whitham’s F function is defined to convert an ideal low sonic boom signature into an equivalent area distribution. A designer can generate an arbitrary target distribution for the design principle by setting parameters such as the mean length, weight, cruising speed, atmospheric density, and temperature.
Using the F function, the target equivalent area distribution can be expressed as shown in Equation (10).
A t ( x ) = 4 0 x F ( y ) · x y d y
Thus, the difference between the equivalent area distribution of the aircraft and the target equivalent distribution d A e is given in Equation (11). Here, l is the aircraft length.
d A e = 0 l A e x A t ( x ) d x
In the Baseline configuration, the target equivalent distribution using the mean length between the nose and the maximum A V position and that using the mean length between the nose and position of the trailing edge of the wing were different. This was because the equivalent area distribution adopted the maximum value in front of the trailing edge of the wing and decreased behind the trailing edge of the wing. In this study, the lift distribution in front of the trailing edge of a wing was designed by modifying the wing planform and amount of lift shared between the main and canard wings. Therefore, to calculate the target equivalent area distribution, we defined the length from the nose to the trailing edge of the wing as the mean length. The equivalent area at this position was reflected by the weight.

4. Results and Discussion

4.1. Equivalent Area Distribution of Baseline Configuration and Target Equivalent Area Distribution

Figure 5 shows the equivalent area distribution of the Baseline configuration owing to the lift, that owing to the volume, the total equivalent area distribution of this configuration, and the target equivalent area distribution. After the total equivalent area distribution increased monotonically to 10 m, because of the increase in the equivalent area distribution owing to the volume, the gradient of the total equivalent area distribution increased from 10 m to 12 m. This occurred because the canard wing generated lift, and the equivalent area distribution owing to the lift increased. Behind the canard wing, the gradient of the total equivalent area distribution decreased, and the total equivalent area distribution increased monotonically.
From the leading edge of the wing, located at a distance of approximately 15 m, the gradient of the equivalent area distribution became large. This occurred because the equivalent area distribution owing to the lift caused by the wing was added to that owing to the volume, which increased monotonically from the nose. From the leading edge of the kink of the wing (located at a distance of approximately 21 m), the equivalent area distribution owing to the volume decreased, and the gradient of the equivalent area distribution owing to the lift increased. At this time, the gradient of the total equivalent area distribution remained constant. This was because of the decrease in that owing to the volume and the increase in that owing to the lift canceled out. The gradient of the equivalent area distribution owing to the lift became gentle at 24 m. Thus, the total equivalent area distribution attained a maximum value near 25 m. This was because the decrease in the equivalent area distribution owing to the volume was larger than the increase in that owing to the lift. Until the trailing edge of the wing was located at approximately 28 m, the decrease in the equivalent area distribution owing to the volume and the increase in that owing to that lift continued to cancel each other out, and the total equivalent area distribution decreased gradually. The total equivalent area distribution decreased linearly with the decrease in the equivalent area distribution owing to the volume. This occurred because the equivalent area distribution owing to lift became constant owing to the absence of a horizontal tail wing behind 28 m. Finally, the total equivalent area distribution became constant behind 30.77 m (the tail end of the aircraft).

4.2. Target Equivalent Area Distribution and Redesign Principle

Figure 6 shows the total equivalent area distribution for the Baseline configuration A e , target equivalent area distribution A t , and the difference between them along the body axis d A e . Four large regions of d A e were identified: regions A (0–9 m), B (9–20 m), C (20–28 m), and D (≥28 m). To reduce the strength of the sonic boom, we redesigned the aircraft geometry by minimizing the area enclosed by the d A e curve and body axis. However, region D was not considered in this study. This was because it was not feasible to reduce the area enclosed by the d A e curve and body axis in this region without redefining the weight of the aircraft and target equivalent area distribution.
Region A was located between the nose and the leading edge of the canard wing. Therefore, we attempted to reduce the area enclosed by the d A e curve and body axis in this region by shifting the lift caused by shifting the canard wing forward. Region B was located between the rear of the canard wing and the leading edge of the kink on the wing. To reduce the area enclosed by the d A e curve and body axis in this region, it was necessary to increase the lift at the front of the wing because A e A t . Region C was located between the leading and trailing edges of the kink in the wing. In this region, we redesigned the geometry of the aircraft to decrease the maximum value of A e because A e A t .
To achieve this, we sequentially redesigned the geometry of the aircraft using the following three steps for each of the three regions:
  • For region A, we shifted the position of the canard wing forward from the Baseline configuration and then shifted the lift distribution of the canard wing forward.
  • For region B, while maintaining the wing span length constant, we moved the position of the kink inside the wing to increase the wing area of the outboard wing with a forward-swept angle and shift the lift distribution forward. Simultaneously, we increased the chord length of the kink and modified that of the root to maintain a constant wing area.
  • For region C, we further stretched the chord length of the kink and expanded the wing area of the forward-swept wing to reduce the maximum value of A e and make its gradient gentler. This also shifted the lift distribution forward and strengthened the effect of Step 2. In addition, we decreased the leading edge swept-back angle and the chord length of the root to decrease the wing area of the inboard wing. This canceled out the increase in the outboard wing area. We also moved the wing-mounted position backward to shift the lift distribution caused by the wing backward and further reduce the maximum value of A e .

4.3. Local Lift Distributions of Redesigned Configurations

To improve regions A, B, and C in Figure 6, the redesign was conducted three times based on Steps 1, 2, and 3 in the previous section. We named these configurations “desC1”, “desW1”, and “desW2”, respectively. In addition, a configuration without the canard wing, “w/o canard”, was evaluated.
Figure 7 shows the top view of each configuration. The w/o canard configuration consists of a fuselage, wing, and vertical tail wing. Meanwhile, the other configurations include a canard wing, fuselage, wing, and vertical tail wing. The canard wing of the desC1 configuration is mounted forward compared with that of the Baseline configuration. The desW1 configuration has a kink in the wing located inboard the Baseline configuration and has a shorter chord length at the root than the Baseline configuration. The desW2 configuration has a longer chord length at the kink and a smaller inboard leading edge swept-back angle than the desW1 configuration. Therefore, the wing mounting position of the desW2 configuration is located further back than that of the desW1 configuration.
The equivalent area distribution owing to the volume of each configuration is shown in Figure 8. According to this figure, the volume difference before and after the redesign is marginal. The wing volumes of the desW1 and desW2 configurations decreased because the chord length of the wing root, which was thick, decreased.
Figure 9 shows the local lift distribution along the body axis of each configuration. In this figure, the first peak of each configuration with the canard wing corresponds to the lift generated by the canard wing. The peak lift for the desC1, desW1, and desW2 configurations was located forward compared with that for the Baseline configuration because the canard wings of these configurations shifted forward. In this study, the area of the main wing, wing planform, and the geometry of the canard wing did not vary. Thus, the lift distribution for all the configurations, which corresponded to the lift generated by the wing, had similar shapes. In the w/o canard configuration, although a lift peak between 7–15 m owing to the canard wing was not observed, the asymmetric fuselage generated a lift of approximately 10 m.
Focusing on the lift distribution at approximately 9–21 m (which corresponds to region B of the d A e distribution in Figure 6), the shapes of the lift distributions of the desW1 and desW2 configurations (which were modified by the wing planform) were different from those of the w/o canard, Baseline, and desC1 configurations. The lift generated by the outboard wing of the desW1 configuration was stronger than that of the desC1 configuration because the kink of the desW1 configuration shifted away from the desC1 configuration. Then, the area of the outboard wing increased. At this time, the wing tip of the desW1 configuration shifted forward from the desC1 configuration and the position at which the lift started shifted from 14 to 13 m. In the desW2 configuration, a local minimum lift was observed at approximately 18 m. Although both desW1 and desW2 configurations had redesigned wings, a local minimum was observed only in the desW2 configuration because of the decreased area of the inboard wing. This occurred because the chord length and leading edge swept-back angle of the inboard wing of desW2 decreased compared with those of desW1. Comparing the local lift distributions of the desC1 and Baseline configurations at approximately 14–17 m, the lift of the desC1 configuration was larger than that of the Baseline configuration. The only difference between the desC1 and Baseline configurations is the mounting position of the canard wing. Thus, the difference in local lift distribution between the desC1 and Baseline configurations could be attributed to the aerodynamic effect of the canard wing on the leading edge of the forward-swept wing.
Focusing on region C of the d A e distribution in Figure 6, the maximum value of the local lift of the desW2 configuration was smaller than those of the configurations. This was particularly so when the position attained a maximum value of 25 m. Considering that the maximum value of the total equivalent area distribution of the Baseline configuration was 25 m (as shown in Figure 6), redesigning desW2 was expected to reduce the maximum value of the total equivalent area distribution. Because the wing-mounting position shifted backward, the local lift of approximately 25–28 m in the desW2 configuration was larger than those of the other configurations. Expanding the area of the outboard wing with a forward-swept angle and shifting the wing backward caused the desW2 configuration to lift more than the other configurations at approximately 13–17 m and 25–28 m. The widely distributed lift along the body axis reduced the maximum value of the local lift distribution.
Figure 10 shows the C M of each configuration. C M is within 0.016 for any configuration. In the future, the sonic boom and aerodynamic performance of SSTs with the forward-swept wing and canard wing should be investigated, considering the pitching trim and the balance between gravity and lift.

4.4. Equivalent Area Distributions of Redesigned Configurations

Figure 11 shows the equivalent area distribution owing to the lift along the body axis of each configuration. The shape of A L was modified by altering the local lift distribution. Focusing on the nose to 9 m corresponding region A in Figure 6, the differences in the total equivalent distribution originated between the Baseline and Target distributions because A L was zero. That is, the total equivalent area consisted of only A V in this range. However, in the desC1, desW1, and desW2 configurations, the canard wings shifted forward. Thus, the equivalent area distribution owing to the lift increased from 7 m to 11 m. This expanded the design range of A L . A comparison of each A L distribution from 13 to 21 m corresponding to region B in Figure 6 revealed that the A L of the desW1 and desW2 configurations expanded the outboard wing area. The forward-swept angle was larger than that of the other configurations. The outboard wing area of the desW2 configuration expanded more than that of the desW1 configuration. Moreover, the A L of the desW2 configuration was higher than that of the desW1 configuration within this range. Focusing on the A L distributions from 21 to 28 m corresponding to region C in Figure 6, no remarkable differences were observed between the configurations in region B. The exception was that the A L distributions of the desW2 configuration were smaller than those of the other configurations in this range.
The total equivalent area distribution for each configuration is shown in Figure 12. Comparing the desW2 and Baseline configurations, the shape of the desW2 configuration was closer to the target equivalent area distribution from 7 to 28 m. To validate the effect of the canard wing, the w/o canard and Baseline configurations were focused on. At approximately 11–22 m, the shape of the Baseline configuration was closer to that of the target distribution than that of w/o canard configuration. This indicated that the application of a canard wing can achieve a good total equivalent area distribution. Comparing the Baseline and desC1 configurations (in which the canard wing was shifted forward from the Baseline configuration to improve region A in Figure 6), the A e shape of the desC1 configuration was closer to that of the target distribution from 7 to 12 m than that of the Baseline configuration because of the increase in A L owing to the movement of the canard wing.
A comparison between the desC1 and desW1 configurations (in which the position of the kink shifted to the root of the wing from the desC1 configuration to improve region B in Figure 6) shows that the A e shape of the desW1 configuration was closer to that of the target distribution from 13 to 21 m than that of the desC1 configuration because of the increase in A L . This variation in A e was smaller than that in A e from w/o canard to the Baseline, Baseline to desC1, and desW1 to desW2.
We compared the A e shapes of the desW1 configuration and desW2 configuration whose outboard wing area was expanded, inboard wing area was contracted, and wing planform was swept back from the desW1 configuration to improve region C (see Figure 6). In the region behind 9 m, the A e shape of the desW2 configuration was closer to that of the target distribution than that of the desW1 configuration because of the increase in the lift from the canard wing. From 14 to 18 m, the A e of the desW2 configuration was larger than that of the desW1 configuration because of the increase in the outboard wing area. From 22 to 27 m, the maximum value of A e for the desW2 configuration decreased because the lift distribution of the desW2-configuration axis widened along the body owing to the redesign of the wing.
The d A e distribution for each configuration is shown in Figure 13. Regions A, B, and C of the d A e distribution decreased in the following order: w/o canard, Baseline, desC1, desW1 and desW2 configurations. Moreover, the A e distribution approached the target distribution. To quantitatively evaluate the difference between A e and A t , the integrated value of d A e along the body axis was calculated (as listed in Table 3). According to this table, the integrated value of the d A e distribution along the body axis decreased with the performance of the redesign process based on the principles described in Section 4.2. If the integrated value of the d A e distribution is used as an objective function of the optimization problem that can be solved using an evolutionary algorithm (EA), a low-boom configuration is likely to be obtained based on the SGD theory.

4.5. Near-Field Pressure Signature of Redesigned Configuration

The near-field pressure signatures of desW2 and Baseline configurations are shown in Figure 14. In both configurations, similar magnitudes of pressure-increase can be observed from 0.8 to 8.0 m (which corresponds to the location of the aircraft’s nose). Intense pressure increases were observed at different positions owing to each canard wing behind the nose.
A pressure increase occurred in the Baseline configuration at approximately 16 m that was similar to the pressure increase owing to the nose. In the desW2 configuration, this type of pressure increase shifted to 18 m. This was caused by a shock wave generated from the leading edge of the wing root. In the desW2 configuration, the leading edge swept-back angle and chord length at the wing root decreased. Thus, the pressure increase in the desW2 configuration shifted backward compared with the Baseline configuration.
The maximum value of p / p was observed at approximately 20 m in the Baseline and desW2 configurations owing to the propagation of the shock wave generated from the wing tip. This peak of the desW2 configuration was located further forward than that of the Baseline configuration. This was because the wing tip of the desW2 configuration was in front of that of the Baseline configuration. Although the desW2 and Baseline configurations had an identical wing span and geometry as the wing tip, the amount of pressure increase owing to the propagation of the shock wave generated from the wing tip differed between these configurations. Thus, this phenomenon was caused by the pressure fluctuations in front of the wing tip. In both configurations, a pressure increase was observed at approximately 23 m owing to the compression wave caused by the kink. The increase in pressure in the desW2 configuration was larger than that in the Baseline configuration. From the discussion regarding the shock wave from the wing tip, the shock wave from the kink of the desW2 configuration was stronger than that of the Baseline configuration because the shock wave from the wing tip of the desW2 configuration became weaker and was located more forward than that of the Baseline configuration. At approximately 26–31 m, the signature of the pressure increase corresponding to the fuselage was different for each configuration because of the difference in the position of the trailing edge of the wing. In the Baseline configuration, the position of the pressure increase on the lower side of the fuselage and the shock wave from the trailing edge of the wing were similar. Thus, after the pressure decrease became gradual at approximately 26 m, the pressure increased dramatically at 29 m. In the desW2 configuration, the trailing edge of the wing shifted backward. Subsequently, the effect of the lower side of the fuselage increased the pressure to approximately 27 m. Moreover, a pressure increase was observed owing to the trailing edge of the wing at approximately 31 m.

4.6. Pressure Signature on the Ground and Perceived Level of Sonic Boom

The pressure signature on the ground for each configuration is shown in Figure 15. Only the w/o configuration had an N-shaped signature. It is a typical sonic boom waveform with a two-stage pressure jump [29]. Canard-wing-mounted configurations such as the Baseline, desC1, desW1, and desW2 configurations had two-stage pressure jumps at the leading boom. Furthermore, their pressure peak values were lower than those of the w/o canard configuration. The only difference between the w/o canard and Baseline (or desC1) configurations was the presence or absence of canard wings. This revealed that the canard wing could divide the leading boom into two stages and reduce its peak pressure.
In the desC1 configuration, where the canard wing shifted forward from the position of the Baseline configuration, the first-stage leading boom also shifted forward. This resulted in a lower peak value compared with the Baseline configuration. The second-stage leading boom of the desC1 configuration was shifted backward and had a lower peak value than that of the Baseline configuration. In the desW1 and desW2 configurations with the modified wings, the position and strength of the first-stage leading boom were almost identical. However, the second-stage leading boom of the desW2 configuration was located behind that of the desW1 configuration and had a lower peak value. This indicates that redesigning the mounted position of the canard wing and wing planform can control the leading boom. This, in turn, can result in an improvement in the sonic boom performance.
Each configuration had a strong single-stage trailing boom. Moreover, the bottom peak values of the trailing boom were almost identical. The bottom peak value of the trailing boom of the w/o canard configuration was the lowest among the evaluated configurations because the pressure jump position of the trailing boom of the w/o canard configuration was located farther back than that of the other configurations. As shown in Figure 15, the Baseline and desW2 configurations exhibited multi-stage trailing booms in the near-field, and their positions and strengths were different. However, these multi-stage trailing booms were unified as single-stage trailing booms during their propagation from the near-field to the ground. To achieve larger sonic boom reduction, it is necessary to convert the single-stage trailing boom at the ground into a multi-stage trailing boom. Several redesigns are required to improve the trailing boom around the near-field. This includes the addition of a horizontal tail to consider the trim stability and the shifting of the position of the trailing edge of the wing.
The perceived level of the sonic boom for each ground configuration is presented in Table 4. The redesign of each configuration resulted in a decrease in d A e and a subsequent decrease in the PLdB value. The Baseline configuration had a PLdB value that was 2.41 PLdB lower than that of the w/o canard configuration because of the multi-stage leading boom that contributed to the effect of the canard wing. The reduction from the w/o canard configuration to the Baseline configuration was larger than that for the other configurations. Dividing a strong leading boom into multiple-stages is important for sonic boom reduction. Thus, the application of a canard wing is effective. The PLdB value of the desW2 configuration was 0.16 PLdB lower than that of the Baseline configuration. Although the d A e value of the desW1 configuration was lower than that of the desC1 configuration; no reduction in the PLdB value was observed. The variation in the total equivalent area distribution (shown in Figure 12) from the desC1 configuration to the desW1 configuration was marginal. Furthermore, no difference in the signature of the first-stage leading boom was observed. In addition, the difference in the signature of the second-stage leading boom was marginal compared with the other redesign processes. Therefore, the variation in wing planform from the desC1 configuration to the desW1 configuration was likely to have been inadequate to affect the PLdB. To reduce the PLdB value further, radical redesign should be required, such as the simultaneous modification of wing and fuselage. If the mounting angle of the canard wing is added to the design variables, the relationship between d A e and PLdB values would vary. This is because the flow field of the canard wing and the equivalent area distribution area owing to the lift around the canard wing affect the flow around the forward-swept wing. At this point, the trim stability should be considered.

5. Conclusions

In this study, to obtain design knowledge for an SST with a forward-swept wing to improve the pressure increase affecting the leading part of the boom signature, the application of a canard wing to an SST with a forward-swept wing was validated. By comparing the equivalent area distribution and target equivalent area distribution based on the SGD theory, the equivalent area distribution was improved by adding the canard wing and redesigning the wing planform and mounted-wing position. Consequently, the sonic boom was reduced by 2.41 PLdB. Three observations were obtained in this study. The canard wing can divide the leading boom on the ground into a multi-stage signature. Considering the mounting position of the canard wing, the pressure jump position of the two-stage leading boom could be altered, and its strength could be controlled. The shape of the equivalent area distribution and its maximum value could be reduced by redesigning the planform of the forward-swept wing. Thus, the position of the second-stage leading boom on the ground and the pressure increase could be improved and reduced. When the lift distribution was redesigned to reduce the difference in the integrated value of the equivalent area distribution along the body axis between the design object and design target, the PLdB values generally decreased with each design.
In this study, for simplicity, the SSTs were redesigned without altering the mounted angle of the canard wing. Next, an integrated design including the fuselage and wing of an SST with a forward-swept wing and canard wing should be carried out using evolutionary computation to minimize sonic boom and aerodynamic drag simultaneously, considering the pitching trim and the balance between the gravity and lift by varying the mounting angle of the canard and the angle of attack of the aircraft. At this time, an equivalent area distribution method that can control leading and trailing booms simultaneously by modifying the geometries of the wing needs to be applied to reduce the sonic boom further.

Author Contributions

Conceptualization, Y.K., R.Y. and M.K.; Formal analysis, R.Y.; investigation, R.Y.; methodology, R.Y.; project administration, M.K.; resource, M.K.; software, R.Y.; visualization, Y.K. and R.Y.; writing—original draft preparation, Y.K.; writing—review and editing, R.Y. and M.K.; supervision, M.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The standard Richardson extrapolation (SRE) [30] was performed to validate the calculation grids. In the SRE, the approximate solution f S R E for the distance between the grids approaches zero and can be estimated using two different calculation grids, namely the fine and coarse grids. The total number of fine grids is, therefore, larger than that of the coarse grids. The approximate solution f S R E can then be expressed as follows:
f S R E = f f i n e + f f i n e f c o a r s e 3
The original SSBJ configuration, developed by JAXA [3], was used to validate the grid. All configurations dealt with in this study are modified from this configuration. The number of fine grids, which had refinement boxes around the aircraft, was approximately 65 million, and the coarse grids, which did not have refinement boxes, was approximately 9.5 million.
The grid dependency of sonic boom performance is shown in Figure A1. The difference between the fine grid and SRE results was approximately 0.25 PLdB, whereas the difference between the coarse grid and the SRE result was more than 1.0 PLdB. Therefore, the fine grid was used for C D P and PL evaluations in this study.
Figure A1. Validation results for grid dependency.
Figure A1. Validation results for grid dependency.
Aerospace 10 00717 g0a1

References

  1. Magee, T.E.; Serrano, L.; Mendoza, R.; Gonzalez, A.; Muscalus, D.R.; Jones, A.M.; Peyton, M.; Shaw, S.G. X-59 CFD and Wind Tunnel Data Comparisons. In Proceedings of the AIAA AVIATION 2023 Forum, San Diego, CA, USA, 12–16 June 2023. [Google Scholar]
  2. Ueno, A.; Kanamori, M.; Makino, Y. Robust Low-Boom Design Based on Near-Field Pressure Signature in Whole Boom Carpet. J. Aircr. 2017, 54, 918–925. [Google Scholar] [CrossRef]
  3. Horinouchi, S. Conceptual Design of a Low Sonic boom SSBJ. In Proceedings of the 43rd AIAA Aerospace Sciences Meeting and Exhibit, Reno, NV, USA, 10–13 January 2005. [Google Scholar]
  4. Weeks, T.M. Advanced Technology Integration for Tommorow’s Fighter Aircraft. In Proceedings of the AIAA Aircraft Systems, Design and Technology Meeting, Dayton, OH, USA, 20–22 October 1986. [Google Scholar]
  5. Moore, M.; Frei, D. X-29 forward Swept Wing Aerodynamic Overview. In Proceedings of the AIAA Applied Aerodynamics Conference, Danvers, MA, USA, 13–15 July 1983. [Google Scholar]
  6. Kruse, M.; Wunderlich, T.; Heinrich, L. A Conceptual Study of a Transonic NLF Transport Aircraft with Forward Swept Wings. In Proceedings of the 30th AIAA Applied Aerodynamics Conference, New Orleans, LA, USA, 25–28 June 2012. [Google Scholar]
  7. Kishi, Y.; Makino, Y.; Kanazaki, M. Numerical Investigation of Sonic Boom for Forward-Swept Wing Based on Euler and Augmented Burgers Equations. Aerosp. Technol. Jpn. Jpn. Soc. Aeronaut. Space Sci. 2020, 19, 1–19. [Google Scholar] [CrossRef]
  8. Kishi, Y.; Kanazaki, M.; Makino, Y. Supersonic Forward-Swept Wing Design Using Multifidelity Efficient Global Optimization. J. Aircr. 2022, 59, 1027–1040. [Google Scholar] [CrossRef]
  9. Darden, C.M. Sonic Boom Theory: Its Status in Prediction and Minimization. J. Aircr. 1977, 14, 569–576. [Google Scholar] [CrossRef]
  10. Darden, C.M. Sonic-Boom Minimization with Nose-Bluntness Relaxation; NASA Technical Paper 1348; NASA Langley Research Center: Hampton, VA, USA, 1979. [Google Scholar]
  11. Whitham, G.B. The Flow Pattern of a Supersonic Projectile. Commun. Pure Appl. Math. 1952, 5, 301–348. [Google Scholar] [CrossRef]
  12. Kasuga, Y.; Yoshida, K.; Ishikawa, H. Wing Planform Optimization Method for Low-Boom and Low-Drag Aircraft. In Proceedings of the AIAA SciTech Forum, Grapevine, TX, USA, 9–13 January 2017. [Google Scholar]
  13. Sasaki, D.; Obayashi, S. Low-Boom Design Optimization for SST Canard-Wing-Fuselage Configuration. In Proceedings of the 16th AIAA Computational Fluid Dynamics Conference, Orlando, FL, USA, 23–26 June 2003. [Google Scholar]
  14. Yoshimoto, M.; Uchiyama, N. Optimization of Canard Surface Positioning of Supersonic Business Jet for Low Boom and Low Drag Design. In Proceedings of the 33rd AIAA Fluid Dynamics Conference and Exhibit, Orlando, FL, USA, 23–26 June 2003. [Google Scholar]
  15. Yetgin, A.; Acar, B. Multidisciolinary Design of Canard. In Proceedings of the 2018 Aviation Technology, Integration, and Operations Conference, Atlanta, GA, USA, 25–29 June 2018. [Google Scholar]
  16. Sasaki, D.; Obayashi, S. Efficient Search for Trade-Offs by Adaptive Range Multi-Objective Genetic Algorithms. J. Aerosp. Inf. Syst. 2005, 2, 44–64. [Google Scholar] [CrossRef]
  17. Shima, E.; Kitamura, K. Parameter-Free Simple Low-Dissipation AUSM-Family Scheme for All Speeds. AIAA J. 2011, 49, 1693–1709. [Google Scholar] [CrossRef]
  18. Sharov, D.; Nakahashi, K. Reordering Hybrid Unstructured Grid for Lower-Upper Symmetric Gauss-Seidel Computations. AIAA J. 1998, 36, 484–486. [Google Scholar] [CrossRef]
  19. Hashimoto, A.; Murakami, K.; Aoyama, T.; Ishiko, K.; Hishida, M.; Sakashita, M. Toward the Fastest Unstructured CFD code “FaSTAR”. In Proceedings of the 50th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, Nashville, TN, USA, 9–12 November 2012. [Google Scholar]
  20. Lahur, P.R. Automatic Hexahedra Grid Generation Method for Component-Based Surface Geometry. In Proceedings of the 17th AIAA Computational Fluid Dynamics Conference, Toronto, ON, Canada, 6–9 January 2005. [Google Scholar]
  21. Hashimoto, A.; Murakami, K.; Aoyama, T. Lift and Drag Prediction Using Automatic Hexahedra Grid Generation Method. In Proceedings of the 47th AIAA Aerospace Sciences Meeting including The New Horizons Forum and Aerospace Exposition, Orlando, FL, USA, 5–8 January 2009. [Google Scholar]
  22. Page, J.A.; Plotkin, K.J. An Efficient Method for Incorporating Computational Fluid Dynamics into Sonic boom Prediction. In Proceedings of the 9th Applied Aerodynamics Conference, Baltimore, MD, USA, 23–25 September 1991. [Google Scholar]
  23. Makino, Y.; Naka, Y.; Hashimoto, A.; Kanamori, M.; Murakami, K.; Aoyama, T. Sonic boom Prediction Tool Development at JAXA. Aeronaut. Space Sci. Jpn. 2013, 7, 237–242. [Google Scholar]
  24. Yamamoto, M.; Hashimoto, A.; Takahashi, T.; Kamakura, T.; Sakai, T. Long-Range Sonic Boom Prediction Considering Atmospheric Effects. In Proceedings of the Inter-Noise 2011, Osaka, Japan, 4–7 September 2011. [Google Scholar]
  25. Stevens, S.S. Perceived Level of Noise by Mark VII and Decibels. J. Acoust. Soc. Am. 1972, 51, 575–601. [Google Scholar] [CrossRef]
  26. Saito, Y.; Ukai, T.; Miyakoshi, K.; Ohtani, K.; Obayashi, S. Sonic Boom Estimation Using the Multipole Method for Free-Flight Experiment. In Proceedings of the AIAA SciTech Forum, National Harbor, MD, USA, 13–17 January 2014. [Google Scholar]
  27. Naka, Y.; Makino, Y.; Hashimoto, A.; Yamamoto, M.; Yamashita, H.; Uchida, T.; Obayashi, S. Validation of Sonic-Boom Propagation Analysis Methods Using D-SEND#1 Data. In Proceedings of the 44th Fluid Dynamics Conference/Aerospace Numerical Simulation Symposium 2012, Toyama, Japan, 5–6 July 2012. [Google Scholar]
  28. Doebler, W.J.; Rathsam, J. How loud is X-59’s shaped sonic boom? In Proceedings of the 177th Meeting of the Acoustical Society of America, Louisville, KY, USA, 13–17 May 2019. [Google Scholar]
  29. Haering, W.A.; Smolka, J.W.; Murray, J.E.; Plotkin, K.J. Flight Demonstration of Low Overpressure N-Wave Sonic Booms and Evanescent Waves. In Proceedings of the 17th International Symposium of Nonlinear Acoustics, State College, PA, USA, 21–22 July 2005. [Google Scholar]
  30. Roy, C.J. Grid Convergence Error Analysis for Mixed-Order Numerical Schemes. AIAA J. 2003, 41, 595–604. [Google Scholar] [CrossRef]
Figure 1. Overview of Baseline configuration.
Figure 1. Overview of Baseline configuration.
Aerospace 10 00717 g001
Figure 2. Design variables for wing planform. (The wing area and aspect ratio were constant.).
Figure 2. Design variables for wing planform. (The wing area and aspect ratio were constant.).
Aerospace 10 00717 g002
Figure 3. Flow chart of the low-boom design method.
Figure 3. Flow chart of the low-boom design method.
Aerospace 10 00717 g003
Figure 4. The Impulsive Noise PL Scale with noise source ranges [28].
Figure 4. The Impulsive Noise PL Scale with noise source ranges [28].
Aerospace 10 00717 g004
Figure 5. Equivalent area distribution of Baseline configuration and target equivalent area distribution.
Figure 5. Equivalent area distribution of Baseline configuration and target equivalent area distribution.
Aerospace 10 00717 g005
Figure 6. Difference between the total equivalent area distribution and the target equivalent area distribution along the body axis.
Figure 6. Difference between the total equivalent area distribution and the target equivalent area distribution along the body axis.
Aerospace 10 00717 g006
Figure 7. Top view of each configuration.
Figure 7. Top view of each configuration.
Aerospace 10 00717 g007
Figure 8. Equivalent area distribution owing to volume of each configuration.
Figure 8. Equivalent area distribution owing to volume of each configuration.
Aerospace 10 00717 g008
Figure 9. Local lift distribution of each configuration.
Figure 9. Local lift distribution of each configuration.
Aerospace 10 00717 g009
Figure 10. Pitching moment coefficient of each configuration.
Figure 10. Pitching moment coefficient of each configuration.
Aerospace 10 00717 g010
Figure 11. Equivalent area distribution owing to lift of each configuration.
Figure 11. Equivalent area distribution owing to lift of each configuration.
Aerospace 10 00717 g011
Figure 12. Total equivalent area distribution of each configuration.
Figure 12. Total equivalent area distribution of each configuration.
Aerospace 10 00717 g012
Figure 13. Difference between A e and A t of each configuration.
Figure 13. Difference between A e and A t of each configuration.
Aerospace 10 00717 g013
Figure 14. Near-field pressure signatures of desW2 and Baseline configurations.
Figure 14. Near-field pressure signatures of desW2 and Baseline configurations.
Aerospace 10 00717 g014
Figure 15. Pressure signature on the ground of each configuration.
Figure 15. Pressure signature on the ground of each configuration.
Aerospace 10 00717 g015
Table 1. Specification of the wing.
Table 1. Specification of the wing.
Wing area [m2] 97.75
Aspect ratio [-] 2.58
Leading edge swept angle [°]Inboard wing76
Outboard wing−30
Trailing edge swept angle [°]Inboard wing0
Outboard wing−65
Table 2. Specification of the canard wing.
Table 2. Specification of the canard wing.
Wing area [m2]6.4
Aspect ratio [-]2.5
Leading edge swept angle [°]45
Trailing edge swept angle [°]5.3
Mounting angle [°]2.0
Table 3. Integrated value of d A e along the body axis of each configuration.
Table 3. Integrated value of d A e along the body axis of each configuration.
Name Integrated   Value   of   d A e ⋅ [m2]
w/o canard20.454
Baseline15.098
desC112.675
desW111.760
desW210.044
Table 4. Perceived level of sonic boom of each configuration.
Table 4. Perceived level of sonic boom of each configuration.
NamePerceived Level [PLdB]
w/o canard102.79
Baseline100.38
desC1100.30
desW1100.30
desW2100.22
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

Kishi, Y.; Yashiro, R.; Kanazaki, M. Low-Boom Design for Supersonic Transport with Canard and Forward-Swept Wings Using Equivalent Area Design Method. Aerospace 2023, 10, 717. https://doi.org/10.3390/aerospace10080717

AMA Style

Kishi Y, Yashiro R, Kanazaki M. Low-Boom Design for Supersonic Transport with Canard and Forward-Swept Wings Using Equivalent Area Design Method. Aerospace. 2023; 10(8):717. https://doi.org/10.3390/aerospace10080717

Chicago/Turabian Style

Kishi, Yuki, Risato Yashiro, and Masahiro Kanazaki. 2023. "Low-Boom Design for Supersonic Transport with Canard and Forward-Swept Wings Using Equivalent Area Design Method" Aerospace 10, no. 8: 717. https://doi.org/10.3390/aerospace10080717

APA Style

Kishi, Y., Yashiro, R., & Kanazaki, M. (2023). Low-Boom Design for Supersonic Transport with Canard and Forward-Swept Wings Using Equivalent Area Design Method. Aerospace, 10(8), 717. https://doi.org/10.3390/aerospace10080717

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