Search NASA⌕ Search

SEARCH · Search NASA

Results for “Layer top”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Climatology of Cloud-Top Radiative Cooling in Marine Shallow Clouds

A one-year’s worth of near-global marine shallow single-layer cloud top radiative cooling (CTRC) is derived from a radiative transfer model with inputs from the satellite cloud retrievals and reanalysis sounding. The mean cloud top radiative flux divergence is 61 Wm -2 , decomposed into the longwave and shortwave components of 73 and -11 W m -2 , respectively. Equatorward of 30°N/S, the CTRC is largely a reflection of free-atmospheric specific humidity distribution: a dry atmosphere enhances CTRC by reducing downward thermal radiation. Consequently, the cooling minimizes in the “wet” tropics and maximizes in the “dry” eastern subtropics. Poleward of 30°N/S, the CTRC decreases slightly due to the colder clouds that emit less effectively. Furthermore, the CTRC exhibits distinctive seasonal cycles with stronger cooling in the winter and has amplitudes of order 10–20 Wm -2 in stratocumulus-rich regions. The datasets were used to train a machine-learning model that substantially speeds up the retrieval.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Interface-induced sign reversal of the anomalous Hall effect in magnetic topological insulator heterostructures

The Berry phase picture provides important insights into the electronic properties of condensed matter systems. The intrinsic anomalous Hall (AH) effect can be understood as the consequence of non-zero Berry curvature in momentum space. Here, we fabricate TI/magnetic TI heterostructures and find that the sign of the AH effect in the magnetic TI layer can be changed from being positive to negative with increasing the thickness of the top TI layer. Our first-principles calculations show that the built-in electric fields at the TI/magnetic TI interface influence the band structure of the magnetic TI layer, and thus lead to a reconstruction of the Berry curvature in the heterostructure samples. Based on the interface-induced AH effect with a negative sign in TI/V-doped TI bilayer structures, we create an artificial “topological Hall effect”-like feature in the Hall trace of the V-doped TI/TI/Cr-doped TI sandwich heterostructures. Our study provides a new route to create the Berry curvature change in magnetic topological materials that may lead to potential technological applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

CMDV (CM)4 Project - University of Washington contribution. Final report

The original goal of this project was to diagnose and improve CLUBB, the turbulence and cloud fraction parameterization used in DOE’s E3SM model. For this purpose, we planned to use large-eddy simulation (LES) and ARM observations from the NE Pacific Ocean and the SGP site, in coordination with ARM’s LASSO program. It was determined that for these cloud regimes, many of the turbulent ‘moment closures’ which underlie CLUBB’s mathematical formulation are inconsistent with LES, which is an appropriate benchmark for testing this. The research assistant found that these closures could be more accurately formulated using machine learning using the LES as a training dataset. However the resulting parameterization proved to quickly drift away from physical plausibility. A novel machine-learning based boundary layer parameterization called MARBLE based on matching the time evolution of a parameterized cloud-topped boundary layer to reanalysis was then developed and published.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Feasibility of main thermal ion heating by ICRF waves using a top launcher in a tokamak with deuterium–tritium plasmas

A scenario of ion cyclotron range of frequency (ICRF) wave injection from a top launcher is proposed as an efficient and direct heating method for thermal deuterium ions in deuterium–tritium tokamak plasmas. Positioned between the tritium cyclotron layer and ion–ion hybrid layer, the top launcher allows effective wave penetration to the ion–ion hybrid layer and enables significant power transfer to thermal deuterium. This is achieved through favorable wave polarization for fundamental cyclotron damping. There is a Doppler broadening around the cyclotron resonance and this overlaps with the ion–ion hybrid layer. Low toroidal mode numbers and ion temperature in the range of 5–20 keV are favorable for enhancing the main ion damping relative to electron damping. In contrast to the neutral beam injection, which penetration strongly depends on machine size and plasma density, the proposed ICRF-based direct ion heating scenario is shown to be scalable and applicable to both larger and smaller tokamak devices within practical constraints.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Opportunistic Short‐Term Water Uptake Dynamics by Subalpine Trees Observed via In Situ Water Isotope Measurements

Abstract Variations in tree water sources are important to understand in semi‐arid ecosystems because climatic shifts towards lower snowpack and increased drought affect water availability in subalpine forests of the western US. Here, we use daily in situ measurements of stable isotopes ( 2 H & 18 O) in soil and tree stem water, soil matric potential and sap flow to study tree water uptake dynamics. We instrumented three soil profiles down to 90 cm, as well as three aspen and engelmann spruce trees near Gothic, Colorado, in the East River watershed. We observed the fate of natural isotopic variations in rainfall, soil, and plants from June to October 2022, and in August 2023 we conducted a 2 H labeled irrigation experiment. Our observations showed that all studied aspen trees compensated for water scarcity in the shallow soil by shifting the dominant water source at 60(±20) cm to ⅔ of uptake from 90 cm within a few days of a dry period. Both species relied on snowmelt stored in the subsoil to sustain transpiration. Intense rainfall caused the plant water uptake to shift partially to top soil layers within 2 days. Spruce transpiration was lower and relied more on snowmelt, because rainfall infiltration was low in the spruce stand due to high canopy interception. Our findings highlight the important role of snowmelt stored in the deep soil layers for subalpine forest drought response and the dominant fate of monsoonal rainfall to become transpiration rather than recharging groundwater and streams in the Upper Colorado River. Plain Language Summary There is a need to understand how trees in mountainous regions respond to dry conditions that lead to water scarcity, because climate projections suggest that such conditions will become more frequent in the future. Here we present a novel data set of measurements of daily stable isotopes of water across soil profiles and in tree stems of aspen and spruce. Our data show that when the upper soil dried out, aspen trees shifted to using water from deeper layers (beneath 60 cm) to keep transpiring. For spruce trees the uptake pattern is less clear, but both types of trees mainly used snowmelt stored in the deeper soil layers to survive the dry summer. After heavy rain, aspen and spruce trees switched to using water from the top 20 cm of soil. However, for spruce, only some rain reached the soil because the dense tree canopy intercepted it, so spruce trees stayed more dependent on snowmelt and used less water overall. This study shows how important deep snowmelt water is for helping forests survive dry periods and suggests that most summer rain is quickly used by trees rather than replenishing streams and groundwater in the headwaters of the Colorado River. Key Points Tree water resources changed within a few days from snow dominated to higher share of rainfall as soils wetted up after a dry period Compensatory plant water uptake by aspen from the deep layer (90 cm), while uptake from soil depths that became drier (60 cm) declined Strong differences between water sources and availability beneath aspen and spruce, respectively

Sprenger, Matthias↗
DOE OSTI2020

Elevated temperature mechanical properties of Inconel 617 surface oxide using nanoindentation

Inconel 617 is a principal candidate material for helium gas cooled very-high-temperature reactors with outlet temperatures of 700–950 °C. Recent findings showed that this alloy develops unique surface oxide layers especially at high temperature (HT) helium environment with distinctive wear, friction and contact properties. This study investigates the elevated temperature mechanical properties of Inconel 617 top surface layers aged in HT helium environment. Nanoindentation technique is used to obtain load-displacement graphs of the alloy top surface oxide in temperatures ranging from 25 °C up to 600 °C. In addition, using finite element analysis along with an iterative regression technique, a semi-numerical method is developed to further measure and quantify the material parameters and, in particular, time-independent and creep characteristics of the oxide. Although Young's modulus of the oxide is found to be relatively close to the bulk for the tested temperatures, the yield strength and hardness, in comparison to the bulk material, increase significantly as the material is oxidized after aging. The oxide exhibits significant softening as the temperature increases to 600 °C. Unlike the bulk material, diffusion through the grains is found to be the dominant creep mechanism for the oxide. Considerable difference between the mechanical properties of the oxide and the bulk material shows the need for accurate measurements of near surface mechanical properties, if reliable predictive contact and tribological models are sought at elevated temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Electrically conductive carbon fiber layers as lightning strike protection for non-conductive epoxy-based CFRP substrate

A large amount of electrically conductive fillers is needed to enhance a Carbon Fiber Reinforced Plastics (CFRP) electrical conductivity enough to withstand lightning strikes of peak currents. However, such large alien constituents hamper the inherent good mechanical properties of CFRP structures. In this work, a solution has been proposed to retain both desired properties in a CFRP laminate. Layer-wise hybrid laminate has been demonstrated as a solution for lightning strike protection of Carbon Fiber Reinforced Plastics (CFRP). Top few layers of a hybrid laminate are prepared using electrically conductive polymer-based resin (CF/C-POLY) to provide effective dissipation of lightning current while epoxy-based CFRP substrate (CF/Epoxy) provides the main structural strength. An insulating adhesive layer is used to bond CF/C-POLY and CF/Epoxy to prepare the laminate. The hybrid laminates were tested for their effectiveness against lightning strikes. Laminates were struck by modified lightning waveform of component A with peak current of -14 kA and -40 kA. The performance of the laminates against lightning strike were evaluated using high speed camera, high-speed and thermal camera. It is found that CF/C-POLY layer successfully defended the main structural component i.e. CF/Epoxy from lightning direct damage.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fabrication of low blaze angle gratings by replication and plasma etch

We suggest a new method of making ultra-low blaze angle gratings for synchrotron application. The method is based on reduction of the blaze angle of a master grating by replication followed by a plasma etch. A master blazed grating with a relatively large blaze angle is fabricated by anisotropic wet etching of a Si single crystal substrate. The surface of the master grating is replicated by a polymer material on top of a quartz substrate by nanoimprinting and then transferred into quartz by a plasma etch. Then a 2nd nanoimprint step is applied to transfer the saw-tooth surface into a resist layer on top of a Si grating substrate. The plasma etch through the patterned resist layer provides transfer of the grooves into the Si substrate and results in reduction of the blaze angle due to the difference in etch rates of the resist and Si. We investigated the impact of the replication process on the groove shape, facet surface roughness, and diffraction efficiency of the fabricated 200 lines/mm low blaze angle grating.

Park, Sooyeon↗
DOE OSTI2023

Field Testing of Self-Healing Metallic Coatings for Internal Corrosion Protection of Natural Gas Pipelines

Internal corrosion occurs in natural gas pipelines primarily due to the presence of water, carbon dioxide, and hydrogen sulfide. Internal corrosion can eventually result in leakage, cracks, and rupture of the pipeline. The objective of this work is to mitigate internal corrosion in steel pipelines transporting natural gas using cold spray coatings. The corrosion behavior of carbon steel coated with self-healing zinc chromium (ZnCr) and zinc niobium (ZnNb) cold spray coatings was investigated in a natural gas environment. For comparison purposes, hot-dip galvanized steel (HDGS) was tested under the same conditions as cold spray coatings. The field test was performed at the NW Natural gas storage facility in Mist, Oregon. The field test was conducted in a 6-inch diameter pipe transporting wet natural gas out of an underground storage well at 500 psi and 4.4 °C. The coupons were tested for 15 and 32 days under stagnant and flow conditions, respectively. Weight loss method was used to measure the corrosion rate of metallic coatings. Post-corrosion surface characterization was performed on the specimen using a scanning electron microscopy (SEM) equipped with energy dispersive X-ray spectroscopy (EDS). Crystalline phases were determined by X-ray diffraction (XRD). The field-test results confirmed that the metallic coatings provided corrosion protection of carbon steel exposed to wet natural gas under stagnant and flow conditions. The formation of the ZnCO 3 layer on top of ZnNb and ZnCr coatings led to passivation of the coatings which helps to reduce the self-corrosion. These layers form a barrier for diffusion of corrosive species to the surface.

03 NATURAL GAS↗
DOE OSTI2022

Origin of Ferroelectric Phase Stabilization via the Clamping Effect in Ferroelectric Hafnium Zirconium Oxide Thin Films

The presence of the top electrode on hafnium oxide-based thin films during processing has been shown to drive an increase in the amount of metastable ferroelectric orthorhombic phase and polarization performance. This “Clamping Effect,” also referred to as the Capping or Confinement Effect, is attributed to the mechanical stress and confinement from the top electrode layer. However, other contributions to orthorhombic phase stabilization have been experimentally reported, which may also be affected by the presence of a top electrode. In this study, it is shown that the presence of the top electrode during thermal processing results in larger tensile biaxial stress magnitudes and concomitant increases in ferroelectric phase fraction and polarization response, whereas film chemistry, microstructure, and crystallization temperature are not affected. Through etching experiments and measurement of stress evolution for each processing step, it is shown that the top electrode locally inhibits out-of-plane expansion in the HZO during crystallization, which prevents equilibrium monoclinic phase formation and stabilizes the orthorhombic phase. This study provides a mechanistic understanding of the clamping effect and orthorhombic phase formation in ferroelectric hafnium oxide-based thin films, which informs the future design of these materials to maximize ferroelectric phase purity and corresponding polarization behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Maximized Hole Trapping in a Polystyrene Transistor Dielectric from a Highly Branched Iminobis(aminoarene) Side Chain

We synthesized highly branched and electron donating side chain subunits and attached them to polystyrene (PS) used as a dielectric layer in a pentacene field effect transistor. The influence of these groups on dielectric function, charge retention and threshold voltage shifts (ΔV th ) depending on their positions in dielectric multilayers was determined. We compared the observations made on an N-perphenylated iminobisaniline side chain with those from the same side chains modified with ZnO nanoparticles and with an adduct formed from tetracyanoethylene (TCNE). We also synthesized an analogue in which six methoxy groups are present instead of two amine nitrogens. At 6 mol % side chain, hopping transport was sufficient to cause shorting of the gate, while at 2 mol %, charge trapping was observable as transistor threshold voltage shifts (ΔV th ). Here, we created three types of devices: with the substituted PS layer as single layer dielectric, on top of a crosslinked PS layer but in contact with the pentacene (bilayers), and sandwiched between two PS layers in trilayers. Particularly large bias stress effects and ΔV th , larger than the hexamethoxy and previously studied dimethoxy analogs, were observed in the second case and the effects increased with the increasing electron donating properties of the modified side chains. The highest ΔV th was consistent with a majority of the side chains stabilizing trapped charge. Trilayer devices showed decreased charge storage capability compared to previous work in which we used less donating side chains but in higher concentrations. The ZnO and TCNE modifications resulted in slightly more and less negative ΔV th , respectively, when the side chain polystyrene was not in contact with the pentacene and isolated from the gate electrode. The results indicate a likely maximum combination of molecular charge stabilizing activity and side chain concentration that still allows gate dielectric function.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Enhanced lightning strike protection using vertically oriented carbon fiber melded with conventional carbon fiber-reinforced composite and its validation through damage analysis

Lightning strike protection (LSP) is one of the prime factors in aerospace and wind industries concerning safety, enhanced service life, and reduced downtime. To ensure better LSP, high electrical conductivity is required to dissipate the current, which is always a challenge for polymer composites due to the inherent insulating nature of the polymer matrix. Conventional carbon fiber-reinforced composites (CFRP) offer electrical conductivity in the planar direction while the vertical i.e. through-thickness conductivity still remains a challenge. Having been motivated by this fact, we have fabricated CFRP interleaved with vertically oriented CF (Z-fiber) for the lightning strike test (100 kA). A paint was applied to mimic the actual service condition of the composite laminates. We have prepared two Z-fiber composites: (a) Z-1 containing one Z-fiber layer on the top along with conventional CF layers and (b) Z-5 containing five interleaved Z-fiber layers and compared the properties with conventional CFRP (Z-0). Interestingly, even a single Z-fiber layer (Z-1) showed lower damage (visual damage diameter 22 mm) than Z-0 (visual damage diameter 26 mm). Z-5 displayed the lowest damage (visual damage diameter 16.7 mm), which was further confirmed by cross-sectional optical microscopy. Improved LSP was attributed to the higher through-thickness electrical conductivity in the case of Z-5 (9 times with respect to Z-0) reflecting a lower electrical anisotropy for Z-fiber composites. The residual mechanical property after the lightning test was analyzed through the flexural test, and the retention of flexural strength and modulus was 66% and 86%, respectively for Z-5 showing a significant improvement in comparison to Z-0 (>40% for both).

Saha, Subhabrata↗
DOE OSTI2023

Advanced low blaze angle x-ray gratings via nanoimprint replication and plasma etch

We developed a new method of making ultra-low blaze angle diffraction gratings for x-ray applications. The method is based on reduction of the blaze angle of a master grating by nanoimprint replication followed by a plasma etch. A master blazed grating with a relatively large blaze angle is fabricated by anisotropic wet etching of a Si single crystal substrate. The surface of the master grating is replicated by a polymer material on top of a quartz substrate by nanoimprinting. Then a second nanoimprinting is performed using the 1st replica as a mold to replicate the saw-tooth surface into a resist layer on top of a Si grating substrate. A reactive ion etch is used to transfer the grating grooves into the Si substrate. The plasma etch provides reduction of the groove depth by a factor defined by the ratio of the etch rates for the resist and Si. We demonstrate reduction of the blaze angle of a master grating by a factor of 5 during fabrication of a 200 lines/mm blazed grating with a blaze angle of 0.2°. We investigated the quality and performance of the fabricated low blaze angle gratings and evaluate process accuracy and reproducibility. The new blaze angle reduction method preserves the planarity of the optical surface of the grating substrate and at the same time provides improvement in the grating groove quality during the reduction process.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Heat treatment to anneal residual stresses during additive manufacturing

The present disclosure relates to a method of producing a product through additive manufacturing with heat treatment. The method involves using a fusing beam to melt powder particles disposed on a substrate, where the fusing beam is impressed with a two dimensional pattern containing image information from a first layer to be printed. The fused powder particles are then heat treated with a beam impressed with an additional two dimensional pattern. The additional two dimensional pattern has image information from the first layer to achieve heat treatment of the product. The heat treatment is completed prior to laying down additional new layers of material. The heat treatment is an annealing operation. The method may further involve providing a new layer of powdered material on top of the layer of fused powder particles subsequent to the heat treatment, and repeating the melting and heat treating operations in a layer-by-layer fashion using the two dimensional pattern and the additional two dimensional pattern, until the part is completed.

DeMuth, James A.↗
DOE OSTI2024

Improving the estimate of higher-order moments from lidar observations near the top of the convective boundary layer

Abstract. Ground-based lidar data have proven extremely useful for profiling the convective boundary layer (CBL). Many groups have derived higher-order moments (e.g., variance, skewness, fluxes) from high-temporal-resolution lidar data using an autocovariance approach. However, these analyses are highly uncertain near the CBL top when the depth of the CBL (zi) is changing during the analysis period. This is because the autocovariance approach is usually applied to constant height levels and the character of the eddies is changing on either side of the changing CBL top. Here, a new approach is presented wherein the autocovariance analysis is performed on a normalized height grid, with a temporally smoothed zi. Output from a large eddy simulation model demonstrates that deriving higher-order moments from time series on a normalized height grid has better agreement with the slab-averaged quantities than the moments derived from the original height grid.

Rosenberger, Tessa E. (ORCID:0000000333205873)↗
DOE OSTI2020

X-ray detection of ultrashort spin current pulses in synthetic antiferromagnets

Here, we explore the ultrafast generation of spin currents in magnetic multilayer samples by applying fs laser pulses to one layer and measuring the magnetic response in the other layer by element-resolved x-ray spectroscopy. In Ni(5 nm)/Ru(2 nm)/Fe(4 nm), the Ni and Fe magnetization directions couple antiferromagnetically due to the Ruderman–Kittel–Kasuya–Yosida interaction but may be oriented parallel through an applied magnetic field. After exciting the top Ni layer with a fs laser pulse, we also find that the Fe layer underneath demagnetizes, with a 4.1±1.9% amplitude difference between parallel and antiparallel orientation of the Ni and Fe magnetizations. We attribute this difference to the influence of a spin current generated by the fs laser pulse that transfers angular momentum from the Ni into the Fe layer. Our results confirm that superdiffusive spin transport plays a role in determining the sub-ps demagnetization dynamics of synthetic antiferromagnetic layers, but also evidence large depolarization effects due to hot electron dynamics, which are independent of the relative alignment of the magnetization in Ni and Fe.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A Model Intercomparison Study to Investigate Mixing Characteristics in Non‐Precipitating Stratocumulus Clouds

Recent aircraft observations of marine stratocumulus clouds consistently showed that cloud microphysical relationships vary with altitude, indicating inhomogeneous mixing characteristics near cloud top and homogeneous mixing characteristics in mid-levels of clouds. Here, we conduct model intercomparison of an idealized, non-precipitating stratocumulus cloud to evaluate model consistency and examine whether simulations can reproduce the observed mixing characteristics. The results show that eleven large-eddy simulations with various dynamics and microphysics schemes show good agreement on the thermodynamical, microphysical, and dynamical properties of the stratocumulus-topped boundary layer in a steady state. The inter-model spread in steady-state liquid water path is significantly reduced compared to previous model intercomparison studies. This improvement might be due to better models and more consistent initial conditions than those used decades ago. In addition, most simulations, including a low-dimensional simulation, capture inhomogeneous mixing characteristics near the cloud top and homogeneous mixing characteristics inside the cloud. Moreover, simulations using Lagrangian microphysics schemes agree better with the observed mixing characteristics compared with those using the bin microphysics schemes. Since most simulations do not fully resolve the entrainment process, the apparent mixing characteristics arise from the variations in the resolved cloud properties. Our results support the vertical circulation mixing hypothesis, which suggests that homogeneous mixing characteristics in mid-levels of clouds are due to the vertical circulation of entrainment-affected and diluted parcels from the cloud top moved to lower levels.

54 ENVIRONMENTAL SCIENCES↗
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↗