Search NASA⌕ Search

SEARCH · Search NASA

Results for “Software dependability”

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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Energy savings modelling of re-tuning energy conservation measures in large office buildings

Today, many large commercial buildings use sophisticated building automation systems (BASs) to manage a wide range of building equipment. While the capabilities of BASs have increased over time, many buildings still do not fully use the BAS’s capabilities and are not properly commissioned, operated or maintained, which leads to inefficient operation, increased energy use, and reduced lifetimes of the equipment. This paper investigates the energy savings potential of several common HVAC system re-tuning measures on a typical large office building, using the Department of Energy’s building energy modeling software, EnergyPlus. The baseline prototype model uses roughly as much energy as an average large office building in existing building stock, but does not utilize any re-tuning measures. Individual re-tuning measures simulated against this baseline include automatic schedule adjustments, damper minimum flow adjustments, thermostat adjustments, as well as dynamic resets (set points that change continuously with building and/or outdoor conditions) to static pressure, supply-air temperature, condenser water temperature, chilled and hot water temperature, and chilled and hot water differential pressure set points. Six combinations of these individual measures have been formulated – each designed to conform to limitations to implementation of certain individual measures that might exist in typical buildings. All the individual measures and combinations were simulated in 16 climate locations representative of specific U.S. climate zones. The modeling results suggest that the most effective energy savings measures are those that affect the demand-side of the building (air-systems and schedules). Many of the demand-side individual measures were capable of reducing annual total HVAC system energy consumption by over 20% in most cities that were modeled. Supply side measures affecting HVAC plant conditions were only modestly successful (less than 5% annual HVAC energy savings for most cities for all measures). Combining many of the re-tuning measures revealed deep savings potential. Some of the more aggressive combinations revealed 35-75% reductions in annual HVAC energy consumption, depending on climate and building vintage.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2018

Mechanism of Low-Temperature Reduction Degradation of Alumina-Containing Hematite Solid Solution Below 550 °C

Low-temperature reduction degradation (LTRD) of sinter has an adverse effect on blast furnace permeability, and it is mainly caused by the stress produced in the reduction process of hematite. This stress is strongly influenced by alumina dissolved in hematite crystal lattice. In this work, the experiments were conducted to investigate the effect of alumina dissolved in hematite solid solution (H{sub ss}) on LTRD by reducing H{sub ss} below 550 °C. Thermogravimetric analysis (TGA), X-ray diffraction (XRD), and optical microscope have been used to characterize the mass change and mineral change of samples, respectively. Jade software has been used to calculate the micro strain in magnetite for quantitatively studying the change of strain in reducing process. The results show that alumina was unfavorable to the reduction of H{sub ss} on thermodynamics, and the starting reduction temperature of H{sub ss} containing 6.0 mol pct alumina was 28 °C higher than that of pure hematite. According to the calculation on kinetics, the generation rate of stress was accelerated by dissolving alumina into hematite crystal lattice. The apparent activation energy of reduction reaction lowered from 47.89 to 28.07 kJ/mol with the increase of alumina content from 0.0 to 6.0 mol pct. The addition of alumina also increased the stress in the reduction products, and this stress was released in the form of LTRD.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Improvements to the Ion Doppler Spectrometer Diagnostic on the HIT-SI Experiments

An Ion Doppler Spectrometer diagnostic system measuring impurity ion temperature and velocity on the HIT-SI and HIT-SI3 spheromak devices has been improved with higher spatiotemporal resolution and lower error than previously described devices. Hardware and software improvements to the established technique have resulted in a record 6.9 µs temporal and <=2.8 cm spatial resolution in the midplane of each device. These allow C III and O II flow, displacement, and temperature profiles to be simultaneously observed. With 72 fused-silica fiber channels in two independent bundles, and an f/8.5 Czerny-Turner spectrometer coupled to video camera, frame-rates of up to ten times the imposed magnetic perturbation frequency of 14.5 kHz were achieved in HIT-SI, viewing the upper 1/2 of the midplane. In HIT-SI3 frame-rates of up to eight times the perturbation frequency were achieved viewing both halves of the midplane. Biorthogonal Decomposition is used as a novel filtering tool, reducing uncertainty in ion temperature from <=13 to <=5 eV (with an instrument temperature of 8-16 eV), and uncertainty in velocity from <=2 to <=1 km/s. Doppler shift and broadening is calculated via the Levenberg-Marquart algorithm, after which errors in velocity and temperature are uniquely specified. Axisymmetric temperature profiles on HIT-SI3 for C III peaked near the inboard current separatrix at approximately 40 eV are observed. Axisymmetric plasma displacement profiles have been measured on HIT-SI3, peaking at approximately 6 cm at the outboard separatrix. Both profiles agree with the upper half of the midplane observable by HIT-SI. With its complete midplane view, HIT-SI3 has unambiguously extracted axisymmetric, toroidal current dependent rotation of up to 3 km/s. Analysis of the temporal phase of the displacement uncovers a coherent structure, locked to the applied perturbation. Previously described diagnostic systems could not achieve such results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Traffic Control via Connected and Automated Vehicles (CAVs): An Open-Road Field Experiment with 100 CAVs

The CIRCLES project aims to reduce instabilities in traffic flow, which are naturally occurring phenomena due to human driving behavior. Also called “phantom jams” or “stop-and-go waves,” these instabilities are a significant source of wasted energy. Toward this goal, the CIRCLES project designed a control system, referred to as the MegaController by the CIRCLES team, that could be deployed in real traffic. Our field experiment, the MegaVanderTest (MVT), leveraged a heterogeneous fleet of 100 longitudinally controlled vehicles as Lagrangian traffic actuators, each of which ran a controller with the architecture described in this article. The MegaController is a hierarchical control architecture that consists of two main layers. The upper layer is called the Speed Planner and is a centralized optimal control algorithm. It assigns speed targets to the vehicles, conveyed through the LTE cellular network. The lower layer is a control layer, running on each vehicle. It performs local actuation by overriding the stock adaptive cruise controller, using the stock onboard sensors. The Speed Planner ingests live data feeds provided by third parties as well as data from our own control vehicles and uses both to perform the speed assignment. The architecture of the Speed Planner allows for the modular use of standard control techniques, such as optimal control, model predictive control (MPC), kernel methods, and others. The architecture of the local controller allows for the flexible implementation of local controllers. Corresponding techniques include deep reinforcement learning (RL), MPC, and explicit controllers. Depending on the vehicle architecture, all onboard sensing data can be accessed by the local controllers or only some. Likewise, control inputs vary across different automakers, with inputs ranging from torque or acceleration requests for some cars to electronic selection of adaptive cruise control (ACC) setpoints in others. The proposed architecture technically allows for the combination of all possible settings proposed previously, that is {Speed Planner algorithms} × {local Vehicle Controller algorithms} × {full or partial sensing} × {torque or speed control}. As a result, most configurations were tested throughout the ramp up to the MegaVandertest (MVT).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2018

Compensation of Axis Errors of Azimuth and Zenith Moving Concentrators in Programmable Solar-Tracking Systems

The solar-energy industry is currently concerned with improving software systems of sun-tracking concentrators. It is a feature of programmable tracking that there is no feedback from the tracked object, unlike in the case of optical tracking. In this regard, the accuracy of programmable tracking is affected by a large number of factors—in particular, axis-alignment errors. An error-checking algorithm is considered for the azimuth and zenith concentrator axes in programmable sun tracking. Effects of these errors on the accuracy of programmable tracking are assessed. It is shown that zenith axis errors have the same amount of influence as the vertical deviation of the azimuth axis, primarily affecting the azimuth angle, while its influence on the solar elevation can be neglected. The cumulative effect of axial errors is most noticeable when they have the same signs and lie in the same plane.

14 SOLAR ENERGY↗
DOE OSTI2018

Improving signal strength in serial crystallography with DIALS geometry refinement

The DIALS diffraction-modeling software package has been applied to serial crystallography data. Diffraction modeling is an exercise in determining the experimental parameters, such as incident beam wavelength, crystal unit cell and orientation, and detector geometry, that are most consistent with the observed positions of Bragg spots. These parameters can be refined by nonlinear least-squares fitting. In previous work, it has been challenging to refine both the positions of the sensors (metrology) on multipanel imaging detectors such as the CSPAD and the orientations of all of the crystals studied. Since the optimal models for metrology and crystal orientation are interdependent, alternate cycles of panel refinement and crystal refinement have been required. To simplify the process, a sparse linear algebra technique for solving the normal equations was implemented, allowing the detector panels to be refined simultaneously against the diffraction from thousands of crystals with excellent computational performance. Separately, it is shown how to refine the metrology of a second CSPAD detector, positioned at a distance of 2.5 m from the crystal, used for recording low-angle reflections. With the ability to jointly refine the detector position against the ensemble of all crystals used for structure determination, it is shown that ensemble refinement greatly reduces the apparent nonisomorphism that is often observed in the unit-cell distributions from still-shot serial crystallography. In addition, it is shown that batching the images by timestamp and re-refining the detector position can realistically model small, time-dependent variations in detector position relative to the sample, and thereby improve the integrated structure-factor intensity signal and heavy-atom anomalous peak heights.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Device Scale Modeling of Solvent Absorption using MFIX-TFM

Recent climate change is largely attributed to greenhouse gases (e.g., carbon dioxide, methane) and fossil fuels account for a large majority of global CO 2 emissions. That said, fossil fuels will continue to play a significant role in the generation of power for the foreseeable future. The extent to which CO 2 is emitted needs to be reduced, however, carbon capture and sequestration are also necessary actions to tackle climate change. Different approaches exist for CO 2 capture including both post-combustion and pre-combustion technologies, oxy-fuel combustion and/or chemical looping combustion. The focus of this effort is on post-combustion solvent-absorption technology. To apply CO 2 technologies at commercial scale, the availability and maturity and the potential for scalability of that technology need to be considered. Solvent absorption is a proven technology but not at the scale needed by typical power plant. The scale up and down and design of laboratory and commercial packed bed reactors depends heavily on the specific knowledge of two-phase pressure drop, liquid holdup, the wetting efficiency and mass transfer efficiency as a function of operating conditions. Simple scaling rules often fail to provide proper design. Conventional reactor design modeling approaches will generally characterize complex non-ideal flow and mixing patterns using simplified and/or mechanistic flow assumptions. While there are varying levels of complexity used within these approaches, none of these models resolve the local velocity fields. Consequently, they are unable to account for important design factors such as flow maldistribution and channeling from a fundamental perspective. Ideally design would be aided by development of predictive models based on truer representation of the physical and chemical processes that occur at different scales. Computational fluid dynamic (CFD) models are based on multidimensional flow equations with first principle foundations. CFD models can include a more accurate physical description of flow processes and be modified to include more complex behavior. Wetting performance and spatial liquid distribution inside the absorber are recognized as weak areas of knowledge requiring further investigation. CFD tools offer a possible method to investigating such topics and gaining a better understanding of their influence on reactor performance. This report focuses first on describing a hydrodynamic model for countercurrent gas-liquid flow through a packed column and then on the chemistry, heat and mass transfer specific to CO 2 absorption using monoethanolamine (MEA). The indicated model is implemented in MFIX, a CFD open source software package. The user defined functions needed to build this model are described in detail along with the keywords for the corresponding input file. A test case is outlined along with a few results. The example serves to briefly illustrate the developed CFD tool and its potential capability to investigate solvent absorption.

09 BIOMASS FUELS↗
DOE OSTI2025

A Climatology and Life‐Cycle Characteristics of Atmospheric Fronts and Their Associated Precipitation

Abstract Atmospheric fronts are one of the main sources of mid‐latitude variability. We employ a novel method for identifying and tracking fronts and frontal precipitation. Thermal and dynamical variables are used to identify fronts as areal objects in space, which are tracked in time using the open‐source TempestExtremes software package. Precipitation objects are co‐located to identify frontal precipitation. The method is subjected to validation and sensitivity tests using manually curated data from the National Weather Service. Climatologies of fronts and frontal precipitation are computed from reanalysis and observations; fronts are present upwards of 14% of the time in the storm tracks, and represent the majority (up to 90%) of total and extreme precipitation. Novel aspects of the method are showcased through the lifetime characteristics of fronts across North America. Three sets of warm and cold fronts were discovered, and their duration, distance‐traveled, and translation velocity are examined. Plain Language Summary Mid‐latitude low‐pressure systems and weather fronts are important for our day‐to‐day experience of weather events, particularly in the mid‐latitudes. This work makes use of standardized atmospheric data and creates a method of automatically tracking these important atmospheric features and their precipitation to quantify their relative role in global precipitation. Weather fronts are persistent in the mid‐latitudes and are associated with the majority of precipitation–particularly the most intense precipitation. Trajectories of fronts over North America are categorized to create a set of archetypal fronts that occur in that region. The differences between these types of fronts are characterized. Key Points An automated, efficient, and skillful frontal detection algorithm is developed and validated Fronts contribute a larger fraction of extreme precipitation than all precipitation in mid‐latitude storm tracks Fronts across North America have substantial variation in characteristics depending on their origin location

extratropical cyclone↗
DOE OSTI2019

Digital Time-Resolved Spot Diagnostic

The digital time-resolved spot diagnostic (D-TRSD) is intended to be a fast-response, time-dependent beam diagnostic for measuring the spot size emitted from radiographic x-ray sources. It will measure spot sizes from submillimeters up to several millimeters for photon energies in the range of 150 keV to 20 MeV with doses as low as 0.03 rad up to 500 rad. The D-TRSD performs fast data acquisition as radiographic x-ray sources typically produce pulses in the tens of nanosecond ranges. The system will be optimized to capture multiple individual pulses in time within the same data record. The data collected will be used to support beam tuning and source optimization with on-the-fly evaluation of radiographic source performance. The diagnostic is based on penumbral imaging techniques. A scintillating fiber array of parallel aligned fibers is positioned in the beam path with a high-Z rolled edge inserted between the source and the array, thereby creating the penumbral x-ray signal that is captured by the D-TRSD system. The light from the scintillation conversion is recorded by individual avalanche photodiode detectors integrated into custom digitizer electronics. The penumbra is used to calculate the spot of the beam on the axis orthogonal to the fiber array based on the magnification from the rolled edge, the resolution of the spaced fibers, and the rotation of the array normal to the beam path. Many arrays can be used to capture additional spot widths about different rotations. The D-TRSD system is modular, software-driven, and easier to use than legacy TRSD systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2013

X-ray phase contrast imaging of the breast: Analysis of tissue simulating materials

Purpose: Phase contrast imaging, particularly of the breast, is being actively investigated. The purpose of this work is to investigate the x-ray phase contrast properties of breast tissues and commonly used breast tissue substitutes or phantom materials with an aim of determining the phantom materials best representative of breast tissues. Methods: Elemental compositions of breast tissues including adipose, fibroglandular, and skin were used to determine the refractive index, n= 1 -{delta}+i {beta}. The real part of the refractive index, specifically the refractive index decrement ({delta}), over the energy range of 5-50 keV were determined using XOP software (version 2.3, European Synchrotron Radiation Facility, France). Calcium oxalate and calcium hydroxyapatite were considered to represent the material compositions of microcalcifications in vivo. Nineteen tissue substitutes were considered as possible candidates to represent adipose tissue, fibroglandular tissue and skin, and four phantom materials were considered as possible candidates to represent microcalcifications. For each material, either the molecular formula, if available, or the elemental composition based on weight fraction, was used to determine {delta}. At each x-ray photon energy, the absolute percent difference in {delta} between the breast tissue and the substitute material was determined, from which three candidates were selected. From these candidate tissue substitutes, the material that minimized the absolute percent difference in linear attenuation coefficient {mu}, and hence {beta}, was considered to be best representative of that breast tissue. Results: Over the energy range of 5-50 keV, while the {delta} of CB3 and fibroglandular tissue-equivalent material were within 1% of that of fibroglandular tissue, the {mu} of fibroglandular tissue-equivalent material better approximated the fibroglandular tissue. While the {delta} of BR10 and adipose tissue-equivalent material were within 1% of that of adipose tissue, the tissue-equivalent material better approximated the adipose tissue in terms of {mu}. Polymethyl methacrylate, a commonly used tissue substitute, exhibited {delta} greater than fibroglandular tissue by {approx}12%. The A-150 plastic closely approximated the skin. Several materials exhibited {delta} between that of adipose and fibroglandular tissue. However, there was an energy-dependent mismatch in terms of equivalent fibroglandular weight fraction between {delta} and {mu} for these materials. For microcalcifications, aluminum and calcium carbonate were observed to straddle the {delta} and {mu} of calcium oxalate and calcium hydroxyapatite. Aluminum oxide, commonly used to represent microcalcifications in the American College of Radiology recommended phantoms for accreditation exhibited {delta} greater than calcium hydroxyapatite by {approx}23%. Conclusions: A breast phantom comprising A-150 plastic to represent the skin, commercially available adipose and fibroglandular tissue-equivalent formulations to represent adipose and fibroglandular tissue, respectively, was found to be best suited for x-ray phase-sensitive imaging of the breast. Calcium carbonate or aluminum can be used to represent microcalcifications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Weatherization Assistant NEAT/MHEA

The software provides a measure selection technique indicating cost effective retrofit activities that can be applied to a home using a standard Savings to Investment Ratio (SIR). Users must provide an input file describing the characteristics of the home to be evaluated. The software takes the input data provided and calculates energy savings and cost savings predicted for a standard set of measures given the input parameters. The Weatherization Assistant computes estimates of pre-retrofit whole building space heating and cooling energy consumptions based on the house description data supplied by the user. The consumptions are computed using a monthly heating and cooling variable base degree-day method by algorithms similar to those developed for the CIRA program [LBL, 1982]. The building consumptions are needed in computing the energy savings from measures affecting the efficiencies of the heating and cooling equipment. Weatherization Assistant then computes the energy savings and costs for each individual measure applicable to the building described as if it were the only measure installed in the house. From these energy savings, a discounted dollar savings over the life of each measure is computed. The ratio of this dollar savings to the cost of installing the measure, the "savings-to-investment ratio" (SIR), is used in an initial ranking of the measures' effectiveness. The "interacted" savings and SIR of measures are then determined assuming the measures are added to the house collectively, in order of their ranking, e.g., the second ranked measure is installed in the house initially described by the user after having been modified by the first ranked measure. If this second-ranked measure's updated SIR is greater than a user-defined limit, the measure is left implemented, else it is removed so that the next measure's effectiveness is not dependent on it. The choice between two mutually exclusive measures (such as different levels of insulation) is made on the basis of their "net present value" (NPV), the difference of life-time savings and installation cost, rather than their SIR. This has been shown to be the more correct criterion on which to base the selection between two measures, both of which cannot be installed. The audit computes and reports to the user the energy savings, discounted dollar savings, installation cost, and SIR for each measure considered cost-effective. For those with SIR greater than the user-designated cutoff, a materials list gives the material name, type, and quantity required for installation of the measure. Weatherization Assistant permits entry of pre-retrofit billing data for gas or electrically heated homes or homes with electric air-conditioning. The user may then make the decision to have the savings of the measures adjusted to reflect the difference in billed consumption and that predicted by the program.

Gettings, Michael↗
DOE OSTI2016

Evaluation of various approaches for assessing dose indicators and patient organ doses resulting from radiotherapy cone-beam CT

Purpose: The aim of this study was to evaluate various approaches for assessing patient organ doses resulting from radiotherapy cone-beam CT (CBCT), by the use of thermoluminescent dosimeter (TLD) measurements in anthropomorphic phantoms, a Monte Carlo based dose calculation software, and different dose indicators as presently defined. Methods: Dose evaluations were performed on a CBCT Elekta XVI (Elekta, Crawley, UK) for different protocols and anatomical regions. The first part of the study focuses on using PCXMC software (PCXMC 2.0, STUK, Helsinki, Finland) for calculating organ doses, adapting the input parameters to simulate the exposure geometry, and beam dose distribution in an appropriate way. The calculated doses were compared to readouts of TLDs placed in an anthropomorphic Rando phantom. After this validation, the software was used for analyzing organ dose variability associated with patients’ differences in size and gender. At the same time, various dose indicators were evaluated: kerma area product (KAP), cumulative air-kerma at the isocenter (K{sub air}), cone-beam dose index, and central cumulative dose. The latter was evaluated in a single phantom and in a stack of three adjacent computed tomography dose index phantoms. Based on the different dose indicators, a set of coefficients was calculated to estimate organ doses for a range of patient morphologies, using their equivalent diameters. Results: Maximum organ doses were about 1 mGy for head and neck and 25 mGy for chest and pelvis protocols. The differences between PCXMC and TLDs doses were generally below 10% for organs within the field of view and approximately 15% for organs at the boundaries of the radiation beam. When considering patient size and gender variability, differences in organ doses up to 40% were observed especially in the pelvic region; for the organs in the thorax, the maximum differences ranged between 20% and 30%. Phantom dose indexes provided better correlation with organ doses than K{sub air} and KAP, with average ratios ranging between 0.9 and 1.1 and variations for different organs and protocols below 20%. The triple phantom setup allowed us to take into account scatter dose contributions, but nonetheless, the correlation with the evaluated organ doses was not improved with this method. Conclusions: The simulation of rotational geometry and of asymmetric beam distribution by means of PCXMC 2.0 enabled us to determine patient organ doses depending on weight, height and gender. Alternatively, the measurement of an in phantom dose indicator combined with proper correction coefficients can be a useful tool for a first dose estimation of in-field organs. The data and coefficients provided in this study can be applied to any patient undergoing a scan by an Elekta XVI equipment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Workshop on Improving Holdup Monitoring in the US

The characterization and quantification of nuclear deposits, or residual nuclear material retained in process equipment, generally referred to as holdup, continues to challenge nuclear processing facilities in DOE program areas such as nuclear criticality safety, material control and accountability, environmental management. The efficiency and effectiveness of nuclear operations in these facilities depends heavily on the results from nondestructive techniques designed to measure nuclear materials in situ. The traditional methods used to measure holdup were established more than 30 years ago and assume unrealistic conditions for most measurement cases. These blind, in situ measurements present unique challenges due to the following attributes: unique, non-ideal geometries; unknown deposit thicknesses; lack of representative calibration standards; large number of holdup deposit locations; and poor accessibility for many measurement locations. These difficulties and the poor assumptions made to handle them have resulted in measurement uncertainties that are often too large for many programs to effectively utilize. In an effort to address and identify solutions to the current holdup challenges facing the US nuclear industry, Oak Ridge National Laboratory hosted the Technical Workshop on Improving Holdup Monitoring in the US. The two-day workshop held August 21–22, 2019, was attended by 35 holdup measurement practitioners, program managers, stakeholders from US Department of Energy sites, the Nuclear Regulatory Commission, and representatives from commercial industry. The primary objectives were to identify critical elements of a successful and defensible holdup measurement program; identify technical challenges, needs, and potential improvements associated with measurement of holdup; and share best practices and lessons learned from recent accumulation events. Fourteen presentations were delivered on holdup topics including necessary elements for effective program management, measurement needs, emerging technologies and other measurement solutions, and lessons learned from recent events. These presentations highlighted best practices related to accumulation monitoring, considerations and proactive steps to facilitate holdup in new facility designs, imaging solutions to improve holdup assumptions, and new software developed in support of holdup measurements. Breakout sessions were also employed to determine additional needs and identify potential solutions to common holdup challenges. The primary discussions focused on causes and prevention of material accumulations, needs for uranium and plutonium measurements, and reduction of measurement uncertainty, all of which were explicitly identified as critical needs during the 2018 Workshop on Technical and Programmatic Needs for a Sustainable NDA Program for the US DOE.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Two-dimensional simulations of internal gravity waves in a $\mathrm{5M_⊙}$ zero-age-main-sequence model

Main-sequence intermediate-mass stars present a radiative envelope that supports internal gravity waves (IGWs). Excited at the boundary with the convective core, IGWs propagate towards the stellar surface and are suspected to impact physical processes such as rotation and chemical mixing. Using the fully compressible time-implicit code MUSIC, we study IGWs in two-dimensional simulations of a zero-age-main-sequence 5 solar mass star model up to 91 percent of the stellar radius with different luminosity and radiative diffusivity enhancements. Our results show that low-frequency waves excited by core convection are strongly impacted by radiative effects as they propagate. This impact depends on the radial profile of radiative diffusivity which increases by almost 5 orders of magnitude between the centre of the star and the top of the simulation domain. In the upper layers of the simulation domain, we observe an increase of the temperature. Our study suggests that this is due to heat added in these layers by IGWs damped by radiative diffusion. We show that non-linear effects linked to large amplitude IGWs may be relevant just above the convective core. Both these effects are intensified by the artificial enhancement of the luminosity and radiative diffusivity, with enhancement factors up to 10 4 times the realistic values. Our results also highlight that direct comparison between numerical simulations with enhanced luminosity and observations must be made with caution. Finally, our work suggests that thermal effects linked to the damping of IGWs could have a non-negligible impact on stellar structure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

waveSZ: A Hardware-Algorithm Co-Design of Efficient Lossy Compression for Scientific Data

Error-bounded lossy compression is critical to the success of extreme-scale scientific research because of ever-increasing volumes of data produced by today's high-performance computing (HPC) applications. Not only can error-controlled lossy compressors significantly reduce the I/O and storage burden but they can retain high data fidelity for post analysis. Existing state-of-the-art lossy compressors, however, generally suffer from relatively low compression and decompression throughput (up to hundreds of megabytes per second on a single CPU core), which considerably restrict the adoption of lossy compression by many HPC applications especially those with a fairly high data production rate. In this paper, we propose a highly efficient lossy compression approach based on field programmable gate arrays (FPGAs) under the state-of-the-art lossy compression model SZ. Our contributions are fourfold. (1) We adopt a wavefront memory layout to alleviate the data dependency during the prediction for higher-dimensional predictors, such as the Lorenzo predictor. (2) We propose a co-design framework named WAVESZ based on the wavefront memory layout and the characteristics of SZ algorithm and carefully implement it by using high-level synthesis. (3) We propose a hardware-algorithm co-optimization method to improve the performance. (4) We evaluate our proposed WAVESZ on three real-world HPC simulation datasets from the Scientific Data Reduction Benchmarks and compare it with other state-of-the-art methods on both CPUs and FPGAs. Experiments show that our WAVESZ can improve SZ's compression throughput by 6.9x similar to 8.7x over the production version running on a state-of-the-art CPU and improve the compression ratio and throughput by 2.1x and 5.8x on average, respectively, respectively, compared with the state-of-the-art FPGA design.

Compression Ratio↗
DOE OSTI2021

Microgrid's Role in Enhancing the Security and Flexibility of City Energy Systems

Smart cities depend on flexible and secure energy systems to ensure resilient power for critical infrastructure; however, recent weather-related events and cyberattacks have highlighted weaknesses in our energy systems, with the potential for widespread economic and security impacts. As stated by the Executive Office of the President, "the resilience of the US electric grid is a key part of the nation's defense against severe weather." To address the energy delivery security challenge, microgrids are rising as a viable solution that enhances the flexibility and resilience of the distribution grid and boosts the reliability of the local supply for the end-user. Traditionally, high capital investment has been a barrier to large-scale adoption of microgrid technology. Understanding the flexibility and resilience benefits of microgrids and accounting for the associated value streams can make the microgrid's proposition economically viable. In this chapter, microgrids' utility and their potential to serve as a flexible and resilient resource for the utility grid by providing capabilities such as peak shaving, demand response, and frequency regulation is presented. Moreover, other value streams, such as (1) their ability to island during a disaster and sustain critical loads which makes them a robust resilience solution for end-users, in the event of the utility grid outage and (2) microgrids also provide a flexible platform for integrating distributed energy resources in conjunction with storage and conventional generation technologies, strengthen microgrid's role in reducing the over-arching goal of emission reduction. Given the myriad of benefits associated with microgrids, we present strategies which can be employed for making microgrid itself resilient against physical and cyberthreats by employing hardware, software, and personnel training solutions to operate the microgrid before, during, and after a potential disaster. This chapter, thus, provides a holistic study of the microgrid as a resilience resource, for the utility grid, and a self-contained end-user for the end-user.

cyber-physical system↗
DOE OSTI2024

VAN-DAMME: GPU-accelerated and symmetry-assisted quantum optimal control of multi-qubit systems

We present an open-source software package, VAN-DAMME (Versatile Approaches to Numerically Design, Accelerate, and Manipulate Magnetic Excitations), for massively-parallelized quantum optimal control (QOC) calculations of multi-qubit systems. To enable large QOC calculations, the VAN-DAMME software package utilizes symmetry-based techniques with custom GPU-enhanced algorithms. This combined approach allows for the simultaneous computation of hundreds of matrix exponential propagators that efficiently leverage the intra-GPU parallelism found in high-performance GPUs. In addition, to maximize the computational efficiency of the VAN-DAMME code, we carried out several extensive tests on data layout, computational complexity, memory requirements, and performance. These extensive analyses allowed us to develop computationally efficient approaches for evaluating complex-valued matrix exponential propagators based on Padé approximants. To assess the computational performance of our GPU-accelerated VAN-DAMME code, we carried out QOC calculations of systems containing 10 - 15 qubits, which showed that our GPU implementation is 18.4× faster than the corresponding CPU implementation. Our GPU-accelerated enhancements allow efficient calculations of multi-qubit systems, which can be used for the efficient implementation of QOC applications across multiple domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Investigating resource-efficient neutron/gamma classification ML models targeting eFPGAs

There has been considerable interest and resulting progress in implementing machine learning (ML) models in hardware over the last several years from the particle and nuclear physics communities. A big driver has been the release of the Python package, hls4ml, which has enabled porting models specified and trained using Python ML libraries to register transfer level (RTL) code. So far, the primary end targets have been commercial field-programmable gate arrays (FPGAs) or synthesized custom blocks on application specific integrated circuits (ASICs). However, recent developments in open-source embedded FPGA (eFPGA) frameworks now provide an alternate, more flexible pathway for implementing ML models in hardware. These customized eFPGA fabrics can be integrated as part of an overall chip design. In general, the decision between a fully custom, eFPGA, or commercial FPGA ML implementation will depend on the details of the end-use application. In this work, we explored the parameter space for eFPGA implementations of fully-connected neural network (fcNN) and boosted decision tree (BDT) models using the task of neutron/gamma classification with a specific focus on resource efficiency. We used data collected using an AmBe sealed source incident on Stilbene, which was optically coupled to an OnSemi J-series silicon photomultiplier (SiPM) to generate training and test data for this study. We investigated relevant input features and the effects of bit-resolution and sampling rate as well as trade-offs in hyperparameters for both ML architectures while tracking total resource usage. The performance metric used to track model performance was the calculated neutron efficiency at a gamma leakage of 10 -3 . The results of the study will be used to aid the specification of an eFPGA fabric, which will be integrated as part of a test chip.

47 OTHER INSTRUMENTATION↗