SEARCH · Engineering Papers
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.
The role of vibration-dissociation coupling in flows representative of AOTV's with afterbodies is discussed. A physical model designed to reflect the behavior of anharmonic diatomic molecules at high temperatures was developed and incorporated into the direct simulation Monte Carlo method of Bird (1976). The model leads to reduced dissociation and increased convective heating for flows with finite catalytic surfaces.
The flow behavior in a vaneless radial turbine scroll was examined experimentally. The data was obtained using the slant sensor technique of hot film anemometry. This method used the unsymmetric heat transfer characteristics of a constant temperature hot film sensor to detect the flow direction and magnitude. This was achieved by obtaining a velocity vector measurement at three sensor positions with respect to the flow. The true magnitude and direction of the velocity vector was then found using these values and a Newton-Raphson numerical technique. The through flow and secondary flow velocity components are measured at various points in three scroll sections.
The thermal history of the earth is investigated using a parameterized model of mantle convection, that includes the effects of volatile exchange between the mantle and the surface reservoir and the softening of the mantle by the dissolved volatiles. The mantle degassing rate is taken to be directly proportional to the rate of seafloor spreading which depends on the mantle heat flow. It is shown that the dependence of the mantle viscosity on the volatile content has important effects on the thermal evolution of planetary interiors and the evolution of planetary atmospheres. Degassing is compensated by an increase in temperature, while regassing is compensated by a decrease in temperature. Reasonable degassing scenarios can account for an early rapid formation of the earth's atmosphere inferred from noble gas abundances.
The Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer (SPHEREx) instrument has been proposed as NASA Small Explorer (SMEX) and Medium Explorer (MIDEX) missions designed to perform an all-sky near-infrared survey. The primary science objective of the SPHEREx instrument is the mapping and cataloging of data on galaxy red shifts and absorption spectra. The SPHEREx instrument utilizes a linear variable band pass filter to provide discrete spectral coverage of the 0.75 to 5.0 μm range. The nominal on-orbit mission lifetime for the instrument is two years with continuous primary science imaging. The SPHEREx instrument will be placed in a near-circular LEO sun-synchronous terminator orbit with a mean altitude of 600 km, and solar beta angles ranging from 59° to 90°. A passive V-groove radiator cooling system was selected for the SPHEREx instrument. This design was chosen for its simplicity and lack of vibration. It includes five radiator stages with the coldest two stages providing temperatures below 80K and 55K for the 2.5μm and 5.3μm detectors respectively. A one quarter scale thermal prototype was constructed to characterize the performance of the SPHEREx thermal control subsystem. This prototype was directly derived and scaled from the current flight design, and accurately replicates all heat flow paths at scale. The prototype was tested in two bounding hot flight-like environments in the Cryogenic Systems Engineering, Advanced Thermal Technology Lab at NASA’s Jet Propulsion Laboratory. Following testing, a previously developed one quarter scale thermal model was correlated to the two conditions seen in test. The resultant correlated model is able to predict the prototype test temperatures to within 1K for all critical radiator stages. An overview of the thermal control design approach, the test configuration, and test and correlation results are presented.
The Mars Sample Return Earth Entry System is a mission concept which would be the first NASA entry vehicle to utilize a woven thermal protection system as well as a 52.5 degree sphere-cone forebody geometry. Due to its novel design, historic experimental data sets were insufficient to validate the models used to characterize the expected aerothermal environment. A wind tunnel test campaign was funded by the Mars Sample Return Earth Entry System project in the NASA Langley 20-Inch Mach 6 Air Tunnel in order to obtain validation data. The December 2023 test entry of this campaign sought to obtain thermographic data to capture the impact of supersonic flow near woven surface roughness elements on convective heating augmentation, which was the largest remaining uncertainty in aerothermal design of the Earth Entry System. Data obtained were in family with prior results obtained on identical woven patterns without supersonic flow present at roughness elements, and it was found that this presence of supersonic flow did not have a noticeable effect on surface convective heating augmentation. Data were directly leveraged to support the NASA Mars Sample Return flight program and this data set is valuable for validating computational solutions on roughness-resolved grids.
The Mars Sample Return Earth Entry System is a mission concept which would be the first NASA entry vehicle to utilize a woven thermal protection system as well as a 52.5 degree sphere-cone forebody geometry. Due to its novel design, historic experimental data sets were insufficient to validate the models used to characterize the expected aerothermal environment. A wind tunnel test campaign was funded by the Mars Sample Return Earth Entry System project in the NASA Langley 20-Inch Mach 6 Air Tunnel in order to obtain validation data. The December 2023 test entry of this campaign sought to obtain thermographic data to capture the impact of supersonic flow near woven surface roughness elements on convective heating augmentation, which was the largest remaining uncertainty in aerothermal design of the Earth Entry System. Data obtained were in family with prior results obtained on identical woven patterns without supersonic flow present at roughness elements, and it was found that this presence of supersonic flow did not have a noticeable effect on surface convective heating augmentation. Data were directly leveraged to support the NASA Mars Sample Return flight program and this data set is valuable for validating computational solutions on roughness-resolved grids.
A computational study is presented for the flow inside an oscillatory cavity. The numerical scheme employs a semiimplicit, time-splitting method to integrate the two-dimensional full Navier-Stokes equations satisfying continuity to machine accuracy. The efficient use of direct solvers for the uncoupled momentum and pressure equations is demonstrated. The oscillatory cavity flow is studied considering the effects of heat transfer, Reynolds number and oscillatory Stokes number.
A vertical flow loop was designed and assembled to determine the local (circumferential and axial) and mean wall temperature distributions for single-phase and two-phase (subcooled and saturated) downward flow in both uniformly-heated and single-side heated vertical channels. Freon-11 was used as the working fluid in order to directly relate and compare the results with a previous experimental campaign which employed this same working fluid. For a given steady-state experiment, the following parameters were held constant: (1) exit pressure, (2) inlet temperature, and (3) mass velocity. For a given configuration of the 2.2 m long cylindrical channel test section, which had a 1.2 m long heated section, the applied heat rate was varied from zero through successive quasi-steady states to a level which corresponded to localized film boiling in the test section. The measurements showed that the boiling curve changes significantly at higher mass velocities with respect to both the circumferential and axial directions. The slope of the boiling curve changes in a non-monotonic fashion with respect to the circumferential directions. The slope of the boiling curve changes in a non-monotonic fashion with respect to the circumferential direction. The measurements point to the existence of a dry-out phenomenon occurring at multiple levels of the applied heat for the single-side heated channel. In comparing the heat transfer for horizontal channel flow with a vertically downward flow, the results show that significantly lower heat transfer occurs in the horizontal flow. However, this trend reverses as both the Reynolds number and the applied heat rate increase. Both the Liu-Winterton and Shah correlations were compared with the experimental data. The Shah correlation predicted the uniformly heated tube data better. When a thermal hydraulic diameter approach was used for the single-side heated case, the data at upstream locations for Z/L less than 0.5 was bounded above by the Liu-Winterton correlation and below by the Shah correlation. At Z/L = 0.5, the Shah correlation bounded the data; and for Z/L greater than 0.5, both correlations overpredicted the data with the Shah correlation being closest to the data. The present results indicate that additional correlational development is needed. In addressing some of the advanced space thermal management objectives concerning accommodating high heat fluxes in non-uniformly heated systems, a large battery of experiments 88 have been completed where local two-dimensional wall temperature variations were measured for both single-phase and two-phase flow in a single-side heated circular tube. As noted above, the results show significant axial and circumferential variations. Accurately accounting for such variations can result in optimized future advanced space, enhanced (high heat flux) thermal management systems.
Spallation is the mechanical removal of small chunks of material gets removed typically due to high shear conditions of the flow field. This reduces the ability of the thermal protection system (TPS) material to protect the spacecraft as well as cause turbulence in the flow causing higher heating rates. In this work, we focus on the material removal through ablation and high shear flow within the boundary layer region of woven TPS material. Woven TPS (WTPS) material is the latest class of material developed by NASA, to be used within the next generation of space flights. They are complex interlocked weaves designed to create a rigid structure that is highly resistant to heat and can be easily designed and tailored for a wide variety of entry environments. Due to material removal resulting from chemical degradation, the structural integrity of TPS material is affected. Spallation occurs when this structurally compromised material is exposed to the high shear flow conditions within the boundary layer. In order to understand the spallation mechanism within WTPS material, we first perform the material removal simulations which occur primarily through oxidation to obtain the microstructure at various stages of degradation. These simulations are performed using the Porous Microstructure Analysis (PuMA) software developed at NASA Ames. The micro-structure geometry used within these simulations were generated artificially to be similar to the 3D weave architecture of MSR-EEV (Mars Sample Return - Earth Entry Vehicle). The various eroded TPS micro-structures are then subjected to the boundary layer flow conditions to obtain critical surface quantities which contribute to the structural failure mechanism such as heat flux, pressure, and shear stress. The direct simulation Monte Carlo (DSMC) methodology is used to perform these simulations in order to accurately capture the strong gradients within the high-temperature boundary layer flow over the intricate geometry of WTPS material. The boundary layer profile is directly taken from the Computational Fluid Dynamics (CFD) simulation and provided as boundary conditions to the DSMC inlet and outlet. Further, the variation of these properties as the microstructure undergoes changes due to oxidation is also investigated. Finally, these quantities are used as input in PuMA to understand the material expansion/compression and strain within the woven TPS geometry and help in developing a comprehensive spallation and structure failure model.
Spallation is the mechanical removal of small chunks of material gets removed typically due to high shear conditions of the flow field. This reduces the ability of the thermal protection system (TPS) material to protect the spacecraft as well as cause turbulence in the flow causing higher heating rates. In this work, we focus on the material removal through ablation and high shear flow within the boundary layer region of woven TPS material. Woven TPS (WTPS) material is the latest class of material developed by NASA, to be used within the next generation of space flights. They are complex interlocked weaves designed to create a rigid structure that is highly resistant to heat and can be easily designed and tailored for a wide variety of entry environments. Due to material removal resulting from chemical degradation, the structural integrity of TPS material is affected. Spallation occurs when this structurally compromised material is exposed to the high shear flow conditions within the boundary layer. In order to understand the spallation mechanism within WTPS material, we first perform the material removal simulations which occur primarily through oxidation to obtain the microstructure at various stages of degradation. These simulations are performed using the Porous Microstructure Analysis (PuMA) software developed at NASA Ames. The micro-structure geometry used within these simulations were generated artificially to be similar to the 3D weave architecture of MSR-EEV (Mars Sample Return - Earth Entry Vehicle). The various eroded TPS micro-structures are then subjected to the boundary layer flow conditions to obtain critical surface quantities which contribute to the structural failure mechanism such as heat flux, pressure, and shear stress. The direct simulation Monte Carlo (DSMC) methodology is used to perform these simulations in order to accurately capture the strong gradients within the high-temperature boundary layer flow over the intricate geometry of WTPS material. The boundary layer profile is directly taken from the Computational Fluid Dynamics (CFD) simulation and provided as boundary conditions to the DSMC inlet and outlet. Further, the variation of these properties as the microstructure undergoes changes due to oxidation is also investigated. Finally, these quantities are used as input in PuMA to understand the material expansion/compression and strain within the woven TPS geometry and help in developing a comprehensive spallation and structure failure model.
The paper presents and compares fluid-flow and heat transfer properties from artificially induced vortices in a flat-plate turbulent boundary layer and naturally occurring vortices due to centrifugal instabilities in a curved-channel laminar flow. Pairs and arrays of vortices are artificially induced by placing half-delta wings on the plate surface. With both arrays and pairs of vortices, streamwise velocities and total pressures are high, and surface heat transfer is locally augmented in vortex downwash regions. In contrast to vortices in the arrays vortices in the pairs tend to move in the streamwise direction with significant divergence (when the common flow between pair is toward the wall) or convergence (when the common flow between pair is away from the wall). The vortices in the arrays cause maximum peak-to-peak heat transfer variations of up to 12 percent of local spanwise-averaged values for initial vortex spacings between 1 to 2.5 generator heights.
The results of an experimental program to evaluate heat transfer and pressure distributions on corrugation roughened flat plates in thick turbulent boundary layers are presented. The experimental program consisted of tests in the tunnel wall boundary layers of the Langley Unitary Plan Wind Tunnel (UPWT) and Continuous Flow Hypersonic Tunnel (CFHT) at freestream Mach numbers of 2.5, 3.5, 4.5, and 10.3. Tests in the UPWT were conducted at a freestream Reynolds number/m of 10,800,000 and in the CFHT, at Reynolds numbers/m of 1,300,000 to 5,800,000. The test configurations consisted of 50.8 x 50.8 cm panels with corrugated beads of two different peak amplitudes, 0.61 and 0.29 cm. The angle of the corrugated beads relative to the flow direction was varied between 0 (aligned) and 90 deg (normal). The measured peak and average heat transfer are analyzed and correlated in terms of the bulk boundary layer, internal boundary layer, and geometric parameters.-
A computational study is presented for the flow inside an oscillatory cavity. The numerical scheme employs a semi-implicit, time-splitting method to integrate the two-dimensional full Navier-Stokes equations satisfying continuity to machine accuracy. The efficient use of direct solvers for the uncoupled momentum and pressure equations is demonstrated. The oscillatory cavity flow is studied considering the effects of heat transfer, Reynolds number, and oscillatory Stokes number.
Structural components in aircraft engines are subjected to multiaxial loads when in service. For such components, life prediction methodologies are dependent on the accuracy of the constitutive models that determine the elastic and inelastic portions of a loading cycle. A threshold surface (such as a yield surface) is customarily used to differentiate between reversible and irreversible flow. For elastoplastic materials, a yield surface can be used to delimit the elastic region in a given stress space. The concept of a yield surface is central to the mathematical formulation of a classical plasticity theory, but at elevated temperatures, material response can be highly time dependent. Thus, viscoplastic theories have been developed to account for this time dependency. Since the key to many of these theories is experimental validation, the objective of this work (refs. 1 and 2) at the NASA Lewis Research Center was to verify that current laboratory techniques and equipment are sufficient to determine flow surfaces at elevated temperatures. By probing many times in the axial-torsional stress space, we could define the yield and flow surfaces. A small offset definition of yield (10 me) was used to delineate the boundary between reversible and irreversible behavior so that the material state remained essentially unchanged and multiple probes could be done on the same specimen. The strain was measured with an off-the-shelf multiaxial extensometer that could measure the axial and torsional strains over a wide range of temperatures. The accuracy and resolution of this extensometer was verified by comparing its data with strain gauge data at room temperature. The extensometer was found to have sufficient resolution for these experiments. In addition, the amount of crosstalk (i.e., the accumulation of apparent strain in one direction when strain in the other direction is applied) was found to be negligible. Tubular specimens were induction heated to determine the flow surfaces at elevated temperatures. The heating system induced a large amount of noise in the data. By reducing thermal fluctuations and using appropriate data averaging schemes, we could render the noise inconsequential. Thus, accurate and reproducible flow surfaces (see the figure) could be obtained.
A set of mathematical formulation is adopted to study vapor deposition from source materials driven by heat transfer process under normal and oblique directions of gravitational acceleration with extremely low pressure environment of 10(exp -2) mm Hg. A series of time animation of the initiation and development of flow and temperature profiles during the course of vapor deposition has been obtained through the numerical computation. Computations show that the process of vapor deposition has been accomplished by the transfer of vapor through a fairly complicated flow pattern of recirculation under normal direction gravitational acceleration. It is obvious that there is no way to produce a homogeneous thin crystalline films with fine grains under such a complicated flow pattern of recirculation with a non-uniform temperature distribution under normal direction gravitational acceleration. There is no vapor deposition due to a stably stratified medium without convection for reverse normal direction gravitational acceleration. Vapor deposition under oblique direction gravitational acceleration introduces a reduced gravitational acceleration in vertical direction which is favorable to produce a homogeneous thin crystalline films. However, oblique direction gravitational acceleration also induces an unfavorable gravitational acceleration along horizontal direction which is responsible to initiate a complicated flow pattern of recirculation. In other words, it is necessary to carry out vapor deposition under a reduced gravity in the future space shuttle experiments with extremely low pressure environment to process vapor deposition with a homogeneous crystalline films with fine grains. Fluid mechanics simulation can be used as a tool to suggest most optimistic way of experiment with best setup to achieve the goal of processing best nonlinear optical materials.
This report presents the modeling progress of a university consortium studying flexible modeling for turbulent heat transfer. In the context of advanced nuclear reactor design, where challenges include non-unity Prandtl fluids, free and mixed convection, and transitional flow, computational fluid dynamics methods are underutilized due to low confidence in modeling approaches and high computational cost. This report evaluates the performance of engineering CFD models in diabatic, buoyant turbulent flow. It finds that all turbulence models including several variants of the k-ε and k-ω models struggle to predict accurate turbulent momentum and heat transfer in such flows. The Nusselt numbers have been compared between the models and the DNS data, where calculations have been performed for each case and trends have shown a good agreement between DNS estimated Nusselt numbers and available correlations as well as experimental data. A novel DNS correlation for high Pr cases as they are transitioning from mixed convection to natural convection has been proposed. While Nusselt number errors relative to DNS range from 20% to 50%, the models capture similar trends to DNS with respect to Nusselt suppression and amplification under varying levels of buoyancy effect. This report also highlights model form error as a significant contributor to CFD predictions and proposes a framework for quantifying model error and improving confidence in CFD calculations. To improve the predictive capability of engineering CFD models, data-driven approaches for turbulence models are investigated. Theoretical frameworks based on the invariant tensor / vector basis neural networks for prediction of Reynolds stresses and turbulent heat fluxes are employed. The models are developed using direct numerical simulations data for forced convection flows of different fluids in vertical planar channel domain. The framework is implemented in spectral element solvers Nek5000 / nekRS and has shown a potential for future development and consideration of mixed convection flows.
The Kelvin–Helmholtz Instability (KHI) is an interface instability that develops between two fluids or plasmas flowing with a common shear layer. KHI occurs in astrophysical jets, solar atmosphere, solar flows, cometary tails, planetary magnetospheres. Two applications of interest, encompassing both space and fusion applications, drive this study: KHI formation at the outer flanks of the Earth’s magnetosphere and KHI growth from non-uniform laser heating in magnetized direct-drive implosion experiments. Here, we study 2D KHI with or without a magnetic field parallel to the flow. We use both the GAMERA code, which solves the compressible Euler equations, and the STRATOSPEC code, which solves the Navier-Stokes equations under the Boussinesq approximation, coupled with the magnetic field dynamics. GAMERA is a global three-dimensional MHD code with high-order reconstruction in arbitrary nonorthogonal curvilinear coordinates, which is developed for a large range of astrophysical applications. STRATOSPEC is a three-dimensional pseudo-spectral code with an accuracy of infinite order (no numerical diffusion). Magnetized KHI is a canonical case for benchmarking hydrocode simulations with extended MHD options. An objective is to assess whether or not, and under which conditions, the incompressibility hypothesis allows to describe a dynamic compressible system. For comparing both codes, we reach the inviscid incompressible regime, by decreasing the Mach number in GAMERA, and viscosity and diffusion in STRATOSPEC. Here, we specifically investigate both single-mode and multi-mode initial perturbations, either with or without magnetic field parallel to the flow. The method relies on comparisons of the density fields, 1D profiles of physical quantities averaged along the flow direction, and scale-by-scale spectral densities. We also address the triggering, formation and damping of filamentary structures under varying Mach number or Atwood number, with or without a parallel magnetic field. Comparisons show very satisfactory results between the two codes. The vortices dynamics is well reproduced, along with the breaking or damping of small-scale structures. We end with the extraction of growth rates of magnetized KHI from the compressible regime to the incompressible limit in the linear regime assessing the effects of compressibility under increasing magnetic field. The observed differences between the two codes are explained either from diffusion or non-Boussinesq effects.
Increased fuel costs, heightened environmental protection requirements, and noise abatement continue to place drag reduction at the forefront of aerospace research priorities. Unfortunately, shortfalls still exist in the fundamental understanding of boundary-layer airflow over aerodynamic surfaces, especially regarding drag arising from skin friction. For example, there is insufficient availability of instrumentation to adequately characterize complex flows with strong pressure gradients, heat transfer, wall mass flux, three-dimensionality, separation, shock waves, and transient phenomena. One example is the acoustic liner efficacy on aircraft engine nacelle walls. Active measurement of shear stress in boundary layer airflow would enable a better understanding of how aircraft structure and flight dynamics affect skin friction. Current shear stress measurement techniques suffer from reliability, complexity, and airflow disruption, thereby compromising resultant shear stress data. The state-of-the-art for shear stress sensing uses indirect or direct measurement techniques. Indirect measurements (e.g., hot-wire, heat flux gages, oil interferometry, laser Doppler anemometry, small scale pressure drag surfaces, i.e., fences) require intricate knowledge of the studied flow, restrictive instrument arrangements, large surface areas, flow disruption, or seeding material; with smaller, higher bandwidth probes under development. Direct measurements involve strain displacement of a sensor element and require no prior knowledge of the flow. Unfortunately, conventional "floating" recessed components for direct measurements are mm to cm in size. Whispering gallery mode devices and Fiber Bragg Gratings are examples of recent additions to this type of sensor with much smaller (m) sensor components. Direct detection techniques are often single point measurements and difficult to calibrate and implement in wind tunnel experiments. In addition, the wiring, packaging, and installation of delicate micro-electromechanical devices impede the use of most direct shear sensors. Similarly, the cavity required for sensing element displacement is sensitive to particulate obstruction. This work was focused on developing a shear stress sensor for use in subsonic wind tunnel test facilities applicable to an array of test configurations. The non-displacement shear sensors described here have minimal packaging requirements resulting in minimal or no disturbance of boundary layer flow. Compared to previous concepts, device installation could be simple with reduced cost and down-time. The novelty lies in the creation of low profile (nanoscale to 100 μm) micropost arrays that stay within the viscous sub-layer of the airflow. Aerodynamic forces, which are related to the surface shear stress, cause post deflection and optical property changes. Ultimately, a reliable, accurate shear stress sensor that does not disrupt the airflow has the potential to provide high value data for flow physics researchers, aerodynamicists, and aircraft manufacturers leading to greater flight efficiency arising from more in-depth knowledge on how aircraft design impacts near surface properties.