Search NASA⌕ Search

SEARCH · Search NASA

Results for “Coupled CFD simulation”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Data-driven model for cross ventilation potential in high-density cities based on coupled CFD simulation and machine learning

Effective urban ventilation through decent urban planning and building design can alleviate the deterioration of the urban built environment. Furthermore, natural ventilation requirements and guidelines in current building codes and standards are either qualitative or quantitative but subject to an absolute indoor airspeed threshold without considering the outdoor wind environment. To fill this gap, this paper develops an urban-scale coupled indoor and outdoor computational fluid dynamics (CFD) model and defines a novel ventilation index to assess natural ventilation potential. The index considers wind environments of both indoor and outdoor spaces.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS2003

Calculation of Turbine Axial Thrust by Coupled CFD Simulations of the Main Flow Path and Secondary Cavity Flow in an SLI LOX Turbine

Each single reusable Space Launch Initiative (SLI) booster rocket is an engine operating at a record vacuum thrust level of over 730,000 Ibf using LOX and LH2. This thrust is more than 10% greater than that of the Delta IV rocket, resulting in relatively large LOX and LH2 turbopumps. Since the SLI rocket employs a staged combustion cycle the level of pressure is very high (thousands of psia). This high pressure creates many engineering challenges, including the balancing of axial-forces on the turbopumps. One of the main parameters in the calculation of the axial force is the cavity pressure upstream of the turbine disk. The flow in this cavity is very complex. The lack of understanding of this flow environment hinders the accurate prediction of axial thrust. In order to narrow down the uncertainty band around the actual turbine axial force, a coupled, unsteady computational methodology has been developed to simulate the interaction between the turbine main flow path and the cavity flow. The CORSAIR solver, an unsteady three- dimensional Navier-Stokes code for turbomachinery applications, was used to solve for both the main and the secondary flow fields. Turbine axial thrust values are presented in conjunction with the CFD simulation, together with several considerations regarding the turbine instrumentation for axial thrust estimations during test.

Dorney, D. J.↗
DOE OSTI2020

An error-controlled adaptive time-stepping method for particle advancement in coupled CFD-DEM simulations

Coupled Computational-Fluid-Dynamics (CFD) and Discrete-Element-Method (DEM) models provide an accurate description of multiphase physical systems where a solid granular particle phase exists in an underlying gaseous continuous medium. The time integration of the granular phase in these simulations is typically handled using an explicit scheme with a constant time-step among all particles that is invariant in time to resolve inter-particle collisions. A locally third-order accurate adaptive time integration technique for particles that employs an embedded locally second-order scheme for error determination is presented in this work. The particle time-step size is dynamically adapted based on solution error, thus leading to significant savings in computational time. The efficacy of our scheme is quantified using four test cases of varying complexity (binary collision, homogeneous cooling system, fluidized bed and hopper discharge). The adaptive time-stepping method exhibits improved performance (~ 2–3 times in most of the cases studied) compared to three commonly used non-adaptive time-step methods (first-order Euler-explicit, second-order Adams-Bashforth and third-order Runge-Kutta schemes), while maintaining the same level of accuracy and parallel scalability.

42 ENGINEERING↗
NASA NTRS2015

Development of an Aerodynamic Analysis Method and Database for the SLS Service Module Panel Jettison Event Utilizing Inviscid CFD and MATLAB

This paper describes the development, testing, and utilization of an aerodynamic force and moment database for the Space Launch System (SLS) Service Module (SM) panel jettison event. The database is a combination of inviscid Computational Fluid Dynamic (CFD) data and MATLAB code written to query the data at input values of vehicle/SM panel parameters and return the aerodynamic force and moment coefficients of the panels as they are jettisoned from the vehicle. The database encompasses over 5000 CFD simulations with the panels either in the initial stages of separation where they are hinged to the vehicle, in close proximity to the vehicle, or far enough from the vehicle that body interference effects are neglected. A series of viscous CFD check cases were performed to assess the accuracy of the Euler solutions for this class of problem and good agreement was obtained. The ultimate goal of the panel jettison database was to create a tool that could be coupled with any 6-Degree-Of-Freedom (DOF) dynamics model to rapidly predict SM panel separation from the SLS vehicle in a quasi-unsteady manner. Results are presented for panel jettison simulations that utilize the database at various SLS flight conditions. These results compare favorably to an approach that directly couples a 6-DOF model with the Cart3D Euler flow solver and obtains solutions for the panels at exact locations. This paper demonstrates a method of using inviscid CFD simulations coupled with a 6-DOF model that provides adequate fidelity to capture the physics of this complex multiple moving-body panel separation event.

Applebaum, Michael P.↗
DOE OSTI2022

Sensitivity study of coupled chemical-CFD simulations for analyzing aluminum-clad spent nuclear fuel storage in sealed canisters

We report the United States Department of Energy (DOE) manages over 50 Metric Tons Heavy Metal (MTHM) of aluminum-clad spent nuclear fuel. One main source for DOE’s Aluminum-clad spent nuclear fuel (ASNF) inventory is the advanced test reactor (ATR) at the INL site, which makes this fuel of particular interest for storage scenarios. Road-ready and final disposition packaging configurations for the ATR fuel dictates storage within helium-backfilled, sealed DOE standard canisters. The conditions within these sealed canisters for extended (greater than50 year dry) storage is of interest. To further this goal, a three-dimensional (3D) multi-physics computational fluid dynamics (CFD) model is developed of the sealed DOE standard canisters. This 3D CFD model is one-way coupled with bulk gas radiolysis reactions considering sealed canisters with inert gas and possible trace amounts of air and water vapor. This study looks at the evolution of the thermal history of the canisters over a 50 year time period with a coupling to the chemical reactions occurring from radiolytic breakdown of residual water. A sensitivity study is then carried out over the parameters of the model including the fuel decay heat, residual water content, sealed pressure, canister external temperature, and canister emissivity. In pure helium, hydrogen generation rates are low, under 10 ppm, but hydrogen generation rates are affected greatly by the presence of even 1% residual air, increasing by 50-plus-fold, and nitric acid generation with residual air also occurs ranging from 500 to 4000 ppm after 50 years. The fuel decay heat and the residual water content show the most importance in the generation of hydrogen gas in pure air, and for nitric acid with a residual air condition. External temperature, canister emissivity and sealed pressure all show minor sensitivity effects to the generation of potentially harmful species.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Forecasting Dynamic Line Rating with Spatial Variation Considerations

Dynamic line rating (DLR) is a technology that allows the ampacity of an electrical conductor to be calculated using real-time or forecasted weather conditions. Historically, the ampacity of a conductor has been determined using a static line rating method which assumes conservative weather assumptions. Therefore, not only can DLR give a more accurate measurement of the true ampacity of a conductor, but it can also increase its ampacity during weather conditions with greater thermal mitigations. The two primary cooling factors in the ampacity calculations are wind speed and direction. In complex terrain, wind speed and direction can have large variations over short distances. Therefore, accurately identifying the limiting span of a transmission line requires high spatial resolution of the wind along its path. One solution is to install dense weather stations along their path, though this can become costly over long distances. Therefore, researchers have investigated the use of Computation Fluid Dynamic (CFD) simulations to accurately compute the wind field along the path of a transmission line and use these results to identify the limiting section of the conductor. This work presents a case study that evaluates the coupling of CFD simulations and forecasted weather simulations using the High-Resolution Rapid Refresh (HRRR) model points over a 2-year span within a region in south eastern Idaho. The primary goal of the work is the evaluation of the number of HRRR model points used, i.e., weather stations, along the path of the line and the accuracy of the resulting ampacity. This was done using 4, 10, 17, 26, and 35 HRRR model points along two transmission line paths. The results indicate that as the number of model points are increased, the DLR ampacity of the lines decrease, yet converge as more points are added and demonstrate little change with additional HRRR points. It is expected that these results can help transmission line operators identify the number of weather stations that must be installed when coupled with CFD simulations and DLR ampacity to ensure accurate ratings and safe operations.

17 WIND ENERGY↗
DOE OSTI2021

Forecasting Dynamic Line Rating with Spatial Variation Considerations

Dynamic line rating (DLR) is a technology that allows the ampacity of an electrical conductor to be calculated using real-time or forecasted weather conditions. Historically, the ampacity of a conductor has been determined using a static line rating method which assumes conservative weather assumptions. Therefore, not only can DLR give a more accurate measurement of the true ampacity of a conductor, but it can also increase its ampacity during weather conditions with greater thermal mitigations. The two primary cooling factors in the ampacity calculations are wind speed and direction. In complex terrain, wind speed and direction can have large variations over short distances. Therefore, accurately identifying the limiting span of a transmission line requires high spatial resolution of the wind along its path. One solution is to install dense weather stations along their path, though this can become costly over long distances. Therefore, researchers have investigated the use of Computation Fluid Dynamic (CFD) simulations to accurately compute the wind field along the path of a transmission line and use these results to identify the limiting section of the conductor. This work presents a case study that evaluates the coupling of CFD simulations and forecasted weather simulations using the High-Resolution Rapid Refresh (HRRR) model points over a 2-year span within a region in south eastern Idaho. The primary goal of the work is the evaluation of the number of HRRR model points used, i.e., weather stations, along the path of the line and the accuracy of the resulting DLR ampacity. This was done using 4, 10, 17, 26, and 35 HRRR model points along two transmission line paths. The results indicate that as the number of model points are increased, the DLR ampacity of the lines decrease, yet converge as more points are added and demonstrate little change with additional HRRR points. It is expected that these results can help transmission line operators identify the number of weather stations that must be installed when coupled with CFD simulations and DLR ampacity to ensure accurate ratings and safe operations.

17 WIND ENERGY↗
NASA NTRS2019

Modeling and Simulation Techniques for the NASA SLS Service Module Panel Separation Event; from Loosely-Coupled Euler to Fully-Coupled 6-DOF, Time-Accurate, Navier-Stokes Methodologies

An aerodynamic database has been generated for use by the Orion Multi-Purpose Crew Vehicle (MPCV) Program to analyze Service Module (SM) panel jettison from the NASA SLS vehicle. The database is a combination of CFD data for the panel aerodynamic coefficients, and MATLAB code written to query the CFD data. The Cart3D inviscid CFD flow solver was used to generate the panel aerodynamic coefficients for static panel orientations and free stream conditions that can occur during the jettison event. The MATLAB code performs the multivariate interpolation to obtain aerodynamic coefficients. The MATLAB code uses input for SM panel parameters and returns the SM panel aerodynamic force and moment coefficients for use with a Six-Degree-of-Freedom (6-DOF) motion solver to model the jettison event. This paper examines the accuracy of the sequential-static database approach by modeling the panel jettison event with a fully-coupled, time-dependent, viscous, moving-body CFD simulation. The fully-coupled simulation is obtained using the Loci/Chem unstructured Navier-Stokes CFD solver. The results show that the fully-coupled approach agrees well with the loosely-coupled database/6-DOF approach, indicating that unsteady effects are minimal for the panel jettison event. These results suggest that the database/6-DOF approach is sufficient. In addition, this paper presents the development of an uncertainty model for use in Monte Carlo analysis of the panel jettison event. Here viscous CFD simulations are obtained with Loci/Chem and compared to the inviscid CFD forces and moments. An uncertainty model based on model-form error and numerical error is presented.

Hall, Leslie H.↗
NASA NTRSDate not supplied

Unified Material-Environment Interaction Model for Binary UHTC Composites

A new model for the surface chemistry of zirconium and hafnium diborides containing silicon-carbide when exposed to high temperature conditions in air is formulated. Oxidation dynamics are simplified by considering the diffusion-limited equilibrium regime within the material at elevated temperatures. This model is first assessed against experimental mass-gain data taken from UHTC samples heated to high temperature in an oxygen environment. Then, further evaluations of the new model are conducted using hypersonic CFD simulations to analyze conditions experienced by HfB2-SiC in arc jet experiments conducted at NASA Ames Research Center. Measurements of stagnation point heat transfer and pressure are then used to calibrate the simulations. Results of coupled CFD-material response simulations are then compared to the surface temperatures measured during the Ames arc jet test. Results support the use of the proposed model while highlighting the need for improved gas-phase data on the additional constituents considered.

ultra high temperature ceramic↗
DOE OSTI2026

Plasma-Assisted Pre-Chamber Ignition System for Highly Dilute Stoichiometric Heavy-Duty Natural Gas Engines (Final Technical Report)

This project explored advanced ignition technologies to significantly enhance efficiency and reduce operating costs for heavy-duty natural gas engines operating at stoichiometric conditions, while meeting ultra-low NOx emission standards. The main goal was to develop and validate a plasma-assisted pre-chamber ignition system that could deliver at least a 2% increase in brake thermal efficiency (BTE) and a 4% decrease in total cost of ownership (TCO) compared to a typical multi-cylinder engine with three-way catalyst aftertreatment, ensuring compatibility with the expected 2027 EPA/CARB regulations. In the first half of the project, the research team concentrated on developing and testing plasma-assisted pre-chamber ignition using nanosecond pulsed discharges. Extensive experiments were conducted in an optically accessible rapid-compression and expansion machine, a constant-volume chamber, and an optical single-cylinder engine. Experiments were coupled with CFD simulations. The work produced unique insights into pre-chamber flame formation, jet ignition, dilution effects, and flame quenching at pressures, temperatures, and dilution levels relevant to engines. Although plasma-assisted ignition showed promise in controlled lab settings, the research also identified fundamental and practical challenges when applying this technology to real engine conditions. Midway through the project, a crucial pivot was made, guided by three key findings. First, the power electronics required for nanosecond plasma discharges were found to be too costly for commercial use, undermining the project’s cost-of-ownership goals. Second, nanosecond plasma ignition was highly sensitive to turbulent flow in the pre-chamber, resulting in lower ignition reliability than traditional spark under engine-like conditions. Third, achieving a truly diffuse low-temperature plasma at high pressures near top dead center was not possible, reducing the anticipated chemical enhancement benefits. These results collectively suggested that continuing with plasma-assisted ignition was unlikely to meet both efficiency and cost objectives. In response, the project shifted focus to a more realistic approach: enhancing traditional spark-based pre-chamber ignition with significantly less spark energy. Using insights gained earlier in the project, the team redesigned the pre-chamber to maintain high dilution tolerance and quick combustion, even with lower ignition energy. Testing confirmed that with optimized pre-chamber design and combustion timing, a lower-energy spark could reliably ignite highly diluted stoichiometric mixtures, reduce burn time, and boost thermal efficiency. Final engine testing and techno-economic analysis verified that this revised approach successfully achieved the project goals. The optimized pre-chamber ignition system provided over a 2% increase in calculated brake thermal efficiency compared to the baseline engine. Notably, the lower ignition energy and simplified hardware reduced component stress, extended maintenance intervals, and lowered the total cost of ownership. When used with stoichiometric operation and traditional three-way aftertreatment, the system remained compatible with near-zero NOx emissions targets without increasing cost or complexity in the emissions control system. In summary, although the project deviated from its initial plasma-assisted ignition idea, the work produced a more practical and commercially viable solution. The results show that precisely optimized, low-energy pre-chamber spark ignition can significantly improve efficiency and reduce overall ownership costs for heavy-duty natural gas engines. This directly aligns with DOE goals for cleaner, more efficient, and cost-effective transportation technologies.

03 NATURAL GAS↗
NASA NTRS2010

Ballistics Model for Particles on a Horizontal Plane in a Vacuum Propelled by a Vertically Impinging Gas Jet

A simple trajectory model has been developed and is presented. The particle trajectory path is estimated by computing the vertical position as a function of the horizontal position using a constant horizontal velocity and a vertical acceleration approximated as a power law. The vertical particle position is then found by solving the differential equation of motion using a double integral of vertical acceleration divided by the square of the horizontal velocity, integrated over the horizontal position. The input parameters are: x(sub 0) and y(sub 0), the initial particle starting point; the derivative of the trajectory at x(sub 0) and y(sub 0), s(sub 0) = s(x(sub 0))= dx(y)/dy conditional expectation y = y((sub 0); and b where bx(sub 0)/y(sub 0) is the final trajectory angle before gravity pulls the particle down. The final parameter v(sub 0) is an approximation to a constant horizontal velocity. This model is time independent, providing vertical position x as a function of horizontal distance y: x(y) = (x(sub 0) + s(sub 0) (y-y(sub 0))) + bx(sub 0) -(s(sub 0)y(sub 0) ((y - y(sub 0)/y(sub 0) - ln((y/y(sub 0)))-((g(y-y(sub 0)(exp 2))/ 2((v(sub 0)(exp 2). The first term on the right in the above equation is due to simple ballistics and a spherically expanding gas so that the trajectory is a straight line intersecting (0,0), which is the point at the center of the gas impingement on the surface. The second term on the right is due to vertical acceleration, which may be positive or negative. The last term on the right is the gravity term, which for a particle with velocities less than escape velocity will eventually bring the particle back to the ground. The parameters b, s(sub 0), and in some cases v(sub 0), are taken from an interpolation of similar parameters determined from a CFD simulation matrix, coupled with complete particle trajectory simulations.

Lane, J. E.↗
NASA NTRS2016

On the Assessment of Acoustic Scattering and Shielding by Time Domain Boundary Integral Equation Solutions

Based on the time domain boundary integral equation formulation of the linear convective wave equation, a computational tool dubbed Time Domain Fast Acoustic Scattering Toolkit (TD-FAST) has recently been under development. The time domain approach has a distinct advantage that the solutions at all frequencies are obtained in a single computation. In this paper, the formulation of the integral equation, as well as its stabilization by the Burton-Miller type reformulation, is extended to cases of a constant mean flow in an arbitrary direction. In addition, a "Source Surface" is also introduced in the formulation that can be employed to encapsulate regions of noise sources and to facilitate coupling with CFD simulations. This is particularly useful for applications where the noise sources are not easily described by analytical source terms. Numerical examples are presented to assess the accuracy of the formulation, including a computation of noise shielding by a thin barrier motivated by recent Historical Baseline F31A31 open rotor noise shielding experiments. Furthermore, spatial resolution requirements of the time domain boundary element method are also assessed using point per wavelength metrics. It is found that, using only constant basis functions and high-order quadrature for surface integration, relative errors of less than 2% may be obtained when the surface spatial resolution is 5 points-per-wavelength (PPW) or 25 points-per-wavelength squared (PPW2).

Hu, Fang Q.↗
DOE OSTI2020

Bridging particle and reactor scales in the simulation of biomass fast pyrolysis by coupling particle resolved simulation and coarse grained CFD-DEM

Efficient utilization of biomass is challenging due to the complexity of the feedstock at molecular, particle, and reactor scales. Direct coupling of particle scale model and reactor scale simulation is computationally infeasible. In this work, the results from a particle resolved simulation were used to calibrate the heat transfer and reaction kinetics of particles. A hybrid drag model was used to consider the different fluidization properties of sands and biomass particles. The computation cost was reduced using a coarse-grained Discrete Element Method (DEM). The simulated conversion rate, char yield, tar yield, and pyrolysis gas yield compare well with experiment. The residence time and axial distributions of biomass were also analyzed. This novel multi-scale method provides an efficient and accurate tool for the modeling of biomass pyrolysis reactors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

CFD Simulation of Floating Body Motion with Mooring Dynamics: Coupling MoorDyn with OpenFOAM

It is increasingly popular to use Computational Fluid Dynamics (CFD) models to study floating structures subjected to ocean waves, especially when it comes to applications of floating offshore wind turbines and Wave Energy Converters (WECs). Mooring dynamics are currently lacking in most of these applications. This paper presents a coupled simulation study of moored floating body motion by coupling two open-source libraries: a finite volume CFD toolbox, OpenFOAM, and a lumped-mass mooring model, MoorDyn. The instantaneous floating body position and velocity are passed from the body motion solver in the CFD model to the mooring model to calculate the fairlead kinematics. The mooring reaction forces, which are calculated by MoorDyn after updating the mooring system states, are then returned to the body motion solver to update the floating body motion. Both mesh deformation and overset mesh methods are used as the mesh motion solver in the CFD model to account for the floating body motion. The coupled model was validated against experimental measurements for a floating box moored with four catenary lines under the action of regular waves, which came from a preliminary test campaign for WECs. Apart from the lumped-mass mooring model, the present work also coupled a quasi-static mooring model and a finite element model with the floating body motion solver in OpenFOAM. The mooring line tensions predicted by these models were compared. The coupled model equipped with three mooring line codes may be further used to carry out survivability studies of FOWTs and WECs subject to severe sea states.

floating body↗
NASA NTRS2019

Extraction of Dynamic Inflow Models for Coaxial and Tandem Rotors from CFD Simulations

The dynamic inflow coupling with rotor/body dynamics is crucial in the analysis of stability and control law design for helicopters. Over the past several decades, finite-state inflow models for single rotor configurations in hover, forward flight, and maneuver have developed (Ref.1-3). By capturing the interference effects between rotors, the extension of pressure potential finite state inflow model has promising result for coaxial rotor configuration (Ref.4-6). Recently, the focus of the dynamic inflow modeling has shifted to tandem rotor configurations (Ref.7, 8). The development of the dynamic inflow models for tandem rotor configuration still have some limitations due to the lack of knowledge of rotor-to-rotor interference, and rotor-wake interference. Experimental methods, and computational fluid dynamics methods are commonly used to understand the rotor performance and rotor airload variations, and measure or predict inflow velocity distributions at the rotor desk. The inflow distributions are subsequently used to improve the dynamic inflow models. Tandem rotor configurations have been studied experimentally and computationally for several decades (Ref.9-12). Sweet (Ref.10) observed that a tandem rotor with 76-percent-radius overlap required 14% more induced power at hovering condition, relative to an isolated rotor of equivalent disk area. Sweet also found that, above a shaft-to-shaft distance of 1.03 diameter, the performance of the tandem rotor was nearly the same as two isolated rotors. The objective of the present study is to apply computational fluid dynamics simulations of tandem rotors for the extraction of dynamic inflow models. The extended methodology is first validated by comparing the computed induced power against test data. Subsequently inflow distributions and wake structures are analyzed.

Extraction of Dynamic Inflow Models↗