Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “example”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Vector Field Decompositions Using Multiscale Poisson Kernel

Extraction of multiscale features using scale-space is one of the fundamental approaches to analyze scalar fields. However, similar techniques for vector fields are much less common, even though it is well known that, for example, turbulent flows contain cascades of nested vortices at different scales. The challenge is that the ideas related to scale-space are based upon iteratively smoothing the data to extract features at progressively larger scale, making it difficult to extract overlapping features. Instead, we consider spatial regions of influence in vector fields as scale, and introduce a new approach for the multiscale analysis of vector fields. Rather than smoothing the flow, we use the natural Helmholtz-Hodge decomposition to split it into small-scale and large-scale components using progressively larger neighborhoods. Our approach creates a natural separation of features by extracting local flow behavior, for example, a small vortex, from large-scale effects, for example, a background flow. We demonstrate our technique on large-scale, turbulent flows, and show multiscale features that cannot be extracted using state-of-the-art techniques.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Sampling in Long-Screened Wells: Issues, Misconceptions, and Solutions

The issues associated with long-screened wells (LSWs) (and open boreholes) at contaminated sites are well documented in the groundwater literature but are still not fully appreciated in practice. As established in seminal and review papers going back over three decades, the interpretation of sampling results from LSWs is challenging in the presence of vertical hydraulic gradients and borehole flow; furthermore, LSWs allow for vertical redistribution of contamination between aquifer layers. Acknowledgment of these issues has led to the development of new technologies and well designs to enable discrete-zone monitoring (DZM), yet LSWs remain common for many reasons, for example, as multipurpose wells, for geophysical logging, and (or) as legacy installations. Despite the literature on LSWs and despite the adoption of DZM at many sites, the use of LSWs persists and the challenges of interpreting sampling results from LSWs remain. In this issue paper, we provide a conceptual overview of the problems posed by LSWs and review existing literature and past work to improve the interpretation of sampling in LSWs. We draw on experience from previous studies at the Hanford Site in eastern WA, USA, and use synthetic examples to illustrate key concepts and challenges for interpretation. A recently published analytical modeling framework is used to develop illustrative synthetic examples and demonstrate a workflow for building scientific intuition to understand issues around interpreting samples from LSWs, which is critical to effective characterization and groundwater remediation at sites with LSWs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A seed-like proteome in oil-rich tubers

There are numerous examples of plant organs or developmental stages that are desiccation-tolerant and can withstand extended periods of severe water loss. One prime example are seeds and pollen of many spermatophytes. However, in some plants, also vegetative organs can be desiccation-tolerant. One example are the tubers of yellow nutsedge ( Cyperus esculentus ), which also store large amounts of lipids similar to seeds. Interestingly, the closest known relative, purple nutsedge ( Cyperus rotundus ), generates tubers that do not accumulate oil and are not desiccation-tolerant. We generated nanoLC-MS/MS-based proteomes of yellow nutsedge in five replicates of four stages of tuber development and compared them to the proteomes of roots and leaves, yielding 2257 distinct protein groups. Our data reveal a striking upregulation of hallmark proteins of seeds in the tubers. A deeper comparison to the tuber proteome of the close relative purple nutsedge ( C. rotundus ) and a previously published proteome of Arabidopsis seeds and seedlings indicates that indeed a seed-like proteome was found in yellow but not purple nutsedge. This was further supported by an analysis of the proteome of a lipid droplet-enriched fraction of yellow nutsedge, which also displayed seed-like characteristics. One reason for the differences between the two nutsedge species might be the expression of certain transcription factors homologous to ABSCISIC ACID INSENSITIVE3, WRINKLED1, and LEAFY COTYLEDON1 that drive gene expression in Arabidopsis seed embryos.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Atomic-Level, Energy-Conversion Heat Transfer

Abstract Heat is stored in quanta of kinetic and potential energies in matter. The temperature represents the equilibrium and excited occupation (boson) of these energy conditions. Temporal and spatial temperature variations and heat transfer are associated with the kinetics of these equilibrium excitations. During energy-conversion (between electron and phonon systems), the occupancies deviate from equilibria, while holding atomic-scale, inelastic spectral energy transfer kinetics. Heat transfer physics reaches nonequilibrium energy excitations and kinetics among the principal carriers, phonon, electron (and holes and ions), fluid particle, and photon. This allows atomic-level tailoring of energetic materials and energy-conversion processes and their efficiencies. For example, modern thermal-electric harvesters have transformed broad-spectrum, high-entropy heat into a narrow spectrum of low-entropy emissions to efficiently generate thermal electricity. Phonoelectricity, in contrast, intervenes before a low-entropy population of nonequilibrium optical phonons becomes a high-entropy heat. In particular, the suggested phonovoltaic cell generates phonoelectricity by employing the nonequilibrium, low-entropy, and elevated temperature optical-phonon produced population—for example, by relaxing electrons, excited by an electric field. A phonovoltaic material has an ultranarrow electronic bandgap, such that the hot optical-phonon population can relax by producing electron-hole pairs (and power) instead of multiple acoustic phonons (and entropy). Examples of these quanta and spectral heat transfer are reviewed, contemplating a prospect for education and research in this field.

Engineering↗
DOE OSTI2022

Surrogate Modeling of Nonlinear Dynamic Systems: A Comparative Study

Surrogate models play a vital role in overcoming the computational challenge in designing and analyzing nonlinear dynamic systems, especially in the presence of uncertainty. This paper presents a comparative study of different surrogate modeling techniques for nonlinear dynamic systems. Four surrogate modeling methods, namely, Gaussian process (GP) regression, a long short-term memory (LSTM) network, a convolutional neural network (CNN) with LSTM (CNN-LSTM), and a CNN with bidirectional LSTM (CNN-BLSTM), are studied and compared. All these model types can predict the future behavior of dynamic systems over long periods based on training data from relatively short periods. The multi-dimensional inputs of surrogate models are organized in a nonlinear autoregressive exogenous model (NARX) scheme to enable recursive prediction over long periods, where current predictions replace inputs from the previous time window. Three numerical examples, including one mathematical example and two nonlinear engineering analysis models, are used to compare the performance of the four surrogate modeling techniques. The results show that the GP-NARX surrogate model tends to have more stable performance than the other three deep learning (DL)-based methods for the three particular examples studied. The tuning effort of GP-NARX is also much lower than its deep learning-based counterparts.

42 ENGINEERING↗
DOE OSTI2022

Critical Shear Stress for Erosion Under Laminar Jet Flow

Particle erosion is an important process in both natural and manmade environments. For example, in natural environments particle or sediment erosion rates determine the course of rivers, depth of streams, and fate of deltas. The interplay between erosion rates and flow patterns can have dramatic impact on shorelines. Where hydrodynamic jet flows or rip currents change direction, sediments may accumulate or disperse altering the shape of coastlines and the ecosystems and economies that rely on them. In man-made systems, jet flows may be used to mobilize particles from surfaces or within industrial scale mixing tanks. For example, in the nuclear industry, impinging jet mixers may be used to mobilize and suspend sludges at the bottom of waste tanks. In each of these examples, the Shields diagram is useful to determine the conditions under which erosion occurs. This diagram plots the scaled critical shear stress for erosion as a function of the particle size as embedded in particle Reynolds numbers, Archimedes numbers, or their equivalents. Although the data, that carries a large spread, is largely in agreement with mathematical models, for the smallest particles, model predictions diverge from each other and from experimental data. Indeed, there is substantial uncertainty about the Shields diagrams at relatively small particle sizes and flow rates. Here we evaluate critical shear stresses for erosion under laminar jet flow. To date critical shear stresses for erosion, as embodied in Shields diagrams, have largely been developed for turbulent flow both experimental and mathematically. However, curves on the Shields diagrams diverge among various models at lower particle Reynolds numbers. Indeed, how critical shear stresses for erosion develop under fully laminar jet flow conditions (not simply laminar boundary layers under turbulent flow) remains unclear. Here we address this gap, which is particularly important for mobilization of the smallest particles. We find the Shields parameter to be inversely proportional to the particle Reynolds numbers under laminar jet flow conditions similar to turbulent conditions.

critical shear stresses, erosion rates, effective ↗
DOE OSTI2020

Role of consistent terminology in XPS reproducibility

The use of inconsistent and incorrect terminology in scientific publications contributes to misunderstanding, confusion, and erroneous results in the scientific literature. This issue is of particular importance in x-ray photoelectron spectroscopy (XPS) that is in widespread use for many different purposes by scientists with diverse backgrounds. A set of surface analysis terminology, approved through consensus by international experts, has been developed by International Organization for Standardization Technical Committee TC201 on Surface Chemical Analysis. To encourage wide use, the terminology is accessible at several web sites free of charge. This short overview provides examples to highlight the importance of agreed terminology in eliminating confusion between similar terms. Examples are provided of terms that are commonly misused or confused in the literature. Other examples highlight terminology that provides a common basis for comparing instrument parameters and performance. As science advances, it is important to clarify terminology for describing evolving concepts and developments important to XPS.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2022

Grassmannian Diffusion Maps--Based Dimension Reduction and Classification for High-Dimensional Data

This work introduces the Grassmannian diffusion maps (GDMaps), a novel nonlinear dimensionality reduction technique that defines the affinity between points through their representation as low-dimensional subspaces corresponding to points on the Grassmann manifold. Here, the method is designed for applications, such as image recognition and data-based classification of constrained high-dimensional data where each data point itself is a high-dimensional object (i.e., a large matrix) that can be compactly represented in a lower-dimensional subspace. The GDMaps is composed of two stages. The first is a pointwise linear dimensionality reduction wherein each high-dimensional object is mapped onto the Grassmann manifold representing the low-dimensional subspace on which it resides. The second stage is a multipoint nonlinear kernel-based dimension reduction using diffusion maps to identify the subspace structure of the points on the Grassmann manifold. To this end, an appropriate Grassmannian kernel is used to construct the transition matrix of a random walk on a graph connecting points on the Grassmann manifold. Spectral analysis of the transition matrix yields low-dimensional Grassmannian diffusion coordinates embedding the data into a low-dimensional reproducing kernel Hilbert space. Further, a novel data classification/recognition technique is developed based on the construction of an overcomplete dictionary of reduced dimension whose atoms are given by the Grassmannian diffusion coordinates. Three examples are considered. First, a "toy" example shows that the GDMaps can identify an appropriate parametrization of structured points on the unit sphere. The second example demonstrates the ability of the GDMaps to revealing the intrinsic subspace structure of high-dimensional random field data. In the last ex- ample, a face recognition problem is solved considering face images subject to varying illumination conditions, changes in face expressions, and occurrence of occlusions. The technique presented high recognition rates (i.e., 95% in the best case) using a fraction of the data required by conventional methods.

42 ENGINEERING↗
DOE OSTI2023

A Probabilistic Scheme for Semilinear Nonlocal Diffusion Equations with Volume Constraints

This work presents a probabilistic scheme for solving semilinear nonlocal diffusion equations with volume constraints and integrable kernels. The nonlocal model of interest is defined by a time-dependent semilinear partial integro-differential equation (PIDE), in which the integro-differential operator consists of both local convection-diffusion and nonlocal diffusion operators. Here, our numerical scheme is based on the direct approximation of the nonlinear Feynman–Kac formula that establishes a link between nonlinear PIDEs and stochastic differential equations. The exploitation of the Feynman–Kac representation avoids solving dense linear systems arising from nonlocal operators. Compared with existing stochastic approaches, our method can achieve first-order convergence after balancing the temporal and spatial discretization errors, which is a significant improvement of existing probabilistic/stochastic methods for nonlocal diffusion problems. Error analysis of our numerical scheme is established. The effectiveness of our approach is shown in two numerical examples. The first example considers a three-dimensional nonlocal diffusion equation to numerically verify the error analysis results. The second example presents a physics problem motivated by the study of heat transport in magnetically confined fusion plasmas.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

BuildingSync® v.2.7.0 (released 9.11.2025) [SWR-18-28]

BuildingSync® is a building data exchange schema to better enable integration between software tools and building data workflows. The schema's original use case was focused on commercial building energy audits; however, several additional use cases have been realized including building energy modeling and more high-level generic building data exchange. Version 2.7.0 adds new elements for file attachment feature and FederalBuilding, and generalizes usage of Optional Elements (e.g. EquipmentCondition, EquipmentID) to all assets/systems. BuildingSync helps streamline the data exchange process, improving the value of the data, minimizing duplication of effort for subsequent building data collection efforts (including audits), and facilitating the achievement of greater energy efficiency. This in done in part by standardizing on (a) reporting audits in an electronic format, (b) tracking proposed, implemented, and discarded energy conservation measures, and (c) storing building characteristics (at multiple levels) for audits, benchmarking, and building energy analysis. BuildingSync has several documents and tools available to help users understand how to best leverage BuildingSync. The list below are only a subset of the resources available. If new resources are discovered, then feel free to create a new pull request with the additions. Generic BuildingSync information is available on the DOE website and the project website. BuildingSync Examples - These examples are kept up to date and show a wide range of implementations. Any new update to BuildingSync is required to pass validation on these example files. BuildingSync Use Case Validator allows for users to determine if their instance complies with a specific use case for BuildingSync by checking if the required elements are implemented in an uploaded instance. An API is also provided for automated integration into other tools. Also, the website contains an easy way to view the entirety of the schema and how elements relate to the Building Exchange Data Exchange Specification. The Validator is open sourced here Use Case TestSuite provides a Python package for easier generation of BuildingSync use cases. BuildingSync use cases depend on the generation of schematron documents, which is time-consuming and difficult to implement well. The TestSuite allows users to define a use case using a more palatable CSV template, which it then turns into a Schematron document. The source code is available here. BuildingSync to OpenStudio/EnergyPlus. The translator is open sourced here. This project will translate a Level 1 (and partial Level 2) ASHRAE Energy Audit to a fully defined OpenStudio and EnergyPlus model. This project is in early Beta testing and any feedback is welcome!

Long, Nicholas [National Renewable Energy Lab. (NR↗
DOE OSTI2026

FAD-Toolset (Floating Array Design Toolset) [SWR-26-056]

The Floating Array Design (FAD) Toolset is a collection of tools for modeling and designing arrays of floating offshore structures. It was originally designed for floating wind systems but has applicability for many offshore applications. A core part of the FAD Toolset is the floating array model, which serves as a high-level library for efficiently modeling a floating array, such as a floating wind array. It combines site condition information and a description of the floating array design, and contains functions for evaluating the array's behavior considering the site conditions. For example, it combines information about site soil conditions, mooring line loads, and an array's anchor characteristics to estimate the holding capacity of each anchor. The library works in conjunction with the tools RAFT, MoorPy, and FLORIS to model floating platforms, wind turbines, mooring systems, power cables, and array wakes respectively. Layered on top of the floating array model is a set of design tools that can be used for algorithmically adjusting or optimizing parts of the a floating array. Specific tools existing for mooring lines, shared mooring systems, dynamic power cables, static power cable routing, and overall array layout. These capabilities work with the design representation and evaluation functions in the floating array model, and they can be applied by users in various combinations to suit different purposes. In addition to standalone uses of the FAD Toolset, a coupling has been made with Ard, (https://github.com/NLRWindSystems/Ard) a sophisticated and flexible wind farm optimization tool. This coupling allows Ard to use certain mooring system capabilities from FAD to perform layout optimization of floating wind farms with Ard's more advanced layout optimization capabilities. The FAD Toolset works with the IEA Wind Task 49 Ontology (https://github.com/IEAWindTask49/Ontology), which provides a standardized format for describing floating wind farm sites and designs. See example use cases in our examples folder (https://github.com/NLRWindSystems/FAD-Toolset/blob/main/examples/README.md) For working with the library, it is important to understand the floating array model structure, which is described more here: https://github.com/NLRWindSystems/FAD-Toolset/blob/main/fad/README.md.

Sirkis, Leah [National Laboratory of the Rockies (↗
DOE OSTI2020

On the vulnerability of data-driven structural health monitoring models to adversarial attack

Many approaches at the forefront of structural health monitoring rely on cutting-edge techniques from the field of machine learning. Recently, much interest has been directed towards the study of so-called adversarial examples; deliberate input perturbations that deceive machine learning models while remaining semantically identical. This article demonstrates that data-driven approaches to structural health monitoring are vulnerable to attacks of this kind. In the perfect information or ‘white-box’ scenario, a transformation is found that maps every example in the Los Alamos National Laboratory three-storey structure dataset to an adversarial example. Also presented is an adversarial threat model specific to structural health monitoring. The threat model is proposed with a view to motivate discussion into ways in which structural health monitoring approaches might be made more robust to the threat of adversarial attack.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

A fast particle-based approach for calibrating a 3-D model of the Antarctic ice sheet

We consider the scientifically challenging and policy-relevant task of understanding the past and projecting the future dynamics of the Antarctic ice sheet. The Antarctic ice sheet has shown a highly nonlinear threshold response to past climate forcings. Triggering such a threshold response through anthropogenic greenhouse gas emissions would drive drastic and potentially fast sea level rise with important implications for coastal flood risks. Previous studies have combined information from ice sheet models and observations to calibrate model parameters. These studies have broken important new ground but have either adopted simple ice sheet models or have limited the number of parameters to allow for the use of more complex models. These limitations are largely due to the computational challenges posed by calibration as models become more computationally intensive or when the number of parameters increases. Here, we propose a method to alleviate this problem: a fast sequential Monte Carlo method that takes advantage of the massive parallelization afforded by modern high-performance computing systems. We use simulated examples to demonstrate how our sample-based approach provides accurate approximations to the posterior distributions of the calibrated parameters. The drastic reduction in computational times enables us to provide new insights into important scientific questions, for example, the impact of Pliocene era data and prior parameter information on sea level projections. These studies would be computationally prohibitive with other computational approaches for calibration such as Markov chain Monte Carlo or emulation-based methods. We also find considerable differences in the distributions of sea level projections when we account for a larger number of uncertain parameters. For example, based on the same ice sheet model and data set, the 99th percentile of the Antarctic ice sheet contribution to sea level rise in 2300 increases from 6.5 m to 13.1 m when we increase the number of calibrated parameters from three to 11. With previous calibration methods, it would be challenging to go beyond five parameters. Here, this work provides an important next step toward improving the uncertainty quantification of complex, computationally intensive and decision-relevant models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Active and tunable nanophotonic metamaterials

Abstract Metamaterials enable subwavelength tailoring of light–matter interactions, driving fundamental discoveries which fuel novel applications in areas ranging from compressed sensing to quantum engineering. Importantly, the metallic and dielectric resonators from which static metamaterials are comprised present an open architecture amenable to materials integration. Thus, incorporating responsive materials such as semiconductors, liquid crystals, phase-change materials, or quantum materials (e.g., superconductors, 2D materials, etc.) imbue metamaterials with dynamic properties, facilitating the development of active and tunable devices harboring enhanced or even entirely novel electromagnetic functionality. Ultimately, active control derives from the ability to craft the local electromagnetic fields; accomplished using a host of external stimuli to modify the electronic or optical properties of the responsive materials embedded into the active regions of the subwavelength resonators. We provide a broad overview of this frontier area of metamaterials research, introducing fundamental concepts and presenting control strategies that include electronic, optical, mechanical, thermal, and magnetic stimuli. The examples presented range from microwave to visible wavelengths, utilizing a wide range of materials to realize spatial light modulators, effective nonlinear media, on-demand optics, and polarimetric imaging as but a few examples. Often, active and tunable nanophotonic metamaterials yield an emergent electromagnetic response that is more than the sum of the parts, providing reconfigurable or real-time control of the amplitude, phase, wavevector, polarization, and frequency of light. The examples to date are impressive, setting the stage for future advances that are likely to impact holography, beyond 5G communications, imaging, and quantum sensing and transduction.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hawaii Wave Surge Energy Converter (HAWSEC) OSU O.H. Hinsdale Basin

The following information and metadata applies to both the Phase I (Hydrodynamics) and Phase II (Full System Power Take-Off) zip folders which contain testing data from the OSU (Oregon State University) O.H. Hinsdale Wave Research Laboratory, from both OSU and the University of Hawaii at Manoa (UH). See zip folders provided further below in the downloads section. For experimental data of the full system, including PTO, see Phase II dataset. There are two main directories in each Phases's zip folder: "OSU_data" and "UH_data". The "OSU_data" directory contains data collected from their DAQ (data acquisition system), which includes all wave gauge observations, as well as body motions derived from their Qualisys motion tracking system. The organization of the directory follows OSU's convention. Detailed information on the instrument setup can be found under "OSU_data/docs/setup/instm_locations". The experiments conducted are documented in the "OSU_data/docs/daq_logs", which provides the trial number to the corresponding data located under "OSU_data/data" in several formats (e.g., ".mat" and ".txt"). Inside the trial directory, data is provided for each of the instruments defined in "OSU_data/docs/setup/instm_locations". The "UH_data" directory contains data collected from their DAQ. The data is stored in a ".tdms" file format. There are free plug-ins for Microsoft Excel and MathWorks MATLAB to read the ".tdms" format. Below are a few links providing methods to read in the data, but a Google search should identify alternatives sources if these no longer exist (valid as of January 2024): Excel: http://www.ni.com/example/27944/en/ MATLAB: https://www.mathworks.com/matlabcentral/fileexchange/30023-tdms-reader The Excel plugin is recommend to get a quick overview of the data. The UH data is organized by directory name, in which the sub-directories for each experiment contains a directory whose name defines the wave height and period for the experimental data within. For example, a directory name "H02_T0275" corresponds to an experiment with wave height 0.1m and a period of 2.75s. For random wave data, the gamma value is also included in the directory name. For example, a directory name "H02_T0225_G18" corresponds to an experiment with a significant wave height of 0.2m, a peak period of 2.25s, and a gamma value of 1.8, with each spectra being a TMA spectrum. For the free decay experiments, the directory name is defined by the initial angular displacement. For example, a directory name "ang05_run01" corresponds to an experiment with an initial angular displacement of 5 degrees. There is a dataset in the UH data for each corresponding experiment defined in the OSU DAQ logs. The ".tdms" data is output from the DAQ at fixed intervals. Therefore, if multiple files are contained within the folder, the data will need to be stitched together. Within the UH dataset, there are two input channels from the OSU DAQ providing a random square wave signal for time synchronization ("ENV-WHT-0010") and a high/low signal ("ENV-WHT-0012") to identify when the wave maker is active (+5V). The UH data is logged as a collection of channel outputs. Channels not in use for the OSU testing (either Phase I or Phase II) are marked "nan" below. If the sensor is disconnected, it will record noise throughout the experiment. Below are the channel definitions in terms of what they measure: GPS Time = time CYL-POS-0001 = position between flap and fixed reference CYL-LCA-0001 = force between flap and hydraulic cylinder REC-LPT-0001 = nan REC-HPT-0001 = nan REC-HPT-0002 = nan REC-HPT-0003 = nan HHT-HPT-0001 = pressure at exhaust ("head" only) REC-FQC-0001 = nan REC-FQC-0002 = nan HHT-FQC-0001 = flow at exhaust ("head" only) ENV-WHT-0001 = nan ENV-WHT-0002 = nan ENV-WHT-0003 = nan ENV-WHT-0010 = random signal from OSU DAQ ENV-WHT-0012 = high/low signal from OSU DAQ Also included is a calibration curve to convert the string pot data to flap pi...

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Revised Kimberlina Groundwater Simulations

Kimberlina Groundwater Qaulity Simulations (Version 3: 2015 October) by Lawrence Livermore National Labratory. Submitted to the EDX team at the National Energy Technology Laboratory for research use. This dataset contains 562 updated groundwater quality simulations based on the Kimberlina site. It is a six-layer model including the aquifers and aquitards inbetween the ground surface and the reservoir. The input files and associated leakage information for these simulations are included in the zip file "nrap_kim_well31_fluxes and inputs.tar.gz". These 562 simulations are the successful results out of 1000 monte carlo simulations, with time length greater than 5 years. Each simulation consists 6 output files, therefore, there are 3372 files (6*562) in total. One single example simulation of full time length (200 year) is also provided in the zip file "single_simulation_200yr_sim0012.tar.gz". These 6 groundwater quality parameters are: - pH (use absolute value) - TDS (mg/kg) - Pressure (Pa) - CO2 Liquid Saturation (fraction range 0.0-1.0) - CO2 Gas Concentration (mol/L) - CO2 Liquid Concentration (mol/kg) Their correponding file names, take the first simulation for example, are: ph0001.ntab (for pH) , tds0001.ntab (for TDS), usnt0001.p.ntab (for pressure), usnt0001.s_liquid.ntab (for liquid saturation), usnt0001.c_co2_gas.ntab (for CO2 gas concentration) and usnt0001.c_co2_liquid.ntab (for CO2 liquid concentration). Each output file contains at least 19 columns (variables). The total number of columns depends on the simulation time length. The column data include the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the jth grid of y-axis 4. k: the kth grid of z-axis 5. element_ref: element type reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The rest (14, 15, 16...) variables are the simulated parameter values over time, take pH as an example, are: 14. ph000y: initial pH value. 15. ph001y: simulated pH value at the end of the 1st year. 16. ph002y: simulated pH value at the end of the 2nd year. ... (and so on)... The zip input files have the input fluxes of brine and CO2 (in kg/s) allocated to model layers 1,2,3 and 5 (layers 4 and 6 are assumed to be no-flow aquitards) through the following nodes: (src1 WELL 7 7 25 25 10 32) #Etchegoin (src2 WELL 7 7 25 25 33 72) #Macoma Chanac (src3 WELL 7 7 25 25 73 83) #Santa Margarita Mclure (src4 WELL 7 7 25 25 84 89) #Fruitvale Round Mountain (aquitard) (src5 WELL 7 7 25 25 90 96) #Olcese Sand (thief zone) (src6 WELL 7 7 25 25 97 99) #tremblor Freeman (aquitard) ex: brn1 is assigned to src 1 from i=(7,7), j=(25,25), k=(10,32) Edited by Ya-Mei Yang at National Energy Technology Laboratory

brine leakage↗
DOE OSTI2020

Field theories with a vector global symmetry

Motivated by recent discussions of fractons, we explore nonrelativistic field theories with a continuous global symmetry, whose charge is a spatial vector. We present several such symmetries and demonstrate them in concrete examples. They differ by the equations their Noether currents satisfy. Simple cases, other than the translation symmetry, are an ordinary (relativistic) one-form global symmetry and its nonrelativistic generalization. In the latter case the conserved charge is associated with a codimension-one spatial manifold, but it is not topological. More general examples involve charges that are integrated over the entire space. We also discuss the coupling of these systems to gauge fields for these symmetries. We relate our examples to known continuum and lattice constructions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Scale Invariance of the Homentropic Inviscid Euler Equations with Application to the Noh Problem

We investigate the inviscid compressible flow (Euler) equations constrained by an ”isentropic” equation of state (EOS), whose functional form in pressure is an arbitrary function of density alone. Under the aforementioned condition, we interrogate using symmetry methods the scale-invariance of the homentropic inviscid Euler equations. We find that under general conditions, we can reduce the inviscid Euler equations into a system of two coupled ordinary differential equations. To exemplify the utility of these results, we formulate two example scale-invariant, self-similar solutions. The first example includes a shock-free expanding bubble scenario, featuring a modified Tait EOS. The second example features the classical Noh problem, coupled to an arbitrary isentropic EOS. In this case, in order to satisfy the conditions set forth in the classical Noh problem, we find that the solution for the flow is given by a transcendental algebraic equation in the shocked density.

97 MATHEMATICS AND COMPUTING↗