Generic Semi-Analytical Simulation Tool (GenSAS)
A generic simulation tool for power grid and general dynamic systems based on semi-analytical solutions (SAS) technology.
SEARCH · Search NASA
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.
A generic simulation tool for power grid and general dynamic systems based on semi-analytical solutions (SAS) technology.
Test-grid-buildouts is a set of notebooks for creating buildouts of synthetic transmission grids with various levels of renewable energy. Workflows can be used to augment various test grid with varying levels of wind/solar generation. Aggregated renewable energy time series (actuals and scenarios) at grid buses can then be used for various transmission grid operation problems. Highly customizable workflow consists of parsing of synthetic grid files, setting up SAM & reV simulations, application of exclusion masks, creating of custom renewable generators on the grid, aggregating of power time series, and writing of modified grid description files.
The South America Affinity Group (SAAG) was established in early 2019 by the National Center for Atmospheric Research (NCAR) Water Systems Program as a community effort focused on improving hydroclimate science over South America. SAAG supports large research efforts such as the ANDEX Regional Hydroclimate Program as well as individual research groups. The group started with a dozen members and quickly grew to over 100 participants from more than ten countries. For the past four years, the SAAG has been meeting online every two weeks and has organized sessions at international conferences such as the American Geophysical Union Fall Meeting and the Convection-Permitting Climate Workshop. At the core of the SAAG effort are two multi-decadal convection permitting (CP) model simulations with 4-km grid spacing for historical and future climates over the South American continent. Additionally, a major observational data collection effort has been undertaken, including in-situ station data from South American meteorological and water services, gridded products, satellite-based observations, and field campaign data. In conclusion, this article discusses the research needs and scientific goals that drive this community of scientists with diverse backgrounds and interests.
This summarizes technical highlights from the DOE Early Career Project entitled “Understanding the Roles of Cloud Microphysics and Land Surface Coupling Feedbacks in Multi-Scale Predictions of Central US Summer Hydroclimate” led by PI Mike Pritchard at the University of California. The project’s original aims were to investigate trade-offs of cloud superparameterization such as the effects of microphysical assumptions on mesoscale convective systems (Section I) and the representation of land-atmosphere coupling under explicit convection (Section II). Under the same broader theme of water cycle dynamics, superparameterization, and land interaction, the project evolved to further investigate the climate dynamics of irrigation, vegetation and streamflow (Section III), the dynamics of the Madden- Julian Oscillation (Section IV), as well as the use of modern machine learning to replace traditional sub grid approaches for climate simulation (Section V). Twenty peer-reviewed journal publications were produced as a result of this project (Section VI).
The Distribution System Operator with Transactive (DSO+T) study investigates the engineering and economic performance of a transactive energy retail market coordinating a high penetration of customer-side flexible energy assets. The study seeks to answer whether such an implementation is cost effective for customers, recovers sufficient revenue for DSOs, and is equally applicable and beneficial to a range of flexible asset types, renewable generation scenarios, and market assumptions. This report volume provides a detailed set of results for the DSO+T study extending results presented in Volumes 1, 2, and 4. The engineering and economic performance of the transactive energy scheme is presented for two separate flexible asset deployments: flexible loads (HVAC units and residential water heaters) and behind-the-meter batteries. The results of each transactive case are compared to a business-as-usual case. These cases are subject to two different renewable generation scenarios, a moderate renewable generation scenario, representative of current levels of renewable generation deployment, and a future high renewables scenario, including the increased deployment of rooftop solar photovoltaic and electric vehicles. The transactive coordination scheme is shown to produce effective and stable control and decrease peak loads 9–15%. The resulting annual demand flexibility provides net economic savings of $3.3–5.0B per year for a region the size of Texas. Detailed analysis shows that net benefits were seen for a range of distribution system operator, customer, and flexible asset types. Both participating customer (with transactive flexible assets) and nonparticipating customers (with nonflexible assets) see reductions in annual utility bills and net annual energy expenses in the range of 10–16%.
Weather models typically operate on averages of 1 km or more. This averages out fine scale turbulence. This results in smooth laminar flows and rapid vertical motions are usually averaged out. In this study, those fine scaled features were simulated for a 2013 field study of transport and diffusion of a tracer gas. Parameterized subgrid turbulence was required to produce dilution representative of observations. Despite this, the simulation with the finest grid and contained parameterized subgrid turbulence produced dilution that was consistent with observations.
Explore the source record for details and available documents.
We present a simple analytic procedure for generating atomic (H i) to molecular (H{sub 2}) density profiles for optically thick hydrogen gas clouds illuminated by far-ultraviolet radiation fields. Our procedure is based on the analytic theory for the structure of one-dimensional H i/H{sub 2} photon-dominated regions, presented by Sternberg et al. Depth-dependent atomic and molecular density fractions may be computed for arbitrary gas density, far-ultraviolet field intensity, and the metallicity-dependent H{sub 2} formation rate coefficient, and dust absorption cross section in the Lyman–Werner photodissociation band. We use our procedure to generate a set of H I-to-H{sub 2} transition profiles for a wide range of conditions, from the weak- to strong-field limits, and from super-solar down to low metallicities. We show that if presented as functions of dust optical depth, the H I and H{sub 2} density profiles depend primarily on the Sternberg “αG parameter” (dimensionless) that determines the dust optical depth associated with the total photodissociated H I column. We derive a universal analytic formula for the H I-to-H{sub 2} transition points as a function of just αG. Our formula will be useful for interpreting emission-line observations of H i/H{sub 2} interfaces, for estimating star formation thresholds, and for sub-grid components in hydrodynamics simulations.
We investigate the evolution of the cosmic web since z = 5 in grid-based cosmological hydrodynamical simulations, focusing on the mass and velocity fields of both baryonic and cold dark matter. The tidal tensor of density is used as the main method for web identification, with λ {sub th} = 0.2–1.2. The evolution trends in baryonic and dark matter are similar, although moderate differences are observed. Sheets appear early, and their large-scale pattern may have been set up by z = 3. In terms of mass, filaments supersede sheets as the primary collapsing structures from z ∼ 2–3. Tenuous filaments assembled with each other to form prominent ones at z < 2. In accordance with the construction of the frame of the sheets, the cosmic divergence velocity, v {sub div}, was already well-developed above 2–3 Mpc by z = 3. Afterwards, the curl velocity, v {sub curl}, grew dramatically along with the rising of filaments, becoming comparable to v {sub div}, for <2–3 Mpc at z = 0. The scaling of v {sub curl} can be described by the hierarchical turbulence model. The alignment between the vorticity and the eigenvectors of the shear tensor in the baryonic matter field resembles that in the dark matter field, and is even moderately stronger between ω and e{sub 1}, and ω and e{sub 3}. Compared with dark matter, there is slightly less baryonic matter found residing in filaments and clusters, and its vorticity developed more significantly below 2–3 Mpc. These differences may be underestimated because of the limited resolution and lack of star formation in our simulation. The impact of the change of dominant structures in overdense regions at z ∼ 2–3 on galaxy formation and evolution is shortly discussed.
We propose a new method for determining the thermal state of the intergalactic medium based on Voigt profile decomposition of the Lyα forest. The distribution of Doppler parameter and column density (b–N H i distribution) is sensitive to the temperature–density relation T = T 0 (ρ/ρ 0 ) γ-1 , and previous work has inferred T 0 and γ by fitting its low-b cutoff. This approach discards the majority of available data and is susceptible to systematics related to cutoff determination. We introduce an approach that exploits all information encoded in the b – N H i distribution by modeling its entire shape. We apply kernel density estimation to discrete absorption lines to generate model probability density functions, and then we use principal component decomposition to create an emulator that can be evaluated anywhere in thermal parameter space. We introduce a Bayesian likelihood based on these models enabling parameter inference via Markov Chain Monte Carlo. The method's robustness is tested by applying it to a large grid of thermal history simulations. By conducting 160 mock measurements, we establish that our approach delivers unbiased estimates and valid uncertainties for a 2D ( T 0 , γ) measurement. Furthermore, we conduct a pilot study applying this methodology to real observational data at z = 2. Using 200 absorbers, equivalent in path length to a single Lya forest spectrum, we measure $\mathrm{log}{T}_{0}={4.092}_{-0.055}^{+0.050}$ and $\gamma ={1.49}_{-0.074}^{+0.073}$ in excellent agreement with cutoff fitting determinations using the same data. Our method is far more sensitive than cutoff fitting, enabling measurements of log T 0 and γ with precision on $\mathrm{log}{T}_{0}$ (γ) nearly two (three) times higher for current data set sizes.
This paper systematically studies the relation between metallicity and mass loss of massive stars. We perform one-dimensional stellar evolution simulations and build a grid of ~2000 models with initial masses ranging between 11 and 60 M ⊙ and absolute metallicities Z between 0.00001 and 0.02. Steady-state winds, comprising hot main-sequence winds and cool supergiant winds, are the main drivers of the mass loss of massive stars in our models. We calculate the total mass loss over the stellar lifetime for each model. Our results reveal the existence of a critical metallicity Z c at Z ~ 10 –3 , where the mass loss exhibits a dramatic jump. If Z > Z c , massive stars tend to evolve into cool supergiants, and a robust cool wind is operational. In contrast, if Z < Z c , massive stars usually remain as blue supergiants, wherein the cool wind is not activated and the mass loss is generally weak. Moreover, we calculate the wind feedback in a 10 5 M ⊙ star cluster with the Salpeter initial mass function. The kinetic energy released by winds does not exhibit any significant transition at Z c because the wind velocity of a cool supergiant wind is low and contributes little to the kinetic energy. The effects of critical metallicity provide implications for the fates of metal-poor stars in the early universe.
The interaction of supernova ejecta with a surrounding circumstellar medium (CSM) generates a strong shock, which can convert ejecta kinetic energy into observable radiation. Given the diversity of potential CSM structures (arising from diverse mass-loss processes such as late-stage stellar outbursts, binary interaction, and winds), the resulting transients can display a wide range of light-curve morphologies. We provide a framework for classifying the transients arising from interaction with a spherical CSM shell. The light curves are decomposed into five consecutive phases, starting from the onset of interaction and extending through shock breakout and subsequent shock cooling. The relative prominence of each phase in the light curve is determined by two dimensionless quantities representing the CSM-to-ejecta mass ratio η, and the breakout parameter ξ. These two parameters define four light-curve morphology classes, where each class is characterized by the location of the shock breakout and the degree of deceleration as the shock sweeps up the CSM. We compile analytic scaling relations connecting the luminosity and duration of each light-curve phase to the physical parameters. We then run a grid of radiation hydrodynamics simulations for a wide range of ejecta and CSM parameters to numerically explore the landscape of interaction light curves, and to calibrate and confirm the analytic scalings. We connect our theoretical framework to several case studies of observed transients, highlighting the relevance in explaining slow-rising and superluminous supernovae, fast blue optical transients, and double-peaked light curves.
Experiments evaluating the uniformity and intensity of rainfall produced by the NGEE Arctic Rainfall Simulator (NARS) were conducted at Los Alamos National Laboratory, New Mexico, over summer 2022. Petri dishes were placed in a grid within the NARS plot. Simulated rainfall was collected in each petri dish and the intensity and uniformity of the simulator was subsequently calculated. This data package contains two .csv files, one that summarizes the rainfall intensity and uniformity for each experiment, the other that contains individual petri dish water volume and intensity for each plot location and experiment. The Python scripts to control NARS are also included. The NGEE Arctic Rainfall Simulator (NARS) is a variable intensity rainfall simulator (RFS) with a frame design based on the Humphry et al. (2002) RFS and a water delivery system based on the Walnut Gulch (Paige et al., 2004) RFS. The NARS uses an aluminum frame that is fully deconstructable for transportation to field locations and a water system that enables variable rain intensity. Rain intensity control and data collection are automated using a Raspberry Pi microcomputer. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a 15-year research effort (2012-2027) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).
This paper discusses the potential of using the Global Positioning System (GPS) to provide autonomous navigation capability to NASA satellites in the 1980 era. Some of the driving forces motivating autonomous navigation are presented. These include such factors as advances in attitude control systems, onboard science annotation, and onboard gridding of imaging data. Simulation results which demonstrate baseline orbit determination accuracies using GPS data on Seasat, Landsat-D, and the Solar Maximum Mission are presented. Emphasis is placed on identifying error sources such as GPS time, GPS ephemeris, user timing biases, and user orbit dynamics, and in a parametric sense on evaluating their contribution to the orbit determination accuracies.
An accurate numerical analysis of the flows associated with rotor-stator configurations in turbomachinery can be very helpful in optimizing the performance of turbomachinery. In this study the unsteady, thin-layer, Navier-Stokes equations are solved using a system of patched and overlaid grids for a rotor-stator configuration of an axial turbine. The equations necessary for an accurate transfer of information between the several grids are briefly described within the framework of an iterative, implicit algorithm. Results in the form of Mach number contours, time-averaged pressures, unsteady pressures, amplitudes, and phase are presented. The numerical results are also compared with experimental data and the agreement is found to be good.
A numerical algorithm based on a finite-volume explicit scheme with Runge-Kutta time integration of the Euler equations is presented for calculating high-speed three-dimensional flow about complex aerospace configurations. The use of enhancing factors such as artificial dissipative terms, enthalpy damping and local time-stepping are described. An algebraic method for generating quasi-three-dimensional computational grids for realistic aerospace configurations is presented. Computed results for various three-dimensional bodies at different Mach numbers and angles of attack have been obtained using the methods for grid-generation and flow simulation. Comparison of computed and experimental data for an advanced tactical aircraft-like configuration is presented, and a reasonable agreement of the data is noticed.
A new numerical algorithm is developed and applied to the direct simulation of compressible, homogeneous turbulent flows at low Mach numbers. A split method in time first solves a subset of the equations explicitly, followed by an implicit treatment of the pressure terms. The definition of an average sound speed over the entire field allows the implicit equations to be solved analytically, while relaxing the severe time limit imposed by the large sound speeds. Results from direct simulations on 96(3) grids provide a data base against which a new subgrid-scale model for compressible homogeneous turbulence is tested. This model reduces to the linear combination model in the absence of compressibility.
We describe a Cartesian grid strategy for the study of three dimensional inviscid flows about arbitrary geometries that uses both conventional and CAD/CAM surface geometry databases. Initial applications of the technique are presented. The elimination of the body-fitted constraint allows the grid generation process to be automated, significantly reducing the time and effort required to develop suitable computational grids for inviscid flowfield simulations.