SEARCH · Search NASA

Results for “GRID MODERNIZATION”

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,009 records · Page 56

All sources NASA NTRS DOE OSTI

DOE OSTI2016

MAGNETO-FRICTIONAL MODELING OF CORONAL NONLINEAR FORCE-FREE FIELDS. I. TESTING WITH ANALYTIC SOLUTIONS

We report our implementation of the magneto-frictional method in the Message Passing Interface Adaptive Mesh Refinement Versatile Advection Code (MPI-AMRVAC). The method aims at applications where local adaptive mesh refinement (AMR) is essential to make follow-up dynamical modeling affordable. We quantify its performance in both domain-decomposed uniform grids and block-adaptive AMR computations, using all frequently employed force-free, divergence-free, and other vector comparison metrics. As test cases, we revisit the semi-analytic solution of Low and Lou in both Cartesian and spherical geometries, along with the topologically challenging Titov–Démoulin model. We compare different combinations of spatial and temporal discretizations, and find that the fourth-order central difference with a local Lax–Friedrichs dissipation term in a single-step marching scheme is an optimal combination. The initial condition is provided by the potential field, which is the potential field source surface model in spherical geometry. Various boundary conditions are adopted, ranging from fully prescribed cases where all boundaries are assigned with the semi-analytic models, to solar-like cases where only the magnetic field at the bottom is known. Our results demonstrate that all the metrics compare favorably to previous works in both Cartesian and spherical coordinates. Cases with several AMR levels perform in accordance with their effective resolutions. The magneto-frictional method in MPI-AMRVAC allows us to model a region of interest with high spatial resolution and large field of view simultaneously, as required by observation-constrained extrapolations using vector data provided with modern instruments. The applications of the magneto-frictional method to observations are shown in an accompanying paper.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

The Lyman α forest in optically thin hydrodynamical simulations

We study the statistics of the Lyα forest in a flat Λ cold dark matter cosmology with the N-body + Eulerian hydrodynamics code NYX.We produce a suite of simulations, covering the observationally relevant redshift range 2 ≤ z ≤ 4. We find that a grid resolution of 20 h-1kpc is required to produce 1 per cent convergence of Lyα forest flux statistics, up to k = 10 h -1 Mpc. In addition to establishing resolution requirements, we study the effects of missing modes in these simulations, and find that box sizes of L > 40h -1 Mpc are needed to suppress numerical errors to a sub-per cent level. Our optically thin simulations with the ionizing background prescription of Haardt & Madau reproduce an intergalactic medium density-temperature relation with T 0 ≈ 10 4 K and γ ≈ 1.55 at z = 2, with a mean transmitted flux close to the observed values. When using the ionizing background prescription of Faucher-Giguère et al., the mean flux is 10-15 per cent below observed values at z = 2, and a factor of 2 too small at z = 4. We show the effects of the common practice of rescaling optical depths to the observed mean flux and how it affects convergence rates. We also investigate the practice of 'splicing' results from a number of different simulations to estimate the 1D flux power spectrum and show it is accurate at the 10 per cent level. Finally, we find that collisional heating of the gas from dark matter particles is negligible in modern cosmological simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2001

Coupled-Flow Simulation of HP-LP Turbines Has Resulted in Significant Fuel Savings

Our objective was to create a high-fidelity Navier-Stokes computer simulation of the flow through the turbines of a modern high-bypass-ratio turbofan engine. The simulation would have to capture the aerodynamic interactions between closely coupled high- and low-pressure turbines. A computer simulation of the flow in the GE90 turbofan engine's high-pressure (HP) and low-pressure (LP) turbines was created at GE Aircraft Engines under contract with the NASA Glenn Research Center. The three-dimensional steady-state computer simulation was performed using Glenn's average-passage approach named APNASA. The areas upstream and downstream of each blade row mutually interact with each other during engine operation. The embedded blade row operating conditions are modeled since the average passage equations in APNASA actively include the effects of the adjacent blade rows. The turbine airfoils, platforms, and casing are actively cooled by compressor bleed air. Hot gas leaks around the tips of rotors through labyrinth seals. The flow exiting the high work HP turbines is partially transonic and, therefore, has a strong shock system in the transition region. The simulation was done using 121 processors of a Silicon Graphics Origin 2000 (NAS 02K) cluster at the NASA Ames Research Center, with a parallel efficiency of 87 percent in 15 hr. The typical average-passage analysis mesh size per blade row was 280 by 45 by 55, or approx.700,000 grid points. The total number of blade rows was 18 for a combined HP and LP turbine system including the struts in the transition duct and exit guide vane, which contain 12.6 million grid points. Design cycle turnaround time requirements ran typically from 24 to 48 hr of wall clock time. The number of iterations for convergence was 10,000 at 8.03x10(exp -5) sec/iteration/grid point (NAS O2K). Parallel processing by up to 40 processors is required to meet the design cycle time constraints. This is the first-ever flow simulation of an HP and LP turbine. In addition, it includes the struts in the transition duct and exit guide vanes.

Veres, Joseph P.↗

DOE OSTI2025

A census of fish passage facilities at US hydropower developments across the conterminous United States

Hydropower provides reliable and secure electricity and contributes significantly to the flexibility and stability of the US electricity grid. Hydropower dams generally are aquatic barriers and hazards to migratory fish in rivers. Fish passage facilities mitigate risks from hydropower to migratory fish, but information on these facilities is incomplete for the conterminous United States (CONUS). Here, we present the first CONUS-scale dataset of fish passage facilities at US hydropower developments in over 30 years. The existence of fish passage facilities (presence or absence) was specified for 1909 hydropower features, 390 of which had at least one facility. Most features had a single passage facility that provided passage in only one direction, with downstream being more common than upstream passage. Bypasses and ladders accounted for 60 % of the fish passage facility types. While we documented 659 fish passage facilities across CONUS, facilities were most common in the New England, Pacific Northwest, Mid-Atlantic, and Great Lakes regions. In general, fish passage facilities were more common at features that were closer to the ocean, at lower elevations, and at shorter dams, but not related to installed electrical generation capacity. Hydrologic sub-basins containing salmonids also contained the largest number of hydropower features, but the proportion of features with passage was generally higher in sub-basins containing multiple migratory taxa. This census provides valuable information on existing fish passage mitigation and is a benchmark to gauge progress toward a modernized hydropower fleet that provides affordable, reliable energy while protecting fishery resources and river ecosystems.

13 HYDRO ENERGY↗

DOE OSTI2016

EMAPS: An Efficient Multiscale Approach to Plasma Systems with Non-MHD Scale Effects

Global interactions of energetic ions with magnetoplasmas and neutral gases lie at the core of many space and laboratory plasma phenomena ranging from solar wind entry into and transport within planetary magnetospheres and exospheres to fast-ion driven instabilities in fusion devices to astrophysics-in-lab experiments. The ability of computational models to properly account for physical effects that underlie such interactions, namely ion kinetic, ion cyclotron, Hall, collisional and ionization processes is important for the success and planning of experimental research in plasma physics. Understanding the physics of energetic ions, in particular their nonlinear resonance interactions with Alfvén waves, is central to improving the heating performance of magnetically confined plasmas for future energy generation. Fluid models are not adequate for high-beta plasmas as they cannot fully capture ion kinetic and cyclotron physics (e.g., ion behavior in the presence of magnetic nulls, shock structures, plasma interpenetration, etc.). Recent results from global reconnection simulations show that even in a MHD-like regime there may be significant differences between kinetic and MHD simulations. Therefore, kinetic modeling becomes essential for meeting modern day challenges in plasma physics. The hybrid approximation is an intermediate approximation between the fluid and fully kinetic approximations. It eliminates light waves, removes the electron inertial temporal and spatial scales from the problem and enables full-orbit ion kinetics. As a result, hybrid codes have become effective tools for exploring ion-scale driven phenomena associated with ion beams, shocks, reconnection and turbulence that control the large-scale behavior of laboratory and space magnetoplasmas. A number of numerical issues, however, make three-dimensional (3D) large-scale hybrid simulations of inhomogeneous magnetized plasmas prohibitively expensive or even impossible. To resolve these difficulties we have developed a novel Event-driven Multiscale Asynchronous Parallel Simulation (EMAPS) technology that replaces time stepping with self-adaptive update events. Local calculations are carried out only on an “as needed basis”. EMAPS (i) guarantees accurate and stable processing of physical variables in time accurate simulations, and (ii) eliminates unnecessary computation. Applying EMAPS to the hybrid model has resulted in the development of a unique parallel code, dimension-independent (compile-time-configurable) HYPERS (Hybrid Parallel Event-Resolved Simulator) that scales to hundreds of thousands of parallel processors. HYPERS advances electromagnetic fields and particles asynchronously on time scales determined by local physical laws and mesh properties. To achieve high computational accuracy in complex device geometries, HYPERS employs high-fidelity Cartesian grids with masked conductive cells. The HYPERS model includes multiple ion species, energy and momentum conserving ion-ion collisions, and provides a number of approximations for plasma resistivity and vacuum regions. Both local and periodic boundary conditions are allowed. The HYPERS solver preserves zero divergence of magnetic field. The project has demonstrated HYPERS capabilities on a number of applications of interest to fusion and astrophysical plasma physics applications listed below. 1. Theta-pinch formation of FRCs The formation, spontaneous spin-up, and stability of theta-pinch formed field-reversed configurations have been studied self-consistently in 3D. The end-to-end hybrid simulations reveal poloidal profiles of implosion-driven fast toroidal plasma rotation and demonstrate three discharge regimes as a function of experimental parameters: the decaying stable configuration, the tilt unstable configuration, and the nonlinear evolution of a fast growing tearing mode. 2. FRC collisions with magnetic mirrors Interactions of fast plasma streams and objects with magnetic obstacles (dipoles, mirrors, etc) lie at the core of many space and laboratory plasma phenomena ranging from magnetoshells and solar wind interactions with planetary magnetospheres to compact fusion plasmas. HYPERS simulations are compared with data from the MSX experiment (LANL) that focuses on the physics of magnetized collisionless shocks through the acceleration and subsequent stagnation of FRC plasmoids against a strong magnetic mirrors and flux-conserving boundaries. 3. Exploding magnetoplasmas Results from hybrid simulations of two experiments at the LAPD and Nevada Terawatt Facility are discussed where short-pulse lasers are used to ablate solid targets to produce plasmas that expand across external magnetic fields. The first simulation recreates flutelike density striations observed at the leading edge of a carbon plasma and predicts an early destruction of the magnetic cavity in agreement with experimental evidence. In the second simulation a polyethylene target is ablated into a mixture of protons and carbon ions. A mechanism is demonstrated that allows protons to penetrate the magnetic field in the form of a collimated flow. The results are compared to experimental data and single-fluid MHD simulations. The EMAPS framework has the potential for wide application in many other engineering and scientific fields, such as climate models, biological systems, electronic devices, seismic events, oil reservation simulators that all involve advancing solutions of partial differential equations in time where the rate of activity can be adapted widely over the spatial domain depending on locally space/time phenomena (“events”).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2013

Regional Climate and Variability of the Summertime Continental United States in Reanalyses

Understanding climate variability at regional scales is an important for research and societal needs. Atmospheric retrospective-analyses (or reanalyses) integrate multitudes of observing systems with numerical models to produce continuous data that include variables not easily observed, if at all. The breadth of variables as well as observational influence included in reanalyses make them ideal for investigating climate variability. In this paper, we assess NASA s Modern Era Retrospective-analysis for Research and Application (MERRA) regional variability in North America, specifically the United States, in conjunction with current satellite data reanalyses. Emphasis is placed on summertime precipitation because 1) it is a difficult parameter to capture in the most difficult season, 2) significant observational resources exist to benchmark comparisons, and 3) accurate assessment of precipitation variability is crucial to a multitude of sectors and applications. Likewise, we have also begun to evaluate surface air temperature. While precipitation biases are identified, year to year variability of the precipitation variations, in many cases, are quite reasonable. However, some spurious long term trends and sudden shifts in the time series are also identified. In surface air temperature, analysis of station observations provides ERA Interim a clear overall advantage. However, in a number of regions, all the reanalyses are quite comparable in variability and trend. In other regions, significant precipitation biases may occur, which has implications for the ancillary process data in a reanalysis, such as surface fluxes. We also characterize the reanalyses ability to capture variability related to ENSO. In general, the summertime variations of precipitation in the reanalyses are more highly correlated (positively) to ENSO (using ENSO34) than are the observations. The Northwestern US shows the largest positive correlations to ENSO34, and reanalyses agree with that, and correlate highly with the gauge observations there. Additional evaluation of the data assimilation using the MERRA Gridded Innovations and Observations (GIO) data, which consists of the assimilated observations as well as forecast an analysis error fields. This work represents an initial evaluation of using the MERRA reanalysis to study regional climate variations of the United States, and identifying its applied limitations

Bosilovich, M G.↗

NASA NTRS2018

Enabling Reanalysis Research Using the Collaborative Reanalysis Technical Environment (CREATE)

Modern atmospheric and oceanic reanalysis are valuable assets for atmospheric research and climate monitoring (Kalnay et al. 1996). Now that most reanalysis records are more than 36 years long, the data have become more useful for climate modeling research (Dole et al. 2008). For investigators who need to use multiple reanalysis, a common challenge is that the data are distributed at various sites and often in different formats. The NASA CREATE system provides access to the data in one location in a standard format (one variable per file and standardized metadata in the CMIP5 style; see Table 1 for a list of the key acronyms used in this paper). The collection includes monthly and 6-hourly data from the seven major atmospheric reanalysis: CFSR (Saha et al. 2010), ERA-Interim (Dee et al. 2011), MERRA (Rienecker et al. 2011), MERRA-2 (Gelaro et al. 2017), JRA-25 (Onogi et al. 2007), JRA-55 (Kobayashi et al. 2014), and 20CRv2c (Compo et al. 2011). An ancillary portion of CREATE includes eight ocean reanalysis: NCEP CFSR, CMCC C-GLORSv5 (Storto et al. 2016), ECMWF ORAS4 (Balmaseda et al.2013), ECMWF ORAP5.0 (Zuo et al. 2015), University of Hamburg GECCO2 (Köhl 2015), GFDL ECDA (Zhang et al. 2007), NOAA GODAS (Saha et al. 2010), and MOVE/MRI.COM-G2i (Toyoda et al. 2016). The ocean state variables were similarly reformatted but were then also regridded onto a common horizontal and vertical grid. This approach facilitated the calculation of an ensemble average and spread that is also published alongside the native gridded data. A third reanalysis product is a global hourly 0.5° land surface air temperature dataset constructed by Wang and Zeng (2013). All three datasets are distributed through the ESGF in a format consistent with the CMIP style described by Cinquini et al. (2014).

Potter, Gerald L.↗

NASA NTRS1993

Algorithm and code development for unsteady three-dimensional Navier-Stokes equations

In the last two decades, there have been extensive developments in computational aerodynamics, which constitutes a major part of the general area of computational fluid dynamics. Such developments are essential to advance the understanding of the physics of complex flows, to complement expensive wind-tunnel tests, and to reduce the overall design cost of an aircraft, particularly in the area of aeroelasticity. Aeroelasticity plays an important role in the design and development of aircraft, particularly modern aircraft, which tend to be more flexible. Several phenomena that can be dangerous and limit the performance of an aircraft occur because of the interaction of the flow with flexible components. For example, an aircraft with highly swept wings may experience vortex-induced aeroelastic oscillations. Also, undesirable aeroelastic phenomena due to the presence and movement of shock waves occur in the transonic range. Aeroelastically critical phenomena, such as a low transonic flutter speed, have been known to occur through limited wind-tunnel tests and flight tests. Aeroelastic tests require extensive cost and risk. An aeroelastic wind-tunnel experiment is an order of magnitude more expensive than a parallel experiment involving only aerodynamics. By complementing the wind-tunnel experiments with numerical simulations the overall cost of the development of aircraft can be considerably reduced. In order to accurately compute aeroelastic phenomenon it is necessary to solve the unsteady Euler/Navier-Stokes equations simultaneously with the structural equations of motion. These equations accurately describe the flow phenomena for aeroelastic applications. At Ames a code, ENSAERO, is being developed for computing the unsteady aerodynamics and aeroelasticity of aircraft and it solves the Euler/Navier-Stokes equations. The purpose of this contract is to continue the algorithm enhancements of ENSAERO and to apply the code to complicated geometries. During the last year, the geometric capability of the code was extended to simulate transonic flows, a wing with oscillating control surface. Single-grid and zonal approaches were tested. For the zonal approach, a new interpolation technique was introduced. The key development of the algorithm was an interface treatment between moving zones for a control surface using the virtual-zone concept. The work performed during the period, 1 Apr. 1992 through 31 Mar. 1993 is summarized. Additional details on the various aspects of the study are given in the Appendices.

Obayashi, Shigeru↗

DOE OSTI2016

Improving Advanced Inverter Control Convergence in Distribution Power Flow

Simulation of modern distribution system powerflow increasingly requires capturing the impact of advanced PV inverter voltage regulation on powerflow. With Volt/var control, the inverter adjusts its reactive power flow as a function of the point of common coupling (PCC) voltage. Similarly, Volt/watt control curtails active power production as a function of PCC voltage. However, with larger systems and higher penetrations of PV, this active/reactive power flow itself can cause significant changes to the PCC voltage potentially introducing oscillations that slow the convergence of system simulations. Improper treatment of these advanced inverter functions could potentially lead to incorrect results. This paper explores a simple approach to speed such convergence by blending in the previous iteration's reactive power estimate to dampen these oscillations. Results with a single large (5MW) PV system and with multiple 500kW advanced inverters show dramatic improvements using this approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2019

Improving the Accuracy of the National Solar Radiation Database (1998-2016): Preprint

The National Solar Radiation Database (NSRDB) provides gridded satellite-based data sets for the United States and other parts of North and South America using cloud properties retrieved from the Geostationary Operational Environmental (GOES) series of satellites. These cloud properties in combination with ancillary atmospheric parameters are used in the Physical Solar Model (PSM) to compute solar radiation. This paper seeks to 1) describe the changes made in the PSM for the recent update to the National Solar Radiation Database (NSRDB 1998-2016) and 2) examine the performance of this update by comparing it to high-quality ground-based measurements. We summarize the changes to the input data for the radiative transfer models as well as updates for the processing algorithms. We also analyze the performance and accuracy of the NSRDB across the United States. Additional details of the data set and model are described in (1), (2), (3), and (4).

14 SOLAR ENERGY↗

DOE OSTI2017

Retrospective Imaging and Characterization of Nuclear Material

Modern techniques for detection of covert nuclear ma-terial requires some combination of real time measurement and/or sampling of the material. More common is real time measure-ment of the ionizing emission caused by radioactive decay or through the materials measured in response to external interroga-tion radiation. One can expose the suspect material with various radiation types, including high energy photons such as x rays or with larger particles such as neutrons and muons, to obtain images or measure nuclear reactions induced in the material. Stand-off detection using imaging modalities similar to those in the medical field can be accomplished, or simple collimated detec-tors can be used to localize radioactive materials. In all such cases, the common feature is that some or all of the nuclear materials have to be present for the measurement, which makes sense; as one might ask, “How you can measure something that is not there?” The current work and results show how to do exactly that: characterize nuclear materials after they have been removed from an area leaving no chemical trace. This new approach is demon-strated to be fully capable of providing both previous source spa-tial distribution and emission energy grouping. The technique uses magnetic resonance for organic insulators and/or lumines-cence techniques on ubiquitous refractory materials similar in theory to the way the nuclear industry carries out worker person-nel dosimetry. Spatial information is obtained by acquiring gridded samples for dosimetric measurements, while energy infor-mation comes through dose depth profile results that are func-tions of the incident radiation energies.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

NASA NTRS2016

The Viability of Trajectory Analysis for Diagnosing Dynamical and Chemical Influences on Ozone Concentrations in the UTLS

The viability of trajectory analysis for diagnosing the interplay between chemistry and dynamics is investigated by comparing ozone mixing ratios modelled using air-parcel pathways to values observed along flight tracks during ATTREX (Airborne Tropical TRopopause EXperiment). Trajectories are initiated at the locations of ozone observations and tracked backward in time to their sources at termini of backward trajectories. The modelled values of ozone utilize 3-dimensional analysis fields from WACCM (Whole Atmosphere Community Climate Model) (a chemical-climate model with dynamical fields nudged towards MERRA (Modern-Era Retrospective Analysis and Research Applications) reanalysis) and ERA-interim (product of ECMWF - the European Centre for Medium-Range Weather Forecasts) to determine source mixing ratios with chemical production and loss terms derived from the ozone chemistry used in WACCM. A statistical base of modelled ozone is constructed with 6 trajectory platforms (adiabatic, diabatic, and kinematic forced by ERA-interim and MERRA), two chemical models (WACCM chemistry and no chemistry), and 4 trajectory lengths (5, 10, 20, and 30 days). Linear regression is employed to separate systematic errors from random errors and to characterize the impact of source mixing ratios, path length, vertical motion, and chemistry on modelled ozone errors. Errors in the analysis ozone fields are large, if not dominant, contributors to model error. Random errors are particularly large for point-by-point comparisons, however averaging over 800 km (75 minutes) flight segments substantially reduces random error and exposes systematic errors. Of the two analysis ozone data sets, WACCM, which incorporates detailed chemistry, provides the smaller systematic errors while ERA-interim, which has crude chemistry but assimilates observational data, has the smaller random errors. Of the different trajectory platforms, adiabatic calculations produce the smaller random errors (irrespective of the use of chemistry) but both vertical motion and chemistry are required to optimally reduce systematic errors. These results suggest that meaningful analysis of dynamical and chemical interactions that control ozone mixing ratios are viable on spatial scales larger than a few reanalysis grid spaces, that errors in the analyzed ozone data sets are large but not prohibitively so, and that vertical velocities and heating rates from reanalysis data, while problematic, contain useful information [on the ozone concentrations in the UTLS (Upper Troposphere/Lower Stratosphere)].

ozone↗

DOE OSTI2019

Reactivity Feedback Modeling in SAM

System Analysis Module (SAM) is being developed at Argonne National Laboratory as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. It utilizes the object-oriented application framework MOOSE to leverage the modern software environment and advanced numerical methods. The capabilities of SAM are being extended to enable the transient modeling, analysis, and design of various advanced nuclear reactor systems. This report presents the development of Point-Kinetics modeling and reactivity feedback mechanisms in SAM. It also presents a new capability for predicting the thermal expansion in various structural components by coupling SAM with an external thermomechanics module (Tensor Mechanics) from the MOOSE framework. The analysis of the transient behavior of a nuclear reactor requires the coupled simulation of reactor kinetics and thermal-hydraulics of the reactor core, especially for those unprotected transients where the reactor scram system may not function properly. The point kinetics model has been widely used for reactor safety analysis due to its simplicity to capture the transient behavior of the reactor. Various reactivity feedback models have been developed and integrated with the Point-Kinetics module, including fuel axial expansion, core radial expansion, fuel Doppler, and coolant density reactivity. The reactivity feedback models in SAM are similar to the respective models used in SAS4A/SASSYS-1. This report first presents the brief theory of the Point-Kinetics module and reactivity feedback models. A number of verification tests have been performed where the code simulations are compared to the analytical model results. The reactivity feedback due to the thermal deformation, such as the fuel axial expansion and core radial expansion, is important for SFR transient analysis. Simplified thermal expansion models for the fuel pin and reactor constrain system (e.g. grid plate) are developed and verified in SAM. Additionally, a coupling interface is developed to couple SAM with external thermomechanical analysis modules for more accurate predictions of the thermal expansion of different components during the transients. The current coupling interface has been tested with the Tensor Mechanics module from MOOSE. These point kinetics and reactivity feedback modeling capabilities have also been demonstrated by simulating the early stage of the unprotected loss-of-flow (ULOF) accident in the Advanced Burner Test Reactor (ABTR). Both the stand-alone SAM and coupled SAM and Tensor Mechanics simulations are performed. It is confirmed that the major physics phenomena in the heat transport system of the ABTR reactor are captured by SAM, and the point kinetics model, reactivity feedback models, and the coupling schemes are working well as expected.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Assessing the suitability of sites near Pine Island Glacier for subglacial bedrock drilling aimed at detecting Holocene retreat–readvance

Abstract. Unambiguous identification of past episodes of ice sheet thinning below the modern surface and grounding line retreat inboard of present requires recovery and exposure dating of subglacial bedrock. Such efforts are needed to understand the significance and potential future reversibility of ongoing and projected change in Antarctica. Here we evaluate the suitability for subglacial bedrock drilling of sites in the Hudson Mountains, which are located in the Amundsen Sea sector of West Antarctica. We use an ice sheet model and field data – geological observations, glaciological observations and bedrock samples from nunataks, and ground-penetrating radar from subglacial ridges – to rate each site against four key criteria: (i) presence of ridges extending below the ice sheet, (ii) likelihood of increased exposure of those ridges if the grounding line was inboard of present, (iii) suitability of bedrock for drilling and geochemical analysis, and (iv) accessibility for aircraft and drilling operations. Our results demonstrate that although no site in the Hudson Mountains is perfect for this study when assessed against all criteria, the accessibility, N–S orientation and basaltic bedrock lithology of Winkie Nunatak's southernmost ridge (74.86° S, 99.77° W) make it a feasible site both for drilling and subsequent cosmogenic nuclide analysis. Furthermore, the ridge is strewn with glacial erratics at all elevations, providing valuable constraints on its early Holocene deglacial history. Based on our experiences during this study, we conclude with a series of recommendations for assessing site suitability for future bedrock drilling campaigns. We emphasise the importance of consulting a range of expertise prior to drilling and ensuring that sufficient field reconnaissance is undertaken (including obtaining detailed grids of radar survey data and bedrock samples).

58 GEOSCIENCES↗

DOE OSTI2012

Radiation modeling using the Uintah heterogeneous CPU/GPU runtime system. In: XSEDE '12 Proceedings of the 1st Conference of the Extreme Science and Engineering Discovery Environment: Bridging from the eXtreme to the campus and beyond, Article No. 4

The Uintah Computational Framework was developed to provide an environment for solving fluid-structure interaction problems on structured adaptive grids on large-scale, long-running, data-intensive problems. Uintah uses a combination of fluid-flow solvers and particle-based methods for solids, together with a novel asynchronou task-based approach with fully automated load balancing. Uintah demonstrates excellent weak and strong scalability at full machine capacity on XSEDE resources such as Ranger and Kraken, and through the use of a hybrid memory approach based on a combination of MPI and Pthreads, Uintah now runs on up to 262k cores on the DOE Jaguar system. In order to extend Uintah to heterogeneous systems, with ever-increasing CPU core counts and additional on-node GPUs, a new dynamic CPU-GPU task scheduler is designed and evaluated in this study. This new scheduler enables Uintah to fully exploit these architectures with support for asynchronous, out-of-order scheduling of both CPU and GPU computational tasks. A new runtime system has also been implemented with an added multi-stage queuing architecture for efficient scheduling of CPU and GPU tasks. This new runtime system automatically handles the details of asynchronous memory copies to and from the GPU and introduces a novel method of pre-fetching and preparing GPU memory prior to GPU task execution. In this study this new design is examined in the context of a developing, hierarchical GPU-based ray tracing radiation transport model that provides Uintah with additional capabilities for heat transfer and electromagnetic wave propagation. The capabilities of this new scheduler design are tested by running at large scale on the modern heterogeneous systems, Keeneland and TitanDev, with up to 360 and 960 GPUs respectively. On these systems, we demonstrate significant speedups per GPU against a standard CPU core for our radiation problem.

Humphrey, Alan↗

DOE OSTI2023

Considerations for Fire Service Response to Residential Battery Energy Storage System Incidents

Renewable energy sources (e.g., rooftop photovoltaics, wind turbines) are capable of addressing our grid-level energy challenges by reducing environmental impacts, increasing resiliency, and increasing the supply of input energy. However, renewable energy sources alone are not a complete solution. The generation of power from major renewable energy sources fluctuates with the weather, creating significant challenges in matching electrical power generation to electrical power consumption. The value of renewable power generation sources is multiplied when paired with battery energy storage systems (ESS), which can store excess power generated from renewable energy sources when it is not consumed and deliver that power later when demand exceeds base power generation. Residential ESS units are frequently installed to support renewable energy initiatives by storing energy from intermittent sources such as photovoltaic panels. Residential ESS units also offer an alternative to gas-powered generators or other backup power means. Lithium-ion batteries are the most common residential ESS technology due to their affordability and energy density. Though lithium-ion systems come with benefits, there is also a risk of thermal runaway within this technology which can result in flammable gas release, fire, and explosion. As the installations of these products increase, the frequency of response to fire incidents involving these products will increase. In response to this new and evolving hazard, UL Solutions (ULS) and UL Fire Safety Research Institute (FSRI) have partnered with the International Association of Fire Fighters (IAFF) to conduct a series of large-scale tests sponsored by the US Department of Energy to characterize the challenges for fire fighters responding to fires involving residential energy storage systems. The project focuses on developing size-up and tactical considerations to support the fire service in navigating the evolving modern fireground.

25 ENERGY STORAGE↗

NASA NTRSDate not supplied

High Performance Access to Archival Data Stored in HDF4 and HDF5 on Cloud Object Stores Without Reformatting the Files

Cloud computing offers numerous advantages for users of extensive Earth science data collections. These benefits encompass direct online access to data files and granules from any location, scalable access supporting parallel computing workflows, and flexible computing tools enabling innovative experimentation with processing techniques. However, older archival file formats designed for distinct computing systems hinder efficient access to decade-long time-series data when compared to data stored in modern cloud-optimized formats like Web Object Stores (WOS), exemplified by Amazon Web Services’ Simple Storage Service (S3). We describe DMR++ (Dataset Metadata Response plus plus), a technology facilitating efficient access to HDF5 (Hierarchical Data Format, version 5) and HDF4 files stored on WOS systems without requiring data reformatting. DMR++ achieves performance comparable to technologies like Zarr while preserving the original file structure, a substantial benefit considering the vast quantity of archival files held by organizations such as NASA. Moreover, DMR++ typically outperforms cloud-optimized versions of HDF5. Essentially an XML (Extensible Markup Language) document usually stored alongside the described data, DMR++ can also be generated on-the-fly but is generally created during data staging to the WOS. Archival files that use HDF4/5 often store large arrays of numerical data. The data in these files is often compressed, typically reducing their size by a factor of four or more. To achieve efficient access to portions of those arrays, they are 'chunked' into smaller sub-arrays, each individually compressed. The chunk size is a compromise, where spinning disks can efficiently access data in smaller chunks while S3 favors larger chunks. A simple optimization of aggregating smaller chunks that are stored adjacently, transferring them in a single access and then individually decompressing them will improve performance. NASA data pose an additional challenge: special Application Programmer Interface (API) libraries are often needed to compute some variables. These libraries are incompatible with WOS environments. Our solution involves storing computed values in the DMR++ document or a companion file, making them accessible like other variables and eliminating the need for specialized APIs. We outline specific optimizations for both satellite grid and swath data stored in HDF4-EOS2 (Earth Observing System).

James Gallagher↗

DOE OSTI2018

Detecting False Data Injection Attacks Against Power System State Estimation With Fast Go-Decomposition Approach

State estimation is a fundamental function in modern energy management system, but its results may be vulnerable to false data injection attacks (FDIAs). FDIA is able to change the estimation results without being detected by the traditional bad data detection algorithms. In this paper, we propose an accurate and computational attractive approach for FDIA detection. In this work, we first rely on the low rank characteristic of the measurement matrix and the sparsity of the attack matrix to reformulate the FDIA detection as a matrix separation problem. Then, four algorithms that solve this problem are presented and compared, including the traditional augmented Lagrange multipliers (ALMs), double-noise-dual-problem (DNDP) ALM, the low rank matrix factorization, and the proposed new “Go Decomposition (GoDec).” Numerical simulation results show that our GoDec algorithm outperforms the other three alternatives and demonstrates a much higher computational efficiency. Furthermore, GoDec is shown to be able to handle measurement noise and applicable for large-scale attacks.

24 POWER TRANSMISSION AND DISTRIBUTION↗