SEARCH · Search NASA
Results for “grid simulation”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
An electrostatic analog for generating cascade grids
Accurate and efficient numerical simulation of flows through turbomachinery blade rows depends on the topology of the computational grids. These grids must reflect the periodic nature of turbomachinery blade row geometries and conform to the blade shapes. Three types of grids can be generated that meet these minimal requirements: through-flow grids, O-type grids, and C-type grids. A procedure which can be used to generate all three types of grids is presented. The resulting grids are orthogonal and can be stretched to capture the essential physics of the flow. A discussion is also presented detailing the extension of the generation procedure to three dimensional geometries.
Grid generation and adaptation for the direct simulation Monte Carlo method
Explore the source record for details and available documents.
Numerical simulation of a hovering rotor using embedded grids
The flow field for a rotor blade in hover was computed by numerically solving the compressible thin-layer Navier-Stokes equations on embedded grids. In this work, three embedded grids were used to discretize the flow field - one for the rotor blade and two to convect the rotor wake. The computations were performed at two hovering test conditions, for a two-bladed rectangular rotor of aspect ratio six. The results compare fairly with experiment and illustrates the use of embedded grids in solving helicopter type flow fields.
An implicit numerical scheme for the simulation of internal viscous flows on unstructured grids
The Navier-Stokes equations are solved numerically for two-dimensional steady viscous laminar flows. The grids are generated based on the method of Delaunay triangulation. A finite-volume approach is used to discretize the conservation law form of the compressible flow equations written in terms of primitive variables. A preconditioning matrix is added to the equations so that low Mach number flows can be solved economically. The equations are time marched using either an implicit Gauss-Seidel iterative procedure or a solver based on a conjugate gradient like method. A four color scheme is employed to vectorize the block Gauss-Seidel relaxation procedure. This increases the memory requirements minimally and decreases the computer time spent solving the resulting system of equations substantially. A factor of 7.6 speed up in the matrix solver is typical for the viscous equations. Numerical results are obtained for inviscid flow over a bump in a channel at subsonic and transonic conditions for validation with structured solvers. Viscous results are computed for developing flow in a channel, a symmetric sudden expansion, periodic tandem cylinders in a cross-flow, and a four-port valve. Comparisons are made with available results obtained by other investigators.
On the metal/ZnO contacts in a sliding-bending piezoelectric nanogenerator
We carry out molecular dynamics (MD) simulations and analyze the piezoelectric and mechanical properties of a ZnO nanowire (NW) under mechanical bending in a sliding-bending piezoelectric nanogenerator. The computational model contains a vertically aligned ZnO NW along the [0001] direction and a Pt (111) metal tip. MD simulations are performed to mimic the interfacial sliding dynamics of an atomic force microscope (AFM) metal tip sliding over a ZnO NW. Ionic charges for Zn and O in the ZnO NW are used to calculate the piezoelectric potential distributions in a grid-point matrix. Simulation results reveal the very detailed distributions of the piezoelectric potential within the ZnO NW, which are different from the continuum field theory predictions. The very similar variations of the lateral friction and contact forces versus the sliding distance of the AFM metal tip over the ZnO NW during repeated forward-backward sliding suggest that the metal/ZnO contact is robust. Moreover, the substantial contact area at the metal/ZnO interface suggests that a Schottky contact should be established at the interface, which needs further investigations.
Quadrature Moments Method for the Simulation of Turbulent Reactive Flows
A sub-filter model for reactive flows, namely the DQMOM model, was formulated for Large Eddy Simulation (LES) using the filtered mass density function. Transport equations required to determine the location and size of the delta-peaks were then formulated for a 2-peak decomposition of the FDF. The DQMOM scheme was implemented in an existing structured-grid LES solver. Simulations of scalar shear layer using an experimental configuration showed that the first and second moments of both reactive and inert scalars are in good agreement with a conventional Lagrangian scheme that evolves the same FDF. Comparisons with LES simulations performed using laminar chemistry assumption for the reactive scalar show that the new method provides vast improvements at minimal computational cost. Currently, the DQMOM model is being implemented for use with the progress variable/mixture fraction model of Pierce. Comparisons with experimental results and LES simulations using a single-environment for the progress-variable are planned. Future studies will aim at understanding the effect of increase in environments on predictions.
Second release of the CoRe database of binary neutron star merger waveforms
Abstract We present the second data release of gravitational waveforms from binary neutron star (BNS) merger simulations performed by the Computational Relativity ( CoRe ) collaboration. The current database consists of 254 different BNS configurations and a total of 590 individual numerical-relativity simulations using various grid resolutions. The released waveform data contain the strain and the Weyl curvature multipoles up to ℓ = m = 4 . They span a significant portion of the mass, mass-ratio, spin and eccentricity parameter space and include targeted configurations to the events GW170817 and GW190425. CoRe simulations are performed with 18 different equations of state, seven of which are finite temperature models, and three of which account for non-hadronic degrees of freedom. About half of the released data are computed with high-order hydrodynamics schemes for tens of orbits to merger; the other half is computed with advanced microphysics. We showcase a standard waveform error analysis and discuss the accuracy of the database in terms of faithfulness. We present ready-to-use fitting formulas for equation of state-insensitive relations at merger (e.g. merger frequency), luminosity peak, and post-merger spectrum.
Battery Performance, Thermal, and Life Modeling for Southern California Edison (Cooperative Research and Development Final Report)
The objective of this project is to develop battery performance models, thermal models and life models for SCE to use in simulation and analysis of grid services involving energy storage. NREL will fit the models to three separate chemistries using test data provided by SCE. The dataset for each chemistry will include electrical, thermal and aging response to different temperatures and cycling conditions, measured under a variety of cell-level and module-level experiments described in SCE test protocol documents. Once tuned to a test dataset, the battery system model software will provide predictions of battery energy and power loss, and thus cycle and calendar life, for any energy storage grid service power profile of interest to SCE. The automated battery life modeling and simulation tool can also be applied to any other grid battery systems as long as necessary input data is available to accelerate battery model development and battery lifetime analysis.
Advanced Inverter Voltage Controls: Simulation and Field Pilot Findings
This report describes work performed by the Hawaiian Electric Companies and the National Renewable Energy Laboratory (NREL) to model and simulate advanced inverter grid-support functions and to validate and expand on those simulations through a field pilot study. This work builds on earlier research, referred to as the Voltage Regulation Operational Strategies study. The objective of both is to investigate functionalities available in most photovoltaic (PV) inverters to modulate active and reactive power autonomously based on local voltage measurements for the purpose of mitigating off-nominal grid voltage conditions. Specifically of interest are volt/volt-ampere reactive (VAR) control and volt/watt control, the effect of those functions on quasi-steady-state feeder voltages, and the impact of the functions on PV energy production. Because volt/VAR and volt/watt control autonomously adjust inverter output based on local conditions without requiring communication with any other devices, they are good candidates to increase PV hosting capacity when it is limited by voltage constraints in a feeder with very large numbers of PV systems.
High fidelity actuator line data from 9 turbine wind farm simulations using ExaWind
This data was generated with the ExaWind code suite (https://github.com/Exawind) to investigate the performance of different Active Wake Mixing turbine control in a wind farm situated in a stable atmospheric boundary layer. All cases correspond to a 3x3 wind farm in a 10km x 10km domain using a total mesh size that varied between 1.6 X 10^9 to 1.85 X 10^9 grid cells. The simulations were run across 1800-2000 GPUs on Frontier. The case description and data generation process is fully documented in Yalla, G. R., Brown, K., Cheung, L., Houck, D., deVelder, N., and Balaji, J. (2025). "Estimating annual energy production of wake mixing control strategies including comparisons to wake steering." Wind Energy Sciences (https://doi.org/10.5194/wes-2025-250).
CRADA Number NFE-19-07846 with Grid Fruit, LLC (CRADA Final Report)
Cooperative Research and Development Agreement (CRADA) NFE-19-07846 between Oak Ridge National Laboratory (ORNL) and Grid Fruit, LLC (Grid Fruit) focused on simulating, monitoring, and adjusting controls of commercial refrigerator and freezer systems to provide load flexibility and demand response services from these machines. Viewing these medium and low-temp coolers as thermal masses, the volumes of chilled air in these coolers make them time-shiftable loads ripe for deployment at optimal times. Grid Fruit is a startup company developing methods to schedule the chilling cycles (e.g., in “build load” and “shed load” grid events) to provide energy efficiency, peak shifting, and other benefits both to the grid and the broader environment. These benefits can be magnified by coordinating chillers with HVAC, lights, computers, and other loads. Grid Fruit completed simulations and then selected and ordered controls hardware to prove the technology were successful, but the diversity of chiller models in the field without digital controls, as well as limited market demand at present, make further exploration necessary before commercialization.
Momentum Transport: 2D and 3D Cloud Resolving Model Simulations
The major objective of this study is to investigate the momentum budgets associated with several convective systems that developed during the TOGA COARE IOP (west Pacific warm pool region) and GATE (east Atlantic region). The tool for this study is the improved Goddard Cumulas Ensemble (GCE) model which includes a 3-class ice-phase microphysical scheme, explicit cloud radiative interactive processes and air-sea interactive surface processes. The model domain contains 256 x 256 grid points (with 2 km resolution) in the horizontal and 38 grid points (to a depth of 22 km) in the vertical. The 2D domain has 1024 grid points. The simulations were performed over a 7-day time period (December 19-26, 1992, for TOGA COARE and September 1-7, 1994 for GATE). Cyclic literal boundary conditions are required for this type of long-term integration. Two well organized squall systems (TOGA, COARE February 22, 1993, and GATE September 12, 1994) were also simulated using the 3D GCE model. Only 9 h simulations were required to cover the life time of the squall systems. the lateral boundary conditions were open for these two squall systems simulations. the following will be examined: (1) the momentum budgets in the convective and stratiform regions, (2) the relationship between momentum transport and cloud organization (i.e., well organized squall lines versus less organized convective), (3) the differences and similarities in momentum transport between 2D and 3D simulated convective systems, and (4) the differences and similarities in momentum budgets between cloud systems simulated with open and cyclic lateral boundary conditions. Preliminary results indicate that there are only small differences between 2D and 3D simulated momentum budgets. Major differences occur, however, between momentum budgets associated with squall systems simulated using different lateral boundary conditions.
Influence of Antarctic and Greenland Continental Shelf Circulation on High‐Latitude Oceans in E3SM
The science objectives of this project are to simulate and understand the impacts of both deep-basin warm-water intrusions and land-ice melt on the continental shelf circulations and sea-ice distributions around the margins of Greenland and Antarctica. As well, the role of subsurface ocean heat from the Atlantic on declining sea-ice cover in the Arctic is explored. Mesoscale processes and fine bathymetry are implicated in cross-shelf property transports around both Greenland and Antarctica. Therefore, we configured and ran an atmospheric reanalysis-forced global ocean/sea-ice simulation on a grid that reduces from 8 km at the Equator to 2 km at the poles (UH8to2) with 60 vertical levels. It was produced using the Energy Exascale Earth System Model ‘‘HiLAT’’ code (E3SMv0-HiLAT) that uses the Parallel Ocean Program (POP) and CICE5 as its ocean and sea-ice components, respectively. Two main UH8to2 simulations were carried out: one for 1975-2009 and the other for July 2016-2020 after it was initialized from a 1/25° data-assimilative ocean/sea-ice prediction system ocean/sea-ice state. The UH8to2 is not coupled to an active land-ice model. Rather, land-ice melt is represented by observationally informed freshwater fluxes (FWFs). Short (multi-year) UH8to2 simulations were conducted to understand sensitivities when Greenland ice sheet (GrIS) melt is released only at the ocean surface or when it is distributed over the upper water column in accordance with fjord melt plume behavior; these cases were compared with a no GrIS melt case. West Greenland continental shelf currents were fastest in the vertical distribution case and an increase in baroclinic conversion at the shelf break associated with increased eddy kinetic energy was found relative to the surface release case. Further, salinity is lower and meltwater volume greater in the eastern Labrador Sea in the vertical distribution case. For the Arctic, the veracity of the UH8to2 was evaluated for 2017-2020 using available observations. Simulated seasonal sea-ice thickness and concentration are realistic, but the ice is unrealistically thin in the central and eastern Arctic in the fall. Comparisons of vertical sections of ocean temperature, salinity, and buoyancy collected from Ice-Tethered Profilers (ITPs) in the eastern Arctic in the fall and winter of 2019/2020 and co-located/concurrent UH8to2 fields show the stratification over the top 100 m of the water column is too low in the model, the simulated mixed layer too deep, and the simulated subsurface Atlantic Water (AW) too warm; these biases may contribute to the sea-ice biases. A model intercomparison study using the UH8to2 and a forced 1/25° regional Arctic ocean/sea-ice (uses the HYbrid Coordinate Ocean Model and CICE5) simulation further investigates the relationship between AW and sea-ice in the eastern Arctic. The models show a mesoscale-rich pulse of Atlantic Water extending into the eastern basin that reaches maximum intensity in late winter of 2018, after which it decreases in strength. Concurrent and co-located sea-ice melt or the inhibition of sea-ice growth is seen and is attributed to halocline mesoscale eddies doming into the mixed layer with convection bringing this heat into the vicinity of the sea-ice.
Changes in Annual Extremes of Daily Temperature and Precipitation in CMIP6 Models
This study presents an analysis of daily temperature and precipitation extremes with return periods ranging from 2 to 50 years in phase 6 of the Coupled Model Intercomparison Project (CMIP6) multimodel ensemble of simulations. Judged by similarity with reanalyses, the new-generation models simulate the present-day temperature and precipitation extremes reasonably well. In line with previous CMIP simulations, the new simulations continue to project a large-scale picture of more frequent and more intense hot temperature extremes and precipitation extremes and vanishing cold extremes under continued global warming. Changes in temperature extremes outpace changes in global annual mean surface air temperature (GSAT) over most landmasses, while changes in precipitation extremes follow changes in GSAT globally at roughly the Clausius–Clapeyron rate of ~7% °C –1 . Changes in temperature and precipitation extremes normalized with respect to GSAT do not depend strongly on the choice of forcing scenario or model climate sensitivity, and do not vary strongly over time, but with notable regional variations. Over the majority of land regions, the projected intensity increases and relative frequency increases tend to be larger for more extreme hot temperature and precipitation events than for weaker events. To obtain robust estimates of these changes at local scales, large initial-condition ensemble simulations are needed. Appropriate spatial pooling of data from neighboring grid cells within individual simulations can, to some extent, reduce the needed ensemble size.
Grid Fin Stabilization of the Orion Launch Abort Vehicle
Wind tunnel tests were conducted by Nielsen Engineering & Research (NEAR) and Rose Engineering & Research (REAR) in conjunction with the NASA Engineering & Safety Center (NESC) on a 6%-scale model of the Orion launch abort vehicle (LAV) configured with four grid fins mounted near the base of the vehicle. The objectives of these tests were to 1) quantify LAV stability augmentation provided by the grid fins from subsonic through supersonic Mach numbers, 2) assess the benefits of swept grid fins versus unswept grid fins on the LAV, 3) determine the effects of the LAV abort motors on grid fin aerodynamics, and 4) generate an aerodynamic database for use in the future application of grid fins to small length-to-diameter ratio vehicles similar to the LAV. The tests were conducted in NASA Ames Research Center's 11x11-foot transonic wind tunnel from Mach 0.5 through Mach 1.3 and in their 9x7-foot supersonic wind tunnel from Mach 1.6 through Mach 2.5. Force- and moment-coefficient data were collected for the complete vehicle and for each individual grid fin as a function of angle of attack and sideslip angle. Tests were conducted with both swept and unswept grid fins with the simulated abort motors (cold jets) off and on. The swept grid fins were designed with a 22.5deg aft sweep angle for both the frame and the internal lattice so that the frontal projection of the swept fins was the same as for the unswept fins. Data from these tests indicate that both unswept and swept grid fins provide significant improvements in pitch stability as compared to the baseline vehicle over the Mach number range investigated. The swept fins typically provide improved stability as compared to the unswept fins, but the performance gap diminished as Mach number was increased. The aerodynamic performance of the fins was not observed to degrade when the abort motors were turned on. Results from these tests indicate that grid fins can be a robust solution for stabilizing the Orion LAV over a wide range of operating conditions.
Validation of Supersonic Film Cooling Modeling for Liquid Rocket Engine Applications
Topics include: upper stage engine key requirements and design drivers; Calspan "stage 1" results, He slot injection into hypersonic flow (air); test articles for shock generator diagram, slot injector details, and instrumentation positions; test conditions; modeling approach; 2-d grid used for film cooling simulations of test article; heat flux profiles from 2-d flat plate simulations (run #4); heat flux profiles from 2-d backward facing step simulations (run #43); isometric sketch of single coolant nozzle, and x-z grid of half-nozzle domain; comparison of 2-d and 3-d simulations of coolant nozzles (run #45); flowfield properties along coolant nozzle centerline (run #45); comparison of 3-d CFD nozzle flow calculations with experimental data; nozzle exit plane reduced to linear profile for use in 2-d film-cooling simulations (run #45); synthetic Schlieren image of coolant injection region (run #45); axial velocity profiles from 2-d film-cooling simulation (run #45); coolant mass fraction profiles from 2-d film-cooling simulation (run #45); heat flux profiles from 2-d film cooling simulations (run #45); heat flux profiles from 2-d film cooling simulations (runs #47, #45, and #47); 3-d grid used for film cooling simulations of test article; heat flux contours from 3-d film-cooling simulation (run #45); and heat flux profiles from 3-d and 2-d film cooling simulations (runs #44, #46, and #47).
Revival of the fittest: exploding core-collapse supernovae from 12 to 25 M⊙
Here, we present results of 2D axisymmetric core-collapse supernova simulations, employing the FORNAX code, of nine progenitor models spanning 12 to 25 M⊙. Four of the models explode with inelastic scattering off electrons and neutrons as well as the many-body correction to neutrino-nucleon scattering opacities. We show that these four models feature sharp Si–O interfaces in their density profiles, and that the corresponding dip in density reduces the accretion rate around the stalled shock and prompts explosion. The non-exploding models lack such a steep feature, highlighting the Si–O interface as one key to explosion. Furthermore, we show that all of the non-exploding models can be nudged to explosion with modest changes to macrophysical inputs, including moderate rotation and perturbations to infall velocities, as well as to microphysical inputs, including reasonable changes to neutrino-nucleon interaction rates, suggesting that all the models are perhaps close to criticality. Exploding models have energies of a few × 10 50 erg at the end of our simulation, and are rising, emphasizing the need to continue these simulations over larger grids and for longer times to reproduce the energies seen in nature. Morphology of the explosion contributes to the explosion energy, with more isotropic ejecta producing larger explosion energies. We do not find evidence for the Lepton-number Emission Self-sustained Asymmetry. Finally, we look at proto-neutron star (PNS) properties and explore the role of dimension in our simulations. We find that convection in the PNS produces larger PNS radii as well as greater ‘ν μ ’ luminosities in 2D compared to 1D.