Helicopter Rotor Flow Analysis Using Mapped Chebyshev Pseudospectral Method and Overset Mesh Topology

Unsteady helicopter rotor flows are solved by a Chebyshev pseudospectral method with overset mesh topology which employs Chebyshev polynomials for solution approximation and a Chebyshev collocation operator to represent the time derivative term of the unsteady flow governing equations. Spatial derivative terms of the flux Jacobians are discretized implicitly while the Chebyshev spectral derivative term is treated in explicit form. Unlike the Fourier spectral method, collocation points of standard Chebyshev polynomials are not evenly distributed and heavily clustered near the extremities of the time interval, which makes the spectral derivative matrix ill-conditioned and deteriorates the stability and convergence of the flow solution. A conformal mapping of an arcsin function is applied to redistribute those points more evenly and thus to improve the numerical stability of the linear system. A parameter study on the condition number of the spectral derivative matrix with respect to the control parameters of the mapping function is also carried out. For the validation of the proposedmethod, both periodic and nonperiodic unsteady flow problemswere solved with two-dimensional problems: an oscillating airfoil with a fixed frequency and a plunging airfoil with constant plunging speed without considering gravitational force. Computation results of the Chebyshev pseudospectral method showed excellent agreements with those of the time-marching computation. Subsequently, helicopter rotor flows in hovering and nonlifting forward flight are solved. Moving boundaries of the rotating rotor blades are efficiently managed by the overset mesh topology. As a set of subgrids are constructed only one time at the beginning of the solution procedure corresponding to the mapped Chebyshev collocation points, computation time for mesh interpolation of hole-cutting between background and near-body grids becomes drastically reduced when compared to the time-marching computation method where subgrid movement and the hole-cutting need to be carried out at each physical time step. The number of the collocation points was varied to investigate the sensitivity of the solution accuracy, computation time, and memory. Computation results are compared with those from the time-marching computation, the Fourier spectral method, and wind-tunnel experimental data. Solution accuracy and computational efficiency are concluded to be great with the Chebyshev pseudospectral method. Further applications to unsteady nonperiodic problems will be left for future work.


Introduction
Challenges in helicopter rotor flow analysis are several.Blade wake interactions need to be solved accurately requiring higher mesh resolution in the wake-trailing regions.Moving boundaries of the rotating blades need to be handled dynamically in time [1][2][3][4][5].A typical local, finite difference method for the time integration of the flow governing equations of the URANS (unsteady Reynolds-averaged Navier-Stokes equations), such as BDF2 (2-step backward differentiation formula) or BDF3 (3-step backward differentiation formula), of the dual-time-stepping method [6,7] requires a small time step for solution stability.Although the pseudo-time integration for the inner iteration and a multigrid technique can make the solution converge faster, the time-marching computation is still expensive even with highly parallel computation as it needs to fully resolve unsteady transient flows before flows reach a periodic steady-state.On the other hand, a spectral method for time integration has been proposed [8][9][10][11][12][13][14][15][16] as an alternative to a time-marching, dual-time stepping method with much reduced computational cost at equivalent solution accuracy.With the flow solution approximation by a discrete Fourier series, a time derivative term of the governing equations reduces to a spectral derivative term and removes time dependence of the flow solutions.The timespectral derivative is the multiplication of the spectral derivative matrix and solution vectors at all time instances.The flow solutions advance in a pseudo-time domain with periodic steady-state remaining throughout the iterations until flows converge.Variants of the time-spectral method have been developed by various researchers.Hall et al. [8] originally developed a harmonic balance method (HBM) to solve internal flows of a two-dimensional compressor and use a discrete Fourier series for the solution approximation.The improved HBM by Thomas et al. [9,10] also approximates flow variables by a discrete Fourier series and represents the time derivative term of the governing equations in Fourier coefficients but solves the FFT-transformed governing equations directly in the time domain by taking the inverse transform of the frequency domain governing equations.This approach is more convenient as the flow solutions are integrated directly in time with residuals checked simultaneously in the time domain.Based on the same principle, McMullen et al. [13,14] developed a nonlinear frequency domain (NLFD) method which is slightly different from the HBM in the details of the implementation.Although the flow approximation by the Fourier series and the formulation of the governing equations in the frequency domain are the same as the HBM, the flow solutions are directly solved in the frequency domain, while the convergence is checked by the residuals reformulated in the time domain by taking the IFFT of the conservative flow variables.Gopinath and Jameson [12] proposed a timespectral method which replaces the time derivative term of the governing equations through the spectral derivative matrix operator using Fourier spectral differentiation.Similar to the HBM developed by Thomas et al. [9,10], it integrates the spectral derivative term of the governing equations directly in the time domain.Ekici et al. [11] solved helicopter rotor flows using the harmonic balance method with periodic boundary conditions both in time and in space, which makes the computation domain more compact than that of the time-spectral method and contributes to further decreasing computation time and memory in comparison to the time-marching computation.Choi et al. [17,18] used the time-spectral CFD computation method for fluidstructure interaction analysis of the rotor flows in unsteady forward flight of the UH60-A helicopter.Adjoint sensitivity analysis has been carried out in the time-spectral form of the URANS equations to compute the sensitivity of aerodynamic performance metrics (drag, lift, torque, etc.) with respect to hundreds of shape design parameters.Despite the rigid blade assumption, aerodynamic shape optimization using timespectral and adjoint sensitivity analysis was validated with the FSI analysis afterward with good performance and shown to be effective in solving unsteady (i.e., periodic steady) flows.A parameter study on the accuracy and computation efficiency of flow solutions with respect to varying number of harmonics has been carried out at various flight conditions.A Chebyshev pseudospectral method [19][20][21][22] uses Chebyshev polynomials for the solution approximation of the boundaryvalue problems, either periodic or nonperiodic, and the spectral derivative matrix is derived correspondingly to continuous differentiation.Standard Chebyshev polynomials are defined in the interval of [−1, 1] and are a member of an orthogonal polynomial set.They can approximate function values with great precision and are applied practically to mathematics, physics, and engineering areas.Dinu et al. [19] applied the Chebyshev pseudospectral method to an aeroelastic problem and solved unsteady aerodynamic forces accurately.Niu et al. [20] solved unsteady cylinder flows with the Chebyshev pseudospectral method applied to the spatial derivative terms of the momentum and energy conservation equations.Flow solutions demonstrated great agreements with those of the time-marching computation which uses the multistage Runge-Kutta method.Khater et al. [21] applied the method to the spatial derivative term of the flow governing equations of the Burgers equations in various forms and excellent agreements were observed in comparison to the exact solutions.Im et al. [23] applied a Chebyshev collocation operator to the time derivative term of the URANS equations and employed a conformal mapping function to improve the numerical stability of the method.They applied the method to the oscillating airfoil with a prespecified frequency and compared the computation results with those of the dual-time stepping method and the harmonic balance method.In the current study, helicopter rotor flows in a hovering condition and forward level-flight are solved using the mapped Chebyshev pseudospectral method applied to the time integration of the URANS equations.To accurately solve the blade-vortex interaction and vortex wake flows trailing from the rotor blades, two key approaches are made: conformal mapping of an arcsin function and the overset mesh topology.The arcsin function is used for the conformal mapping of the uneven collocation points of the standard Chebyshev polynomials and yields points more evenly distributed over the physical time interval of [ initial ,  final ].Improved numerical stability is shown with respect to varying control parameters of the function.The application of the overset mesh topology is particularly advantageous with the mapped Chebyshev pseudospectral method, as the overset mesh construction and a hole-cutting procedure is carried out once and for all for subtime intervals corresponding to mapped collocation points of the Chebyshev polynomials.Cost for identifying holecutting topology between the near-body and background grids is reduced considerably when compared to the timemarching computation where it is carried out at each physical time step with the size often limited to a small value due to the numerical stability.The ratio of the number of subtime intervals between the time-spectral and the finite difference method is about (10) : (10 2 ), and the corresponding cost saving is more than an order of magnitude.Despite some disadvantage of increase in computation memory per solution iteration due to the storage of meshes and flow solutions at all subtime intervals, the mapped Chebyshev pseudospectral method has shown cost effectiveness in the present study.In the current work, the overset mesh topology will be used for the mapped Chebyshev pseudospectral method to solve helicopter flows of Caradonna and Tung's rotors in both hovering and unsteady nonlifting level flights.Aerodynamic results including section Cp distribution, thrust, and torque show excellent agreements in all computations and with wind-tunnel experimental data.

Mapped Chebyshev Pseudospectral Method
A Chebyshev pseudospectral method uses orthogonal Lagrange polynomials to evaluate state and control parameters of a given function to approximate or reconstruct.Either Legendre-Gauss-Lobatto (LGL) or Chebyshev-Gauss-Lobatto (CGL) points are used to evaluate function derivatives.However, the LGL or CGL points are heavily clustered around the end points of the interval, and the corresponding spectral derivative matrix becomes illconditioned with a large condition number.If the original function shows strong nonlinearity around the center of the interval where only a few CGL points are available for the spectral interpolation, the prediction of function value becomes very inaccurate around the region.The uneven distribution of the collocation points also causes numerical instability of the Chebyshev pseudospectral method [24][25][26][27][28]. Conformal mapping based on the arcsin function is applied to the CGL points to reduce the condition number of the spectral derivative matrix and improve the stability.[29].An arbitrary function () defined in time over the interval of [, ] can be represented as Chebyshev polynomials with  + 1 GCL points:

Chebyshev Collocation Operator
for  = 0, 1, 2, . . ., A vector of polynomial coefficients is written as A vector of function derivatives can be written similarly as then the coefficient vector of the function derivative is represented as A variable of  , in ( 7) is defined only when the value of (+) becomes an even number and is zero for the odd number.
For the interval of [, ], not [−1, 1], the elements of  , are divided by ( − )/2.Therefore, the time derivative of the arbitrary function is summarized as where  ch represents the Chebyshev collocation operator and is a matrix of size [ + 1,  + 1].

Conformal Mapping.
The purpose of the conformal mapping is to mitigate the problem of ill-conditioning of the spectral derivative matrix.Collocation points or nodes in a given time interval are rewritten by the mapping function of (  ), where   is the standard Chebyshev collocation points represented as a cosine mesh.
To have an interval bounded by [0, ], not by [−1, 1], a linear transformation is applied as Then, the function derivative of   is written by the chain rule, For the conformal mapping function of (), the arcsin function is used which was originally suggested by Kosloff and Tal-Ezer [25,26].A distribution characteristic is controlled by the constants of ,  with the value of  close to one resulting in more even distribution whereas smaller values push the collocation points towards the two ends of the interval.Analytic form of () is Variables of  and , shown in (11), are the functions of , , and time and can be written as below over the time interval of [ 0 ,  1 ], where  0 is typically set to be zero for the interval of time.
The interval size of () is shown in Figure 1.The graph shows the distribution of Δ() of inverse sine function.
The number of mapping points is 37, where 0 < () <  max .The distribution of the inverse sine mapping function is consistent and uniform and is less sensitive to the () values.A derivative function of the mapping function,   , is as follows: It is a diagonal matrix of size ( + 1) with zero off-diagonal terms.Therefore, the mapped Chebyshev collocation operator with the mapping function applied is represented as The size of the matrix shown in ( 16) is ( + 1) × ( + 1), and it represents the mapped Chebyshev pseudospectral collocation operator mapped to [0, ] by a tangent function.
The value of  ch corresponds to the interval of [−1, 1] as 1/  takes into account the [0,  1 ].Table 1 shows the values of the condition number for the standard and the mapped Chebyshev spectral derivative matrices with the total number of collocation points increasing from 24 to 96 by a factor of 2, where the values of  and  are optimized for the lowest condition number for each case.Reduction of the condition number by an order of magnitude is shown in all cases in Table 1 with the application of the conformal mapping of the arcsin function shown in (16).

Implementation for Unsteady Reynolds-Averaged Navier-Stokes Equations
Unsteady Reynolds-averaged Navier-Stokes (URANS) equations in three dimensions are written as where  represents conservative variables and the solutions of RANS equations.Flux vectors in , , and  directions are , , and , respectively, and include both viscous and inviscid fluxes.To apply the spectral method, vectors of solutions and fluxes of , , , and  are approximated by Chebyshev polynomials with  + 1 CGL points.
Unlike the Fourier spectral method which uses a discrete Fourier series with equally distributed collocation points, the cosine mesh points of the standard Chebyshev polynomials have uneven distribution and those are redistributed by the mapping function of () as After we replace the time derivative term in (17) by the mapped Chebyshev pseudospectral derivative operator, then (17) can be rewritten as where  → ,  →  ,  → , and  →  are the vectors of size ( + 1), and their elements correspond to mapped collocation points by the arcsin function.
The flux terms of  → , and  → are linearized with flux Jacobian matrices of , , and ,  → where , , and  are the flux Jacobian matrices of size ( + 1) × ( + 1) with each element to be a 5 × 5 submatrix and can be expressed as Substituting ( 24) in ( 23) can be expressed as The flux term in ( 23) is a function of variables of each subtime level, and (25) for each element except the diagonal term can be considered 0. Therefore, simplified implicit flux Jacobian matrix is represented as The element of the flux Jacobians of , , and  is a 5 × 5 submatrix and corresponds to the mapped collocation point.A noteworthy point here is that a simplified diagonal Jacobian matrix is the same as the flux Jacobian matrix in time domain solver.Therefore, the additional Jacobian matrix does not create in this method to be able to apply the Chebyshev pseudospectral method easily.With (29) substituted and multiplied by infinitesimally small time step Δ, (26) is rewritten as where The final form of the mapped Chebyshev pseudospectral method implemented in the URANS equations has a similar discretized structure of the steady form of RANS equations with Chebyshev derivative source terms added to the right hand side.

Numerical Validation
4.1.Pitching Airfoil with a Given Frequency.Helicopter rotor blades are in pitching and flapping motions during the forward flight to overcome the imbalance of aerodynamic forces and moments between the advancing and the retreating side.Flows on the advancing side have higher speed due to the positive contribution from rotation speed, whereas flows on the retreating side experience lower speed due to negative contribution.This causes imbalance of the aerodynamic forces and moments, and to overcome the issue a pitch angle is lowered on the advancing side and increased on the retreating side.In earlier research on the helicopter rotor blade design, aerodynamic analyses of pitching airfoil have been carried out in various wind-tunnel facilities to investigate the associated flow characteristics.Experimental data from AGARD Report 702 [30] are used in the present study to validate the mapped Chebyshev pseudospectral method.A C-type mesh of size 309 × 81 is constructed and shown in Figure 2 to solve two-dimensional URANS equations.Mach number is 0.755 and a pitching motion is described by where the damping frequency is   = 0.1628 ∘ , the mean angle of attack is  0 = 0.016 ∘ , and the pitching amplitude is   = 2.51 ∘ . * is the nondimensional time and Reynolds number based on the chord length is 5.5 × 10 6 .For the pseudo time integration, the mapped Chebyshev spectral derivative term, a DADI (diagonal alternate direction implicit), method is applied [15,23,31].A Roe's upwind flux-difference splitting (FDS) method [32] is employed for the flux term discretization with a third order MUSCL scheme to improve solution accuracy.A two-equation turbulence model of - WD+ [33]  through the Fourier spectral interpolation.With slight discrepancies in the lift and the drag coefficient in computations with 4 and 6 subtime intervals, the present method with subtime intervals more than 6 seems to converge with good agreements with those of the HBM computation.

Plunging Airfoil. One of the biggest advantages of the
Chebyshev pseudospectral method is the capability to solve nonperiodic unsteady problems if the positions at the two end points of the time interval are known a priori.The Fourier spectral method, or the nonlinear frequency domain method, uses the Fourier transform which requires the specification of the motion frequency a priori with the periodic condition imposed at the ends of the time interval; however, the Chebyshev pseudospectral method has flexibility in having different values at the two boundary points of the time interval and  is applicable to nonperiodic problems.Thus, the second validation problem is chosen as a plunging airfoil which moves downwards at a constant plunging speed prespecified in (34), where external gravity force is ignored.Freestream Mach number is 0.8 in a horizontal direction with the initial angle of attack of the airfoil to be zero.With the end point of the time interval arbitrarily set, Chebyshev collocation points are distributed based on the plunging speed and distance (Figure 5).
Moreover, it is reported in the reference paper [23] that unsteady flows around an airfoil plunging at constant speed with nonzero freestream velocity are equivalent to steady flows of an stationary airfoil fixed at a certain angle of attack with the same freestream velocity, which is explained graphically in Figures 6(a) and 6(b).An overset mesh topology for the plunging airfoil is constructed to facilitate the dynamic motion of the airfoil with a near-body grid of 257 × 45 and a background grid of 301 × 151 as shown in Figure 7. Figure 8(a) shows the schematics of the conventional time-marching method with the overset grid technique, which uses a local, finite difference-based approach.The overset mesh technique consists of two major steps, hole-cutting and donor cell identification.Wall boundaries are used for the hole-cutting and the cutting boundaries are constituted with a Cartesian mesh.The cutting surfaces are approximately made up with a collection of cubes.A zone of interface scheme and a hole-map algorithm are used for the chimera hole-cutting.For the donor cell identification, a stencil walk and gradient searching method is applied [34].The stencil walk finds the nearest point to a given interpolation point by comparing the distance between the interpolation point and the donor candidate cells.The gradient search method checks whether the interpolation point actually lies inside the cube that is made with the eight donor cells or not [35].For the time-marching method, time integration is carried out by a finite difference method, and therefore the application of the body movement, subgrid movement, hole-cutting and donor cell identification, and solving the RANS equations per dual-time step are made iteratively over the entire interval until flow solutions converge.An initial transient flow state needs to be fully resolved before flow reaches a periodic steady-state.Therefore, the overset grid technique needs to cut surfaces with a zone of interface scheme and hole-map algorithm and to interpolate with donor cell identification every physical time step.On the other hand, the mapped Chebyshev pseudospectral method constructs and stores a set of subgrid meshes, hole-cutting, and donor cell information corresponding to the mapped CGL points once at the beginning of the solution process, and the information is not modified throughout the iterations until flow solutions converge in a pseudo-time domain.This fact is greatly favorable in efficiency by reducing computation time considerably, which is more so when the size of grid becomes larger.Although it should be noted that computation memory increases accordingly from storing the meshes and flow solutions at all subgrid intervals, the overall computation efficiency has been shown higher in the Chebyshev pseudospectral computation for this plunging airfoil problem and the helicopter rotor flow analysis discussed in Section 5; Figures 9 and 10 show the contours of pressure and Mach number, respectively, around the airfoils.Contours shown in Figures 9(a) and 10(a) are computed for the plunging airfoil by the mapped Chebyshev pseudospectral method, whereas those shown in Figures 9(b) and 10(b) are solved by the steady RANS computation for the stationary airfoil with the angle of attack fixed at 8 ∘ .A single block, C-type structured mesh is used for the steady RANS computation.Comparisons shown in Figures 9 and 10 show excellent agreements between the two computations, although slightly nonsmooth contours are observed in the Chebyshev pseudospectral computation especially in the overlapping mesh regions between the near-body and the background grid.Figure 11 shows a direct comparison of the surface pressure coefficients with good agreements.condition of the helicopter rotor flow analysis [36].The rotor system has two blades which have a rectangular-shape planform with no twist and taper in the span direction, and their cross section is NACA0012 airfoil.The span length is 3.75 ft (1.143 m) with an aspect ratio of 7. A collective pitch angle of 8 ∘ is applied to the blades at all azimuth angles with Mach numbers of 0.439 and 0.877.All computations are carried out in parallel in a Linux cluster using a total of 56 CPUs of i7 Core processors with 24 GB RAM per node, which is effective in handling the increased memory requirement for the mapped Chebyshev pseudospectral method with meshes and solutions at all subtime intervals.Figure 12 shows the multiblock, overset mesh topology used for the Caradonna and Tung's rotor.A near-body mesh is H-type with 155 × 239 × 239 points and the background mesh has the grid of 97 × 106 × 217.A total of 11 Chebyshev collocation points, corresponding to 11 subtime intervals due to the periodicity of the rotor motion, are used for the analysis of the hovering flight.Computations of the time-marching method, the HBM, and the mapped Chebyshev pseudospectral method are carried out independently and resultant aerodynamic forces and moments are compared with the wind-tunnel

Helicopter Rotor Flow Analysis
for  = 0, 1, 2, . . ., ,    Chebyshev points.Computation results with a total of 11 points predict relatively well the shock structure in terms of the location and magnitude.Figure 18 shows the contours of surface pressure distribution on the blades at the mapped Chebyshev points which vary from 7, 9, to 11. Figure 18(d) shows the contours of surface pressure distribution on the blades reconstructed from the 11 points to 18 points.We checked the wall clock time with increasing interval numbers of HBM and present method in Figure 19.For time accurate method the dual-time stepping is applied with 360 intervals for one cycle and 20 subiteration for each interval, and three cycles for the periodic convergence of the unsteady flow.The computational time of present method is increasing linearly with respect to the number of intervals as like HBM, and Figure 19 shows present method is as efficient as HBM below 19 intervals.

Conclusions
The mapped Chebyshev pseudospectral method was developed and implemented for URANS solvers with the overset mesh topology.It uses Chebyshev polynomials for solution approximation with the CGL points for the collocation nodes.
The time derivative term of the governing equations is replaced by mapped Chebyshev spectral derivative term and discretized in explicit form for time integration in a pseudo time domain.Flux terms are discretized in implicit form and this different discretization strategy makes the integration of the PDE governing equations particularly efficient in both time and space.The governing equations in discrete form are similar to those of the steady RANS equations.Unlike a discrete Fourier series, Chebyshev polynomials have uneven distribution of the collocation points, which often leads to the ill-conditioned spectral derivative matrix and deteriorates the numerical stability of the linear system.An arcsin conformal mapping function is applied for redistributing the Chebyshev points more evenly with optimum value of the control parameters of  and  to be 0.99.For the validation of the present method, airfoils in oscillation and plunging motions were solved separately to investigate the solution accuracy of the method in periodic and nonperiodic flows.The number of subtime intervals was also varied to demonstrate the sensitivity of the solution accuracy to the number of Chebyshev points.Convergence rate did not change to the varying number of subtime intervals, although computation time and memory to compute one pseudo-time step increase linearly.Only a small number of Chebyshev points were needed to accurately predict the airfoil solution problems.Subsequently, the Chebyshev pseudospectral method was applied to the helicopter rotor flow problems in hovering and nonlifting forward flight conditions using the Caradonna and Tung rotor.As a set of subgrids are constructed once at the beginning of the solution procedure only at given mapped Chebyshev points along with the hole-cutting method for solution interpolation, computation time is saved dramatically in comparison to the time-marching computation which carries out the mesh movement with the hole-cutting at each physical time step which is often limited to less than 10 due to the numerical stability limit.The time saving was on the orders of magnitude.Aerodynamic performance of thrust and torque as well as surface pressure coefficients at several radial locations around azimuth angles shows excellent agreements when compared to those of the timemarching and the HBM method.All computation results show good agreements with wind-tunnel experimental data for both flight conditions.The number of subtime interval more than 11 is needed for the solution accuracy of the forward flight, while the corresponding computation time increases linearly with respect to the number of the subtime intervals.Overall, the present approach of the mapped Chebyshev pseudospectral method is efficient in solving both periodic and nonperiodic problems at greatly reduced cost

Figure 6 :
Figure 6: An airfoil plunging at constant speed ẏ and a stationary airfoil at the AOA of .
(b) Results by the steady computation for the stationary airfoil with  = 0
(a) Results by the Chebyshev spectral method for the plung-Results by the steady computation for the stationary airfoil with  = 0

Figure 11 :Figure 12 :
Figure 11: Comparison of the surface pressure coefficients.

Figure 13 :
Figure 13: Schematic for solution reconstruction using the mapped Chebyshev points.

Figure 16 :
Figure 16: Contour plots of the blade surface pressure at Mach = 0.877.

Cheby 13 ＭFigure 17 :
Figure 17: Pressure coefficients at the 0.89 radial location during the forward flight with an azimuth angle varying from 30 deg to 180 deg by 30 deg.

Figure 18 :
Figure 18: Surface pressure distribution on the blades during the forward flight at the mapped Chebyshev points varying from seven, nine, and eleven in total.Surface pressure distribution on the blades at the reconstructed at 11 mapped Chebyshev points.

Figure 19 :
Figure 19: Wall clock time (in seconds) along with the number of interval points.

Table 1 :
Comparison of condition number for the differential operator.
) term and the source term can give the implicit convergence characteristic.In other words, the convergence rate of the existing time domain solver is maintained.First, the combination of explicit source and implicit residual value is derived as pseudospectral approach.Chebyshev pseudospectral source term is considered explicitly, so the influence of the source term in nondiagonal matrix is ignored.For this reason the source term is added in flux term at subtime level.Though the Chebyshev pseudospectral source term is used explicitly, Diagonal ADI technique that treats together with the flux 5.1.Hovering Condition.A geometry of the Caradonna andTung's rotor is popularly used to validate the hovering Results by the Chebyshev spectral method for the plunging airfoil (a)

Table 2 :
Aerodynamic coefficients at Mach number of 0.439.and,aswell as the figure of merit (FM) is shown in Table2where the computations are carried with the Mach number of 0.439.The same Caradonna and Tung rotors are solved at a nonlifting forward flight condition.Blades are not in a pitching or flapping motion and a collective pitch angle is set to be zero.Only forward speed is considered with corresponding tip Mach number of 0.5 and the advance ratio of 0.2.The CPU time is shown in Table3with respect to the number of Chebyshev collocation points from 11 to 23 and is normalized by the computational time applied with

Table 3 :
[37]imensionalized time with respect to the varying number of Chebyshev collocation points and corresponding grid module time for mesh interpolation and hole cutting.Table3also shows overset grid module time which is defined as the time for hole-cutting and donor cell identification with respect to the number of Chebyshev collocation points.For the time-marching computation for the unsteady analysis with this overset grid technique, the grid module time is about 3.4 seconds per physical time step which means 1 ∘ azimuth angle, which is inferred to be more than 20 minutes for one revolution which means 360 ∘ azimuth angle.The rotor flow usually converges after five revolutions and it only takes 1.7 hours for overset grid module time.However, by the present mapped Chebyshev pseudospectral method, the grid module time for the holecutting and mesh interpolation needs is needed only once at the beginning of the solution procedure and takes less than 2 minutes even with 23 Chebyshev points.Figure17shows the plots of the pressure coefficients at the 0.89 radial location of the blade when the azimuth angle varies from 30 deg to 180 deg by 30 deg.Results of the mapped Chebyshev pseudospectral method are computed with a total of 11, 13, and 15 collocation points and compared with the wind-tunnel data[37].Strong shock appears in the advancing side of the rotor, approximately from 60 deg to 150 deg of the azimuth angle, demonstrating solution accuracy higher with more