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.
Research conducted at the Institute for Computer Applications in Science and Engineering in applied mathematics, numerical analysis, fluid mechanics including fluid dynamics, acoustics, and combustion, aerodynamics, and computer science during the period 1 Apr. 1992 - 30 Sep. 1992 is summarized.
High-temperature gas-cooled reactors (HTGRs) are rapidly approaching deployment. Confidence in transient analysis of these systems for design, optimization, and licensing calculations requires modeling and simulation tools that have been validated against data relevant to HTGR conditions. The High Temperature Test Facility (HTTF) is an integral effects thermal hydraulics test facility for prismatic HTGRs. In spring and summer of 2019, HTTF was used for a series of experiments that now serve as the basis for the OECD/NEA Thermal Hydraulic Code Validation Benchmark for High Temperature Gas-Cooled Reactors using HTTF Data (HTGR T/H Benchmark). This benchmark contains problems for systems code, computational fluid dynamics (CFD), and coupled systems code/CFD modeling representing lower plenum mixing and both the depressurized and pressurized conduction cooldown (DCC and PCC respectively) transients. Benchmark problems include exercises for code-to-code and code-to-data comparisons as well as an exercise for error scaling between HTTF and the Modular High Temperature Gas-Cooled Reactor, which serves as the basis for the HTTF design. Previous analysis as part of the HTGR T/H benchmark used a RELAP5-3D model developed at Idaho National Laboratory (INL) and demonstrated an ability to reproduce trends in the measured data but difficulties reproducing experimental values within their uncertainty. These difficulties were largely attributed to assumptions made during the development of the initial RELAP5-3D model, which predated the HTTF experiments. A significant cause of difficulty reproducing the measured temperatures may be the radial nodalization of the previous RELAP5-3D model. The new model provides a finer nodalization to assess the impact of radial nodalization and allows for asymmetric heating within the core, which was a feature of multiple HTTF experiments. In this paper, we present the new RELAP5-3D model of HTTF. In addition to describing the new model, this paper compares the new and old models and provides results for a full-power steady state, a DCC, and a PCC in HTTF. These analyses are based on the code-to-code comparison exercises for the DCC and PCC problems of the HTGR T/H benchmark. We present the results of these exercises from the new model and compare them to the results of the old model.
Solar thermal power tower systems are the primary technology being proposed for solar electricity from thermal energy. Operational limits on these towers are often driven by mechanical properties under significant thermal loads, particularly at the receiver where incoming flux is converted to thermal energy. While the solar receiver’s efficiency is largely driven by its optical properties, thermomechanical stresses on the receiver limits the operational envelope. One pathway to higher efficiency is greater allowable solar fluxes on the receiver but novel materials are required. The present study uses computational fluid dynamics to describe a parametric design space for a microvascular carbon/carbon composite solar receiver as a new material option for high flux solar receivers. Simulations are conducted for different microvascular geometries considering the role of material properties and heat transfer fluids, for the impact on thermal efficiency, and allowable strain. Results show that microscale receiver modules made of the proposed carbon/carbon composite could achieve thermal efficiencies over 90% and full-scale receivers can achieve up to 85% thermal efficiency for the design explored considering realistic strain limits, flux levels, and material properties. These values are highly dependent on the heat transfer fluid pairing, the through plane thermal conductivity of the carbon/carbon composite, the path architecture of the microscale receiver, and the incident solar flux profiles.
In laser powder bed fusion (L-PBF) additive manufacturing, the mechanical performance, microstructure and defects of fabricated parts are closely associated with the melt pool morphology, e.g., its dimension and shape through the building process. Past studies have largely focused on how the process parameters such as laser power and scan speed affect melt pool characteristics. In this study, the melt pool morphology variation as a function of preheating temperature in the conduction, transition, and keyhole regimes and the underlying mechanisms in each regime are investigated through ex-situ sample characterization and computation thermal fluid dynamics (CtFD) simulation. Single tracks with different combinations of laser power and scan speed are deposited on an Inconel 718 bare plate preheated to a temperature range of 100–500 °C in the experiment. Significant changes are observed in melt pool morphology as a function of preheating temperature from optical measurements of melt track cross sections. The depth of melt pool in the three regimes increases monotonically with preheating temperature, e.g., at 500 °C, the experimental melt pool depth is increased by 49% in conduction regime, 34% in transition regime and 33% in keyhole regime, respectively, while the variation of melt pool width in each regime does not all follow an increasing trend but depends on the melt pool regimes. Melt pool width variation in the conduction and transition regimes is found to depend on the enhanced heat conduction directly related to temperature dependent thermal properties. Through validated CtFD simulations, it is found that in the keyhole regime the evaporation mass, recoil pressure, and laser drilling effect is enhanced with higher preheating temperature, which gives rise to a deeper melt pool. The simulations also reveal that preheating temperature significantly elongates the melt track length due to the increased flow rate and strong recoil pressure that accelerates the backward flow.
An enhanced performance solid rocket booster concept for the space shuttle system has been proposed. The concept booster will have strong commonality with the existing, proven, reliable four-segment Space Shuttle Reusable Solid Rocket Motors (RSRM) with individual component design (nozzle, insulator, etc.) optimized for a five-segment configuration. Increased performance is desirable to further enhance safety/reliability and/or increase payload capability. Performance increase will be achieved by adding a fifth propellant segment to the current four-segment booster and opening the throat to accommodate the increased mass flow while maintaining current pressure levels. One development concept under consideration is the static test of a "standard" RSRM with a fifth propellant segment inserted and appropriate minimum motor modifications. Feasibility studies are being conducted to assess the potential for any significant departure in component performance/loading from the well-characterized RSRM. An area of concern is the aft motor (submerged nozzle inlet, aft dome, etc.) where the altered internal flow resulting from the performance enhancing features (25% increase in mass flow rate, higher Mach numbers, modified subsonic nozzle contour) may result in increased component erosion and char. To assess this issue and to define the minimum design changes required to successfully static test a fifth segment RSRM engineering test motor, internal flow studies have been initiated. Internal aero-thermal environments were quantified in terms of conventional convective heating and discrete phase alumina particle impact/concentration and accretion calculations via Computational Fluid Dynamics (CFD) simulation. Two sets of comparative CFD simulations of the RSRM and the five-segment (IBM) concept motor were conducted with CFD commercial code FLUENT. The first simulation involved a two-dimensional axi-symmetric model of the full motor, initial grain RSRM. The second set of analyses included three-dimensional models of the RSRM and FSM aft motors with four-degree vectored nozzles.
Numerical simulations of a freely bubbling cylindrical fluidized bed are carried out using a coupled computational fluid dynamics and discrete element method (CFD-DEM) model and compared to recent experimental data. The experiments were conducted using high-resolution and high-frequency magnetic resonance imaging providing high-fidelity data of the bubbling within a central 10 mm slice of the bed. Qualitatively, we nd more regular (i.e., less chaotic) structures in the simulated beds than observed experimentally. Quantitatively, however, the bubble diameter and number of bubbles as a function of height within the bed is predicted well by the base model. Unfortunately, the regularity in the simulations manifests as a considerable discrepancy in the speed of the (dense) emulsion phase. The simulated velocity probability distribution functions show an accumulation of low-speed regions and deficiency of high-speed regions. A simple parametric study of the base model is also carried out considering many of the most common CFD-DEM modeling parameters. It is found that the fluid grid size, geometry resolution, transfer kernel and drag law did not have a significant effect on bubble or particle dynamics.
Due to the challenging design requirements, elbows are often unavoidable in duct configuration, and these 90-degree bends introduce swirl, velocity variations, recirculation, and secondary flow. These disturbances make it difficult for nuclear facilities to meet particle sampling standards. A series of numerical analyses are conducted to track aerosols in a nuclear facility duct having a 90-degree elbow with the assistance of computational fluid dynamics (CFD). A turbulence model, a continuity, and a momentum, a discrete phase model, and species transport equations are solved simultaneously to track aerosols in the duct. The effect of turbulence models, turbulent dispersion models, droplet drag model, aerosol amount, aerosol spray configuration, guide vanes, and mixers are investigated. Simulation results are analyzed per relevant testing codes such as DOE-HDBK-1169, ASME AG1, ISO 14644-3, ACGIH, and ANSI/HPS N13.1.
Critical technologies are being developed to support the goals of the NASA Office of Aeronautics and Space Transportation Technology Access to Space initiative for next-generation reusable space transportation systems. From the perspective of aerothermodynamic performance throughout the flight trajectory, the Reusable Launch Vehicle program incorporates conceptual analysis, ground-based testing, and computational fluid dynamics to provide flyable suborbital flight demonstrator vehicles. This report provides an overview of the hypersonic aeroheating wind tunnel test program conducted at the NASA Langley Research Center in support of one of these vehicles, the X-34 small reusable technology demonstrator program. Global surface heat transfer images, surface streamline patterns, and shock shapes were measured on 0.0153- and 0.0183-scale models of proposed X-34 flight vehicles at Mach 6 and 10 in air. The primary parametrics that were investigated include angles-of-attack from 0 to 35 deg. and freestream unit Reynolds numbers from 0.5 to 8 million per foot (which was sufficient to produce laminar, transitional, and turbulent heating data), both with and without control surface deflections. Comparisons of the experimental data to computational predictions are included, along with a discussion of the implications of some of the experimental flow features for the flight vehicle.
A computational model based on the axisymmetric, thin-layer Navier-Stokes equations is developed to predict the convective, radiation and conductive heat transfer in high temperature space nuclear reactors. An implicit-explicit, finite volume, MacCormack method in conjunction with the Gauss-Seidel line iteration procedure is utilized to solve the thermal and fluid governing equations. Simulation of coolant and propellant flows in these reactors involves the subsonic and supersonic flows of hydrogen, helium and uranium tetrafluoride under variable boundary conditions. An enthalpy-rebalancing scheme is developed and implemented to enhance and accelerate the rate of convergence when a wall heat flux boundary condition is used. The model also incorporated the Baldwin and Lomax two-layer algebraic turbulence scheme for the calculation of the turbulent kinetic energy and eddy diffusivity of energy. The Rosseland diffusion approximation is used to simulate the radiative energy transfer in the optically thick environment of gas core reactors. The computational model is benchmarked with experimental data on flow separation angle and drag force acting on a suspended sphere in a cylindrical tube. The heat transfer is validated by comparing the computed results with the standard heat transfer correlations predictions. The model is used to simulate flow and heat transfer under a variety of design conditions. The effect of internal heat generation on the heat transfer in the gas core reactors is examined for a variety of power densities, 100 W/cc, 500 W/cc and 1000 W/cc. The maximum temperature, corresponding with the heat generation rates, are 2150 K, 2750 K and 3550 K, respectively. This analysis shows that the maximum temperature is strongly dependent on the value of heat generation rate. It also indicates that a heat generation rate higher than 1000 W/cc is necessary to maintain the gas temperature at about 3500 K, which is typical design temperature required to achieve high efficiency in the gas core reactors. The model is also used to predict the convective and radiation heat fluxes for the gas core reactors. The maximum value of heat flux occurs at the exit of the reactor core. Radiation heat flux increases with higher wall temperature. This behavior is due to the fact that the radiative heat flux is strongly dependent on wall temperature. This study also found that at temperature close to 3500 K the radiative heat flux is comparable with the convective heat flux in a uranium fluoride failed gas core reactor.
This keynote presentation explores the behavior of headers—essential components of pipeline systems—using ANSYS simulation software and machine learning techniques. The study aims to predict the thermal and mechanical performance of headers under diverse conditions through both steady-state and transient simulations. We investigate critical parameters such as heat transfer coefficient, fluid velocity, and temperature to optimize header design. Conducted as part of a DOE project led by NCAT in collaboration with UNC Charlotte, this research encompasses multiple key topics. The initial section focuses on the behavior of header systems under steady-state conditions using ANSYS simulation. It underscores the importance of headers in industrial infrastructure, especially in the energy sector, and examines the implications of material selection and flow direction on heat transfer dynamics. Methodologically, we employ Computational Fluid Dynamics (CFD) analysis through ANSYS, detailing the development of models, material properties, geometry specifications, boundary conditions, and meshing strategies. Our simulations explore various operational parameters, including temperature and mass flow rates, crucial for predicting heat transfer coefficients and enhancing header design. Results from the study include parametric investigations into mesh sensitivity, viscosity model evaluations, and the effects of heat transfer locations, all validated against theoretical calculations. We conclude with insights on mesh optimization, the suitability of viscosity models, and recommendations for future research aimed at improving header system efficiency and sustainability in industrial applications.
A computational study focused on leading-edge radius effects and associated Reynolds number sensitivity for a High Speed Civil Transport configuration at transonic conditions was conducted as part of NASA's High Speed Research Program. The primary purposes were to assess the capabilities of computational fluid dynamics to predict Reynolds number effects for a range of leading-edge radius distributions on a second-generation supersonic transport configuration, and to evaluate the potential performance benefits of each at the transonic cruise condition. Five leading-edge radius distributions are described, and the potential performance benefit including the Reynolds number sensitivity for each is presented. Computational results for two leading-edge radius distributions are compared with experimental results acquired in the National Transonic Facility over a broad Reynolds number range.
Ensuring the efficient removal of decay heat from the reactor vessel is essential for the safety of advanced reactor technologies. Several Generation-IV concepts incorporate variations in the reactor vessel cooling systems to achieve this objective. High-temperature gas-cooled reactors utilize a reactor cavity cooling system (RCCS), a passive ex-vessel system designed to operate without active components or external power during accident conditions. The RCCS removes decay heat primarily through radiative and convective heat transfer mechanisms. Here, this study presents a comprehensive validation of a computational fluid dynamics Reynolds-averaged Navier-Stokes model for the University of Wisconsin-Madison air-cooled RCCS facility. Validation was conducted for both high- and low-power natural convection cases under a uniform heating profile. Near-wall resolution was found to be critical for accurately modeling natural convection in the RCCS; employing an all-𝑦 + wall treatment resulted in wall temperature discrepancies exceeding 50 °𝐶 compared to a wall-resolved mesh. Thermal-hydraulic behaviors under natural and forced convection conditions were compared within the heated cavity and RCCS. A turbulence model sensitivity analysis indicated that low-Reynolds number k-ɛ, k-ω shear stress transport (SST), and Reynolds stress transport models produce similar wall temperature predictions. A buoyancy modeling sensitivity study revealed that the Boussinesq approximation significantly underpredicted thermal-hydraulic behavior in the RCCS. Based on these findings, modeling recommendations are provided. The validated data set along with identified sensitivities refine the modeling of natural convection in the RCCS. The information produced by this study supports RCCS design, optimization, and safety evaluations, enabling the calibration and verification of reduced-order thermal-hydraulic models.
On-orbit fluid transfer such as refueling of propellant tanks and life-support systems can enable long-duration space missions. For safe and efficient liquid transfer operations, prior knowledge of liquid positioning and liquid-vapor interface behavior while draining in a low-gravity environment is required. Numerical models capable of predicting vapor ingestion (or vapor pull-through) can be used to design liquid transfer operations while reducing liquid residuals, mission risk and settling thrust required to prevent vapor ingestion. An experimental program conducted in the 2.2 Second Drop Tower facility at NASA Glenn Research Center investigated the vapor ingestion phenomenon for a range of outflow rates and tank sizes providing a database for validation. This study presents a Computational Fluid Dynamics model capable of accurately predicting the vapor ingestion using the Volume-of-Fluid multiphase solver in commercial code STAR-CCM+. A description of the experimental setup and general trends from similar studies are presented. Comparisons of the numerical prediction and test data in normal and low gravity show good agreement and give confidence in pursuing design of full-scale propellant transfer systems.
The lack of a reliable method to evaluate the convergence of molecular dynamics simulations has contributed to discrepancies in different areas of molecular dynamics. Here, the method of information entropy is introduced to molecular dynamics for stationarity assessment. The Shannon information entropy formalism is used to monitor the convergence of the atom motion to a steady state in a continuous spatial domain and is also used to assess the stationarity of calculated multidimensional fields such as the temperature field in a discrete spatial domain. It is demonstrated in this work that monitoring the information entropy of the atom position matrix provides a clear indicator of reaching steady state in radiation damage simulations, non-equilibrium molecular dynamics thermal conductivity computations, and simulations of Poiseuille and Couette flow in nanochannels. A main advantage of the present technique is that it is non-local and relies on fundamental quantities available in all molecular dynamics simulations. Unlike monitoring average temperature, the technique is applicable to simulations that conserve total energy such as reverse non-equilibrium molecular dynamics thermal conductivity computations and to simulations where energy dissipates through a boundary as in radiation damage simulations. The method is applied to simulations of iron using the Tersoff/ZBL splined potential, silicon using the Stillinger–Weber potential, and to Lennard–Jones fluid. Its applicability to both solids and fluids shows that the technique has potential for generalization to other areas in molecular dynamics.
On-orbit fluid transfer such as refueling of propellant tanks and life-support systems can enable long-duration space missions. For safe and efficient liquid transfer operations, prior knowledge of liquid positioning and liquid-vapor interface behavior while draining in a reduced-gravity environment is required. Numerical models capable of predicting vapor ingestion (or vapor pull-through) can be used to design liquid transfer operations while reducing liquid residuals, mission risk and settling thrust required to prevent vapor ingestion. An experimental program conducted in the 2.2 Second Drop Tower facility at NASA Lewis Research Center in 1969 investigated the vapor ingestion phenomenon for a range of outflow rates and tank sizes providing a database for validation. This study presents a Computational Fluid Dynamics model capable of accurately predicting the vapor ingestion using the Volume-of-Fluid multiphase solver in commercial code STAR-CCM+. A description of the experimental setup and general trends from similar studies are presented. Comparisons of the numerical prediction and test data in normal and low gravity show good agreement and give confidence in pursuing design of full-scale propellant transfer systems.
On-orbit fluid transfer such as refueling of propellant tanks and life-support systems can enable long-duration space missions. For safe and efficient liquid transfer operations, prior knowledge of liquid positioning and liquid-vapor interface behavior while draining in a reduced-gravity environment is required. Numerical models capable of predicting vapor ingestion (or vapor pull-through) can be used to design liquid transfer operations while reducing liquid residuals, mission risk and settling thrust required to prevent vapor ingestion. An experimental program conducted in the 2.2 Second Drop Tower facility at NASA Lewis Research Center in 1969 investigated the vapor ingestion phenomenon for a range of outflow rates and tank sizes providing a database for validation. This study presents a Computational Fluid Dynamics model capable of accurately predicting the vapor ingestion using the Volume-of-Fluid multiphase solver in commercial code STAR-CCM+. A description of the experimental setup and general trends from similar studies are presented. Comparisons of the numerical prediction and test data in normal and low gravity show good agreement and give confidence in pursuing design of full-scale propellant transfer systems.
The use of a flow augmented thermal management system for entry and re-entr environments is one method for reducing heat and drag loads. This concept relies on jet penetration from supersonic and hypersonic counterflowing jets that could significantly weaken and disperse the shock-wave system of the spacecraft flow field. The objective of this research effort is to conduct parametric studies of the supersonic flow over a 2.6% scale model of the Apollo capsule, with and without the counterflowing jet, using time-accurate and steady-state computational fluid dynamics simulations. The numerical studies, including different freestream Mach number angle of attack counterflowing jet mass flow rate, and nozzle configurations, were performed to examine their effect on the drag and beat loads and to explore the counternowing jet condition. The numerical results were compared with the test data obtained from transonic blow-down wind-tunnel experiments conducted independently at NASA MSFC.
This article examines the use of computational fluid dynamics (CFD) simulations to predict the response of a floating wind platform to irregular-wave excitation. This work was conducted as part of the Reproducible CFD Joint Industry Project for Floating Offshore Wind Applications and involved verification across several participants and modeling tools and validation against experiments. The authors pay special attention to the uncertainties in both CFD results and experimental measurements. In this work, we perform detailed comparisons of the incident waves and the motion of a moored structure. The nonlinear, low-frequency resonance motion is of particular interest because it potentially drives the mooring and tower-base loads. The verification and validation study is partially successful in that the CFD simulations capture the low-frequency slow-drift motion well but underpredict the low-frequency pitch resonance. This underprediction can be attributed in part to the CFD incident waves, which showed some discrepancies with the experimental waves, especially around extreme events. The effects of the wave discrepancies are also estimated and investigated using a mid-fidelity OpenFAST model. Overall, the present study increases our confidence in using CFD simulations to predict the global performance of offshore wind platforms in irregular waves and to produce data for the calibration of lower-fidelity models.