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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

An Investigation of the Effects of Volume Fraction on Drag Coefficient of Non-Spherical Particles Using PR-DNS

Prediction of the drag coefficient is required in gas-particle multiphase flow modeling and simulation. Experimental data and correlations on the fixed-bed system of spherical particles with high volume fractions for various possible arrangements are available in the literature. However, the effect of volume fraction on the drag coefficient of non-spherical particles is not well studied. In solving the momentum equation, the volume fraction plays a vital role in determining the flow resistances. In this paper, we study the impact of volume fraction in the range of 0.069 to 0.65 on the drag coefficient using the computational fluid dynamics (CFD) simulation of air for Reynold number in the range of 10 to 10000 using particle resolved direct numerical solution (PR-DNS). Regular non-spherical particles such as a cube, tetrahedron, and spheroids are used in this study since their single particle’s drag coefficient data are available in the literature for comparison. For this work, the simulations are carried out in the Ansys Fluent using polyhedral mesh, which consumes significantly less computational time and power. The study showed the sphericity and volume fraction have significant impact on the bed pressure drop and average drag coefficient of the particles in the bed especially in high Reynolds number regime. The bed of the spheroid experiences the lowest drag being the most streamlined particle, and the particles with the edges result in a large drag coefficient due to flow separation at the discontinuity. The vector plots verify this behavior where large wake regions are observed behind the tetrahedron particle.

Mahyawansi, Pratik↗
DOE OSTI2021

A two-equation soot-in-flamelet modeling approach applied under Spray A conditions

Soot production (including formation and oxidation) is studied in the transient, high-pressure and turbulent n-dodecane Spray A flames from the Engine Combustion Network (ECN) using computational fluid dynamics (CFD) simulations. A two-equation soot-in-flamelet modeling approach is applied within the framework of the Unsteady Flamelet Progress Variable (UFPV) model and results are validated against experimental data. Equations for soot mass fraction and soot number density derived in the mixture fraction space are solved in the context of detailed flamelet calculations. Source terms for the different steps in the soot chemistry are tabulated and incorporated in the flamelet manifold. For the reference condition, the modeling approach based on the tabulated flamelet manifold reduces the computational cost of a CFD calculation by approximately 40 times compared to a non-tabulated well-mixed (WM) modeling approach. The soot-in-flamelet approach is then extended to study the effect of ambient oxygen concentration, ambient mixture composition and ambient temperature on soot production. Results show that the modeling approach is able to capture the experimental trends for the soot volume fraction (SVF) with good quantitative agreement, especially in the soot ramp-up region.

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

Reduced-order CFD modeling of cryogenic hydrogen isotope extrusion for pellet fueling

This study presents a reduced-order model (ROM) for computational fluid dynamics (CFD) simulations of cryogenic hydrogen isotope extrusions, focusing on protium (H₂) and deuterium (D₂) piston extruders. Using a 2D axisymmetric ROM in ANSYS-Polyflow, significant computational savings were achieved (runtime reduced from 9∼24 h to 3∼5 min), with extrusion force discrepancies between the 2D ROM and 3D models being on the order of 1%. Parametric studies identified optimal cutoff shear rates in the viscosity model (0.01/s for H₂ and 0.001/s for D₂), providing recommendations for future simulations. Finally, a comprehensive comparison of ROM results with experimental data was performed across varying geometries, cryogenic materials, temperatures, extrusion lengths, and piston velocities. Predictions at low extrusion temperatures met the objective of providing quick and efficient solutions with an acceptable extrusion force error of approximately 10% or less, validating the effectiveness of the 2D ROM approach. However, at high temperatures closer to the triple point, extrusion force error grows, which necessitates developing an improved model that accounts for temperature effects, e.g. melting. Nevertheless, the findings still represent a significant improvement in efficiency of CFD modeling of cryogenic hydrogenic extrusion. The ROM framework can also be extended to tritium (T2) and screw extruders, which will ultimately provide a fast and effective tool for optimizing pellet injector design for ITER and future reactor systems.

Fan, Joy [ORNL] (ORCID:0000000229751735)↗
NASA NTRS2002

Load Balancing Strategies for Multi-Block Overset Grid Applications

The multi-block overset grid method is a powerful technique for high-fidelity computational fluid dynamics (CFD) simulations about complex aerospace configurations. The solution process uses a grid system that discretizes the problem domain by using separately generated but overlapping structured grids that periodically update and exchange boundary information through interpolation. For efficient high performance computations of large-scale realistic applications using this methodology, the individual grids must be properly partitioned among the parallel processors. Overall performance, therefore, largely depends on the quality of load balancing. In this paper, we present three different load balancing strategies far overset grids and analyze their effects on the parallel efficiency of a Navier-Stokes CFD application running on an SGI Origin2000 machine.

Djomehri, M. Jahed↗
NASA NTRS2017

Progress Towards Understanding Fan Inlet Implications of Top-Mounted Propulsion

Computational fluid dynamic (CFD) simulations were performed on an N+2 commercial supersonic transport aircraft design that featured top-mounted propulsion. The simulations were run at take-off conditions at both 0 degrees and 8 degrees angle of attack. The results showed little separation in around the inlets with inlet performance on par with an under-the-wing configuration. The next step will be to take these results and determine the acoustic impact of the top-mounted propulsion system.

Computational Fluid Dynamics (CFD)↗
NASA NTRS2009

Lessons Learned from Numerical Simulations of the F-16XL Aircraft at Flight Conditions

Nine groups participating in the Cranked Arrow Wing Aerodynamics Project International (CAWAPI) project have contributed steady and unsteady viscous simulations of a full-scale, semi-span model of the F-16XL aircraft. Three different categories of flight Reynolds/Mach number combinations were computed and compared with flight-test measurements for the purpose of code validation and improved understanding of the flight physics. Steady-state simulations are done with several turbulence models of different complexity with no topology information required and which overcome Boussinesq-assumption problems in vortical flows. Detached-eddy simulation (DES) and its successor delayed detached-eddy simulation (DDES) have been used to compute the time accurate flow development. Common structured and unstructured grids as well as individually-adapted unstructured grids were used. Although discrepancies are observed in the comparisons, overall reasonable agreement is demonstrated for surface pressure distribution, local skin friction and boundary velocity profiles at subsonic speeds. The physical modeling, steady or unsteady, and the grid resolution both contribute to the discrepancies observed in the comparisons with flight data, but at this time it cannot be determined how much each part contributes to the whole. Overall it can be said that the technology readiness of CFD-simulation technology for the study of vehicle performance has matured since 2001 such that it can be used today with a reasonable level of confidence for complex configurations.

Rizzi, Arthur↗
NASA NTRS2018

Performance of Supersonic Parachutes Behind Slender Bodies

NASA’s ASPIRE (Advanced Supersonic Parachute Inflation Research Experiments) project is investigating the supersonic deployment, inflation and aerodynamics of full-scale disk-gap-band (DGB) parachutes. The first two flight tests were carried out in October 2017 and March 2018, while a third test is planned for the fall of 2018. In these tests, Mars-relevant conditions are achieved by deploying the parachutes at high altitudes over Earth using a sounding rocket test platform. As a result, the parachute is deployed behind a slender body (roughly 1/6-th the diameter of the capsule that will use this parachute for descent at Mars). Because there is limited flight and experimental data for supersonic DGBs behind slender bodies, the development of the parachute aerodynamic models was informed by CFD simulations of both the leading body wake and the parachute canopy. This presentation will describe the development of the pre-flight parachute aerodynamic models and compare pre-flight predictions with the reconstructed performance of the parachute during the flight tests. Specific attention will be paid to the differences in parachute performance behind blunt and slender bodies.

CFD↗
NASA NTRS2009

What was Learned from Numerical Simulations of F-16XL (CAWAPI) at Flight Conditions

Nine groups participating in the CAWAPI project have contributed steady and unsteady viscous simulations of a full-scale, semi-span model of the F-16XL aircraft. Three different categories of flight Reynolds/Mach number combinations were computed and compared with flight-test measurements for the purpose of code validation and improved understanding of the flight physics. Steady-state simulations are done with several turbulence models of different complexity with no topology information required and which overcome Boussinesq-assumption problems in vortical flows. Detached-eddy simulation (DES) and its successor delayed detached-eddy simulation (DDES) have been used to compute the time accurate flow development. Common structured and unstructured grids as well as individually-adapted unstructured grids were used. Although discrepancies are observed in the comparisons, overall reasonable agreement is demonstrated for surface pressure distribution, local skin friction and boundary velocity profiles at subsonic speeds. The physical modeling, be it steady or unsteady flow, and the grid resolution both contribute to the discrepancies observed in the comparisons with flight data, but at this time it cannot be determined how much each part contributes to the whole. Overall it can be said that the technology readiness of CFD-simulation technology for the study of vehicle performance has matured since 2001 such that it can be used today with a reasonable level of confidence for complex configurations.

Rizzi, Arthur↗
NASA NTRS2022

Recommended Direct Simulation Monte Carlo Collision Model Parameters for Reacting Methane Flows

While hydrocarbon combustion and gas-surface interactions have been simulated for a variety of applications in continuum solvers, such processes have been rarely investigated in Direct Simulation Monte Carlo (DSMC). Thus Variable Hard/Soft Sphere (VHS/VSS)collision parameters are not available for most molecules of interest in the decomposition of methane, a significant pyrolysis gas product and potential atmospheric species for Titan re-entries. As these properties are essential for accurate modeling of gas transport and thermochemical effects, a method must be devised to compute them for use in DSMC. In this study, the methodology outlined in Stephani et al is utilized to compute the VHS/VSS parameters by fitting the Ω(1,1)and Ω(2,2) collision integrals [1]. These collision integrals are typically computed from intermolecular potentials and are used in CFD simulations to calculate transport properties such as viscosity, thermal conductivity, and diffusion. This approach ensures consistency between micro-scale collisions and macro-scale transport properties. To produce the necessary collision integrals needed for the fitting routine, we utilize the method described by Laricchiuta et al, where collision integrals are computed based on the constituent particles’ polarizability and effective electron number [2].Preliminary computations of pure-species Ω(1,1)and Ω(2,2)collision integrals produced from the coupled approach of Stephani et al and Laricchiuta et al for methane show good agreement with those computed via available ab-initio Potential Energy Surface Computations, as seen in Figure 1. Above 300K, no more than 12% deviation is observed between the two results over the temperature range of the fitted PES data [3].Current results of this study include full collision-specificVSS transport properties for a complex-species mixture involving major byproducts of pure methane decomposition. Targeted final results will include a comprehensive list of VSS parameters relevant toDSMC simulations of Titan atmospheric entry, based on the reduced kinetic model of Savajano et al [4].

M R Gosma↗
NASA NTRS2017

Lessons Learned from Inlet Integration Analysis of NASA's Low Boom Flight Demonstrator

In 2016, NASA's Aeronautics Research Mission Directorate announced the New Aviation Horizons Initiative with a goal of designing/building several X-Planes, including a Low Boom Flight Demonstrator (LBFD). That same year, NASA awarded a contract to Lockheed Martin (LM) to advance the LBFD concept through preliminary design. Several configurations of the LBFD aircraft were analyzed by both LM engineers and NASA researchers. This presentation focuses on some of the CFD simulations that were run by NASA Glenn researchers. NASA's FUN3D V13.1 code was used for all adjoint-based grid refinement studies and Spalart-Allmaras turbulence model was used during adaptation. It was found that adjoint-based grid adaptation did not accurately capture inlet performance for high speed top-aft-mounted propulsion.

Adjoint-based Grid Adaptatio↗
NASA NTRS2019

Space Launch System Booster Separation Supersonic Powered Testing with Surface and Off-Body Measurements

A wind tunnel test was run in the NASA Langley Unitary Plan Wind Tunnel simulating the separation of the two solid rocket boosters (SRB) from the core stage of the NASA Space Launch System (SLS). The test was run on a 0.9% scale model of the SLS Block 1B Cargo (27005) configuration and the SLS Block 1B Crew (28005) configuration at a Mach of 4.0. High pressure air was used to simulate plumes from the booster separation motors located at the nose and aft skirt of the two boosters. Force and moment data were taken on both SRBs and on the core stage. Schlieren still photos and video were recorded throughout testing. A set of points were acquired using Cross-correlation Doppler Global Velocimetry (CCDGV) readings to get 3 component velocity measurements between the core and the left-hand SRB. The CCDGV laser was utilized to record flow visualization in the same location, between the core and the left-hand SRB. Pressure Sensitive Paint data were taken on a separate set of runs. Computational Fluid Dynamics (CFD) runs were computed on a subset of the wind tunnel data points for comparison. A combination of the force/moment, CCDGV and Pressure Sensitive Paint (PSP) data (as well as schlieren images) at the CFD-specified test conditions will be used te the CFD simulations that will be used to build an SLS booster separation database flight conditions.

Winski, Courtney S.↗
NASA NTRSDate not supplied

Space Launch System Booster Separation Supersonic Powered Testing with Surface and Off-body Measurements

A wind tunnel test was run in the NASA Langley Unitary Plan Wind Tunnel simulating the separation of the two solid rocket boosters (SRB) from the core stage of the NASA Space LaunchSystem (SLS). The test was run on a 0.9% scale model of the SLS Block 1B Cargo (27005) configuration and the SLS Block 1B Crew (28005) configuration at a Mach of 4.0. The SRBs weremounted separately from the core, and could be pitched, yawed, and translated in the X, Y, andZ directions relative to the core. High pressure air was used to simulate plumes from the booster separation motors located at the nose and aft skirt of the two boosters. Force and moment data were taken on both SRBs and on the core stage. Schlieren still photos and video were recorded throughout testing. A set of points were acquired using Cross-correlation Doppler Global Velocimetry (CCDGV) readings to get 3 component velocity measurements between the core and the left-hand SRB. The CCDGV laser was utilized to record flow visualization in the same location, between the core and the left-hand SRB. Pressure Sensitive Paint (PSP) data were taken on a separate set of runs. Computational Fluid Dynamics (CFD) runs were computed on a subset of the wind tunnel data points for comparison. A combinationof the force and moment, CCDGV, and PSP data (as well as schlieren images) at the CFD-specified test conditions will be used to validate the CFD simulations that will be used to build an SLS booster separation database at flight conditions.

SLS↗
DOE OSTI2020

Data-driven selection of stiff chemistry ODE solver in operator-splitting schemes

Most computational fluid dynamics simulations of practical combustion applications employ operator-splitting schemes, where chemistry and transport are separated and integrated with distinct numerical methods. The changes in composition due to chemistry are evaluated by solving ordinary differential equations (ODE) in each cell of the computational domain, which typically dominates the computational cost when detailed chemistry is considered. In this work, a data-driven approach for the selection of chemistry ODE solvers in operator-splitting schemes is presented. Neural networks are used to predict the ODE solvers CPU times and errors for a given thermochemical state. This allows the selection of an optimal ODE solver on a cell-by-cell, timestep-by-timestep basis. The models are trained using a wide set of thermochemical states generated through partially-stirred reactors and flames simulations. The methodology is validated by quantifying the prediction errors, the classification accuracy, and the computational speedup. The model predicts the optimal ODE solver for 70 to 95% of the validation cases and decreases the computional cost by a factor of 3 or more. The generalizability of the methodology to different chemical mechanisms and different fuels is assessed and it is shown that the model’s performance is only slightly degraded and its applicability is significantly enhanced if the inputs to the neural networks are restricted to a small set of thermochemical state variables present in most chemical mechanisms. In conclusion, the models are used in an homogeneous reactor case and a multi-dimensional CFD simulation of a diesel spray at high pressure where a speedup of more than 3 is achieved.

42 ENGINEERING↗
NASA NTRS2011

J-2X Fuel Pump Impeller Seal Simulations

The J-2X engine was originally designed for the upper stage of the previously cancelled Crew Launch Vehicle. Although the Crew Launch Vehicle was cancelled the J-2X engine, which is currently undergoing hot-fire testing, may be used on future programs. The J-2X engine is a direct descendent of the J-2 engine which powered the upper stage during the Apollo program. Many changes including a thrust increase from 230K to 294K lbf have been implemented in this engine. The rotor-dynamic stability of the fuel turbopump is highly dependent on the tangential velocity of the fluid as it enters the the front face impeller seal. Rotor-dynamic analysis predicts that a much lower tangential velocity will be required for stability than was needed for previous engines. The geometry at the seal entrance for this engine is very complex and vastly different than previous engines. In order to better determine the fluid dynamics and tangential velocity in this seal several CFD simulations were performed. The results of these simulations show that for this seal geometry a great reduction in the tangential velocity is to be expected. The simulations also provided insight into methods that could be employed to drive the swirl velocity to near zero. Unsteady and time-averaged results of several simulations will be presented.

Schmauch, Preston B.↗
NASA NTRS2005

Using CFD as a Rocket Injector Design Tool: Recent Progress at Marshall Space Flight Center

New programs are forcing American propulsion system designers into unfamiliar territory. For instance, industry s answer to the cost and reliability goals set out by the Next Generation Launch Technology Program are engine concepts based on the Oxygen- Rich Staged Combustion Cycle. Historical injector design tools are not well suited for this new task. The empirical correlations do not apply directly to the injector concepts associated with the ORSC cycle. These legacy tools focus primarily on performance with environment evaluation a secondary objective. Additionally, the environmental capability of these tools is usually one-dimensional while the actual environments are at least two- and often three-dimensional. CFD has the potential to calculate performance and multi-dimensional environments but its use in the injector design process has been retarded by long solution turnaround times and insufficient demonstrated accuracy. This paper has documented the parallel paths of program support and technology development currently employed at Marshall Space Flight Center in an effort to move CFD to the forefront of injector design. MSFC has established a long-term goal for use of CFD for combustion devices design. The work on injector design is the heart of that vision and the Combustion Devices CFD Simulation Capability Roadmap that focuses the vision. The SRL concept, combining solution fidelity, robustness and accuracy, has been established as a quantitative gauge of current and desired capability. Three examples of current injector analysis for program support have been presented and discussed. These examples are used to establish the current capability at MSFC for these problems. Shortcomings identified from this experience are being used as inputs to the Roadmap process. The SRL evaluation identified lack of demonstrated solution accuracy as a major issue. Accordingly, the MSFC view of code validation and current MSFC-funded validation efforts were discussed in some detail. The objectives of each effort were noted. Issues relative to code validation for injector design were discussed in some detail. The requirement for CFD support during the design of the experiment was noted and discussed in terms of instrumentation placement and experimental rig uncertainty. In conclusion, MSFC has made significant progress in the last two years in advancing CFD toward the goal of application to injector design. A parallel effort focused on program support and technology development via the SCIT Task have enabled the progress.

Tucker, Kevin↗
DOE OSTI2024

Design, Analysis, and Experimental Testing of Hydrogen Lean Direct Injection Nozzles at Elevated Pressure

Abstract There are many challenges of commissioning a hydrogen combustor into future gas turbine engines; especially regarding achieving emissions goals. Previously, Escudero et al. and Tran et al. conducted a study to adapt the liquid fuel Lean Direct Injection (LDI) concept from Jet-A to gaseous natural gas-hydrogen blends and pure hydrogen [1], [2]. Experimental data was collected at atmospheric conditions using a Box Behnken design of experiments. The design of experiments suggested that biasing the air split in favor of the inner air circuit and increasing the swirl strength of this inner air passage resulted in improved NOx emissions, while the inverse was true for stability, which was quantified by studying the lean blowoff point (LBO) [1], [2]. The trends revealed by the original experiment [1], [2] provided a design direction for further iterations of the experimental hardware. The study presented herein describes the further investigation of such LDI injectors through experimental methods and computational fluid dynamic (CFD) simulations at atmospheric conditions, which were used to identify potential flow behaviors driving enhanced emissions performance. Further evaluation of select injectors from both studies was then conducted at elevated pressures up to 6 atmospheres. The results from both experiments are presented in this study, which include flame observations, emissions measurements, and operational challenges. NOx emissions results are reported on a volume basis in ppmvd corrected to 15% O2 and corrected for fuel. A predictive model for relating NOx emissions to test conditions at atmospheric conditions show high significance to adiabatic flame temperature while little to no significance to fuel composition for the best performing configurations. The results illustrate the connection between atmospheric testing and testing elevated pressures. The design direction indicated by the initial tests and CFD results in promising configurations for implementation into a Multi-point LDI array.

08 HYDROGEN↗
DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗
NASA NTRS2019

A Multidisciplinary Modeling Approach of Plant Gas Exchange in Reduced Gravity Environments

In-situ food production is a necessary step for human exploration of the solar system and requires a deep understanding of plant growth in reduced gravity environments. In particular, the lack of buoyancy-driven convection changes the gas exchange at the leaf surface, which decreases photosynthesis and transpiration rates, and ultimately biomass production. To understand the intricate relations between physical, chemical, and biochemical processes, the following methodology combines the development of a mechanistic model of plant growth in reduced gravity environments, computational fluid dynamics (CFD) simulations, and experiments in different time frames.The model presented here is a coupled mass and energy balance using the single round leaf assumption, including gravity as an entry parameter, and the leaf surface temperature as an output variable. Measures of the leaf surface temperature using infra-red cameras allow for a computation of the transpiration rate. This approach was followed to design a parabolic flight experiment, which performed 7 flights, and enabled data collection for model validation in different gravity and ventilation settings on a short time frame. Current measures of carbon assimilation and transpiration rate at the leaf and canopy level using an infra-red gas analyzer (Li-6800) in 1g lab conditions on several species will enable a validation on longer time frames and further calibration of the model. CFD studies both on the parabolic flight and on the lab experimental set-up allow the precise assessment of ventilation above the canopy and plants' leaves.Ultimately, this work will provide recommendations for the design of future plant growth hardware, especially on the lowest adequate ventilation for optimal plant growth in reduced gravity environments, as well as assessing biomass and oxygen production rates on planetary surfaces and space stations. This work was funded by CNES, CNRS, Clermont Auvergne Metropole, and NASA Space Biology through NASA postdoctoral program / USRA.

Poulet, Lucie↗