Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “temporal computing”

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 541 records · Page 30

Analytical Models of Frequency and Voltage in Large-Scale All-Inverter Power Systems

Low-order frequency response models for power systems have a decades-long history in optimization and control problems such as unit commitment, economic dispatch, and wide-area control. With a few exceptions, these models are built upon the Newtonian mechanics of synchronous generators, assuming that the frequency dynamics across a system are approximately homogeneous, and assume the dynamics of nodal voltages for most operating conditions are negligible, and thus are not directly computed at all buses. As a result, the use of system frequency models results in the systematic underestimation of frequency minimum nadir and maximum RoCoF, and provides no insight into the reactive power-voltage dynamics. This paper proposes a low-order model of both frequency and voltage response in grid-forming inverter-dominated power systems. The proposed model accounts for spatial-temporal variations in frequency and voltage behavior across a system and as a result, demonstrates the heterogeneity of frequency response in future renewable power systems. Electromagnetic transient (EMT) simulations are used to validate the utility, accuracy, and computational efficiency of these models, setting the basis for them to serve as fast, scalable alternatives to EMT simulation, especially when dealing with very large-scale systems, for both planning and operational studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Integrated Lift- and Wall-Force Closures in Computational Fluid Dynamics (CMFD)

In this report, we provide the necessary background to put this work in context and for a more detailed background, the reader is referred to CASL 2013 report. Central constraints on the design of LWRs are thermal-structural limits, e.g., maximum fuel pin temperature and maximum allowable working pressure in a pressure vessel. In both of these cases, the knowledge of the spatially and temporally varying two-phase flow is required to make accurate predictions as the distribution of both the phases have an impact on heat extraction as well as neutron moderation. Engineered two-phase flow systems have complex geometry and often high Reynolds numbers which increase the computational cost of an analysis, and generally make direct numerical simulations of the flow field impossible even on the fastest computers today. Consequently, a number of models which seek to reduce this computational cost have been created over the years. One of these called the two-fluid model has become the standard in three dimensional CFD (Computational Fluid Dynamics) codes and is and is also referred to as CMFD (Computational Multiphase Fluid Dynamics). The two-fluid model splits each conserved quantity into two fields, one for each phase that is treated as interspersed continua. Splitting one field into two requires additional interfacial closure relationships. In the two-fluid model, the momentum flux between the phases is governed by a number of these closure relationships. Generally, the net force on the dispersed phase is decomposed into inertial, added mass, buoyancy, drag, lift, and wall forces. In addition, there is a time-dependent Basset force, as well as turbulent dispersion effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Relative Hausdorff distance for network analysis

Similarity measures are used extensively in machine learning and data science algorithms. The newly proposed graph Relative Hausdorff (RH) distance is a lightweight yet nuanced similarity measure for quantifying the closeness of two graphs. In this work we study the effectiveness of RH distance as a tool for detecting anomalies in time-evolving graph sequences. We apply RH to cyber data with given red team events, as well to synthetically generated sequences of graphs with planted attacks. In our experiments, the performance of RH distance is at times comparable, and sometimes superior, to graph edit distance in detecting anomalous phenomena. Furthermore, our results suggest that in appropriate contexts, RH distance has advantages over more computationally intensive similarity measures.

97 MATHEMATICS AND COMPUTING↗

Noise correction of turbulent spectra obtained from Acoustic Doppler Velocimeters

Accurately estimated auto-spectral density functions are essential for characterization of turbulent flows, and they also have applications in computational fluid dynamics modeling, site and inflow characterization for hydrokinetic turbines, and inflow turbulence generation. The Acoustic Doppler Velocimeter (ADV) provides single-point temporally resolved data, that are used to characterize turbulent flows in rivers, seas, and oceans. However, ADV data are susceptible to contamination from various sources, including instrument noise, which is the intrinsic limit to the accuracy of acoustic velocity measurements. Due to the presence of instrument noise, the spectra obtained are altered at high frequencies. The focus of this study is to develop a robust and effective method for accurately estimating auto-spectral density functions from ADV data by reducing or removing the spectral contribution derived from instrument noise. For this purpose, the “Noise Auto-Correlation” (NAC) approach was developed, which exploits the correlation properties of instrument noise to identify and remove its contribution from spectra. The spectra estimated using the NAC approach exhibit increased fidelity and a slope of -5/3 in the inertial range, which is typically observed for turbulent flows. Finally, this study also compares the effectiveness of low-pass Gaussian filters in removing instrument noise with that of the NAC approach. For the data used in this study, both the NAC and Gaussian filter approaches are observed to be capable of removing instrument noise at higher frequencies from the spectra. However, the NAC results are closer to the expected frequency power of -5/3 in the inertial sub-range.

Acoustic Doppler Velocimeter↗

Plasticity of irradiated materials at the nano and micro-scales

Here, we review here our recent work on plastic deformation in irradiated materials at the nano- and micro-scales, as revealed by Discrete Dislocation Dynamics (DDD) simulations. Two methods of including irradiation effects in the DDD framework are presented. The first directly captures the atomistic interaction mechanisms, while the second can effectively study high-dose irradiation. Computer simulations lead to new understanding of the dynamics of collective dislocation-irradiation defect interactions, as well as the quantitative analysis of the temporal and spatial characteristics associated with plastic instabilities. Based on these insights, theoretical models are developed to predict the critical conditions for dislocation channel formation. A simple probability model is proposed and demonstrated to predict the width of dislocation channels in bulk irradiated materials with good agreement with experimental data. The fundamental understanding of the origins of plastic flow localization in irradiated materials sheds light on the design of future generations of radiation-resistant materials.

36 MATERIALS SCIENCE↗

A simple three-component mixing problem for the evaluation of a new reaction rate model

A simple computational mixing problem is presented which can be utilized to assess the behavior of Reynolds-averaged reaction rate models in a problem with temporally varying mixedness. In this problem, three mixing components are homogeneously distributed but initially separated in a triply periodic domain. Further, these components are initialized within a Taylor–Green-like velocity field, which creates a mixing history evolving from the so-called “no-mix limit” to a well-mixed state. Large-eddy simulation results from this problem in configurations involving both premixed and nonpremixed reactants are then compared with zero-dimensional Reynolds-averaged Navier–Stokes results utilizing a new model for multicomponent reacting mixtures. The new model is shown to appropriately respect the no-mix limit and outperforms an earlier model (Morgan, 2022), particularly at early times when components are near the no-mix limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Time-coded aperture for x-ray imaging

Changing conditions and relative motion of an experiment’s components during a single exposure causes motion blur for X-ray imaging techniques. In this letter, we describe a computational imaging system, where image capture and post-processing are co-designed to improve final image quality. The idea is based on temporally encoding the motion by fluttering the detector shutter open and closed with a binary pseudo-random sequence for guaranteeing an invertible motion blur. This work is especially relevant to high-speed X-ray imaging applications where time-resolution is limited by the detector and/or available photon flux.

47 OTHER INSTRUMENTATION↗

Evaluation of Methods for Comparison of Spatiotemporal and Time Series Datasets

Data used by the National Renewable Energy Laboratory (NREL) in energy analysis are often produced by industry and licensed or purchased for analysis. While this practice provides needed flexibility in selecting data for analysis it presents challenges in understanding the differences among multiple, ostensibly similar, datasets. As options for source data become more varied, it is important to be able to articulate why certain datasets were chosen and to ensure those include the data that best meet the boundaries and/or limitations of a particular analysis. This report represents the first of three phases of research intended to develop methods to quantitatively assess and compare both input datasets and the results of analyses performed at NREL. This capability is critical to identifying tipping points in the costs or benefits of achieving high spatial and temporal resolution of input data.

97 MATHEMATICS AND COMPUTING↗

A harmonic balance method for PDEs

In this report, we describe some approaches to calculate the non-linear system of equations prescribed by the harmonic balance method (HB), a frequency domain analysis technique for modelling a non-linear system of partial differential equations (PDEs). The approach which we ultimately pursue can be seen as a time-collocation approach, except that the harmonic balance equations are obtained weakly (in the sense used in the calculus of variations). This weak formulation allows us to adapt existing transient or stationary PDEs models in the Panzer/Trilinos framework for frequency domain analysis via the harmonic balance method. We begin with a motivatation for the harmonic balance method and outline its mathematical formulation. We then describe some approaches to calculate the harmonic balance formulae, and their means of implementation through the modification of a Panzer tutorial problem - a stationary Helmholtz equation with a constant Dirichlet boundary condition and a non-linear source. For each of these approaches, we outline the necessary adaptations to solve the corresponding (periodically) transient Helmholtz equation with a (temporally) periodic Dirichlet boundary condition and non-linear source.

97 MATHEMATICS AND COMPUTING↗

Data-Driven Method for Groundwater-Level Mapping and Monitoring-Well Network Optimization at Hanford

This report summarizes the initial results and outcomes of a physics-informed, data-driven groundwater level (GWL) mapping capability for the Hanford Site. GWL mapping at Hanford is typically conducted annually and requires a significant amount of computational and expert resources, and it does not allow assessment of the informational value of specific monitoring wells. The proposed method produces spatially and temporally resolved fields consistent with sparse, irregularly sampled, and nonuniformly distributed well measurements. Implemented successfully, this capability will allow rapid mapping of groundwater levels and provide an opportunity to optimize monitoring activities (both location and sampling frequency) based on data information value evaluation. The approach integrates a diffusion-based generative model – trained on MODFLOW simulation data from the Plateau-to-River (P2R) model – with score-based data assimilation (SDA), allowing observation-conditioned mapping without retraining for each monitoring-network layout.

54 ENVIRONMENTAL SCIENCES↗

RISK MONITORING CAPABILITES FROM DYNAMIC PRA DATA

Risk Monitors (RMs) are software tools that are employed in nuclear power plants to monitor the safety properties of the plant itself during operation. These safety properties are generated using classical PRA methods which are typically based on static Boolean structures such as Fault-Trees (FTs) and Event-Trees (ETs). The outcomes of these methods are Core Damage Frequency (CDF) and risk importance values of plant components (e.g., valves and pumps). Examples of Risk Importance Measures (RIMs) commonly used are: Risk Achievement Worth (RAW), Risk Reduction Worth (RRW), Birnbaum (B) and Fussell-Vesely (FV) [1,2]. In a RM these values are updated in real time once plant components status is updated (e.g., under maintenance, off-line, testing). In addition, RMs can be used to optimize maintenance schedules so those safety features of the plant are below plant standard values. This paper presents an approach to create temporal profiles of RIMs (i.e., RM capabilities) given a temporal change of the plant configuration. Instead of employing data generated by ETs and FTs, the proposed approach will employ data generated by Dynamic PRA methods. In contrast to classical PRA methods, Dynamic PRA methods couple stochastic methods with safety analysis codes to determine risk associate to complex systems such as nuclear plants. Compared to classical PRA methods, they can evaluate with higher resolution the safety impact of timing and sequencing of events on the accident progression. We will describe how classical RIMs can be generated by Dynamic PRA data and then we will show how the temporal profiles of the RIMs can be generated given a temporal profile of the system configuration.

97 MATHEMATICS AND COMPUTING↗

Probabilistic Groundwater Modeling of the RDX Plume at Los Alamos National Laboratory to Support Risk Assessment - 20359

A plume of the contaminant hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) with concentrations greater than the New Mexico tap water drinking standard (9.66 ppb) is present in the regional aquifer near the southwestern boundary of Los Alamos National Laboratory (LANL). A risk assessment is performed for exposure to regional aquifer groundwater, with long-term predictions of RDX concentrations provided by a calibrated, probabilistic, numerical fate and transport model. The structure of the model is hierarchical, with the RDX Regional Aquifer (RA) groundwater model acting as the primary tool for analysis of downgradient RDX concentrations. The RA model is deeply informed by the conceptual site model (CSM) and is calibrated using site RDX concentration data and hydraulic head measurements, along with other analyses. The model is calibrated using data through December 2019. Where data are scarce other lines of evidence are used to inform inputs, including the multiphase RDX Vadose Zone (VZ) model and Pipe and Disk analytical screening tool. Model inputs are described with informative prior distributions using a robust approach to development that incorporates all available lines of evidence for every parameter used as an input to the RA model. Calibration is performed using a classical nonlinear optimization routine, which is then used to initialize a Bayesian calibration. The Bayesian calibration constrains the uncertainty in the classical calibration, ultimately providing posterior distributions for all model parameters. The challenges of the calibration include high-dimensional parameter space, including spatially heterogeneous hydraulic conductivities, comparatively sparse data, and low RDX concentrations. Posterior parameter distributions developed in the calibration process are then used for stochastic predictive model runs into the future. The result of the forward runs is spatially and temporally explicit estimates of head and concentration with uncertainty at all points in space and time. The probabilistic modeling approach presented here includes innovative computational and statistical methods that leverage high-performance computing (HPC) resources. It makes use of prior modeling work performed at the chromium plume site in the central LANL area, with extensive updates. The risk assessment will ultimately be used to support decision-making at the site, using multiple appropriately weighted sources of information, as well as uncertainty, sensitivity, and value-of-information analyses. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Response of Subsurface Nitrogen-Cycling Microbial Communities to Environmental Fluctuations (Final Technical Report)

Riparian floodplains are dynamic ecosystems linking terrestrial and riverine systems. These floodplains experience hydrological shifts such as changes in water table height, flooding, and drought and can be ‘hotspots’ of biogeochemical cycling due to shifting sediment moisture (and saturation) and subsurface exchanges of water, nutrients, and other compounds across different sediment layers. Subsurface microbial communities are the primary drivers of biogeochemical processes in floodplains, and thus their structure and function can directly influence both surface and groundwater quality. The microbial nitrogen (N) cycle is particularly important in floodplains as it affects nutrient availability and removal. Two functional guilds of chemoautotrophic (i.e. CO2-fixing) microorganisms are responsible for the first oxidative step of the N cycle, nitrification: ammonia-oxidizing archaea (AOA) and bacteria (AOB) catalyze the oxidation of ammonia to nitrite, while nitrite-oxidizing bacteria (NOB) oxidize nitrite to nitrate. Despite the critical role nitrification plays in N-cycling in both terrestrial and aquatic ecosystems, our understanding of the diversity, ecophysiology, and activity of nitrifying organisms in subsurface floodplain soils/sediments is extremely limited. To help address this critical knowledge gap, the overarching goal of this project was to determine how shifts in key environmental parameters and gradients impact microbial N-cycling communities/processes, with particular emphasis on nitrification, within hydrologically-variable floodplain sediments in the Wind River Basin near Riverton, Wyoming. The three specific objectives of this project were to: (1) to associate in situ environmental drivers of N cycling with distinct functional guilds; (2) determine the guild response to variation in key ecosystem drivers; and (3) develop a dynamic ecosystem model of the microbial N cycle with the Riverton subsurface using community genomic and biogeochemical data collected in the first two objectives. Over the course of this project, we employed both 16S rRNA gene amplicon sequencing and genome-resolved metagenomics to examine the phylogenetic diversity and metabolic potential of subsurface nitrifier communities within 68 samples collected across multiple sites, depths, and time points within the Riverton floodplain, allowing for both spatial and temporal investigations at different scales. This project benefitted tremendously from recent advances in high-throughput sequencing technologies coupled with dramatic improvements in the computational tools and algorithms available for analyzing such large, complex genomic datasets. By pairing these cutting-edge genomic approaches with depth-resolved sampling and detailed geochemical analyses of the Riverton floodplain, we have gained novel insights into the structure and function of subsurface nitrifier communities in relation to both hydrology and biogeochemistry. This project resulted in the most detailed and comprehensive characterization of N-cycling floodplain microbial communities to date and will hopefully inspire and pave the way for future studies using similar approaches in other floodplains. Indeed, such information is critical for understanding subsurface biogeochemical cycling and how elemental stores are altered from perturbations initiated by the water cycle within floodplains. Finally, because of the terrestrial-aquatic nature of the Riverton floodplain, results from this project are also of relevance to disciplines such as soil science, estuarine science, limnology & oceanography, biogeochemistry, geobiology, environmental engineering, as well as genomics and data science.

54 ENVIRONMENTAL SCIENCES↗

Model-Based Approaches to Generate Knowledge from Data in a Plant Reliability Context

One challenge that nuclear power plant system engineers are facing is continuous generation of an extremely large amount of equipment reliability (ER) data. These data elements come in textual (e.g., condition reports) and numeric (e.g., generated by monitoring systems) forms. They provide system engineers with valuable insights and information by discovering anomalous behaviors or degradation trends, identifying possible causes behind such behaviors and trends, and predicting their direct consequences. This paper directly targets the knowledge generation from ER data by putting “data into context.” We employ model-based system engineering (MBSE) of systems and assets to represent and capture their architecture and functional (i.e., cause-effect) relations. ER data elements are processed by first identifying which of the developed MBSE elements they are referring to. This task is harder for textual data since the information contained in issue or maintenance reports needs to be “understood” by a computational tool. We called this process “knowledge extraction” since our methods extract knowledge from textual data. Last, once numeric and textual ER data elements have been processed and “understood,” we discover possible cause-effect relations among them. This is performed by observing whether a logical connection through the MBSE models exists, and if there is a temporal relationship among them. The logic and temporal are the two main ingredients to perform “machine reasoning” from ER data.

97 - MATHEMATICS AND COMPUTING↗

Using gap-filled MAIAC AOD and WRF-Chem to estimate daily PM 2.5 concentrations at 1 km resolution in the Eastern United States

To link short-term exposures of air pollutants to health outcomes, scientists must use high temporal and spatial resolution estimates of PM 2.5 concentrations. We develop a daily PM 2.5 product at 1 × 1 km 2 spatial resolution across the eastern United States (east of 90° W) with the aid of 1 × 1 km 2 MAIAC aerosol optical depth (AOD) data, 36 × 36 km 2 WRF-Chem output, 1 × 1 km 2 land-use type from the National Land Cover Database, and 0.125° × 0.125° ERA-Interim re-analysis meteorology. A gap-filling technique is applied to MAIAC AOD data to construct robust daily estimates of AOD when the satellite data are missing (e.g., areas obstructed by clouds or snow cover). The input data are incorporated into a multiple-linear regression model trained to surface observations of PM 2.5 from the EPA Air Quality System (AQS) monitoring network. The model generates a high-fidelity estimate (r 2 = 0.75 using a 10-fold random cross-validation) of daily PM 2.5 throughout the eastern United States. Of the inputs to the statistical model, WRF-Chem output (r 2 = 0.66) is the most important contributor to the skill of the model. MAIAC AOD is also a strong contributor (r 2 = 0.52). Daily PM 2.5 output from our statistical model can be easily integrated into county-level epidemiological studies. The novelty of this project is that we are able to simulate PM 2.5 in a computationally efficient manner that is constrained to ground monitors, satellite data, and chemical transport model output at high spatial resolution (1 × 1 km 2 ) without sacrificing the temporal resolution (daily) or spatial coverage (>2,000,000 km 2 ).

47 OTHER INSTRUMENTATION↗

Model-Based Approaches to Generate Knowledge from Data in a Plant Reliability Context

One challenge that nuclear power plant system engineers are facing is that the amount of equipment reliability (ER) data being continuously generated are extremely large. These data elements come in different forms: textual (e.g., condition reports) and numeric (e.g., generated by monitoring systems) and they provide system engineers with valuable insights and information regarding the discovery of anomalous behaviors or degradation trends, the identification of the possible causes behind such behaviors and trends, and the prediction of their direct consequences. This paper directly targets the generation of knowledge from ER data by putting “data into context”. Here, we employ model-based system engineering (MBSE) models of systems and assets to represent and capture their architecture and functional (i.e., cause-effect) relations. ER data elements are processed by identifying first which elements of the developed MBSE elements they are referring to. This task is much harder for textual data since the information contained in issue or maintenance reports needs to “be understood” by a computational tool. Here we called this process “knowledge extraction” where our methods to extract knowledge from textual data. Lastly, once numeric and textual ER data elements have been processed and “understood”, we discover possible cause-effect relations among them. This is performed by observing if a logical connection through the MBSE models exists, and if there is a temporal relation among them. The logic and temporal are the two main ingredients to perform “machine reasoning” from ER data.

97 MATHEMATICS AND COMPUTING↗

Parallel-In-Time Multigrid with Adaptive Spatial Coarsening for The Linear Advection and Inviscid Burgers Equations

We apply a multigrid reduction-in-time (MGRIT) algorithm to hyperbolic partial differential equations in one spatial dimension. This study is motivated by the observation that sequential time-stepping is a computational bottleneck when attempting to implement highly concurrent algorithms; thus parallel-in-time methods are desirable. MGRIT adds parallelism by using a hierarchy of successively coarser temporal levels to accelerate the solution on the finest level. In the case of explicit time-stepping, spatial coarsening is a suitable approach to ensure that stability conditions are satisfied on all levels, and it may be useful for implicit time-stepping by producing cheaper multigrid cycles. Unfortunately, uniform spatial coarsening results in extremely slow convergence when the wave speed is near zero, even if only locally. We present an adaptive spatial coarsening strategy that addresses this issue for the variable coefficient linear advection equation and the inviscid Burgers equation using first-order explicit or implicit time-stepping methods. Serial numerical results show this method offers significant improvements over uniform coarsening and is convergent for the inviscid Burgers equation with and without shocks. Parallel scaling tests on up to 128K cores indicate that run-time improvements over serial time-stepping strategies are possible when spatial parallelism alone saturates, and that scalability is robust for oscillatory solutions which change on the scale of the grid spacing.

97 MATHEMATICS AND COMPUTING↗

DiffESM: Conditional Emulation of Temperature and Precipitation in Earth System Models With 3D Diffusion Models

Earth system models (ESMs) are essential for understanding the interaction between human activities and the Earth's climate. However, the computational demands of ESMs often limit the number of simulations that can be run, hindering the robust analysis of risks associated with extreme weather events. While low-cost climate emulators have emerged as an alternative to emulate ESMs and enable rapid analysis of future climate, many of these emulators only provide output on at most a monthly frequency. This temporal resolution is insufficient for analyzing events that require daily characterization, such as heat waves or heavy precipitation. We propose using diffusion models, a class of generative deep learning models, to effectively downscale ESM output from a monthly to a daily frequency. Trained on a handful of ESM realizations, reflecting a wide range of radiative forcings, our DiffESM model takes monthly mean precipitation or temperature as input, and is capable of producing daily values with statistical characteristics close to ESM output. Combined with a low-cost emulator providing monthly means, this approach requires only a small fraction of the computational resources needed to run a large ensemble. We evaluate model behavior using a number of extreme metrics, showing that DiffESM closely matches the spatio-temporal behavior of the ESM output it emulates in terms of the frequency and spatial characteristics of phenomena such as heat waves, dry spells, or rainfall intensity.

54 ENVIRONMENTAL SCIENCES↗