Engineering Papers⌕ Search

Engineering topics

Candler, Graham V.

Publications and source records attributed to Candler, Graham V..

At least 19 records

Full Facility Shock Frame Simulations of the Electric Arc Shock Tube

Radiative heating computations are performed for high speed lunar return experiments conducted in the Electric Arc Shock Tube (EAST) facility at NASA Ames Research Center. The nonequilibrium radiative transport equations are solved via NASA's in-house radiation code NEQAIR using flow field input from US3D flow solver. The post-shock flow properties for the 10 km/s Earth entry conditions are computed using the stagnation line of a blunt-body and a full facility CFD (Computational Fluid Dynamics) simulation of the EAST shock tube. The shocked gas in the blunt-body flow achieves a thermochemical equilibrium away from the shock front whereas EAST flow exhibits a nonequilibrium behavior due to strong viscous dissipation of the shock by boundary layer. The full-tube flow calculations capture the influence of the boundary layer on the shocked gas state and provide a realistic fluid dynamic input for the radiative predictions. The integrated radiance behind the shock is calculated in NEQAIR for wavelength regimes from Vacuum-UltraViolet (VUV) to InfraRed (IR), which are pertinent to the emission characteristics of high enthalpy shock waves in air. These radiance profiles are validated against corresponding EAST shots. The full-tube simulations successfully predict a sharp radiance peak at the shock front which gets smeared in the test data due to the spatial resolution in the measurements. The full facility based radiance behind the shock shows a slightly better match with the test data in the VUV and Red spectral regions, as compared to that from a blunt-body based predictions. The UV radiance is very similar for both geometries and under-predicts the test behavior. The IR test data matches better with the blunt-body based predictions where the full-tube simulations show a significant over-prediction.

Radiation↗

Estimation of inflow uncertainties in laminar hypersonic double-cone experiments

We propose herein a probabilistic framework for assessing the consistency of an experimental dataset, i.e., whether the stated experimental conditions are consistent with the measurements provided. In case the dataset is inconsistent, our framework allows one to hypothesize and test sources of inconsistencies. This is crucial in model validation efforts. The framework relies on Bayesian inference to estimate experimental settings deemed uncertain, from measurements deemed accurate. The quality of the inferred variables is gauged by its ability to reproduce held-out experimental measurements. We test the correctness of the framework on three double-cone experiments conducted in the CUBRC Inc.'s LENS-I shock tunnel, which have also been numerically simulated successfully. Thereafter, we use the framework to investigate two double-cone experiments (executed in the LENS-XX shock tunnel) which have encountered difficulties when used in model validation exercises. We detect an inconsistency with one of the LENS-XX experiments. In addition, we hypothesize two causes for our inability to simulate LEXS-XX experiments accurately and test them using our framework. We find that there is no single cause that explains all the discrepancies between model predictions and experimental data, but different causes explain different discrepancies, to larger or smaller extent. We end by proposing that uncertainty quantification methods be used more widely to understand experiments and characterize facilities, and we cite three different methods to do so, the third of which we present in this paper.

Thermophysics↗

Chapter 15 - Characterization of Freestream Disturbances in Conventional Hypersonic Wind Tunnels

Prediction of boundary-layer transition is a critical part of the design of hypersonic vehicles because of the large increase in skin-friction drag and surface heating associated with the onset of transition. Testing in conventional (noisy) wind tunnels has been an important means of characterizing and understanding the boundary-layer transition (BLT) behavior of hypersonic vehicles. Because the existing low disturbance, i.e., quiet, facilities operate only at Mach 6, moderate Reynolds numbers, fairly small sizes, and low freestream enthalpy, conventional facilities will continue to be employed for testing and evaluation of hypersonic vehicles, especially for ground testing involving other Mach numbers, higher freestream enthalpies, and larger models. To enable better use of transition data from conventional facilities and more accurate extrapolation of wind-tunnel results to flight, one needs an in-depth knowledge of the broadband disturbance environment in those facilities as well as of the interaction between the freestream disturbances with laminar boundary layers.

42 ENGINEERING↗

Subgrid-scale effects in compressible variable-density decaying turbulence

We present that many turbulent flows are characterized by complex scale interactions and vorticity generation caused by compressibility and variable-density effects. In the large-eddy simulation of variable-density flows, these processes manifest themselves as subgrid-scale (SGS) terms that interact with the resolved-scale flow. This paper studies the effect of the variable-density SGS terms and quantifies their relative importance. We consider the SGS terms appearing in the density-weighted Favre-filtered equations and in the unweighted Reynolds-filtered equations. The conventional form of the Reynolds-filtered momentum equation is complicated by a temporal SGS term; therefore, we derive a new form of the Reynolds-filtered governing equations that does not contain this term and has only double-correlation SGS terms. The new form of the filtered equations has terms that represent the SGS mass flux, pressure-gradient acceleration and velocity-dilatation correlation. To evaluate the dynamical significance of the variable-density SGS effects, we carry out direct numerical simulations of compressible decaying turbulence at a turbulent Mach number of 0.3. Two different initial thermodynamic conditions are investigated: homentropic and a thermally inhomogeneous gas with regions of differing densities. The simulated flow fields are explicitly filtered to evaluate the SGS terms. The importance of the variable-density SGS terms is quantified relative to the SGS specific stress, which is the only SGS term active in incompressible constant-density turbulence. It is found that while the variable-density SGS terms in the homentropic case are negligible, they are dynamically significant in the thermally inhomogeneous flows. Investigation of the variable-density SGS terms is therefore important, not only to develop variable-density closures but also to improve the understanding of scale interactions in variable-density flows.

42 ENGINEERING↗

Characterization of Freestream Disturbances in Conventional Hypersonic Wind Tunnels

While low disturbance ("quiet") hypersonic wind tunnels are believed to provide more reliable extrapolation of boundary layer transition behavior from ground to flight, the presently available quiet facilities are limited to Mach 6, moderate Reynolds numbers, low freestream enthalpy, and subscale models. As a result, only conventional ("noisy") wind tunnels can reproduce both Reynolds numbers and enthalpies of hypersonic flight configurations, and must therefore be used for flight vehicle test and evaluation involving high Mach number, high enthalpy, and larger models. This article outlines the recent progress and achievements in the characterization of tunnel noise that have resulted from the coordinated effort within the AVT-240 specialists group on hypersonic boundary layer transition prediction. New Direct Numerical Simulation (DNS) datasets elucidate the physics of noise generation inside the turbulent nozzle wall boundary layer, characterize the spatiotemporal structure of the freestream noise, and account for the propagation and transfer of the freestream disturbances to a pitot-mounted sensor. The new experimental measurements cover a range of conventional wind tunnels with different sizes and Mach numbers from 6 to 14 and extend the database of freestream fluctuations within the spectral range of boundary layer instability waves over commonly tested models. Prospects for applying the computational and measurement datasets for developing mechanism-based transition prediction models are discussed.

Duan, Lian↗

Scalar conservation and boundedness in simulations of compressible flow

With the proper combination of high-order, low-dissipation numerical methods, physics-based subgrid-scale models, and boundary conditions it is becoming possible to simulate many combustion flows at relevant conditions. However, non-premixed flows are a particular challenge because the thickness of the fuel/oxidizer interface scales inversely with Reynolds number. Sharp interfaces can also be present in the initial or boundary conditions. When higher-order numerical methods are used, there are often aphysical undershoots and overshoots in the scalar variables (e.g.passive scalars, species mass fractions or progress variable). These numerical issues are especially prominent when low-dissipation methods are used, since sharp jumps in flow variables are not always coincident with regions of strong variation in the scalar fields: consequently, special detection mechanisms and dissipative fluxes are needed. Most numerical methods diffuse the interface, resulting in artificial mixing and spurious reactions. In this paper, we propose a numerical method that mitigates this issue. As a result, we present methods for passive and active scalars, and demonstrate their effectiveness with several examples.

97 MATHEMATICS AND COMPUTING↗

Development of the US3D Code for Advanced Compressible and Reacting Flow Simulations

Aerothermodynamics and hypersonic flows involve complex multi-disciplinary physics, including finite-rate gas-phase kinetics, finite-rate internal energy relaxation, gas-surface interactions with finite-rate oxidation and sublimation, transition to turbulence, large-scale unsteadiness, shock-boundary layer interactions, fluid-structure interactions, and thermal protection system ablation and thermal response. Many of the flows have a large range of length and time scales, requiring large computational grids, implicit time integration, and large solution run times. The University of Minnesota NASA US3D code was designed for the simulation of these complex, highly-coupled flows. It has many of the features of the well-established DPLR code, but uses unstructured grids and has many advanced numerical capabilities and physical models for multi-physics problems. The main capabilities of the code are described, the physical modeling approaches are discussed, the different types of numerical flux functions and time integration approaches are outlined, and the parallelization strategy is overviewed. Comparisons between US3D and the NASA DPLR code are presented, and several advanced simulations are presented to illustrate some of novel features of the code.

CFD↗

Parallelization of Unsteady Adaptive Mesh Refinement for Unstructured Navier-Stokes Solvers

This paper explores the implementation of the MPI parallelization in a Navier-Stokes solver using adaptive mesh re nement. Viscous and inviscid test problems are considered for the purpose of benchmarking, as are implicit and explicit time advancement methods. The main test problem for comparison includes e ects from boundary layers and other viscous features and requires a large number of grid points for accurate computation. Ex- perimental validation against double cone experiments in hypersonic ow are shown. The adaptive mesh re nement shows promise for a staple test problem in the hypersonic com- munity. Extension to more advanced techniques for more complicated ows is described.

Schwing, Alan M.↗

Implementation of Implicit Adaptive Mesh Refinement in an Unstructured Finite-Volume Flow Solver

This paper explores the implementation of adaptive mesh refinement in an unstructured, finite-volume solver. Unsteady and steady problems are considered. The effect on the recovery of high-order numerics is explored and the results are favorable. Important to this work is the ability to provide a path for efficient, implicit time advancement. A method using a simple refinement sensor based on undivided differences is discussed and applied to a practical problem: a shock-shock interaction on a hypersonic, inviscid double-wedge. Cases are compared to uniform grids without the use of adapted meshes in order to assess error and computational expense. Discussion of difficulties, advances, and future work prepare this method for additional research. The potential for this method in more complicated flows is described.

Schwing, Alan M.↗

Geometric Effects on the Amplification of First Mode Instability Waves

The effects of geometric changes on the amplification of first mode instability waves in an external supersonic boundary layer were investigated using numerical techniques. Boundary layer stability was analyzed at Mach 6 conditions similar to freestream conditions obtained in quiet ground test facilities so that results obtained in this study may be applied to future test article design to measure first mode instability waves. The DAKOTA optimization software package was used to optimize an axisymmetric geometry to maximize the amplification of the waves at first mode frequencies as computed by the 2D STABL hypersonic boundary layer stability analysis tool. First, geometric parameters such as nose radius, cone half angle, vehicle length, and surface curvature were examined separately to determine the individual effects on the first mode amplification. Finally, all geometric parameters were allowed to vary to produce a shape optimized to maximize the amplification of first mode instability waves while minimizing the amplification of second mode instability waves. Since first mode waves are known to be most unstable in the form of oblique wave, the geometries were optimized using a broad range of wave frequencies as well as a wide range of oblique wave angles to determine the geometry that most amplifies the first mode waves. Since first mode waves are seen most often in flows with low Mach numbers at the edge of the boundary layer, the edge Mach number for each geometry was recorded to determine any relationship between edge Mach number and the stability of first mode waves. Results indicate that an axisymmetric cone with a sharp nose and a slight flare at the aft end under the Mach 6 freestream conditions used here will lower the Mach number at the edge of the boundary layer to less than 4, and the corresponding stability analysis showed maximum first mode N factors of 3.

Kirk, Lindsay C.↗

Hypersonic Navier-Stokes Comparisons to Orbiter Flight Data

During the STS-119 flight of Space Shuttle Discovery, two sets of surface temperature measurements were made. Under the HYTHIRM program3 quantitative thermal images of the windward side of the Orbiter with a were taken. In addition, the Boundary Layer Transition Flight Experiment 4 made thermocouple measurements at discrete locations on the Orbiter wind side. Most of these measurements were made downstream of a surface protuberance designed to trip the boundary layer to turbulent flow. In this paper, we use the US3D computational fluid dynamics code to simulate the Orbiter flow field at conditions corresponding to the STS-119 re-entry. We employ a standard two-temperature, five-species finite-rate model for high-temperature air, and the surface catalysis model of Stewart.1 This work is similar to the analysis of Wood et al . 2 except that we use a different approach for modeling turbulent flow. We use the one-equation Spalart-Allmaras turbulence model8 with compressibility corrections 9 and an approach for tripping the boundary layer at discrete locations. In general, the comparison between the simulations and flight data is remarkably good

Candler, Graham V.↗

Transition Analysis for the HIFiRE-5 Vehicle

The Hypersonic International Flight Research and Experimentation (HIFiRE) 5 flight experiment by Air Force Research Laboratories and Australian Defense Science and Technology Organization is designed to provide in-flight boundary-layer transition data for a canonical 3D configuration at hypersonic Mach numbers. This paper outlines the progress, to date, on boundary layer stability analysis for the HIFiRE-5 flight configuration, as well as for selected test conditions from the wind tunnel experiments supporting the flight test. At flow conditions corresponding to the end of the test window, rather large values of linear amplification factor are predicted for both second mode (N>40) and crossflow (N>20) instabilities, strongly supporting the feasibility of first in-flight measurements of natural transition on a fully three-dimensional hypersonic configuration. Additional results highlight the rich mixture of instability mechanisms relevant to a large segment of the flight trajectory, as well as the effects of angle of attack and yaw angle on the predicted transition fronts for ground facility experiments at Mach 6.

Choudhari, Meelan M.↗

Validation of Afterbody Aeroheating Predictions for Planetary Probes: Status and Future Work

A review of the relevant flight conditions and physical models for planetary probe afterbody aeroheating calculations is given. Readily available sources of afterbody flight data and published attempts to computationally simulate those flights are summarized. A current status of the application of turbulence models to afterbody flows is presented. Finally, recommendations for additional analysis and testing that would reduce our uncertainties in our ability to accurately predict base heating levels are given.

Wright, Michael J.↗

A perspective on computational aerothermodynamics: Current modeling difficulties

Future hypersonic vehicles are going to be designed largely with computational fluid dynamic methods based on appropriate physical models. The question on how much of this design process can be completed with the present state of computational aerothermodynamics is addressed. Some limitations of current models are discussed. It is shown that much more research is required before it will be possible to accurately design a hypersonic vehicle for all of its flight conditions. The quantities that must be computed accurately so that a minimum weight hypersonic vehicle can be designed are discussed. The use of computational fluid dynamics methods coupled with current thermochemical models in order to compute the quantities under specific flow conditions is considered.

Candler, Graham V.↗

Data-parallel lower-upper relaxation method for reacting flows

The implicit lower-upper symmetric Gauss-Seidel (LU-SGS) method of Yoon and Jameson is modified for use on massively parallel computers. The method has been implemented on the Thinking Machines CM-5 and the MasPar MP-1 and MP-2, where large percentages of the theoretical peak floating point performance are obtained. It is shown that the new data-parallel LU relaxation method has better convergence properties than the original method for two different inviscid compressible flow simulations. The convergence is also improved for five-species reacting air computations. The performance of the method on various partitions of the CM-5 and on the MasPar computers is discussed. The new method shows promise for the efficient simulation of very large perfect gas and reacting flows.

Candler, Graham V.↗