Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “US DOE”

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

Characterization of Spent Purolite A530E Resin with Implications for Long-Term Radioactive Contaminant Removal

Direct removal of pertechnetate anion (TcO4-) from the subsurface contaminated plumes at legacy nuclear weapons production sites, is effective remediation strategy to prevent their spread to adjacent aquifers but is also challenging because TcO4- is a trace component of the contaminated groundwater typically containing large access of other anionic constituents, such nitrate, sulfate, chloride and others. It can be achieved through ion exchange treatment using Purolite A530E resin which to date however has only been evaluated under laboratory conditions and many questions regarding its long-term operational performance remain. To address this need, here we report comprehensive characterization of the spent Purolite A530E resin which processed over 5.38x109 L of contaminated groundwater and successfully removed about 3.78 Ci of Tc-99 during four years of operation at the 200 West Pump & Treat facility at the US DOE Hanford site. This Tc-99 loading constitutes however only about 0.034% of the theoretical capacity of the resin which retained significant amounts of sulfate. Among other radioactive contaminants, small quantities of U-238, Co-60, and I-129 were retained by the resin. The total loading of iodine (combined 1-127 and I-129) on the spent resin exceeded that of Tc-99. To elucidate the mechanism of iodine retention, ion exchange behavior of iodide and iodate was investigated. Purolite A530E resin exhibited highly efficient uptake of iodide and only moderate affinity for iodate in accord with their Gibbs energy of hydration. Sorption isotherms for both anions obeyed Freundlich model.

Levitskaia, Tatiana G.↗
DOE OSTI2024

Development of multi-scale computational frameworks to solve fusion materials science challenges

Over the past two decades, the US-DOE has funded multiple projects that rely on high-performance computing and exascale computing platforms to accelerate scientific discoveries and address grand scientific challenges, such as harnessing fusion energy. In this article, we review in detail one of these efforts aimed at enhancing our capability to model plasma-facing materials subject to plasma and high-energy ion/neutron irradiation. The plasma surface interactions project has built a multi-scale modeling framework where many of the plasma- and high-energy ion/neutron irradiation-induced effects occurring in tungsten are explored. Here, this knowledge is used to develop atomistically-informed, high-fidelity continuum and meso-scale models that can be validated against experiments. We review the developments within this project, with attention to experimental validation efforts, and specifically highlight activities associated with: helium bubble bursting and equation of state, and hydrogen-helium interactions in tungsten; atomistically-informed model development for beryllium-tungsten material mixing; coupling of scrape-of-layer plasma, sheath and material models; and coupling of stochastic cluster-dynamics and crystal plasticity models to address radiation effects in tungsten under stress. Finally, we present how the project is preparing for future computational architectures, for instance through efforts to adapt atomistic methods to exascale computing.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Bipolar plate flow channel designs for vanadium redox flow battery: a review

Vanadium redox flow battery is one of the preferred systems for grid scale energy storage due to long service life (>20000 cycles), higher efficiency (>85 %), deep discharge capability (>95 % DoD), and inherent scalability and safety. Among system components, flow field is most critical as it governs electrolyte distribution, mass transport and hydraulic performance. Here, this review examines emerging flow channel geometries, highlighting the impact of channel width (0.66 to 1.5 mm), depth (1 to 1.5 mm) and land width (0.5 to 1.5 mm) can reduce pressure drop to <10 kPa while enabling power densities above 600 mW.cm -2 . A key finding is that low channel width to depth ratio (<1) enhances voltage and energy efficiencies by improving under rib convection. The review also provides an overview of progress and perspective in bipolar plate materials, manufacturability, and shunt current mitigation strategies for stack scaling up. Cost analysis emphasizes the influence of flow field designs on the levelized cost of storage and pathways towards the US DOE's $\$$0.05 per kWh energy generation cost target. In addition, opportunities for AI/ML/DT tools assisted design and data driven optimization strategies are also outlined to accelerate next generation flow field development.

Efficiency optimization↗
DOE OSTI2023

Advanced silicon tracking detector developments for the future Electron-Ion Collider

The proposed Electron-Ion Collider (EIC) will operate high-luminosity high-energy electron+proton and electron+nucleus collisions at the collision energies from 20 GeV to 141 GeV to solve several fundamental questions in the high energy and nuclear physics fields. Its instantaneous luminosity can reach 10 33-34 cm -2 s -1 and the bunching crossing rate is around 10 ns. The EIC project has received CD1 approval from the US DOE in 2021 and moves toward the machine design and preparation for construction. To realize various particle measurements with high precision at the future EIC, a low material-budget and high-granularity silicon vertex and tracking detector with fine spatial and momentum resolutions and nearly 4π solid angle coverage is desired. The Monolithic Active Pixel Sensor (MAPS) and AC Coupled Low Gain Avalanche Diode (AC-LGAD) technologies stand out of several advanced technology options for the EIC silicon vertex and tracking detector subsystems. The MAPS technology has advanced features of low material budget, low power consumption, good radiation resistance and fine spatial resolution. The AC-LGAD technology can achieve fast timing resolution. Latest studies and progress of the EIC silicon vertex and tracking detector conceptual design, performance validations in simulation and ongoing MAPS and AC-LGAD R&D will be shown. Furthermore, schedule and plan of the EIC project detector development will be discussed as well.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Measurement and analysis of GPU-accelerated applications with HPCToolkit

To address the challenge of performance analysis on the US DOE’s forthcoming exascale supercomputers, Rice University has been extending its HPCToolkit performance tools to support measurement and analysis of GPU-accelerated applications. To help developers understand the performance of accelerated applications as a whole, HPCToolkit’s measurement and analysis tools attribute metrics to calling contexts that span both CPUs and GPUs. To measure GPU-accelerated applications efficiently, HPCToolkit employs a novel wait-free data structure to coordinate monitoring and attribution of GPU performance. To help developers understand the performance of complex GPU code generated from high-level programming models, HPCToolkit constructs sophisticated approximations of call path profiles for GPU computations. To support fine-grained analysis and tuning, HPCToolkit uses PC sampling and instrumentation to measure and attribute GPU performance metrics to source lines, loops, and inlined code. To supplement fine-grained measurements, HPCToolkit can measure GPU kernel executions using hardware performance counters. To provide a view of how an execution evolves over time, HPCToolkit can collect, analyze, and visualize call path traces within and across nodes. Finally, on NVIDIA GPUs, HPCToolkit can derive and attribute a collection of useful performance metrics based on measurements using GPU PC samples. Here, we illustrate HPCToolkit’s new capabilities for analyzing GPU-accelerated applications with several codes developed as part of the Exascale Computing Project.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Nitrogen-Doped PtNi Catalysts on Polybenzimidazole-Functionalized Carbon Support for the Oxygen Reduction Reaction in Polymer Electrolyte Membrane Fuel Cells

PtM (M = 3d transition metals) alloys are known as the promising oxygen reduction reaction catalysts and have been considered as the replacement of pure Pt catalysts for the commercialization of proton exchange membrane fuel cells. Although great progress has been made in the past three decades, the performance and durability of PtM catalysts still face stringent challenges from practical applications. Functionalization of a catalyst carbon support with nitrogen-contained groups can add charges onto its surface, which can be utilized to build a more complete ionomer/catalyst interface, to reduce the catalyst particle size, and to improve particle size distribution. Nitriding of PtNi catalysts can effectively improve the catalyst activity and stability by the modification of lattice strain. Hereby, we propose a synergistic approach of combining polybenzimidazole-grafted Vulcan XC72 carbon as the catalyst carbon support and the nitriding of PtNi to develop PtNiN/XC72-polybenzimidazole catalysts. Such PtNiN/XC72-PBI catalysts exhibit the excellent performance of fuel cell membrane electrode assembly (i.e., mass activity, 440 mA mgPt –1 ; electrochemical surface area, 51 m 2 gPt –1 ; and rated power density, 836 mW cm –2 ) as well as promising catalyst stability. The developed PtNiN/XC72-PBI meets the US DOE 2020 targets of mass activity for the fuel cell catalysts. This work provides a novel approach and a promising pathway on the development of the catalyst using such a synergistic approach-modification of the catalyst structure by nitrogen doping and functionalization of carbon support by polybenzimidazole for both high performance and high durability.

25 ENERGY STORAGE↗
DOE OSTI2022

A solar energy desalination analysis tool, sedat, with data and models for selecting technologies and regions

Abstract There is interest for desalination technologies powered by solar energy as arid areas are typically bestowed with good solar potential. In response to a US DOE call for solar desalination analysis tools, we developed an open-source solar energy desalination analysis tool, sedat , for techno-economical evaluation of desalination technologies and selection of regions with the highest potential for using solar energy to power desalination plants. It is expected that this software will simplify the planning, design, and valuation of solar desalination systems in the U.S. and worldwide. Sedat uses Dash for integrating various layers of large volumes of GIS data with Python-based models of solar energy generation and desalination technologies. It derives time-series of energy generation and water production, with details of plant performance and suggestions for improving the solar-desalination coupling. This paper summarizes the various phases of the tool’s development, presents example results showing the potential, under multiple objectives, of solar desalination in parts of the U.S. southwest, and discusses method details that would be useful for future model development.

14 SOLAR ENERGY↗
DOE OSTI2020

Risk-Informed Safety Analysis for Accident Tolerant Fuels

Accident Tolerant Fuels (ATF) are being tested by different nuclear vendors and research organizations and their introduction in the US light water reactors fleet is planned for the second half of the 2020`s. In the framework of the US-DOE Light-Water Reactor Sustainability program, Risk-Informed Systems Analysis pathway, researches are being conducted at the Idaho National Laboratory for developing tools and methods that can help the industry quantifying the benefits from the adoption of ATF technology. In this paper we describe the developed risk-informed methodology including the safety analysis codes improvements, and we present some results for selected accident scenarios. The developed methodology combines INL state-of-the-art deterministic Best Estimate code RELAP5-3D, and Probabilistic Risk Analysis tools RAVEN and SAPHIRE. The analyses are performed on a three-loops pressurized water reactor, simulating Station Blackout and LB-LOCA accidents and considering near-term ATFs, or FeCrAl and Chromium-coated clads. Finally, we show how, applying our methodology, the new Core Damage Frequency can be assessed. The results indicated that the main benefits in introducing near-term ATFs concern a significant reduction in the hydrogen production during accident conditions. No significant Core Damage Frequency reduction was found.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Mesoscale to Microscale Coupling for Wind Energy Applications: Addressing the Challenges

The purpose of the US DOE's Mesoscale to Microscale Coupling (MMC) Project is to develop, verify, and validate physical models and modeling techniques that bridge the most important atmospheric scales that determine wind plant performance and reliability. The project seeks to create a new predictive numerical simulation capability that represents a range of dynamic atmospheric flow conditions impacting wind plant performance.

17 WIND ENERGY↗
DOE OSTI2023

High average power ultrafast laser technologies for driving future advanced accelerators

Large scale laser facilities are needed to advance the energy frontier in high energy physics and accelerator physics. Laser plasma accelerators are core to advanced accelerator concepts aimed at reaching TeV electron electron colliders. In these facilities, intense laser pulses drive plasmas and are used to accelerate electrons to high energies in remarkably short distances. A laser plasma accelerator could in principle reach high energies with an accelerating length that is 1000 times shorter than in conventional RF based accelerators. Notionally, laser driven particle beam energies could scale beyond state of the art conventional accelerators. LPAs have produced multi GeV electron beams in about 20 cm with relative energy spread of about 2 percent, supported by highly developed laser technology. This validates key elements of the US DOE strategy for such accelerators to enable future colliders but extending best results to date to a TeV collider will require lasers with higher average power. While the per pulse energies envisioned for laser driven colliders are achievable with current lasers, low laser repetition rates limit potential collider luminosity. Applications will require rates of kHz to tens of kHz at Joules of energy and high efficiency, and a collider would require about 100 such stages, a leap from current Hz class LPAs. This represents a challenging 1000 fold increase in laser repetition rates beyond current state of the art. Here, this whitepaper describes current research and outlook for candidate laser systems as well as the accompanying broadband and high damage threshold optics needed for driving future advanced accelerators.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Cryogenic and safety design of the future high field cable test facility at Fermilab

The HFVMTF (High Field Vertical Magnet Test Facility) is a new experimental facility under development at Fermi National Accelerator Laboratory (FNAL) to test large superconducting magnets (up to 20 tons weight and 1.3 m diameter) in a double bath superfluid helium cryostat (1.9 K and 1.2 bar). Coupled with a superconducting dipole magnet fabricated by Lawrence Berkeley National Laboratory (LBNL), this facility will be able to test future high-temperature superconductor (HTS) cables under a background magnetic field of 15 T for fusion magnets. This paper describes the design of the cryostat and its 1.4-meter diameter lambda plate, as well as the different components for a safe operation of the facility, even during critical events such a magnet quench or a vacuum breaking situation. The project is funded by US DOE Offices of Science, High Energy Physics (HEP), and Fusion Energy Sciences (FES).

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

A Perspective on Scalable AI on High-Performance Computing and Leadership Class Supercomputing Facilities [Industrial and Governmental Activities]

Many scientific applications that support the mission of the US Department of Energy (US-DoE) require modeling complex engineering and/or physical systems. Here, examples of such complex systems arise from: (a) materials science to develop new compounds with exceptional mechanical and thermodynamical properties (e.g., resistance to mechanical stresses and high temperatures), (b) structural and nuclear engineering to model the temporal evolution of the structural damage of concrete shields exposed to continuous neutron and gamma radiations emitted by the nuclear reactor core, (c) urban sciences (e.g., transportation and smart buildings), and (d) power grid systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Investigation of Membrane Chemical Degradation as a Function of Catalyst Platinum Loading

Membrane chemical degradation is one of many factors that can impact fuel cell durability. Additionally, the fuel cell’s lifetime heavily depends on the membrane and its ability to maintain chemical and mechanical integrity. Previous studies indicate that chemical degradation is due to the formation of hydroxyl radicals that attack the polymer structure resulting in membrane thinning, pinhole formation, and the release of fluoride and sulfate ions. Membrane durability was investigated using ultra-low Pt electrode loadings (≤ 0.1 mg Pt cm -2 ). Accelerated stress testing (US-DOE protocols) demonstrated that the degradation rate was found to increase with higher Pt loadings. This is most likely due to more heterogeneous sites for radical formation due to hydrogen crossover to the cathode. We also explored membrane degradation rates while varying catalyst layer thickness, ionomer to carbon ratio, and types of carbon support. All of the aforementioned variables impact the membrane degradation rates.

25 ENERGY STORAGE↗
DOE OSTI2023

HydraGNN 2.0

HydraGNN is an Oak Ridge National Laboratory (ORNL)-branded implementation of distributed multi-tasking graph neural networks that supports several scientific applications within the ORNL portfolio to support the US-DOE mission. New or improved capabilities included in v2.0.0 release are as follows: 1) Enhancement in message passing layers through class inheritance 2) Adding transformation to ensure translation and rotation invariance 3) Supporting various optimizers 4) Atomic descriptors 5) Integration with continuous CI test 6) Distributed printouts and timers 7) Profiling 8) Support of ADIOS2 for scalable data loading 9) Large-scale system support, including Summit (ORNL) and Perlmutter (NERSC)

Lupo Pasini, Massimiliano [Oak Ridge National Labo↗
DOE OSTI2026

JACC.jl: Julia for Accelerators

CPU/GPU portable array, parallel_for/parallel_reduce in Julia for productive science. Funded by the US DOE Advanced Scientific Computing Research (ASCR).

Fackler, Philip [Oak Ridge National Laboratory (OR↗
DOE OSTI2020

Kraken2 Metagenomic Virus Database

The Database: Kraken2 [1] database built from a classification tree containing over 700k metagenomic viruses from JGI IMG/VR [2]. (1) Wood, D. E., Lu, J., and Langmead, B. (2019). Improved metagenomic analysis with Kraken 2. Genome Biol., 20(1), 1–13. doi: 10.1186/s13059-019-1891-0 (2) Paez-Espino D, Chen I-MA, Palaniappan K, Ratner A, Chu K, Szeto E, et al. IMG/VR: a database of cultured and uncultured DNA Viruses and retroviruses. Nucleic Acids Res. 2017;45:D457-65. For Paper: Title: A k-mer based approach for virus classification in metatranscriptomic and metagenomic samples identifies viral associations in the Populus phytobiome and autism brains Abstract Background Viruses are an underrepresented taxa in the study and identification of microbiome constituents; however, they play an important role in health, microbiome regulation, and transfer of genetic material. Only a few thousand viruses have been isolated, sequenced, and assigned a taxonomy, which further limits the ability to identify and quantify viruses in the microbiome. Additionally, the vast diversity of viruses represents a challenge for classification, not only in constructing a viral taxonomy, but also in identifying similarities between a virus' genotype and its phenotype. However, the diversity of viral sequences can be leveraged to classify their sequences in metagenomic and metatranscriptomic samples. Methods To identify viruses in transcriptomic and genomic samples, we developed a dynamic programming algorithm for creating a classification tree out of 715,672 metagenome viruses. To create the classification tree, we clustered proportional similarity scores generated from the k-mer profiles of each of the metagenome viruses. We then integrated the viral classification tree with the NCBI taxonomy for use with ParaKraken, a metagenomic/transcriptomic classifier. Results To illustrate the breadth of our utility for classifying viruses with ParaKraken, we analyzed data from a plant metagenome study identifying the differences between two Populus genotypes in three different compartments and on a human metatranscriptome study identifying the differences between Autism Spectrum Disorder patients and controls in post mortem brain biopsies. In the Populus study, we identified genotype and compartment specific viral signatures, while in the Autism study we identified a significant increased abundance of eight viral sequences in Autism brain biopsies. Conclusion Viruses represent an important aspect of the microbiome. The ability to classify viruses represents the first step in being able to better understand their role in the microbiome. The viral classification method presented here allows for more complete identification of viral sequences for use in identifying associations between viruses and the host and viruses and other microbiome members. Acknowledgements and Funding This research used resources of the Oak Ridge Leadership Computing Facility at the Oak Ridge National Laboratory, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725. This research was also supported by the Plant-Microbe Interfaces Scientific Focus Area in the Genomic Science Program, the Office of Biological and Environmental Research (BER) in the U.S. Department of Energy Office of Science, and by the Department of Energy, Laboratory Directed Research and Development funding (ProjectID 8321), at the Oak Ridge National Laboratory. Oak Ridge National Laboratory is managed by UT-Battelle, LLC, for the US DOE under contract DE-AC05-00OR22725. This research used resources of the Compute and Data Environment for Science (CADES).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Micromechanical Properties of the SiC and Pyrolytic Carbon Layers in Tristructural-Isotropic Coated Particles

Tristructural isotropic (TRISO) coated particle fuel was initially developed for high temperature gas cooled reactors (HTGR) and has been proposed for several other advanced reactor concepts. The design of TRISO particles focuses on preventing the release of fission products in normal and off-normal reactor conditions. The particle design features an actinide bearing fuel kernel that is surrounded by three pyrolytic carbon (PyC) layers and a silicon carbide layer (SiC). The mechanical stability of the particle and fission product retention for both metallic and gaseous fission products depend on the SiC layer. Post irradiation examination (PIE) of TRISO fuel from the first two US DOE Advanced Gas Reactor Fuel Development and Qualification Program irradiation campaigns, AGR-1 and AGR-2, had identified a low rate of particles exhibiting cracking in the SiC layer that did not propagate across the SiC layer. While cracking in the SiC is rare for test conditions and particles associated with the AGR program, understanding the stress state and mechanical properties of the SiC and PyC layers related to particle architecture can aid predicting thermomechanical response of TRISO fuel under the prescribed operation envelope and beyond as well as aiding in the development of similar fuel concepts for other advanced reactors. The presented investigation shows the relationship of the mechanical properties and mechanical response (e.g., understanding crack propagation) of the SiC and PyC layers relative to position within the particle. Testing was conducted on the inner and outer PyC layers of TRISO particles to quantify differences in mechanical behavior.

Montoya, Katherine [ORNL] (ORCID:0000000326955086)↗
DOE OSTI2020

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including 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 ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element 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 remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗