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,693 records · Page 94

All sources NASA NTRS DOE OSTI

NASA NTRS2014

A Walsh Function Module Users' Manual

The solution of partial differential equations (PDEs) with Walsh functions offers new opportunities to simulate many challenging problems in mathematical physics. The approach was developed to better simulate hypersonic flows with shocks on unstructured grids. It is unique in that integrals and derivatives are computed using simple matrix multiplication of series representations of functions without the need for divided differences. The product of any two Walsh functions is another Walsh function - a feature that radically changes an algorithm for solving PDEs. A FORTRAN module for supporting Walsh function simulations is documented. A FORTRAN code is also documented with options for solving time-dependent problems: an advection equation, a Burgers equation, and a Riemann problem. The sample problems demonstrate the usage of the Walsh function module including such features as operator overloading, Fast Walsh Transforms in multi-dimensions, and a Fast Walsh reciprocal.

Gnoffo, Peter A.↗

DOE OSTI2012

CoCoMANS Ocean App [Slides]

Thermohaline circulation is that part of ocean circulation which is driven by density gradients. It is widely believed that changes in the ocean thermohaline circulation on the 10 to 100 year time scale can instigate rapid climate change. Scientists world wide are currently using simulation tools based on semi-implicit and split methods that are first-order in time. From a numerical error perspective, the quality of such simulations is suspect. As a result of numerical stability issues resulting from semi-implicit and split methods they are forced to use low resolution grids; ie high resolution grids will require prohibitive numbers of timesteps. Having the ability to perform such simulation with reduced splitting errors, using high resolution grids, and using second-order accurate time integration methods will be a quantum leap forward in predictive ocean simulation.

36 MATERIALS SCIENCE↗

DOE OSTI2015

Large-eddy and unsteady RANS simulations of a shock-accelerated heavy gas cylinder

Two-dimensional numerical simulations of the Richtmyer–Meshkov unstable “shock-jet” problem are conducted using both large-eddy simulation (LES) and unsteady Reynolds-averaged Navier–Stokes (URANS) approaches in an arbitrary Lagrangian–Eulerian hydrodynamics code. Turbulence statistics are extracted from LES by running an ensemble of simulations with multimode perturbations to the initial conditions. Detailed grid convergence studies are conducted, and LES results are found to agree well with both experiment and high-order simulations conducted by Shankar et al. (Phys Fluids 23, 024102, 2011). URANS results using a k–L approach are found to be highly sensitive to initialization of the turbulence lengthscale L and to the time at which L becomes resolved on the computational mesh. As a result, it is observed that a gradient diffusion closure for turbulent species flux is a poor approximation at early times, and a new closure based on the mass-flux velocity is proposed for low-Reynolds-number mixing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2012

Radiation Coupling with the FUN3D Unstructured-Grid CFD Code

The HARA radiation code is fully-coupled to the FUN3D unstructured-grid CFD code for the purpose of simulating high-energy hypersonic flows. The radiation energy source terms and surface heat transfer, under the tangent slab approximation, are included within the fluid dynamic ow solver. The Fire II flight test, at the Mach-31 1643-second trajectory point, is used as a demonstration case. Comparisons are made with an existing structured-grid capability, the LAURA/HARA coupling. The radiative surface heat transfer rates from the present approach match the benchmark values within 6%. Although radiation coupling is the focus of the present work, convective surface heat transfer rates are also reported, and are seen to vary depending upon the choice of mesh connectivity and FUN3D ux reconstruction algorithm. On a tetrahedral-element mesh the convective heating matches the benchmark at the stagnation point, but under-predicts by 15% on the Fire II shoulder. Conversely, on a mixed-element mesh the convective heating over-predicts at the stagnation point by 20%, but matches the benchmark away from the stagnation region.

Wood, William A.↗

DOE OSTI2016

Investigation of boundary-layer wind predictions during nocturnal low-level jet events using the Weather Research and Forecasting model

Simulations of two periods featuring three consecutive low level jet (LLJ) events in the US Upper Great Plains during the autumn of 2011 were conducted to explore the impacts of various setup configurations and physical process models on simulated flow parameters within the lowest 200 m above the surface, using the Weather Research and Forecasting (WRF) model. Sensitivities of simulated flow parameters to the horizontal and vertical grid spacing, planetary boundary layer (PBL) and land surface model (LSM) physics options, were assessed. Data from a Light Detection and Ranging (lidar) system, deployed to the Weather Forecast Improvement Project (WFIP; Finley et al. 2013) were used to evaluate the accuracy of simulated wind speed and direction at 80 m above the surface, as well as their vertical distributions between 120 and 40 m, covering the typical span of contemporary tall wind turbines. All of the simulations qualitatively captured the overall diurnal cycle of wind speed and stratification, producing LLJs during each overnight period, however large discrepancies occurred at certain times for each simulation in relation to the observations. 54-member ensembles encompassing changes of the above discussed configuration parameters displayed a wide range of simulated vertical distributions of wind speed and direction, and potential temperature, reflecting highly variable representations of stratification during the weakly stable overnight conditions. Root mean square error (RMSE) statistics show that different ensemble members performed better and worse in various simulated parameters at different times, with no clearly superior configuration . Simulations using a PBL parameterization designed specifically for the stable conditions investigated herein provided superior overall simulations of wind speed at 80 m, demonstrating the efficacy of targeting improvements of physical process models in areas of known deficiencies. However, the considerable magnitudes of the RMSE values of even the best performing simulations indicate ample opportunities for further improvements.

Mirocha, Jeff D.↗

DOE OSTI2022

Performance Evaluation of the Meteorology and Air Quality Conditions From Multiscale WRF-CMAQ Simulations for the Long Island Sound Tropospheric Ozone Study (LISTOS)

The Long Island Sound (LIS) Tropospheric Ozone Study was a multi-agency collaborative field campaign conducted during the summer of 2018 to improve the understanding of ozone chemistry and transport from New York City to areas downstream, especially the LIS and adjacent Connecticut coastline. Measurements made during this campaign were leveraged to test and evaluate the coupled WRF-CMAQ model at 12 km, 4 and 1.33 km horizontal grid spacing. Special attention was placed on the model's representation of sea breeze circulations, low level jets, and boundary layer evolution. The evaluation suggests using higher resolutions resulted in improved surface meteorology statistics throughout the whole summer, with temperature biases seeing the biggest statistical improvements when using 1.33-km grid spacing, going from –0.12 to 0.08 K. Additionally, 4-km grid spacing provided the biggest advantage when simulating ozone over the region of interest, with biases being reduced from 2.40 to 0.57 to 0.37 ppbV with increased resolution. Case studies of two high ozone concentration events (July 10 and August 6) revealed that sound breezes and low-level jets had a critical role in transporting pollutant-rich, shallow marine air masses from the LIS inland over the Connecticut coast. Modifications were made to the representation of sea surface temperatures, which subsequently improved the simulation of surface ozone predictions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2016

The fan of influence of streams and channel feedbacks to simulated land surface water and carbon dynamics

Large-scale land models assume unidirectional land-to-river hydrological interactions, without considering feedbacks between channels and land. Using a tested, physically based model with explicit multiway interactions between overland, channel, wetland, and groundwater flows, we assessed how the representation and properties of channels influence simulated land surface hydrologic, biogeochemical, and ecosystem dynamics. A zone near the channels where various fluxes and states are significantly influenced by the channels, referred to as the fan of influence (FoI) of channels, has been identified. We elucidated two mechanisms inducing the model-derived FoI: the base flow mechanism, in which incised, gaining streams lower the water table and induce more base flow, and the relatively more efficient conveyance of the channel network compared to overland flow. We systematically varied drainage density and grid resolution to quantify the size of the FoI, which is found to span a large fraction of the watershed (25-50%) for hydrologic variables including depth to water table and recharge, etc. The FoI is more pronounced with low-resolution simulations but remains noticeable in hyperresolution (25 m) subbasin simulations. The FoI and the channel influence on basin-average fluxes are also similar in simulations with alternative parameter sets. We found that high-order, entrenched streams cause larger FoI. In addition, removing the simulated channels has disproportionally large influence on modeled wetland areas and inundation duration, which has implications for coupled biogeochemical or ecological modeling. Our results suggest that explicit channel representation provides important feedbacks to land surface dynamics which should be considered in meso or large-scale simulations. Since grid refinement incurs prohibitive computational cost, subgrid channel parameterization has advantages in efficiency over grid-based representations that do not distinguish between overland flow and channel flow.

58 GEOSCIENCES↗

NASA NTRSDate not supplied

FUN3D Analyses of the 5th AIAA Propulsion Aerodynamics Workshop Inlet Test Case

This work presents the results of FUN3D analyses that were performed for the 5th AIAA Propulsion Aerodynamics Workshop. The workshop was separated into two sections that focused on inlet and nozzle flows. This paper focuses on the inlet section of the workshop, which considered an axisymmetric inlet in close-proximity to the ground and subjected to a crosswind. Four levels of grid refinement were provided by the workshop committee. Participants were asked to perform simulations for three crosswind velocities on a minimum of three levels of grid refinement, for a total of nine required cases. This work employed an updated set of grids provided by the committee after completion of the workshop based on feedback from participants. The updated grids featured a larger computational domain and increased refinement near the ground vortex region. This work employed the three coarsest grid levels from the updated grid set. Note that an additional extension downstream of the inlet was added to the provided grids for this work to mitigate observed convergence issues. Additionally, a limited turbulence model study was performed for the highest value of crosswind velocity and the finest grid. The results show that grid convergence was only achieved for the case with the lowest value of crosswind velocity. However, the results were generally observed to improve, relative to the experimental data, with increasing grid refinement. The FUN3D predictions were observed to agree qualitatively with the provided experimental data. However, FUN3D exhibited a tendency to underpredict the minimum value of the radially averaged total pressure at the Aerodynamic Interface Plane, which is consistent with the general findings from the workshop. Overall, FUN3D was found to perform well relative to its peers for the PAW5 workshop inlet test case.

PAI↗

NASA NTRSDate not supplied

FUN3D Analyses of the 5th AIAA Propulsion Aerodynamics Workshop Inlet Test Case

PAI↗

DOE OSTI2025

Synchronous Machine Governor Upgrade

Conventional generation sources play a critical role in the stability and reliability of the electrical grid, particularly as we transition towards more renewable energy sources. To understand and accurately emulate their behavior for optimizing grid operations and ensuring seamless integration with renewable technologies, it is essential to better emulate the grid- and plant-level impacts of conventional generation sources, such as natural gas (NG) driven heat recovery steam generators (HRSGs) and combustion turbines (CTs). Therefore, a governor model is developed in a programmable logic controller (PLC) to investigate the performance of the conventional generator under various dynamic operating conditions and to identify the impact on grid stability in a controlled environment. The governor model aims to enable the hardware-in-the-loop (HIL) based emulation of these conventional generation sources using the existing 2 MVA synchronous machine/generator that is driven by a flexible 2.5 MW variable speed drive. This setup will allow us to replicate the dynamic characteristics and response behaviors of NG-driven HRSGs and CTs. The controls for the emulated conventional plants follow the industry standard and are adjustable, ensuring they accurately reflect the operational capabilities and limitations of real-world systems. These controls include load-following capabilities, ramp rates, startup and shutdown sequences, and emissions characteristics. By incorporating these adjustable controls, we aim to capture the nuanced impacts of conventional generation, such as their ability to provide ancillary services like frequency regulation, voltage support, and spinning reserve. In this report, we simulate two types of dynamic operations: grid-connected and islanding. For each dynamic operation, representative starting sequences are tested, including turbine purge, ignition, speed ramping up, generator excitation and synchronizing, and breaker close. The HIL based tests provides insights for field deployment, specifically the high-fidelity governor model provides results to predict the potential stability and reliability risk and suggest possible integration measures (e.g., generation and load balancing, tuning of governor control parameters). Ultimately, this enhanced emulation capability will be integrated into our Advanced Research on Integrated Energy Systems (ARIES), enabling us to conduct comprehensive studies on the interactions between conventional and renewable energy sources. By better understanding these interactions, we can develop strategies to optimize the overall performance and reliability of the grid. This will support the deployment of advanced grid management techniques, such as demand response, grid-forming inverters, and energy storage systems. The main contributions are summarized as follows: (1) This report introduces a PLC-based governor model for gas turbines. This model accurately simulates the dynamic behavior of conventional generation sources under various operational scenarios; (2) The model is integrated with an HIL testbed that includes a 2.5 MW variable speed drive and a 2 MVA synchronous machine. This setup enables realistic, real-time emulation of conventional power plants, particularly NG driven HRSGs and CTs; (3) The developed model is adaptable to various gas turbine configurations and allows for precise control over parameters such as MW ramp rates. This flexibility makes it a valuable tool for future research and industry collaboration; and (4) By incorporating the model into the National Renewable Energy Laboratory's Advanced Research on Integrated Energy Systems, the report lays the groundwork for future studies on interactions between conventional and renewable energy sources, enhancing the ability to develop advanced grid management strategies.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2023

EMT and Phasor-Domain Co-simulation of a Low Frequency AC Transmission System for Offshore Wind Integration

This work investigates system level interactions of multi-terminal low frequency high-voltage alternating current (LF-HVac) system using back-to back converters and grid-level 60 Hz system. A co-simulation framework by integrating high- fidelity dynamic model of multi-terminal LF-HVac system in electromagnetic transient (EMT) simulation platform and phasor model of IEEE 39 bus system in transient-stability (phasor) program. The co-simulation framework for multi-terminal LF- HVac system with offshore wind integrated in a bulk power system is developed successfully and its steady state and dynamic operation is demonstrated and compared with EMT simulation. The co-simulation results show accurate representation of the dynamic responses with reduced computation burden and time compared to EMT platform.

low-frequency transmission, co-simulation, electro↗

DOE OSTI2024

Newton-Raphson AC Power Flow Convergence Based on Deep Learning Initialization and Homotopy Continuation

Power flow forms the basis of many power system studies. With the increased penetration of renewable energy, grid planners tend to perform multiple power flow simulations under various operating conditions and not just selected snapshots at peak or light load conditions. Getting a converged AC power flow (ACPF) case remains a significant challenge for grid planners especially in large power grid networks. This paper proposes a two-stage approach to improve Newton-Raphson ACPF convergence and was applied to a 6102 bus Electric Reliability Council of Texas (ERCOT) system. The first stage utilizes a deep learning-based initializer with data re-training. Here a deep neural network (DNN) initializer is developed to provide better initial voltage magnitude and angle guesses to aid in power flow convergence. This is because Newton-Raphson ACPF is quite sensitive to the initial conditions and bad initialization could lead to divergence. The DNN initializer includes a data re-training framework that improves the initializer's performance when faced with limited training data. The DNN initializer successfully solved 3,285 cases out of 3,899 non-converging dispatch and performed better than random forest and DC power flow initialization methods. ACPF cases not solved in this first stage are then passed through a hot-starting algorithm based on homotopy continuation with switched shunt control. The hot-starting algorithm successfully converged 416 cases out of the remaining 614 non-converging ACPF dispatch. In conclusion, the combined two-stage approach achieved a 94.9% success rate, by converging a total of 3,701 cases out of the initial 3,899 unsolved cases.

Deep learning↗

DOE OSTI2022

Improving Predictive Capability in REHEDS Simulations with Fast, Accurate, and Consistent Non-Equilibrium Material Properties

Predictive design of REHEDS experiments with radiation-hydrodynamic simulations requires knowledge of material properties (e.g. equations of state (EOS), transport coefficients, and radiation physics). Interpreting experimental results requires accurate models of diagnostic observables (e.g. detailed emission, absorption, and scattering spectra). In conditions of Local Thermodynamic Equilibrium (LTE), these material properties and observables can be pre-computed with relatively high accuracy and subsequently tabulated on simple temperature-density grids for fast look-up by simulations. When radiation and electron temperatures fall out of equilibrium, however, non-LTE effects can profoundly change material properties and diagnostic signatures. Accurately and efficiently incorporating these non-LTE effects has been a longstanding challenge for simulations. At present, most simulations include non-LTE effects by invoking highly simplified inline models. These inline non-LTE models are both much slower than table look-up and significantly less accurate than the detailed models used to populate LTE tables and diagnose experimental data through post-processing or inversion. Because inline non-LTE models are slow, designers avoid them whenever possible, which leads to known inaccuracies from using tabular LTE. Because inline models are simple, they are inconsistent with tabular data from detailed models, leading to ill-known inaccuracies, and they cannot generate detailed synthetic diagnostics suitable for direct comparisons with experimental data. This project addresses the challenge of generating and utilizing efficient, accurate, and consistent non-equilibrium material data along three complementary but relatively independent research lines. First, we have developed a relatively fast and accurate non-LTE average-atom model based on density functional theory (DFT) that provides a complete set of EOS, transport, and radiative data, and have rigorously tested it against more sophisticated first-principles multi-atom DFT models, including time-dependent DFT. Next, we have developed a tabular scheme and interpolation methods that compactly capture non-LTE effects for use in simulations and have implemented these tables in the GORGON magneto-hydrodynamic (MHD) code. Finally, we have developed post-processing tools that use detailed tabulated non-LTE data to directly predict experimental observables from simulation output.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Object-oriented analysis as a foundation for building climate storylines of compounding short-term drought and crop heat stress

Introduction: Crops are vulnerable to precipitation and heat extremes during late spring through summer. Methods: We analyzed for a north-central U.S. region short-term drought and agricultural heat stress during April-May-June-July. We used the 4-km Parameter Elevation Regression on Independent Slopes Model (PRISM) for observations, aggregated to a 25-km grid, and two 25-km Regional Climate Model version 4 (RegCM4) simulns used either GFDL- or MPI-GCM boundary conditions. We chose 1981-2000 as our contemporary time period, and 2041- 2060 as our scenario time period, which used the Representative Concentration Pathway 8.5 emissions scenario. We used object-oriented analysis to identify events of interest in observations and simulations by identifying objects in a space-time domain that meet specified criteria, such as exceeding a heat-stress temperature threshold. The event diagnosis allowed analysis of compound events, occurring when temperature and drought objects overlap. Results: Identified objects yielded events that can undermine agricultural productivity and which are thus relevant to decision makers, making them building blocks for possible climate storylines. The observations and simulations showed similar spatial distributions of event frequencies across the analysis region. However, the simulations attained this distribution by having fewer events that tend to cover larger areas compared to observed events, suggesting that the effective resolution of the simulations was coarser than their 25-km grids. Short-term drought frequency increased and heat-stress frequency decreased in transitioning to the scenario climate. When compounding occurred heat-stress events generally preceded the short-term drought events. The overlapping, compound events tended to be more extreme compared to non-overlapping events of either type. Discussion: The information yielded projected changes in these agriculturally motivated events. One prominent conditional behavior emerging from the work was that a heat-stress event should be a warning to watch for potential drought, as both could compound each other to more intense levels.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Reconfigurable Real-Time Power Grid Emulator for Systems With High Penetration of Renewables

Novel power system control and new utility devices need to be tested before their actual deployment to the power grid. To assist with such a testing need, real-time digital emulators such as RTDS and Opal-RT can be used to connect to the physical world and form a hardware in the loop (HIL) emulation. However, due to the limitations of today's computational resources, the accuracy and fidelity suffer from different levels of model reductions in purely digital simulations. CURENT has developed a reconfigurable electric grid hardware testbed (HTB) to overcome the limitations of digital emulators. The HTB has been used to develop measurement, control, modeling, and actuation techniques for a national grid with a high penetration of renewables. The power electronic-based system includes emulators for synchronous generators; photovoltaics with grid-interfacing inverter; wind turbines; induction motor loads, ZIP loads, power electronic loads; batteries; ac and dc transmission lines; short circuit faults and grid relay protection; and a multiterminal HVDC overlay including power electronics interfaces. The system contains real elements of power flow, measurement, communication, protection, and control that mimic what would be seen in an actual electric grid. This paper presents an overview of the HTB and several scenarios that have been run to determine control and actions needed for the future power grid.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2026

LaueMatching: an approach for rapid and robust indexing of Laue diffraction patterns

Traditional Laue diffraction pattern indexing often struggles with noisy data, weak signals, peak overlap and missing reflections, particularly from complex or deformed microstructures. Here, we introduce LaueMatching, a high-throughput indexing algorithm designed to overcome these limitations. LaueMatching utilizes a fundamentally different approach based on direct pattern correlation: experimentally pre-processed images are compared against a comprehensive pre-computed library of simulated diffraction patterns corresponding to a dense grid of possible orientations. This approach bypasses the need for explicit peak identification and fitting, steps that are often a failure point for traditional methods. The algorithm rapidly and robustly indexes multiple crystallographic orientations and crystal systems simultaneously, even from challenging patterns. LaueMatching's effectiveness and accuracy have been rigorously tested and validated on diverse experimental (Ni, Al, EuAl 2 O 4 ) and simulated diffraction patterns, demonstrating high-fidelity orientation refinement. Code to implement this approach on both CPU and GPU resources can be downloaded from https://github.com/AdvancedPhotonSource/LaueMatching.

36 MATERIALS SCIENCE↗

DOE OSTI2015

SU-E-T-315: Dosimetric Effects of Couch Top Shift On VMAT Delivery in Absence of Indexing

Purpose: To investigate dosimetric effects of couch top shift for volumetric-modulated arc therapy (VMAT) in absence of indexing of immobilization devices. Methods: A total of twelve VMAT treatment plans were selected from three regions (lung, abdomen, and pelvis) to account for the variation of the patient position relative to the couch top. The treatment plans were generated using the Varian Eclipse system. A pinpoint ionization chamber (PTW TN31014) was placed at the center of 16-cm solid water phantom and the dose was delivered using the Varian TrueBeam STx with BrainLAB ExacTrac couch top. To simulate the day-to-day variation of the patient position relative to couch top, the couch top was laterally shifted up to 50 mm, with an increment of 5 mm from 0 to 20 mm; and of 10 mm afterwards, and the phantom was moved back to 0 cm shift for measurement. The dose was also delivered using a Varian tennis racket grid insert at 0 cm shift to simulate the absence of couch top. The treatment plans were delivered with 6, 10, and 15 MV photons using the same leaf sequencing to investigate the energy dependence of couch top shift. The dose difference was normalized to 0 cm shift for the regular couch top for comparison. Results: The percent difference of dose was found to increase with lateral shift for all energies; however, the average differences were close to 0% and the maximum difference was within 1% along the lateral shifts. The differences with the absence of couch top were 2.2±0.5% (6MV), 1.7±0.3% (10MV), and 1.6±0.2% (15MV), respectively. Conclusion: The inclusion of couch top is recommended in treatment planning to minimize the dosimetric uncertainty between calculated and delivered dose even in absence of indexing of immobilization devices in VMAT delivery.

60 APPLIED LIFE SCIENCES↗

NASA NTRS2018

Unstructured Grid Adaptation and Solver Technology for Turbulent Flows

Unstructured grid adaptation is a tool to control Computational Fluid Dynamics (CFD) discretization error. However, adaptive grid techniques have made limited impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Issues that prevent the use of adaptive grid methods are identified by applying unstructured grid adaptation methods to a series of benchmark cases. Once identified, these challenges to existing adaptive workflows can be addressed. Unstructured grid adaptation is evaluated for test cases described on the Turbulence Modeling Resource (TMR) web site, which documents uniform grid refinement of multiple schemes. The cases are turbulent flow over a Hemisphere Cylinder and an ONERA M6Wing. Adaptive grid force and moment trajectories are shown for three integrated grid adaptation processes with Mach interpolation control and output error based metrics. The integrated grid adaptation process with a finite element (FE) discretization produced results consistent with uniform grid refinement of fixed grids. The integrated grid adaptation processes with finite volume schemes were slower to converge to the reference solution than the FE method. Metric conformity is documented on grid/metric snapshots for five grid adaptation mechanics implementations. These tools produce anisotropic boundary conforming grids requested by the adaptation process.

Park, Michael A.↗