Engineering PapersSearch

Engineering topics

Lillard, Randolph P.

Publications and source records attributed to Lillard, Randolph P..

Turbulent Axial Odometer Model

A Galilean-invariant field equation is proposed and tested on standard turbulence model test cases. The field equation provides an additional non-dimensional outer scale which allows the turbulence models to reproduce the axial normal stress increase with Re seen in high Reynolds numbers experiments. The field equation provides a Reynolds number for every point based on the length of turbulent flow upstream of that point in the domain. This outer scale equation can be considered an odometer that gives a length scale conjectured to be related to the large stream wise structures that are seen in turbulent flow and that require run length to develop. A new RANS model using this additional scale is able to match the Reynolds number variation of the normal stresses seen at high Reynolds number. Furthermore, the good attached flow prediction capabilities of current RANS models appears to be attained. Using this scale equation, the entire Reynolds-stress state appears to be predicted correctly, over a large run length Reynolds number range such as experienced in aircraft design.

Olsen, Michael E.

Turbulent Axial Odometer Model

A Galilean-invariant field equation is proposed and tested on standard turbulence model test cases. The field equation provides an additional non-dimensional outer scale which allows the turbulence models to reproduce the axial normal stress increase with Re seen in high Reynolds numbers experiments. The field equation provides a Reynolds number for every point based on the length of turbulent flow upstream of that point in the domain. This outer scale equation can be considered an odometer that gives a length scale conjectured to be related to the large streamwise structures that are seen in turbulent flow and that require run length to develop. A new RANS model using this additional scale is able to match the Reynolds number variation of the normal stresses seen at high Reynolds number. Furthermore, the good attached flow prediction capabilities of current RANS models appears to be attained. Using this scale equation, the entire Reynolds-stress state appears to be predicted correctly, over a large run length Reynolds number range such as experienced in aircraft design.

Olsen, Michael E.

Using Adaptive Mesh Refinement to Study Grid Resolution Effects for Shock-Boundary Layer Interactions

Adaptive Mesh Refinement (AMR) promises a much more computationally efficient means to obtain a discrete approximation to a continuous boundary value problem of a specified accuracy than classic isotropic grid refinement. The AMR capability of OVERFLOW (a computational fluid dynamics (CFD) code) is utilized to provide estimates of the exact analytical solutions to problems of interest to turbulence modeling. Predictions of surface pressure and skin friction, essentially the state of stress at the surface, shows little difference with grids believed to be "grid resolved." Velocity profiles, on the other hand, show marked differences in flows with shocks. The AMR method, as implemented in OVERFLOW 2.2k, appears to provide the ability to produce arbitrarily accurate solutions at a predictable cost much smaller than classic uniform mesh refinement.

Turbulence Modeling

Using Adaptive Mesh Refinement to Study Grid Resolution Effects for Shock/Boundary-Layer Interactions

Adaptive Mesh Refinement (AMR) promises a much more computationally efficient meansto obtain a discrete approximation to a continuous boundary value problem of a specifiedaccuracy than classic isotropic grid refinement. The AMR capability of OVERFLOW is utilizedto provide estimates of the exact analytical solutions to problems of interest to turbulencemodeling. Predictions of surface pressure and skin friction, essentially the state of stress at thesurface, shows little difference with grids believed to be "grid resolved." Velocity profiles, on theother hand, show marked differences in flows with shocks. The AMR method, as implementedin OVERFLOW2.2k, appears to provide the ability to produce arbitrarily accurate solutionsat a predictable cost much smaller than classic uniform mesh refinement.

Adaptive Mesh Refinement

Improvement of Reynolds-Stress and Triple-Product Lag Models

The Reynolds-stress and triple product Lag models were created with a normal stress distribution which was defined by the accepted 4:3:2 distribution of streamwise, spanwise and wall normal stresses, and a ratio of τ(sub w) = 0.3k in the log layer region of high Reynolds number flat plate flow, which implies R11(+) = 4/(9/2)∗.3 ≈ 2.96. More recent measurements show a more complex picture of the log layer region at high Reynolds numbers. The first cut at improving these models along with the direction for future refinements is described. Comparison with recent high Reynolds number data shows areas where further work is needed, but also shows inclusion of the modeled turbulent transport terms improve the prediction where they influence the solution. Additional work is needed to develop a model that better matches experiments, but there is significant improvement in many of the details of the log layer behavior.

Turbulence

Revised Reynolds Stress and Triple Product Models

Revised versions of Lag methodology Reynolds-stress and triple product models are applied to accepted test cases to assess the improvement, or lack thereof, in the prediction capability of the models. The Bachalo-Johnson bump flow is shown as an example for this abstract submission.

fluid dynamics

Improvement of Reynolds-Stress and Triple-Product Lag Models

The Reynolds-stress and triple product Lag models were created with a normal stress distribution which was defined by the accepted 4:3:2 distribution of streamwise, spanwise and wall normal stresses, and a ratio of τ(sub w) = 0.3k in the log layer region of high Reynolds number flat plate flow, which implies R11(+) = 4/(9/2)∗.3 ≈ 2.96. More recent measurements show a more complex picture of the log layer region at high Reynolds numbers. The first cut at improving these models along with the direction for future refinements is described. Comparison with recent high Reynolds number data shows areas where further work is needed, but also shows inclusion of the modeled turbulent transport terms improve the prediction where they influence the solution. Additional work is needed to develop a model that better matches experiments, but there is significant improvement in many of the details of the log layer behavior.

Fluid dynamics

Boundary Layer Transition and Trip Effectiveness on an Apollo Capsule in the JAXA High Enthalpy Shock Tunnel (HIEST) Facility

Computational assessments were performed to size boundary layer trips for a scaled Apollo capsule model in the High Enthalpy Shock Tunnel (HIEST) facility at the JAXA Kakuda Space Center in Japan. For stagnation conditions between 2 MJ/kg and 20 MJ/kg and between 10 MPa and 60 MPa, the appropriate trips were determined to be between 0.2 mm and 1.3 mm high, which provided kappa/delta values on the heatshield from 0.15 to 2.25. The tripped configuration consisted of an insert with a series of diamond shaped trips along the heatshield downstream of the stagnation point. Surface heat flux measurements were obtained on a capsule with a 250 mm diameter, 6.4% scale model, and pressure measurements were taken at axial stations along the nozzle walls. At low enthalpy conditions, the computational predictions agree favorably to the test data along the heatshield centerline. However, agreement becomes less favorable as the enthalpy increases conditions. The measured surface heat flux on the heatshield from the HIEST facility was under-predicted by the computations in these cases. Both smooth and tripped configurations were tested for comparison, and a post-test computational analysis showed that kappa/delta values based on the as-measured stagnation conditions ranged between 0.5 and 1.2. Tripped configurations for both 0.6 mm and 0.8 mm trip heights were able to effectively trip the flow to fully turbulent for a range of freestream conditions.

Kirk, Lindsay C.

Technology Needs to Support Future Mars Exploration

The Mars Program Planning Group (MPPG) under the direction of Dr. Orlando Figueroa, was chartered to develop options for a program-level architecture for robotic exploration of Mars consistent with the objective to send humans to Mars in the 2030's. Scientific pathways were defined for future exploration, and multiple architectural options were developed that meet current science goals and support the future human exploration objectives. Integral to the process was the identification of critical technologies which enable the future scientific and human exploration goals. This paper describes the process for technology capabilities identification and examines the critical capability needs identified in the MPPG process. Several critical enabling technologies that have been identified to support the robotic exploration goals and with potential feedforward application to human exploration goals. Potential roadmaps for the development and validation of these technologies are discussed, including options for subscale technology demonstrations of future human exploration technologies on robotic missions.

Mars Program Planning Group (MPPG)

Aeroheating Measurement of Apollo Shaped Capsule with Boundary Layer Trip in the Free-piston Shock Tunnel HIEST

An aeroheating measurement test campaign of an Apollo capsule model with laminar and turbulent boundary layer was performed in the free-piston shock tunnel HIEST at JAXA Kakuda Space Center. A 250mm-diameter 6.4%-scaled Apollo CM capsule model made of SUS-304 stainless steel was applied in this study. To measure heat flux distribution, the model was equipped with 88 miniature co-axial Chromel-Constantan thermocouples on the heat shield surface of the model. In order to promote boundary layer transition, a boundary layer trip insert with 13 "pizza-box" isolated roughness elements, which have 1.27mm square, were placed at 17mm below of the model geometric center. Three boundary layer trip inserts with roughness height of k=0.3mm, 0.6mm and 0.8mm were used to identify the appropriate height to induce transition. Heat flux records with or without roughness elements were obtained for model angles of attack 28º under stagnation enthalpy between H(sub 0)=3.5MJ/kg to 21MJ/kg and stagnation pressure between P(sub 0)=14MPa to 60MPa. Under the condition above, Reynolds number based on the model diameter was varied from 0.2 to 1.3 million. With roughness elements, boundary layer became fully turbulent less than H(sub 0)=9MJ/kg condition. However, boundary layer was still laminar over H(sub 0)=13MJ/kg condition even with the highest roughness elements. An additional experiment was also performed to correct unexpected heat flux augmentation observed over H(sub 0)=9MJ/kg condition.

Hideyuki, TANNO

Turbulence Modeling for Shock Wave/Turbulent Boundary Layer Interactions

Accurate aerodynamic computational predictions are essential for the safety of space vehicles, but these computations are of limited accuracy when large pressure gradients are present in the flow. The goal of the current project is to improve the state of compressible turbulence modeling for high speed flows with shock wave / turbulent boundary layer interactions (SWTBLI). Emphasis will be placed on models that can accurately predict the separated region caused by the SWTBLI. These flows are classified as nonequilibrium boundary layers because of the very large and variable adverse pressure gradients caused by the shock waves. The lag model was designed to model these nonequilibrium flows by incorporating history effects. Standard one- and two-equation models (Spalart Allmaras and SST) and the lag model will be run and compared to a new lag model. This new model, the Reynolds stress tensor lag model (lagRST), will be assessed against multiple wind tunnel tests and correlations. The basis of the lag and lagRST models are to preserve the accuracy of the standard turbulence models in equilibrium turbulence, when the Reynolds stresses are linearly related to the mean strain rates, but create a lag between mean strain rate effects and turbulence when nonequilibrium effects become important, such as in large pressure gradients. The affect this lag has on the results for SWBLI and massively separated flows will be determined. These computations will be done with a modified version of the OVERFLOW code. This code solves the RANS equations on overset grids. It was used for this study for its ability to input very complex geometries into the flow solver, such as the Space Shuttle in the full stack configuration. The model was successfully implemented within two versions of the OVERFLOW code. Results show a substantial improvement over the baseline models for transonic separated flows. The results are mixed for the SWBLI assessed. Separation predictions are not as good as the baseline models, but the over prediction of the peak heat flux downstream of the reattachment shock that plagues many models is reduced.

Lillard, Randolph P.

The lagRST Model: A Turbulence Model for Non-Equilibrium Flows

This study presents a new class of turbulence model designed for wall bounded, high Reynolds number flows with separation. The model addresses deficiencies seen in the modeling of nonequilibrium turbulent flows. These flows generally have variable adverse pressure gradients which cause the turbulent quantities to react at a finite rate to changes in the mean flow quantities. This "lag" in the response of the turbulent quantities can t be modeled by most standard turbulence models, which are designed to model equilibrium turbulent boundary layers. The model presented uses a standard 2-equation model as the baseline for turbulent equilibrium calculations, but adds transport equations to account directly for non-equilibrium effects in the Reynolds Stress Tensor (RST) that are seen in large pressure gradients involving shock waves and separation. Comparisons are made to several standard turbulence modeling validation cases, including an incompressible boundary layer (both neutral and adverse pressure gradients), an incompressible mixing layer and a transonic bump flow. In addition, a hypersonic Shock Wave Turbulent Boundary Layer Interaction with separation is assessed along with a transonic capsule flow. Results show a substantial improvement over the baseline models for transonic separated flows. The results are mixed for the SWTBLI flows assessed. Separation predictions are not as good as the baseline models, but the over prediction of the peak heat flux downstream of the reattachment shock that plagues many models is reduced.

Lillard, Randolph P.

Aerothermal Testing for Project Orion Crew Exploration Vehicle

The Project Orion Crew Exploration Vehicle aerothermodynamic experimentation strategy, as it relates to flight database development, is reviewed. Experimental data has been obtained to both validate the computational predictions utilized as part of the database and support the development of engineering models for issues not adequately addressed with computations. An outline is provided of the working groups formed to address the key deficiencies in data and knowledge for blunt reentry vehicles. The facilities utilized to address these deficiencies are reviewed, along with some of the important results obtained thus far. For smooth wall comparisons of computational convective heating predictions against experimental data from several facilities, confidence was gained with the use of algebraic turbulence model solutions as part of the database. For cavities and protuberances, experimental data is being used for screening various designs, plus providing support to the development of engineering models. With the reaction-control system testing, experimental data were acquired on the surface in combination with off-body flow visualization of the jet plumes and interactions. These results are being compared against predictions for improved understanding of aftbody thermal environments and uncertainties.

Berry, Scott A.

Aeroheating Testing and Predictions for Project Orion CEV at Turbulent Conditions

An investigation of the aeroheating environment of the Project Orion Crew Exploration Vehicle was performed in the Arnold Engineering Development Center Hypervelocity Wind Tunnel No. 9 Mach 8 and Mach 10 nozzles and in the NASA Langley Research Center 20 - Inch Mach 6 Air Tunnel. Heating data were obtained using a thermocouple-instrumented approx.0.035-scale model (0.1778-m/7-inch diameter) of the flight vehicle. Runs were performed in the Tunnel 9 Mach 10 nozzle at free stream unit Reynolds numbers of 1x10(exp 6)/ft to 20x10(exp 6)/ft, in the Tunnel 9 Mach 8 nozzle at free stream unit Reynolds numbers of 8 x 10(exp 6)/ft to 48x10(exp 6)/ft, and in the 20-Inch Mach 6 Air Tunnel at free stream unit Reynolds numbers of 1x10(exp 6)/ft to 7x10(exp 6)/ft. In both facilities, enthalpy levels were low and the test gas (N2 in Tunnel 9 and air in the 20-Inch Mach 6) behaved as a perfect-gas. These test conditions produced laminar, transitional and turbulent data in the Tunnel 9 Mach 10 nozzle, transitional and turbulent data in the Tunnel 9 Mach 8 nozzle, and laminar and transitional data in the 20- Inch Mach 6 Air Tunnel. Laminar and turbulent predictions were generated for all wind tunnel test conditions and comparisons were performed with the experimental data to help define the accuracy of computational method. In general, it was found that both laminar data and predictions, and turbulent data and predictions, agreed to within less than the estimated 12% experimental uncertainty estimate. Laminar heating distributions from all three data sets were shown to correlate well and demonstrated Reynolds numbers independence when expressed in terms of the Stanton number based on adiabatic wall-recovery enthalpy. Transition onset locations on the leeside centerline were determined from the data and correlated in terms of boundary-layer parameters. Finally turbulent heating augmentation ratios were determined for several body-point locations and correlated in terms of the boundary-layer momentum Reynolds number.

Hollis, Brian R.

Experimental Investigation of Project Orion Crew Exploration Vehicle Aeroheating in AEDC Tunnel 9

An investigation of the aeroheating environment of the Project Orion Crew Entry Vehicle has been performed in the Arnold Engineering Development Center Tunnel 9. The goals of this test were to measure turbulent heating augmentation levels on the heat shield and to obtain high-fidelity heating data for assessment of computational fluid dynamics methods. Laminar and turbulent predictions were generated for all wind tunnel test conditions and comparisons were performed with the data for the purpose of helping to define uncertainty margins for the computational method. Data from both the wind tunnel test and the computational study are presented herein.

Hollis, Brian R.

Aeroheating Test of CEV Entry Vehicle at Turbulent Conditions

An investigation of the aeroheating environment of the Project Orion Crew Entry Vehicle has been performed in the Arnold Engineering Development Center Tunnel 9. Data were measured on a approx. 3.5% scale model (0.1778m/7-inch diam.) of the vehicle using coaxial thermocouples in the Mach 8 and Mach 10 nozzles of Tunnel 9. Runs were performed at free stream Reynolds numbers of 1 106/ft to 20 10(exp 6)/ft in the Mach 10 nozzle and 8 10(exp 6)/ft to 48 10(exp 6)/ft in the Mach 8 nozzle. The test gas in Tunnel 9 is pure N2, which at these operating conditions remains un-dissociated and may be treated as a perfect gas. At these conditions, laminar, transitional, and turbulent flow was produced on the model at Mach 10, and transitional and turbulent conditions were produced on the model at Mach 8. The majority of runs were made on a clean, smooth-surface model configuration and a limited number of runs were made in which inserts with varying boundary-layer trips configurations were used to force the occurrence of transition. Laminar and turbulent predictions were generated for all wind tunnel test conditions and comparisons were performed with the data for the purpose of helping to define uncertainty margins for the computational method. Data from both the wind tunnel test and the computations are presented herein. Figure 1 shows a schematic of the thermocouple locations on the model and figures 2 and 3 show a photo and schematic of the AEDC Hypervelocity Tunnel 9. Figure 4 shows a typical grid used in the computations. From the comparisons shown in figures 5 through 8 it was concluded that for perfect-gas conditions, the computations could predict either fully-laminar or full-turbulent flow to within +/-10% of the experimental data. The experimental data showed that transition began on the leeside of the heatshield at a free stream Reynolds number of 9 10(exp 6)/ft in the Mach 10 nozzle and fully-developed turbulent flow was produced at 20 10(exp 6)/ft. In the Mach 8 nozzle, transition on the leeside of the heat-shield was observed for all test conditions, and full-developed turbulent flow occurred at a free stream Reynolds number of 18 10(exp 6)/ft. On the aftbody of the vehicle no evidence of turbulence was detected at Mach 10 conditions, and at Mach 8 conditions, transition appeared to begin on the windside of the aftbody at free stream Reynolds number of 18x10(exp 6)/ft with fully-developed turbulent flow occurring only at the highest test condition of 48x10(exp 6)/ft.

Hollis, Brian R.

The Lag Model Applied to High Speed Flows

The Lag model has shown great promise in prediction of low speed and transonic separations. The predictions of the model, along with other models (Spalart-Allmaras and Menter SST) are assessed for various high speed flowfields. In addition to skin friction and separation predictions, the prediction of heat transfer are compared among these models, and some fundamental building block flowfields, are investigated.

Olsen, Michael E.

Laminar Heating Validation of the OVERFLOW Code

OVERFLOW, a structured finite difference code, was applied to the solution of hypersonic laminar flow over several configurations assuming perfect gas chemistry. By testing OVERFLOW's capabilities over several configurations encompassing a variety of flow physics a validated laminar heating was produced. Configurations tested were a flat plate at 0 degrees incidence, a sphere, a compression ramp, and the X-38 re-entry vehicle. This variety of test cases shows the ability of the code to predict boundary layer flow, stagnation heating, laminar separation with re-attachment heating, and complex flow over a three-dimensional body. In addition, grid resolutions studies were done to give recommendations for the correct number of off-body points to be applied to generic problems and for wall-spacing values to capture heat transfer and skin friction. Numerical results show good comparison to the test data for all the configurations.

Lillard, Randolph P.