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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS2003

Pneumatically Actuated Miniature Peristaltic Vacuum Pumps

Pneumatically actuated miniature peristaltic vacuum pumps have been proposed for incorporation into advanced miniature versions of scientific instruments that depend on vacuum for proper operation. These pumps are expected to be capable of reaching vacuum-side pressures in the torr to millitorr range (from .133 down to .0.13 Pa). Vacuum pumps that operate in this range are often denoted roughing pumps. In comparison with previously available roughing pumps, these pumps are expected to be an order of magnitude less massive and less power-hungry. In addition, they would be extremely robust, and would operate with little or no maintenance and without need for oil or other lubricants. Portable mass spectrometers are typical examples of instruments that could incorporate the proposed pumps. In addition, the proposed pumps could be used as roughing pumps in general laboratory applications in which low pumping rates could be tolerated. The proposed pumps could be designed and fabricated in conventionally machined and micromachined versions. A typical micromachined version (see figure) would include a rigid glass, metal, or plastic substrate and two layers of silicone rubber. The bottom silicone layer would contain shallow pump channels covered by silicone arches that could be pushed down pneumatically to block the channels. The bottom silicone layer would be covered with a thin layer of material with very low gas permeability, and would be bonded to the substrate everywhere except in the channel areas. The top silicone layer would be attached to the bottom silicone layer and would contain pneumatic- actuation channels that would lie crosswise to the pump channels. This version is said to be micromachined because the two silicone layers containing the channels would be fabricated by casting silicone rubber on micromachined silicon molds. The pneumatic-actuation channels would be alternately connected to a compressed gas and (depending on pump design) either to atmospheric pressure or to a partial vacuum source. The design would be such that the higher pneumatic pressure would be sufficient to push the silicone arches down onto the substrates, blocking the channels. Thus, by connecting pneumatic- actuation channels to the two pneumatic sources in spatial and temporal alternation, waves of opening and closing, equivalent to peristalsis, could be made to move along the pump channels. A pump according to this concept could be manufactured inexpensively. Pneumatic sources (compressors and partial vacuum sources) similar those needed for actuation are commercially available; they typically have masses of .100 g and power demands of the order of several W. In a design-optimization effort, it should be possible to reduce masses and power demands below even these low levels and to integrate pneumatic sources along with the proposed pumps into miniature units with overall dimensions of no more than a few centimeters per side.

Feldman, Sabrina↗
NASA NTRS1977

Mathematical models for the reflection coefficients of lossy dielectric half-spaces with application to transient responses of chirped pulses

Reflection coefficients are found at normal incidence for a large class of homogeneous lossy half-spaces with a one-dimensionally inhomogeneous or stratified lossy layer on top. Solutions are in terms of Hankel functions of complex argument to decrease cancellation error at high frequencies. One special case is that of layers on a homogeneous half-space where the dielectric constant in each layer may vary in a quite general manner. A Wronskian is used to insure the critical computations are correct. The reflection of chirped pulses is considered. Solutions are obtained by applying the fast Fourier transform. It is found that for a typical relatively long normalized 'long' pulse the power reflected as a function of time is essentially the power reflection coefficient for the frequencies swept out, whereas for a relatively short 'long' pulse, with the same relative change in frequency and the same number of oscillations there is only the uniform attenuation by the power reflection coefficient of the center frequency. By a 'long' pulse we mean a pulse whose spatial length is long compared to the thickness of the reflecting layer.

Evans, D. D.↗
NASA NTRS1980

Entrainment into a stratocumulus layer with distributed radiative cooling

It is shown that the radiative cooling of a cloud layer strongly influences the turbulent flux profiles and the entrainment rate, and that the radiative cooling should be modeled as acting inside the turbulent layer. Numerical experiments demonstrate that a cloud-topped mixed-layer model, similar to that of Lilly (1968), is quite sensitive to delta (p) sub R, the depth of the radiatively cooled layer near cloud top. As delta (p) sub R increases, the model's sensitivity to the entrainment assumption is markedly heightened; for large delta (p) sub R the cloud top and cloud base rise as the entrainment parameter k is increased, while for small delta (p) sub R an increase in k has almost no effect. The model is most sensitive to delta (p) sub R for the cold-water, strong-divergence regime of greatest interest.

Randall, D. A.↗
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↗
NASA NTRS1980

Deep convection on day 261 of GATE

The structural features of the deep convection observed on September 18, 1974, day 261 of the GARP Atlantic Tropical Experiment (GATE), as the ridge axis of a 700 mb wave passed the center of the GATE B-scale network are reported. Satellite and aircraft maps indicate the presence of clouds penetrating above 2.5 km into the middle troposphere organized in bands about 9 km apart and aligned roughly along the direction of the wind shear in the cloud layer. Radar echoes corresponding to cumulus convection of lifetime, peak height and peak rainfall rates on the orders of 30 min, 6 km and 1.3 mm/h, respectively, were observed to triple in number density as convergence at 950 hPa increased from 1.5 to 3 x 10 to the -5th/sec. The structural features of the radar echoes indicate that the day was similar to a mesoscale precipitation feature of Leary and Houe (1979), with the cluster consisting of many echoes appearing in succession. Data from aircraft penetrations of the deep convection reveal downdrafts accompanying the precipitation and updrafts immediately to their south. Shipboard and rawinsonde observations show that the convective downdrafts brought down air of low pseudo-equivalent potential temperature, with local surface convergence of up to 0.001/sec. Mean wind shears through the cloud layer to the top of a main cloud layer are found to be only 75% greater than those of Malkus (1958) for the Caribbean, with shears just above the cloud base several factors larger.

Warner, C.↗
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)↗
NASA NTRS2002

TRMM Radar Observations of Cloud Tops in the Tropical Tropopause Layer

Air dehydrates to stratospheric abundances in the tropical tropopause layer (TTL). The role of overshooting convection in the dehydration process is not well understood. To study this effect, we use the TRMM (Tropical Rainfall Measuring Mission) precipitation radar (PR) to measure the altitudes of cloud tops forming in the TTL. Because the radar signal is dominated by scatter from large particles, these cloud observations imply the presence of strong convective systems with large updraft. Both winter and summer data from two different years are examined to study both interseasonal and interannual variability. The global distribution of these clouds is in good agreement with those of the surface precipitation rates. In addition, the altitude distributions of these clouds follow an exponential dependence. However, clouds over continental regions typically extend to higher altitudes in the tropics. Almost no cloud tops were observed above 20 km. Comparison between the radar cloud tops and colocated IR brightness temperature measurements reveal a large difference in both the diurnal cycle and intensity between continental and oceanic convection.

Alcala, C. M.↗
NASA NTRS2024

Optimizing Thin Cover Layers for Optically Clear Electrodynamic Dust Shields

This project looked at optimizing the dielectric cover layer on top of an optically clear glass electrodynamic dust shield (EDS). The EDS is a device that uses high voltage electric fields to remove dust from surfaces. Currently, the EDS uses 2000 volts to remove dust, but by optimizing the cover layer and the geometry of the EDS, this voltage can be lowered to improve safety and reduce electromagnetic interference (EMI) and the complexity of supporting systems. This project was able to demonstrate successful dust clearing on an EDS with 750 volts in high vacuum (~1E-5 torr) conditions. The EDS had a trace spacing of 500 μm, and the Parylene C cover layer had a thickness of 25 μm. There is still a substantial amount of room for improvement of these devices, which may be able to push the operating voltage even lower.

Particulate Contamination↗
NASA NTRS2001

Toward a Unified Parameterization of the Boundary Layer and Moist Convection: Simulations of Clear and Cloudy Convection - Part 3

A model that employs a new form of mass-flux closure (described in Part I of this paper) is applied to a variety of clear and cloudy planetary boundary layers (PBLs) including dry convection from the Wangara Experiment, trade wind cumulus from the Barbados Oceanographic and Meteorological Experiment (BOMEX), and marine stratocumulus from the Atlantic Stratocumulus Experiment (ASTEX). For Wangara, the simulated variances and fluxes match that expected from similarity arguments, while the mean state is a little less mixed than the observations. In the BOMEX simulation. the shape and magnitude of the fluxes and the turbulence kinetic energy budget agree with LES results and observations. However, the liquid water mixing ratio is too large. This is attributed to an underprediction of the skewness. In agreement with observations from the ASTEX experiment, many of the model-simulated fields distinctly reflect a regime in transition between the trade wind cumulus and the classic stratocumulus-topped boundary layers. In general, the simulated entrainment rate tends to be a little underpredicted in regimes where there is little cloud-top radiative cooling (Wangara and BOMEX), while it is overpredicted in regimes where this process is more critical (e.g., ASTEX). Prior work suggests that this may be related to the manner in which the pressure terms are parameterized in the model. Overall, the model is able to capture some key physical features of these PBL regimes, and appears to have the potential to represent both cloud and boundary layer processes. Thus, this approach is a first step toward unifying these processes in large-scale models.

Lappen, Cara-Lyn↗
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↗
NASA NTRS1985

Three-layer Generic Ganymedes: Structure and Evolution

Accurate, three layered structural models were generated for Ganymede and Callisto. Three layered satellites consist of a rock core, a region of mixed rock and ice, and an outer shell of pure ice. This structure would result from either accretional melting of the outermost region of the satellite or differentiation subsequent to localized melting. A completely undifferentiated ice rock satellite possesses only the single, top boundary layer and, in this case, melting is also relatively difficult to initiate; this suggests that, if accretional melting results in a small degree of differentiation initially, subsequent melting may readily occur as the satellite warms due to radiogenic heating. Alternatively, if accretion is not accompanied by a small degree of differentiation, it may prove too difficult to initiate nonaccretional melting and the satellite might remain undifferentiated indefinitely.

Mueller, S. W.↗
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↗