Open Access

ARTICLE

# Numerical Simulation of Aerodynamic Interaction Effects in Coaxial Compound Helicopters

School of Aeronautic Science and Engineering, Beihang University, Beijing, 100191, China

* Corresponding Author: Yihua Cao. Email:

*Fluid Dynamics & Materials Processing* **2023**, *19*(5), 1301-1315. https://doi.org/10.32604/fdmp.2023.023435

**Received** 26 April 2022; **Accepted** 22 August 2022; **Issue published** 30 November 2022

## Abstract

The so-called coaxial compound helicopter features two rigid coaxial rotors, and possesses high-speed capabilities. Nevertheless, the small separation of the coaxial rotors causes severe aerodynamic interactions, which require careful analysis. In the present work, the aerodynamic interaction between the various helicopter components is investigated by means of a numerical method considering both hover and forward flight conditions. While a sliding mesh method is used to deal with the rotating coaxial rotors, the Reynolds-Averaged Navier-Stokes (RANS) equations are solved for the flow field. The Caradonna & Tung (CT) rotor and Harrington-2 coaxial rotor are considered to validate the numerical method. The results show that the aerodynamic interaction of the two rigid coaxial rotors significantly influences hover’s induced velocity and pressure distribution. In addition, the average thrust of an isolated coaxial rotor is smaller than that of the corresponding isolated single rotor. Compared with the isolated coaxial rotor, the existence of the fuselage results in an increment in the thrust of the rotors. Furthermore, these interactions between the components of the considered coaxial compound helicopter decay with an increase in the advance ratio.## Keywords

Nomenclature

| Chord of the rotor blade section |

| Fuselage lift |

| Air pressure |

| Freestream pressure |

| Rotor torque |

| Radial distance to the blade section |

| Rotor radius |

| Rotor thrust |

| Forward flight speed |

| Air density |

| Rotor angular velocity |

| Advance ratio |

| Lift coefficient of the fuselage |

| Pressure coefficient of the blade section |

| Rotor thrust coefficient |

| Rotor torque coefficient |

Although conventional helicopters possess excellent hover-ability and low-speed maneuverability, the maximum forward speed of conventional helicopters is restricted by the adverse effects of the compressibility on the advancing blades and stalling on the retreating blades of the main rotor [1]. In recent decades, a coaxial configuration with a tail-mounted propeller, referred to as the coaxial compound helicopter, has been developed to increase the maximum forward speed. The coaxial compound helicopter features two rigid rotors located on one axis and rotating in opposite directions. The rigid rotor originates from the Advancing Blade Concept (ABC) with the increased lift ability of the advancing blades [2,3]. However, the coaxial compound helicopter suffers severe aerodynamic interactions because of the small separation between the two rigid coaxial rotors. In addition, the downwash flow of the main rotor is blocked by the fuselage in hover and low-speed forward flights, which complicates the aerodynamic characteristics of the whole helicopter.

The methods to investigate the aerodynamic characteristics of coaxial configurations include experimental tests, vortex methods, and computational fluid dynamics (CFD) simulations. Experimental tests focus on the static-thrust performance in hover [4], the impact of axial rotor spacing on the performance characteristics [5,6], and wake structures [7,8]. The vortex methods, including the free-wake method [9,10], the vorticity transport model [11], and the vortex particle method [12], are utilized to analyze the aerodynamic interference effects and unsteady aerodynamic loads. However, the experimental tests are high-costs and the simplifications limit the fidelity of the vortex methods.

CFD simulations have been widely employed to investigate the aerodynamic characteristics of coaxial configurations with acceptable costs and high fidelity. Based on the structured overset-mesh method and Reynolds-Averaged Navier-Stokes (RANS) solver, Lakshminarayan et al. [13] investigated the hover aerodynamic characteristics of a microscale coaxial configuration varying from a range of rpm and rotor spacing. The overset mesh system contained a stationary background mesh region and several rotating near-body mesh regions. A hole-cutting technique was used to find the connectivity of the two parts. Based on the unstructured overset-mesh method, Xu et al. [14,15] investigated the aerodynamics of a coaxial configuration with a Robin fuselage in hover and low-speed forward flight by solving the unsteady Euler equations. Compared with the RANS method, the Euler method did not include viscous effects. Zhao et al. [16] utilized the sliding mesh method and RANS method to predict the aerodynamic performance of the rigid coaxial rotor in hover and high-speed forward flight. The interface boundary was employed to connect the stationary mesh and rotating mesh in the sliding mesh system. Deng et al. [17] analyzed the rotor lift-drag ratio and lateral lift-offset of the rigid coaxial rotor through numerical simulations and experiments. The lateral lift-offset is a unique feature of the rigid coaxial rotor. Based on mixed meshes and solvers, Park et al. [18] investigated the influence of the coaxial rotor spacing on the aerodynamic interactions in hover and forward flights. The inviscid fluxes in the off-body region were calculated with the 7th order Weighted Essentially Non-Oscillatory (WENO) scheme. Utilizing high-fidelity CFD and CSD (Computational Structural Dynamics) loose coupling approach, Jia et al. [19] investigated the impact of the pitch attitude on the aerodynamic interactions and acoustics of a rigid coaxial rotor in high-speed forward flight. In that study, the influence of rotor hubs was also investigated. The presence of rotor hubs results in hub-wake interactions, while the magnitude of the root-induced blade-vortex interaction is reduced in the presence of rotor hubs.

Most of the literature focuses on the rigid coaxial rotor, and the fuselage is not involved. In this paper, a numerical simulation for a compound helicopter is constructed based on the sliding mesh method and the RANS method. The cases of the isolated single rotor, isolated coaxial rotor, and whole helicopter are calculated. Then, the aerodynamic interactions between the components of the compound helicopter are analyzed. The simulated states include hover and forward flights.

2.1 Mathematic & Numerical Methods

Although the Navier-Stokes (NS) equations are applicable for all turbulence flow, the direct numerical simulation of the turbulence flow is prohibitively expensive. The RANS method provides a way to reduce the cost, focusing on the time-averaged flow and the effects of turbulence on time-averaged flow properties. The RANS equations contain the continuity equation, the momentum equations, and the scalar transport equation, as follows [20]:

where

where

In this paper, the compressible RANS equations are solved by the Finite Volume Method (FVM) [22]. The pressure-velocity coupling is dealt with the Coupled scheme, and the spatial discretization method is the second-up wind scheme. The turbulence model used is the Shear-Stress Transport (SST)

where

The sliding mesh method is a special case of general dynamic mesh motion wherein the nodes move rigidly in a given dynamic mesh zone. There are at least two cell zones that connect through non-conformal interfaces. This means that one boundary face has at least two neighboring cells. In practice, the two interfaces of the neighboring zones are replaced with a new set of faces formed by their intersections [24]. For the two-dimension mesh, the interfaces are shown in Fig. 1. The A-B, B-C, D-E, and E-F are the interfaces, and the a-b, b-c, c-d, d-e, and e-f are the intersections. The utilization of the intersection avoids interpolation among these neighbors. For instance, the flux crossing the D-E interface is the sum of the flux crossing the b-c and c-d intersections, and the flux of the b-c and c-d intersections is just from the A-B interface and B-C interface, respectively.

A typical coaxial compound helicopter is the X2, which utilizes a similar rotor system to the XH-59A [25]. However, the detailed geometry of X2 is difficult to find in the existing literature. Therefore, a simplified coaxial compound helicopter, combined with the X2 and XH-59A, is constructed in this paper. This coaxial compound helicopter comprises two rigid coaxial rotors and a fuselage with a tailplane and a vertical fin, as shown in Fig. 2. It should be noted that the rotor hub and tail-mounted propeller are ignored. The geometric parameters of this coaxial compound helicopter are shown in Tables 1 and 2 [26].

Three unstructured poly-hex-core mesh regions, including two rotating cylindrical regions and a stationary region, are generated for this numerical simulation. The poly-hex-core mesh contains the polyhedral mesh and the hex-core mesh. The polyhedral mesh is used in the field close to the boundaries, while the hex-core mesh is used inside the core of the domain. The poly-hex-core mesh has a lower mesh count than the equivalent tetrahedral mesh, and it can handle complex geometries.

The rotating cylindrical mesh regions cover the upper and lower rotor blades. Mesh refinement is performed on the blades’ leading edge and trailing edge. The mesh size of the rotating region interface is set to march the length when the rotor rotates

To speed up the convergence and reduce computation time, the steady-state simulation of the flow field is firstly calculated by the moving reference frames method, and then the steady-state results are used as the initial values of transient analysis of the sliding mesh method. The time step of transient analysis is set as the time when the main rotor rotates

The Caradonna & Tung (CT) rotor with experimental data [27] is utilized to verify the validity of the constructed numerical method. This rotor has two rectangular blades without torsion, and the blades use a NACA0012 profile. The diameter is 2.286 m, and the cord is 0.1905 m. The hover state of

The Harrington 2 rotor is a coaxial rotor from the reference [4], which is widely used to verify the validity of numerical simulations for coaxial rotors. This coaxial rotor has two rectangular blades for each rotor, and the thickness ratio of the blades is tapered. The diameter is 7.62 m, and the cord is 0.4572 m. Then, the cases of 0.293 Mach number and

In this section, cases of the isolated single rotor, isolated coaxial rotor, and the whole helicopter are calculated, and the aerodynamic interactions between these components are analyzed in hover and forward flights. The simulation conditions are 500 m height, 0.566 Mach number of the blade tip, and

The hover state of the isolated single rotor is calculated with three kinds of grid sizes. Each grid is divided into two regions, the rotating region, and the stationary region. The grids of both regions are refined. The average thrust coefficient (

4.2 Aerodynamic Interactions in Hover

4.2.1 Rotor-Rotor Aerodynamic Interactions

Fig. 6 presents the Z velocity contours of the isolated single rotor and the isolated coaxial rotor in hover. It can be seen that the induced velocity of the isolated coaxial rotor is larger than that of the isolated single rotor. In addition, the induced velocity of the lower rotor is affected by the upper rotor, which reduces the effective angle of attack. Then, the pressure distributions of the blade section changed. Fig. 7 presents the pressure distributions of the upper and lower rotors when the upper and lower blades encounter each other. The negative peak pressure of the lower rotor is smaller than that of the upper rotor, and the discrepancy increases at the outboard of the blade. This demonstrates the velocity distribution above.

Fig. 8 presents the variations of transient thrusts for the isolated single rotor and the isolated coaxial rotor in one revolution. The

4.2.2 Rotor-Fuselage Aerodynamic Interactions

Fig. 9 presents the Z velocity distributions of the plane 0.4 m below the lower rotor disc at the azimuth of

Fig. 10 compares the transient thrusts between the isolated coaxial rotor and the coaxial rotor in the whole helicopter. The thrusts of the upper and lower rotors in the whole helicopter are larger than that of the isolated coaxial rotor. This indicates that the existence of the fuselage has a positive effect on the rotor lift capability.

4.3 Aerodynamic Interactions in Forward Flight

4.3.1 Rotor-Rotor Aerodynamic Interactions

The cases of the isolated single rotor, the isolated coaxial rotor, and the whole helicopter are calculated at advance ratios of 0.1, 0.2, and 0.3.

Fig. 11 presents the transient thrusts of the rotors at these advance ratios. It can be seen that the thrust of the coaxial rotor is still smaller than that of the single rotor due to the interactions between the coaxial rotors. Moreover, the periodic thrust variation of the isolated single rotor remains three-per-revolution. However, the addition of the forward speed changes the periodic thrust variation of the coaxial rotors. Compared with the hover state, the periodic thrust variation of the upper and lower rotors changes to three-per-revolution in forward flight. Furthermore, the thrust of the upper rotor oscillates at the azimuths of

To identify the interactions of the coaxial rotors in forward flights, the pressure distributions of the isolated single rotor and coaxial rotors at the advancing and retreating sides are compared, as shown in Fig. 12. The disc azimuth here is

4.3.2 Rotor-Fuselage Aerodynamic Interactions

Fig. 13 presents the velocity contours in the longitudinal cross-section at different advance ratios. At the advance ratio of 0.1, the wake of the main rotors is blocked by the fuselage on a large scale, which causes severe aerodynamic interactions. As the advance ratio increases, the wake of the main rotor tilts in the backward direction. The impact of the main rotor downwash on the fuselage is weakened at the advance ratio of 0.2. Moreover, the downwash of the main rotor has little influence on the fuselage at the advance ratio of 0.3. This phenomenon is consistent with the variation of fuselage lift, as shown in Fig. 14. The fuselage lift coefficient is negative because of the rotor downwash. However, the value of the fuselage lift coefficient decreases with the increment of the advance ratio. These indicate that the interactions between the main rotor and the fuselage decay with the increase of the advance ratio.

Table 4 illustrates the fuselage’s effect on the coaxial rotor’s average thrusts in one revolution. Because the downwash of the main rotor is blocked by the fuselage, the average thrust of the main rotor increases. However, this effect is weakened when the advance ratio increases, which is in accordance with Fig. 13.

Based on the sliding mesh method, a numerical method for coaxial rotor helicopter aerodynamic interaction simulation is constructed. The rotor-rotor interactions and rotor-fuselage interactions are analyzed in hover and forward flights. The following conclusions are obtained:

1. Compared with the isolated single rotor, the aerodynamic interaction of the two rigid coaxial rotors significantly influences the induced velocity and pressure distribution in hover. This interaction also reduces the thrust of the isolated coaxial rotor.

2. Compared with the hover state, the periodic thrust variation of the upper rotor and lower rotor changes from six-per-revolution to three-per-revolution in forward flights. The encounter of the upper and lower rotors results in the oscillation of the transient thrusts.

3. The thrust of the coaxial rotor increases due to the fuselage’s blocked downwash of the coaxial rotor. This indicated that the fuselage positively affects the lift capability of the coaxial rotor.

4. The interactions between the components of the coaxial helicopter decay with the advance ratio increasing.

Funding Statement: This work was supported by Rotor Aerodynamics Key Laboratory [Grant No. RAL202102-4].

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.

## References

- Cheney, M. C. (1969). The ABC helicopter. VTOL Research, Design, and Operations Meeting, Geogia.
- Ruddell, A. J. (1981). Advancing blade concept (ABC) development test program. 1st Flight Test ConferenceNevada.
- Ruddell, A. J. (1981). Advancing blade concept (ABC) technology demonstrator. NASA-ADA-100181.
- Harrington, D. (1951). Full-scale-tunnel investigation of the static-thrust performance of a coaxial helicopter rotor. NACA-TN-2318.
- Landgrebe, A. J., Bellinger, E. D. (1974). Experimental investigation of model variable-geometry and ogee tip rotors. NASA-CR-2275.
- Ramasamy, M. (2013). Measurements comparing hover performance of single, coaxial, tandem, and tilt-rotor configurations. AHS 69th Annual Forum, Phoenix, AZ.
- Konus, M. F. (2017). Vortex wake of coaxial rotors in hover (Ph.D. Thesis). Mechanical Engineering, University of California, Berkeley, USA.
- Silwal, L., Raghav, V. (2020). Preliminary study of the near wake vortex interactions of a coaxial rotor in hover. AIAA Scitech 2020 Forum.
- Lyu, W., & Xu, G. (2018). Interactional effect of propulsive propeller location on counter-rotating coaxial main rotor.
*Journal of Aircraft*,*55(6)*, 2538-2545. [Google Scholar] [CrossRef] - Syal, M., & Leishman, J. G. (2012). Aerodynamic optimization study of a coaxial rotor in hovering flight.
*Journal of the American Helicopter Society*,*57(4)*, 1-15. [Google Scholar] [CrossRef] - Kim, H. W., & Brown, R. E. (2010). A rational approach to comparing the performance of coaxial and conventional rotors.
*Journal of the American Helicopter Society*,*55(1)*, 012003. [Google Scholar] [CrossRef] - Tan, J., Sun, Y., & Barakos, G. N. (2018). Unsteady loads for coaxial rotors in forward flight computed using a vortex particle method.
*The Aeronautical Journal*,*122(1251)*, 693-714. [Google Scholar] [CrossRef] - Lakshminarayan, V. K., & Baeder, J. D. (2010). Computational investigation of microscale coaxial-rotor aerodynamics in hover.
*Journal of Aircraft*,*47(3)*, 940-955. [Google Scholar] [CrossRef] - Xu, H., & Ye, Z. (2010). Coaxial rotor helicopter in hover based on unstructured dynamic overset grids.
*Journal of Aircraft*,*47(5)*, 1820-1824. [Google Scholar] [CrossRef] - Xu, H., & Ye, Z. (2011). Numerical simulation of unsteady flow around forward flight helicopter with coaxial rotors.
*Chinese Journal of Aeronautics*,*24(1)*, 1-7. [Google Scholar] [CrossRef] - Zhao, X., Zhou, P., Yan, X., Lin, Y. (2017). Numerical study of aerodynamic performance and flow interaction for coaxial rigid rotor in hover and forward flight. 35th AIAA Applied Aerodynamics ConferenceColorado.
- Deng, J. H., Fan, F., Liu, P. A., Huang, S. L., & Lin, Y. F. (2019). Aerodynamic characteristics of rigid coaxial rotor by wind tunnel test and numerical calculation.
*Chinese Journal of Aeronautics*,*32(3)*, 568-576. [Google Scholar] [CrossRef] - Park, S. H., & Kwon, O. J. (2020). Numerical study about aerodynamic interaction for coaxial rotor blades.
*International Journal of Aeronautical and Space Sciences*,*22(2)*, 277-286. [Google Scholar] [CrossRef] - Jia, Z., & Lee, S. (2021). Aerodynamically induced noise of a lift-offset coaxial rotor with pitch attitude in high-speed forward flight.
*Journal of Sound and Vibration*,*491*, 1-34. [Google Scholar] [CrossRef] - Versteeg, H. K., Malalasekera, W. (2007). An introduction to computational fluid dynamics: The finite volume method. England: Pearson Education.
- Hinze, J. O. (1975). Turbulence. New York: McGraw-Hill.
- Patankar, S. V. (1980). Numerical heat transfer and fluid flow. New York: Hemisphere Publishing Corporation.
- Menter, F. R. (1994). Two-equation eddy-viscosity turbulence models for engineering applications.
*AIAA Journal*,*32(8)*, 1598-1605. [Google Scholar] [CrossRef] - Mathur, S. (1994). Unsteady flow simulations using unstructured sliding meshes. Fluid Dynamics Conference, Colorado Springs, CO.
- Walsh, D., Weiner, S., Arifian, K., Bagai, A., Lawrence, T. et al. (2009). Development testing of the sikorsky X2 technologyTM demonstrator. The 65th Annual Forum of the American Helicopter Society InternationalGrapevine, TX.
- Kim, H. W., Kenyon, A. R., Duraisamy, K., Brown, R. E. (2008). Interactional aerodynamics and acoustics of hingeless coaxial helicopter with an auxiliary propeller in forward flight. 9th International Powered Lift ConferenceLondon.
- Caradonna, F. X., Tung, C. (1981). Experimental and analytical studies of a model helicopter rotor in hover. NASA-TM-81232.

## Cite This Article

Wang, M., Wang, Y., Cao, Y., Zhang, Q. (2023). Numerical Simulation of Aerodynamic Interaction Effects in Coaxial Compound Helicopters.*FDMP-Fluid Dynamics & Materials Processing, 19(5)*, 1301–1315.