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 1,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Parallel-in-Time Solution of Hyperbolic PDE Systems via Characteristic-Variable Block Preconditioning

We consider the parallel-in-time solution of both linear and nonlinear hyperbolic partial differential equation (PDE) systems in one spatial dimension. In the nonlinear setting, the discretized equations are solved with a preconditioned residual iteration based on a global linearization. The linear(ized) equation systems are approximately solved parallel-in-time using a block preconditioner applied in the characteristic variables of the underlying linear(ized) hyperbolic PDE. This change of variables is motivated by the observation that intervariable coupling between characteristic variables is weak, at least locally where spatio-temporal variations in the eigenvectors of the associated flux Jacobian are sufficiently small, while that between the original variables is not. For an ℓ-dimensional system of PDEs, applying the preconditioner consists of solving a sequence of ℓ scalar linear(ized)-advection-like problems, each associated with a different characteristic wave-speed in the underlying linear(ized) PDE. Furthermore, we approximately solve these linear advection problems using multigrid reduction-in-time (MGRIT); however, any other suitable parallel-in-time method could be used. Numerical examples are shown for the (linear) acoustics equations in heterogeneous media and for the (nonlinear) shallow water equations and Euler equations of gas dynamics with shocks and rarefactions. For many test problems, the solver converges in just a handful of iterations and with mesh-independent convergence rates.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Temporal evolution of soil moisture statistical fractal and controls by soil texture and regional groundwater flow

Soil moisture statistical fractal is an important tool for downscaling remotely-sensed observations and has the potential to play a key role in multi-scale hydrologic modeling. The fractal was first introduced two decades ago, but relatively little is known regarding how its scaling exponents evolve in time in response to climatic forcings. Previous studies have neglected the process of moisture re-distribution due to regional groundwater flow. Here, we used a physically-based surface-subsurface processes model and numerical experiments to elucidate the patterns and controls of fractal temporal evolution in two U.S. Midwest basins. Groundwater flow was found to introduce large-scale spatial structure, thereby reducing the scaling exponents (τ), which has implications for the transferability of calibrated parameters to predict τ. However, the groundwater effects depend on complex interactions with other physical controls such as soil texture and land use. The fractal scaling exponents, while in general showing a seasonal mode that correlates with mean moisture content, display hysteresis after storm events that can be divided into three phases, consistent with literature findings: (a) wetting, (b) re-organizing, and (c) dry-down. Modeling experiments clearly show that the hysteresis is attributed to soil texture, whose "patchiness" is the primary contributing factor. We generalized phenomenological rules for the impacts of rainfall, soil texture, groundwater flow, and land use on τevolution. Grid resolution has a mild influence on the results and there is a strong correlation between predictions of τ from different resolutions. Overall, our results suggest that groundwater flow should be given more consideration in studies of the soil moisture statistical fractal, especially in regions with a shallow water table.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A dictionary learning algorithm for compression and reconstruction of streaming data in preset order

There has been an emerging interest in developing and applying dictionary learning (DL) to process massive datasets in the last decade. Many of these efforts, however, focus on employing DL to compress and extract a set of important features from data, while considering restoring the original data from this set a secondary goal. On the other hand, although several methods are able to process streaming data by updating the dictionary incrementally as new snapshots pass by, most of those algorithms are designed for the setting where the snapshots are randomly drawn from a probability distribution. In this paper, we present a new DL approach to compress and denoise massive dataset in real time, in which the data are streamed through in a preset order (instances are videos and temporal experimental data), so at any time, we can only observe a biased sample set of the whole data. Here, our approach incrementally builds up the dictionary in a relatively simple manner: if the new snapshot is adequately explained by the current dictionary, we perform a sparse coding to find its sparse representation; otherwise, we add the new snapshot to the dictionary, with a Gram-Schmidt process to maintain the orthogonality. To compress and denoise noisy datasets, we apply the denoising to the snapshot directly before sparse coding, which deviates from traditional dictionary learning approach that achieves denoising via sparse coding. Compared to full-batch matrix decomposition methods, where the whole data is kept in memory, and other mini-batch approaches, where unbiased sampling is often assumed, our approach has minimal requirement in data sampling and storage: i) each snapshot is only seen once then discarded, and ii) the snapshots are drawn in a preset order, so can be highly biased. Through experiments on climate simulations and scanning transmission electron microscopy (STEM) data, we demonstrate that the proposed approach performs competitively to those methods in data reconstruction and denoising.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Stride search: A general algorithm for storm detection in high-resolution climate data

This study discusses the problem of identifying extreme climate events such as intense storms within large climate data sets. The basic storm detection algorithm is reviewed, which splits the problem into two parts: a spatial search followed by a temporal correlation problem. Two specific implementations of the spatial search algorithm are compared: the commonly used grid point search algorithm is reviewed, and a new algorithm called Stride Search is introduced. The Stride Search algorithm is defined independently of the spatial discretization associated with a particular data set. Results from the two algorithms are compared for the application of tropical cyclone detection, and shown to produce similar results for the same set of storm identification criteria. Differences between the two algorithms arise for some storms due to their different definition of search regions in physical space. The physical space associated with each Stride Search region is constant, regardless of data resolution or latitude, and Stride Search is therefore capable of searching all regions of the globe in the same manner. Stride Search's ability to search high latitudes is demonstrated for the case of polar low detection. Wall clock time required for Stride Search is shown to be smaller than a grid point search of the same data, and the relative speed up associated with Stride Search increases as resolution increases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

The RISMC Approach to Perform Advanced PRA Analyses

The existing fleet of nuclear power plants is in the process of extending its lifetime and increasing the power generated from these plants via power uprates. In order to evaluate the impact of these two factors on the safety of the plant, the RISMC Pathway aims to develop simulation-based tools and methods to assess risks for existing nuclear power plants in order to optimize safety. This pathway, by developing new methods, is extending the state-of-the-practice methods that have been traditionally based on logic structures such as Event-Trees and Fault-Trees. These static types of models mimic system response in an inductive and deductive way respectively, yet are restrictive in the ways they can represent spatial and temporal constructs. RISMC analyses are performed by using a combination of thermal-hydraulic codes and a stochastic analysis tool currently under development at the Idaho National Laboratory, i.e. RAVEN. This paper presents a case study in order to show the capabilities of the RISMC methodology to assess impact of power uprate of a boiling water reactor system during a station blackout accident scenario. We employ the system simulator code, RELAP5-3D, coupled with RAVEN which perform the stochastic analysis. Our analysis is in fact performed by: 1) sampling values of a set of parameters from the uncertainty space of interest, 2) simulating the system behavior for that specific set of parameter values and 3) analyzing the set of simulation runs. Results obtained give a detailed investigation of the issues associated with a plant power uprate including the effects of station blackout accident scenarios. We are able to quantify how the timing of specific events was impacted by a higher nominal reactor core power. Such safety insights can provide useful information to the decision makers to perform risk informed margins management.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Advancing spatiotemporal forecasts of CO 2 plume migration using deep learning networks with transfer learning and interpretation analysis

Accurate and timely forecasts of CO 2 plume distribution throughout the injection and post-injection phases are crucial for detecting plume migration, assessing leakage risks, and supporting operational decisions in geologic carbon storage (GCS). Current convolutional neural network-based approaches primarily focus on spatial information and overlook temporal dependencies in plume distributions, thus limiting their ability to capture dynamic movement effects and provide accurate predictions of plume migration. In this work, we propose two deep learning models, Auto-Encoder (AE)-LSTM and Encoder-Decoder (ED)-ConvLSTM, each uniquely designed to capture both spatial and temporal features. We apply the proposed methods to forecast the dynamic distribution of CO 2 plumes based on 108 reservoir simulations over a 30-year injection and a 30-year post-injection period. The results indicate that the ED-ConvLSTM model outperforms the AE-LSTM model in accurately predicting the spatiotemporal dynamics of CO 2 plume migration, achieving R 2 values above 0.99. To provide a deeper understanding of these model predictions, we employ a gradient-based explanation method on the trained models. This approach provides insights into the influence of input variables on plume migration forecasts and uncovers the underlying prediction mechanisms of the proposed models. Furthermore, we introduce a transfer learning technique, enabling fast and accurate plume migration forecasting in the post-injection phase by leveraging the trained model during the injection phase. This reduces the necessity for extensive data collection or re-training. In conclusion, the methods proposed in our work enhances the performance and interpretability of CO 2 plume migration forecasts, thereby facilitating informed decision-making throughout the entire lifecycle of GCS applications.

58 GEOSCIENCES↗
DOE OSTI2024

Dashboard for Marine Energy Site Assessment and Monitoring

The marine energy (ME) industry presently relies upon fragmented site assessment solutions that require high resource expenditure for deployment at each site and do not leverage the wealth of readily available tools and information. A wave energy resource assessment dashboard, currently in development, will substantially improve siting, permitting, operations, and maintenance of ME projects by providing an integrated solution that is a one-stop-shop for a developer’s needs. The Site Energy Assessment and MOnitoring Dashboard (SEAMOD) will be of commercial interest to anyone seeking to deploy an ME project and is easily expandable to include tidal and wind energy site assessments. The integrated dashboard is being developed using state-of-the-art database and cloud computing methods and data-assimilative modeling tools that can be coupled with low-cost, rapidly deployable wave buoys and environmental sensing hardware. The combined software and hardware dashboard will reduce wave energy site characterization and wave climate monitoring costs by more than 60 percent and provide assessments that meet international industry standards. To realize a thriving global ME industry, the physical environment at a potential deployment site must be understood, not only for resource characterization, but also for optimization of device and power conversion performance. SEAMOD directly addresses these needs with a commercially marketable product. SEAMOD is a low-cost solution that provides comprehensive ME resource assessments, baseline environmental monitoring, and offshore characterizations required for successful ME development. The key technical objectives for Phase I were a series of software development goals, which when implemented with monitoring solutions, produced an initial proof-of-concept low-cost wave energy resources dashboard. In Phase II, the development of the prototype SEAMOD continued. The basic framework employed was the development of a revised dashboard and monitoring tool customized for ME applications by focusing on IEC site assessment and method requirements. Development was focused on the integration of full hindcast metocean products to provide hindcast resource characterization and environmental information. The final integrated dashboard provides a low-cost solution that delivers comprehensive, scalable, industry-standard energy resource assessments and offshore characterizations required for successful ME development. The integrated dashboard offers visibility of the most recent site modeling, measurements, and historical data. The application and integration of consensus-based standards for wave energy resource assessment, as determined by the International Electrotechnical Commission (IEC), are crucial for the impact and value of SEAMOD. SEAMOD includes monthly, seasonal, and yearly statistics, as well as the total 30-year record, offering temporal resolution of the IEC parameters to aid potential developers in determining the available wave energy resources in their area of interest.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Verification and validation of the open-source plasma fluid code: Zapdos

Zapdos is an open-source finite element plasma fluid solver based on the MOOSE multiphysics framework. This paper outlines Zapdos verification, benchmarking, and validation efforts for 1D and 2D RF capacitively coupled plasma discharge models for mid-range pressures (0.1 - 1 Torr). The verification process involved using the method of manufactured solutions to assess Zapdos spatial and temporal error convergence. L2 errors ranged from 10 -2 to 10 -4 , while the convergence’s slope were in agreement with the predicted slopes for the tested variable and time integration orders. The benchmarking process involved comparisons to previously results from the validated finite element code, LSODI. These works included 1D and 2D simulations for a range of plasma parameters (densities, temperatures, voltage, etc.). For the 1D cases, Zapdos and LSODI results were in very good agreement. In the 2D cases, variable behaviors matched, with slight discrepancies in peak values. The validation process involved comparisons to experimental works including electron density measurements by microwave interferometry and metastable density measurements by planar laser-induced fluorescence imaging. Results shown reasonable agreement at higher pressure, with results starting to diverge at low pressures. Furthermore, probable causes for this diverges are the limitation of the fluid assumption for plasmas at low pressure, or the need for more robust boundary conditions. Overall, Zapdos shown reasonable results for the verification, benchmarking, and validation efforts, and Zapdos can be downloaded at https://github.com/shannon-lab/zapdos.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Maximum bound principle preserving integrating factor Runge–Kutta methods for semilinear parabolic equations

A large class of semilinear parabolic equations satisfy the maximum bound principle (MBP) in the sense that the time-dependent solution preserves for any time a uniform pointwise bound imposed by its initial and boundary conditions. Here, the MBP plays a crucial role in understanding the physical meaning and the well-posedness of the mathematical model. Investigation on numerical algorithms with preservation of the MBP has attracted increasingly attentions in recent years, especially for the temporal discretizations, since the violation of MBP may lead to nonphysical solutions or even blow-ups of the algorithms. In this paper, we study high-order MBP-preserving time integration schemes by means of the integrating factor Runge-Kutta (IFRK) method. Beginning with the space-discrete system of semilinear parabolic equations, we present the IFRK method in general form and derive the sufficient conditions for the method to preserve the MBP. In particular, we show that the classic four-stage, fourth-order IFRK scheme is MBP preserving for some typical semilinear systems although not strong stability preserving, which can be instantly applied to the Allen-Cahn type of equations. To our best knowledge, this is the first time to present a fourth-order linear numerical method preserving the MBP. In addition, convergence of these numerical schemes is proved theoretically and verified numerically, as well as their efficiency by simulations of 2D and 3D long-time evolutional behaviors. Numerical experiments are also carried out for a model which is not a typical gradient flow as the Allen-Cahn type of equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Analysis of the Temporal Response of Coupled Asymmetrical Zero-Power Subcritical Bare Metal Reactor Systems

The behavior of symmetrical coupled-core systems has been extensively studied, yet there is a dearth of research on asymmetrical systems due to the increased complexity of the analysis of such systems. In this research, the multipoint kinetics method is applied to asymmetrical zeropower, subcritical, bare metal reactor systems. Existing research on asymmetrical reactor systems assumes symmetry in the neutronic coupling; however, it will be shown that this cannot always be assumed. Deep subcriticality adds another layer of complexity and requires modification of the multipoint kinetics equations to account for the effect of the external neutron source. A modified set of multipoint kinetics equations is derived with this in mind. Subsequently, the Rossi-alpha equations are derived for a two-region asymmetrical reactor system. The predictive capabilities of the radiation transport code MCNP6 for neutron noise experiments are shown in a comparison to the results of a series of Rossi-alpha measurements performed by J. Mihalczo utilizing a coupled set of symmetrical bare highly-enriched uranium (HEU) cylinders. The ptrac option within MCNP6 can generate time-tagged counts in a cell (list-mode data). The list-mode data can then be processed similarly to measured data to obtain values for system parameters such as the dual prompt neutron decay constants observable in a coupled system. The results from the ptrac simulations agree well with the historical measured values. A series of case studies are conducted to study the effects of geometrical asymmetry in the coupling between two bare metal HEU cylinders. While the coupling behavior of symmetrical systems has been reported on extensively, that of asymmetrical systems remains sparse. In particular, it appears that there has been no previous research in obtaining the coupling time constants for asymmetrically-coupled systems. The difficulty in observing such systems is due in part to the inability to determine the individual coupling coefficients from measurement: unlike the symmetrical cases, only the product of the values can be obtained. A method is proposed utilizing MCNP6 tally ratios to separate the coupling coefficients for such systems. This work provides insight into the behavior of asymmetrically-coupled systems as the separation distance between the two cores is changed and also as the asymmetry is increased. As the asymmetry increases, both the slower and the faster observable prompt neutron decay constants increase in magnitude. The coupling time constants are determined from the measured decay constants. As the separation distance increases, both coupling coefficients decrease as expected. Based on these findings, an effective computational method utilizing MCNP6 and the Rossialpha technique can be applied to the prediction of asymmetrical coupled system measurements.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2019

Estimation of cerebral blood flow velocity during breath-hold challenge using artificial neural networks

The effect of untreated Obstructive Sleep Apnoea (OSA) on cerebral haemodynamics and CA impairment is an active field of research interest. A breath-hold challenge is usually used in clinical and research settings to simulate cardiovascular and cerebrovascular changes that mimic OSA events. This work utilises temporal arterial oxygen saturation (SpO 2 ) and photoplethysmography (PPG) signals to estimate the temporal cerebral blood flow velocity (CBFv) waveform. Measurements of CBFv, SpO 2 , and PPG, were acquired concurrently from volunteers performing two different protocols of breath-hold challenge in the supine position. Past values of the SpO 2 and PPG signals were used to estimate the current values of CBFv using different permutations and topologies of supervised learning with shallow artificial neural networks (ANNs). Measurements from one protocol were used to train the ANNs and find the optimum topologies, which in turn were tested using the other protocol. Data collected from 10 normotensive, healthy subjects (four females, age 28.5 ± 6.1 years, Body Mass Index (BMI) 24.0 ± 4.7 kg/m 2 ) were used in this study. The results show that different subjects have different optimum topologies for ANNs, thus indicating the effects of inter-subject variability on ANNs. Successfully reconstructed blind waveforms for the same subject group in the second protocol showed a reasonable accuracy of 60–80% estimation compared to the measured waveforms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Modeling and observations of North Atlantic cyclones: Implications for U.S. Offshore wind energy

To meet the Biden-Harris administration's goal of deploying 30 GW of offshore wind power by 2030 and 110 GW by 2050, expansion of wind energy into U.S. territorial waters prone to tropical cyclones (TCs) and extratropical cyclones (ETCs) is essential. This requires a deeper understanding of cyclone-related risks and the development of robust, resilient offshore wind energy systems. Here, this paper provides a comprehensive review of state-of-the-science measurement and modeling capabilities for studying TCs and ETCs, and their impacts across various spatial and temporal scales. We explore measurement capabilities for environments influenced by TCs and ETCs, including near-surface and vertical profiles of critical variables that characterize these cyclones. The capabilities and limitations of Earth system and mesoscale models are assessed for their effectiveness in capturing atmosphere–ocean–wave interactions that influence TC/ETC-induced risks under a changing climate. Additionally, we discuss microscale modeling capabilities designed to bridge scale gaps from the weather scale (a few kilometers) to the turbine scale (dozens to a few meters). We also review machine learning (ML)-based, data-driven models for simulating TC/ETC events at both weather and wind turbine scales. Special attention is given to extreme metocean conditions like extreme wind gusts, rapid wind direction changes, and high waves, which pose threats to offshore wind energy infrastructure. Finally, the paper outlines the research challenges and future directions needed to enhance the resilience and design of next-generation offshore wind turbines against extreme weather conditions.

17 WIND ENERGY↗
DOE OSTI2018

Frequency–wavenumber spectral analysis of spatio-temporal flows

We propose a fully spatio-temporal approach for identifying spatially varying modes of oscillation in fluid dynamics simulation output by means of multitaper frequency–wavenumber spectral analysis. One-dimensional spectrum estimation has proven to be a valuable tool in the analysis of turbulence data applied spatially to determine the rate of energy transport between spatial scales, or temporally to determine frequencies of oscillatory flows. It also allows for the quantitative comparison of flow characteristics between two scenarios using a standard basis. It has the limitation, however, that it neglects coupling between spatial and temporal structures. Two-dimensional frequency–wavenumber spectral analysis allows one to decompose waveforms into standing or travelling variety. The extended higher-dimensional multitaper method proposed here is shown to have improved statistical properties over conventional non-parametric spectral estimators, and is accompanied by confidence intervals which estimate their uncertainty. Multitaper frequency–wavenumber analysis is applied to a canonical benchmark problem, namely, a direct numerical simulation of von Kármán vortex shedding off a square wall-mounted cylinder with two inflow scenarios with matching momentum-thickness Reynolds numbers $Re_{\unicode[STIX]{x1D703}}\approx 1000$ at the obstacle. Frequency–wavenumber analysis of a two-dimensional section of these data reveals that although both the laminar and turbulent inflow scenarios show a turbulent $-5/3$ cascade in wavenumber ( $\unicode[STIX]{x1D708}$ ) and frequency ( $f$ ), the flow characteristics differ in that there is a significantly more prominent discrete harmonic oscillation near $(f,\unicode[STIX]{x1D708})=(0.2,0.21)$ in wavenumber and frequency in the laminar inflow scenario than the turbulent scenario. This frequency–wavenumber pair corresponds to a travelling wave with velocity near one near the centre path of the vortex street.

Geoga, Christopher J. (ORCID:0000000250604799)↗
DOE OSTI2013

Respiratory effort correction strategies to improve the reproducibility of lung expansion measurements

Purpose: Four-dimensional computed tomography (4DCT) can be used to make measurements of pulmonary function longitudinally. The sensitivity of such measurements to identify change depends on measurement uncertainty. Previously, intrasubject reproducibility of Jacobian-based measures of lung tissue expansion was studied in two repeat prior-RT 4DCT human acquisitions. Difference in respiratory effort such as breathing amplitude and frequency may affect longitudinal function assessment. In this study, the authors present normalization schemes that correct ventilation images for variations in respiratory effort and assess the reproducibility improvement after effort correction.Methods: Repeat 4DCT image data acquired within a short time interval from 24 patients prior to radiation therapy (RT) were used for this analysis. Using a tissue volume preserving deformable image registration algorithm, Jacobian ventilation maps in two scanning sessions were computed and compared on the same coordinate for reproducibility analysis. In addition to computing the ventilation maps from end expiration to end inspiration, the authors investigated the effort normalization strategies using other intermediated inspiration phases upon the principles of equivalent tidal volume (ETV) and equivalent lung volume (ELV). Scatter plots and mean square error of the repeat ventilation maps and the Jacobian ratio map were generated for four conditions: no effort correction, global normalization, ETV, and ELV. In addition, gamma pass rate was calculated from a modified gamma index evaluation between two ventilation maps, using acceptance criterions of 2 mm distance-to-agreement and 5% ventilation difference.Results: The pattern of regional pulmonary ventilation changes as lung volume changes. All effort correction strategies improved reproducibility when changes in respiratory effort were greater than 150 cc (p < 0.005 with regard to the gamma pass rate). Improvement of reproducibility was correlated with respiratory effort difference (R = 0.744 for ELV in the cohort with tidal volume difference greater than 100 cc). In general for all subjects, global normalization, ETV and ELV significantly improved reproducibility compared to no effort correction (p = 0.009, 0.002, 0.005 respectively). When tidal volume difference was small (less than 100 cc), none of the three effort correction strategies improved reproducibility significantly (p = 0.52, 0.46, 0.46 respectively). For the cohort (N = 13) with tidal volume difference greater than 100 cc, the average gamma pass rate improves from 57.3% before correction to 66.3% after global normalization, and 76.3% after ELV. ELV was found to be significantly better than global normalization (p = 0.04 for all subjects, and p = 0.003 for the cohort with tidal volume difference greater than 100 cc).Conclusions: All effort correction strategies improve the reproducibility of the authors' pulmonary ventilation measures, and the improvement of reproducibility is highly correlated with the changes in respiratory effort. ELV gives better results as effort difference increase, followed by ETV, then global. However, based on the spatial and temporal heterogeneity in the lung expansion rate, a single scaling factor (e.g., global normalization) appears to be less accurate to correct the ventilation map when changes in respiratory effort are large.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2017

An agent-based stochastic Occupancy Simulator

Occupancy has significant impacts on building performance. However, in current building performance simulation programs, occupancy inputs are static and lack diversity, contributing to discrepancies between the simulated and actual building performance. This work presents an Occupancy Simulator that simulates the stochastic behavior of occupant presence and movement in buildings, capturing the spatial and temporal occupancy diversity. Each occupant and each space in the building are explicitly simulated as an agent with their profiles of stochastic behaviors. The occupancy behaviors are represented with three types of models: (1) the status transition events (e.g., first arrival in office) simulated with probability distribution model, (2) the random moving events (e.g., from one office to another) simulated with a homogeneous Markov chain model, and (3) the meeting events simulated with a new stochastic model. A hierarchical data model was developed for the Occupancy Simulator, which reduces the amount of data input by using the concepts of occupant types and space types. Finally, a case study of a small office building is presented to demonstrate the use of the Simulator to generate detailed annual sub-hourly occupant schedules for individual spaces and the whole building. The Simulator is a web application freely available to the public and capable of performing a detailed stochastic simulation of occupant presence and movement in buildings. Future work includes enhancements in the meeting event model, consideration of personal absent days, verification and validation of the simulated occupancy results, and expansion for use with residential buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Using X-ray radiography to study oxygen flow in a proton exchange membrane electrolyzer operating under balanced pressure conditions

Of the various water electrolyzer technologies, the proton exchange membrane electrolyzer (PEMWE) is one of the best solutions for producing clean hydrogen without releasing CO 2 . In order to allow for widespread use of clean hydrogen, it is necessary to decrease its cost, which is intrinsically related to system operation. Current PEMWE plants operate in differential mode, directly pressurizing hydrogen and benefiting from thermodynamic compression, which increases overall system efficiency. However, high differential pressure above 30 bar can cause membrane stress, resulting in membrane creeping and failure. Pressurizing the water and operating at balanced pressure allows hydrogen to be produced at higher pressures while preserving the integrity of the membrane and porous layers. Nevertheless, the impact of pressurizing water on PEMWE performance must be better understood to maximize performance under balanced pressure conditions. Here, this study examined the impact of water pressure on two-phase flow. A high-pressure electrolyzer setup was developed to perform operando X-ray radiography and examine oxygen transport with high temporal resolution. The imaging segmentation process, developed to capture bubble properties in the channel, was applied to a specific experiment. The results clearly showed that as pressure increased up to 30 bars, the initial bubbly flow transitioned to slug flow, which led to channel saturation with oxygen. This work demonstrates that two-phase flow in an electrolyzer can be studied using X-ray radiography, which has the advantages of fast measurements and the ability to probe dense materials, such as those required for pressurized electrolyzers.

Balanced high pressure operation↗
DOE OSTI2021

Displacement Analysis of Geothermal Field Based on PSInSAR And SOM Clustering Algorithms A Case Study of Brady Field, Nevada—USA

The availability of free and high temporal resolution satellite data and advanced SAR techniques allows us to analyze ground displacement cost-effectively. Our aim was to properly define subsidence and uplift areas to delineate a geothermal field and perform time-series analysis to identify temporal trends. A Persistent Scatterer Interferometry (PSI) algorithm was used to estimate vertical displacement in the Brady geothermal field located in Nevada by analyzing 70 Sentinel-1A Synthetic-Aperture Radar (SAR) images, between January 2017 and December 2019. To classify zones affected by displacement, an unsupervised Self-Organizing Map (SOM) algorithm was applied to classify points based on their behavior in time, and those clusters were used to determine subsidence, uplift, and stable regions automatically. Finally, time-series analysis was applied to the clustered data to understand the inflection dates. The maximum subsidence is –19 mm/yr with an average value of –6 mm/yr within the geothermal field. The maximum uplift is 14 mm/yr with an average value of 4 mm/yr within the geothermal field. The uplift occurred on the NE of the field, where the injection wells are located. On the other hand, subsidence is concentrated on the SW of the field where the production wells are located. The coupling of the PSInSAR and the SOM algorithms was shown to be effective in analyzing the direction and pattern of the displacements observed in the field.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

A split imaging spectrometer for temporally and spatially resolved titanium absorption spectroscopy

We present a temporally and a spatially resolved spectrometer for titanium x-ray absorption spectroscopy along 2 axial symmetric lines-of-sight. Each line-of-sight of the instrument uses an elliptical crystal to acquire both the 2p and 3p Ti absorption lines on a single, time gated channel of the instrument. The 2 axial symmetric lines-of-sight allow the 2p and 3p absorption features to be measured through the same point in space using both channels of the instrument. The spatially dependent material temperature can be inferred by observing the 2p and the 3p Ti absorption features. The data are recorded on a two strip framing camera with each strip collecting data from a single line-of-sight. The design is compatible for use at both the OMEGA laser and the National Ignition Facility. The spectrometer is intended to measure the material temperature behind a Marshak wave in a radiatively driven SiO{sub 2} foam with a Ti foam tracer. In this configuration, a broad band CsI backlighter will be used for a source and the Ti absorption spectrum measured.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗