Search NASA⌕ Search

SEARCH · Search NASA

Results for “Flow Simulator”

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 397 records · Page 22

Direct Simulation of Reentry Flows with Ionization

The Direct Simulation Monte Carlo (DSMC) method is applied in this paper to the study of rarefied, hypersonic, reentry flows. The assumptions and simplifications involved with the treatment of ionization, free electrons and the electric field are investigated. A new method is presented for the calculation of the electric field and handling of charged particles with DSMC. In addition, a two-step model for electron impact ionization is implemented. The flow field representing a 10 km/sec shock at an altitude of 65 km is calculated. The effects of the new modeling techniques on the calculation results are presented and discussed.

Carlson, Ann B.↗

Numerical simulation of viscous supersonic flow over a generic fighter configuration

A procedure is presented, as well as some results, to calculate the flow over a generic fighter configuration. A parabolized marching Navier-Stokes code is used to obtain the solution over a wing-canopy body. The flow conditions simulate supersonic cruise with a freestream Mach number of 2.169 and angles of attack of 4 and 10 deg. The body surface is considered to be an adiabatic wall and the flow is assumed to be turbulent for the given Reynolds number.

Chaussee, D. S.↗

Gaussian FLOWERS: Wind-rose-based analytical integration of Gaussian wake model for extremely fast AEP estimation

A major cost in the study of wind farm layout optimization is the repeated evaluation of the annual energy production (AEP). The current approach to estimating AEP requires a large set of flow simulations to be performed that cover each discrete wind speed and direction combination contained within the wind rose, followed by a probability-weighted sum of the power production resulting from each simulation. Even with inexpensive engineering wake models, this numerical integration scheme can lead to high computational costs. In this paper, we derive an analytical formulation for estimating farm AEP across every wind direction, based on a Gaussian wake velocity model, which reduces the number of wind farm simulations to a single function evaluation. As a result, we find that the Gaussian-FLOWERS approach reduces the time for AEP calculations by more than two orders of magnitude with a small trade-off in accuracy when compared to a conventional approach. This massive reduction in computation cost is useful to reduce overall costs in wind farm layout optimization studies.

17 WIND ENERGY↗

Numerical solution of axisymmetric boattail flow fields with plume simulators

Turbulent separating flows over axisymmetric afterbody-boattail configurations with solid sting plume simulators are computed with a time-dependent finite-difference method to solve the compressible Navier-Stokes equations. The Reynolds stress terms are replaced with a two-layer eddy viscosity model including a relaxation formula to model the nonequilibrium effects of the separated flow. The mesh is alined with the boattail body through an analytic transformation which accommodates a wide variety of boattail geometries. Numerical results for a series of boattail geometries over a wide range of Reynolds number (140,000 to 140 million) are presented and, when possible, compared with experimental data or independent numerical results.

Holst, T. L.↗

Flexible AI Models for Grid Resilience

The rapid growth in size and complexity of artificial intelligence (AI) and machine learning (ML) models has led to increased energy demands, posing a threat to the reliability of the existing power grid. This project addresses the challenge of highly intermittent and energy-intensive inference workloads by (1) developing fidelity-adaptive neural networks capable of dynamic response to grid conditions and (2) integrating these networks with power flow simulations to assess their impact on power grid reliability. We will explore both top-down and bottom-up approaches to create hierarchies of submodels that provide a controlled trade-off between power draw and prediction accuracy. The top-down method utilizes NN pruning to reduce a flagship model into progressively smaller, energy-efficient variants. The bottom-up approach employs geometrically principled weight setting strategies to construct depth-efficient models from the ground up. A real-time hardware-in-the-loop (HIL) platform will be developed to simulate a scaled AC power grid, integrating live AI workload power draw and enabling dynamic model switching in response to grid feedback. This work will provide a novel framework for evaluating the impact of flexible AI/ML workloads on grid performance and establish new methodologies for energy-aware computing in data centers. The outcomes will demonstrate that adaptive AI/ML can play a critical role in improving grid stability while advancing NREL's leadership in energy-efficient computing research.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Large-eddy simulations of turbulent wake flows behind helical- and straight-bladed vertical axis wind turbines rotating at low tip speed ratios

Turbulent wake flows behind helical- and straight-bladed vertical axis wind turbines (VAWTs) rotating at low tip speed ratios (TSRs) are studied numerically. The turbulent flows are simulated using the large-eddy simulation (LES) model, and the rotating turbine blades are modeled using the actuator line method. The helical VAWT has identical key parameters as the straight VAWT except for the 135° helical twist of the blades over the 0.3m vertical span. A set of LES runs are performed for two TSRs, 0.6 and 0.4, and the results are reported and analyzed. At these low TSRs, the wake behind the straight-bladed VAWT exhibits two-dimensional dominant flow motions (in the horizontal plane perpendicular to the straight blades) in the near-wake region that cause considerable spanwise expansion of the wake as it extends downstream. In contrast, the helical-bladed VAWT generates highly three-dimensional (3D) wake flow structures and upward/downward mean flow motions within the wake that cause the wake to expand mainly in the vertical direction. Turbulence statistical analyses also show that the 3D wake flow features induced by the helical blades accelerate the wake transition to turbulence and enhance the small-scale turbulent dissipation (as shown by the subgrid-scale turbulent dissipation in the LES), which leads to a more rapid decay of the wake turbulence intensity than that in the straight-bladed VAWT case at the same TSR. Finally, compared with the straight-bladed VAWT, the helical-bladed VAWT also exhibits much smaller temporal variations for the torque and power coefficients during the rotation cycle, which can be beneficial for wind power generation.

17 WIND ENERGY↗

CRADA Number NFE-24-10110 with Qubit Engineering Inc. (CRADA Final Report)

Over the past year, the Qubit Engineering team has pushed the frontiers of power‑grid optimization, working in close collaboration with Oak Ridge National Laboratory (ORNL) and the Tennessee Valley Authority (TVA). Their progress is reflected in three newly submitted conference papers, “Unified Relational GNN Architecture for AC Optimal Power Flow Calculations in Electric Grids,” “Graph‑Based Attention Mechanisms for Solving the AC Optimal Power Flow Problem in Electrical‑Power Networks,” and “Enhanced Power‑Grid Maintenance Planning and Quantum‑Inspired Combinatorial Prospects.” These publications showcase state‑of‑the‑art graph‑neural‑network methods for AC‑OPF and novel quantum‑inspired heuristics for maintenance scheduling. Beyond the academic results, the Qubit team has converted the research into two production‑grade tools built on TVA data: Neuro‑Grid, an AI‑driven power‑flow simulator that provides instant, interactive full‑grid load‑flow visualizations, and Quanta‑Grid, a quantum‑inspired maintenance‑scheduling engine to support logistics optimization for power utilities. Together, these advances demonstrate how Qubit’s partnership with ORNL and TVA is delivering practical, physics‑grounded analytics for next‑generation grid management.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Numerical aerodynamic simulation of the space shuttle ascent environment

After the STS 51-L accident, an extensive review of the Space Shuttle Orbiter's ascent aerodynamic loads uncovered several questionable areas that required further analysis. The insight gained by comparing the Shuttle ascent CFD numerical simulations, obtained by the NASA Ames Space Shuttle Flow Simulation Group, to the current IVBC-3 aerodynamic loads database was instrumental in resolving uncertainties on the Orbiter payload bay doors and fuselage. Initial confidence in the numerical simulations was gained by comparing them with the limited flight data that had been obtained during the Orbiter Flight Test (OFT) program. Current CFD results exist for Mach numbers 0.6, 0.9, 1.05, 1.55, 2.0, and 2.5. Since the pre STS-1 wind tunnel test program (IA-105) often yields considerable differences when compared to STS-5 flight data, the M(sub infinity) = 1.05 transonic case is the most investigated. The IA308 mated-vehicle hot gas plume wind tunnel test, recently completed at AEDC 16T (transonic) and Lewis (hypersonic), is also used to compare with the computation where applicable.

Slotnick, Jeff P.↗

Viscous/Inviscid Interaction On An Oscillating Airfoil

Viscous/inviscid-interaction procedure developed for computation of two-dimensional flow of air about airfoil, angle of attack of which oscillates at high amplitude. Purpose of procedure to study flow effects in dynamic stall. Bulk of flow treated as inviscid, while effects of viscosity treated only on surface of airfoil and in boundary layer. Inviscid part of flow simulated via time-marching, unsteady-subsonic-flow panel method incorporating appropriate distributions of unsteady surface and wake singularities.

Carr, L. W.↗

Data-parallel lower-upper relaxation method for reacting flows

The implicit lower-upper symmetric Gauss-Seidel (LU-SGS) method of Yoon and Jameson is modified for use on massively parallel computers. The method has been implemented on the Thinking Machines CM-5 and the MasPar MP-1 and MP-2, where large percentages of the theoretical peak floating point performance are obtained. It is shown that the new data-parallel LU relaxation method has better convergence properties than the original method for two different inviscid compressible flow simulations. The convergence is also improved for five-species reacting air computations. The performance of the method on various partitions of the CM-5 and on the MasPar computers is discussed. The new method shows promise for the efficient simulation of very large perfect gas and reacting flows.

Candler, Graham V.↗

A Parametric Study of a Plug Nozzle, Using the Liquid Propellant Program (LPP) Code

The Liquid Propellant Program (LPP) computer code is a super-set of the industry standard Two Dimensional Kinetics (TDK) computer code. The TDK code uses a two dimensional method of characteristics solution with fully coupled finite rate kinetics for axially symmetric nozzles. The chemical reactions are modeled with a generalized reaction package that include three dimensional body efficiencies and four reaction rate forms. The code performs optional solutions for frozen or equilibrium flow. TDK evaluates discrete shocks, both attached or induced. The Transonic module models variable mixture ratio profiles from the combustion chamber injector. The Mass Addition Boundary Layer module (MABL) calculates the boundary parameters with the same chemistry options, and includes transpiration or tangential slot injection of gas at the wall. The LPP upgrades include: planar nozzle, scarfed nozzles, plug nozzles, and scramjet nozzle configurations. The code evaluates both upper and lower wall flow simulation, and includes the interaction with the external flow. The MABL module evaluates equilibrium radiation heat transfer for both upper and lower walls. In addition, LPP code models combustion effects due to injector inefficiencies with the Spray Combustion Analysis Program (SCAP) module. The LPP package provides extensive post plotting capabilities for flow visualization. The LPP is sufficiently fast and robust to provide performance predictions for extensive parametric studies and sufficiently accurate to provide flow field and performance solutions for detailed studies.

Dunn, Stuart S.↗

Direct numerical simulations of fluid flow, heat transfer and phase changes

Direct numerical simulations of fluid flow, heat transfer, and phase changes are presented. The simulations are made possible by a recently developed finite difference/front tracking method based on the one-field formulation of the governing equations where a single set of conservation equations is written for all the phases involved. The conservation equations are solved on a fixed rectangular grid, but the phase boundaries are kept sharp by tracking them explicitly by a moving grid of lower dimension. The method is discussed and applications to boiling heat transfer and the solidification of drops colliding with a wall are shown.

Phase Transformations↗

On the Numerical Convergence to Steady State of Hypersonic Flows Over Bodies with Concavities

Two recent numerical studies of hypersonic flows over bodies with concavities revealed problems with convergence to a steady state with an oft-used application of local-time-stepping. Both simulated flows showed a time-like, periodic shedding of vortices in a subsonic domain bounded by supersonic external flow although the simulations, using local-time-stepping, were not time accurate. Simple modifications to the numerical algorithm were implemented to enable implicit, first-order accurate in time simulations. Subsequent time-accurate simulations of the two test problems converged to a steady state. The baseline algorithm and modifications for temporal accuracy are described. The requirement for sub-iterations to achieve convergence is demonstrated. Failure to achieve convergence without time accuracy is conjectured to arise from temporal errors being continuously refocused into a subsonic domain.

Gnoffo, Peter A.↗

Monte Carlo simulation of radiating reentry flows

The Direct Simulation Monte Carlo (DSMC) method is applied to a radiating, hypersonic, axisymmetric flow over a blunt body in the near continuum regime. The ability of the method to predict the flowfield radiation and the radiative heating is investigated for flow over the Project Fire II configuration at 11.36 kilometers per second at an altitude of 76.42 kilometers. Two methods that differ in the manner in which they treat ionization and estimate electronic excitation are employed. The calculated results are presented and compared with both experimental data and solutions where radiation effects were not included. Differences in the results are discussed. Both methods ignore self absorption and, as a result, overpredict measured radiative heating.

Taylor, Jeff C.↗

Computer simulation of hypersonic flow over the Space Shuttle Orbiter

Computer simulations of the flow field around the Space Shuttle Orbiter are described. Results of inviscid calculations are presented for the shock wave pattern and bottom centerline pressure distribution at 30 deg angle of attack. Results of viscous calculations are presented for wall pressure and heat transfer distributions for simple configurations representative of regions where shock wave-boundary layer interactions occur. The computer codes are verified by comparisons with wind-tunnel data and can be applied to flight conditions.

Inouye, M.↗

Numerical simulation of unsteady flow in an axisymmetric shear layer

Numerical simulation of unsteady flow in an axisymmetric subsonic shear layer is accomplished by solving the time-dependent compressible Navier-Stokes equations. The objective of the effort is to investigate by numerical means the influence of various flow parameters on the shear layer behavior. The parameters investigated include the velocity ratio of two streams, total temperature, and nozzle lip thickness. The computations are performed on a CRAY-IS computer using McCormack's explicit finite difference scheme. The computed results generally show qualitative agreement with experimental data.

Scott, J. N.↗

Methodology for simulation of unsteady flows to determine the time-dependent interference between stationary and moving boundaries

A new methodology is developed to simulate unsteady flows about prescribed and aerodynamically determined moving boundary problems. The method couples the fluid dynamics and rigid-body dynamics equations to capture the time-dependent interference between stationary and moving boundaries. The unsteady, compressible, inviscid (Euler) equations are solved on dynamic, unstructured grids by an explicit, finite-volume, upwind method. For efficiency, the grid adaptation is performed within a window around the moving object. The Eulerian equations of the rigid-body dynamics are solved by a Runge-Kutta method in a non-inertial frame of reference. The two-dimensional flow solver is validated by computing the flow past a sinusoidally-pitching airfoil and comparing these results with the experimental data. The overall methodology is used for two two-dimensional examples: the flow past an airfoil which is performing a three-degrees-of-freedom motion in a transonic freestream, and the free-fall of a store after separation from a wing-section. Then the unstructured mesh methodology is extended to three-dimensions to simulate unsteady flow past bodies in relative motion, where the trajectory is determined from the instantaneous aerodynamics. The flow solver and the adaptation scheme in three dimensions are validated by simulating the transonic, unsteady flow around a wing undergoing a forced, periodic, pitching motion, and comparing the results with the experimental data. To validate the trajectory code, the six-degrees-of-freedom motion of a store separating from a wing was computed using the experimentally determined force and moment fields, then comparing with an independently generated trajectory. Finally, the overall methodology was demonstrated by simulating the unsteady flowfield and the trajectory of a store dropped from a wing. The methodology, its computational cost notwithstanding, has proven to be accurate, automated, easy for dynamic gridding, and relatively efficient for the required man-hours.

Singh, Kamakhya Prasad↗