Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CFD simulations”

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 721 records · Page 40

Photophoretic Propulsion Enabling Mesosphere Exploration NIAC Phase I Final Report

This Phase I report presents a comprehensive study on photophoretic flyers—innovative, ultralight, solar-powered vehicles that harness photophoretic forces generated via Knudsen pumping to achieve sustained flight in the mesosphere (50–80 km altitude). By integrating advanced materials such as nanocardboard— characterized by its extremely low areal density (~1 g/m²) and high bending stiffness—with ultrathin light-absorbing coatings, the project converts incident solar radiation directly into a directed thrust. Extensive experimental investigations, coupled with high-fidelity computational fluid dynamics (CFD) simulations using ANSYS Fluent, validate the concept across various three-dimensional geometries, including spherical, conical, and rocket-shaped configurations. These simulations bridge the gap between free-molecular and continuum flow regimes, demonstrating that optimized designs can generate lift forces sufficient to support kilogram-scale payloads even in low-pressure environments. At the heart of this innovation is the use of Knudsen pumping, where temperature gradients across porous surfaces induce directional gas flow, creating a modest overpressure that provides lift. The report introduces an analytical framework that interpolates between the well-known low-Reynolds number drag regime and the high-Reynolds number momentum theory. This model accurately predicts lift based on design parameters such as microchannel dimensions, porous wall geometry, areal density, and nozzle exit area. For instance, simulations indicate that 10-meter-scale structures with carefully engineered porous walls can achieve the necessary pressure differential to support scientifically significant payloads (~1 kg). The study also explores a hybrid propulsion approach that combines solar buoyancy with photophoretic lift. Initially, solar heating creates a buoyant force that elevates the flyer to mesospheric altitudes. Once in the optimal pressure range, the photophoretic mechanism—powered by Knudsen pumping—takes over as the primary source of lift, ensuring stable, long-duration flight. This dual-mode operation not only facilitates the deployment of photophoretic flyers but also broadens the potential applications for mesospheric exploration. In addition to propulsion, the report investigates the integration of photophoretic thrusters for trajectory control of existing research balloons in the upper stratosphere. By dynamically adjusting the nozzle orientation and controlling the flow-through velocity, these thrusters provide precise maneuverability, enabling the flyers to counteract atmospheric disturbances and adjust their flight paths in real time. For example, a photophoretic thruster approximately 7.5 by 7.5 meters in size could be unfolded below a payload gondola of a 60 million-cubic-foot zero-pressure balloon. Such a thruster can provide horizontal speed control of approximately 1 m/s using only sunlight and no moving parts (except those needed to track the Sun and control the jet direction). Importantly, photophoretic thrusters operate more efficiently at higher altitudes, which is complementary to known trajectory control techniques, such as propellers and tethered wings, which are more effective at lower altitudes. Finally, the report identifies three scientific research thrusts where mesospheric aircraft technology can have a profound impact: atmospheric tides, characterization of gravity waves, and investigation of mesospheric instabilities. Overall, the findings of this Phase I project represent a significant advancement in photophoretic propulsion technology. By demonstrating that large-scale, ultralight structures can be powered solely by solar radiation—via carefully engineered Knudsen pumping—this work lays a robust foundation for scalable, near-space flight architectures. Future refinements in material fabrication, structural optimization, and integrated trajectory control are expected to further enhance performance, paving the way for operational demonstrations that could revolutionize atmospheric science, remote sensing, and communication networks.

Knudsen Pump↗

Boundary Layer Analysis of Shock Tube Flows

Shock tubes offer a controlled environment to reproduce kinetic and radiative phenomena characteristic of atmospheric entry flows under ground-test conditions. The boundary layer developing behind the incident shock wave determines the available test time, influences particle residence times important for similarity scaling, and can affect radiative energy transport. In this study, we couple a quasi-1D space marcher with the compressible boundary-layer equations to numerically compute the post-shock flow in a shock frame of reference for various test gas mixtures representative of different planetary atmospheres. The solvers are individually verified against CFD simulations and analytical correlations available in the literature. Coupled solutions are computed for a finite-rate chemistry in the boundary layer and a non-catalytic isothermal wall. Results yield refined estimates of the maximum separation distance, as well as insights into concentration profiles of relevant species within the boundary layer.

Andrea Fagnani↗

A Transient Hydrodynamic Model of Screen Channel Liquid Acquisition Devices for In-Space Cryogenic Propellant Management

Screen channel liquid acquisition devices (LADs) will play a crucial role in future deep space travel. It is essential that vapor-free delivery of propellants during tank-to-tank transfer is ensured to maximize yield from storage tanks and prevent potential combustion instabilities. The screen channel LAD utilizes a fine screen wire mesh that can separate phases in a low Bond number (i.e. microgravity) environment using surface tension forces. This study presents the development and verification of a new model for transient screen compliance, one of the influential factors for screen channel LAD design. Screen compliance is crucial during LAD channel outflow transients because the slight deflection of the screen can provide needed mass to satisfy rapid outflow demands and reduce the pressure difference across the screen. The model is successfully verified against CFD simulations. In addition, the characteristic speed for the governing screen compliance equations is derived which allows for numerical stability criteria to be established. As shown in this study, the transient maximum pressure difference across the screen can greatly exceed the steady state maximum pressure difference across the screen in many cases.

Hydrodynamics Simulations↗

Thermal Stress Modeling and Analysis of Packed-bed Thermocline Energy Storage Tank for INL Thermal Energy Distribution System (TEDS)

The Thermal Energy Distribution System (TEDS) at Idaho National Laboratory (INL) is a thermal-hydraulic flow loop to support the integration of co-located multiple experimental systems, where a packed-bed thermal energy storage (TES) is installed as a thermal buffer and storage unit for TEDS. The packed-bed TES is adopted in TEDS because of its benefit as a low-cost single-tank storage option compared to the traditional two-tank storage. However, thermal ratcheting is one potential design concern which is caused by the rearrangement of granular filler inside a packed-bed tank during continuous thermal cycling operation of the packed-bed TES tank. If the thermally induced stress exceeds yield strength of the tank wall, it may cause catastrophic consequences like rupture of the thermal storage tank. Thus, it is crucial to understand the phenomenon to ensure the robust operation. Based on the temperature boundary conditions given by transient thermal analyses with computational fluid dynamics (CFD) simulations, the thermal ratcheting analysis is then conducted to evaluate the hoop stress and resultant thermal ratcheting potential of the TES tanks with two different modeling approaches: (1) infinite rigidity model and (2) Drucker-Prager (DP) model. The validity of each modeling method was examined by comparing the numerical simulation with the experimental data obtained from the packed-bed TES tank for Solar One Plant and evaluate the thermal ratcheting potential of the TEDS TES tank.

25 ENERGY STORAGE↗

Technology Enablers for Advanced Gasoline Compression Ignition Engines

To meet current emission regulations and increasingly demanding global fleet CO2 standards on fuel economy and future trends towards life cycle GHG emissions, advanced combustion engines remain significant in the passenger vehicle sector to achieve high efficiency and low emissions over the full operating range. Lean burn gasoline compression ignition (GCI) technology has shown to have the most potential in reaching these goals, although it faces challenges in the operating range. In this study, an advanced GCI engine is considered with the capability to operate under two combustion modes, namely low temperature combustion (LTC) and GCI. They are enabled with the use of two Hyundai in-house developed technologies; an advanced valve control mechanism known as continuously variable valve duration (CVVD) and a high-pressure gasoline injection system. At low load, the engine utilizes dual CVVD and dual CVVT (continuously variable valve timing) mechanisms for both intake and exhaust valvetrains to enable NVO (negative valve overlap) to trap hot residuals. The hot residuals enable low load auto-ignition, and the phasing control of that auto-ignition is achieved by varying the degrees of NVO, and in turn, the amount of hot residuals. This is done in real time with respect to the changes in engine loads or operating conditions, significantly improving cylinder-phasing control in LTC mode, and subsequently is able to realize a fuel economy gain and reduced emission benefits. Early direct injection during NVO or intake stroke is used to form homogenous or quasi-homogenous air-fuel mixture for LTC. At mid to high loads, late direct injection of gasoline or GCI mode is employed to take advantage of the high volatility and high resistance to auto-ignition of gasoline fuel, so that diesel-engine like high combustion efficiency can be achieved with lower soot emissions for the same NOx level as diesel combustion. In this paper, the key control and fueling technologies that enable this multiple mode combustion are introduced, and the resultant analysis on engine test and CFD simulation at LTC and GCI operating conditions are presented in detail.

Gasoline compression ignition, low temperature com↗

Study of effect of PWR cold leg temperature gradient on reactor core condition

Effects of temperature and flow gradients in Westinghouse designed three-loop Pressurized Water Reactor (PWR) cold legs, the piping between the main coolant pump and the reactor vessel, were evaluated using Computational Fluid Dynamics (CFD) code STAR CCM+ and coupled neutronic and thermal-hydraulic (T/H) code system VERA. In the parametric study, several symmetric and asymmetric temperature gradients that were significantly larger than those observed from plant measurements were applied to the cold leg inlets for comparison with the base case without any temperature gradient. A CFD model using the STAR-CCM+ code was developed for a portion of the RCS region between the Reactor Coolant Pump (RCP) and the core inlet based on previously validated modeling approach. The CFD simulation results were processed for the temperature and flow rate distributions at the core inlet as input to the VERA calculations. The VERA code system consists of COBRA-TF (CTF) for thermal-hydraulics, MPACT for reactor physics and neutron transport, and ORIGEN for isotopic depletion. The VERA model was for depletion calculations of a high-burnup loading pattern with the reactor core in pin-by-pin and subchannel resolution. The results of the study indicate that the postulated temperature gradients within the PWR cold legs do not result in any significant changes in the core inlet temperature distributions and the core power distributions during the reactor operation. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Impacts of Biomass Feedstock Pre-Processing on Heat and Mass Transfer During Pyrolysis Using X-Ray Computed Tomography and Multiscale Modeling

Knowledge of the transport properties of biomass particles such as porosity, tortuosity, and permeability is paramount for high-fidelity modeling of biomass pyrolysis due to the heat and mass transfer limitations imposed by particle microstructure. X-ray computed tomography (XCT) is a non-destructive imaging method that enables full 3D reconstructions of the biomass particle microstructure with high resolution, permitting direct calculation of porosity, tortuosity, and permeability from real particle geometries. In this study, XCT imaging revealed the 3D microstructures of particles and chars from pyrolytic conversion of cylindrically cut or milled/pelletized loblolly pine samples. The porosity, tortuosity, and permeability were calculated directly from the XCT geometries via open-source microstructural analysis tool MATBOX+TauFactor (https://github.com/NREL/MATBOX_Microstructure_analysis_toolbox) and computational fluid dynamics (CFD) simulations using our solver, Mesoflow (https://github.com/NREL/mesoflow). These properties were used in a reactor scale model developed in COMSOL of the single particle reactor at NREL to investigate the impact of feedstock pre-processing on biomass conversion during pyrolysis with rigorous experimental validation.

biomass↗

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↗

Use, Assessment, and Improvement of the Loci-CHEM CFD Code for Simulation of Combustion in a Single Element GO2/GH2 Injector and Chamber

This document is a viewgraph presentation of a paper that documents a continuing effort at Marshall Space Flight Center (MSFC) to use, assess, and continually improve CFD codes to the point of material utility in the design of rocket engine combustion devices. This paper describes how the code is presently being used to simulate combustion in a single element combustion chamber with shear coaxial injectors using gaseous oxygen and gaseous hydrogen propellants. The ultimate purpose of the efforts documented is to assess and further improve the Loci-CHEM code and the implementation of it. Single element shear coaxial injectors were tested as part of the Staged Combustion Injector Technology (SCIT) program, where detailed chamber wall heat fluxes were measured. Data was taken over a range of chamber pressures for propellants injected at both ambient and elevated temperatures. Several test cases are simulated as part of the effort to demonstrate use of the Loci-CHEM CFD code and to enable us to make improvements in the code as needed. The simulations presented also include a grid independence study on hybrid grids. Several two-equation eddy viscosity low Reynolds number turbulence models are also evaluated as part of the study. All calculations are presented with a comparison to the experimental data. Weaknesses of the code relative to test data are discussed and continuing efforts to improve the code are presented.

Westra, Douglas G.↗

CFD Technology for Rotorcraft Gearbox Windage Aerodynamics Simulation

A computational fluid dynamics (CFD) method is adapted, validated and applied to spinning gear systems with emphasis on predicting windage losses. Several spur gears and a disc are studied. The CFD simulations return good agreement with measured windage power loss. Turbulence modeling choices, the relative importance of viscous and pressure torques with gear speed and the physics of the complex 3-D unsteady flow field in the vicinity of the gear teeth are studied.

Handschuh, Robert↗

Controls/CFD Interdisciplinary Research Software Generates Low-Order Linear Models for Control Design From Steady-State CFD Results

The NASA Lewis Research Center is developing analytical methods and software tools to create a bridge between the controls and computational fluid dynamics (CFD) disciplines. Traditionally, control design engineers have used coarse nonlinear simulations to generate information for the design of new propulsion system controls. However, such traditional methods are not adequate for modeling the propulsion systems of complex, high-speed vehicles like the High Speed Civil Transport. To properly model the relevant flow physics of high-speed propulsion systems, one must use simulations based on CFD methods. Such CFD simulations have become useful tools for engineers that are designing propulsion system components. The analysis techniques and software being developed as part of this effort are an attempt to evolve CFD into a useful tool for control design as well. One major aspect of this research is the generation of linear models from steady-state CFD results. CFD simulations, often used during the design of high-speed inlets, yield high resolution operating point data. Under a NASA grant, the University of Akron has developed analytical techniques and software tools that use these data to generate linear models for control design. The resulting linear models have the same number of states as the original CFD simulation, so they are still very large and computationally cumbersome. Model reduction techniques have been successfully applied to reduce these large linear models by several orders of magnitude without significantly changing the dynamic response. The result is an accurate, easy to use, low-order linear model that takes less time to generate than those generated by traditional means. The development of methods for generating low-order linear models from steady-state CFD is most complete at the one-dimensional level, where software is available to generate models with different kinds of input and output variables. One-dimensional methods have been extended somewhat so that linear models can also be generated from two- and three-dimensional steady-state results. Standard techniques are adequate for reducing the order of one-dimensional CFD-based linear models. However, reduction of linear models based on two- and three-dimensional CFD results is complicated by very sparse, ill-conditioned matrices. Some novel approaches are being investigated to solve this problem.

Melcher, Kevin J.↗

CFD modeling and simulation for corrosive wear of refractory in molten slag

The development of high-wear resistant refractories having minimal production costs is facilitated by characterizing the wear mechanisms associated with their corrosive wear. Static cup testing is a commonly used method for comparing the corrosion resistance performance of two or more refractory materials. Although the static cup test conditions are not as severe as dynamic tests, this study shows that the thermal gradient present within the system during heating and cooling stages serves to generate movement of the slag leading to mechanical wear. The thermal gradient within the refractory, and between the slag and the refractory, occurs during the ramp stage of the test and lasts until the soaking stage is reached bringing the system to a thermal equilibrium. Using computational fluid dynamics (CFD) capabilities embedded within ANSYS software, this study modelled and quantified the convection currents within the slag and associated shear stresses generated on the refractory walls due to the thermal gradient. A traditional ladle furnace was employed as a case study to verify the results of the studied CFD model. The corrosion rate of the refractory lining was found to depend on the mass transfer coefficient of the refractory dissolution into the slag, and a velocity term which governs the extent of corrosion at any given location. This velocity term is a function of slag viscosity, as well as the concentration gradient and/or temperature gradient at the triple points. In this study, wall shear stress was used as a reliable proxy for identifying high-velocity regions prone to excessive corrosive wear. Elevated wall shear stress near the slag/air and slag/molten steel interfaces align with observed corrosion grooves, which reflects the intensified corrosive wear at these locations.

Ramteke, Rajat Rajat Durgesh Ramteke [University o↗

Studying Rate of Fire Spread Using Landscape Statistics and FIRETEC

Wildfires are multi-system and multi-scale processes that are influenced by the land and atmosphere. Here, we study the influence of land, in terms of topography and vegetation, on uphill and downhill fire spread using data of the state of California as a testbed. The California landscape is first quantified by reclassifying its dominant vegetation types into grass, trees, and shrubs. Next, statistics of slope, vegetation height, and live fuel moisture percentage for each dominant vegetation type are gathered to create ensembles of idealized computational fluid dynamics (CFD) model simulations representative of the landscape. Rate of spread values for each CFD simulation are recorded to observe the differences in perturbing slope and fuel moisture for each vegetation type. We have found that all vegetation types prefer slopes of less than 30 degrees. Trees in the California data are 10 to 40 meters high, whereas shrub heights are less than 3 meters and grass heights are less than 1 meter. The live fuel moisture percentage of both trees and shrubs fluctuates around approximately 100% over one year (2019-2020). Using these parameters as inputs into the CFD model, fire spread rate is found to be fastest in the grass ensemble (maximum 12 m s -1 ), followed by shrubs (maximum 10 m s -1 ) and lastly by trees (maximum 7 m s -1 ). Increasing slope increases the spread rate of the uphill fire but the downhill spread decreases for slopes greater than 15 degrees. Increasing fuel moisture decreases the spread rate both up and down the hill. This study highlights the impacts that topography and vegetation have on fire spread.

54 ENVIRONMENTAL SCIENCES↗

User Interface Developed for Controls/CFD Interdisciplinary Research

The NASA Lewis Research Center, in conjunction with the University of Akron, is developing analytical methods and software tools to create a cross-discipline "bridge" between controls and computational fluid dynamics (CFD) technologies. Traditionally, the controls analyst has used simulations based on large lumping techniques to generate low-order linear models convenient for designing propulsion system controls. For complex, high-speed vehicles such as the High Speed Civil Transport (HSCT), simulations based on CFD methods are required to capture the relevant flow physics. The use of CFD should also help reduce the development time and costs associated with experimentally tuning the control system. The initial application for this research is the High Speed Civil Transport inlet control problem. A major aspect of this research is the development of a controls/CFD interface for non-CFD experts, to facilitate the interactive operation of CFD simulations and the extraction of reduced-order, time-accurate models from CFD results. A distributed computing approach for implementing the interface is being explored. Software being developed as part of the Integrated CFD and Experiments (ICE) project provides the basis for the operating environment, including run-time displays and information (data base) management. Message-passing software is used to communicate between the ICE system and the CFD simulation, which can reside on distributed, parallel computing systems. Initially, the one-dimensional Large-Perturbation Inlet (LAPIN) code is being used to simulate a High Speed Civil Transport type inlet. LAPIN can model real supersonic inlet features, including bleeds, bypasses, and variable geometry, such as translating or variable-ramp-angle centerbodies. Work is in progress to use parallel versions of the multidimensional NPARC code.

Source record↗

Thermal-Hydraulics Modeling and Simulations of Hot Pool Using the SAS-CFD Coupled Code

The main goal of this activity is to test the dynamic coupling of the SAS4A/SASSYS-1 (SAS) and CFD models, using a recently patched version of the SAS code intended to address an undocumented limitation that hindered the Versatile Test Reactor (VTR) simulation efforts in FY21. As described in previous VTR calculation reports, the undocumented limitation in SAS v5.4 does not allow the user to activate the CFD coupling option during restart calculations. Since the analysts were unaware of this limitation, prior SAS-CFD simulation results for the protected station blackout (PSBO) transient were erroneous. Root-cause analysis was performed to determine the cause of this undocumented limitation in SAS v5.4, the SAS software was updated in a new patch, and the SAS-CFD simulations were repeated with this patched software. The results of the SAS-CFD simulations documented in this report show that the software patch does address the cited issue, and that the patched software indeed supports the activation of the CFD coupling model in restart calculations. The SAS development team will determine the schedule for implementing the patch in an official software release. This report documents updated SAS-CFD simulations of the PSBO transient response in the VTR. The hot pool is modeled with the CFD code STAR-CCM+, which is coupled at the flow boundaries to the SAS model of the primary heat transport system. SAS computes the mass flow rate and temperature at each core subassembly outlet, the thermal insulation cavity bypass, and the IHX inlet windows. CFD in turn computes the absolute pressure and temperature at each of these boundaries. The SAS code will ignore the temperature data at flow boundaries where flow is directed into the hot pool, i.e., at the core subassembly outlets unless flow reversal occurs. Similarly, CFD will ignore temperature data at boundaries where the flow is directed out of the hot pool, i.e., at the IHX inlets except under flow reversal. The focus of this work is to ensure that the SAS software patch addresses the undocumented limitation described in prior VTR calculation reports, rather than the accurate assessment of thermal stratification in the VTR during protected transients. This motivates the development of a new, simplified CFD model with a coarser mesh to accelerate the testing process. The updated model, and simplifying assumptions, are documented in this report. In future work, the thermal stratification assessment should be performed in more detail. The simplified CFD model can be improved by performing grid convergence studies sensitivity studies of turbulence parameters (e.g., Prandtl number, turbulence production and dissipation parameters) on temperature distributions and thermal stratification.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

CFD Approaches for Simulation of Wing-Body Stage Separation

A collection of computational fluid dynamics tools and techniques are being developed and tested for application to stage separation and abort simulation for next-generation launch vehicles. In this work, an overset grid Navier-Stokes flow solver has been enhanced and demonstrated on a matrix of proximity cases and on a dynamic separation simulation of a belly-to-belly wing-body configuration. Steady cases show excellent agreement between Navier-Stokes results, Cartesian grid Euler solutions, and wind tunnel data at Mach 3. Good agreement has been obtained between Navier-Stokes, Euler, and wind tunnel results at Mach 6. An analysis of a dynamic separation at Mach 3 demonstrates that unsteady aerodynamic effects are not important for this scenario. Results provide an illustration of the relative applicability of Euler and Navier-Stokes methods to these types of problems.

Buning, Pieter G.↗

Application of CFD codes for the simulation of scramjet combustor flowfields

An overview of CFD activities in the Hypersonic Propulsion Branch is given. Elliptic and PNS codes that are being used for the simulation of hydrogen-air combusting flowfields for scramjet applications are discussed. Results of the computer codes are shown in comparison with those of the experiments where applicable. Two classes of experiments will be presented: parallel injection of hydrogen into vitiated supersonic air flow; and normal injection of hydrogen into supersonic crossflow of vitiated air.

Chitsomboom, Tawit↗