Engineering Papers⌕ Search

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.

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

At least 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A comparative study of machine learning models for predicting the state of reactive mixing

Mixing phenomena are important mechanisms controlling flow, species transport, and reaction processes in fluids and porous media. Accurate predictions of reactive mixing are critical for many Earth and environmental science problems such as contaminant fate and remediation, macroalgae growth, and plankton biomass evolution. Here, to investigate the evolution of mixing dynamics under different scenarios (e.g., anisotropy, fluctuating velocity fields), a finite-element-based numerical model was built to solve the fast, irreversible bimolecular reaction-diffusion equations to simulate a range of reactive-mixing scenarios. A total of 2,315 simulations were performed using different sets of model input parameters comprising various spatial scales of vortex structures in the velocity field, time-scales associated with velocity oscillations, the perturbation parameter for the vortex-based velocity, anisotropic dispersion contrast (i.e., ratio of longitudinal-to-transverse dispersion), and molecular diffusion. The outputs comprised concentration profiles of reactants and products. The inputs to and outputs from these simulations were concatenated into feature and label matrices, respectively, to train 20 different machine learning (ML) models intended to emulate system behavior. These 20 ML emulators, based on linear methods, Bayesian methods, ensemble learning methods, and multilayer perceptrons (MLPs), were trained to classify the state of mixing and predict three quantities of interest (QoIs) characterizing species production, decay (i.e., average concentration, square of average concentration), and degree of mixing (i.e., variances of species concentration). Unsurprisingly, linear classifiers and regressors failed to reproduce the QoIs; however, ensemble methods (classifiers and regressors) and the MLP model accurately classified the state of reactive mixing and the QoIs. Among ensemble methods, random forest and decision-tree-based AdaBoost faithfully predicted the QoIs. At run time, trained ML emulators produced results times faster than the finite-element simulations. Due to their low computational expense and high accuracy, ensemble and MLP models are excellent emulators for these numerical simulations and great utilities in uncertainty quantification exercises, which can require 1,000s of forward model runs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Historic and modern approaches for discovery of abandoned wells for methane emissions mitigation in Oil Creek State Park, Pennsylvania

Hundreds of oil wells were drilled along Oil Creek in Pennsylvania in the mid-16 1800s, birthing the modern oil industry. No longer in operation, many wells are now classified as abandoned, and, due to their age, their locations are either unknown or inaccurately recorded. These historic well sites present environmental, safety, and economic concerns in the form of possible methane leaks and physical hazards. Airborne magnetic and LiDAR surveys were conducted in the Pioneer Run watershed in Oil Creek State Park to find abandoned wells in a historically significant but physically challenging location. Wells were drilled in this area prior to modern geolocation and legal documentation. Although a large number of old wells were abandoned summarily without remediation of the site, much of the land area within Oil Creek State Park is now covered in trees and dense underbrush, which can obscure wellheads. The thick vegetation and steep terrain limited the possibility of ground-based surveys to easily find well sites for methane emissions studies. The data from remote sensing surveys were used to corroborate potential well locations from historic maps and photographs. Potential well sites were verified in a ground-based field survey and monitored for methane emissions. Two historic photographs documenting oil activity in the late 1800s were georeferenced using a combination of magnetic and LiDAR data. LiDAR data, which were more useful in georeferencing and in field verification, identified 290 field locations in the Pioneer Run watershed, 86% of which were possible well sites. Sixty-two percent of the ground-verified wells remained unplugged and comprised the majority of leaking wells. The mean methane emissions factor for unplugged wells was 0.027 ± 0.099 kg/day, lower than other Appalachian Basin methane emissions estimates. LiDAR was used for the first time, in combination with an airborne magnetic survey, to reveal underground oil industry features and inform well identification and remediation efforts in difficult-to-navigate regions. In the oldest oil fields, where well casing has been removed or wood conductor casing was installed, historic photographs provide additional lines of evidence for oil wells where ground disturbances have concealed surface features. Identification of well sites is necessary for mitigation efforts, as unplugged wells emit methane, a potent greenhouse gas.

02 PETROLEUM↗
DOE OSTI2025

Binding and Translocation of Substrate Allosterically Promotes Functional Interactions Within the AlkB–AlkG Electron Transfer Complex

The alkane monooxygenase AlkB and rubredoxin AlkG form an electron transfer complex that hydroxylates terminal alkanes to produce alcohols. The recent cryoEM study of Fontimonas thermophila AlkB-AlkG complex revealed its architecture, including a dodecane (D12) substrate at the active site. However, FtAlkBG molecular mechanism of action of remains unknown. Here, in this study, we examined its dynamics and interactions by multiscale computations, including molecular dynamics simulations, elastic network models, and QM/MM of the oxygen activation mechanism at the AlkB catalytic site. D12 maintained stable interactions within the catalytic site during two MD runs, coordinated by hydrophobic residues L263-L264, I267, I133. A third extended run revealed that D12 could translocate to a membrane-exposed site near S49/F46 along a hydrophobic channel gated by I54. During this translocation, D12 was temporarily stabilized at intermediate sites IS1 (lined by I27/L30-G31/G50/L53-I54/P59/S124/A127-V128) and IS2 (I33-G34/L37/L45-F46/S49) before nearly exiting the protein, and diffused back to the active site, assisted by L30. Substrate binding and translocation across those intermediate sites affects the coupling between the iron centers in AlkBG, and interfacial interactions between AlkB-AlkG. The channel was further connected to the cytosol, near two surface-exposed arginines, potentially allowing for O 2 passage. The allosteric effects between D12 putative entry site, catalytic site and AlkB-AlkG interface were analyzed by ENM-based methods which confirmed the cooperative perturbation-responses and strongly correlated movements of residues belonging to those distal regions. Our study provides new mechanistic insights into key sites and their interactions that could be targeted for developing AlkB-variants with desirable alkane conversion functions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Synthesis and structural characterization of the new Zintl phases Ba 3 Cd 2 P 4 and Ba 2 Cd 2 P 3 . Rare example of small gap semiconducting behavior with negative thermopower within the range 300 K-700 K

The new Zintl phases Ba 3 Cd 2 P 4 and Ba 2 Cd 2 P 3 have been synthesized using Pb flux, which allowed for the growth of 4-5 mm large crystals. The structures were determined utilizing single-crystal X-ray diffraction methods. Both compounds crystalize in the monoclinic crystal system (space group C2/m (No. 12)) and their structures are closely related. The structure of Ba 3 Cd 2 P 4 can be seen as being comprised of divalent Ba atoms and conjoined CdP 4 tetrahedra in the form of [Cd 2 P 4 ] 6- layers. Within the layers, homoatomic P–P bonds are present, which if cleaved, leave two infinite [CdP 3 ] 7- chains running along the crystallographic b-axis. The other structure, that of Ba 2 Cd 2 P 3 , can be rationalized as also having divalent Ba atoms and conjoined CdP4 tetrahedra in the form of [Cd 2 P 3 ] 6- layers. These layers, again, can be visualized as chains that run down the crystallographic b-axis, which are further connected by P-P dimers. Electronic band structure calculations show that each structure has an optimal number of valence electrons, and therefore conform to the Zintl-Klemm concept. Accordingly, the two compounds can be considered small band gap semiconductors, with band gaps of ca. 0.1 eV and 0.6 eV for Ba3Cd2P4 and Ba 2 Cd 2 P 3 , respectively. Electrical resistivity measurements show that Ba3Cd2P4 displays a large resistivity value at room temperature and an experimental band gap of ca. 0.05 eV, which fits reasonably well with the theoretical predictions. Thermopower measurements show that throughout the temperature range 300 K-700 K, Ba 3 Cd 2 P 4 displays a negative Seebeck coefficient. Here, the extremum value of -84 μV is reached at 630 K, suggestive of an n-type semiconductor, a rarity among Zintl phases.

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

A carbon molecular sieve membrane-based reactive separation process for pre-combustion CO 2 capture

In this study, we discuss a hybrid system combining a membrane reactor (MR) and an adsorptive reactor (AR), with the MR's reject stream serving as the AR's feed. We apply this system for the water gas shift (WGS) reaction for H 2 generation and simultaneous CO 2 capture in the context of the Integrated Gas Combined Cycle (IGCC) process for power generation from coal and biomass. This MR-AR system attains a high conversion exceeding equilibrium, produces a pure H 2 product for power generation, and delivers a high-pressure CO 2 stream ready for sequestration. Specifically, in our study we use carbon molecular sieve membranes (CMSMs) and a commercial sour-shift WGS catalyst. Lab experiments were carried-out to determine the membrane characteristics, and the MR performance under IGCC-relevant conditions, i.e., for temperatures up to 250 °C and pressures up to 25 bar, employing a model coal gasifier syngas. The CMSM and the catalyst have displayed robust and stable performance during a long-term run (~750 h of syngas exposure). We evaluated the MR-AR system in multi-cycle runs and it has demonstrated superior performance to that of a conventional packed-bed reactor, producing a high-purity H 2 product directly useable in a turbine for power generation. We conclude from the study, that the CMSM-based MR-AR system is a good candidate technology for environmentally-benign power generation. We are currently constructing a pilot-scale system for field demonstration of the technology.

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

Dense Ru single-atoms integrated with sulfoacids for cellulose valorization to isosorbide

Metal-acid bifunctional catalysts are the unity of two opposites (metal-acid repelling) for cellulosic biomass valorization to isosorbide. It is challenging to improve the selectivity of metal hydrogenation catalysts and their synergism with acids for catalytic hydrolysis and dehydration to achieve substantial isosorbide. Herein, dense Ru single-atoms (10.1 wt% of Ru SAs) are anchored on sulfoacid-functionalized hollow mesoporous carbon shells, designed by assembling silica and 8-hydroxyquinoline-modified chitosan (HQ-CTS) through in situ Stöber templating strategy before pyrolysis and acid treatment. Based on X-ray absorption fine structure and computational modeling results, the structure of Ru SAs is determined as RuN 4 , which is more selective for a transitional glucose hydrogenation to sorbitol than Ru001 of Ru clusters. A lower-energy barrier of 1.21 (0.72) eV is delivered over RuN 4 (Ru001) for glucose hydrogenation (isomerization). These Ru SAs are integrated with sulfoacids (SO 3 H) but resistent against acids, rendering enhanced isosorbide yield in water as compared to Ru clusters, via a one-pot cascade reaction under harsh conditions (220 °C, 6 MPa H 2 ). The elaborately fabricated dense Ru SAs and sulfoacids, achieved by varying the addition time of HQ-CTS during the in situ Stöber templating process, improve the synergism of glucose hydrogenation with cellulose hydrolysis and sorbitol dehydration. In conclusion, this study provides a new idea for rational design of high-performance metal-acid bifunctional catalysts toward one-pot conversion of cellulose to isosorbide.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Benchmarking the Performance of Neuromorphic and Spiking Neural Network Simulators

Software simulators play a critical role in the development of new algorithms and system architectures in any field of engineering. Neuromorphic computing, which has shown potential in building brain-inspired energy-efficient hardware, suffers a slow-down in the development cycle due to a lack of flexible and easy-to-use simulators of either neuromorphic hardware itself or of spiking neural networks (SNNs), the type of neural network computation executed on most neuromorphic systems. While there are several openly available neuromorphic or SNN simulation packages developed by a variety of research groups, they have mostly targeted computational neuroscience simulations, and only a few have targeted small-scale machine learning tasks with SNNs. Evaluations or comparisons of these simulators have often targeted computational neuroscience-style workloads. In this work, we seek to evaluate the performance of several publicly available SNN simulators with respect to non-computational neuroscience workloads, in terms of speed, flexibility, and scalability. We evaluate the performance of the NEST, Brian2, Brian2GeNN, BindsNET and Nengo packages under a common front-end neuromorphic framework. Our evaluation tasks include a variety of different network architectures and workload types to mimic the computation common in different algorithms, including feed-forward network inference, genetic algorithms, and reservoir computing. We also study the scalability of each of these simulators when running on different computing hardware, from single core CPU workstations to multi-node supercomputers. Our results show that the BindsNET simulator has the best speed and scalability for most of the SNN workloads (sparse, dense, and layered SNN architectures) on a single core CPU. However, when comparing the simulators leveraging the GPU capabilities, Brian2GeNN outperforms the others for these workloads in terms of scalability. NEST performs the best for small sparse networks and is also the most flexible simulator in terms of reconfiguration capability NEST shows a speedup of at least 2x compared to the other packages when running evolutionary algorithms for SNNs. The multi-node and multi-thread capabilities of NEST show at least 2x speedup compared to the rest of the simulators (single core CPU or GPU based simulators) for large and sparse networks. We conclude our work by providing a set of recommendations on the suitability of employing these simulators for different tasks and scales of operations. We also present the characteristics for a future generic ideal SNN simulator for different neuromorphic computing workloads.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Differences in the response of two light guide technologies and two readout technologies after an exchange of liquid argon in the dewar

In this investigation the response to the scintillation light generated by through-going cosmic muons in liquid argon (LAr) was measured by two light guide technologies and two readout technologies after five weeks of running in the TallBo dewar at Fermilab. The response was remeasured after the dewar was drained of LAr, refilled, and then run again for an additional four weeks. After the dewar was refilled, there was clear evidence that the scintillation signal had dropped significantly. The two light guide technologies were developed at Indiana University and MIT/Fermilab. The two readout technologies were boards that passively or actively ganged 12 Hamamatsu MPPCs. Two possible explanations were identified for the degraded signal: the response of the two light guide technologies degraded due to damage caused by thermal cycling, and/or unknown differences in the trace residual Xe contamination in the fills of LAr led to the observed drop in scintillation light. Neither absorption nor quenching by N 2 , O 2 , and H 2 O contamination can account for the degradation. Neither the individual Hamamatsu MPPCs nor the passive/active ganging boards appear to have been affected by the thermal cycling. The path length distributions of the cosmics traversing the dewar appear quite similar in both event samples.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Genetic algorithm optimization of tin–copper graded shielding for improved plutonium safeguards measurements

Plutonium nuclear material safeguards is essential to ensure that states do not divert plutonium to clandestine nuclear weapons programs. Traditional safeguards methods rely on correlated neutron measurements. Detection systems based on organic scintillators, which are sensitive to fast neutrons and photons, are gaining support for use in safeguards applications due to increased capabilities and shorter measurement times compared to the current state-of-the-art. A drawback of organic scintillators in plutonium measurements is data throughput limitations; 241 Am builds up in plutonium samples and emits 60-keV gamma rays with a high specific activity. In this work, a genetic algorithm is used to design graded shielding made of tin and copper to attenuate the 60-keV gamma rays. The shield is designed to limit the attenuation of the higher-energy, fission photons, because fission photons can be used as a signature. The genetic algorithm optimization takes days to run, while a complete search space exploration would take years. The genetic algorithm was run for different levels of 60-keV gamma-ray attenuation, and the thickness of the required shielding layers was determined at each level. Tin–copper shielding was placed in front of an array of organic scintillation detectors for a measurement of approximately 4.5 kg of alpha-phase, weapons-grade plutonium, an AmLi source, and a 252 Cf source. We show that the shield successfully attenuated low-energy, non-fission gamma rays, which significantly improved the organic scintillator performance.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

First 3D vector tracking of helium recoils for fast neutron measurements at Super$\mathrm{KEKB}$

Here we present results from the first deployment of novel, high definition, compact gas Time Projection Chambers (TPCs) with pixel chip readout as part of the BEAST II beam background measurement project at SuperKEKB. The TPCs provide detailed 3D imaging of ionization from neutron-induced nuclear recoils in a helium and carbon dioxide target gas mixture at standard temperature and pressure. We present the TPC performance and the neutron backgrounds observed during the initial stage of collider commissioning. We find excellent electron background rejection, leading to background-free nuclear recoil measurements above 50 keV ee , despite the extreme high-background environment. We measure an angular resolution better than 20° for recoil tracks longer than 1.7 mm, corresponding to an average ionization energy of approximately 100 keV ee . We also obtain the full 3D vector direction of helium recoils by utilizing charge profile measurements along the recoil axis, with a correct head/tail assignment efficiency of approximately 80%. With this performance, we present comparisons between measured and simulated event rates, recoil energy spectra, and directional distributions originating from beam–gas and Touschek beam losses at SuperKEKB. We utilize head/tail recognition to distinguish neutron components traveling with positive radial velocity in the Belle II coordinate system from those traveling in the opposite direction. Finally, we present a novel method of discriminating beam–gas interactions from Touschek beam losses that can eliminate the need for dedicated accelerator runs for background measurements. This method is still statistics-limited. However, future studies should be able to verify this method, which in turn could lead to neutron background analysis runs symbiotic with normal Belle II operation. The capabilities demonstrated here also suggest that high definition recoil imaging in gas TPCs is applicable to low energy, low-background experiments, such as directional dark matter searches.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Experiment-modeling studies comparing energy dissipation in the DIII-D SAS and SAS-VW divertors

Recent DIII-D experiments on Small Angle Slot (SAS) divertors have confirmed that a combination of divertor closure and target shaping can enhance cooling across the divertor target and increase energy dissipation, but with significant dependence on B T (toroidal magnetic field) direction. In these novel divertors, the roles of closure, target shaping, drifts, and scale lengths are all interconnected in optimizing dissipation, with the separatrix electron density n eSEP being the key parameter associated with the level of dissipation/detachment. After modifying the original flat-targeted graphite SAS to include a V shape with a tungsten coating on the outer side of the divertor (SAS-VW), matched series of discharges were run to compare to detailed SOLPS-ITER modeling. Experimentally, when run as designed with the outer strike point at the slot vertex, SAS-VW requires nearly identical n eSEP for detachment as the original SAS, with little difference in dissipation for the new geometry. This is in contrast to (1) earlier modeling predictions that a small change of the SAS geometry to a V shape should enhance dissipation at the same n eSEP for magnetic configurations having better H-mode access (ion B × ∇B drift directed into the divertor), and (2) despite the achievement of significantly higher (2-7x) neutral pressures and compression in the SAS-VW slot. Comparisons of experimental density scans to the most recent SOLPS-ITER modeling with ExB drifts show reasonable agreement for dissipation/detachment onset when using separatrix density as the independent parameter. In order to help understand the discrepancy in modeled vs actual performance for the new configuration, additional measurements varying gas injection location and impurity injection were undertaken. In-slot D 2 gas fueling is more effective (5–22 %) in promoting detachment, in accord with modeling. In-slot impurity injection (N 2 or Ne) can yield 30 % lower core Z eff and 15 % less confinement degradation after detachment compared to main chamber puffing, as well as relatively lower tungsten leakage from the divertor. Modeling can also reproduce the improved detachment seen as the strike point moves inboard of the slot vertex. While we can explain the effects of the most important parameters causing energy dissipation in these slot divertors, it remains that many aspects of their behavior cannot be accurately modeled using state-of-art codes such as SOLPS-ITER. This is of concern for future model-driven designs utilizing similar V-shaped geometries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Advanced two-phase subchannel method via non-linear iteration

A fast-running, robust two-phase flow, sub-channel model is presented based on non-linear solution of the steady-state subchannel fluid flow equations. The drift-flux model solves for conservation of liquid and vapor mass, mixture energy, and axial and transverse mixture momentum as part of an efficient planar marching scheme and nonlinear, nested outer and inner iteration. Here, models based on mechanistic subcooled boiling, two-phase turbulent void mixing, and drift are included. Solution verification and mesh convergence studies were performed for modern GE 10 × 10 fuel geometry and are shown to have excellent convergence behavior. Run time performance for a 50 axial mesh model showed 2.2 seconds on a single CPU core to tightly converge all 3D distributions (flow, void, pressure) for the GE 10 × 10 fuel geometry, supporting its efficient use within the Virtual Environment for Reactor Applications boiling water reactor framework.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Irradiation performance of a U-7Mo in Al-Si matrix dispersion full-size fuel plate assembly

The Korea Atomic Energy Research Institute (KAERI) is leading the Ki-Jang Research Reactor (KJRR) project with the intent to develop a new reactor for medical isotope production and other nuclear research purposes. The KJRR core is designed to use high density fuel system where uranium alloyed with 7 wt% molybdenum (U-7Mo) particles are dispersed in a matrix of aluminum alloyed with 5 wt% silicon (Al-5Si) and clad in aluminum alloy 6061 (Al-6061) to form fuel plates. KAERI developed a fabrication facility to construct KJRR fuel assemblies and partnered with the Idaho National Laboratory (INL) to irradiate a full-size fuel assembly, with 21 total fuel plates, in the Advanced Test Reactor (ATR). Irradiation testing and subsequent Post Irradiation Exam (PIE) campaigns were performed successfully over a multi-year project. Monte Carlo neutronic calculations, coupling with a depletion code, were performed based on ATR’s as-run power history which showed that the highest power plate (plate 20) reached 83.1 % end-of-life (EOL) local burnup based on initial 235 U content. Additionally, finite element thermal modeling was performed based as-run power history which showed a beginning-of-life (BOL) peak local heat flux of 184 W/cm 2 . No anomalous fuel performance was observed during the irradiation and target test conditions were achieved. PIE showed favorable performance of the fuel assembly regarding all important phenomena. This paper describes the KJRR fuel assembly irradiation conditions and PIE data to support the conclusion that it performed well, without evidence of unexpected or problematic fuel performance, within an irradiation test designed to bound the KJRR design environment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Immortal rays: Rethinking random ray neutron transport on GPU architectures

The Random Ray Method (TRRM) is a recently developed adaptation of the Method of Characteristics for neutral particle transport simulations. TRRM has demonstrated excellent performance on 3D nuclear reactor benchmark problems using CPU-based compute systems. When porting to GPU-based systems, however, new performance challenges arise that are unique to processors targeting massive fine-grained parallelism. For smaller problems, or for large problems that are domain decomposed across many computational nodes, the problem size per node has insufficient parallelism to saturate GPU node resources, thus greatly limiting speedup. In this study, we report on a newly developed “immortal ray” variant of TRRM. Here, the immortal ray technique exposes significantly more fine-grained parallelism by fundamentally reformulating the numerical details of ray discretization, resulting in performance tradeoffs with significant overall benefit on GPUs. For very small 2D simulation problems we found the new immortal ray variant allowed for up to a 4.4x speedup when run on a single GPU. For larger 3D simulation problems we found the new variant improved strong scaling by 3x when run on the Summit supercomputer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Machine learning based rate optimization under geologic uncertainty

We propose a novel approach for rate optimization during a waterflood under geologic uncertainty in reservoir properties such as permeability and porosity. The traditional approach typically involves several runs of the forward simulator. This may not scale well when the optimization is to be performed at the full field-level and over multiple geologic realizations. A machine-learning (ML) based approach which is quick and scalable for rate optimization over multiple geologic realizations is proposed instead. The training data for the model is generated by running the forward simulator with randomly assigned well rates using multiple geologic realizations. A reduced order representation of the permeability heterogeneity in each of the realizations is derived using a grid connectivity transformation (GCT). This step involves finding basis functions corresponding to the different modal frequencies of the grid connectivity represented by the grid Laplacian. The projection of the heterogeneous property field along these basis functions gives the basis coefficients that form the reduced order representation. Subsequently, for each training datapoint, streamlines are traced and the minimum time of flight (TOF) representing the tracer breakthrough time at each producer is recorded. The basis coefficients and well rates are fed to a machine learning model as input and the minimum TOF at the producers forms the output of the model. This trained model can then be used along with an optimizer for computing the optimal injection rates to maximize the injection sweep efficiency. This corresponds to minimizing the variance in the minimum TOF within each well group. Different architectures of neural network are tested using 5-fold cross validation to decide the best ML model to compute the streamline time of flight. The trained model is used to perform well rate optimization over multiple realizations of geology by using a risk tolerance penalty. The optimal well rates thus obtained are compared with two cases: a) equal well rates assigned to all injectors and producers and b) well rates obtained by optimizing over a single realization without considering the uncertainty in geology. The optimal well rates are seen to offer better oil recovery and sweep efficiency than both cases.

02 PETROLEUM↗
DOE OSTI2023

CLC: A cross-level program characterization method

Characterization of program execution plays a key role in performance improvement. There are numerous transformations applied to each step that a program takes on its lowering from source code to a compiler intermediate representation to machine language to microarchitecture-specific execution. The unpredictable benefit of each transformation step could lead a notionally superior algorithm to exhibit inferior performance once actually run, and it can be hard to discern which step in the transformation path contradicted the code developer’s assumptions. Conventional approaches to program-execution characterization consider the behavior after only a single one of those steps, which limits the information that can be provided to the user. To help address the issue of myopic views of program execution, this paper presents a novel cross-level characterization approach for understanding the behavior of program execution at different levels in the process of writing, compiling, and running a program. We show that this approach provides a richer view of the sources of performance gains and losses and helps identify program execution in a more accurate manner.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Turbulence-parameter estimation for current-energy converters using surrogate model optimization

Surrogate models maximize information utility by building predictive models in place of computational or experimentally expensive model runs. Marine hydrokinetic current energy converters require large-domain simulations to estimate array efficiencies and environmental impacts. Meso-scale models typically represent turbines as actuator discs that act as momentum sinks and sources of turbulence and its dissipation. An OpenFOAM model was developed where actuator disc k-ε turbulence was characterized using an approach developed for flows through vegetative canopies. Turbine-wake data from laboratory flume experiments collected at two influent turbulence intensities were used to calibrate parameters in the turbulence-source terms in the k-ε equations. Additionally, parameter influences on longitudinal wake profiles were estimated using Gaussian process regression with subsequent optimization minimizing the objective function within 3.1% of those obtained using the full model representation, but for 74% of the computational cost (far fewer model runs). This framework facilitates more efficient parameterization of the turbulence-source equations using turbine-wake data.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Multi-decade high-resolution regional hindcasts for wave energy resource characterization in U.S. coastal waters

Long-term, high-resolution, regional wave hindcast datasets were generated using unstructured-grid Simulating WAves Nearshore (SWAN) models for the U.S. coastal waters to support nearshore wave energy development in the U.S. including those bordering U.S. territorial islands. The model domains resolved the entire U.S. exclusive economic zones, with a spatial resolution of approximately 200 m nearshore. The regional SWAN models were driven by the global WAVEWATCH III® model outputs and run for a 42-year period from 1979 to 2020. Extensive model validations were performed using buoy observations and altimeter data. Regional resource characterization was performed based on hindcast data points at 2 km from shore and along the 100 m isobath. Aggregations of wave resource parameters were produced, and spatial and seasonal variations were analyzed for all the regions. Wave resource metrics recommended by international standards, including a 3-hour time series of six resource parameters, hourly frequency- and directionally resolved wave spectra at selected “virtual buoy” locations, and average-annual values of omni-directional wave power, significant wave height, and energy period are publicly disseminated through an Amazon Web Service and a Marine Energy Atlas web application tool to facilitate wave energy research and a wide range of coastal ocean applications.

16 TIDAL AND WAVE POWER↗