Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CFD modeling”

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.

At least 397 records · Page 22

Modeling Droplet Heat and Mass Transfer during Spray Bar Pressure Control of the Multipurpose Hydrogen Test Bed (MHTB) Tank in Normal Gravity

A CFD model for simulating pressure control in cryogenic storage tanks through the injection of a subcooled liquid into the ullage is presented and applied to the 1g MHTB spray bar cooling experiments. An Eulerian-Lagrangian approach is utilized to track the spray droplets and capture the interaction between the discrete droplets and the continuous ullage phase. The spray model is coupled with the VOF model by performing particle tracking in the ullage, removing particles from the ullage when they reach the interface, and then adding their contributions to the liquid. A new model for calculating the droplet-ullage heat and mass transfer is developed. In this model, a droplet is allowed to warm up to the saturation temperature corresponding to the ullage vapor pressure, after which it evaporates while remaining at the saturation temperature. The droplet model is validated against the results of the MHTB spray-bar cooling experiments with 50% and 90% tank fill ratios. The predictions of the present T-sat based model are compared with those of a previously developed kinetic-based droplet mass transfer model. The predictions of the two models regarding the evolving tank pressure and temperature distributions, as well as the droplets' trajectories and temperatures, are examined and compared in detail. Finally, the ullage pressure and local vapor and liquid temperature evolutions are validated against the corresponding data provided by the MHTB spray bar mixing experiment.

mass transfer↗

Modeling of Slosh Dynamics in Cryogenic Propellant Tanks in Microgravity Environments

The slosh dynamics in cryogenic fuel tanks under microgravity is a pressing problem that severely affects the reliability of launching spacecraft. After reaching low Earth orbit, the propellant in a multistage rocket experiences large and cyclic changes in temperature as a result of solar heating. Tank wall heating can induce thermal stratification and propellant boiloff, particularly during slosh-inducing vehicle maneuvers. Precise understanding of the dynamic and thermodynamic effects of propellant slosh caused by these maneuvers is critical to mission performance and success. Computational fluid dynamics (CFD) analysis is used extensively within the space vehicle industry in an attempt to characterize the behavior of liquids in microgravity, yet experimental data to quantify these predictions is very limited and reduces confidence in the analytical predictions. A novel approach designed to produce high-fidelity data for correlation to CFD model predictions is being developed with the assistance of Florida Institute of Technology (FIT) and Sierra Lobo, Inc. With few exceptions, previous work in slosh dynamics was theoretical or treated the mass of fuel as a variable of inertia only; such models did not consider the viscosity, surface tension, or other important fluid effects. The challenges in this research are in the development of instrumentation able to measure the required parameters, the computational ability to quantify the fluid behaviors, and the means to assess both the measurements and predictions. The design of this experiment bridges the understanding of slosh dynamics in microgravity by a comprehensive approach that combines CFD tools, dynamic simulation tools, semianalytical models of the predominant fluid effects, and an experimental framework that includes measurement and characterization of liquid slosh in one-degree-of-freedom (DOF) and two-DOF experiments, and ultimately experiments in a NASA low-gravity aircraft.

Source record↗

Preliminary Analysis of a Water Shield for a Surface Power Reactor

A water based shielding system is being investigated for use on initial lunar surface power systems. The use of water may lower overall cost (as compared to development cost for other materials) and simplify operations in the setup and handling. The thermal hydraulic performance of the shield is of significant interest. The mechanism for transferring heat through the shield is natural convection. A simple 1-D thermal model indicates the necessity of natural convection to maintain acceptable temperatures and pressures in the water shield. CFD analysis is done to quantify the natural convection in the shield, and predicts sufficient natural convection to transfer heat through the shield with small temperature gradients. A test program will he designed to experimentally verify the thermal hydraulic performance of the shield, and to anchor the CFD models to experimental results.

Pearson, J. Boise↗

Numerical Investigation of LO2 and LCH4 Storage Tanks on the Lunar Surface

Currently NASA is developing technologies to enable human exploration of the lunar surface for duration of up to 210 days. While trade studies are still underway, a cryogenic ascent stage using liquid oxygen (LO2) and liquid methane (LCH4) is being considered for the Altair lunar lander. For a representative Altair cryogenic ascent stage, we present a detailed storage analysis of the LO2 and LCH4 propellant tanks on the lunar surface for durations of up to 210 days. Both the LO2 and LCH4 propellant tanks are assumed to be pressurized with gaseous helium at launch. A two-phase lumped-vapor computational fluid dynamics model has been developed to account for the presence of a noncondensable gas in the ullage. The CFD model is used to simulate the initial pressure response of the propellant tanks while they are subjected to representative heat leak rates on the lunar surface. Once a near stationary state is achieved within the liquid phase, multizone model is used to extrapolate the solution farther in time. For fixed propellant mass and tank size, the long-term pressure response for different helium mass fractions in both the LO2 and LCH4 tanks is examined.

Moder, Jeff↗

Computational Fluid Dynamics Simulation of the Swordfish Oxygen/Methane Rotating Detonation Rocket Engine

Rotating Detonation Rocket Engines (RDREs) provide several theoretical performance benefits over traditional deflagration combustion rocket engines. RDREs can potentially generate higher specific impulse (ISP) for a given propellant combination enabling significant consumable mass savings for a spaceflight vehicle. Conversely, RDREs are a relatively new technology with a shallower analytical modeling base compared to traditional combustion engine architectures. In this work, we utilize time-accurate, Computational Fluid Dynamics (CFD) modeling to simulate the complex fluid environment inside an RDRE combustion chamber and propellant manifolds. The selected engine of study is the Swordfish RDRE which was developed and tested and at the National Aeronautics and Space Administration (NASA) Marshall Space Flight Center (MSFC). CFD simulations were conducted using the as-tested Swordfish geometry and conditions to facilitate comparison between the analytical predictions and experimental results.

Brian R. Richardson↗

Reduced-Order Aerodynamic Modeling Based on CFD Frequency Responses from Multisine Inputs

A system identification analysis was performed to determine reduced-order models of a computational fluid dynamics (CFD) solver for linear aeroelastic analysis and control design. The application was to the FUN3D code and the flexible half-span wind tunnel test article, in transonic flow conditions, used in the NASA-Boeing collaboration called the Integrated Adaptive Wing Technology Maturation (IAWTM) project. Multiple inputs (structural mode displacements and control surface deflections) were simultaneously excited with orthogonal phase-optimized multisines and multiple outputs (generalized aerodynamic forces) were recorded, from which the matrix of frequency responses were computed using a single CFD run. A state-space model was then fit to the frequency response data using a maximum-likelihood estimator.

System identification↗

Coupling particle scale model and SuperDEM–CFD for multiscale simulation of biomass pyrolysis in a packed bed pyrolyzer

An efficient biomass pyrolysis process requires a comprehensive understanding of the chemical and physical phenomena that occur at multi-length and time scales. In this study, a multiscale computational approach was developed and validated for biomass pyrolysis in a packed-bed reactor by integrating pyrolysis kinetics, a particle scale model, and SuperDEM-CFD in open-source code MFiX. A one-dimensional particle-scale model that discretizes the characteristic length of biomass particle into layers was developed to predict the intraparticle phenomena inside a single particle. The 1D model was validated by comparing it with a single biomass particle pyrolysis experiment. Here, a recently developed SuperDEM-CFD model was employed to simulate the non-spherical particle-particle contact and fluid-particle interaction. The coupled model was applied to simulate the pyrolysis of cubic biomass particles in a packed bed and validated by comparing with experimental data. Simulation with and without particle-scale model was compared, and the effect of the gas-solid heat transfer models was also investigated.

09 BIOMASS FUELS↗

Modeling and analysis of heat emissions from buildings to ambient air

Heat emissions from buildings is a significant source of anthropogenic heat influencing the urban microclimate; however, they are usually oversimplified in urban climate and microclimate modeling. This study developed a bottom-up physics-based approach to calculate heat emissions from buildings to the ambient air and implemented the approach in EnergyPlus. In this work, a simple result verification was conducted by comparing the EnergyPlus simulated results against the spreadsheet calculations. Simulations covering 16 commercial building types, four climates, and two energy efficiency levels were conducted to understand and evaluate the building heat emissions and their temporal patterns as well as three major components: (1) building envelope (convective heat transfer to ambient air), (2) zones (air exfiltration and exhaust air), and (3) HVAC systems (relief air and heat rejection from condensers or cooling towers). The main findings are: (1) heat emissions are usually higher than the site energy use (about 2.5 times), and their dynamics should be considered; (2) building characteristics and their energy systems lead to differences in heat emission contributions from the three components, and their dynamics, for example, in the warehouse models, the envelope component accounts for 90.4%, while it is 12.7% for the large office models; (3) for most building typologies, the climate has a strong impact on heat emissions, for example, buildings with dominant heat emissions from the zone exhaust air and/or the HVAC reject heat, a general decrease in heat emissions in hotter climates is observed, while envelope-dominated buildings show the opposite; and (4) building technologies that reduce energy use in buildings may perform differently in reducing heat emissions. The developed heat emissions calculation method can be adopted in EnergyPlus and most other building energy modeling programs. It can provide dynamic building heat emissions as an input to urban climate computational fluid dynamics (CFD) models at a higher spatial and temporal resolution than is currently available, to improve the simulation accuracy of the urban microclimate and capture the urban heat island effect and urban overheating.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A Computational Fluid Dynamic Study on Polymer Heat Exchangers

Polymer heat exchangers have been developed for the applications involving weight restrictions or chemical compatibility and fouling issues owing to the low density, anticorrosive properties, and low thermal expansion of polymers. Recently development of additive manufacturing also brings new opportunities to make polymer heat exchangers with desired design that was hard to realize before. However, due to the lower thermal conductivity polymer composite, the overall heat transfer performance is still a challenge in the polymer heat exchanger technology. In present work, a computational fluid dynamic (CFD) model has been developed to study the overall heat transfer performance of additively manufactured polymer heat exchangers. The CFD results offer an insight of fluid flow and temperature distribution in the polymer heat exchangers. This study provides a guidance not only on the polymer material selection but also on the design of polymer heat exchangers. The conclusions will be helpful to design a polymer heat exchanger whose overall heat transfer performance is comparable to a metal heat exchanger.

Zhang, Mingkan↗

Modeling and Analysis of Stirling Power Convertors

Modeling and Analysis of Stirling Power Convertors Luis A. Rodriguez1 Steven M. Geng, Terry V. Reid, Scott D. Wilson NASA Glenn Research Center, Cleveland, OH, 44135, USA NASA Glenn’s Thermal Energy Conversion Branch is supporting the development of the next generation free-piston Stirling power convertors. American Superconductor (AMSC) and Sunpower Inc. are the two firms under contract to develop the Flexure Isotope Stirling Convertor (FISC) and the Sunpower Robust Stirling Convertor (SRSC), respectively. To comprehend and forecast convertor performance, Sage, ANSYS® Maxwell, and ANSYS® Fluent were used to model the Stirling thermodynamic cycle, alternator electromagnetics, and piston and displacer dynamics. I. Introduction Stirling convertors are being developed by NASA as a potential steady source of electrical power for NASA’s future scientific space missions. Currently, NASA Glenn Research Center has two corporations under contract, American Superconductor (AMSC) and Sunpower Inc., for the development of the next generation of free-piston Stirling convertors for dynamic radioisotope power systems. AMSC is developing the Flexure Isotope Stirling Convertor (FISC), which uses flexures to prevent side motion and rubbing of the piston. Similarly, Sunpower Inc, is developing the Sunpower Robust Stirling Convertor (SRSC). The SRSC uses gas bearings to prevent radial contact of the moving piston. As convertor development continues, it is increasingly important to understand and predict the interactions of components in the system, how they respond to one another, and how they perform as a response to changes in operating conditions. A suitable and enlightening way to demonstrate and foresee these interactions is with the use of accurate modeling software. Sage, ANSYS® Maxwell, and ANSYS® Fluent are the current modeling tools used by NASA to analytically determine convertor performance. Sage is a one-dimensional object-oriented commercial software package used for modeling and optimizing Stirling convertors for Dynamic Radioisotope Power Systems (DRPS) and it is one of the most accurate Stirling convertor codes in use by NASA. This code is the successor to GLIMPS (Globally-Implicit Stirling Cycle Simulation) and GLOP (GLIMPS Optimization) software created by Gedeon Associates [1]. Model input parameters are typically material/gas type, component physical dimensions, temperatures, frequency, charge pressure, and number of time/space nodes. Sage is used to model both the FISC’s and SRSC’s Stirling cycle thermodynamics and piston/displacer dynamics. Performance maps were created and analyzed for both power systems to better understand the relationship between the following conditions: cold-end temperature, hot-end temperature, piston/displacer amplitudes, pressure drop, and thermal input power. The synergy between these conditions will help determine parameter sensitivity. ANSYS® Maxwell was used to create a three-dimensional (3-D) axisymmetric model for both FISC and SRSC alternators. The significant physical components included in each model are the magnets, magnet carrier, outer/inner laminations, and the coil. Inputs to the model are piston amplitude, piston frequency, alternator load, coil resistance, tuning capacitance, and specific material properties. The alternator models calculate terminal voltage, current, piston/current phase, voltage/current phase, coil inductance, terminal power and efficiency. The RI2 losses, core (hysteresis and eddy) losses, and magnet/can eddy losses are also a part of the final results. ANSYS® Fluent is used to build 3-D computational fluid dynamic (CFD) models to examine the Stirling cycle thermodynamics for both the FISC and SRSC systems. Three-dimensional Computer Aided Design (CAD) models were used to create the physical components of each convertor. Steady-state simulations were conducted for hardware testing, prediction of environmental losses during testing, and generation of radiation look-up tables. The model inputs to the aforementioned analysis are the material properties and boundary thermal conditions. The steady-state model calculates temperature and heat flow distributions. Transient 3-D calculations were also part of the CFD analysis. In this study a physically reduced version of the FISC is used to obtain a prediction of available engine power. For the gas bearing SRSC, the transient effort is used to obtain a prediction of bearing pad performance and its sensitivity to micro-channel geometric variation. The model inputs to the transient simulations are the piston amplitude, displacer amplitude, frequency, displacer/piston phase angle, dynamic deforming CFD grid, temperature boundary conditions, and user defined files describing motion profile of piston/displacer. The results of the model are temperature distributions, heat distributions, and PV power produced at pre-determined conditions.

Luis A Rodriguez↗

A dynamic 2D Borehole Thermal Energy Storage (BTES) model for enhanced computational efficiency

Progressing toward a future increasingly reliant on renewable energy sources, the development of effective, durable energy storage solutions becomes essential to balance supply and demand fluctuations. Borehole Thermal Energy Storage (BTES) is a long-duration thermal energy storage technology that captures excess heat generated from renewable energy sources and stores it underground for later use, enabling the efficient utilization of sustainable energy. This approach is particularly valuable in district energy networks when integrated with Ground Source Heat Pumps (GSHP) to provide stable heating and cooling. However, traditional three-dimensional (3D) numerical models of BTES systems demand extensive computational resources, limiting their practicality for real-time and large-scale applications. This study introduces a novel two-dimensional (2D) modeling approach that reduces computational costs while maintaining high accuracy. By employing a radial ring-based discretization method, the model simulates heat injection, retention, and retrieval dynamics over seasonal cycles. A new thermal-mass weighted-average temperature parameter is introduced to evaluate the performance of BTES systems. Model validation against FEFLOW simulations demonstrates a 17-fold improvement in computational speed compared to traditional Computational Fluid Dynamics (CFD) models while achieving a mean absolute percentage error (MAPE) of 2 % during charging and 4 % during discharging. Additionally, a trade-off analysis between computational efficiency and accuracy is conducted, ensuring the model's applicability for real-world scenarios. The findings of this research contribute to the development of computationally efficient BTES models, facilitating better optimization, control, and integration into renewable energy systems. This work provides a foundation for further studies in techno-economic analysis, multi-year performance evaluation, and real-time operational strategies for BTES applications, supporting a more sustainable energy future.

2D modeling↗

Design of high-deflection foils MHK applications - FEA models - Helical turbines

The Ocean Renewable Power Company's (ORPC's) goal is to design, develop, and test hydrofoils with large deflections. The effects of the deflections on cross-flow turbine performance would be evaluated in order to inform design considerations for full-scale water turbines and other marine hydrokinetic devices. FEA models - NASTRAN Helical foil turbines tested at UNH tow tank Glass and carbon composite material properties Loads derived from CFD models

16 TIDAL AND WAVE POWER↗

Characterizing Performance of Additively Manufacturing Regenerative Cooled Combustion Chambers Through Hot Fire Testing

NASA Marshall Space Flight Center (MSFC) has developed several additive manufactured regenerative cooled combustion chambers from 2014 to present. Chambers have been constructed using the Selective Laser Melting (SLM) powder bed additive manufacturing (AM) technique. The materials used for the chambers has varied from Inconel 718, Inconel 625, GRCop-84, and C18150 metal alloys. The surface finish and build techniques used to successfully manufacture combustion chambers using the SLM AM result in as built surface finish and off nominal geometric features. This paper discusses and reviews the effects of surface roughness, flow area, and hydraulic performance of regenerative cooled combustion chambers. The thrust chambers designs tested ranged from 200 to 1,400 psia producing 1,000 to 35,000 lbf thrust. Empirical data from four test campaigns will be reviewed and related to design modeling parameters such as pressure drop, surface roughness, and actual fluid flow geometry. Correlation results with simple 2-D CFD modeling will also be presented. The goal of the work is to better understand how the finished additive manufactured component characteristics inside internal features can effect regenerative cooling chamber design. Adjustments to design methodology will be discussed and recommended based on the above results.

Garcia, C. P.↗

Multi-Node Modeling of Cryogenic Tank Pressurization System Using Generalized Fluid System Simulation Program

Cryogenic Tanks are pressurized by inert gas such as Helium or Nitrogen to maintain the required pressure of the propellant delivered to the turbo-pump of a liquid rocket engine. Thermo-fluid system simulation tools are used to analyze the pressurization process of a cryogenic tank. Most system level codes (GFSSP and ROCETS) use single node1 to represent ullage which is the gaseous space in the tank. Ullage space in a cryogenic tank is highly stratified because the entering inert gas is at ambient temperature whereas the liquid propellant is at a cryogenic temperature. A single node model does not account for the effect of temperature gradient in the ullage. High fidelity Navier-Stokes based CFD model of Tank Pressurization is not practical for running a long duration transient model with thousands and millions of nodes. A possible recourse is to construct a multi-node model with system level code that can account for ullage stratification.

Bandyopadhyay, Alak↗

Liquid Piston with Spray Cooling Near-Isothermal Compressor

The goal of this project was to prototype and characterize the performance of a liquid-piston spray-cooled gas compressor. The working principle of the compressor enables optimized high-efficiency operation over a very wide range of operating conditions, unlike conventional compressors that are optimized for a narrow range of operating conditions. The compressor technology is suitable for many applications, such as gas pipeline transport, gas storage, and commercial and residential heat pumps. Both physical testing and computational fluid dynamics (CFD) modeling of the processes using the Oak Ridge National Laboratory high-performance computing center were completed. The experimental and CFD studies focused on a near-isothermal liquid-piston compressor (LPC) that uses propylene glycol to compress CO 2 . The first prototype demonstrated isothermal operation during several sequentially executed cycles of CO 2 compression and raised the temperature of the compressed CO 2 by only 2 K, compared with approximately 6 K when the gas was compressed non-isothermally. Isothermal operation was demonstrated at CO 2 flow rates of up to 2 L/min. The second prototype was designed with two compression chambers to allow continuous flow of high-pressure CO 2 . However, the design of the valve train to direct flow between the compression chambers was not sufficient to allow demonstration of CO 2 compression. Numerical simulations of the LPC in which the compression chamber was filled with propylene glycol injected from the bottom inlet were performed using large eddy simulation (LES) with the wall-adapting local eddy-viscosity subgrid-scale model coupled with the multiphase volume of fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfers within the compression chamber. In this effort, the effects of boundary conditions applied to the LES-VOF calculations (i.e., no wall, an adiabatic wall, and a wall with a heat flux subscribed) on the overall pressure and temperature of the CO 2 gas as well as the transient evolution of flow and heat transfer within the compression chamber were investigated and are discussed in this report. The LES calculation with no wall showed no dynamical flow patterns, and the volume-averaged temperature of CO 2 increased from 305 to 392.7 K, whereas LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO 2 temperatures. The results of the LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO 2 gas layers and circulation flow patterns within the chamber that contributed to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor. Though an industrial partnership for commercializing the compressor was not secured, the technology was attractive for an industrial partner to use in two research proposals in response to US Department of Energy funding opportunity announcements.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Numerical Study of Ammonia Leak and Dispersion in the International Space Station

Release of ammonia into the International Space Station (ISS) cabin atmosphere can occur if the water/ammonia barrier breach of the active thermal control system (ATCS) interface heat exchanger (IFHX) happens. After IFHX breach liquid ammonia is introduced into the water-filled internal thermal control system (ITCS) and then to the cabin environment through a ruptured gas trap. Once the liquid water/ammonia mixture exits ITCS, it instantly vaporizes and mixes with the U.S. Laboratory cabin air that results in rapid deterioration of the cabin conditions. The goal of the study is to assess ammonia propagation in the Station after IFHX breach to plan the operation procedure. A Computational Fluid Dynamics (CFD) model for accurate prediction of airflow and ammonia transport within each of the modules in the ISS cabin was developed. CFD data on ammonia content in the cabin aisle way of the ISS and, in particular, in the Russian On- Orbit Segment during the period of 15 minutes after gas trap rupture are presented for four scenarios of rupture response. Localized effects of ammonia dispersion and risk mitigation are discussed.

Son, Chang H.↗

Using RotCFD to Predict Isolated XV-15 Rotor Performance

Experimental techniques to measure rotorcraft aerodynamic performance are widely used. However, the need exists to understand the limitations of ground based testing by augmenting the analysis of experimental test results with Computational Fluid Dynamics (CFD) modeling. The immediate objective of the present research is to develop an XV-15 Tilt Rotor Research Aircraft rotor model for investigation of wind tunnel wall interference. The predicted performance of the XV-15 during various flight modes is compared to theoretical and experimental data. This research is performed to support wind tunnel tests scheduled for 2016. A mid-fidelity RANS solver, RotCFD, is used with an unsteady, incompressible flow model and a realizable k- turbulence model. The rotor is modeled using an actuator disk model or blade element model with a momentum source approach. In RotCFD the setup, grid generation and running of cases is faster than many CFD codes which makes it a useful engineering tool. Performance predictions need not be as accurate as high-fidelity CFD codes, as long as wall effects can be properly simulated. Being able to accurately predict unsteady rotorcraft performance on desktop-class computers provides a quicker analysis of highly complex flows during the initial design phase.

CFD↗

Microstructural and infiltration properties of woven preforms during chemical vapor infiltration

Interface-resolved direct numerical simulations (DNSs) of chemical vapor infiltration (CVI) have been performed over a range of furnace-operating conditions (Thiele moduli) and for practical woven preform geometries. A level-set method is used to resolve the geometry of the initial preform at tow scale. The interface between the vapor and solid phase is then evolved in time through the entire CVI densification cycle, fully resolving the time-varying topology between the two phases. In contrast to previous level-set methods for CVI simulation, the physical reaction and diffusion processes govern the level-set movement in the current approach. The surface deposition kinetics is described by the usual one-step model. In this paper, the DNS data are used to study the evolving porosity, surface-to-volume ratio, and flow infiltration properties (permeability and effective diffusivities). Additionally, comparisons are made to popularly-assumed structure functions and the standard, Kozeny–Carmen porous media model commonly employed in modeled CFD simulations of CVI. The virtual DNS experiments reveal a Thiele modulus and preform geometry (fabric layup) dependence which the existing microstructural and infiltration models are not able to describe throughout the entire densification process. The DNS-based, woven geometry-specific correlations can be applied directly to mean-field, furnace-scale CFD simulations.

36 MATERIALS SCIENCE↗