Search NASA⌕ Search

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.

At least 901 records · Page 50

On the Dependence of Simulated Convection on Domain Size in CRMs

Abstract We present a heuristic model to explain the suppression of deep convection in convection‐resolving models (CRMs) with a small number of grid columns, such as those used in super‐parameterized or multi‐scale modeling framework (MMF) general circulation models (GCM) of the atmosphere. Domains with few grid columns require greater instability to sustain convection because they force a large convective fraction, driving strong compensating subsidence warming. Updraft dilution, which is stronger for reduced horizontal grid spacing, enhances this effect. Thus, suppression of deep convection in CRMs with few grid columns can be reduced by increasing grid spacing. Radiative‐convective equilibrium simulations using standalone CRM simulations with the System for Atmospheric Modeling (SAM) and using GCM‐coupled CRM simulations with the Energy Exascale Earth System Model (E3SM)‐MMF confirm the heuristic model results.

CRM↗

Simulating image coaddition with the Nancy Grace Roman Space Telescope – I. Simulation methodology and general results

The upcoming Nancy Grace Roman Space Telescope will carry out a wide-area survey in the near-infrared. A key science objective is the measurement of cosmic structure via weak gravitational lensing. Roman data will be undersampled, which introduces new challenges in the measurement of source galaxy shapes; a potential solution is to use linear algebra-based coaddition techniques such as imcom that combine multiple undersampled images to produce a single oversampled output mosaic with a desired ‘target’ point spread function (PSF). We present here an initial application of imcom to 0.64 square degrees of simulated Roman data, based on the Roman branch of the Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC) Data Challenge 2 (DC2) simulation. We show that imcom runs successfully on simulated data that includes features such as plate scale distortions, chip gaps, detector defects, and cosmic ray masks. We simultaneously propagate grids of injected sources and simulated noise fields as well as the full simulation. We quantify the residual deviations of the PSF from the target (the ‘leakage’), as well as noise properties of the output images; we discuss how the overall tiling pattern as well as Moiré patterns appear in the final leakage and noise maps. We include appendices on interpolation algorithms and the interaction of undersampling with image processing operations that may be of broader applicability. The companion paper (‘Paper II’) explores the implications for weak lensing analyses.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluating the resilience benefits of marine energy in microgrids

Marine energy resources could promote clean energy and resilience of coastal and island microgrids, and thus, these applications are a key future market for marine energy development. To demonstrate these benefits, this paper illustrates how inclusion of wave resources into energy resilience solutions can improve overall grid efficiency and sustainability, as well as maintain electricity supply during grid outages. The paper describes a case study evaluation of the potential to add wave energy to the Moloka'i grid as Hawaii strives to meet a 100% clean energy target. The Microgrid Component Optimization for Resilience tool is used to simulate operation in off-grid conditions and size different combinations of wave, solar photovoltaic (PV), wind, storage, and fuel resources required to meet resilience objectives. This research investigates how including wave resources in a microgrid contributes to reducing or eliminating biofuel generation, producing a zero-greenhouse gas emission profile in the latter case, and avoiding the over-sizing of PV and battery systems to accommodate periods of unavailability or high demand. Insight from this paper supports the value proposition of wave resources for future markets and informs the relationship between marine generators and microgrids or isolated grids.

16 TIDAL AND WAVE POWER↗

Forward Modeling of Atmospheric Carbon Dioxide in GEOS-5: Uncertainties Related to Surface Fluxes and Sub-Grid Transport

Forward GEOS-5 AGCM simulations of CO2, with transport constrained by analyzed meteorology for 2009-2010, are examined. The CO2 distributions are evaluated using AIRS upper tropospheric CO2 and ACOS-GOSAT total column CO2 observations. Different combinations of surface C02 fluxes are used to generate ensembles of runs that span some uncertainty in surface emissions and uptake. The fluxes are specified in GEOS-5 from different inventories (fossil and biofuel), different data-constrained estimates of land biological emissions, and different data-constrained ocean-biology estimates. One set of fluxes is based on the established "Transcom" database and others are constructed using contemporary satellite observations to constrain land and ocean process models. Likewise, different approximations to sub-grid transport are employed, to construct an ensemble of CO2 distributions related to transport variability. This work is part of NASA's "Carbon Monitoring System Flux Pilot Project,"

Pawson, Steven↗

Advanced Research on Integrated Energy Systems Cyber Range

As digital technologies expand to meet the needs of a more autonomous, interconnected, and advanced power system, new cybersecurity complexities and vulnerabilities arise. The ARIES Cyber Range enables the energy sector to evaluate these evolutions and validate cybersecurity solutions without impacting live systems. Combining power grid-scale hardware with emulation and simulation approaches, the ARIES Cyber Range can faithfully replicate modern energy systems - from grid physics to communication networks, and everything in between - with real-world fidelity. At NLR, researchers and partners are answering complex power system cybersecurity questions, examining emerging threats to the electric sector, and de risking new security technologies, all at a mission-relevant speed that keeps pace with rapidly evolving systems and hazards.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Power Hardware-in-the-Loop Interfaces for Inverter-Based Microgrid Experiments Including Transitions: Preprint

Microgrids continue to proliferate, particularly those that rely on inverter-based resources as the microgrid's voltage and frequency leader. Hardware-in-the-loop (HIL) simulations can de-risk microgrid deployments, and therefore we need power-hardware-in-the-loop (PHIL) interfaces for grid-forming (GFM) inverters to enable simulations of seamless transitions between islanded and grid-connected operation. Such microgrid interfaces are particularly challenging to model when the inverters need to switch modes, i.e., between grid-following (GFL) and GFM as the microgrid transitions between grid-connected and islanded operation. This paper presents two PHIL interfaces for such microgrids and the preliminary experimental results: One interface is suitable for those who have access to a power amplifier that can operate as a controllable current source, and the other is for those who need a solution that uses a voltage-mode power amplifier.

battery energy storage system↗

wrfhydro

This dataset contains a continuous WRF-Hydro simulation for the CACTI field campaign. The simulated period is 1-Aug-2018 through 21-Mar-2019 with hourly model output. The simulation is made with two grids that match the two outer LASSO-CACTI grids, with 7.5 km and 2.5 km grid spacings. Forcing data for WRF-Hydro was obtained from the ERA5-Land product. Values from this WRF-Hydro simulation were used to initialize the soil state for each LASSO-CACTI simulation.

54 ENVIRONMENTAL SCIENCES↗

Electromagnetic Transient (EMT) Simulation Algorithms for Evaluation of Large-Scale Extreme Fast Charging Systems (T&D Models)

Simulation of high-fidelity models of extreme fast charging (XFC) systems and large-area power grids with many XFCs can be time consuming in traditional simulators. Traditional simulators use a single method of discretization for all the components that results in imposing a large computational burden of inverting a large matrix as well as increased computations related to single method of discretization (that is typically a trapezoidal method). To overcome the problem of simulating large-area power grids with many XFCs, in this paper, advanced numerical simulation algorithms are applied for the first time together to reduce the dimension of matrix inversion. Here, the algorithms include numerical stiffness-based segregation, time constant-based segregation, clustering and aggregation on differential algebraic equations (DAEs), and multi-order integration approaches. These algorithms apply multiple discretization algorithms rather than a single discretization algorithm that further reduces the computational burden. The approaches mentioned here have resulted in speed-up of up to 18x in the simulation of a single distribution system with 15 XFCs and of up to 271x in the simulation of a transmission-distribution system with 300 XFCs in multiple distribution feeders with respect to conventional simulators (like power systems computer aided design [PSCAD]).

42 ENGINEERING↗

Distributed Real-Time Simulation and its Applications to Wind Energy Research

This paper presents a novel approach - distributed Real-Time Simulations (RTS) - for innovative wind energy and grid integration research. Distributed RTS is performed with Digital Real-Time Simulators (DRTS) in different geographical locations and connected over the Internet. Geographically distributed RTS supports - higher computation capability for RTS, leveraging distributed assets, novel testing techniques and hybridized smart grid testbeds. This paper discusses some of the basic aspects of distributed RTS with the communication delay as a probabilistic variable. To leverage hardware resources and enable flexible collaboration of national laboratories, the distributed RTS connection between the simulators at Idaho National Laboratory (INL) and National Wind Test Center (NWTC) at National Renewable Energy Laboratory (NREL) is established. The main objective of this distributed lab setup is to perform wind turbine characterization and assessment of wind integration studies in large scale grids. Several systematic steps based on integrating models are carried out to assess the impact of communication delay on distributed RTS.

17 WIND ENERGY↗

Hybrid Reynolds-Averaged/Large Eddy Simulation of a Cavity Flameholder; Assessment of Modeling Sensitivities

Steady-state and scale-resolving simulations have been performed for flow in and around a model scramjet combustor flameholder. The cases simulated corresponded to those used to examine this flowfield experimentally using particle image velocimetry. A variety of turbulence models were used for the steady-state Reynolds-averaged simulations which included both linear and non-linear eddy viscosity models. The scale-resolving simulations used a hybrid Reynolds-averaged / large eddy simulation strategy that is designed to be a large eddy simulation everywhere except in the inner portion (log layer and below) of the boundary layer. Hence, this formulation can be regarded as a wall-modeled large eddy simulation. This effort was undertaken to formally assess the performance of the hybrid Reynolds-averaged / large eddy simulation modeling approach in a flowfield of interest to the scramjet research community. The numerical errors were quantified for both the steady-state and scale-resolving simulations prior to making any claims of predictive accuracy relative to the measurements. The steady-state Reynolds-averaged results showed a high degree of variability when comparing the predictions obtained from each turbulence model, with the non-linear eddy viscosity model (an explicit algebraic stress model) providing the most accurate prediction of the measured values. The hybrid Reynolds-averaged/large eddy simulation results were carefully scrutinized to ensure that even the coarsest grid had an acceptable level of resolution for large eddy simulation, and that the time-averaged statistics were acceptably accurate. The autocorrelation and its Fourier transform were the primary tools used for this assessment. The statistics extracted from the hybrid simulation strategy proved to be more accurate than the Reynolds-averaged results obtained using the linear eddy viscosity models. However, there was no predictive improvement noted over the results obtained from the explicit Reynolds stress model. Fortunately, the numerical error assessment at most of the axial stations used to compare with measurements clearly indicated that the scale-resolving simulations were improving (i.e. approaching the measured values) as the grid was refined. Hence, unlike a Reynolds-averaged simulation, the hybrid approach provides a mechanism to the end-user for reducing model-form errors.

Baurle, R. A.↗

Segmented Domain Decomposition Multigrid For 3-D Turbomachinery Flows

A Segmented Domain Decomposition Multigrid (SDDMG) procedure was developed for three-dimensional viscous flow problems as they apply to turbomachinery flows. The procedure divides the computational domain into a coarse mesh comprised of uniformly spaced cells. To resolve smaller length scales such as the viscous layer near a surface, segments of the coarse mesh are subdivided into a finer mesh. This is repeated until adequate resolution of the smallest relevant length scale is obtained. Multigrid is used to communicate information between the different grid levels. To test the procedure, simulation results will be presented for a compressor and turbine cascade. These simulations are intended to show the ability of the present method to generate grid independent solutions. Comparisons with data will also be presented. These comparisons will further demonstrate the usefulness of the present work for they allow an estimate of the accuracy of the flow modeling equations independent of error attributed to numerical discretization.

Celestina, M. L.↗

Model of Atmospheric Links on Optical Communications from High Altitude

Optical communication links have the potential to solve many of the problems of current radio and microwave links to satellites and high-altitude aircraft. The higher frequency involved in optical systems allows for significantly greater signal bandwidth, and thus information transfer rate, in excess of 10 Gbps, and the highly directional nature of laser-based signals eliminates the need for frequency-division multiplexing seen in radio and microwave links today. The atmosphere, however, distorts an optical signal differently than a microwave signal. While the ionosphere is one of the most significant sources of noise and distortion in a microwave or radio signal, the lower atmosphere affects an optical signal more significantly. Refractive index fluctuations, primarily caused by changes in atmospheric temperature and density, distort the incoming signal in both deterministic and nondeterministic ways. Additionally, suspended particles, such as those in haze or rain, further corrupt the transmitted signal. To model many of the atmospheric effects on the propagating beam, we use simulations based on the beam-propagation method. This method, developed both for simulation of signals in waveguides and propagation in atmospheric turbulence, separates the propagation into a diffraction and refraction problem. The diffraction step is an exact solution, within the limits of numerical precision, to the problem of propagation in free space, and the refraction step models the refractive index variances over a segment of the propagation path. By applying refraction for a segment of the propagation path, then diffracting over that same segment, this method forms a good approximation to true propagation through the atmospheric medium. Iterating over small segments of the total propagation path gives a good approximation to the problem of propagation over the entire path. Parameters in this model, such as initial beam profile and atmospheric constants, are easily modified in a simulation such as this, which allows for the rapid analysis of different propagation scenarios. Therefore, this method allows the development of a near-optimal system design for a wide range of situations, typical of what would be seen in different atmospheric conditions over a receiving ground station. A simulation framework based upon this model was developed in FORTRAN, and for moderate grid sizes and propagation distances these simulations are computable in reasonable time on a standard workstation. This presentation will discuss results thus far.

Subich, Christopher↗

The Role of Dredge-up in Double White Dwarf Mergers

We present the results of an investigation of the dredge-up and mixing during the merger of two white dwarfs (WDs) with different chemical compositions by conducting hydrodynamic simulations of binary mergers for three representative mass ratios. In all the simulations, the total mass of the two WDs is $\lesssim$1.0 M ⊙ . Mergers involving a CO and a He WD have been suggested as a possible formation channel for R Coronae Borealis (RCB)–type stars, and we are interested in testing if such mergers lead to conditions and outcomes in agreement with observations. Even if the conditions during the merger and subsequent nucleosynthesis favor the production of 18 O, the merger must avoid dredging up large amounts of 16 O, or else it will be difficult to produce sufficient 18 O to explain the oxygen ratio observed to be of order unity. Here, we performed a total of nine simulations using two different grid-based hydrodynamics codes using fixed and adaptive meshes and one smooth particle hydrodynamics (SPH) code. We find that in most of the simulations, >10 -2 M ⊙ of 16 O is indeed dredged up during the merger. However, in SPH simulations where the accretor is a hybrid He/CO WD with a ~0.1 M ⊙ layer of helium on top, we find that no 16 O is being dredged up, while in the q = 0.8 simulation <10 -4 M ⊙ of 16 O has been brought up, making a WD binary consisting of a hybrid CO/He WD and a companion He WD an excellent candidate for the progenitor of RCB stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Decentralized Distribution System Restoration with Grid-Forming/Following Inverter-Based Resources

The high penetration of distributed energy resources (DERs) in active distribution systems has posed challenges to the centralized distribution system restoration (DSR) strategies in current practice. On the other hand, the advancement in smart inverter technologies enables the bottom-up restoration capability. This paper is motivated to develop a 3-layered hierarchical framework for decentralized DSR, based on the grid-forming (GFM) and grid-following (GFL) grid-edge inverters. The first layer presents the tertiary control, which determines the load pickup schedule and generation dispatch of DERs, using the alternating direction method of the multipliers algorithm. The second layer consists of two control functions: GFM control, which regulates voltage and frequency, establishing a stable grid for GFL inverters to follow; and GFL control, which regulates the real and reactive power. In the third layer, the primary control is proposed to regulate the inverter voltage and current, which is developed based on the virtual oscillator control (VOC). Furthermore, the developed framework is tested in the modified IEEE 13-node test feeder. Two scenarios of grid-connected and islanded operating modes are designed, and simulation results demonstrate the effectiveness of decentralized DSR strategies for controlling grid-edge inverters to enhance the distribution system resilience.

24 POWER TRANSMISSION AND DISTRIBUTION↗

ASGarD: Adaptive Sparse Grid Discretization

Many areas of science exhibit physical processes that are described by high dimensional partial differential equations (PDEs), e.g., the 4D, 5D and 6D models describing magnetized fusion plasmas, models describing quantum chemistry, or derivatives pricing. Such problems are affected by the so-called “curse of dimensionality” where the number of degrees of freedom (or unknowns) required to be solved for scales as N D where N is the number of grid points in any given dimension D. A simple, albeit naive, 6D example is demonstrated in the left panel of Figure 1. With N = 1000 grid points in each dimension, the memory required just to store the solution vector, not to mention forming the matrix required to advance such a system in time, would exceed an exabyte - and also the available memory on the largest of supercomputers available today. The right panel of Figure 1 demonstrates potential savings for a range of problem dimensionalities and grid resolution. While there are methods to simulate such high-dimensional systems, they are mostly based on Monte-Carlo methods, which rely on a statistical sampling such that the resulting solutions include noise. Since the noise in such methods can only be reduced at a rate proportional to $\sqrt{N_p}$ where N p is the number of Monte-Carlo samples, there is a need for continuum, or grid/mesh-based methods for high-dimensional problems, which both do not suffer from noise and bypass the curse of dimensionality. We present a simulation framework that provides such a method using adaptive sparse grids.

97 MATHEMATICS AND COMPUTING↗

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work have been limited. Here, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. Boiling involves both fluid flow and heat transfer and thus requires the solution of the Navier-Stokes and the energy equations. The numerical method is based on writing one set of governing transport equations which is valid in both the liquid and vapor phases. This local, single-field formulation incorporates the effect of the interface in the governing equations as source terms acting only at the interface. These sources account for surface tension and latent heat in the equations for conservation of momentum and energy as well as mass transfer across the interface due to phase change. The single-field formulation naturally incorporates the correct mass, momentum and energy balances across the interface. Integration of the conservation equations across the interface directly yields the jump conditions derived in the local instant formulation for two-phase systems. In the numerical implementation, the conservation equations for the whole computational domain (both vapor and liquid) are solved using a stationary grid and the phase boundary is followed by a moving unstructured two-dimensional grid. While two-dimensional simulations have been used for preliminary studies and to examine the resolution requirement, the focus is on fully three-dimensional simulations. The numerical methodology, including the parallelization and grid refinement strategy is discussed, and preliminary results shown. For buoyancy driven flow, the heat transfer is in good agreement with experimental correlations. The changes when gravity is turned off and/or fluid shear is added are discussed, as well as the difference between simulations of a layer freely releasing bubbles versus simulations using only one wavelength initial perturbation. Figure 1 shows the early stages of the formation of a three-dimensional bubble from a thin vapor layer. The boundary conditions are periodic in the x and y direction, the bottom is a hot and the top allows a free outflow. The jagged edge of the surface close to the bottom of the computational domain is due to some of the surface elements being on the other side of the domain and some elements not plotted by our plotting routine. In the second figure, we show the temperature distribution through two perpendicular planes.

Tryggvason, G.↗

Performance of explicit and IMEX MRI multirate methods on complex reactive flow problems within modern parallel adaptive structured grid frameworks

Large-scale multiphysics simulations are computationally challenging due to the coupling of multiple processes with widely disparate time scales. The advent of exascale computing systems exacerbates these challenges since these systems enable ever-increasing size and complexity. In recent years, there has been renewed interest in developing multirate methods as a means to handle the large range of time scales, as these methods may afford greater accuracy and efficiency than more traditional approaches of using implicit-explicit (IMEX) and low-order operator splitting schemes. However, to date there have been few performance studies that compare different classes of multirate integrators on complex application problems. In this work, we study the performance of several newly developed multirate infinitesimal (MRI) methods, implemented in the SUNDIALS solver package, on two reacting flow model problems built on structured mesh frameworks. The first model revisits prior work on a compressible reacting flow problem with complex chemistry that is implemented using BoxLib but where we now include comparisons between a new explicit MRI scheme with the multirate spectral deferred correction (SDC) methods in the original paper. The second problem uses the same complex chemistry as the first problem, combined with a simplified flow model, but runs at a large spatial scale where explicit methods become infeasible due to stability constraints. Two recently developed IMEX MRI multirate methods are tested. These methods rely on advanced features of the AMReX framework on which the model is built, such as multilevel grids and multilevel preconditioners. The results from these two problems show that MRI multirate methods can offer significant performance benefits on complex multiphysics application problems and that these methods may be combined with advanced spatial discretization to compound the advantages of both.

97 MATHEMATICS AND COMPUTING↗

Arctic ice-ocean interactions in an 8-to-2 kilometer resolution global model

In the last decades, the Arctic climate has changed dramatically, with the loss of multiyear sea ice one of the clearest consequences. These changes have occurred on relatively rapid timescales, and both accurate short-term Arctic prediction (e.g., 10 days to three months) and climate projection of future Arctic scenarios present ongoing challenges. Here we describe a representation of the Arctic ocean and sea ice in a ultrahigh resolution simulation in which the horizontal grid mesh reduces from 8 km at the equator to 2 km at the poles (UH8to2) for the years 2017-2020. We find the simulation reproduces observed distributions of seasonal sea-ice thickness and concentration realistically, although concentration is biased low in the spring and summer and low biases in thickness are found in the central and eastern basins in the fall. Volume, fresh water, and heat transports through key passages are realistic, lying within observationally determined ranges. Climatological comparisons reveal that the UH8to2 Atlantic Water is shallower, warmer, and saltier than the World Ocean Atlas 2018 climatology for 2005-2017 in the eastern basin. Our analysis suggests that these biases, combined with a lack of stratification in the upper 100 m of the simulated ocean, contribute to the winter biases in modeled sea ice thickness. This relationship between biases in the sea ice and ocean points to a potential positive feedback within the model, illuminating challenges for long term model predictive power in a changing Arctic climate.

58 GEOSCIENCES↗