SEARCH · Search NASA

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,621 records · Page 90

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Development and evaluation of a skeletal mechanism for EHN additized gasoline mixtures in large Eddy simulations of HCCI combustion

Advanced Low Temperature Combustion modes, such as the Sandia proposed Additive-Mixing Fuel Injection (AMFI), can unlock significant potential to boost fuel conversion efficiency and ultimately improve the energy conversion of internal combustion engines. This is a novel improved combustion process that is enabled by supplying small (<5%) variable amounts of autoignition improver to the fuel to enhance the engine operation and control. Common, diesel-fuel ignition-quality enhancing additive, 2-ethylexyl nitrate (EHN), is doped into gasoline to enable Sandia LTGC + AMFI combustion. This manuscript focuses on the development of a reduced sub-mechanism for EHN chemical kinetics at engine relevant conditions that is implemented into a skeletal mechanism for chemical kinetic studies of gasoline surrogate fuels. The mechanism validation utilized zero-dimensional numerical simulations and comparison to shock tube ignition-delay data of pure and EHN-doped n-heptane. Additional validation is presented with Homogeneous Charge Compression-Ignition (HCCI) engine data of pure and EHN-doped research-grade E10 gasoline. Then, the mechanism was deployed in a 3-D computational fluid dynamics (CFD) using Large Eddy Simulations (LES) to model the HCCI engine experiments of 0.4% vol EHN additized E10 gasoline at several equivalence ratios. Simulations showed a very good performance of the mechanism, and the model accurately reproduced (a) the ignition point, (b) combustion phasing, (c) combustion duration, and (d) the peak of the heat release rates of the engine experiments. The results show that EHN promotes Low-Temperature Heat Release, ultimately driving the gasoline to autoignite at thermodynamic conditions where the fuel would not otherwise ignite. Overall, this work demonstrates a viable reduced chemical-kinetic mechanism for EHN and shows that it can be combined with a skeletal gasoline mechanism for CFD-LES analysis of well-mixed LTGC that matches well with experimental results. The CFD-LES analysis also shows the spatial distribution of EHN-fuel interactions that control the autoignition throughout the combustion chamber.

Guleria, Gaurav↗

DOE OSTI2024

MFiX Development Updates

This presentation discusses recent developments of the Multiphase Flow with Interphase eXchanges (MFiX) software. A brief review of all modeling approaches along with their cost versus accuracy is provided to guide users when selecting a model for a given application. The main new features of the past five official releases of MFiX are described. Improvement in chemistry management allow for easier simulation setup, and faster simulation speed. Progress in the implementation of a thin-wall boundary condition is presented. The major new model development over the past year is the release of the Glued Sphere Particle model (GSP), where component spheres are combined together to represent non-spherical particles. The integration of Machine Learning (ML) workflow in the CFD process is discussed with two applications: a surrogate model for stiff chemistry and the development of a PIC stress model using ML.

Dietiker, Jeff↗

NASA NTRS2019

A Streamlined Approach to Spatial Mapping of Complex 3D Thermal Boundary Condition Data

Spatial mapping of 1D, 2D, and particularly 3D thermal boundary condition data is essential for accurate thermal simulation of spacecraft and launch vehicles. This information typically is available in simple ASCII file format but can be problematic to use for thermal simulations. Difficulties that can be encountered include the ability to directly use the available file format, transformation of the data into the desired reference coordinate system in 3D space, and accurate interpolation of the data onto the thermal calculation domain. This paper presents a streamlined approach to spatial mapping of complex unsymmetrical 3D cold wall heat flux data using the fields’ capability in Simcenter3D. A recent example of 3D thermal plume simulation at MSFC will be employed to illustrate the use of text, csv, and Excel files, specifying coordinate systems for transformation of data in 3D space, and use of the various interpolation schemes available for spatial mapping. Spatial boundary condition verification is also very important, and the spacecraft thermal analyst needs effective visualization tools to develop confidence in the boundary condition definitions. Visualization of the spacecraft plume raw field data, the calculated field data, and the interpolated field data onto the thermal calculation domain will be discussed, with a particular focus on real-time viewing of the interpolated data for available interpolation schemes. The objective of this paper is to outline the modeling techniques used with the integrated thermal analysis of a terrestrial lunar lander demonstrator. Specifically, the techniques outlined herein focus on spatial mapping of convective and radiative exhaust plume results from a computational fluid dynamics (CFD) analysis. Two cases were investigated, however, the modeling techniques that were utilized translate to both cases seamlessly. Those two cases were as follows: • Lunar lander demonstrator firing on-pad (worst case). • Lunar lander demonstrator hovering at an elevated position above the launch pad. For the purposes of information control, this paper will focus on the modeling techniques only, using altered or dimensionless data where necessary. BACKGROUND The thermal model examined in this paper is of the lunar lander demonstrator known as XL-1T (terrestrial), born from a collaborative effort between Masten Space Systems (MSS) and NASA, under NASA’s Lunar CATALYST (also known as Lunar Cargo Transportation and Landing by Soft Touchdown) initiative. The lander is a reusable vertical takeoff/vertical landing (VTVL) test bed which is controlled by four throttleable main engines utilizing green hypergolic propellants.

Hawkins, Robert F.↗

DOE OSTI2019

Utilization of Hydroturbines in Wastewater Treatment Plants

Wastewater treatment plants (WWTPs) are a significant energy consumer, yet there are several opportunities for implementing on-site power generation systems. Within the treatment process, the high flow rate of effluent is produced and discharged to a nearby water body by gravity. Thus, hydroturbines can be utilized to generate power in such an application due to a difference in elevation and high flow rate. This paper presents a case study of introducing a hydroturbine in a WWTP in Wisconsin and evaluating the power output in addition to determining the energy savings. The WWTP considered in this study has an effluent flow rate of 190 MGD (million gallons per day) and elevation difference of 3 m (10 ft) between the final stage of treatment and the discharge point. Based on the parameters above; hubless rim-drive Kaplan type hydroturbine (RDT) is the optimal choice to be used in such an application. The RDT is designed and optimized by using in-house code. A computational fluid dynamics (CFD) software is applied to evaluate the performance of the proposed model, and the system is simulated through homer software to validate the results generated by the CFD. The expected savings is estimated to be 1564 MWh/yr.

Energy & Fuels↗

DOE OSTI2025

Summary of Multiphysics Modeling for Sublimation Mass Transfer

This report summarizes the effort to develop a multiphysics modeling framework for sublimation mass transfer processes. The aim of the project was to develop a multiphysics code capable of predicting phase change between solid and gas in a closed container, and the movement of material within such a container when exposed to various exterior environmental conditions. Two modeling frameworks were developed towards this aim: one with a high-fidelity computational fluid dynamics (CFD) structure, and the other as a fast reduced-order model. Ultimately, both attempts were unsuccessful due to instabilities in the code and physical processes for which there is no adequate numerical representation. Both modeling attempts are briefly detailed before providing a brief survey of recent updates in the literature, and a recommendation for future work on this subject.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2015

Numerical Calibration of Mass Flow Plug for Inlet Testing

A simple control volume model has been developed to calculate the discharge coefficient through a mass flow plug (MFP) and validated with a calibration experiment. The maximum error of the model within the operating region of the MFP is 0.54%. The control volume analysis developed work is comprised of a sequence of flow calculations through the MFP. The model uses the MFP geometry and operating pressure and temperature to couple continuity, momentum, energy, an equation of state, and wall shear. The discharge coefficient calculation also includes the effects of boundary layer growth, including the reduction in cross-sectional flow area as characterized by the boundary layer displacement thickness. The last calculation in the sequence uses an integral method to calculate the growth of the boundary layer, from which the displacement thickness is then determined. The result of these successive calculations is an accurate one-dimension model of the velocity, pressure, and temperature through the MFP. For comparison, a computational fluid dynamic (CFD) calibration is shown, which when compared to the presented numerical model, had a lower accuracy with a maximum error of 1.35% in addition to being slower by a factor of 100."

Flow Measurement↗

NASA NTRS2008

Use of HART-II Measured Motion in CFD

This presentation examines the use of HART-II measured rotor blade motion in computational fluid dynamics (CFD). Historically, comprehensive analyses were used for input to acoustic calculations. These analyses focused on lifting line aerodynamics and beam models. However, there is a a need to evolve lifting line aerodynamics to first principles, notably the use of CFD instead of lifting line. The current analysis focuses on CFD and computational structural dynamics (CSD) coupling. Beam models are still very good (CSD is typically from comprehensive analysis), but generally CFD replaced aerodynamics in comprehensive analysis. This presentation examines both CFD and CSD individually and includes predictions using measured motion as well as predictions using measured motion versus coupled motion and calculations of "correct" airloads, noise and vibration.

Boyd, D. Douglas, Jr.↗

DOE OSTI2023

Turbulence Modeling with Nek5000/RS, SOD2D and Alya

We present Validation and Verification (V & V) study of two high-order spectral element method (SEM) based Computational Fluid Dyanimcs (CFD) codes that we will utilize for work on turbulence modeling: Nek5000 and SOD2D. While the former solves the incompressible form of Navier-Stokes equation, the latter works with compressible set of equations and uses an entropy-viscosity formulation to account for the discontinuities for high Mach number flows. We demonstrate the accuracy of these codes for two benchmark sub-sonic turbulent flows: periodic channel and pipe flow, by carrying out first and second order statistical analysis, grid convergence and turbulent structure analysis using wall-resolved large eddy simulations (WRLES). Later, we report on the implementation and testing of various wall-modeling strategies for large eddy simulation of turbulent flows in Nek5000. These include both classical log-law based and decision-tree based machine-learning models. Accuracy of these closure strategies are analyzed and necessary future work is outlined.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Potassium carbonate decomposition modeling within rotating detonation engines for direct power extraction applications

We report that recently, there has been a significant interest in detonation-based combustion systems, such as rotating detonation engines (RDE), due to potential performance advantages in propulsion and energy applications. This includes improvements in thermodynamic efficiency, inherently high gas velocities, and the possibility of achieving an increase in total pressure (i.e. pressure gain) through the combustion process. Due to approximating a constant volume combustion process, these devices also tend to exhibit extremely high local gas temperatures relative to a comparable constant pressure combustor. These advantages overlap well with desirable performance characteristics of direct power extraction (DPE) technologies, such as a magnetohydrodynamic (MHD) generator. Typically, in DPE systems hot combustion products are seeded with an easily ionizable material such as potassium carbonate (K 2 CO 3 ) in order to boost the electrical conductivity. However, due to the short gas residence times within an RDE, it was unclear whether forming an electrically conductive combustion plasma would be feasible for integration with a downstream MHD generator. A model is presented which describes the heating, decomposition, and ionization of solid K 2 CO 3 particles and aqueous solutions of K 2 CO 3 in water, for a given initial particle diameter. This model was combined with available computational fluid dynamics (CFD) data for an oxygen-methane RDE in a one-way coupled Eulerian-Lagrangian framework to predict particle trajectories and the corresponding heating, decomposition, and ionization histories. Electrical conductivities were computed using a previously developed model, and a method was proposed to determine an equivalent average electrical conductivity. Results show that particle sizes below ~30 µm are able to fully decompose before reaching the exit of the RDE. While the one-way coupled nature of the simulations precluded rigorous evaluation of the effects of seed material on detonability, a substantial temperature reduction is expected at the detonation wave front due to heating and decomposition. A preliminary comparison is presented between the RDE and an equivalent constant pressure adiabatic combustor, showing a potential performance advantage for the RDE.

42 ENGINEERING↗

NASA NTRS2018

PIV and Rotational Raman-Based Temperature Measurements for CFD Validation in a Single Injector Cooling Flow

Film cooling is used in a wide variety of engineering applications for protection of surfaces from hot or combusting gases. The design of more efficient thin film cooling geometries/configurations could be facilitated by an ability to accurately model and predict the effectiveness of current designs using computational fluid dynamics (CFD) code predictions. Hence, a benchmark set of flow field property data were obtained for use in assessing current CFD capabilities and for development of better turbulence models. Both Particle Image Velocimetry (PIV) and spontaneous rotational Raman scattering (SRS) spectroscopy were used to acquire high quality, spatially-resolved measurements of the mean velocity, turbulence intensity and also the mean temperature and normalized root mean square (rms) temperatures in a single injector cooling flow arrangement. In addition to flowfield measurements, thermocouple measurements on the plate surface enabled estimates of the film effectiveness. Raman spectra in air were obtained across a matrix of radial and axial locations downstream from a 68.07 mm square nozzle blowing heated air over a range of temperatures and Mach numbers, across a 30.48 cm long plate equipped with a single injector cooling hole. In addition, both centerline streamwise 2-component PIV and cross-stream 3-component Stereo PIV data at 15 axial stations were collected in the same flows. The velocity and temperature data were then compared against Wind-US CFD code predictions for the same flow conditions. The results of this and planned follow-on studies will support NASA's development and assessment of turbulence models for heated flows.

Raman Scattering↗

NASA NTRS2018

PIV and Rotational Raman-Based Temperature Measurements for CFD Validation in a Single Injector Cooling Flow

Film cooling is used in a wide variety of engineering applications for protection of surfaces from hot or combusting gases. The design of more efficient thin film cooling geometries/configurations could be facilitated by an ability to accurately model and predict the effectiveness of current designs using computational fluid dynamics (CFD) code predictions. Hence, a benchmark set of flow field property data were obtained for use in assessing current CFD capabilities and for development of better turbulence models. Both Particle Image Velocimetry (PIV) and spontaneous rotational Raman scattering (SRS) spectroscopy were used to acquire high quality, spatially-resolved measurements of the mean velocity, turbulence intensity and also the mean temperature and normalized root mean square (rms) temperatures in a single injector cooling flow arrangement. In addition to flowfield measurements, thermocouple measurements on the plate surface enabled estimates of the film effectiveness. Raman spectra in air were obtained across a matrix of radial and axial locations downstream from a 68.07 mm square nozzle blowing heated air over a range of temperatures and Mach numbers, across a 30.48cm long plate equipped with a single injector cooling hole. In addition, both centerline streamwise 2-component PIV and cross-stream 3-component Stereo PIV data at 15 axial stations were collected in the same flows. The velocity and temperature data were then compared against Wind-US CFD code predictions for the same flow conditions. The results of this and planned follow-on studies will support NASA's development and assessment of turbulence models for heated flows.

Wernet, Mark P.↗

NASA NTRS2018

PIV and Rotational Raman-Based Temperature Measurements for CFD Validation in a Single Injector Cooling Flow

Film cooling is used in a wide variety of engineering applications for protection of surfaces from hot or combusting gases. The design of more efficient thin film cooling geometries/configurations could be facilitated by an ability to accurately model and predict the effectiveness of current designs using computational fluid dynamics (CFD) code predictions. Hence, a benchmark set of flow field property data were obtained for use in assessing current CFD capabilities and for development of better turbulence models. Both Particle Image Velocimetry (PIV) and spontaneous rotational Raman scattering (SRS) spectroscopy were used to acquire high quality, spatially-resolved measurements of the mean velocity, turbulence intensity and also the mean temperature and normalized root mean square (rms) temperatures in a single injector cooling flow arrangement. In addition to flowfield measurements, thermocouple measurements on the plate surface enabled estimates of the film effectiveness. Raman spectra in air were obtained across a matrix of radial and axial locations downstream from a 68.07 mm square nozzle blowing heated air over a range of temperatures and Mach numbers, across a 30.48cm long plate equipped with a single injector cooling hole. In addition, both centerline streamwise 2-component PIV and cross-stream 3-component Stereo PIV data at 15 axial stations were collected in the same flows. The velocity and temperature data were then compared against Wind-US CFD code predictions for the same flow conditions. The results of this and planned follow-on studies will support NASA's development and assessment of turbulence models for heated flows.

Turbulence Modeling↗

NASA NTRS2019

PIV and Rotational Raman-Based Temperature Measurements for CFD Validation of a Perforated Plate Cooling Flow: Part I

Film cooling is used in a wide variety of engineering applications for protection of surfaces from hot or combusting gases. The design of more efficient film cooling geometries/configurations could be facilitated by an ability to accurately model and predict the effectiveness of current designs using computational fluid dynamics (CFD) code predictions. Hence, a benchmark set of flow field property data were obtained for use in assessing current CFD capabilities and for development of better modeling approaches for these turbulent flow fields where accurate calculation of turbulent heat flux is important. Both Particle Image Velocimetry (PIV) and spontaneous rotational Raman scattering (SRS) spectroscopy were used to acquire high quality, spatially-resolved measurements of the mean velocity, turbulence intensity as well as the mean temperature and root mean square (rms) temperatures in a film cooling flow field. In addition to off-body flow field measurements, infrared thermography (IR) and thermocouple measurements on the plate surface enabled estimates of the film effectiveness. Raman spectra in air were obtained across a matrix of axial locations downstream from a 68.07 mm square nozzle blowing heated air over a range of temperatures (up to TR = 2.7) and Mach numbers (up to M0.9), across a 30.48 cm long plate equipped with three patches of 45 small (~1 mm) diameter cooling holes arranged in a staggered configuration. In addition, both centerline streamwise 2-component PIV and cross-stream 3-component Stereo PIV data at 14 axial stations were collected in the same flows. Only a subset of the data collected in the test program is included in this Part I report and are available from the NASA STI office. The final portion of the data will be published in a future report, Part II, along with CFD predictions of the complex cooling film flow.

IR camera Measurments↗

NASA NTRS2019

PIV and Rotational Raman-Based Temperature Measurements for CFD Validation of a Perforated Plate Cooling Flow: Part 1

Raman Scattering↗

NASA NTRS2014

Combined Uncertainty and A-Posteriori Error Bound Estimates for General CFD Calculations: Theory and Software Implementation

This workshop presentation discusses the design and implementation of numerical methods for the quantification of statistical uncertainty, including a-posteriori error bounds, for output quantities computed using CFD methods. Hydrodynamic realizations often contain numerical error arising from finite-dimensional approximation (e.g. numerical methods using grids, basis functions, particles) and statistical uncertainty arising from incomplete information and/or statistical characterization of model parameters and random fields. The first task at hand is to derive formal error bounds for statistics given realizations containing finite-dimensional numerical error [1]. The error in computed output statistics contains contributions from both realization error and the error resulting from the calculation of statistics integrals using a numerical method. A second task is to devise computable a-posteriori error bounds by numerically approximating all terms arising in the error bound estimates. For the same reason that CFD calculations including error bounds but omitting uncertainty modeling are only of limited value, CFD calculations including uncertainty modeling but omitting error bounds are only of limited value. To gain maximum value from CFD calculations, a general software package for uncertainty quantification with quantified error bounds has been developed at NASA. The package provides implementations for a suite of numerical methods used in uncertainty quantification: Dense tensorization basis methods [3] and a subscale recovery variant [1] for non-smooth data, Sparse tensorization methods[2] utilizing node-nested hierarchies, Sampling methods[4] for high-dimensional random variable spaces.

CFD↗

DOE OSTI2021

A Scoping Study to Develop a Computational Fluid Dynamics Based Model to Predict Radiological Materials Packaging Temperatures within a Generic Staging Building

The objective of this work is to perform computational fluid dynamics simulations of a ventilated radiological-material-package staging building to determine the effect of including or excluding a variety of physical effects and computational methods on both the predicted package temperatures and required computational resources. The generic building contains 640 drum-packages containing heat-generating, radiological material supported on four racks that are eight levels tall. Additionally, there is a forced ventilation system, lighting, and insulated walls. Computational models were constructed that included or excluded (a) shelving, (b) effects of unsteadiness, and (c) radiation heat transfer. Simulations with each drum modeled separately were compared to simpler simulations with sets of four drums represented by an equivalent box-package. These models employed between 106 to 107 elements. Steady state simulations predicted package temperatures that were within 0.1°C of simulations that included transient effects and required only one-eighth the computational resources. Calculations that excluded shelving predicted temperatures within 0.6°C of simulations that included shelving and required one-fourth the computational resources. Excluding radiation heat transfer systematically increased temperatures by around 1.5°C but reduced computational resources by a factor of four. The box-package model reduced the computational resources by a factor of 3, but systematically predicted higher temperatures by around 1.1°C. These results will be used to develop an efficient computational fluid dynamics model to assess the ability of different staging building designs to prevent the temperature of package components from exceeding specified limits.

42 ENGINEERING↗

NASA NTRS2015

Characterizing Droplet Formation from Non-Linear Slosh in a Propellant Tank

The Fluid Dynamics Branch (ER42) at the Marshall Space Flight Center (MSFC) was tasked with characterizing the formation and evolution of liquid droplets resulting from nonlinear propellant slosh in a storage tank. Lateral excitation of propellant tanks can produce high amplitude nonlinear slosh waves through large amplitude excitations and or excitation frequencies near a resonance frequency of the tank. The high amplitude slosh waves become breaking waves upon attaining a certain amplitude or encountering a contracting geometry such as the upper dome section of a spherical tank. Inherent perturbations in the thinning regions of breaking waves result in alternating regions of high and low pressure within the fluid. Droplets form once the force from the local pressure differential becomes larger than the force maintaining the fluid interface shape due to surface tension. Droplets released from breaking waves in a pressurized tank may lead to ullage collapse given the appropriate conditions due to the increased liquid surface area and thus heat transfer between the fluids. The goal of this project is to create an engineering model that describes droplet formation as a function of propellant slosh for use in the evaluation of ullage collapse during a sloshing event. The Volume of Fluid (VOF) model in the production level Computational Fluid Dynamics (CFD) code Loci-Stream was used to predict droplet formation from breaking waves with realistic surface tension characteristics. Various excitation frequencies and amplitudes were investigated at multiple fill levels for a single storage tank to create the engineering model of droplet formation from lateral propellant slosh.

Brodnick, Jacob↗

DOE OSTI2022

Prediction of local concentration fields in porous media with chemical reaction using a multi scale convolutional neural network

The study of solute transport in porous media is of interest in many chemical engineering systems. Some example applications include packed bed catalytic reactors, filtration devices, and batteries. The pore scale modeling of these systems is time consuming and may require large computing resources, for this reason computational fluid dynamics (CFD) simulations are not practical if a large number of simulations is required, like in multiscale modeling, where a model at a large scale calls for pore scale simulations. It has been shown that neural networks can be trained with a dataset of flow simulations and then predict fields orders of magnitude faster, and with less computational resources, in new domains. However, it is crucial to provide the neural network with an effective description of the domain and the undergoing operating conditions to be able to train models that generalize accurately in unseen samples. Therefore, research is needed to employ neural networks in new complex systems. The appropriate training of a network for predicting coupled flow and solute transport processes is an outstanding problem due to the complex interplay between geometry and operating conditions. In this work, we train a multi scale convolutional neural network (MSNet) with a diverse dataset of simulations of transport and chemical reaction in porous media to predict the local concentration fields in images of porous media. Our dataset contains a wide diversity of sphere pack arrangements under different operating conditions (Péclet and Reynolds numbers). Further, we train a robust model by employing different input descriptors that represent the medium and the different operating conditions of each system. Our trained model is able to provide nearly instantaneous predictions, compared to around twenty hours of the CFD workflow, with less than 3.5% error on new geometries and transport conditions. Thus the model could be easily integrated in a multiscale workflow where fast response is needed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗