#### Abstract

The linearized Euler equations (LEEs) solver for aeroacoustic problems has been developed on block-structured Cartesian mesh to address complex geometry. Taking advantage of the benefits of Cartesian mesh, we employ high-order schemes for spatial derivatives and for time integration. On the other hand, the difficulty of accommodating curved wall boundaries is addressed by the immersed boundary method. The resulting LEEs solver is robust to complex geometry and numerically efficient in a parallel environment. The accuracy and effectiveness of the present solver are validated by one-dimensional and three-dimensional test cases. Acoustic scattering around a sphere and noise propagation from the JT15D nacelle are computed. The results show good agreement with analytical, computational, and experimental results. Finally, noise propagation around fuselage-wing-nacelle configurations is computed as a practical example. The results show that the sound pressure level below the over-the-wing nacelle (OWN) configuration is much lower than that of the conventional DLR-F6 aircraft configuration due to the shielding effect of the OWN configuration.

#### 1. Introduction

To satisfy the severe requirements of environmentally friendly aircraft, the applications of significant technological advances to conventional aircraft are needed. Presently, some of the most urgent targets are noise, emission, and fuel consumption. Green aviation, which seeks to meet the demands of the noise regulations of the next generation, has been proposed and intensively studied worldwide [1–9]. The continual reduction of noise has been steadily advanced for about 40 years mainly due to improvements in engine performance [4]. Examples of these advances are an increase of the bypass ratio, the use of geared fan machinery, and the improvement of liners and vanes. However, airport noise regulations are becoming increasingly stringent. Aircraft of the next generation must meet these regulations with sufficient margin for future, tougher limits.

The noise generated from aircraft is classified into two main groups, airframe noise and engine noise. Airframe noise is mainly turbulent noise generated from high-lift devices and landing gears. Engine noise is the sum of fan and jet noises. During takeoff, fan and jet noises are dominant because the engine is at full throttle. On the other hand, the airframe and fan noises are dominant during landing. Under both conditions, fan noise is one of the dominant factors of community noise generated from aircraft. To realize a dramatic reduction of engine noise, shielding the fan noise from the engine nacelle by the fuselage or wings is one of the reasonable approaches that have been applied [10–12]. Various unique aircrafts have been proposed and evaluated from this point of view using computational and experimental approaches. To design such unconventional aircraft, it is imperative that noise propagation be reliably analyzed to estimate the shielding effect accurately.

A common computational approach used to analyze noise propagation is based on the linearized Euler equations (LEEs), which is one of the accurate methods available for capturing the propagation of sound [13]. LEEs can address noise propagation from known noncompact sound sources in a nonuniform mean flow field. The method can also address scattering, diffraction, and reflection by objects. LEEs are often solved on body-fitted structured and unstructured meshes [14–17]. From a practical point of view, these meshes have some drawbacks in the computation of realistic complex geometries. Body-fitted structured mesh is fit to the object and is able to adapt to detailed components precisely; however, the computational cost of mesh generation around an entire aircraft or more complex geometry is quite high. Unstructured mesh is applicable to complex geometry; however, it generally achieves low-order accuracy in space. In addition, high-order unstructured methods incur considerable computational cost. Due to these reasons, a simple Cartesian mesh is the focus of the present research. Generation of Cartesian mesh is based on the simple division of the computational domain, and, thus, it can easily accommodate complex geometries. Furthermore, a high-order scheme is easily employed by the extension of stencils. Rapid computation is achieved owing to the simple computational structure of the Cartesian mesh. However, the main problem associated with Cartesian mesh is numerical error caused by the staircase wall boundary. In this research, a block-structured Cartesian mesh method, denoted as the building-cube method (BCM) [18–22], and the immersed boundary method (IBM) are applied to solve the abovementioned problem. Using the BCM, an appropriate mesh resolution can be locally assigned over the entire computational domain. IBM helps to maintain a practical minimum mesh size and to achieve highly accurate wall boundary conditions.

The purpose of this research is to develop the LEEs solver to compute noise propagation around complex geometry easily and robustly. To achieve this purpose, the BCM is employed in conjunction with IBM. First, a one-dimensional wave propagation problem is evaluated. Second, acoustic scattering around a sphere is computed. Third, noise propagation from the JT15D nacelle is computed and compared with results obtained from other sources. At last, noise propagation around fuselage-wing-nacelle configurations is computed, and the noise shielding effect of the over-the-wing nacelle (OWN) configuration is estimated as a realistic computation around complex geometry.

#### 2. Computational Method

##### 2.1. Linearized Euler Equations

In the Euler equations, the physical quantities are divided into a mean flow field vector and fluctuation vector . By linearization of the Euler equations around , the time evolution equations of are obtained in (1). These are the LEEs. In the LEEs, the sound source term is added to incorporate a realistic sound source. For computation of the LEEs,* S* and the mean flow field vector are introduced. The time evolution of is then computed. The governing equations are nondimensionalized by the mean flow density , speed of sound* c*, and reference length* D*:Here, is the specific heat ratio. is the Kronecker delta. In the computations with mean flow field, , , and are neglected to avoid the instability waves caused by the shear mean flow field. The effect of neglecting the terms about the gradient of the mean flow field is discussed in [23].

##### 2.2. Computational Mesh of the BCM

The computational mesh of BCM is generated by the following procedures in two dimensions [24]. The computational domain is divided into an aggregation of square areas (cubic areas in three dimensions), where each area is denoted as a “cube,” as shown in Figure 1(a). Each cube is then divided by an equispaced Cartesian mesh, as shown in Figure 1(b). Cells located outside the wall boundary are defined as fluid cells. On the other hand, cells located inside the wall boundary are defined as wall cells. In this method, all cubes have the same number of cells so that the computational effort required for all cubes is essentially equivalent in parallel computation, and excellent parallel efficiency is achieved. Each cube has three overlap cells, as shown by the hatched cells in Figure 1(b) for data exchange. When a cube is locally refined, the cube selected for refinement is divided into four cubes (eight cubes in three dimensions), and each cube is subdivided by prescribed cells. After the refinement, the size of the cube is smoothed so that the sizes of adjacent cubes are restricted to the same, double, or half that size. The sizes of adjacent cubes are checked for all cubes. If the size of an adjacent cube is larger than double the size of a focusing cube, the adjacent cube is divided into four cubes (eight cubes in three dimensions), and each cube is subdivided by prescribed cells.

**(a) Computational domain and cube boundary**

**(b) Computational cells in a cube (15 × 15 cells, 3 overlap cells)**

##### 2.3. Computational Algorithm of the LEEs Solver

Figure 2 shows the computational algorithm of the LEEs solver. At the beginning, mesh information, object shape, initial conditions, and the mean flow field are given. The mean flow field is computed by the compressible Euler solver in this research. The details of the compressible Euler solver are provided in [25]. The computational scheme of inviscid flux is changed from that of the reference to the simple low-dissipative advection upstream splitting (SLAU) [26] scheme.

The spatial derivative of the LEEs solver is computed by the dispersion relation preserving (DRP) scheme [27]. Through the use of the overlap cells, a fourth-order DRP scheme using seven-point stencils can be implemented over the entire area of the computational domain. Time integration is computed by a six-stage fourth-order Runge-Kutta scheme [28]. In the time integral, six subiterations constitute one time step. In addition, artificial selective damping [27] of the seven-point stencils is applied at each iteration to eliminate the nonphysical oscillation. The above computations, including boundary treatments, are parallelized for all cubes using Open Multi-Processing (OpenMP).

##### 2.4. Wall Boundary

For computation by Cartesian mesh, it is important to apply mesh to wall boundaries realistically because real surfaces have curvature. Therefore, various IBMs have been proposed for diverse implementations. In the present solver, an IBM that employs a ghost cell (GC) approach using an image point (IP) [29] is applied. Wall cells adjacent to fluid cells are defined as GCs. The IP is defined from the GC in the direction normal to the closest wall boundary, as shown in Figure 3(a). In this process, surface stereo lithography (STL) data is used to determine the intersection point of the normal vector with the wall boundary. The physical quantities of the IP are interpolated by the inverse distance weighting method of (2) using 3 × 3 × 3 = 27 stencils, as illustrated in Figure 3(b):Here, is the physical quantities of the IP, isthe physical quantities of stencils, is the weight function based on the distance between the IP and stencils, and mask is a value that reflects whether the stencil is a fluid cell or wall cell. Fluid cells have , and wall cells have .

**(a) Definition of GC and IP**

**(b) Interpolation to IP**

The physical quantities of a GC are computed by the following equations so that the slip condition is satisfied using the physical quantities of the IP:Here, is the distance from the IP to the wall surface; is the distance from the GC to the wall surface; is the unit normal vector of the wall surface; and , , and are the velocity vector, pressure, and density at the IP and , , and are those at the GC, respectively. In this study, to avoid the instability caused by the wall boundary condition. The spatial derivative at the fluid cell adjacent to the GC is conducted using the second-order central difference based on the physical quantities of the GC.

##### 2.5. Cube Boundary

At the boundaries between cubes with different sizes, the three overlap cells are hanging nodes. Therefore, data exchange with interpolation is needed. In the solver, Lagrange interpolation [30] is employed for data exchange at the boundary, as given by Here, represents the physical quantities of an overlap cell; represents the physical quantities of stencils; , , and are computed weight functions based on the distance between an overlap cell and stencils; and , and are the coordinates of an overlap cell.

Figure 4(a) illustrates the stencils used for data exchange from smaller to larger cubes in two dimensions. In this figure, black points are the fluid cells of a larger cube. The shaded fluid cells of the larger cube are the overlap cells. The white circles are stencils of the smaller cube and are used for the interpolation. Interpolation to three overlap cells per one row is needed. The stencils of the overlap cell that is closest to the boundary are 2 × 2 × 2 = 8 points and those of the other two cells are 4 × 4 × 4 = 64 points so that the location of stencils is symmetrical. For interpolation from larger to smaller cubes, as shown in Figure 4(b), stencils of 3 × 3 × 3 = 27 points are used for four rows in three dimensions, which achieves third-order accuracy.

**(a) From a smaller to larger cube**

**(b) From a larger to smaller cube**

##### 2.6. Outer Boundary

In LEEs computation, outgoing waves should be damped so that the inner sound field is not disturbed by reflected waves. To this end, the buffer zone boundary condition [31] is implemented in the present solver. To damp the outgoing waves gradually, the magnitude of the damping coefficients varies as a quadratic function of the coordinates toward the outer boundary. Equations (5) provide the damping equation in the buffer zone. The fluctuation vector is damped in each direction individually:Here, is the fluctuation vector at each iteration and with respect to each direction; , , and are the damped fluctuation vectors; , and are the damping coefficients in the buffer zone; , and* Δz *are the mesh spaces; , , and are the widths of the buffer zone; and , and are the distances from the inner end of buffer zone. Figure 5 shows the buffer zone (illustrated as the hatched area) and an expanded view of the bottom left corner in a two-dimensional case.

**(a) Buffer zone in the computational domain**

**(b) Expanded view of the bottom left corner of (a)**

#### 3. Evaluation of the Order of Accuracy

The order of accuracy of the present computational method is evaluated by the one-dimensional wave propagation problem. The IBM and Lagrange interpolation are employed in their three-dimensional forms. Equations (6) are the one-dimensional wave equations. The initial pressure distribution is given byThe root mean square error (RMSE) between the analytical solution and the computed solution is computed byThe RMSE is computed for two meshes. Figure 6 illustrates the computational cubes with and without cube refinement. The size of the computational domain is . The periodic boundary condition is set at the end of the -directional boundaries. The RMSE is also computed in the case where a wall boundary exists. When the wall boundary exists, the wave reflects at the end of the -directional boundaries.

**(a) Cubes without cube refinement**

**(b) Cubes with cube refinement**

The RMSE and the order of accuracy are listed in Table 1. The order of accuracy of each mesh is computed using the RMSE of one mesh and that of a coarser mesh. Case 1 employs no wall boundary and is conducted without cube refinement. Therefore, the prescribed order of accuracy is as shown. In Case 2, the order of accuracy is slightly reduced because the spatial derivative at the fluid cell adjacent to the GC is conducted using the second-order central difference. The order of accuracy of Case 3 is approximately two because the accuracy of the cube boundary is dominated by the interpolation from a larger to smaller cube. Case 4 provides a slightly smaller order of accuracy than that of Case 3 because of the existence of the wall boundary. These results indicate that the general order of accuracy of the present computational method is approximately two.

#### 4. Acoustic Scattering around a Sphere

Acoustic scattering around a sphere [32] is computed under various mesh conditions in this section. The magnitude of error relative to an analytical solution and the computational time are compared. A sphere is located at the origin of a three-dimensional domain. The reference length is the diameter of the sphere. A monopole Gaussian sound source is given by the following equation as a function of time* t*:

The influence of three important parameters on the error is investigated: the minimum cell size around the sphere, points per wavelength (PPW) in the computational domain, and the size of the buffer zone. The minimum cell size around the sphere has an effect on the error generated from the wall boundary. A small cell size not only suppresses the error but also restricts the time step size and results in greater computational time. The PPW is an important factor for acoustic simulations because insufficient PPW introduces dissipation and dispersion errors in wave propagation. The size of the buffer domain is important for setting the outer boundary condition. If the width of the buffer zone is short, the reflected wave is generated at the outer boundary, and the reflected wave disturbs the inner sound field.

Figure 7 shows computational domains of the Coarse, BCM_Middle2, and BCM_Fine3 meshes listed in Table 2. The gray lines indicate the cube boundaries. Table 2 summarizes the mesh information of all the cases used for the parametric study. Coarse mesh forms the base mesh of the other meshes employed. The size of the computational domain is , referring to [33] where 128 cubes of a size constitute the domain, as shown in Figure 7(a). The number of cells in one cube is 48 × 48 × 48 and the PPW is eight. In the Uniform_Middle and Uniform_Fine meshes, the domain size and the number of cubes are the same as those of the Coarse mesh. The number of cells in one cube is 72 × 72 × 72 in the Uniform_Middle mesh and 96 × 96 × 96 in the Uniform_Fine mesh; thus, the PPW is 12 and 16, respectively. In the BCM_Middle1 mesh, the eight cubes in the Coarse mesh surrounding the sphere are subdivided into half size cubes. Similarly, the eight cubes surrounding the sphere are subdivided from the BCM_Middle1 mesh in the BCM_Fine1 mesh. The number of cubes is 184 in the BCM_Middle1 mesh and 240 in the BCM_Fine1 mesh, and the maximum PPW is 16 and 32, respectively. In the BCM_Middle2 mesh shown in Figure 7(b) and BCM_Fine2 meshes, the buffer zone is expanded to add 112 cubes to the BCM_Middle1 and BCM_Fine1 meshes. In the BCM_Middle3 and BCM_Fine3 (shown in Figure 7(c)) meshes, eight cubes located at 2 * x * 6, −2 * y * 2, and −2 * z * 2, and eight cubes located at −6 * x * −2, −2 * y * 2, and −2 * z * 2 are subdivided from the BCM_Middle2 and BCM_Fine2 meshes to increase the PPW of the computational domain. The inner end of the buffer zone is set at* x *= [−7*D*, 7*D*],* y *= [−3*D*, 3*D*], and* z *= [−3*D*, 3*D*] in all meshes to compare the sizes of the buffer zones.

**(a) Coarse**

**(b) BCM_Middle2**

**(c) BCM_Fine3**

Figure 8 shows the instantaneous pressure distribution around the sphere computed using the Coarse mesh. In this figure, the buffer zone is not drawn. Waves generated from the sound source are scattered by the sphere, and interference fringes appear. The computation is executed in nondimensional time , and a periodic steady state is observed after . Figure 9 shows the pressure distribution near the inner end of the buffer zone and the time history of the pressure near the sphere for all meshes. The pressure distribution is sampled at on the -axis. The time history of the pressure is sampled at (−2*D*, 0, 0) and .

**(a) Pressure distributions of uniform meshes**

**(b) Time histories of uniform meshes**

**(c) Pressure distributions of BCM_Middle meshes**

**(d) Time histories of BCM_Middle meshes**

**(e) Pressure distributions of BCM_Fine meshes**

**(f) Time histories of BCM_Fine meshes**

In Figures 9(a) and 9(b), the uniform meshes are compared. The pressure distributions do not follow the decreasing trend of the analytical solution in Figure 9(a). Reflection occurs at the outer boundary because of the insufficient width of the buffer zone. Additionally, the Uniform_Middle and Uniform_Fine meshes show a higher pressure amplitude compared with the pressure distribution of the Coarse mesh. It is assumed that the sound wave is damped at the Coarse mesh because of the insufficient PPW. The time histories of the uniform meshes provide lower amplitudes than the analytical solution.

In Figures 9(c) and 9(d), the BCM_Middle meshes are compared. The BCM_Middle2 and BCM_Middle3 meshes follow the analytical solution in Figure 9(c). These meshes have a buffer zone of width 9*D* in the direction, and outgoing waves are damped appropriately. However, the BCM_Middle2 mesh provides a lower amplitude. This is owing to the insufficient PPW at the sampling points. On the other hand, phase error is not observed in the BCM_Middle2 mesh. In Figure 9(d), the amplitude of the pressure converges to the analytical solution for the BCM_Middle2 mesh. Only the BCM_Middle1 mesh provides a lower amplitude. This amplitude is the same as that of the uniform meshes. Therefore, it is confirmed that the amplitude of the inner sound field is diminished when using a buffer zone of short width.

In Figures 9(e) and 9(f), the BCM_Fine meshes are compared. The pressure distributions and the time histories of the BCM_Fine meshes exhibit similar tendencies as those of the BCM_Middle meshes.

To compare the magnitude of errors more quantitatively, the root mean square pressure, , of the scattered sound is computed and the normalized RMSE (NRMSE) and the normalized maximum error (MME) between the analytical solution and the solutions of each mesh are compared, as given by

The number of sampling points are set clockwise at a radius *D *from the origin. The points (−2*D*, 0, 0) and (2*D*, 0, 0) correspond to 0 deg. and 180 deg. The errors are summarized in Table 3. As shown in the table, the refined meshes provide lower errors than the coarser mesh, and the meshes that utilize a large buffer zone also provide lower errors. The minimum cell size around the sphere also affects the errors. The smallest errors are provided by the BCM_Fine3 mesh. A comparison of the uniform and BCM meshes indicates that the BCM_Middle3 and BCM_Fine3 meshes, which use cubes of different sizes at the appropriate locations, provide lower errors than Uniform_Middle meshes that utilize a similar number of mesh points. Moreover, the errors of the BCM_Fine3 mesh are about one-half the errors of the Uniform_Fine mesh, which has twice the number of mesh points as the BCM_Fine3 mesh. These results indicate the effectiveness of the BCM arrangement.

The computational wall-clock times for a single time step per cell and the total wall-clock time required to reach a nondimensional time are also listed in Table 3. All meshes were computed using 128 CPUs of an SGI AltixUV1000 computer. Regarding the wall-clock time per one time step per cell, the computational times of the uniform and BCM meshes are comparable. However, the total wall-clock times of the BCM meshes are slightly longer than those of the uniform meshes. The BCM meshes employ a larger number of cubes than the number of CPUs used in the evaluation. The mismatch between the numbers of cubes and CPUs can degrade the load balance. Computations using an equivalent number of CPUs and cubes are more effective for BCM meshes.

#### 5. Noise Propagation from the JT15D Nacelle

Noise propagation from the JT15D [34] nacelle is computed and the far-field sound pressure level (SPL) is compared with data derived from experiment and the computational results of others. The JT15D is a small Pratt and Whitney turbofan engine with a bellmouth inlet. The sound source is introduced via a discrete frequency spinning mode. The spinning mode input [35] is a typical sound source generated by the fan or rotor-stator intersections in the engine. The sound source of a single spinning mode is defined byHere, and are the* m*th order Bessel functions of the first and second kind, respectively, is the angular frequency, is the axial wavenumber, is the radial wavenumber, and is the phase. The th solution of the following equation is determined by the hard-wall boundary condition of the duct:Here, and are the bypass duct inner wall radius and the inner hub radius, respectively. The axial wave number is computed from where is the Mach number at the fan face. The selection of plus or minus signs in the parentheses is determined by the propagation direction of the sound wave, where plus (+) represents the positive propagation direction along the axial coordinate, and vice versa. The constant satisfies the following relation:

Parameters are set from experimental data: , m, and = 0.0998 m. The reference length* D* = 0.6223 mis the length from the fan face to the leading edge of the nacelle. The blade passing frequency (BPF) is set to 3,150 Hz, which is derived from the fan rotating speed of 6,750 rpm. The experiments were conducted in the static condition, and the fan face axial Mach number is 0.175. The present computation is also conducted in the static condition. In the present computation, the axial flow velocity of the ghost cell at the fan face boundary is directly modified to this Mach number. The pressure and the density are the same values as those of the adjacent fluid cells. The SPL is measured at a radius of 49*D*. The integration method of Ffowcs-Williams and Hawkings (FW-H) [36] is employed to estimate the far-field pressure from the near-field pressure.

Figure 10 shows the computational domain and Table 4 lists the mesh information. The computations of the compressible Euler solver and the LEEs solver are conducted using the same mesh. The PPW in Table 4 is the PPW reduced by the Doppler effect of the fan face axial Mach number. The computational domain is set at , , and . The fan face is at . The inner end of the buffer zone is set at , , and . The FW-H surfaces are located at the cube boundaries between the smallest size cubes and the larger size cubes.

Figure 11 shows the Mach number distribution around the JT15D nacelle. In Figure 11, the acceleration region near the internal wall of the nacelle lip and the stagnation at the spike are shown. Because the experiments and computations were conducted under the static condition, the mean flow field exists only near the nacelle.

Figure 12 shows the instantaneous pressure distributions at the planes defined at and . In the figure, only the distributions of the smallest sized cubes are drawn. The fan noise is diffracted at the inlet of the nacelle and propagates along the radial direction. The 13 peaks and valleys of the swirling fan noise at the plane are clearly captured, which is equivalent to the number of spinning modes. The time history of the pressure is sampled at the established points on the FW-H surface, and the periodic steady state is observed after .

**(a)**plane

**(b)**planeFigure 13 shows a comparison of the predicted SPL with the experimental results conducted by Heidmann et al. [37] and the computational results computed by Lan et al. [34]. The horizontal axis in the figure is the angle from the +-axial direction. In the figure, the predicted SPL using the present solver shows good agreement with the computational result of Lan et al. (within 1 dB). Compared with the experimental data, a difference of about 3 dB is found at 60 deg. However, this discrepancy is also shown between the experimental data and the result of Lan et al. The SPL peak, obtained near 55 deg., is equivalent to that obtained by Lan et al. and that obtained experimentally.

#### 6. Noise Propagation around a Fuselage-Wing-Nacelle Configuration

The noise propagation around a fuselage-wing-nacelle configuration is computed as a realistic example. The focus of the example was that of the OWN configuration [38], where the engine nacelle is mounted over the main wing, which serves as one of the configurations designed to achieve substantial airport noise reduction. The main feature of this configuration is that the main wing serves as a shield between the ground and the noise generated by the engine. To investigate the effect of the OWN configuration, the propagation of fan noise of both the conventional DLR-F6 aircraft configuration [39] and the OWN configuration is computed and the SPL below the aircraft for both configurations is compared. These configurations have complex geometry with large curvatures near the wing-body and the nacelle-pylon junctions. In the OWN configuration, the nacelle is moved 1.1*D* in the + direction and 0.6*D* in the +*y* direction relative to its location in the DLR-F6 configuration. Here, the reference length m is the longitudinal length of the nacelle. The cross-sectional geometry of the pylon is used without any changes and has a sweepback angle of 40 deg. A mean flow field of Mach number 0.3 with zero angle of attack is computed using the compressible Euler solver. The sound source is introduced via a discrete frequency spinning mode. As an example of a modern turbofan engine, the CFM56-7B engine is considered. The bypass ratio is 5.5, the fan diameter is 1.54 m, the number of blades is 24, the maximum rotational rate is 5,382 rpm, and the fundamental BPF is 2,152.8 Hz. The BPF is used as the input frequency of the spinning mode. The spinning mode is set to (−24, 0). The SPL is measured at a radius of 10*D* using the FW-H method.

Figure 14 shows the computational domain of each configuration. The black dot shows the origin of coordinate system. The computations of the compressible Euler and LEEs solvers are conducted using an equivalent mesh. The size of the computational domains is . The symmetric boundary condition is set at the −*z* boundary. The computational domains are set to ensure that sufficient computational domains surround the nacelle. It is assumed that the geometry out of the computational domain has little influence on the noise directivity below the aircraft. The minimum size cubes are located within the domain wherein the reflection and diffraction by the aircraft must be resolved. The inner end of the buffer zone is set at , , and in the DLR-F6 calculation and at , , and in the OWN calculation. FW-H surfaces are located at the inner end of the buffer zone and the symmetrical plane. The size of the buffer zones of the DLR-F6 calculation is 0.75*D *in the and directions, 0.5*D* in the −*y* and* +y* directions, and 0.25*D* in the* +z* direction. The size of the buffer zones of the OWN calculation is 0.5*D *in the −*x *direction, 0.75*D *in the* +x* direction, 0.25*D* in the −*y* direction*, *0.75*D *in* +y* direction, and 0.25*D *in the* +z* direction. Table 5 lists the mesh information of the computations. The PPW in Table 5 is the PPW reduced by the Doppler effect of the mean flow field Mach number.

**(a) DLR-F6 configuration**

**(b) OWN configuration**

Figure 15 shows the Mach number distributions at the nacelle center for the two configurations. In Figure 15, the acceleration region over the wing and the stagnation at the inlet of the nacelle are shown. A Mach number greater than 0.3 is shown in the acceleration region. Figure 16 shows the SPL distributions on the aircraft surface for the two configurations. Noise from the nacelle propagates in the radial direction and the noise is shielded by the main wing in the OWN configuration. On the other hand, noise from the inlet reaches the fuselage and generates a high SPL in the OWN configuration compared to the DLR-F6 configuration.

**(a) DLR-F6 configuration**

**(b) OWN configuration**

**(a) DLR-F6 configuration**

**(b) OWN configuration**

To check the propriety of the computations, the SPL distribution computed for the OWN configuration on the suction side of the main wing is compared with data captured during flight and the computations of others [40]. The flight data of a twin-engine business jet aircraft was acquired at a flight Mach number of 0.3 at an altitude of 1,500 m. This aircraft has the aft fuselage nacelle configuration that the inlet of the engine nacelle is located over the main wing. Therefore, qualitative comparisons can be conducted. The SPL on the suction side of the main wing was measured using Kulite microphones. The computation of Xu et al. was conducted using the discrete frequency spinning mode (20, 0) without a mean flow field. The computed SPL is normalized so that the peak SPL is equal to that of the flight data. Figure 17 shows the SPL versus the angle from the inlet axis. The computed SPL agrees with the flight data from 50 to 70 deg. within 3 dB. The peak SPL is at 63 deg. This is similar to the peak given by the flight data, which has a peak at 60 deg. The discrepancy at an angle of 40 deg. is caused by the different clearances between the engine and the main wing, where the clearance of the aircraft employed in the experiment is smaller than that of the present computation.

Figure 18 shows the estimated SPL distribution at a radius of 10*D* based on the center of the front face of nacelle. Based on International Civil Aviation Organization rules, aircraft noise is determined by the noise levels at the side and the bottom locations relative to the fuselage. Therefore, the SPL distribution below the aircraft is estimated. From Figure 18, the SPL of the OWN configuration is lower than that of the DLR-F6 configuration by about 10 dB over the range of 40 deg. to 70 deg. This represents significant noise reduction relative to conventional aircraft design. Using the proposed method, it would be possible to optimize the position of the nacelle to obtain maximum noise shielding performance.

#### 7. Conclusion

The linearized Euler equationssolver on block-structured Cartesian mesh of the building-cube method** (**BCM) coupled with the immersed boundary method (IBM) has been developed to compute noise propagation around complex geometry. The accuracy and effectiveness of the solver for practical problems were validated by application to one-dimensional and three-dimensional noise propagation problems and the application of fan noise propagation around fuselage-wing-nacelle configurations.

For the one-dimensional noise propagation problem, the IBM slightly decreased the order of accuracy of the system because the second-order central difference is employed at fluid cells adjacent to ghost cells. However, a fourth order of accuracy is nearly retained even when employing the IBM. Interpolation between cubes of different sizes causes degradation to a second order of accuracy.

For the problem of acoustic scattering around a sphere, the error was suppressed by mesh refinement near the sphere and by expansion of the buffer zone. The results computed on the BCM meshes provided quantitative agreement with the analytical solution. A normalized root mean square error of 1.62% and a normalized maximum error of 5% for the results computed by the BCM_Fine3 mesh were about one-half of the errors computed by the Uniform_Fine mesh, which has twice as many mesh points as the BCM_Fine3 mesh. The computational time required for the BCM mesh was slightly longer than that required for the uniform mesh because of the mismatch between the number of CPUs and cubes.

The estimated sound pressure level (SPL) from the JT15D nacelle provided good agreement with experimental data and with other computational results. The present solver demonstrated accurate computations for a curved object like an axisymmetric nacelle.

Noise propagation around two fuselage-wing-nacelle configurations was computed as a realistic example to verify the capability of the present solver and to estimate the noise shielding effect of the OWN configuration. The aircraft configuration has complicated geometries, especially near the wing-body and the nacelle-pylon junctions, and the present solver demonstrated a robust treatment of these geometries. For the OWN configuration, the noise from the nacelle inlet was shielded effectively by the main wing because the noise diffracted strongly in the radial direction. The computed SPL at the suction side of the main wing agrees with data obtained during flight. The SPL distribution of the OWN configuration below the aircraft was lower by about 10 dB than that of the conventional DLR-F6 configuration.

Through application to simple test cases and realistic cases, the effectiveness of the solver was confirmed. Using the present solver, noise propagation around next-generation unconventional aircraft can be estimated, and the method is expected to contribute to the future of aircraft design.

#### Nomenclature

: | Cell index coordinates |

: | Coordinates in the computational domain |

: | Mesh spacing |

: | Physical quantities vector |

: | Fluctuation vector |

: | Velocity fluctuation of , , and direction |

: | Mean flow field vector |

: | Mean velocity of , and direction |

: | Sound source vector |

: | Specific heat ratio |

: | Kronecker delta |

: | Unit normal vector |

: | Velocity vector at the image point |

: | Velocity vector at the ghost cell |

: | Distance from the image point to the wall surface |

: | Distance from the ghost cell to the wall surface |

: | Physical quantities at an overlap cell |

: | Physical quantities at a surrounding cell |

: | coordinates of an overlap cell |

: | Damping coefficient in the buffer zone |

: | Distance from the inner end of the buffer zone |

: | Width of the buffer zone |

: | Reference length |

: | Initial pressure |

: | Computed pressure |

: | Analytical solution of the pressure |

: | Nondimensional time |

: | Root mean square pressure |

: | Speed of sound |

: | Input spinning mode |

: | th order Bessel functions of the first and second kind. |

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This work was supported by JSPS KAKENHI Grant nos. 21226018 and 259225. The computation in this research was conducted using an SGI AltixUV1000 of the Institute of Fluid Science, Tohoku University. The Ffowcs-Williams and Hawkings code used in this research was provided by the Aviation Program Group of Japan Aerospace Exploration Agency (JAXA).