SEARCH · Engineering Papers

Results for “RuN”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 721 records · Page 40

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Accelerating Machine Learning Inference with GPUs in ProtoDUNE Data Processing

Abstract We study the performance of a cloud-based GPU-accelerated inference server to speed up event reconstruction in neutrino data batch jobs. Using detector data from the ProtoDUNE experiment and employing the standard DUNE grid job submission tools, we attempt to reprocess the data by running several thousand concurrent grid jobs, a rate we expect to be typical of current and future neutrino physics experiments. We process most of the dataset with the GPU version of our processing algorithm and the remainder with the CPU version for timing comparisons. We find that a 100-GPU cloud-based server is able to easily meet the processing demand, and that using the GPU version of the event processing algorithm is two times faster than processing these data with the CPU version when comparing to the newest CPUs in our sample. The amount of data transferred to the inference server during the GPU runs can overwhelm even the highest-bandwidth network switches, however, unless care is taken to observe network facility limits or otherwise distribute the jobs to multiple sites. We discuss the lessons learned from this processing campaign and several avenues for future improvements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Lax-Oleinik-Type Formulas and Efficient Algorithms for Certain High-Dimensional Optimal Control Problems

Two of the main challenges in optimal control are solving problems with state-dependent running costs and developing efficient numerical solvers that are computationally tractable in high dimension. In this paper, we provide analytical solutions to certain optimal control problems whose running cost depends on the state variable and with constraints on the control. We also provide Lax-Oleinik-type representation formulas for the corresponding Hamilton-Jacobi partial differential equations with state-dependent Hamiltonians. Additionally, we present an efficient, grid-free numerical solver based on our representation formulas, which is shown to scale linearly with the state dimension, and thus, to overcome the curse of dimensionality. Using existing optimization methods and the min-plus technique, we extend our numerical solvers to address more general classes of convex and nonconvex initial costs. We demonstrate the capabilities of our numerical solvers using implementations on a central processing unit (CPU) and a field-programmable gate array (FPGA). In several cases, our FPGA implementation obtains over a 10 times speedup compared to the CPU, which demonstrates the promising performance boosts FPGAs can achieve. Furthermore, our numerical results show that our solvers have the potential to serve as a building block for solving broader classes of high-dimensional optimal control problems in real-time.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

DISCOVR strain pipeline screening – Part II: Winter and summer season areal productivities and biomass compositional shifts in climate-simulation photobioreactor cultures

Assessing the seasonal biomass productivity and compositional shift dynamics under simulated outdoor culture conditions of the top 22 algae strains selected during Tier I flask screening is an important step in the further prioritization of strains with regard to outdoor pond cultivation. These top 22 strains were subjected to Tier II testing in the PNNL Laboratory Environmental Algae Pond Simulator (LEAPS) photobioreactors, simulating light and temperature conditions of 20 cm deep outdoor ponds during the Arizona winter and summer season. All strains were grown in two consecutive nutrient-replete batch culture experiments at their particular optimal medium salinity to quantify their respective seasonal linear-phase areal biomass productivities. To determine biomass compositional shifts in response to nutrient-depletion, the LEAPS cultures were allowed to enter a 9-day nutrient depletion phase at the end of the second batch run. The following strains were evaluated in winter-season climate-simulated cultures and are listed in the order from highest (7.9 g m -2 day -1 ) to lowest (2.3 g m -2 day -1 ) areal N-replete biomass productivity: Monoraphidium minutum 26B-AM, Tetraselmis striata LANL1001, Chlorella vulgaris LRB AZ-1201, Micractinium reisseri NREL14-F2, Monoraphidium sp. MONOR1, Chlorella vulgaris NREL4-C12, Scenedesmus obliquus UTEX393, Scenedesmus acutus LRB-AP-0401, Nannochloropsis oceanica CCAP849/10, and Stichococcus minutus CCALA727. The following strains were evaluated in summer-season climate-simulated cultures and are listed in the order from highest (31.8 g m -2 day -1 ) to lowest (14.2 g m -2 day -1 ) areal N-replete biomass productivity: Picochlorum renovo NREL39-A8, Scenedesmus obliquus UTEX393, Porphyridium cruentum CCMP675, Picochlorum celeri TG2-WT-CSM/EMRE, Chlorella sorokiniana DOE1116, Stichococcus minor CCMP819, Picochlorum oklahomensis CCMP2329, Chlorella sorokiniana DOE1412 (UTEXB3016), Scenedesmus rubescens NREL46B-D3, Picochlorum soloecismus DOE101, Tetraselmis striata LANL1001, Scenedesmus obliquus DOE 0152.z, and Agmenellum quadruplicatum UTEX2268. All LEAPS cultures experienced a significant reduction in areal biomass productivity in response to nutrient-depletion, from 7-16% in the winter season simulation and up to 1-60% in the summer season simulation. Finally, for 10 of the strains tested, the carbohydrate content more than doubled upon nutrient depletion, and for 9 strains, the lipid content increased by over 50% of the initial content.

09 BIOMASS FUELS↗

DOE OSTI2020

Reliability of emergency and standby diesel generators: Impact on energy resiliency solutions

Emergency diesel generators are the most common form of backup power for critical loads when the grid fails and are most often deployed as stand-alone generators (<2000 kW) tied to individual buildings for hospitals, emergency services, military bases, ports, airports, industries, and commercial facilities. Understanding the finite reliability of emergency diesel generators during continuous operation is critical for energy planners, managers, and end-users. A new analysis of two large non-public emergency diesel generator operational data sets shows that commonly used reliability metrics are inadequate to predict the performance during a grid outage. Additional metrics addressing the likelihood of an emergency diesel generator to be available at the time an outage starts, successfully start and carry load, and constantly run for the duration of the outage are needed to predict the continuous run-time reliability of emergency diesel generators during a grid outage. This is the first analysis of the reliability of modern generators that follow standard maintenance protocols commonly used for backup power and these new metrics are shown to be dependent on the level of emergency diesel generator maintenance with values that can vary by more than an order of magnitude. Even using the new metrics for well-maintained emergency diesel generators, the analysis shows that the single emergency diesel generators configurations are only 80% likely to provide power for the duration of a two-week grid outage.

ENERGY PLANNING, POLICY, AND ECONOMY,POWER TRANSMI↗

DOE OSTI2022

Effect of PVD-coated chromium on the subcooled flow boiling performance of nuclear reactor cladding materials

Here we elucidate the separate effect of a thin Cr coating deposited by physical vapor deposition (PVD) on the subcooled flow boiling performance of zircaloy-4. First, we run flow boiling experiments on prototypical zircaloy-4 surfaces mimicking the scratch pattern and surface roughness of nuclear reactor claddings. Then, we PVD-coat a 0.3 µm thick chromium layer on the same exact surface and repeat the same flow boiling investigations. All experiments are run using deionized water at atmospheric pressure, flowing on a 1 × 3 cm 2 rectangular cross section channel at a rate of 1000 kg/m 2 /s and a subcooling of 10 K. We measure the average temperature of the boiling surface at increasing surface heat fluxes, covering a wide range of heat transfer regimes, from single-phase forced convection to the boiling crisis. We also record high-speed videos of the boiling process, which we postprocess to measure bubble nucleation site density, growth time, departure diameter and frequency. The surface analysis reveals that, while the chromium coating does not seem change the surface roughness and morphology, it improves surface wettability. However, it decreases the critical heat flux. The chromium coating causes an increase of nucleation temperature, bubble departure diameter and growth time, and a reduction of the nucleation site density. The concurrence of these observations indicates that a size reduction of the nucleation sites, conformally covered by the chromium coating, may be the cause of the boiling performance deterioration. We confirm this hypothesis repeating the same analysis on a FeCrAl sample prepared and tested using the same protocol as the zircaloy-4 sample, but with a different initial surface texture.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2021

Search for low-energy neutrinos from astrophysical sources with Borexino

We report on searches for neutrinos and antineutrinos from astrophysical sources performed with the Borexino detector at the Laboratori Nazionali del Gran Sasso in Italy. Electron antineutrinos ( $\bar{ν}_e$ ) are detected in an organic liquid scintillator through the inverse β -decay reaction. In the present work we set model-independent upper limits in the energy range 1.8-16.8 MeV on neutrino fluxes from unknown sources that improve our previous results, on average, by a factor 2.5. Using the same data set, we first obtain experimental constraints on the diffuse supernova $\bar{ν}_e$ fluxes in the previously unexplored region below 8 MeV. A search for $\bar{ν}_e$ in the solar neutrino flux is also presented: the presence of $\bar{ν}_e$ would be a manifestation of a non-zero anomalous magnetic moment of the neutrino, making possible its conversion to antineutrinos in the strong magnetic field of the Sun. We obtain a limit for a solar $\bar{ν}_e$ flux of 384 cm -2 s -1 (90% C.L.), assuming an undistorted solar 8 B neutrinos energy spectrum, that corresponds to a transition probability pνe→$\bar{ν}_e$< 7.2 ×10-5 (90% C.L.) for E $\bar{ν}_e$ > 1.8 MeV.At lower energies, by investigating the spectral shape of elastic scattering events, we obtain a new limit on solar 7 Be- ν e conversion into $\bar{ν}_e$ of pνe→$\bar{ν}_e$< 0.14 (90% C.L.) at 0.862 MeV. Last, we investigate solar flares as possible neutrino sources and obtain the strongest up-to-date limits on the fluence of neutrinos of all flavor neutrino below 3–7 MeV. Assuming the neutrino flux to be proportional to the flare’s intensity, we exclude an intense solar flare as the cause of the observed excess of events in run 117 of the Cl-Ar Homestake experiment.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2022

Rare earth element separations by high-speed counter-current chromatography

We report that following the initial development of High-Speed Counter-Current Chromatography (HSCCC) in the 1960s, several studies have explored its applicability in the separation of rare earth elements (REEs). More recently, however, HSCCC publications have transitioned towards the separation of natural products or pharmaceuticals, leaving the application for REEs largely unexplored from a practical standpoint. Herein, we expand upon prior work in this field by evaluating the suitability of HSCCC to separation of a subset of non-radioactive REEs (Nd, Sm, Eu, Tb, and Y) at 10 -4 mol levels using di-(2-ethylhexyl)phosphoric acid (HDEHP) in n-heptane as the stationary phase and hydrochloric acid as the mobile phase. First, the effect of flow rate on the stationary phase volume retention ratio and resolution of Nd/Sm/Eu subgroup was evaluated followed by optimization of step-gradient elution profiles resulting in additional recovery of Tb and Y within a seven-hour window. The five REEs were separated at the baseline resolution level or above. Elution profiles obtained from multiple runs across two independently operated columns and across independent runs were cross analyzed. Reproducibility in elution profiles point to future applications in radioelement separation chemistry, where both chemical and radiochemical purity are of importance.

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

Physics-Informed Machine Learning Models for Predicting the Progress of Reactive-Mixing

This paper presents a physics-informed machine learning (ML) framework to construct reduced-order models (ROMs) for reactive-transport quantities of interest (QoIs) based on high-fidelity numerical simu-lations. QoIs include species decay, product yield, and degree of mixing. The ROMs for QoIs are applied to quantify and understand how the chemical species evolve over time. First, high-resolution datasets for constructing ROMs are generated by solving anisotropic reaction-di?usion equations using a non-negative finite element formulation for di?erent input parameters. The reactive-mixing model input parameters are: time-scale associated with flipping of velocity, spatial-scale controlling small/large vortex structures of velocity, perturbation parameter of the vortex-based velocity, anisotropic dispersion strength/contrast, and molecular diffusion. Second, random forests, F-test, and mutual information criterion are used to evaluate the importance of model inputs/features with respect to QoIs. We observed that anisotropic dispersion strength/contrast is the most important feature and time-scale associated with flipping of velocity is the least important feature. Third, Support Vector Machines (SVM) and Support Vector Regression (SVR) are used to construct ROMs based on the model inputs. The constructed SVR-ROMs are then used to predict scaling of QoIs. We also present estimates and inequalities on the QoIs, which inform that the species decay, mix, and produce in an exponential fashion. These inequalities also inform that a radial basis function is the most suitable kernel for the SVM/SVR models for QoIs. It is observed that R2-score for SVR-ROMs on unseen data is greater than 0.9, implying that the SVR-ROMs are able to predict the reaction-diffusion system state reasonably well. Finally, in terms of the computational cost, the proposed SVM-ROMs are O(107) times faster than running a high-fidelity finite element simulation for evaluating QoIs. This makes the proposed ML-based ROMs attractive for reactive-transport sensing and real-time monitoring applications as they are significantly faster yet reasonably accurate.

Mudunuru, Maruti K.↗

DOE OSTI2022

ExaCA: A performance portable exascale cellular automata application for alloy solidification modeling

Modeling the as-solidified grain structures that form during alloy processing is a critical component in understanding process-property relationships, particularly for additive manufacturing (AM) where grain structure is very sensitive to processing conditions. While cellular automata (CA)-based models have proven able to predict aspects of microstructure for several alloys and AM process conditions, long run times and large resource sets required limit the utility and the problem size to which existing CA models can be applied. As part of the ExaAM project, an initiative within the Exascale Computing Project (ECP) to develop, test, and optimize an exascale-capable coupled and self-consistent model of AM parts, we developed ExaCA (https://github.com/LLNL/ExaCA) for the liquid–solid phase transformation in the wake of AM melt pools. The CA-based code is parallelized using MPI and the Kokkos programming model, the latter enabling simulation on both CPUs and GPUs within a single-source implementation. Here, we detail the steps taken to transform a baseline, MPI-based CA code into one that is performant on CPUs and GPUs. Performance testing of ExaCA on Summit (a pre-exascale machine at Oak Ridge National Laboratory) was used to quantify CPU–GPU speedup comparing with equal numbers of nodes. Testing showed comparable CPU performance to the MPI-only CA code and a 5-20x speedup when running AM-based test problems using GPUs. The improved performance of CA through GPU utilization and the performance portable nature of ExaCA will enable accurate part-scale modeling by harnessing the power of current and future generations of high performance computing resources. Future work will include improving the strong scaling of ExaCA on GPUs by reducing load imbalance associated with the locality of the problem, and continuing performance optimization across exascale hardware.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Woven ceramic matrix composite surrogate model based on physics-informed recurrent neural network

A recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear constitutive behavior of woven ceramic matrix composites (CMCs) driven by matrix damage at multiple length scales. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the multiscale generalized method of cells (MSGMC) approach coupled with a matrix damage model. This coupling permits simulating the nonlinear behavior of woven CMCs based on constituent response at the micro-, meso-, and macroscales. The multiscale repeating unit cell is loaded under non-monotonic conditions including multiple load / unload cycles and tension / compression. The fiber volume fraction as well as the intra- and intertow void volume fractions are also varied in the generation of training data. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input strain sequence and fiber and void volume fractions, the resulting stress versus strain response while satisfying physical constraints such as positive semi-definiteness of the tangent stiffness matrix and linear elastic unloading. Further, the trained surrogate model effectively matches the stress versus strain response and successfully predicts the tangent modulus throughout the loading regime. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex architectures, nonlinear multiaxial material response, and under non-monotonic loading conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

SOMAFOAM: An OpenFOAM based solver for continuum simulations of low-temperature plasmas

Here, we report the development of SOMAFOAM, a finite volume framework for performing continuum simulations of low-temperature plasmas. The primary goal of this work is to discuss the features of SOMAFOAM along with representative results provided as examples for a range of operating conditions and geometries. This includes plasma and plasma–dielectric systems operating in direct current, radio frequency, and microwave regimes from pressures as low as 100 mTorr to atmospheric pressure. The code has several useful features including the ability to run massively parallel simulations using arbitrary geometries, structured/unstructured meshes, choice of models such as drift–diffusion/full-momentum at runtime, and species-dependent timesteps to name a few. The verification/validation studies presented include comparison with previously published continuum simulations (low-pressure direct current and radio frequency plasma), with experiments (Gaseous Electronics Conference Reference Cell and microwave microplasma ignited in a split ring resonator), and previously published kinetic simulations (low-pressure radio frequency plasma). Other examples provided include a direct current atmospheric pressure microplasma bounded by dielectric sidewalls and a helium–nitrogen plasma ignited using a needle electrode facing a dielectric. The performance of the code is also discussed with serial and distributed memory parallel runs demonstrated up to 512 cores. The design and implementation of the code in a modular object-oriented framework allows for easy extension and seamless coupling with other codes and can be expected to play an important role in both academia and industry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

TUMME: Tsinghua University Minnesota Master Equation program

We report that TUMME is a program for assembling and solving master equations for gas-phase chemical kinetics based on chemically significant eigenmodes. TUMME has interfaces to the Gaussian, Polyrate, and/or MSTor output files that allow the master equation code to obtain the microcanonical flux coefficients needed for the coefficient matrix of the master equation. The flux coefficients for reactions with barriers can be calculated by multi-structural variational transition state theory with small-curvature tunneling (MS-VTST/SCT) or by simpler approximations to this such as conventional transition state theory without tunneling (also called RRKM theory). The flux coefficients for barrierless reactions are provided by a hard-sphere model. TUMME is written in double precision with Python 3; quadruple and octuple precision are also available for some subtasks in C++. The Python code can run in serial or parallel (MP or MPI), and the C++ code can run on a single processor or on multiple processors with OpenMP.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

$\mathrm{DRAGON}$: A multi-GPU orbital-free density functional theory molecular dynamics simulation package for modeling of warm dense matter

As progress in electronic structure theoretical methods is made, ab initio molecular dynamics (MD) based on orbital-free density functional theory (OF-DFT) is becoming increasingly more successful at substituting the traditional, very accurate but computationally costly Kohn–Sham (KS) approach for simulations of matter at the challenging warm dense matter (WDM) regime. However, despite the significant cost alleviation of eliminating the dependence on the KS orbitals, OF-DFT MD runs require ~10 2 to 10 3 CPU cores running for days, or even weeks, for simulations of systems comprised of 10 2 to 10 3 atoms, depending on thermodynamic conditions. Here we present DRAGON, a multi-GPU OF-DFT MD code for fast and efficient simulations of WDM. With a relatively small allocation of resources (4 to 8 GPU devices) it can provide an order of magnitude speedup for simulations containing $\mathscr{O}$(10 4 ) atoms and target systems composed of $\mathscr{O}$(10 5 ) atoms at conditions within the WDM regime, which is currently outside the capabilities of CPU codes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

$\mathrm{D}$ark$\mathrm{F}$lux: A new tool to analyze indirect-detection spectra of next-generation dark matter models

Here we present DarkFlux, a software tool designed to analyze indirect-detection signatures for next-generation models of dark matter (DM) with multiple annihilation channels. Version 1.0 of this tool accepts user-generated models with 2→2 tree-level dark matter annihilation to pairs of Standard Model (SM) particles and analyzes DM annihilation to γ rays. The tool consists of three modules, which can be run in a loop in order to scan over DM mass if desired: The annihilation fraction module calls an internal installation of MadDM, a dark matter phenomenology plugin for the Monte Carlo event generator MadGraph5 _ aMC@NLO, to compute the thermally averaged cross section $\langle σv\rangle_i$ for each annihilation channel $χχ$ ($\bar{χ}$,$χ^†$)→ i $\∈{SM, SM}. . The module then computes the fractional annihilation rate (annihilation fraction) into each channel. The flux module combines the flux spectrum from each annihilation channel, weighted by the appropriate annihilation fractions, to compute the total flux at Earth due to DM annihilation. In DarkFlux v1.0, this module specifically computes the γ-ray flux for each channel using the publicly available PPPC4DMID tables. The analysis module compares the total flux to observational data and computes the upper limit at 95% confidence level (CL) on the total thermally averaged DM annihilation cross section. In DarkFlux v1.0, this module compares the total γ-ray flux to a joint-likelihood analysis of fifteen dwarf spheroidal galaxies (dSphs) analyzed by the Fermi-LAT collaboration. DarkFlux v1.0 automatically provides data tables and can plot the output of these three modules. In this manual, we briefly motivate this indirect-detection computer tool and review the essential DM physics. We then describe the several modules of DarkFlux in greater detail. Finally, we show how to install and run DarkFlux and provide two worked examples demonstrating its capabilities.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

In situ deformation of antigorite-olivine two-phase mixtures: Implications for dynamics and seismic anisotropy in the mantle wedge

Water released from hydrous minerals in subducting slabs reacts with the overlying plate, resulting in widespread serpentinization in the mantle wedge. Deformation of serpentinized peridotites has been invoked to explain forearc seismic anisotropy, yet studies of the mechanical properties and deformation behaviors of serpentine-bearing multiphase aggregates remain limited. Here we deformed olivine-antigorite mixtures containing 70, 50, and 20 vol.% of antigorite at 2.5 – 7.6 GPa, 673 K and strain rates of ∼10 –5 –10 –4 s –1 . Elasto-viscoplastic self-consistent simulations, constrained by synchrotron X-ray diffraction (XRD) data, were used to estimate lattice strain, stress–strain partitioning, crystallographic preferred orientations (CPO) development, and aggregate strength. Selected run products were also analyzed by electron backscatter diffraction for comparison with the CPO results obtained from XRD experiments. We found olivine transitions from A- or B-type to C-type when antigorite fraction drops to 20 vol.%, coinciding with a microstructural change from interconnected weak layers to a load-bearing framework (LBF). An additional run on a sample Atg50/Ol50 with preexisting microstructures suggested the formation of LBF was promoted by these microstructures, although the preexisting antigorite CPO has been overprinted at 20.8 % strain and could be erased completely by subsequent deformation in nature. Estimated viscosity of the two-phase mixtures suggests that low-degree serpentinization (≤20 %) in the mantle wedge may increase the strength of olivine-rich peridotite and hinder slab-mantle decoupling, whereas high-degree serpentinization (≥50–70 %) weakens the peridotite and favors decoupling if sufficient viscosity contrast (>10) develops. Seismic anisotropy shows a nonlinear dependence on antigorite fraction: antigorite CPO governs the anisotropy of the mixtures with ≥50 vol.% antigorite, whereas olivine CPO dominates at low fractions (∼20 vol.%). The presence of pre-existing microstructures reduces seismic anisotropy of the deformed mixtures, however the persistence of pre-existing CPO in actively subducting slabs remains uncertain, making their significance over geological timescales questionable.

Crystallographic preferred orientation↗

DOE OSTI2021

Use of a Lewis acid, a Brønsted acid, and their binary mixtures for the hydrothermal liquefaction of lignocellulose

The main objective of the present study is to investigate the effect of a Lewis acid, Brønsted acid, and their combined use on the hydrothermal liquefaction of lignocellulosic biomass. Hydrothermal liquefaction of teak wood was conducted at 250, 300 and 350 °C for 15, 30 and 60 min. Hydrothermal liquefaction of teak wood was carried out at 300 °C for 30 min (the best optimum conditions) without and with the use of Mg(ClO 4 ) 2 , HClO 4 , and HClO 4 /Mg(ClO 4 ) 2 at various loadings (2–10 mmol/15 g wood). The highest bio-oil yield was obtained with the non-catalytic run. All tested catalysts have negative effect on bio-oil yields. The bio-oil yields generally decreased with increasing the catalyst loadings. The deoxygenation degree in bio-oils changed depending on the type of catalyst and loading. A high degree of de-oxygenation took place with Mg(ClO 4 ) 2 catalysts. An increased catalyst loading led to decreased aromatic contents of bio-oils catalyzed by either Mg(ClO 4 ) 2 or HClO 4 . The use of a catalyst increased total naphtha fractions in bio-oils. The highest heating value of the bio-oil was estimated to be approximately 30 MJ/kg. Gas chromatography–mass spectrometry analysis revealed that the bio-oils from the non-catalytic and catalytic runs contained aldehydes, ketones, phenols, acids, esters and alcohols. Furthermore, the relative yields of the oxygenated compounds were affected by catalyst type.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Adaptive elasticity policies for staging-based in situ visualization

In situ processing aims to alleviate the growing gap between computation and I/O capabilities by performing data processing close to the data source. In situ processing is widely used to process data generated by multiple data sources, including observation data from edge devices or scientific observational facilities and the simulation data generated by scientific computation on a high-performance computing (HPC) platform. For a scientific workflow that is run on an HPC platform and composed of a simulation program and an in situ data analytics or visualization (abbreviated as ana/vis) task, there is an implicit assumption that the computing resources assigned to the workflow keep static during the workflow execution. However, with the converging trend between the HPC and cloud computing platform, running the in situ ana/vis task in an elastic way is promising to decrease its overhead and improve its resource utilization rate. Resource elasticity represents the ability to change resource configurations such as the number of computing nodes/processes during workflow execution. An elastic job may dynamically adjust resource configurations; it may use a few resources at the beginning and more resources toward the end of the job when interesting data appear. However, it is hard to predict a priori how many computing nodes/processes need to be added/removed during the workflow execution to adapt to changing workflow needs. How to efficiently guide elasticity operations, such as growing or shrinking the number of processes used for in situ analysis during workflow execution, is an open-ended research question. In this article, we present adaptive elasticity policies that adopt workflow runtime information collected during workflow execution to predict how to trigger the addition/removal of processes in order to minimize in situ processing overhead. Taking in situ visualization tasks as an example, we integrate the presented elasticity policies into a staging-based elastic workflow and evaluate its efficiency in multiple elasticity scenarios. Compared with the situation without elasticity or with a static elasticity policy that uses a fixed number of processes for each rescaling operation, the adaptive elasticity policy can save overhead in finding a proper resource configuration and improve resource utilization efficiency. Furthermore, one experiment illustrates that the adaptive elasticity policy saves 41% of core-hours compared with the situation without the resource elasticity.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Assessing Impacts on Pressure Stabilization and Leasing Acreage for CO 2 Storage Utilizing Oil Migration Concepts

Favorable geological storage for CO 2 has long been pictured as large anticlines with thick sandstones, similar to oil reservoirs in the petroleum system. Unlike oil, however, stored CO 2 does not need to be recoverable, which raises the possibility of using dissolution and residual trapping to augment the capacity of buoyant traps and tap more of the bulk rock volume. The work presented builds on that idea, asking the following question: If we inject CO 2 down to a syncline – analogous to the carrier bed in the petroleum system – how would this injection mechanism impact storage capacity and plume shape, migration, and stabilization? To address this question, we built a reservoir model, based on seismic interpretation of Middle Miocene strata, offshore Galveston, Texas. 3-D seismic and well logs were used to characterize key intervals. Reservoirs chosen for modeling are progradational-aggradational sands with mud intercalation. They have a higher degree of heterogeneity than the more conventional reservoirs commonly targeted for CO 2 storage. Modeling investigated how far the CO 2 plume would migrate under two scenarios: (1) injecting CO 2 at the base of the salt withdrawal basin (syncline scenario) and (2) injecting CO 2 at the base of the structural closure, similar to a common injection well location for EOR purposes (base scenario). For each scenario, we separately simulated injection of 30 MT of CO 2 and 60 MT of CO 2 continuously for 30 years and observe the plume and pressure evolution for 100 years after the injection stops. The simulation shows that injecting the CO 2 into a syncline limits the vertical migration of CO 2 , thus making synclinal injection more secure. In the syncline scenario, the geological layer around the injection point is more heterogeneous than the layer in the base scenario; thus, the CO 2 tends to migrate laterally. Additionally, in the syncline scenario, the plume does not even reach the upper part of the anticline, allowing us to safely store an additional amount of CO 2 into the reservoir. Furthermore, the simulation also shows that in the syncline scenario, the times needed for the reservoir to reach its stabilized pressure after the end of injections are faster. To summarize, CO 2 injection at the base of a syncline could provide additional storage, increase the safety of the project from the limited vertical plume migration, and expedite plume stabilization, which could result in the decrease of monitoring frequency as the project runs, thus lowering the operating cost of the project in the long run.

58 GEOSCIENCES↗