Results for “computational fluid dynamics simulation” — Page 42

NASA NTRSDate not supplied

A Green’s Function Sensor Fusion Approach for Evaluating Spacecraft Entry Heating From on-Board Thermal Instrumentation

During atmospheric entry, distributed thermal measurements are critical to enable the evaluation of heat loads on spacecraft thermal protection systems (TPS). In recent space exploration missions, Schmidt-Boelter-type heat flux gauges have been integrated into the TPS alongside conventional temperature measurement instrumentation to measure total (convective and radiative) and radiative heat transfer rates1,2. While direct heat flux sensors (HFS) are able to provide valuable information detailing the thermal loads experienced by spacecraft, the interpretation of these measurements in unsteady, convective environments requires a correction factor to account for local heating augmentations at the cold wall HFS surface2,3. Current efforts to estimate cold wall correction factors, and thus recover the hot wall TPS heat flux, rely on time-marching computational fluid dynamics (CFD) simulations2. Simulation-based methods are susceptible to large uncertainties, however, as they require estimations of vehicle trajectory, gas kinetics, wall catalysis models, and other flight conditions as input parameters2,4. Furthermore, CFD simulations are computationally expensive and cannot efficiently survey all possible entry scenarios, exacerbating the uncertainty of reconstructed hot wall heat flux values. These drawbacks motivate the development of alternative hot wall heat flux reconstruction methods that are not reliant on CFD-based correction factors.

Kenneth McAfee↗

NASA NTRSDate not supplied

A Green’s Function Sensor Fusion Approach for Evaluating Spacecraft Entry Heating From on-Board Thermal Instrumentation

During atmospheric entry, distributed thermal measurements are critical to enable the evaluation of heat loads on spacecraft thermal protection systems (TPS). In recent space exploration missions, Schmidt-Boelter-type heat flux gauges have been integrated into the TPS alongside conventional temperature measurement instrumentation to measure total (convective and radiative) and radiative heat transfer rates1,2. While direct heat flux sensors (HFS) are able to provide valuable information detailing the thermal loads experienced by spacecraft, the interpretation of these measurements in unsteady, convective environments requires a correction factor to account for local heating augmentations at the cold wall HFS surface2,3. Current efforts to estimate cold wall correction factors, and thus recover the hot wall TPS heat flux, rely on time-marching computational fluid dynamics (CFD) simulations2. Simulation-based methods are susceptible to large uncertainties, however, as they require estimations of vehicle trajectory, gas kinetics, wall catalysis models, and other flight conditions as input parameters2,4. Furthermore, CFD simulations are computationally expensive and cannot efficiently survey all possible entry scenarios, exacerbating the uncertainty of reconstructed hot wall heat flux values. These drawbacks motivate the development of alternative hot wall heat flux reconstruction methods that are not reliant on CFD-based correction factors.

Kenneth McAfee↗

NASA NTRS1989

A zonal Navier-Stokes methodology for flow simulation about a complete aircraft

The thin layer, Reynolds-averaged, Navier-Stokes equations are used to simulate the transonic viscous flow about the complete F-16A fighter aircraft. These computations demonstrate how computational fluid dynamics can be used to simulate turbulent viscous flow about realistic aircraft geometries. A zonal grid approach is used to provide adequate viscous grid clustering on all aircraft surfaces. Zonal grids extend inside the F-16A inlet and up to the compressor face while power on conditions are modeled by employing a zonal grid extending from the exhaust nozzle to the far field. Computations are compared with existing experimental data and are in fair agreement. Computations for the F-16A in side slip are also presented.

Flores, Jolen↗

NASA NTRS1992

Improved Depiction Of Computed Unsteady Flows

PLOT4D and STREAKER programs increase ease of production and utility of animated depictions of unsteady flows simulated numerically by advanced computational fluid dynamics. Programs developed to maximize information displayed while minimizing amount stored. PLOT4D, an extended version of PLOT3D, used to depict steady three-dimensional flows, incorporates features making it more computationally efficient than to visualize unsteady flow as succession of steady flows, animating frame by frame. PLOT4D postprocesses numerical output of fluid-dynamical computation into maps of flow quantities on two-coordinate surfaces depicting cut through flow field or conform to body surface. STREAKER generates streaklines, close approximations of patterns traced by streams of dye or smoke released into simulated flows. Streaklines defined via spatial succession of endpoints of pathlines of massless particles released into flow. STREAKER computes endpoints by fourth-order Runge-Kutta integration of velocity vectors.

Smith, Merritt H.↗

NASA NTRS2017

DSMC simulations of the Shuttle Plume Impingement Flight EXperiment(SPIFEX)

During orbital maneuvers and proximity operations, a spacecraft fires its thrusters inducing plume impingement loads, heating and contamination to itself and to any other nearby spacecraft. These thruster firings are generally modeled using a combination of Computational Fluid Dynamics (CFD) and DSMC simulations. The Shuttle Plume Impingement Flight EXperiment(SPIFEX) produced data that can be compared to a high fidelity simulation. Due to the size of the Shuttle thrusters this problem was too resource intensive to be solved with DSMC when the experiment flew in 1994.

Stewart, Benedicte↗

NASA NTRSDate not supplied

Validation of a Two-Phase CFD Model for Autogenous Pressurization and Expulsion

This paper presents a two-phase computational fluid dynamics (CFD) model for simulating autogenous pressurization and expulsion in a cryogenic liquid hydrogen storage tank, utilizing a volume-of-fluid (VOF) approach combined with the kinetics-based Schrage equation to accurately capture the phase boundary and interfacial heat, mass, and momentum transfer between vapor and liquid phases. The model is validated against experimental data from NASA's K-site tank facility, specifically experimental case 225, which involved autogenous pressurization and controlled expulsion of liquid hydrogen. Various turbulence models are evaluated to assess their influence on model accuracy. The CFD simulations successfully replicate key thermodynamic behaviors observed during the experiments, including pressure evolution, temperature profiles, and phase-change dynamics at the vapor-liquid interface. The predicted tank pressures and temperatures agree well with experimental data, whereas the pressurant mass predictions are within 16% of observed values. This study emphasizes the importance of selecting appropriate turbulence models to accurately simulate the complex flow and heat transfer phenomena during tank autogenous pressurization and expulsion. By improving the accuracy and reliability of CFD models for these processes, this research contributes to developing efficient cryogenic propellant management strategies for future space missions.

Computational Fluid Dynamics↗

NASA NTRSDate not supplied

Validation of a Two-Phase CFD Model for Autogenous Pressurization and Expulsion

This paper presents a two-phase computational fluid dynamics (CFD) model for simulating autogenous pressurization and expulsion in a cryogenic liquid hydrogen storage tank, utilizing a volume-of-fluid (VOF) approach combined with the kinetics-based Schrage equation to accurately capture the phase boundary and interfacial heat, mass, and momentum transfer between liquid and vapor phases. The model is validated against experimental data from NASA's K-site tank facility, specifically experimental case 225, which involved pressurization and controlled expulsion of liquid hydrogen. Various turbulence models are evaluated to assess their influence on model accuracy. The CFD simulations successfully replicate key thermodynamic behaviors observed during the experiments, including pressure evolution, temperature profiles, and phase-change dynamics at the vapor-liquid interface. The predicted tank pressures, temperatures, and pressurant requirements agree with experimental data, with pressurant mass predictions within 16\% of observed values. This study highlights the importance of selecting appropriate turbulence models to accurately simulate complex flow and heat transfer phenomena during tank pressurization and expulsion. By enhancing the accuracy and reliability of CFD models for liquid hydrogen under cryogenic conditions, this research contributes to developing efficient cryogenic propellant management strategies for future space missions.

Computational Fluid Dynamics↗

DOE OSTI2025

Verification and Validation of Open-Source Computational Fluid Dynamics Capabilities in MOOSE

The Multiphysics Object Oriented Simulation Environment (MOOSE) advances predictive modeling in nuclear research by simplifying the application design process and making modeling more efficient and accessible to researchers and the broader community. Because of MOOSE’s critical role in the scientific computation field, verification and validation are important to evaluate the reliability and credibility of MOOSE’s code, simulations, and models. To assess MOOSE’s Navier Stokes module and turbulence models, the following verification and validation University of Manchester ERCOFTAC cases were simulated in MOOSE, and the results were compared: Case 032 Turbulent Flow in a 2D Channel Case 079 Natural Convection in an Enclosed Cavity

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

NASA NTRS1993

Mathematical description of complex chemical kinetics and application to CFD modeling codes

A major effort in combustion research at the present time is devoted to the theoretical modeling of practical combustion systems. These include turbojet and ramjet air-breathing engines as well as ground-based gas-turbine power generating systems. The ability to use computational modeling extensively in designing these products not only saves time and money, but also helps designers meet the quite rigorous environmental standards that have been imposed on all combustion devices. The goal is to combine the very complex solution of the Navier-Stokes flow equations with realistic turbulence and heat-release models into a single computer code. Such a computational fluid-dynamic (CFD) code simulates the coupling of fluid mechanics with the chemistry of combustion to describe the practical devices. This paper will focus on the task of developing a simplified chemical model which can predict realistic heat-release rates as well as species composition profiles, and is also computationally rapid. We first discuss the mathematical techniques used to describe a complex, multistep fuel oxidation chemical reaction and develop a detailed mechanism for the process. We then show how this mechanism may be reduced and simplified to give an approximate model which adequately predicts heat release rates and a limited number of species composition profiles, but is computationally much faster than the original one. Only such a model can be incorporated into a CFD code without adding significantly to long computation times. Finally, we present some of the recent advances in the development of these simplified chemical mechanisms.

Bittker, D. A.↗

NASA NTRS1993

Mathematical Description of Complex Chemical Kinetics and Application to CFD Modeling Codes

A major effort in combustion research at the present time is devoted to the theoretical modeling of practical combustion systems. These include turbojet and ramjet air-breathing engines as well as ground-based gas-turbine power generating systems. The ability to use computational modeling extensively in designing these products not only saves time and money, but also helps designers meet the quite rigorous environmental standards that have been imposed on all combustion devices. The goal is to combine the very complex solution of the Navier-Stokes flow equations with realistic turbulence and heat-release models into a single computer code. Such a computational fluid-dynamic (CFD) code simulates the coupling of fluid mechanics with the chemistry of combustion to describe the practical devices. This paper will focus on the task of developing a simplified chemical model which can predict realistic heat-release rates as well as species composition profiles, and is also computationally rapid. We first discuss the mathematical techniques used to describe a complex, multistep fuel oxidation chemical reaction and develop a detailed mechanism for the process. We then show how this mechanism may be reduced and simplified to give an approximate model which adequately predicts heat release rates and a limited number of species composition profiles, but is computationally much faster than the original one. Only such a model can be incorporated into a CFD code without adding significantly to long computation times. Finally, we present some of the recent advances in the development of these simplified chemical mechanisms.

Bittker, D. A.↗

NASA NTRSDate not supplied

Modeling Self-Pressurization and Spray Bar Pressure Control of A Cryogenic Storage Tank in Normal Gravity

This paper presents computational fluid dynamics (CFD) models for simulating self-pressurization and spray-bar pressure control processes in a large-scale liquid hydrogen storage tank under normal gravity conditions. For self-pressurization, the model employs the kinetics-based Schrage equation alongside the volume-of-fluid (VOF) method to account for interfacial mass transfer. The CFD predictions of pressure and temperature are compared with experimental data from the Multipurpose Hydrogen Test Bed (MHTB) experiment, and the predicted interfacial mass transfer rates are also presented. A CFD model simulating pressure control using a spray bar has also been developed. An Eulerian-Lagrangian approach models the interactions between discrete droplets and the continuous ullage (vapor) phase. The spray model is coupled with the VOF method by tracking droplets in the ullage and removing them when they reach the liquid interface. The T-sat model calculates droplet-ullage heat and mass transfer, where droplets warm up to the saturation temperature corresponding to the ullage vapor pressure before evaporating while remaining at the saturation temperature. The evolution of tank pressure, vapor temperature, and liquid temperature predicted by the CFD model is validated against data from the MHTB spray-bar mixing experiment. Overall, the CFD models agree with experimental data, demonstrating their capability to simulate self-pressurization and pressure control processes in large-scale cryogenic storage tanks. These models can be valuable tools for designing and optimizing cryogenic fluid management systems in future applications.

Computational Fluid Dynamics↗

NASA NTRSDate not supplied

OVERFLOW Analysis of Supersonic Retropropulsion Testing on a Blunt Mars Entry Vehicle Concept

Supersonic Retropropulsion (SRP) flow over a Hypersonic Inflatable Aerodynamic Decelerator (HIAD) blunt-body vehicle was simulated using the OVERFLOW Computational Fluid Dynamics (CFD) solver. Two nozzle configurations were tested (1E and 1F) as a subset of the seven total configurations in the Descent Systems Study (DSS) testing campaign. A generalized, Adaptive Mesh Refinement (AMR) shock capturing and plume refinement technique was developed and calibrated for producing automatic and optimized grid systems unique to any given testing condition. Simulations were subjected to successive levels of grid refinement in order to determine solution independence from grid resolution. OVERFLOW Reynolds Averaged Navier-Stokes (RANS) solutions produced realistic SRP flow phenomena, including Mach disk normal shocks contained in high-thrust plumes and bow shock triple-points. Solutions for the 1E nozzle configuration were steady. Conversely, a subset of the 1F test conditions were unsteady and exhibited periodic, non-sinsoidal expansions and contractions of the streamwise-oriented plume, resulting in unsteady aerodynamic loads on the vehicle.

Supersonic Retropropulsion↗

NASA NTRS1995

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗

NASA NTRSDate not supplied

Modeling Self-Pressurization and Spray Bar Pressure Control of A Cryogenic Storage Tank in Normal Gravity

This paper presents computational fluid dynamics (CFD) models for simulating self-pressurization and spray-bar pressure control processes in a large-scale liquid hydrogen storage tank under normal gravity conditions. For self-pressurization, the model employs the kinetics-based Schrage equation alongside the volume-of-fluid (VOF) method to account for interfacial mass transfer. The CFD predictions of pressure and temperature are compared with experimental data from the Multipurpose Hydrogen Test Bed (MHTB) experiment, and the predicted interfacial mass transfer rates are also presented. A CFD model simulating pressure control using a spray bar has also been developed. An Eulerian-Lagrangian approach models the interactions between discrete droplets and the continuous ullage (vapor) phase. The spray model is coupled with the VOF method by tracking droplets in the ullage and removing them when they reach the liquid interface. The T-sat model calculates droplet-ullage heat and mass transfer, wherein droplets warm up to the saturation temperature corresponding to the ullage vapor pressure before evaporating while remaining at the saturation temperature. The evolution of tank pressure, vapor temperature, and liquid temperature predicted by the CFD model is validated against data from the MHTB spray-bar mixing experiment. Overall, the CFD models agree with experimental data, demonstrating their capability to simulate self-pressurization and pressure control processes in large-scale cryogenic storage tanks. These models can be valuable tools for designing and optimizing cryogenic fluid management systems in future applications.

Self-Pressurization↗

DOE OSTI2022

Enhancing $\phi$-sensitivity of ignition delay times through dilution of fuel-air mixture

The high $\phi$-sensitivity (η) of ignition delay time (τ IDT ) is one of the desirable fuel properties for the high-load extension of advanced compression ignition engine. Recent studies revealed the effectiveness of a high dilution rate (x D ) for enhancing η at the engine-relevant conditions. This study aims to quantify the effect of dilution on the η of isooctane using a combined experiment-simulation approach. The ignition delay of the isooctane/air mixture was measured with an Advanced Fuel Ignition Delay Analyzer (AFIDA) over the temperature range of 623 - 923 K at 10 bar pressure with global $\phi$ = 0.3 - 0.6, with and without 28.6% of additional N 2 dilution. For precise evaluation of experimental η, the facility effect of the AFIDA experiment was characterized with three-dimensional computational fluid dynamics (3-D CFD) simulation. The temperature in the combustion chamber from 3-D CFD indicated a substantial temporal dependency, varying up to ~52 K by charge-cooling of fuel injection and heat transfer from the wall. We introduced the dimensionless number θ(t) for characterizing the temporal profile of chamber temperature. Consideration of the facility effect using θ(t) resulted in better agreement between the experimental ..eta.. and zero-dimensional (0-D) kinetics simulation. The refined η were then further utilized to quantify the effectiveness of dilution to η The extent of η enhancement with dilution strategy was maximized at the low-temperature chemistry regime, increasing η by 77% with a 28.6% dilution rate. Further analysis on the dilution effect was carried out using 0-D kinetics simulation, revealing the critical dimensionless numbers relevant to the effectiveness of dilution to η enhancement. Here, this study is the first experiment-simulation combined research to quantify the effect of dilution on η, facilitating the kinetics model refinement for better reproduction of $\phi$-sensitivity.

33 ADVANCED PROPULSION SYSTEMS↗

NASA NTRS2003

Evaluation of a Second-Order Accurate Navier-Stokes Code for Detached Eddy Simulation Past a Circular Cylinder

We evaluate the applicability of a production computational fluid dynamics code for conducting detached eddy simulation for unsteady flows. A second-order accurate Navier-Stokes code developed at NASA Langley Research Center, known as TLNS3D, is used for these simulations. We focus our attention on high Reynolds number flow (Re = 5 x 10(sup 4) - 1.4 x 10(sup 5)) past a circular cylinder to simulate flows with large-scale separations. We consider two types of flow situations: one in which the flow at the separation point is laminar, and the other in which the flow is already turbulent when it detaches from the surface of the cylinder. Solutions are presented for two- and three-dimensional calculations using both the unsteady Reynolds-averaged Navier-Stokes paradigm and the detached eddy simulation treatment. All calculations use the standard Spalart-Allmaras turbulence model as the base model.

Vatsa, Veer N.↗

NASA NTRSDate not supplied

Numerical Simulations of an Axisymmetric Shock-Wave/Boundary-Layer Interaction

A series of Reynolds-averaged Navier-Stokes (RANS) simulations were performed of two Mach 2.5 axisymmetric shock-wave/boundary-layer interactions (SWBLIs) based off of the NASA Turbulent Computational Fluid Dynamics Validation Experiment (TCFDVE). The goals of the simulations were to (1) determine how well RANS can predict the flow physics of an axisymmetric SWBLI and (2) generate a new validation case for the NASA Turbulence Modeling Resource (TMR). The simulations were run with two different turbulence models and three different computational fluid dynamics (CFD) codes while looking at two different conical shock generator angles, one of which was shown to cause flow separation in the experiment. The results in this paper show that, as expected, the CFD simulations agreed well with each other and the experiment when SWBLI-induced separation was not present. When separation was present, differences in the predicted separation characteristics were noted. These included a stronger interaction and subsequently longer separation extent for the Menter Shear-Stress Transport k-ω (SST) turbulence model simulations compared to both the Spalart-Allmaras (SA) turbulence model simulations and the experiment.

Computational Fluid Dynamics↗

DOE OSTI2021

Design Optimization of A-Type Heat Exchangers based on High Performance, Non-Round Tubes

Air-to-refrigerant heat exchangers (HXs) have been the topic of considerable research as they are the fundamental heat transfer components of HVAC&R systems. For residential and commercial applications, indoor heat exchangers are often in A-type configurations to reduce air-conditioner system footprint while still delivering the required cooling or heating. Classical HX design practices typically assume a uniform frontal air velocity profile to facilitate the use of conventional tube-fin airside heat transfer and pressure drop correlations for performance prediction. However, it is well-reported throughout the literature that A-type HXs experience significant airflow maldistribution which can severely degrade HX performance. Additionally, recent advancements in simulation software such as Computational Fluid Dynamics (CFD), Finite Element Analysis (FEA) and optimization algorithms have led researchers to consider primary heat transfer surface optimization to achieve highly compact HXs which do not require fins. In this paper, an A-type HX which utilizes shape-optimized non-round tubes is optimized for minimum duct size and airside pressure drop. This is achieved through coupling of an experimentally validated finite volume HX model with automated 2D CFD simulations of airflow through the HX for airside thermal-hydraulic performance evaluation. Preliminary results show that the optimal designs deliver similar capacity while achieving 20% reductions in airside pressure drop and duct cross-section area. Additionally, the optimal design apex angles in this study were similar to those from previous A-Type HX optimization studies for conventional tube-fin HX geometries, suggesting that the apex angle for optimal A-Type configuration performance may be independent of tube geometry.

42 ENGINEERING↗