Experimental and analytical studies of flow through a ventral and axial exhaust nozzle system for STOVL aircraft

Flow through a combined ventral and axial exhaust nozzle system was studied experimentally and analytically. The work is part of an ongoing propulsion technology effort at NASA Lewis Research Center for short takeoff, vertical landing (STOVL) aircraft. The experimental investigation was done on the NASA Lewis Powered Lift Facility. The experiment consisted of performance testing over a range of tailpipe pressure ratios from 1 to 3.2 and flow visualization. The analytical investigation consisted of modeling the same configuration and solving for the flow using the PARC3D computational fluid dynamics program. The comparison of experimental and analytical results was very good. The ventral nozzle performance coefficients obtained from both the experimental and analytical studies agreed within 1.2 percent. The net horizontal thrust of the nozzle system contained a significant reverse thrust component created by the flow overturning in the ventral duct. This component resulted in a low net horizontal thrust coefficient. The experimental and analytical studies showed very good agreement in the internal flow patterns

Full Navier-Stokes analysis of a two-dimensional mixer/ejector nozzle for noise suppression

A three-dimensional full Navier-Stokes (FNS) analysis was performed on a mixer/ejector nozzle designed to reduce the jet noise created at takeoff by a future supersonic transport. The PARC3D computational fluid dynamics (CFD) code was used to study the flow field of the nozzle. The grid that was used in the analysis consisted of approximately 900,000 node points contained in eight grid blocks. Two nozzle configurations were studied: a constant area mixing section and a diverging mixing section. Data are presented for predictions of pressure, velocity, and total temperature distributions and for evaluations of internal performance and mixing effectiveness. The analysis provided good insight into the behavior of the flow

A High-Resolution Capability for Large-Eddy Simulation of Jet Flows

A large-eddy simulation (LES) code that utilizes high-resolution numerical schemes is described and applied to a compressible jet flow. The code is written in a general manner such that the accuracy/resolution of the simulation can be selected by the user. Time discretization is performed using a family of low-dispersion Runge-Kutta schemes, selectable from first- to fourth-order. Spatial discretization is performed using central differencing schemes. Both standard schemes, second- to twelfth-order (3 to 13 point stencils) and Dispersion Relation Preserving schemes from 7 to 13 point stencils are available. The code is written in Fortran 90 and uses hybrid MPI/OpenMP parallelization. The code is applied to the simulation of a Mach 0.9 jet flow. Four-stage third-order Runge-Kutta time stepping and the 13 point DRP spatial discretization scheme of Bogey and Bailly are used. The high resolution numerics used allows for the use of relatively sparse grids. Three levels of grid resolution are examined, 3.5, 6.5, and 9.2 million points. Mean flow, first-order turbulent statistics and turbulent spectra are reported. Good agreement with experimental data for mean flow and first-order turbulent statistics is shown

Evaluation of Industry Standard Turbulence Models on an Axisymmetric Supersonic Compression Corner

Reynolds-averaged Navier-Stokes computations of a shock-wave/boundary-layer interaction (SWBLI) created by a Mach 2.85 flow over an axisymmetric 30-degree compression corner were carried out. The objectives were to evaluate four turbulence models commonly used in industry, for SWBLIs, and to evaluate the suitability of this test case for use in further turbulence model benchmarking. The Spalart-Allmaras model, Menter's Baseline and Shear Stress Transport models, and a low-Reynolds number k- model were evaluated. Results indicate that the models do not accurately predict the separation location; with the SST model predicting the separation onset too early and the other models predicting the onset too late. Overall the Spalart-Allmaras model did the best job in matching the experimental data. However there is significant room for improvement, most notably in the prediction of the turbulent shear stress. Density data showed that the simulations did not accurately predict the thermal boundary layer upstream of the SWBLI. The effect of turbulent Prandtl number and wall temperature were studied in an attempt to improve this prediction and understand their effects on the interaction. The data showed that both parameters can significantly affect the separation size and location, but did not improve the agreement with the experiment. This case proved challenging to compute and should provide a good test for future turbulence modeling work

Prediction of Turbulent Temperature Fluctuations in Hot Jets

Large-eddy simulations were used to investigate turbulent temperature fluctuations and turbulent heat flux in hot jets. A high-resolution finite-difference Navier-Stokes solver, WRLES, was used to compute the flow from a 2-inch round nozzle. Several different flow conditions, consisting of different jet Mach numbers and temperature ratios, were examined. Predictions of mean and fluctuating velocities were compared to previously obtained particle image velocimetry data. Predictions of mean and fluctuating temperature were compared to new data obtained using Raman spectroscopy. Based on the good agreement with experimental data for the individual quantities, the combined quantity turbulent heat flux was examined

Solutions of the Taylor-Green Vortex Problem Using High-Resolution Explicit Finite Difference Methods

A computational fluid dynamics code that solves the compressible Navier-Stokes equations was applied to the Taylor-Green vortex problem to examine the code s ability to accurately simulate the vortex decay and subsequent turbulence. The code, WRLES (Wave Resolving Large-Eddy Simulation), uses explicit central-differencing to compute the spatial derivatives and explicit Low Dispersion Runge-Kutta methods for the temporal discretization. The flow was first studied and characterized using Bogey & Bailley s 13-point dispersion relation preserving (DRP) scheme. The kinetic energy dissipation rate, computed both directly and from the enstrophy field, vorticity contours, and the energy spectra are examined. Results are in excellent agreement with a reference solution obtained using a spectral method and provide insight into computations of turbulent flows. In addition the following studies were performed: a comparison of 4th-, 8th-, 12th- and DRP spatial differencing schemes, the effect of the solution filtering on the results, the effect of large-eddy simulation sub-grid scale models, and the effect of high-order discretization of the viscous terms

Three-dimensional compressible turbulent computations for a diffusing S-duct

The purpose of the present study was to evaluate the capability of the computational fluid dynamics computer program PARC3D to model flow in a typical diffusing subsonic S-duct, with strong secondary flows. This evaluation is needed to provide confidence in the analysis of aircraft inlets, which have similar geometries. The performance predictions include total pressure profiles, static pressures, velocity profiles, boundary layer data, and skin friction data. Flow in the S-duct is subsonic, and the boundary layers are assumed to be turbulent. The results for both H and O grid solutions, are compared with existing test data

Crises and bail-outs of banks and countries: Interconnections, analogies, differences

During the last decade or so a number of banking crises has occurred in both industrial and developing countries. Bank unsoundness is the focus of deep concerns for its possible interactions with macroeconomic instability. The potential size and severity of sovereign debt crises in integrated financial markets were highlighted by the outbreak of the Mexican crisis at the end of 1994. Subsequently, in a number of international bodies policy makers have undertaken to consider various avenues and arrangements to prevent, anticipate and resolve sovereign debt crises. In the second half of 1997, East Asia's turmoil testified to the interdependence between currency devaluations, stockmarket crashes, and banks' insolvencies. After a general discussion of financial crises and bank runs, the paper discusses the prototype case of banks for which both established economic doctrine and the practical experience of policy makers and regulators offer clearer indications. Then it turns to countries exploring interactions, analogies and differences with the case of banks. Fundamental differences remain between the cases of financial intermediaries and nations. Ultimately the decision of a sovereign state to default or suspend its debt-service payments is largely a voluntary one and the safeguards against moral hazard built into domestic bankruptcy codes cannot be applied to it. Important consequences follow in terms of international co-ordination and the function of the "lender of last resort"

Three-dimensional compressible turbulent computations for a nondiffusing S-duct

The PARC3D code was used to compute the compressible turbulent flow within a three dimensional, nondiffusing S-duct. A frame of reference is provided for future computational fluid dynamics studies of internal flows with strong secondary flows and provides an understanding of the performance characteristics of a typical S-duct with attached flow. The predicted results, obtained with both H- and O-grids, are compared with the experimental wall pressure, static and total pressure fields, and velocity vectors. Additionally, computed boundary layer thickness, velocity profiles in wall coordinates, and skin friction values are presented

Comparison of High-Order and Low-Order Methods for Large-Eddy Simulation of a Compressible Shear Layer

The objective of this work is to compare a high-order solver with a low-order solver for performing large-eddy simulations (LES) of a compressible mixing layer. The high-order method is the Wave-Resolving LES (WRLES) solver employing a Dispersion Relation Preserving (DRP) scheme. The low-order solver is the Wind-US code, which employs the second-order Roe Physical scheme. Both solvers are used to perform LES of the turbulent mixing between two supersonic streams at a convective Mach number of 0.46. The high-order and low-order methods are evaluated at two different levels of grid resolution. For a fine grid resolution, the low-order method produces a very similar solution to the high-order method. At this fine resolution the effects of numerical scheme, subgrid scale modeling, and filtering were found to be negligible. Both methods predict turbulent stresses that are in reasonable agreement with experimental data. However, when the grid resolution is coarsened, the difference between the two solvers becomes apparent. The low-order method deviates from experimental results when the resolution is no longer adequate. The high-order DRP solution shows minimal grid dependence. The effects of subgrid scale modeling and spatial filtering were found to be negligible at both resolutions. For the high-order solver on the fine mesh, a parametric study of the spanwise width was conducted to determine its effect on solution accuracy. An insufficient spanwise width was found to impose an artificial spanwise mode and limit the resolved spanwise modes. We estimate that the spanwise depth needs to be 2.5 times larger than the largest coherent structures to capture the largest spanwise mode and accurately predict turbulent mixing