Computational aspects of zonal algorithms for solving the compressible Navier-Stokes equations in three dimensions

Transonic flow fields about wing geometries are computed using an Euler/Navier-Stokes approach in which the flow field is divided into several zones. The flow field immediately adjacent to the wing surface is resolved with fine grid zones and solved using a Navier-Stokes algorithm. Flow field regions removed from the wing are resolved with less finely clustered grid zones and are solved with an Euler algorithm. Computational issues associated with this zonal approach, including data base management aspects, are discussed. Solutions are obtained that are in good agreement with experiment, including cases with significant wind tunnel wall effects. Additional cases with significant shock induced separation on the upper wing surface are also presented

Natural Rolling Responses of a Delta Wing in Transonic and Subsonic Flows

The unsteady, three-dimensional, full Navier-Stokes (NS) equations and the Euler equations of rigid-body dynamics are sequentially solved to simulate the natural rolling response of slender delta wings of zero thickness at moderate to high angles of attack, to transonic and subsonic flows. The governing equations of fluid flow and dynamics of the present multi-disciplinary problem are solved using the time-accurate solution of the NS equations with the implicit, upwind, Roe flux-difference splitting, finite-volume scheme and a four-stage Runge-Kutta scheme, respectively. The main focus is to analyze the effect of Mach number and angle of attack on the leading edge vortices and their breakdown, the resultant rolling motion, and overall aerodynamic response of the wing. Three cases demonstrate the natural response of a 65 deg swept, cropped delta wing in a transonic flow with breakdown of the leading edge vortices and an 80 deg swept delta wing in a subsonic flow undergoing either damped or self-excited limit-cycle rolling oscillations as a function of angle of attack. Comparisons with an experimental investigation completes this study, validating the analysis and illustrating the complex details afforded by computational investigations

Time-Spectral Rotorcraft Simulations on Overset Grids

The Time-Spectral method is derived as a Fourier collocation scheme and applied to NASA's overset Reynolds-averaged Navier-Stokes (RANS) solver OVERFLOW. The paper outlines the Time-Spectral OVERFLOWimplementation. Successful low-speed laminar plunging NACA 0012 airfoil simulations demonstrate the capability of the Time-Spectral method to resolve the highly-vortical wakes typical of more expensive three-dimensional rotorcraft configurations. Dealiasing, in the form of spectral vanishing viscosity (SVV), facilitates the convergence of Time-Spectral calculations of high-frequency flows. Finally, simulations of the isolated V-22 Osprey tiltrotor for both hover and forward (edgewise) flight validate the three-dimensional Time-Spectral OVERFLOW implementation. The Time-Spectral hover simulation matches the time-accurate calculation using a single harmonic. Significantly more temporal modes and SVV are required to accurately compute the forward flight case because of its more active, high-frequency wake

Assessment of the Vortex Particle-Mesh Method for Efficient LES of Hovering Rotors and their Wakes

Midway between a vortex method and a grid-based CFD, the Vortex Particle-Mesh method with Immersed Lifting Lines is intended to provide medium-fidelity results on rotor loadings, together with a realistic representation of the 3-D vortical wakes and their dynamics over long distances. We assess the potential of this hybrid approach for the Large-Eddy Simulation of helicopter rotors in hover. A novel Poisson solver with mixed unbounded-outflow boundary conditions here further enables the computation of turbulent hovering scenarios in tight domains. Considering the Knight and Hefner experiment and the S-76 test case as references, we present and compare blade integrated and distributed loads, induction velocities, and wake characteristics. While the quality of the performance predictions highly depends on external polar data, the obtained wakes exhibit similar characteristics to those recently identified in other CFD analyses, here at a moderate computational cost. Based on these results, we further investigate the secondary vortex structures forming between the main tip vortices, and we bring additional clues on their relation to the phenomenon of wake breakdown. We finally discuss the strengths and difficulties of hybrid vortex methods for the challenging analysis of hovering rotors