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.

SourceResearch connectionFrom yearThrough year

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

At least 1,171 records · Page 65

All sources NASA NTRS DOE OSTI

DOE OSTI2023

RANS Simulation of Variable Density Turbulent Round Jets with Coflow using xRAGE Hydrodynamic Code and BHR Turbulence Models

This work for the fiscal year 2023 (FY23) is a continuation of previous efforts to evaluate the BHR turbulence models for their ability to accurately simulate variable density turbulent round jets with coflow. As before, RANS simulations are carried out using the xRAGE hydrodynamic code. The following are some of the previous findings. Israel showed that i) the symmetry boundary conditions for the BHR models in axisymmetric simulations were in error, ii) three grids of different resolutions did not lead to converging solutions, and iii) BHR 3.1 simulation exhibited instabilities and did not reach a steady state. Saenz and Rauenzahn derived and implemented into xRAGE the correct BHR boundary conditions at the symmetry axis. Cline conducted sensitivity studies with various parameters including the BHR models (versions 2, 3.1 and 4), gravity, material pressure, specific heat, initial turbulent kinetic energy and initial turbulent length scale, and found that the largest factor impacting on simulation results was the BHR model version, followed by the initial turbulent length scale. In addition, freeze boundary conditions at the exit and the side wall of the computational domain were used to remove anomalous flow behavior. Cline adjusted the inlet jet velocity, initial turbulent kinetic energy and initial turbulent length scale to obtain the best reasonable match with the experimental data of Charonko and Prestridge. It was found that the BHR 2 and 3.1 models performed in a similar manner, but the BHR 2 model produced much lower levels of density-specific-volume covariance and turbulent kinetic energy. The main focus for the FY23 is to study the effects of computational parameters related to the boundary conditions, mesh, domain size and timestep size. The reasoning behind this is that, unless simulation results can be shown to be reasonably independent from the aforementioned computational parameters, it would be difficult to attribute any discrepancies between simulation and experimental results to turbulence models. This important aspect has largely been overlooked in the previous years, and therefore will be studied comprehensively here. Additionally, effects of varying the initial turbulent length scale will be examined because it was previously identified as a major factor affecting the flow fields.

42 ENGINEERING↗

DOE OSTI2020

Application of Intelligent Load Control to Manage Building Loads to Support Rapid Growth of Distributed Renewable Generation

Electricity utilities are faced with the mounting challenge of providing a stable supply of power to meet the growing demand while also integrating rapid growth in distributed variable renewable generation. Traditional means of balancing short- and long-term supply and demand imbalance will be expensive. Alternative approaches of using flexible loads in buildings are needed to mitigate the imbalance at a lower cost. This paper shows how the intelligent load control (ILC) process can be used to manage loads in buildings by dynamically prioritizing loads for curtailment using both quantitative and qualitative criteria. The ILC process can be deployed on low-cost computing platforms without the need for any additional sensing. ILC was first validated in a simulation environment to provide two grid service use cases: (1) managing monthly peak electricity demand and (2) managing buildings’ electricity consumption during a capacity bidding event. After being successfully tested in a simulation environment, ILC was deployed on real buildings to manage electricity consumption to provide two different use cases under different outdoor operating conditions. Both the simulation tests and the real building experiments were deployed using VOLTTRON™, a distributed sensing and control platform. The results from the tests and experiments showed that ILC was able to manage the controllable loads (heat pumps) in the building to maintain the electricity consumption at the desired level without a significant impact on occupant comfort. Overall, the results demonstrate that the ILC allows coordination of the controllable loads and provides a more intelligent means of load management than the traditional duty-cycling approach.

Kim, Woohyun↗

DOE OSTI2024

Multi-State Transportation Electrification Impact Study: Preparing the Grid for Light-, Medium-, and Heavy-Duty Electric Vehicles

Recent U.S. Environmental Protection Agency (EPA) notices of proposed rulemakings for GHG emissions standards for light-, medium-, and heavy-duty on-road vehicles would accelerate ongoing advancements already happening in the industry because of private investment, consumer demand, state-level policies, and federal incentives. As the EPA finalizes these regulations, questions persist regarding the cost of the requisite charging infrastructure and associated upgrades to the nation's electric grid. With support from the U.S. Department of Energy, U.S. Joint Office of Energy and Transportation, and the EPA, a multidisciplinary team conducted a Multi-State Transportation Electrification Impact Study that quantitatively assesses the incremental investment necessary to enable the levels of vehicle electrification expected to be induced by pending EPA regulations and to estimate the potential value of deferred investments in electric distribution infrastructure stemming from proactive vehicle-grid integration planning and deployment. This study finds the simulated incremental capital cost of charging infrastructure (including grid upgrades) to be at least 2.5 times smaller than the lifetime net benefits of vehicle electrification (including fuel savings but excluding the value of avoided emissions). Additionally, the incremental distribution grid upgrade cost of the EPA Action-Unmanaged scenario was found to be approximately 3% of existing utility distribution system investments (on an annual basis). Finally, the potential for managed charging to defer distribution grid upgrades was found to be significant with costs found to decrease from $2.3 billion to an incremental cost of $1 billion across five states in the Action-Managed scenario (relative to the No Action-Unmanaged scenario).

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2019

Generalized, Energy-conserving Numerical Simulations of Particles in General Relativity. II. Test Particles in Electromagnetic Fields and GRMHD

Observations of compact objects, in the form of radiation spectra, gravitational waves by LIGO/Virgo, and direct imaging with the Event Horizon Telescope, are currently the main information sources on plasma physics in extreme gravity. Modeling such physical phenomena Requires numerical methods that allow for the simulation of microscopic plasma dynamics in the presence of both strong gravity and electromagnetic fields. In Bacchini et al. (2018), we presented a detailed study of numerical techniques for the integration of free geodesic motion. In this paper, we extend the study by introducing electromagnetic forces in the simulation of charged particles in curved spacetimes. We extend the Hamiltonian energy-conserving method presented in Bacchini et al. to include the Lorentz force, and we test its performance compared to that of standard explicit Runge–Kutta and implicit midpoint rule schemes against analytic solutions. Then, we show the application of the numerical schemes to the integration of test particle trajectories in general relativistic magnetohydrodynamic (GRMHD) simulations by modifying the algorithms to handle grid-based electromagnetic fields. We test this approach by simulating ensembles of charged particles in a static GRMHD configuration obtained with the black hole accretion code (BHAC).

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Transportation Electrification Impact Study (TEIS)

ADVANCED PROPULSION SYSTEMS,POWER TRANSMISSION AND↗

DOE OSTI2017

A low-rank control variate for multilevel Monte Carlo simulation of high-dimensional uncertain systems

Multilevel Monte Carlo (MLMC) is a recently proposed variation of Monte Carlo (MC) simulation that achieves variance reduction by simulating the governing equations on a series of spatial (or temporal) grids with increasing resolution. Instead of directly employing the fine grid solutions, MLMC estimates the expectation of the quantity of interest from the coarsest grid solutions as well as differences between each two consecutive grid solutions. When the differences corresponding to finer grids become smaller, hence less variable, fewer MC realizations of finer grid solutions are needed to compute the difference expectations, thus leading to a reduction in the overall work. Here this paper presents an extension of MLMC, referred to as multilevel control variates (MLCV), where a low-rank approximation to the solution on each grid, obtained primarily based on coarser grid solutions, is used as a control variate for estimating the expectations involved in MLMC. Cost estimates as well as numerical examples are presented to demonstrate the advantage of this new MLCV approach over the standard MLMC when the solution of interest admits a low-rank approximation and the cost of simulating finer grids grows fast.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

Generating one-column grids with fractal flow dimension

The grid generation capability built into the numerical simulator TOUGH for multi-phase fluid and heat flow through geologic media can create one-column grids with linear or radial geometry, corresponding to one-dimensional or two-dimensional radial flow, respectively. The integral-finite-difference-method that TOUGH employs for spatial discretization makes it very simple to generalize the grid-generation algorithm from integer to non-integer (fractal) flow dimension. Here the grid-generation algorithm is generalized to create one-column grids with fractal flow dimension ranging from less than 1 to 3. The fractal grid generation method is verified by comparing numerical simulation results to an analytical solution for a generalized Theis solution for integer and non-integer flow dimensions between 0.4 and 3. It is then applied to examine gas production decline curves from hydraulically fractured shale that is modeled as a fractal-dimensioned fracture network with flow dimensions between 0.25 and 3. Grids with fractal flow dimension are useful for representing flow through fracture networks or highly heterogeneous geologic media with fractal geometry, and may be particularly useful for inverse methods.

58 GEOSCIENCES↗

NASA NTRS2012

The Influence of Pickup Protons, from Interstellar Neutral Hydrogen, on the Propagation of Interplanetary Shocks from the Halloween 2003 Solar Events to ACE and Ulysses: A 3-D MHD Modeling Study

We describe our 3-D, time ]dependent, MHD solar wind model that we recently modified to include the physics of pickup protons from interstellar neutral hydrogen. The model has a time-dependent lower boundary condition, at 0.1 AU, that is driven by source surface map files through an empirical interface module. We describe the empirical interface and its parameter tuning to maximize model agreement with background (quiet) solar wind observations at ACE. We then give results of a simulation study of the famous Halloween 2003 series of solar events. We began with shock inputs from the Fearless Forecast real ]time shock arrival prediction study, and then we iteratively adjusted input shock speeds to obtain agreement between observed and simulated shock arrival times at ACE. We then extended the model grid to 5.5 AU and compared those simulation results with Ulysses observations at 5.2 AU. Next we undertook the more difficult tuning of shock speeds and locations to get matching shock arrival times at both ACE and Ulysses. Then we ran this last case again with neutral hydrogen density set to zero, to identify the effect of pickup ions. We show that the speed of interplanetary shocks propagating from the Sun to Ulysses is reduced by the effects of pickup protons. We plan to make further improvements to the model as we continue our benchmarking process to 10 AU, comparing our results with Cassini observations, and eventually on to 100 AU, comparing our results with Voyager 1 and 2 observations.

Detman, T. R.↗

DOE OSTI2021

A Decision Support Framework for Grid-Aware Electric Bus Charge Scheduling

While there are many advantages to electric public transit vehicles, they also pose new challenges for fleet operators. One key challenge is defining a charge scheduling policy that minimizes operating costs and power grid disruptions while maintaining schedule adherence. An uncoordinated policy could result in some buses running out of charge before completing their trip, while a grid agnostic policy might incur higher energy costs or cause an adverse impact on the grid’s distribution system. We present a grid aware decision-theoretic framework for electric bus charge scheduling that accounts for energy price, grid load, and bus status. The framework connects environmental models for traffic (SUMO) and the electric grid (GridLAB-D) to a planner, which evaluates bus charging decisions with regard to their long-term effect on grid reliability and cost. The framework has been evaluated on a simulation of Richland, WA’s bus and grid network, and compared to greedy methods.

EV, machine leaning, EV-Grid, Transportation, Grid↗

DOE OSTI2025

Dominant balance-based adaptive mesh refinement for incompressible fluid flows

This work introduces a novel adaptive mesh refinement (AMR) method that utilizes dominant balance analysis (DBA) for efficient and accurate grid adaptation in computational fluid dynamics (CFD) simulations. The proposed method leverages a Gaussian mixture model (GMM) to classify grid cells into active and passive regions based on the dominant physical interactions within the equation space. By modeling truncation error probabilistically from discretized terms, the method identifies regions of high interaction where numerical accuracy is most sensitive to resolution. Unlike traditional AMR strategies, this approach does not rely on heuristic-based sensors or user-defined thresholds, providing a fully automated and problem-independent framework for AMR. Applied to the incompressible Navier-Stokes equations for steady and unsteady flow past a cylinder, the DBA-based AMR method achieves comparable accuracy to high-resolution grids while reducing computational costs by up to 70 %. The validation highlights the method’s effectiveness in capturing complex flow features while minimizing grid cells, directing computational resources toward regions with the most critical dynamics. This modular and scalable strategy is adaptable to a wide range of applications, presenting a promising tool for efficient high-fidelity simulations in CFD and other multiphysics domains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

High-resolution Navier-Stokes simulations of Richtmyer-Meshkov instability with reshock

The interaction of a Mach-1.45 shock wave with a perturbed planar interface between sulfur hexafluoride and air is studied through high-resolution two-dimensional (2D) and three-dimensional (3D) shock-capturing adaptive mesh refinement simulations of multispecies Navier-Stokes equations. The sensitivities of time-dependent statistics on grid resolution for 2D and 3D simulations are evaluated to ensure the accuracy of the results. Here, the numerical results are used to examine the differences between the development of 2D and 3D Richtmyer-Meshkov instability during two different stages: (1) initial growth of hydrodynamic instability from multimode perturbations after first shock and (2) transition to chaotic or turbulent state after reshock. The effects of the Reynolds number on the mixing in 3D simulations are also studied through varying the transport coefficients.

42 ENGINEERING↗

DOE OSTI2026

Data and scripts associated with a manuscript analyzing ELM-FATES parameter sensitivity under pre-fire and postfire scenarios using machine learning

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript “Fire Severity-Dependent Shifts in Vegetation Parameter Sensitivity: A Pre- and Post-Fire Analysis Using ELM-FATES and Explainable AI” submitted to Journal of Advances in Modeling Earth Systems (Zahura et al. 2026). The study examines vegetation physiological parameters controlling pre-fire and post-fire vegetation dynamics. To support this analysis, 73 vegetation parameters in Functionally Assembled Terrestrial Ecosystem Simulator (FATES) (Fisher et al., 2018) , which is coupled with E3SM (Energy Exascale Earth System Model) land model (ELM, ELM-FATES), were perturbed using a Sobol sequence to generate 1,024 ensemble members for two plant functional types: needleleaf evergreen extratropical trees (NEET) and C3 grass. Simulations were conducted for the pre-fire period (2016) and post-fire period (2018–2023). Burn severity was represented by modifying the Nesterov index in FATES to 75,000, 150,000, and 300,000 for low, moderate, and high severity, respectively. A no-fire scenario was also included. Simulations were performed for 16 grid cells in the American River Watershed across different burn severities and plant functional types. XGBoost (eXtreme Gradient Boosting) models were trained using the parameter ensembles and ELM-FATES-simulated outputs, including leaf area index (LAI), gross primary productivity (GPP), aboveground biomass, vegetation evaporation, transpiration, and soil evaporation. Models were trained separately for each year and burn severity, followed by SHAP (SHapley Additive exPlanations) analysis to identify changes in dominant parameters after fire disturbance. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. The data package contains the ELM-FATES simulation data. The scripts and data related to the analysis will be added later. The inputs and outputs from ELM-FATES are inside the “FATES” folder. “FATES_domain_surface” contains the domain and surface netcdfs that were used to run ELM-FATES in the study area. “FATES_parameters” contains the 1024 ensembles that were generated using Sobol sequence. “FATES_outputs” folder contains ELM-FATES simulated variables. All files are .csv and .nc (NetCDF).

Aboveground biomass↗

DOE OSTI2026

Reliable Protection for an Inverter-Based Resources Dominant Grid: Technology Development and Field Demonstration

The project aimed to address the challenges posed by the rapid growth of inverter-based resources (IBR) such as solar, wind, and battery storage, which have fundamentally changed fault behavior, system strength, and protection performance in bulk power systems. Traditional protection schemes, designed for synchronous generator-dominated grids, are inadequate under high IBR penetration. Overall, the project materially advances protection modeling and simulation capabilities needed to maintain reliable grid protection under high IBR penetration. The results build utility confidence in operating power systems safely and reliably across a wide range of generation mixes, supporting grid modernization goals.

24 POWER TRANSMISSION AND DISTRIBUTION↗

NASA NTRS2016

Dynamic Mesh CFD Simulations of Orion Parachute Pendulum Motion During Atmospheric Entry

This paper demonstrates the usage of computational fluid dynamics to study the effects of pendulum motion dynamics of the NASAs Orion Multi-Purpose Crew Vehicle parachute system on the stability of the vehicles atmospheric entry and decent. Significant computational fluid dynamics testing has already been performed at NASAs Johnson Space Center, but this study sought to investigate the effect of bulk motion of the parachute, such as pitching, on the induced aerodynamic forces. Simulations were performed with a moving grid geometry oscillating according to the parameters observed in flight tests. As with the previous simulations, OVERFLOW computational fluid dynamics tool is used with the assumption of rigid, non-permeable geometry. Comparison to parachute wind tunnel tests is included for a preliminary validation of the dynamic mesh model. Results show qualitative differences in the flow fields of the static and dynamic simulations and quantitative differences in the induced aerodynamic forces, suggesting that dynamic mesh modeling of the parachute pendulum motion may uncover additional dynamic effects.

Halstrom, Logan D.↗

NASA NTRSDate not supplied

Automated Unstructured Grid Adaptation on a Strut Fuel Injector at Hypervelocity Flow Conditions

Computational fluid dynamics (CFD) analysis is presented with the use of an automated unstructured grid adaptation tool on a strut fuel injector at hypervelocity flow conditions. The analysis was carried out with the VULCAN-CFD solver using Reynolds-averaged simulations (RAS). The hypervelocity flow conditions match the high Mach number flow of the experiments conducted as part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center (LaRC). The current work utilizes an automated grid adaptation tool recently implemented into VULCAN-CFD, and explores this tool’s ability to solve high-speed mixing problems. Simulation results obtained using the unstructured adaptive grid approach are compared to those on a user generated structured grid. These results are evaluated by analyzing how efficiently comparable fidelity results are obtained from both adapted and structured simulations. In addition, two adaptation strategies were used to explore the impact on the final solution. In the current work, the unstructured grid adaptation tool automatically generates unstructured grids and performs adaptation of the grid based on a Hessian error estimate of a specified flow field parameter. Multiple adaptations were executed using each run strategy with the one-dimensional values of the mixing efficiency used to determine grid convergence and for comparison with the structured grid simulation results. It was found that the unstructured adaptive grid simulations were able to produce results that matched closely with those on structured grids using far fewer grid cells, and thus, requiring far less computational time to reach the solution. It was also discovered that the adaptation run strategy influenced the total number of grid cells and the efficiency with which a final grid-adapted solution was reached. Overall, the investigation demonstrated that the automated unstructured grid adaptation tool implemented in VULCAN-CFD is capable of accurately and efficiently solving complex highspeed mixing problems using only a fraction of the grid cells required to obtain comparable results using a user-generated structured grid.

grid adaptation↗

NASA NTRS1995

Computations of Viscous Flows in Complex Geometries Using Multiblock Grid Systems

Generating high quality, structured, continuous, body-fitted grid systems (multiblock grid systems) for complicated geometries has long been a most labor-intensive and frustrating part of simulating flows in complicated geometries. Recently, new methodologies and software have emerged that greatly reduce the human effort required to generate high quality multiblock grid systems for complicated geometries. These methods and software require minimal input form the user-typically, only information about the topology of the block structure and number of grid points. This paper demonstrates the use of the new breed of multiblock grid systems in simulations of internal flows in complicated geometries. The geometry used in this study is a duct with a sudden expansion, a partition, and an array of cylindrical pins. This geometry has many of the features typical of internal coolant passages in turbine blades. The grid system used in this study was generated using a commercially available grid generator. The simulations were done using a recently developed flow solver, TRAF3D.MB, that was specially designed to use multiblock grid systems.

Steinthorsson, Erlendur↗

DOE OSTI2016

An approach for drag correction based on the local heterogeneity for gas-solid flows

The drag models typically used for gas-solids interaction are mainly developed based on homogeneous systems of flow passing fixed particle assembly. It has been shown that the heterogeneous structures, i.e., clusters and bubbles in fluidized beds, need to be resolved to account for their effect in the numerical simulations. Since the heterogeneity is essentially captured through the local concentration gradient in the computational cells, this study proposes a simple approach to account for the non-uniformity of solids spatial distribution inside a computational cell and its effect on the interaction between gas and solid phases. Finally, to validate this approach, the predicted drag coefficient has been compared to the results from direct numerical simulations. In addition, the need to account for this type of heterogeneity is discussed for a periodic riser flow simulation with highly resolved numerical grids and the impact of the proposed correction for drag is demonstrated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Application Note: Power Grid Modeling With Xyce™

This application note describes how to model steady-state power flows and transient events in electric power grids with the SPICE-compatible Xyce™ Parallel Electronic Simulator developed at Sandia National Labs. This application notes provides a brief tutorial on the basic devices (branches, bus shunts, transformers and generators) found in power grids. The focus is on the features supported and assumptions made by the Xyce models for power grid elements. It then provides a detailed explanation, including working Xyce netlists, for simulating some simple power grid examples such as the IEEE 14-bus test case.

24 POWER TRANSMISSION AND DISTRIBUTION↗