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

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

Rocket-Based Combined Cycle Engine Technology Development: Inlet CFD Validation and Application

A CFD methodology has been developed for inlet analyses of Rocket-Based Combined Cycle (RBCC) Engines. A full Navier-Stokes analysis code, NPARC, was used in conjunction with pre- and post-processing tools to obtain a complete description of the flow field and integrated inlet performance. This methodology was developed and validated using results from a subscale test of the inlet to a RBCC 'Strut-Jet' engine performed in the NASA Lewis 1 x 1 ft. supersonic wind tunnel. Results obtained from this study include analyses at flight Mach numbers of 5 and 6 for super-critical operating conditions. These results showed excellent agreement with experimental data. The analysis tools were also used to obtain pre-test performance and operability predictions for the RBCC demonstrator engine planned for testing in the NASA Lewis Hypersonic Test Facility. This analysis calculated the baseline fuel-off internal force of the engine which is needed to determine the net thrust with fuel on

A Monte Carlo study of the three-dimensional Coulomb frustrated Ising ferromagnet

We have investigated by Monte-Carlo simulation the phase diagram of a three-dimensional Ising model with nearest-neighbor ferromagnetic interactions and small, but long-range (Coulombic) antiferromagnetic interactions. We have developed an efficient cluster algorithm and used different lattice sizes and geometries, which allows us to obtain the main characteristics of the temperature-frustration phase diagram. Our finite-size scaling analysis confirms that the melting of the lamellar phases into the paramgnetic phase is driven first-order by the fluctuations. Transitions between ordered phases with different modulation patterns is observed in some regions of the diagram, in agreement with a recent mean-field analysis.Comment: 14 pages, 10 figures, submitted to Phys. Rev.

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

Modeling of Turbulent Free Shear Flows

The modeling of turbulent free shear flows is crucial to the simulation of many aerospace applications, yet often receives less attention than the modeling of wall boundary layers. Thus, while turbulence model development in general has proceeded very slowly in the past twenty years, progress for free shear flows has been even more so. This paper highlights some of the fundamental issues in modeling free shear flows for propulsion applications, presents a review of past modeling efforts, and identifies areas where further research is needed. Among the topics discussed are differences between planar and axisymmetric flows, development versus self-similar regions, the effect of compressibility and the evolution of compressibility corrections, the effect of temperature on jets, and the significance of turbulent Prandtl and Schmidt numbers for reacting shear flows. Large eddy simulation greatly reduces the amount of empiricism in the physical modeling, but is sensitive to a number of numerical issues. This paper includes an overview of the importance of numerical scheme, mesh resolution, boundary treatment, sub-grid modeling, and filtering in conducting a successful simulation

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 highorder 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

Inlet Development for a Rocket Based Combined Cycle, Single Stage to Orbit Vehicle Using Computational Fluid Dynamics

Design and analysis of the inlet for a rocket based combined cycle engine is discussed. Computational fluid dynamics was used in both the design and subsequent analysis. Reynolds averaged Navier-Stokes simulations were performed using both perfect gas and real gas assumptions. An inlet design that operates over the required Mach number range from 0 to 12 was produced. Performance data for cycle analysis was post processed using a stream thrust averaging technique. A detailed performance database for cycle analysis is presented. The effect ot vehicle forebody compression on air capture is also examined

A Comparison of Three Navier-Stokes Solvers for Exhaust Nozzle Flowfields

A comparison of the NPARC, PAB, and WIND (previously known as NASTD) Navier-Stokes solvers is made for two flow cases with turbulent mixing as the dominant flow characteristic, a two-dimensional ejector nozzle and a Mach 1.5 elliptic jet. The objective of the work is to determine if comparable predictions of nozzle flows can be obtained from different Navier-Stokes codes employed in a multiple site research program. A single computational grid was constructed for each of the two flows and used for all of the Navier-Stokes solvers. In addition, similar k-e based turbulence models were employed in each code, and boundary conditions were specified as similarly as possible across the codes. Comparisons of mass flow rates, velocity profiles, and turbulence model quantities are made between the computations and experimental data. The computational cost of obtaining converged solutions with each of the codes is also documented. Results indicate that all of the codes provided similar predictions for the two nozzle flows. Agreement of the Navier-Stokes calculations with experimental data was good for the ejector nozzle. However, for the Mach 1.5 elliptic jet, the calculations were unable to accurately capture the development of the three dimensional elliptic mixing layer

Geometry parametrization and aerodynamic characteristics of axisymmetric afterbodies

A key aspect of the preliminary design process for a new generation combat aircraft is the prediction of afterbody aerodynamic drag. Current prediction methods for preliminary design are constrained in terms of number of independent geometric degrees of freedom that can be studied due to the classic circular arc or conical afterbody geometry parametrization. In addition, the amount of data available for the construction of the reliable performance correlations is too sparse. This paper presents a methodology for the generation of aerodynamic performance maps for transonic axisymmetric afterbody and exhaust systems. It uses a novel parametric geometry definition along with a compressible flow solver to conduct an extensive design space exploration. The proposed geometry parametrization is based on the Class Shape Transformation method and it enables the assessment of the aerodynamic performance of a wider range of afterbodies at the expense of one additional geometric degree of freedom. Relative to the conventional approach, this enables the exploration of a wider design space and the construction of more complete aerodynamic performance maps. This research quantifies the impact of a number of geometric degrees of freedom on the aerodynamic performance of transonic afterbody and exhaust systems at different operating conditions