Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sparse sampling”

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 703 records · Page 39

Assessing entropy for catalytic processes at complex reactive interfaces

When chemical reactions are accelerated by a catalyst, entropy differences between reactants and their transient intermediates can be the driving force behind the promotion or inhibition of desired and parasitic chemical pathways. Understanding and controlling catalytic processes therefore requires both a fundamental and practicable understanding of entropy in addition to enthalpy. In unstructured media such as the vapor phase equilibrated with sparsely covered surfaces, entropy can be adequately accounted for by well-established approaches based on translational, rotational, and harmonic vibrational partition functions. However, these approximations become inadequate in more complex condensed phase environments, e.g., solid liquid interfaces of confined reaction spaces. In this chapter, we provide an overview of the state-of-art in the computational quantification of entropy and its known ramifications on catalysis. The fundamental roles of thermodynamics and kinetics in catalysis are covered in enough detail to appreciate and contextualize the computational methods employed to compute chemically accurate estimates of entropy. These methods are discussed in appropriate detail and range from the ubiquitous harmonic oscillator approximation where entropy unrelated to high frequency oscillations is typically underestimated, to enhanced free energy sampling with molecular dynamics where the desired accuracy must be weighed against the associated computational cost of obtaining it. The rising importance of machine learning and artificial intelligence in accelerating methodological progress in this field is touched upon, as well. Finally, applications, successes, and pitfalls of using these methods are provided to showcase past and present accomplishments while clarifying where improvements in both understanding and methodology are still needed.

Kollias, Loukas↗

Joint Retrieval of Surface BRDF from Geostationary and Polar-Orbiting Satellite Sensors

The latest geostationary sensors like GOES 16/17 ABI and Himawari 8/9 AHI provide high frequent observations of the Earth surface with continuously changing solar illumination geometries, which allow us to retrieve the surface Bidirectional Reflectance Distribution Function (BRDF) with dedicated atmospheric correction algorithms like MAIAC (Multi-Angle Implementation of Atmospheric Correction). However, because the viewing geometry of a specific location from the geostationary satellites are fixed, the angular sampling of surface BRDF by GEO (Geostationary Earth Orbit) sensors is far from comprehensive. This study tries to address this issue by exploring a GEO-LEO (Low-Earth-Orbit) synergy, in particular, jointly retrieving surface BRDF parameters with concurrent ABI/AHI and VIIRS top-of-atmosphere (TOA) reflectance for the near-infrared (NIR) band. The NIR band is chosen because the ABI, AHI, and VIIRS instruments have very similar spectral response functions in this band and therefore simplifies the requirements for cross-sensor radiometric calibration. We compile ABI/AHI and VIIRS TOA data with the best quality (cloud-free and low aerosol loading) at chosen AERONET sites. We then run the GeoNEX MAIAC algorithm to retrieve the Ross-Thick-Li-Sparse (RTLS) surface BRDF parameters with or without the AEORNET measured atmospheric aerosol optical depth (AOD) as inputs. The joint retrieval results are considered the best estimate of surface BRDF. We compare the joint BRDF retrievals with the corresponding MAIAC BRDF products, retrieved with ABI/AHI or VIIRS separately, to evaluate their differences. We expect that the jointly retrieved BRDF data are more robust than the standard products and may help us reduce uncertainties in higher-level earth observation satellite products.

Remote Sensing↗

Sparse signal reconstruction from polychromatic X-ray CT measurements via mass attenuation discretization

We propose a method for reconstructing sparse images from polychromatic x-ray computed tomography (ct) measurements via mass attenuation coefficient discretization. The material of the inspected object and the incident spectrum are assumed to be unknown. We rewrite the Lambert-Beer’s law in terms of integral expressions of mass attenuation and discretize the resulting integrals. We then present a penalized constrained least-squares optimization approach for reconstructing the underlying object from log-domain measurements, where an active set approach is employed to estimate incident energy density parameters and the nonnegativity and sparsity of the image density map are imposed using negative-energy and smooth ℓ{sub 1}-norm penalty terms. We propose a two-step scheme for refining the mass attenuation discretization grid by using higher sampling rate over the range with higher photon energy, and eliminating the discretization points that have little effect on accuracy of the forward projection model. This refinement allows us to successfully handle the characteristic lines (Dirac impulses) in the incident energy density spectrum. We compare the proposed method with the standard filtered backprojection, which ignores the polychromatic nature of the measurements and sparsity of the image density map. Numerical simulations using both realistic simulated and real x-ray ct data are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

High temperature oxidation behavior of Ni-based superalloy GH202

The isothermal oxidation kinetic and oxidation behavior of Ni-based superalloy GH202 were investigated at temperatures from 800 °C to 1100 °C. The microstructure of GH202 was characterized by scanning electron microscope (SEM) and electron back-scattered diffraction (EBSD). The composition and morphologies of oxide scales of the surface and cross-section samples were examined by X-ray diffraction (XRD), electron probe microanalysis (EPMA), and SEM combined with energy-dispersive spectroscopy (EDS). The results show that the growth of the oxide scales on GH202 obeyed a parabolic law with an activation energy of 247.8 KJ·mol{sup −1} from 800 °C to 1100 °C. The continuous and dense oxide scales that formed at 800 °C and 900 °C were composed of Cr{sub 2}O{sub 3}, NiCr{sub 2}O{sub 4} and a small amount of TiO{sub 2}. Beneath the oxide scales, aluminum-rich oxides formed along the grain boundaries due to internal oxidation. After treatment at 1000 °C for 150 h, the oxide scale was separated into two layers, the outer layer with sparse holes that was mainly composed of Cr{sub 2}O{sub 3} and TiO{sub 2} and an inner layer of dense Cr{sub 2}O{sub 3}. At 1100 °C for 100 h the outer layer was spalled, due to the looser oxide scale with large blocks of TiO{sub 2} and NiCr{sub 2}O{sub 4} which provided channels for oxygen diffusion. - Highlights: •Isothermal oxidation kinetic of GH202 obeys a parabolic law from 800 °C to 1100 °C. •The activation energy of GH202 is closed to the diffusion activation of Cr{sup 3+}. •The evolution process of oxide scales from 800 °C to 1100 °C is displayed by EPMA. •The oxidation behavior of GH202 is mainly controlled by Cr{sup 3+} and Ti{sup 4+}.

36 MATERIALS SCIENCE↗

Statistics of Wide Pre-main-sequence Binaries in the Orion OB1 Association

Statistics of low-mass pre-main-sequence binaries in the Orion OB1 association with separations ranging from 0."6 to 20" (220 to 7400 au at 370 pc) are studied using images from the VISTA Orion mini survey and astrometry from Gaia. The input sample based on the CVSO catalog contains 1137 stars of K and M spectral types (masses between 0.3 and 0.9 M ⊙ ), 1021 of which are considered to be association members. There are 135 physical binary companions to these stars with mass ratios above ~0.13. The average companion fraction is 0.09 ± 0.01 over 1.2 decades in separation, slightly less than, but still consistent with, the field. We found a difference between the Ori OB1a and OB1b groups, the latter being richer in binaries by a factor of 1.6 ± 0.3. No overall dependence of the wide-binary frequency on the observed underlying stellar density is found, although in the Ori OB1a off-cloud population, these binaries seem to avoid dense clusters. The multiplicity rates in Ori OB1 and in sparse regions like Taurus differ significantly, hinting that binaries in the field may originate from a mixture of diverse populations.

79 ASTRONOMY AND ASTROPHYSICS↗

Characterizing the Influence of Fracture Density on Network Scale Transport

The topology of natural fracture networks is inherently linked to the structure of the fluid velocity field and transport therein. Here we study the impact of network density on flow and transport behaviors. We stochastically generate fracture networks of varying density and simulate flow and transport with a discrete fracture network model, which fully resolves network topology at the fracture scale. We study conservative solute trajectories with Lagrangian particle tracking and find that as fracture density decreases, solute channelization to large local fractures increases, thereby reducing plume spreading. Furthermore, in sparse networks mean particle travel distance increases and local network features, such as velocity zones where flow is counter to the primary pressure gradient, become increasingly important for transport. As the network density increases, network statistics homogenize and such local features have a reduced impact. We quantify local topological influence on transport behavior with an effective tortuosity parameter, which measures the ratio of total advective distance to linear distance at the fracture scale; large tortuosity values are correlated to slow-velocity regions. Additionally, these large tortuosity, slow-velocity regions delay downstream transport and enhance tailing on particle breakthrough curves. Finally, we predict transport with an upscaled, Bernoulli spatial Markov random walk model and parameterize local topological influences with a novel tortuosity parameter. Furthermore, Bernoulli model predictions improve when sampling from a tortuosity distribution, as opposed to a fixed value as has previously been done, suggesting that local network topological features must be carefully considered in upscaled modeling efforts of fracture network systems.

58 GEOSCIENCES↗

The role of topography in geodetic gravity field modelling

Masses associated with the topography, bathymetry, and its isostatic compensation are a dominant source of gravity field variations, especially at shorter wavelengths. On global scales the topographic/isostatic effects are also significant, except for the lowest harmonics. In practice, though, global effects need not be taken into account as such effects are included in the coefficients of the geopotential reference fields. On local scales, the short-wavelength gravity variations due to the topography may, in rugged terrain, be an order of magnitude larger than other effects. In such cases, explicit or implicit terrain reduction procedures are mandatory in order to obtain good prediction results. Such effects may be computed by space-domain integration or by fast Fourier transformation (FFT) methods. Numerical examples are given for areas of the Canadian Rockies. In principle, good knowledge of the topographic densities is required to produce the smoothest residual field. Densities may be determined from sample measurements or by gravimetric means, but both are somewhat troublesome methods in practice. The use of a standard density, e.g., 2.67 g/cu cm, may often yield satisfactory results and may be put within a consistent theoretical framework. The independence of density assumptions is the key point of the classical Molodensky approach to the geodetic boundary value problem. The Molodensky solutions take into account that land gravity field observations are done on a non-level surface. Molodensky's problem may be solved by integral expansions or more effective FFT methods, but the solution should not be intermixed with the use of terrain reductions. The methods are actually complimentary and may both be required in order to obtain the smoothest possible signal, least prone to aliasing and other effects coming from sparse data coverage, typical of rugged topography.

Forsberg, R.↗

Inspection design using 2D phased array, TFM and cueMAP software

A simulation suite, cueMAP, has been developed to facilitate the design of inspection processes and sparse 2D array configurations. At the core of cueMAP is a Total Focusing Method (TFM) imaging algorithm that enables computer assisted design of ultrasonic inspection scenarios, including the design of bespoke array configurations to match the inspection criteria. This in-house developed TFM code allows for interactive evaluation of image quality indicators of ultrasonic imaging performance when utilizing a 2D phased array working in FMC/TFM mode. The cueMAP software uses a series of TFM images to build a map of resolution, contrast and sensitivity of imaging performance of a simulated reflector, swept across the inspection volume. The software takes into account probe properties, wedge or water standoff, and effects of specimen curvature. In the validation process of this new software package, two 2D arrays have been evaluated on 304n stainless steel samples, typical of the primary circuit in nuclear plants. Thick section samples have been inspected using a 1MHz 2D matrix array. Due to the processing efficiency of the software, the data collected from these array configurations has been used to investigate the influence sub-aperture operation on inspection performance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

CyberGAN: Generating High-fidelity Cybersecurity Data With Generative Adversarial Networks

Machine learning for cyber defense offers the promise of detecting adversarial activity against the ground data systems managing critical space assets. A fundamental challenge facing machine learning research in cybersecurity is the lack of high-fidelity, shareable datasets for robust evaluation and testing of machine learning-based solutions. High-fidelity, real-world datasets are necessary for reliable benchmarking of nominal system behavior and malicious activity. Unfortunately, such realistic datasets of both nominal and adversarial activity are rarely shared publicly by data owners due to security and privacy concerns. Besides, the available adversarial data is sparse, which makes training models on malicious activity much harder. This situation has impeded and continues to impede the research and successful adoption of machine learning methods for cyber defense. Researchers have dealt with this problem by generating data within a low-fidelity lab environment, using classified and thus unshareable datasets, or downloading low-fidelity public datasets made available by others. We propose an innovative solution to the problem by employing machine learning methods to generate high-fidelity data. Specifically, we propose the use of Generative Adversarial Networks (GANs) to generate high-fidelity data for cybersecurity purposes. GANs have found successful image processing and natural language applications, but have not yet been investigated for cyber data generation. Our proposed approach first involves training the `discriminator' network of the GAN with a sample of real-world data consisting of malicious and nominal samples. We then use the `generator' network to generate new high-fidelity data samples consisting of an appropriate mix of malicious and nominal activity. We demonstrate applications of our architecture by generating high-fidelity cybersecurity data containing both malicious and nominal samples. We thoroughly evaluate the fidelity of our generated data using heuristics and evaluate its usefulness for machine learning applications using three different datasets. Overall, our approach results in high-fidelity, shareable datasets.

Zhang, Yuening↗

Improving signal strength in serial crystallography with DIALS geometry refinement

The DIALS diffraction-modeling software package has been applied to serial crystallography data. Diffraction modeling is an exercise in determining the experimental parameters, such as incident beam wavelength, crystal unit cell and orientation, and detector geometry, that are most consistent with the observed positions of Bragg spots. These parameters can be refined by nonlinear least-squares fitting. In previous work, it has been challenging to refine both the positions of the sensors (metrology) on multipanel imaging detectors such as the CSPAD and the orientations of all of the crystals studied. Since the optimal models for metrology and crystal orientation are interdependent, alternate cycles of panel refinement and crystal refinement have been required. To simplify the process, a sparse linear algebra technique for solving the normal equations was implemented, allowing the detector panels to be refined simultaneously against the diffraction from thousands of crystals with excellent computational performance. Separately, it is shown how to refine the metrology of a second CSPAD detector, positioned at a distance of 2.5 m from the crystal, used for recording low-angle reflections. With the ability to jointly refine the detector position against the ensemble of all crystals used for structure determination, it is shown that ensemble refinement greatly reduces the apparent nonisomorphism that is often observed in the unit-cell distributions from still-shot serial crystallography. In addition, it is shown that batching the images by timestamp and re-refining the detector position can realistically model small, time-dependent variations in detector position relative to the sample, and thereby improve the integrated structure-factor intensity signal and heavy-atom anomalous peak heights.

36 MATERIALS SCIENCE↗

NASA Tech Briefs, May 2010

Topics covered include: Instrument for Analysis of Greenland's Glacier Mills Cryogenic Moisture Apparatus; A Transportable Gravity Gradiometer Based on Atom Interferometry; Three Methods of Detection of Hydrazines; Crossed, Small-Deflection Energy Analyzer for Wind/Temperature Spectrometer; Wavefront Correction for Large, Flexible Antenna Reflector; Novel Micro Strip-to-Waveguide Feed Employing a Double-Y Junction; Thin-Film Ferro Electric-Coupled Microstripline Phase Shifters With Reduced Device Hysteresis; Two-Stage, 90-GHz, Low-Noise Amplifier; A 311-GHz Fundamental Oscillator Using InP HBT Technology; FPGA Coprocessor Design for an Onboard Multi-Angle Spectro-Polarimetric Imager; Serrating Nozzle Surfaces for Complete Transfer of Droplets; Turbomolecular Pumps for Holding Gases in Open Containers; Triaxial Swirl Injector Element for Liquid-Fueled Engines; Integrated Budget Office Toolbox; PLOT3D Export Tool for Tecplot; Math Description Engine Software Development Kit; Astronaut Office Scheduling System Software; ISS Solar Array Management; Probabilistic Structural Analysis Program; SPOT Program; Integrated Hybrid System Architecture for Risk Analysis; System for Packaging Planetary Samples for Return to Earth; Offset Compound Gear Drive; Low-Dead-Volume Inlet for Vacuum Chamber; Simple Check Valves for Microfluidic Devices; A Capillary-Based Static Phase Separator for Highly Variable Wetting Conditions; Gimballing Spacecraft Thruster; Finned Carbon-Carbon Heat Pipe with Potassium Working Fluid; Lightweight Heat Pipes Made from Magnesium; Ceramic Rail-Race Ball Bearings; Improved OTEC System for a Submarine Robot; Reflector Surface Error Compensation in Dual-Reflector Antennas; Enriched Storable Oxidizers for Rocket Engines; Planar Submillimeter-Wave Mixer Technology with Integrated Antenna; Widely Tunable Mode-Hop-Free External-Cavity Quantum Cascade Laser; Non-Geiger-Mode Single-Photon Avalanche Detector with Low Excess Noise; Using Whispering-Gallery-Mode Resonators for Refractometry; RF Device for Acquiring Images of the Human Body; Reactive Collision Avoidance Algorithm; Fast Solution in Sparse LDA for Binary Classification; Modeling Common-Sense Decisions in Artificial Intelligence; Graph-Based Path-Planning for Titan Balloons; Nanolaminate Membranes as Cylindrical Telescope Reflectors; Air-Sea Spray Airborne Radar Profiler Characterizes Energy Fluxes in Hurricanes; Large Telescope Segmented Primary Mirror Alignment; and Simplified Night Sky Display System.

Source record↗

UQ Toolkit v. 3.0

The Uncertainty Quantification (UQ) Toolkit is a software library for the characterization and propagation of uncertainties in computational models. This library provides Bayesian inference tools to infer uncertain model parameters, Bayesian compressive sensing methods for discovering sparse representations of high-dimensional input-output response surfaces, methods for constructing Karhunen-Loève representations of stochastic processes, and global sensitivity analysis tools used to compute Sobol indices in order to characterize the importance of uncertain inputs or the interactions between them. The basis for many of these methods relies on representing random variables with Polynomial Chaos expansions (PCEs). This library implements several spectral basis function types (e.g. Hermite basis functions in terms of Gaussian random variables or Legendre basis functions in terms of uniform random variables) that can be used to represent random variables with PCEs. For the propagation of uncertainty, tools are provided to propagate PCEs that describe the input uncertainty through the computational model using either intrusive methods (Galerkin projection of equations onto basis functions) or non-intrusive methods (perform deterministic operation at sampled values of the random values and project the obtained results onto basis functions).

Sargsyan, Khachik↗

THE RUNAWAYS AND ISOLATED O-TYPE STAR SPECTROSCOPIC SURVEY OF THE SMC (RIOTS4)

We present the Runaways and Isolated O-Type Star Spectroscopic Survey of the SMC (RIOTS4), a spatially complete survey of uniformly selected field OB stars that covers the entire star-forming body of the Small Magellanic Cloud (SMC). Using the IMACS (Inamori-Magellan Areal Camera and Spectrograph) multislit spectrograph and MIKE (Magellan Inamori Kyocera Echelle) echelle spectrograph on the Magellan telescopes, we obtained spectra of 374 early-type field stars that are at least 28 pc from any other OB candidates. We also obtained spectra of an additional 23 field stars in the SMC bar identified from slightly different photometric criteria. Here, we present the observational catalog of stars in the RIOTS4 survey, including spectral classifications and radial velocities. For three multi-slit fields covering 8% of our sample, we carried out monitoring observations over 9–16 epochs to study binarity, finding a spectroscopic, massive binary frequency of at least ∼60% in this subsample. Classical Oe/Be stars represent a large fraction of RIOTS4 (42%), occurring at much higher frequency than in the Galaxy, consistent with expectation at low metallicity. RIOTS4 confirmed a steep upper initial mass function in the field, apparently caused by the inability of the most massive stars to form in the smallest clusters. Our survey also yields evidence for in situ field OB star formation, and properties of field emission-line star populations, including sgB[e] stars and classical Oe/Be stars. We also discuss the radial velocity distribution and its relation to SMC kinematics and runaway stars. RIOTS4 presents a first quantitative characterization of field OB stars in an external galaxy, including the contributions of sparse, but normal, star formation; runaway stars; and candidate isolated star formation.

79 ASTRONOMY AND ASTROPHYSICS↗

Methodology for sensitivity analysis, approximate analysis, and design optimization in CFD for multidisciplinary applications

In this study involving advanced fluid flow codes, an incremental iterative formulation (also known as the delta or correction form) together with the well-known spatially-split approximate factorization algorithm, is presented for solving the very large sparse systems of linear equations which are associated with aerodynamic sensitivity analysis. For smaller 2D problems, a direct method can be applied to solve these linear equations in either the standard or the incremental form, in which case the two are equivalent. Iterative methods are needed for larger 2D and future 3D applications, however, because direct methods require much more computer memory than is currently available. Iterative methods for solving these equations in the standard form are generally unsatisfactory due to an ill-conditioning of the coefficient matrix; this problem can be overcome when these equations are cast in the incremental form. These and other benefits are discussed. The methodology is successfully implemented and tested in 2D using an upwind, cell-centered, finite volume formulation applied to the thin-layer Navier-Stokes equations. Results are presented for two sample airfoil problems: (1) subsonic low Reynolds number laminar flow; and (2) transonic high Reynolds number turbulent flow.

Taylor, Arthur C., III↗

Assessments of epistemic uncertainty using Gaussian stochastic weight averaging for fluid-flow regression

Here, we use Gaussian stochastic weight averaging (SWAG) to assess the epistemic uncertainty associated with neural-network-based function approximation relevant to fluid flows. SWAG approximates a posterior Gaussian distribution of each weight, given training data, and a constant learning rate. Having access to this distribution, it is able to create multiple models with various combinations of sampled weights, which can be used to obtain ensemble predictions. The average of such an ensemble can be regarded as the 'mean estimation', whereas its standard deviation can be used to construct 'confidence intervals', which enable us to perform uncertainty quantification (UQ) with regard to the training process of neural networks. We utilize representative neural-network-based function approximation tasks for the following cases: (i) a two-dimensional circular-cylinder wake; (ii) the DayMET dataset (maximum daily temperature in North America); (iii) a three-dimensional square-cylinder wake; and (iv) urban flow, to assess the generalizability of the present idea for a wide range of complex datasets. SWAG-based UQ can be applied regardless of the network architecture, and therefore, we demonstrate the applicability of the method for two types of neural networks: (i) global field reconstruction from sparse sensors by combining convolutional neural network (CNN) and multi-layer perceptron (MLP); and (ii) far-field state estimation from sectional data with two-dimensional CNN. We find that SWAG can obtain physically-interpretable confidence-interval estimates from the perspective of epistemic uncertainty. This capability supports its use for a wide range of problems in science and engineering.

97 MATHEMATICS AND COMPUTING↗

Construction of ground-state preserving sparse lattice models for predictive materials simulations

First-principles based cluster expansion models are the dominant approach in ab initio thermodynamics of crystalline mixtures enabling the prediction of phase diagrams and novel ground states. However, despite recent advances, the construction of accurate models still requires a careful and time-consuming manual parameter tuning process for ground-state preservation, since this property is not guaranteed by default. In this paper, we present a systematic and mathematically sound method to obtain cluster expansion models that are guaranteed to preserve the ground states of their reference data. The method builds on the recently introduced compressive sensing paradigm for cluster expansion and employs quadratic programming to impose constraints on the model parameters. The robustness of our methodology is illustrated for two lithium transition metal oxides with relevance for Li-ion battery cathodes, i.e., Li 2x Fe 2(1-x) O 2 and Li 2x Ti 2(1-x) O 2 , for which the construction of cluster expansion models with compressive sensing alone has proven to be challenging. We demonstrate that our method not only guarantees ground-state preservation on the set of reference structures used for the model construction, but also show that out-of-sample ground-state preservation up to relatively large supercell size is achievable through a rapidly converging iterative refinement. This method provides a general tool for building robust, compressed and constrained physical models with predictive power.

36 MATERIALS SCIENCE↗

A Probabilistic Autoencoder for Type Ia Supernova Spectral Time Series

We construct a physically parameterized probabilistic autoencoder (PAE) to learn the intrinsic diversity of Type Ia supernovae (SNe Ia) from a sparse set of spectral time series. The PAE is a two-stage generative model, composed of an autoencoder that is interpreted probabilistically after training using a normalizing flow. We demonstrate that the PAE learns a low-dimensional latent space that captures the nonlinear range of features that exists within the population and can accurately model the spectral evolution of SNe Ia across the full range of wavelength and observation times directly from the data. By introducing a correlation penalty term and multistage training setup alongside our physically parameterized network, we show that intrinsic and extrinsic modes of variability can be separated during training, removing the need for the additional models to perform magnitude standardization. We then use our PAE in a number of downstream tasks on SNe Ia for increasingly precise cosmological analyses, including the automatic detection of SN outliers, the generation of samples consistent with the data distribution, and solving the inverse problem in the presence of noisy and incomplete data to constrain cosmological distance measurements. We find that the optimal number of intrinsic model parameters appears to be three, in line with previous studies, and show that we can standardize our test sample of SNe Ia with an rms of 0.091 ± 0.010 mag, which corresponds to 0.074 ± 0.010 mag if peculiar velocity contributions are removed.

79 ASTRONOMY AND ASTROPHYSICS↗

Basalt Reactivity Variability with Reservoir Depth in Supercritical CO2 and Aqueous Phases

Long term storage of CO{sub 2} in geologic formations is currently considered the most attractive option to reduce greenhouse gas emissions while continuing to utilize fossil fuels for energy production. Injected CO{sub 2} is expected to reside as a buoyant water-saturated supercritical fluid in contact with reservoir rock, the caprock system, and related formation waters. As was reported for the first time at the GHGT-9 conference, experiments with basalts demonstrated surprisingly rapid carbonate mineral formation occurring with samples suspended in the scCO{sub 2} phase. Those experiments were limited to a few temperatures and CO{sub 2} pressures representing relatively shallow (1 km) reservoir depths. Because continental flood basalts can extend to depths of 5 km or more, in this paper we extend the earlier results across a pressure-temperature range representative of these greater depths. Different basalt samples, including well cuttings from the borehole used in a pilot-scale basalt sequestration project (Eastern Washington, U.S.) and core samples from the Central Atlantic Magmatic Province (CAMP), were exposed to aqueous solutions in equilibrium with scCO{sub 2} and water-rich scCO{sub 2} at six different pressures and temperatures for select periods of time (30 to 180 days). Conditions corresponding to a shallow injection of CO{sub 2} (7.4 MPa, 34 C) indicate limited reactivity with basalt; surface carbonate precipitates were not easily identified on post-reacted basalt grains. Basalts exposed under identical times appeared increasingly more reacted with simulated depths. Tests, conducted at higher pressures (12.0 MPa) and temperatures (55 C), reveal a wide variety of surface precipitates forming in both fluid phases. Under shallow conditions tiny clusters of aragonite needles began forming in the wet scCO{sub 2} fluid, whereas in the CO{sub 2} saturated water, cation substituted calcite developed thin radiating coatings. Although these types of coatings were sparse, conditions corresponding to deeper depths showed increasing carbonate precipitation. Basalts exposed to aqueous dissolved CO{sub 2} (25.5 MPa, 116 C) for 30 days were coated in tiny nodules of precipitate ({approx}100 {micro}m in diameter) that were identified by micro x-ray diffraction as ankerite, [Ca(Fe,Mg)(CO{sub 3}){sub 2}], a variety of dolomite commonly associated with hydrothermal and metamorphic environments. Surface characterization by SEM revealed well-developed round nodules composed of discrete individual platelets. In contrast, reaction products forming on the basalt in the corresponding wet scCO{sub 2} phase had completely different morphology, appearing in an optical microscope as a surface coating instead of discrete nodules. Examination by SEM revealed layers of discrete platelets forming a cover over a few discrete nodules. Longer test durations (180 days) produced severe iron staining along with minerals structures similar to rhodochrosite and kutnohorite. These preliminary experiments show strong evidence of the faster rate of increase in mineralization reactions taking place in the scCO{sub 2} phase, transformation reactions that are just beginning to be explored in detail.

54 ENVIRONMENTAL SCIENCES↗