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.
For the 2020 Summer Internship period, optimization and analysis was performed of a new iteration of the LIFQWIK PLIF model, used to evaluate the accuracy of CFD simulations of fuel injector performance. This presentation details the accuracy and performance optimization steps taken and presents the results of that optimization.
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.
Digital twin models are virtual representations of physical systems that use real-time data to simulate and optimize performance. This study presents the development and initial implementation of a digital twin (DT) for the electrolyzer-based hydrogen production facility at the National Renewable Energy Laboratory (NREL)'s Advanced Research on Integrated Energy Systems (ARIES), focused on enhancing safety and optimizing sensor placement through physics-based simulations and metadata integration. The DT incorporates detailed facility-specific information, including component layout, leak locations, and controlled release parameters, to model hydrogen dispersion under varying environmental conditions. Using steady-state computational fluid dynamics (CFD) simulations informed by real meteorological data, such as wind speed, direction, and vertical wind profiles, the DT enables visualization of hydrogen plume behavior and spatial concentration distributions. Comparative analysis between high and low wind speed scenarios illustrates the significant influence of wind dynamics on plume shape and extent, with horizontal momentum dominating dispersion at higher speeds, while buoyancy effects become more prominent under low wind conditions. These simulations generate a rich dataset embedded within the DT, allowing users to assess potential leak outcomes and identify optimal sensor locations based on concentration thresholds. The model supports scenario-based analysis to guide safety strategies and equipment deployment for open-area hydrogen infrastructure. The digital twin thus serves as a dynamic platform for virtual prototyping, providing predictive insight into hydrogen behavior and enhancing risk-informed decision-making. This initial phase establishes a validated foundation for future integration of transient, uncontrolled leak scenarios and real-time sensor feedback, positioning the DT as a critical tool for safety design, operational planning, and adaptive monitoring in hydrogen systems. Overall, the approach demonstrates the value of combining environmental data with digital simulations to inform safer and more efficient deployment of hydrogen technologies.
This paper presents an analysis and simulation of evaporation and condensation at a motionless liquid/vapor interface. A 1-D model equation, emphasizing heat and mass transfer at the interface, is solved in two ways, and incorporated into a subgrid interface model within a CFD simulation. Simulation predictions are compared with experimental data from the CPST Engineering Design Unit tank, a cryogenic fluid management test tank in 1-g. The numerical challenge here is the physics of the liquid/vapor interface; pressurizing the ullage heats it by several degrees, and sets up an interfacial temperature gradient that transfers heat to the liquid phase-the rate limiting step of condensation is heat conducted through the liquid and vapor. This physics occurs in thin thermal layers O(1 mm) on either side of the interface which is resolved by the subgrid interface model. An accommodation coefficient of 1.0 is used in the simulations which is consistent with theory and measurements. This model is predictive of evaporation/condensation rates, that is, there is no parameter tuning.
An anti-vortex baffle is a liquid propellant management device placed adjacent to an outlet of the propellant tank. Its purpose is to substantially reduce or eliminate the formation of free surface dip and vortex, as well as prevent vapor ingestion into the outlet, as the liquid drains out through the flight. To design an effective anti-vortex baffle, Computational Fluid Dynamic (CFD) simulations were undertaken for the NASA Ares I vehicle LOX tank subjected to the simulated flight loads with and without the anti-vortex baffle. The Six Degree-Of-Freedom (6- DOF) dynamics experienced by the Crew Launch Vehicle (CLV) during ascent were modeled by modifying the momentum equations in a CFD code to accommodate the extra body forces from the maneuvering in a non-inertial frame. The present analysis found that due to large moments, the CLV maneuvering has significant impact on the vortical flow generation inside the tank. Roll maneuvering and side loading due to pitch and yaw are shown to induce swirling flow. The vortical flow due to roll is symmetrical with respect to the tank centerline, while those induced by pitch and yaw maneuverings showed two vortices side by side. The study found that without the anti-vortex baffle, the swirling flow caused surface dip during the late stage of drainage and hence early vapor ingestion. The flow can also be non-uniform in the drainage pipe as the secondary swirling flow velocity component can be as high as 10% of the draining velocity. An analysis of the vortex dynamics shows that the swirling flow in the drainage pipe during the Upper Stage burn is mainly the result of residual vortices inside the tank due to conservation of angular momentum. The study demonstrated that the swirling flow in the drainage pipe can be effectively suppressed by employing the anti-vortex baffle.
An anti-vortex baffle is a liquid propellant management device placed adjacent to an outlet of the propellant tank. Its purpose is to substantially reduce or eliminate the formation of free surface dip and vortex, as well as prevent vapor ingestion into the outlet, as the liquid drains out through the flight. To design an effective anti-vortex baffle, Computational Fluid Dynamic (CFD) simulations were undertaken for the NASA Ares I vehicle LOX tank subjected to the simulated flight loads with and without the anti-vortex baffle. The Six Degree-Of-Freedom (6-DOF) dynamics experienced by the Crew Launch Vehicle (CLV) during ascent were modeled by modifying the momentum equations in a CFD code to accommodate the extra body forces from the maneuvering in a non-inertial frame. The present analysis found that due to large moments, the CLV maneuvering has a significant impact on the vortical flow generation inside the tank. Roll maneuvering and side loading due to pitch and yaw are shown to induce swirling flow. The vortical flow due to roll is symmetrical with respect to the tank centerline, while those induced by pitch and yaw maneuverings showed two vortices side by side. The study found that without the anti-vortex baffle, the swirling flow caused surface dip during the late stage of drainage and hence early vapor ingestion. The flow can also be non-uniform in the drainage pipe as the secondary swirling flow velocity component can be as high as 10% of the draining velocity. An analysis of the vortex dynamics shows that the swirling flow in the drainage pipe during the Upper Stage burn is mainly the result of residual vortices inside the tank due to the conservation of angular momentum. The study demonstrated that the swirling flow in the drainage pipe can be effectively suppressed by employing the anti-vortex baffle.
An anti-vortex baffle is a liquid propellant management device placed adjacent to an outlet of the propellant tank. Its purpose is to substantially reduce or eliminate the formation of free surface dip and vortex, as well as prevent vapor ingestion into the outlet, as the liquid drains out through the flight. To design an effective anti-vortex baffle, Computational Fluid Dynamic (CFD) simulations were undertaken for the NASA Ares I vehicle LOX tank subjected to the simulated flight loads with and without the anti-vortex baffle. The Six Degree-Of-Freedom (6-DOF) dynamics experienced by the Crew Launch Vehicle (CLV) during ascent were modeled by modifying the momentum equations in a CFD code to accommodate the extra body forces from the maneuvering in a non-inertial frame. The present analysis found that due to large moments, the CLV maneuvering has a significant impact on the vortical flow generation inside the tank. Roll maneuvering and side loading due to pitch and yaw are shown to induce swirling flow. The vortical flow due to roll is symmetrical with respect to the tank centerline, while those induced by pitch and yaw maneuverings showed two vortices side by side. The study found that without the anti-vortex baffle, the swirling flow caused surface dip during the late stage of drainage and hence early vapor ingestion. The flow can also be non-uniform in the drainage pipe as the secondary swirling flow velocity component can be as high as 10% of the draining velocity. An analysis of the vortex dynamics shows that the swirling flow in the drainage pipe during the Upper Stage burn is mainly the result of residual vortices inside the tank due to the conservation of angular momentum. The study demonstrated that the swirling flow in the drainage pipe can be effectively suppressed by employing the anti-vortex baffle.
This chapter describes a new intelligent platform for learning optimal designs of morphing wings based on Variable Camber Continuous Trailing Edge Flaps (VCCTEF) in conjunction with a leading edge flap called the Variable Camber Krueger (VCK). The new platform consists of a Computational Fluid Dynamics (CFD) methodology coupled with a semi-supervised learning methodology. The CFD component of the intelligent platform comprises of a full Navier-Stokes solution capability (NASA OVERFLOW solver with Spalart-Allmaras turbulence model) that computes flow over a tri-element inboard NASA Generic Transport Model (GTM) wing section. Various VCCTEF/VCK settings and configurations were considered to explore optimal design for high-lift flight during take-off and landing. To determine globally optimal design of such a system, an extremely large set of CFD simulations is needed. This is not feasible to achieve in practice. To alleviate this problem, a recourse was taken to a semi-supervised learning (SSL) methodology, which is based on manifold regularization techniques. A reasonable space of CFD solutions was populated and then the SSL methodology was used to fit this manifold in its entirety, including the gaps in the manifold where there were no CFD solutions available. The SSL methodology in conjunction with an elastodynamic solver (FiDDLE) was demonstrated in an earlier study involving structural health monitoring. These CFD-SSL methodologies define the new intelligent platform that forms the basis for our search for optimal design of wings. Although the present platform can be used in various other design and operational problems in engineering, this chapter focuses on the high-lift study of the VCK-VCCTEF system. Top few candidate design configurations were identified by solving the CFD problem in a small subset of the design space. The SSL component was trained on the design space, and was then used in a predictive mode to populate a selected set of test points outside of the given design space. The new design test space thus populated was evaluated by using the CFD component by determining the error between the SSL predictions and the true (CFD) solutions, which was found to be small. This demonstrates the proposed CFD-SSL methodologies for isolating the best design of the VCK-VCCTEF system, and it holds promise for quantitatively identifying best designs of flight systems, in general.
The accurate prediction of pressure drop in the fuel rod bundle is of importance for the design and thermal-hydraulic analysis of nuclear reactor system. Due to very complex geometry of wire-wrapped fuel rod bundle geometry for the Sodium Fast Reactor (SFR), a simple method using hydraulic (equivalent) diameter is not sufficient to predict the pressure drop of rod bundle accurately. For several decades, many experimental works have been performed to develop an accurate correlation for wire-wrapped fuel bundle friction factor. In these works, the friction factor correlation of wire-wrapped fuel bundle has been developed as a function of Reynolds number of the flow and the design parameters such as fuel pin diameter, wire diameter, rod pitch, and wire lead length. The Cheng and Todreas correlations are well-known and used in many subchannel analysis code, such as CADET, TRIO-U, MATRA-LMW, SE2-ANL, etc. Novendstern model, Rehme model, Baxi and Dalle Donne model are also well-known and widely referred. Computational Fluid Dynamics is another typical approach to predict the pressure drop of wire-wrapped fuel rod bundle, but it requires the model validation because the accuracy of pressure drop prediction by the CFD simulation depends on what turbulence model is employed. The present work aims at evaluating the existing pressure drop correlations for the wire-wrapped fuel rod bundle and evaluating the pressure drop of 217-pin rod bundle by the CFD analysis.
It is desirable for fourth-generation Small Modular Reactors to be passively cooled in standard and accident operations. Passive Containment Cooling Systems can reject heat from the containment structure, without using pumps or blowers. The targeted design containment structure is a large, domed, stainless steel, cylindrical vessel. In a postulated Design Basis Accident, steam will flash inside containment. Steam condensation occurs on the inner containment wall and transfers heat through the steel containment into a large body of water known as the annular reservoir (AR) surrounding the vessel serving as the ultimate heat sink. Natural circulation drives the flow in the AR and heat will be released to the environment by evaporation of water. Unique containment geometry requires a separate effects test (SET) facility for the verification and validation of the computer code and evaluation model development and assessment for reactor licensing efforts. In this study, STAR-CCM+, a computational fluid dynamics (CFD) code was used to inform the decision-making process on the design of the SET. The CFD simulation modeled, a two-phase turbulent flow with fluid film development and heat transfer for different containment geometries. The Reactor Excursion and Leak Analysis Program will also be used in a code-to-code verification against the CFD results.
The Waste Treatment and Immobilization Plant (WTP) at the Hanford site is nearing the start of the Direct-Feed Low-Activity Waste (DFLAW) operations. DFLAW is destined to convert a pretreated low activity waste portion of the 56 million gallons of tank waste into a stable solid glass. In the subsequent decade completion of the high-level waste (HLW) facility is anticipated. Sustained operational missions of both LAW and HLW melter facilities are expected over multiple decades. In high-temperature glass melters, the refractory lining corrodes over time, which could potentially be an issue for longer term operations, this refractory corrosion is higher at the level of the glass-air interface due to surface tension driven flow. The glass viscosity, melt pool temperature, and glass chemical composition can impact the rate at which the refractory corrodes. This rate is important to quantify for the various waste glasses to be produced at the WTP since the integrity of the refractory should not be a limiting factor affecting the lifetime of the melter. To this end, a series of glasses representative of the first batches of waste glass produced by the WTP will be melted in small-scale crucibles with Monofrax® K-3 coupons inserted. The corrosion of the K-3 will be measured in the melt and at the meltline (or neckline). A model for the corrosion rate will be constructed and implemented into a previously developed framework for a computational fluid dynamics (CFD) model of the full-scale WTP. To assist with experimental design and validate the implementation of the model in the full-scale melter, CFD simulations of the small-scale crucible tests were performed. The bubbling that occurs in the small-scale crucible is initially validated here with a model that uses silicone oil at room temperature. The viscosity of the oil ranges from 1 to 100 Pa•s, which corresponds to operating glass pool temperatures near 1150 °C down to idling temperatures near 950 °C. The simulation results show good agreement with the bubble sizes that form during experiments. CFD modeling of the crucible setup was used to determine bubbling characteristics to match the range of near-wall velocities expected in the full-scale WTP. This study presents the initial CFD modeling results, corrosion testing plan, and some preliminary corrosion samples with an outline for the next steps for the development of the corrosion model.
High-fidelity Computational Fluid Dynamics (CFD) simulations for multi-rotor vehicles have been carried out. The three-dimensional unsteady Navier-Stokes equations are solved on overset grids employing high order accurate schemes, dual-time stepping, and a hybrid turbulence model using NASA's CFD code Over- flow. The vehicles studied consist of small to medium sized drones, and bigger vehicles for future Urban Air Mobility (UAM) applications. The performances for different configurations and rotor mounting are calculated in hover and in forward flight. Understanding the complex flows and the interactions between rotors and with other elements will help design the future multi-rotor vehicles to be quieter, safer, and more efficient.
Much effort has been made to enhance exploration on Mars. In addition to a rover and Mars-orbiting satellites, a Mars helicopter (MH) was proposed in order to augment planetary research. Computational Fluid Dynamics (CFD) simulations have been performed to have a better understanding of the behavior and performance of vertical lift Planetary Aerial Vehicles (PAV). Due to the large differences in atmospheric conditions between Mars and Earth, predicting and testing rotorcraft performance is a complex task. The goal of this project is to understand the capability of the mid-fidelity CFD software RotCFD to predict rotor performance in terms of thrust at 1013.25 milibar and 14 milibar corresponding to Terrestrial and Martian conditions, respectively. Also, in order to characterize the wind tunnel wall effects free field and wind tunnel simulations were performed, analyzed and compared. Different analytical tools have been used in order to aid with the design process for the future vertical lift planetary aerial vehicles. One of them includes experimental tests performed on a rotor in the Aeolian Wind Tunnel (AWT) facility at NASA Ames Research Center under different pressure conditions ranging from Terrestrial to Martian atmospheric conditions. Other software was used as well in order to capture the aerodynamic coefficients of the corresponding rotor sections based on the Mach and Reynolds numbers used for the experimental tests. The aerodynamic coefficients were input into RotCFD, and various simulations were performed under Terrestrial and Martian conditions in order to mimic the experimental test. Then, the obtained results from RotCFD were compared with the AWT collected data.
A detailed aerothermal model requires the application of the heat flux profile over all exposed aerodynamic surfaces. These loads are created using both engineering methods (e.g., FORTRAN codes implementing correlations) and CFD simulations (e.g., solutions to Navier-Stokes equations). The goal is to define the spatial and temporal dependency of the heat flux upon the model.
A family of high aspect ratio nozzles were designed to provide a parametric database of canonical embedded propulsion concepts. Nozzle throat geometries with aspect ratios of 2:1, 4:1, and 8:1 were chosen, all with convergent nozzle areas. The transition from the typical round duct to the rectangular nozzle was designed very carefully to produce a flow at the nozzle exit that was uniform and free from swirl. Once the basic rectangular nozzles were designed, external features common to embedded propulsion systems were added: extended lower lip (a.k.a. bevel, aft deck), differing sidewalls, and chevrons. For the latter detailed Reynolds-averaged Navier-Stokes (RANS) computational fluid dynamics (CFD) simulations were made to predict the thrust performance and to optimize parameters such as bevel length, and chevron penetration and azimuthal curvature. Seventeen of these nozzles were fabricated at a scale providing a 2.13 inch diameter equivalent area throat." ! The seventeen nozzles were tested for far-field noise and a few data were presented here on the effect of aspect ratio, bevel length, and chevron count and penetration. The sound field of the 2:1 aspect ratio rectangular jet was very nearly axisymmetric, but the 4:1 and 8:1 were not, the noise on their minor axes being louder than the major axes. Adding bevel length increased the noise of these nozzles, especially on their minor axes, both toward the long and short sides of the beveled nozzle. Chevrons were only added to the 2:1 rectangular jet. Adding 4 chevrons per wide side produced some decrease at aft angles, but increased the high frequency noise at right angles to the jet flow. This trend increased with increasing chevron penetration. Doubling the number of chevrons while maintaining their penetration decreased these effects. Empirical models of the parametric effect of these nozzles were constructed and quantify the trends stated above." Because it is the objective of the Supersonics Project that future design work be done more by physics-based computations and less by experiments, several codes under development were evaluated against these test cases. Preliminary results show that the RANS-based code JeNo predicts the spectral directivity of the low aspect ratio jets well, but has no capability to predict the non-axisymmetry. An effort to address this limitations, used in the RANS-based code of Leib and Goldstein, overpredicted the impact of aspect ratio. The broadband shock noise code RISN, also limited to axisymmetric assumptions, did a good job of predicting the spectral directivity of underexpanded 2:1 cold jet case but was not as successful on high aspect ratio jets, particularly when they are hot. All results are preliminary because the underlying CFD has not been validated yet. An effort using a Large Eddy Simulation code by Stanford University predicted noise that agreed with experiments to within a few dB.
The feasibility of using a two-stage unshrouded impeller turbopump to replace the current three-stage reusable launch vehicle engine shrouded impeller hydrogen pump has been evaluated from the standpoint of turbopump weight reduction and overall payload improvement. These advantages are a by-product of the higher tip speeds that an unshrouded impeller can sustain. The issues associated with the effect of unshrouded impeller tip clearance on pump efficiency and head have been evaluated with one-dimensional tools and full three-dimensional rotordynamic fluid reaction forces and coefficients have been established through time dependent computational fluid dynamics (CFD) simulation of the whole 360 degree impeller with different rotor eccentricities and whirling ratios. Unlike the shrouded impeller, the unshrouded impeller forces are evaluated as the sum of the pressure forces on the blade and the pressure forces on the hub using the CFD results. The turbopump axial thrust control has been optimized by adjusting the first stage impeller backend wear ring seal diameter and diverting the second stage backend balance piston flow to the proper location. The structural integrity associated with the high tip speed has been checked by analyzing a 3D-Finite Element Model at maximum design conditions (6% higher than the design speed). This impeller was fabricated and tested in the NASA Marshall Space Flight Center water-test rig. The experimental data will be compared with the analytical predictions and presented in another paper. The experimental data provides validation data for the numerical design and analysis methodology. The validated numerical methodology can be used to help design different unshrouded impeller configurations.
Due to the aft, upper surface engine location on the Hybrid Wing Body (HWB) planform, there is potential to shed vorticity and separated wakes into the engine when the vehicle is operated at off-design conditions and corners of the envelope required for engine and airplane certification. CFD simulations were performed of the full-scale reference propulsion system, operating at a range of inlet flow rates, flight speeds, altitudes, angles of attack, and angles of sideslip to identify the conditions which produce the largest distortion and lowest pressure recovery. Pretest CFD was performed by NASA and Boeing, using multiple CFD codes, with various turbulence models. These data were used to make decisions regarding model integration, characterize inlet flow distortion patterns, and help define the wind tunnel test matrix. CFD was also performed post-test; when compared with test data, it was possible to make comparisons between measured model-scale and predicted full-scale distortion levels. This paper summarizes these CFD analyses.