Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shear Layers”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,117 records · Page 62

Initiation of Long-Wave Instability of Vortex Pairs at Cruise Altitudes

Previous studies have usually attributed the initiation of the long-wave instability of a vortex pair to turbulence in the atmosphere or in the wake of the aircraft. The purpose here is to show by use of observations and photographs of condensation trails shed by aircraft at cruise altitudes that another initiating mechanism is not only possible but is usually the mechanism that initiates the long-wave instability at cruise altitudes. The alternate initiating mechanism comes about when engine thrust is robust enough to form an array of circumferential vortices around each jet-engine-exhaust stream. In those cases, initiation begins when the vortex sheet shed by the wing has rolled up into a vortex pair and descended to the vicinity of the inside bottom of the combined shear-layer vortex arrays. It is the in-and-out (up and down) velocity field between sequential circumferential vortices near the bottom of the array that then impresses disturbance waves on the lift-generated vortex pair that initiate the long-wave instability. A time adjustment to the Crow and Bate estimate for vortex linking is then derived for cases when thrust-based linking occurs.

Rossow, Vernon J.↗

The Prediction of Broadband Shock-Associated Noise Including Propagation Effects

An acoustic analogy is developed based on the Euler equations for broadband shock- associated noise (BBSAN) that directly incorporates the vector Green's function of the linearized Euler equations and a steady Reynolds-Averaged Navier-Stokes solution (SRANS) as the mean flow. The vector Green's function allows the BBSAN propagation through the jet shear layer to be determined. The large-scale coherent turbulence is modeled by two-point second order velocity cross-correlations. Turbulent length and time scales are related to the turbulent kinetic energy and dissipation. An adjoint vector Green's function solver is implemented to determine the vector Green's function based on a locally parallel mean flow at streamwise locations of the SRANS solution. However, the developed acoustic analogy could easily be based on any adjoint vector Green's function solver, such as one that makes no assumptions about the mean flow. The newly developed acoustic analogy can be simplified to one that uses the Green's function associated with the Helmholtz equation, which is consistent with the formulation of Morris and Miller (AIAAJ 2010). A large number of predictions are generated using three different nozzles over a wide range of fully expanded Mach numbers and jet stagnation temperatures. These predictions are compared with experimental data from multiple jet noise labs. In addition, two models for the so-called 'fine-scale' mixing noise are included in the comparisons. Improved BBSAN predictions are obtained relative to other models that do not include the propagation effects, especially in the upstream direction of the jet.

Miller, Steven↗

Time Resolved PIV for Space-Time Correlations in Hot Jets

Temporally Resolved Particle Image Velocimetry (TR-PIV) is the newest and most exciting tool recently developed to support our continuing efforts to characterize and improve our understanding of the decay of turbulence in jet flows -- a critical element for understanding the acoustic properties of the flow. A new TR-PIV system has been developed at the NASA Glenn Research Center which is capable of acquiring planar PIV image frame pairs at up to 25 kHz. The data reported here were collected at Mach numbers of 0.5 and 0.9 and at temperature ratios of 0.89 and 1.76. The field of view of the TR-PIV system covered 6 nozzle diameters along the lip line of the 50.8 mm diameter jet. The cold flow data at Mach 0.5 were compared with hotwire anemometry measurements in order to validate the new TR-PIV technique. The axial turbulence profiles measured across the shear layer using TR-PIV were thinner than those measured using hotwire anemometry and remained centered along the nozzle lip line. The collected TR-PIV data illustrate the differences in the single point statistical flow properties of cold and hot jet flows. The planar, time-resolved velocity records were then used to compute two-point space-time correlations of the flow at the Mach 0.9 flow condition. The TR-PIV results show that there are differences in the convective velocity and growth rate of the turbulent structures between cold and hot flows at the same Mach number

Wernet, Mark P.↗

Underexpanded Supersonic Plume Surface Interactions: Applications for Spacecraft Landings on Planetary Bodies

Numerical and experimental investigations of both far-field and near-field supersonic steady jet interactions with a flat surface at various atmospheric pressures are presented in this paper. These studies were done in assessing the landing hazards of both the NASA Mars Science Laboratory and Phoenix Mars spacecrafts. Temporal and spatial ground pressure measurements in conjunction with numerical solutions at altitudes of approx.35 nozzle exit diameters and jet expansion ratios (e) between 0.02 and 100 are used. Data from steady nitrogen jets are compared to both pulsed jets and rocket exhaust plumes at Mach approx.5. Due to engine cycling, overpressures and the plate shock dynamics are different between pulsed and steady supersonic impinging jets. In contrast to highly over-expanded (e <1) and underexpanded exhaust plumes, results show that there is a relative ground pressure load maximum for moderately underexpanded (e approx.2-5) jets which demonstrate a long collimated plume shock structure. For plumes with e much >5 (lunar atmospheric regime), the ground pressure is minimal due to the development of a highly expansive shock structure. We show this is dependent on the stability of the plate shock, the length of the supersonic core and plume decay due to shear layer instability which are all a function of the jet expansion ratio. Asymmetry and large gradients in the spatial ground pressure profile and large transient overpressures are predominantly linked to the dynamics of the plate shock. More importantly, this study shows that thruster plumes exhausting into martian environments possess the largest surface pressure loads and can occur at high spacecraft altitudes in contrast to the jet interactions at terrestrial and lunar atmospheres. Theoretical and analytical results also show that subscale supersonic cold gas jets adequately simulate the flow field and loads due to rocket plume impingement provided important scaling parameters are in agreement. These studies indicate the critical importance of testing and modeling plume-surface interactions for descent and ascent of spacecraft and launch vehicles.

Mehta, M.↗

Acoustic Performance of an Advanced Model Turbofan in Three Aeroacoustic Test Facilities

A model advanced turbofan was acoustically tested in the NASA Glenn 9- by 15-Foot-Low-Speed Wind Tunnel (LSWT), and in two other aeroacoustic facilities. The Universal Propulsion Simulator (UPS) fan was designed and manufactured by the General Electric Aircraft Engines (GEAE) Company, and featured active core, as well as bypass, flow paths. The reference test configurations were with the metal, M4, rotor with hardwall and treated bypass flow ducts. The UPS fan was tested within an airflow at a Mach number of 0.20 (limited flow data were also acquired at a Mach number of 0.25) which is representative of aircraft takeoff and approach conditions. Comparisons were made between data acquired within the airflow (9x15 LSWT and German-Dutch Wind Tunnel (DNW)) and outside of a free jet (Boeing Low Speed Aero acoustic Facility (LSAF) and DNW). Sideline data were acquired on an 89-in. (nominal 4 fan diameters) sideline using the same microphone assembly and holder in the 9x15 LSWT and DNW facilities. These data showed good agreement for similar UPS operating conditions and configurations. Distortion of fan spectra tonal content through a free jet shear layer was documented, suggesting that in-flow acoustic measurements are required for comprehensive fan noise diagnostics. However, there was good agreement for overall sound power level (PWL) fan noise measurements made both within and outside of the test facility airflow.

Woodward, Richard P.↗

Sofia Observatory Performance and Characterization

The Stratospheric Observatory for Infrared Astronomy (SOFIA) has recently concluded a set of engineering flights for Observatory performance evaluation. These in-flight opportunities have been viewed as a first comprehensive assessment of the Observatory's performance and will be used to address the development activity that is planned for 2012, as well as to identify additional Observatory upgrades. A series of 8 SOFIA Characterization And Integration (SCAI) flights have been conducted from June to December 2011. The HIPO science instrument in conjunction with the DSI Super Fast Diagnostic Camera (SFDC) have been used to evaluate pointing stability, including the image motion due to rigid-body and flexible-body telescope modes as well as possible aero-optical image motion. We report on recent improvements in pointing stability by using an Active Mass Damper system installed on Telescope Assembly. Measurements and characterization of the shear layer and cavity seeing, as well as image quality evaluation as a function of wavelength have been performed using the HIPO+FLITECAM Science Instrument configuration (FLIPO). A number of additional tests and measurements have targeted basic Observatory capabilities and requirements including, but not limited to, pointing accuracy, chopper evaluation and imager sensitivity. SCAI activities included in-flight partial Science Instrument commissioning prior to the use of the instruments as measuring engines. This paper reports on the data collected during the SCAI flights and presents current SOFIA Observatory performance and characterization.

Temi, Pasquale↗

Prediction of Broadband Shock-Associated Noise Including Propagation Effects Originating NASA

An acoustic analogy is developed based on the Euler equations for broadband shock-associated noise (BBSAN) that directly incorporates the vector Green s function of the linearized Euler equations and a steady Reynolds-Averaged Navier-Stokes solution (SRANS) to describe the mean flow. The vector Green s function allows the BBSAN propagation through the jet shear layer to be determined. The large-scale coherent turbulence is modeled by two-point second order velocity cross-correlations. Turbulent length and time scales are related to the turbulent kinetic energy and dissipation rate. An adjoint vector Green s function solver is implemented to determine the vector Green s function based on a locally parallel mean flow at different streamwise locations. The newly developed acoustic analogy can be simplified to one that uses the Green s function associated with the Helmholtz equation, which is consistent with a previous formulation by the authors. A large number of predictions are generated using three different nozzles over a wide range of fully-expanded jet Mach numbers and jet stagnation temperatures. These predictions are compared with experimental data from multiple jet noise experimental facilities. In addition, two models for the so-called fine-scale mixing noise are included in the comparisons. Improved BBSAN predictions are obtained relative to other models that do not include propagation effects.

Miller, Steven↗

The Prediction of Jet Noise Ground Effects Using an Acoustic Analogy and a Tailored Green's Function

An assessment of an acoustic analogy for the mixing noise component of jet noise in the presence of an infinite surface is presented. The reflection of jet noise by the ground changes the distribution of acoustic energy and is characterized by constructive and destructive interference patterns. The equivalent sources are modeled based on the two-point cross- correlation of the turbulent velocity fluctuations and a steady Reynolds-Averaged Navier-Stokes (RANS) solution. Propagation effects, due to reflection by the surface and refaction by the jet shear layer, are taken into account by calculating the vector Green's function of the linearized Euler equations (LEE). The vector Green's function of the LEE is written in relation to Lilley's equation; that is, approximated with matched asymptotic solutions and the Green's function of the convective Helmholtz equation. The Green's function of the convective Helmholtz equation for an infinite flat plane with impedance is the Weyl-van der Pol equation. Predictions are compared with an unheated Mach 0.95 jet produced by a nozzle with an exit diameter of 0.3302 meters. Microphones are placed at various heights and distances from the nozzle exit in the peak jet noise direction above an acoustically hard and an asphalt surface. The predictions are shown to accurately capture jet noise ground effects that are characterized by constructive and destructive interference patterns in the mid- and far-field and capture overall trends in the near-field.

Miller, Steven A. E.↗

DNS, Enstrophy Balance, and the Dissipation Equation in a Separated Turbulent Channel Flow

The turbulent flows through a plane channel and a channel with a constriction (2-D hill) are numerically simulated using DNS and RANS calculations. The Navier-Stokes equations in the DNS are solved using a higher order kinetic energy preserving central schemes and a fifth order accurate upwind biased WENO scheme for the space discretization. RANS calculations are performed using the NASA code CFL3D with the komega SST two-equation model and a full Reynolds stress model. Using DNS, the magnitudes of different terms that appear in the enstrophy equation are evaluated. The results show that the dissipation and the diffusion terms reach large values at the wall. All the vortex stretching terms have similar magnitudes within the buffer region. Beyond that the triple correlation among the vorticity and strain rate fluctuations becomes the important kinematic term in the enstrophy equation. This term is balanced by the viscous dissipation. In the separated flow, the triple correlation term and the viscous dissipation term peak locally and balance each other near the separated shear layer region. These findings concur with the analysis of Tennekes and Lumley, confirming that the energy transfer terms associated with the small-scale dissipation and the fluctuations of the vortex stretching essentially cancel each other, leaving an equation for the dissipation that is governed by the large-scale motion.

Balakumar, Ponnampalam↗

Iodine Tagging Velocimetry and Mechanism in the Hypersonic Near Wake of a MultiPurpose Crew Vehicle

This study demonstrates a new molecular tagging velocimetry (MTV) method for velocity measurements of high speed flow. It demonstrates offbody Iodine Tagging Velocimetry (ITV) in the hypersonic near wake of a MultiPurpose Crew Vehicle (MPCV) model. Experiments are performed in the NASA-Langley 31-inch Mach 10 air wind tunnel. A 0.5% I2 / N2 mixture is seeded on the leeward backshell of the model using a pressure tap. I2 laser-induced fluorescence is excited along a 5.5 mm line using an ArF excimer laser near 193 nm. Results indicate I2 absorbs at least 2 photons to produce iodine ions and electrons. These recombine as the tagged region is displaced downstream to produce I (2P3/2) whose emission is monitored at 206 nm. Results at P0 = 2.41 MPa (350 psi), T0 = 990K, and 10 micro-sec transit times produce velocities from 630-820 m/sec across the I2 seeded jet at a distance of 38.2 mm (25.5 jet diameters) downstream from the jet orifice. Maximum wake jet velocities near the shear layer are 59% of freestream velocity.

Balla, R. Jeffrey↗

Acoustic Data Processing and Transient Signal Analysis for the Hybrid Wing Body 14- by 22-Foot Subsonic Wind Tunnel Test

An advanced vehicle concept, the HWB N2A-EXTE aircraft design, was tested in NASA Langley's 14- by 22-Foot Subsonic Wind Tunnel to study its acoustic characteristics for var- ious propulsion system installation and airframe con gurations. A signi cant upgrade to existing data processing systems was implemented, with a focus on portability and a re- duction in turnaround time. These requirements were met by updating codes originally written for a cluster environment and transferring them to a local workstation while en- abling GPU computing. Post-test, additional processing of the time series was required to remove transient hydrodynamic gusts from some of the microphone time series. A novel automated procedure was developed to analyze and reject contaminated blocks of data, under the assumption that the desired acoustic signal of interest was a band-limited sta- tionary random process, and of lower variance than the hydrodynamic contamination. The procedure is shown to successfully identify and remove contaminated blocks of data and retain the desired acoustic signal. Additional corrections to the data, mainly background subtraction, shear layer refraction calculations, atmospheric attenuation and microphone directivity corrections, were all necessary for initial analysis and noise assessments. These were implemented for the post-processing of spectral data, and are shown to behave as expected.

Bahr, Christopher J.↗

Direct Numerical Simulations of Transitional/Turbulent Wakes

The interest in transitional/turbulent wakes spans the spectrum from an intellectual pursuit to understand the complex underlying physics to a critical need in aeronautical engineering and other disciplines to predict component/system performance and reliability. Cylinder wakes have been studied extensively over several decades to gain a better understanding of the basic flow phenomena that are encountered in such flows. Experimental, computational and theoretical means have been employed in this effort. While much has been accomplished there are many important issues that need to be resolved. The physics of the very near wake of the cylinder (less than three diameters downstream) is perhaps the most challenging of them all. This region comprises the two detached shear layers, the recirculation region and wake flow. The interaction amongst these three components is to some extent still a matter of conjecture. Experimental techniques have generated a large percentage of the data that have provided us with the current state of understanding of the subject. More recently computational techniques have been used to simulate cylinder wakes, and the data from such simulations are being used to both refine our understanding of such flows as well as provide new insights. A few large eddy and direct numerical simulations (LES and DNS) of cylinder wakes have appeared in the literature in the recent past. These investigations focus on the low Reynolds number range where the cylinder boundary layer is laminar (sub-critical range). However, from an engineering point of view, there is considerable interest in the situation where the upper and/or lower boundary layer of an airfoil is turbulent, and these turbulent boundary layers separate from the airfoil to contribute to the formation of the wake downstream. In the case of cylinders, this only occurs at relatively large unit Reynolds numbers. However, in the case of airfoils, the boundary layer has the opportunity to transition to turbulence on the airfoil surface at a relatively lower unit Reynolds number because the characteristic length of the airfoil is typically one to two orders of magnitude larger than the trailing edge diameter. This transition to turbulence would occur unless there is a strong favorable pressure gradient that results in the boundary layer remaining laminar or transitional over the surface of the airfoil. This presentation will focus on two direct numerical simulations that have been performed at NASA ARC. The first is of a cylinder wake with laminar separating boundary layers. The second is the wake of a flat plate with a circular trailing edge. The upper and lower plate surface boundary layers are both turbulent and statistically identical. Thus the computed wake is symmetric in a statistical sense. This flow is more representative of airfoil wakes than cylinder wakes. Results from the two simulations including flow visualization and turbulence statistics in the near wake will be presented at the seminar.

wakes↗

Characterization of the NASA Langley Arc Heated Scramjet Test Facility Using NO PLIF

The nitric oxide planar laser-induced fluorescence (NO PLIF) imaging was used to characterize the air flow of the NASA Langley Arc Heated Scramjet Test Facility (AHSTF) configured with a Mach 6 nozzle. The arc raises the enthalpy of the test gas in AHSTF, producing nitric oxide. Nitric oxide persists as the temperature drops through the nozzle into the test section. NO PLIF was used to qualitatively visualize the flowfield at different experimental conditions, measure the temperature of the gas flow exiting the facility nozzle, and visualize the wave structure downstream of the nozzle at different operating conditions. Uniformity and repeatability of the nozzle flow were assessed. Expansion and compression waves on the free-jet shear layer as the nozzle flow expands into the test section were visualized. The main purpose of these experiments was to assess the uniformity of the NO in the freestream gas for planned experiments, in which NO PLIF will be used for qualitative fuel-mole-fraction sensitive imaging. The shot-to-shot fluctuations in the PLIF signal, caused by variations in the overall laser intensity as well as NO concentration and temperature variations in the flow was 20-25% of the mean signal, as determined by taking the standard deviation of a set of images obtained at constant conditions and dividing by the mean. The fluctuations within individual images, caused by laser sheet spatial variations as well as NO concentration and temperature variations in the flow, were about 28% of the mean in images, determined by taking standard deviation within individual images, dividing by the mean in the same image and averaged over the set of images. Applying an averaged laser sheet intensity correction reduced the within-image intensity fluctuations to about 10% suggesting that the NO concentration is uniform to within 10%. There was no significant difference in flow uniformity between the low and high enthalpy settings. While not strictly quantitative, the temperature maps show qualitative agreement with the computations of the flow.

Kidd, F. Gray, III↗

The Prediction of Scattered Broadband Shock-Associated Noise

A mathematical model is developed for the prediction of scattered broadband shock-associated noise. Model arguments are dependent on the vector Green's function of the linearized Euler equations, steady Reynolds-averaged Navier-Stokes solutions, and the two-point cross-correlation of the equivalent source. The equivalent source is dependent on steady Reynolds-averaged Navier-Stokes solutions of the jet flow, that capture the nozzle geometry and airframe surface. Contours of the time-averaged streamwise velocity component and turbulent kinetic energy are examined with varying airframe position relative to the nozzle exit. Propagation effects are incorporated by approximating the vector Green's function of the linearized Euler equations. This approximation involves the use of ray theory and an assumption that broadband shock-associated noise is relatively unaffected by the refraction of the jet shear layer. A non-dimensional parameter is proposed that quantifies the changes of the broadband shock-associated noise source with varying jet operating condition and airframe position. Scattered broadband shock-associated noise possesses a second set of broadband lobes that are due to the effect of scattering. Presented predictions demonstrate relatively good agreement compared to a wide variety of measurements.

Miller, Steven A. E.↗

Advanced Model for Extreme Lift and Improved Aeroacoustics (AMELIA)

With the very recent advent of NASA's Environmentally Responsible Aviation Project (ERA), which is dedicated to designing aircraft that will reduce the impact of aviation on the environment, there is a need for research and development of methodologies to minimize fuel burn, emissions, and reduce community noise produced by regional airliners. ERA tackles airframe technology, propulsion technology, and vehicle systems integration to meet performance objectives in the time frame for the aircraft to be at a Technology Readiness Level (TRL) of 4-6 by the year of 2020 (deemed N+2). The proceeding project that investigated similar goals to ERA was NASA's Subsonic Fixed Wing (SFW). SFW focused on conducting research to improve prediction methods and technologies that will produce lower noise, lower emissions, and higher performing subsonic aircraft for the Next Generation Air Transportation System. The work provided in this investigation was a NASA Research Announcement (NRA) contract #NNL07AA55C funded by Subsonic Fixed Wing. The project started in 2007 with a specific goal of conducting a large-scale wind tunnel test along with the development of new and improved predictive codes for the advanced powered-lift concepts. Many of the predictive codes were incorporated to refine the wind tunnel model outer mold line design. The large scale wind tunnel test goal was to investigate powered lift technologies and provide an experimental database to validate current and future modeling techniques. Powered-lift concepts investigated were Circulation Control (CC) wing in conjunction with over-the-wing mounted engines to entrain the exhaust to further increase the lift generated by CC technologies alone. The NRA was a five-year effort; during the first year the objective was to select and refine CESTOL concepts and then to complete a preliminary design of a large-scale wind tunnel model for the large scale test. During the second, third, and fourth years the large-scale wind tunnel model design would be completed, manufactured, and calibrated. During the fifth year the large scale wind tunnel test was conducted. This technical memo will describe all phases of the Advanced Model for Extreme Lift and Improved Aeroacoustics (AMELIA) project and provide a brief summary of the background and modeling efforts involved in the NRA. The conceptual designs considered for this project and the decision process for the selected configuration adapted for a wind tunnel model will be briefly discussed. The internal configuration of AMELIA, and the internal measurements chosen in order to satisfy the requirements of obtaining a database of experimental data to be used for future computational model validations. The external experimental techniques that were employed during the test, along with the large-scale wind tunnel test facility are covered in great detail. Experimental measurements in the database include forces and moments, and surface pressure distributions, local skin friction measurements, boundary and shear layer velocity profiles, far-field acoustic data and noise signatures from turbofan propulsion simulators. Results and discussion of the circulation control performance, over-the-wing mounted engines, and the combined performance are also discussed in great detail.

Subsonic Fixed Wing↗

Characterization of Flap Edge Noise Radiation from a High-Fidelity Airframe Model

The results of an experimental study of the noise generated by a baseline high-fidelity airframe model are presented. The test campaign was conducted in the open-jet test section of the NASA Langley 14- by 22-foot Subsonic Tunnel on an 18%-scale, semi-span Gulfstream airframe model incorporating a trailing edge flap and main landing gear. Unsteady surface pressure measurements were obtained from a series of sensors positioned along the two flap edges, and far field acoustic measurements were obtained using a 97-microphone phased array that viewed the pressure side of the airframe. The DAMAS array deconvolution method was employed to determine the locations and strengths of relevant noise sources in the vicinity of the flap edges and the landing gear. A Coherent Output Power (COP) spectral method was used to couple the unsteady surface pressures measured along the flap edges with the phased array output. The results indicate that outboard flap edge noise is dominated by the flap bulb seal cavity with very strong COP coherence over an approximate model-scale frequency range of 1 to 5 kHz observed between the array output and those unsteady pressure sensors nearest the aft end of the cavity. An examination of experimental COP spectra for the inboard flap proved inconclusive, most likely due to a combination of coherence loss caused by decorrelation of acoustic waves propagating through the thick wind tunnel shear layer and contamination of the spectra by tunnel background noise at lower frequencies. Directivity measurements obtained from integration of DAMAS pressure-squared values over defined geometric zones around the model show that the baseline flap and landing gear are only moderately directional as a function of polar emission angle.

Humphreys, William M., Jr.↗

DNS/LES Simulations of Separated Flows at High Reynolds Numbers

Direct numerical simulations (DNS) and large-eddy simulations (LES) simulations of flow through a periodic channel with a constriction are performed using the dynamic Smagorinsky model at two Reynolds numbers of 2800 and 10595. The LES equations are solved using higher order compact schemes. DNS are performed for the lower Reynolds number case using a fine grid and the data are used to validate the LES results obtained with a coarse and a medium size grid. LES simulations are also performed for the higher Reynolds number case using a coarse and a medium size grid. The results are compared with an existing reference data set. The DNS and LES results agreed well with the reference data. Reynolds stresses, sub-grid eddy viscosity, and the budgets for the turbulent kinetic energy are also presented. It is found that the turbulent fluctuations in the normal and spanwise directions have the same magnitude. The turbulent kinetic energy budget shows that the production peaks near the separation point region and the production to dissipation ratio is very high on the order of five in this region. It is also observed that the production is balanced by the advection, diffusion, and dissipation in the shear layer region. The dominant term is the turbulent diffusion that is about two times the molecular dissipation.

Balakumar, P.↗

Improving Fidelity of Launch Vehicle Liftoff Acoustic Simulations

Launch vehicles experience high acoustic loads during ignition and liftoff affected by the interaction of rocket plume generated acoustic waves with launch pad structures. Application of highly parallelized Computational Fluid Dynamics (CFD) analysis tools optimized for application on the NAS computer systems such as the Loci/CHEM program now enable simulation of time-accurate, turbulent, multi-species plume formation and interaction with launch pad geometry and capture the generation of acoustic noise at the source regions in the plume shear layers and impingement regions. These CFD solvers are robust in capturing the acoustic fluctuations, but they are too dissipative to accurately resolve the propagation of the acoustic waves throughout the launch environment domain along the vehicle. A hybrid Computational Fluid Dynamics and Computational Aero-Acoustics (CFD/CAA) modeling framework has been developed to improve such liftoff acoustic environment predictions. The framework combines the existing highly-scalable NASA production CFD code, Loci/CHEM, with a high-order accurate discontinuous Galerkin (DG) solver, Loci/THRUST, developed in the same computational framework. Loci/THRUST employs a low dissipation, high-order, unstructured DG method to accurately propagate acoustic waves away from the source regions across large distances. The DG solver is currently capable of solving up to 4th order solutions for non-linear, conservative acoustic field propagation. Higher order boundary conditions are implemented to accurately model the reflection and refraction of acoustic waves on launch pad components. The DG solver accepts generalized unstructured meshes, enabling efficient application of common mesh generation tools for CHEM and THRUST simulations. The DG solution is coupled with the CFD solution at interface boundaries placed near the CFD acoustic source regions. Both simulations are executed simultaneously with coordinated boundary condition data exchange.

Liever, Peter↗