Search NASA⌕ Search

SEARCH · Search NASA

Results for “reactive transport model”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Using Reactive Transport Modeling to Understand Formation of the Stimson Sedimentary Unit and Altered Fracture Zones at Gale Crater, Mars

Water flowing through sediments at Gale Crater, Mars created environments that were likely habitable, and sampled basin-wide hydrological systems. However, many questions remain about these environments and the fluids that generated them. Measurements taken by the Mars Science Laboratory Curiosity of multiple fracture zones can help constrain the environments that formed them because they can be compared to nearby associated parent material (Figure 1). For example, measurements of altered fracture zones from the target Greenhorn in the Stimson sandstone can be compared to parent material measured in the nearby Big Sky target, allowing constraints to be placed on the alteration conditions that formed the Greenhorn target from the Big Sky target. Similarly, CheMin measurements of the powdered < 150 micron fraction from the drillhole at Big Sky and sample from the Rocknest eolian deposit indicate that the mineralogies are strikingly similar. The main differences are the presence of olivine in the Rocknest eolian deposit, which is absent in the Big Sky target, and the presence of far more abundant Fe oxides in the Big Sky target. Quantifying the changes between the Big Sky target and the Rocknest eolian deposit can therefore help us understand the diagenetic changes that occurred forming the Stimson sedimentary unit. In order to interpret these aqueous changes, we performed reactive transport modeling of 1) the formation of the Big Sky target from a Rocknest eolian deposit-like parent material, and 2) the formation of the Greenhorn target from the Big Sky target. This work allows us to test the relationships between the targets and the characteristics of the aqueous conditions that formed the Greenhorn target from the Big Sky target, and the Big Sky target from a Rocknest eolian deposit-like parent material.

Hausrath, E. M.↗
NASA NTRS2019

Interpreting Aqueous Alteration in the Murray Formation Using Reactive Transport Modeling

Abundant evidence for liquid water exists at Gale crater, Mars. However, the characteristics of past water remain an area of active research. The first exposures of the Murray formation in Gale crater, Mars (Fig. 1) were studied with four samples analyzed using CheMin: Buckskin, Telegraph Peak, Mojave, and Confidence Hills. Analyses indicate differences in mineralogy and chemistry between the samples which have been attributed to changes in pH and oxidation state of depositional and diagenetic environments. Recent work also suggests that hydrothermal fluids may have been present based on the presence of Se, Zn, Pb, and other elements.

Hausrath, E. M.↗
NASA NTRS2017

Investigating Weathering of Basaltic Materials in Gale Crater, Mars: A Combined Laboratory, Modeling and Terrestrial Field Approach

Recent observations from Gale Crater, Mars document past aqueous alteration both in the formation of the Stimson sandstone unit, as well as in the formation of altered fractures within that unit. Geochemical and mineralogical data from Curiosity also suggest Fe-rich amorphous weathering products are present in most samples measured to date. Here we interpret conditions of possible past weathering in Gale Crater using a combination of field, laboratory, and modeling work. In order to better understand secondary Fe-rich phases on Mars, we are examining formation of weathering products in high Fe and Mg and low Al serpentine soils in the Klamath Mountains, CA. We have isolated potential weathering products from these soils, and are analyzing them using synchrotron μXRF and μXRD as well as FullPat for a direct comparison to analyses from Gale Crater. In order to interpret the implications of the persistence of potential secondary Fe-containing phases on Mars, we are also measuring the dissolution rates of the secondary weathering products allophane, Fe-rich allophane, and hisingerite. Ongoing dissolution experiments of these materials suggest that they dissolve significantly more rapidly than more crystalline secondary minerals with similar chemical compositions. Finally, to quantify the specific conditions of past aqueous alteration in Gale Crater we are performing reactive transport modeling of a range of possible past environmental conditions. Specifically, we are testing the conditions under which a Stimson unit-like material forms from a parent material similar to Rocknest or Bagnold eolian deposits, and the conditions under which observed altered fracture zones form from a Stimson unit-like parent material. Our modeling results indicate that the formation of the Stimson unit is consistent with leaching of an eolian deposit with a solution of pH = 6-8, and that formation of the altered fracture zones is consistent with leaching with a very acidic (pH = 2-3) high sulfate solution containing Ca. These results suggest circumneutral pH conditions during authigenesis or early diagenesis in the Stimson formation sediments followed by diagenetic alteration by very acidic solutions along fracture zones.

Hausrath, Elisabeth↗
NASA NTRS2016

Biological Potential in Serpentinizing Systems

Generation of the microbial substrate hydrogen during serpentinization, the aqueous alteration of ultramafic rocks, has focused interest on the potential of serpentinizing systems to support biological communities or even the origin of life. However the process also generates considerable alkalinity, a challenge to life, and both pH and hydrogen concentrations vary widely across natural systems as a result of different host rock and fluid composition and differing physical and hydrogeologic conditions. Biological potential is expected to vary in concert. We examined the impact of such variability on the bioenergetics of an example metabolism, methanogenesis, using a cell-scale reactive transport model to compare rates of metabolic energy generation as a function of physicochemical environment. Potential rates vary over more than 5 orders of magnitude, including bioenergetically non-viable conditions, across the range of naturally occurring conditions. In parallel, we assayed rates of hydrogen metabolism in wells associated with the actively serpentinizing Coast Range Ophiolite, which includes conditions more alkaline and considerably less reducing than is typical of serpentinizing systems. Hydrogen metabolism is observed at pH approaching 12 but, consistent with the model predictions, biological methanogenesis is not observed.

the potential of serpentinizing systems to support↗
NASA NTRS2017

Energy in a Planetary Context

The potential present day habitability of solar system bodies beyond Earth is limited to subsurface environments, where the availability of energy in biologically useful form is a paramount consideration. Energy availability is commonly quantified in terms of molar Gibbs energy changes for metabolisms of interest, but this can provide an incomplete and even misleading picture. A second aspect of life's requirement for energy - the rate of delivery, or power - strongly influences habitability, biomass abundance, growth rates, and, ultimately, rates of evolution. We are developing an approach to quantify metabolic power, using a cell-scale reactive transport model in which physical and chemical environmental parameters are varied. Simultaneously, we evaluate cell-specific energy flux requirements and their dependence on environmental "extremes". Comparison of metabolic power supply and demand provides a constraint on how biomass abundance varies across a range of environmental parameters, and thereby a prediction of the relative habitability of different environments. We are evaluating the predictive capability of this approach through comparison to observed distributions of microbial abundance in a range of subsurface (predominantly serpentinizing) systems.

MARINE GEOMICROBIOLOGY↗
NASA NTRS2018

Reactive Transport and Mass Balance Modeling of the Stimson Sedimentary Formation and Altered Fracture Zones Constrain Diagenetic Conditions at Gale Crater, Mars

On a planet as cold and dry as present-day Mars, evidence of multiple aqueous episodes offers an intriguing view into very different past environments. Fluvial, lacustrine, and eolian depositional environments are being investigated by the Mars Science Laboratory Curiosityin Gale crater, Mars. Geochemical and mineralogical observations of these sedimentary rocks suggest diagenetic processes affected the sediments. Here, we analyze diagenesis of the Stimson formation eolian parent material, which caused loss of olivine and formation of magnetite. Additional, later alteration in fracture zones resulted in preferential dissolution of pyroxene and precipitation of secondary amorphous silica and Ca sulfate. The ability to compare the unaltered parent material with the reacted material allows constraints to be placed on the characteristics of the altering solutions. In this work we use a combination of a mass balance approach calculating the fraction of a mobile element lost or gained, τ, with fundamental geochemical kinetics and thermodynamics in the reactive transport code CrunchFlow to examine the characteristics of multiple stages of aqueous alteration at Gale crater, Mars. Our model results indicate that early diagenesis of the Stimson sedimentary formation is consistent with leaching of an eolian deposit by a near-neutral solution, and that formation of the altered fracture zones is consistent with a very acidic, high sulfate solution containing Ca, P and Si. These results indicate a range of past aqueous conditions occurring at Gale crater, Mars, with important implications for past Martian climate and environments.

Gale crater↗
NASA NTRS2018

Reactive Transport and Mass Balance Modeling of the Stimson Sedimentary Formation and Altered Fracture Zones Constrain Diagenetic Conditions at Gale Crater, Mars

On a planet as cold and dry as present-day Mars, evidence of multiple aqueous episodes offers an intriguing view into very different past environments. Fluvial, lacustrine, and eolian depositional environments are being investigated by the Mars Science Laboratory Curiosity in Gale crater, Mars. Geochemical and mineralogical observations of these sedimentary rocks suggest diagenetic processes affected the sediments. Here, we analyze diagenesis of the Stimson formation eolian parent material, which caused loss of olivine and formation of magnetite. Additional, later alteration in fracture zones resulted in preferential dissolution of pyroxene and precipitation of secondary amorphous silica and Ca sulfate. The ability to compare the unaltered parent material with the reacted material allows constraints to be placed on the characteristics of the altering solutions. In this work we use a combination of a mass balance approach calculating the fraction of a mobile element lost or gained, tau, with fundamental geochemical kinetics and thermodynamics in the reactive transport code CrunchFlow to examine the characteristics of multiple stages of aqueous alteration at Gale crater, Mars. Our model results indicate that early diagenesis of the Stimson sedimentary formation is consistent with leaching of an eolian deposit by a near-neutral solution, and that formation of the altered fracture zones is consistent with a very acidic, high sulfate solution containing Ca, P and Si. These results indicate a range of past aqueous conditions occurring at Gale crater, Mars, with important implications for past martian climate and environments.

Gale↗
NASA NTRS2012

Progress Toward Affordable High Fidelity Combustion Simulations Using Filtered Density Functions for Hypersonic Flows in Complex Geometries

Significant progress has been made in the development of subgrid scale (SGS) closures based on a filtered density function (FDF) for large eddy simulations (LES) of turbulent reacting flows. The FDF is the counterpart of the probability density function (PDF) method, which has proven effective in Reynolds averaged simulations (RAS). However, while systematic progress is being made advancing the FDF models for relatively simple flows and lab-scale flames, the application of these methods in complex geometries and high speed, wall-bounded flows with shocks remains a challenge. The key difficulties are the significant computational cost associated with solving the FDF transport equation and numerically stiff finite rate chemistry. For LES/FDF methods to make a more significant impact in practical applications a pragmatic approach must be taken that significantly reduces the computational cost while maintaining high modeling fidelity. An example of one such ongoing effort is at the NASA Langley Research Center, where the first generation FDF models, namely the scalar filtered mass density function (SFMDF) are being implemented into VULCAN, a production-quality RAS and LES solver widely used for design of high speed propulsion flowpaths. This effort leverages internal and external collaborations to reduce the overall computational cost of high fidelity simulations in VULCAN by: implementing high order methods that allow reduction in the total number of computational cells without loss in accuracy; implementing first generation of high fidelity scalar PDF/FDF models applicable to high-speed compressible flows; coupling RAS/PDF and LES/FDF into a hybrid framework to efficiently and accurately model the effects of combustion in the vicinity of the walls; developing efficient Lagrangian particle tracking algorithms to support robust solutions of the FDF equations for high speed flows; and utilizing finite rate chemistry parametrization, such as flamelet models, to reduce the number of transported reactive species and remove numerical stiffness. This paper briefly introduces the SFMDF model (highlighting key benefits and challenges), and discusses particle tracking for flows with shocks, the hybrid coupled RAS/PDF and LES/FDF model, flamelet generated manifolds (FGM) model, and the Irregularly Portioned Lagrangian Monte Carlo Finite Difference (IPLMCFD) methodology for scalable simulation of high-speed reacting compressible flows.

Drozda, Tomasz G.↗
NASA NTRS2018

Chemical Modeling of the Reactivity of Short-Lived Greenhouse Gases: A Model Inter-Comparison Prescribing a Well-Measured, Remote Troposphere

We develop a new protocol for merging in situ measurements with 3-D model simulations of atmospheric chemistry with the goal of integrating over the data to identify the most reactive air parcels in terms of tropospheric production and loss of the greenhouse gases ozone and methane. Presupposing that we can accurately measure atmospheric composition, we examine whether models constrained by such measurements agree on the chemical budgets for ozone and methane. In applying our technique to a synthetic data stream of 14,880 parcels along 180W, we are able to isolate the performance of the photochemical modules operating within their global chemistry-climate and chemistry-transport models, removing the effects of modules controlling tracer transport, emissions, and scavenging. Differences in reactivity across models are driven only by the chemical mechanism and the diurnal cycle of photolysis rates, which are driven in turn by temperature, water vapor, solar zenith angle, clouds, and possibly aerosols and overhead ozone, which are calculated in each model. We evaluate six global models and identify their differences and similarities in simulating the chemistry through a range of innovative diagnostics. All models agree that the more highly reactive parcels dominate the chemistry (e.g., the hottest 10% of parcels control 25-30% of the total reactivities), but do not fully agree on which parcels comprise the top 10%. Distinct differences in specific features occur, including the regions of maximum ozone production and methane loss, as well as in the relationship between photolysis and these reactivities. Unique, possibly aberrant, features are identified for each model, providing a benchmark for photochemical module development. Among the 6 models tested here, 3 are almost indistinguishable based on the inherent variability caused by clouds, and thus we identify 4, effectively distinct, chemical models. Based on this work, we suggest that water vapor differences in model simulations of past and future atmospheres may be a cause of the different evolution of tropospheric O3 and CH4, and lead to different chemistry-climate feedbacks across the models.

greenhouse gases↗
NASA NTRS1991

The influence of polar heterogeneous processes on reactive chlorine at middle latitudes - Three dimensional model implications

Three-dimensional model calculations with the NASA/GSFC chemistry and transport model have been designed to consider the impact of heterogeneous processes occurring on polar stratospheric clouds (PSCs) in the Arctic vortex on the HCl distribution. By examining the HCl concentration for a calculation with PSCs relative to a calculation with gas phase chemistry only, the impact of polar processing on reactive chlorine species at middle latitudes is inferred. Results from the chemistry and transport model reproduce basic features of the ClO measurements (Toohey et al., 1991), which were made on the ferry flights of the ER-2 from Stavanger, Norway to Moffett Field, California via Wallops Island, Virginia on February 20 and 21, 1989. The model indicates that perturbed air which is contained within the polar vortex during winter is not homogeneously mixed, and that the ferry flights were made through air with the largest conversion of HCl to reactive chlorine that is seen at middle latitudes.

Douglass, Anne R.↗
NASA NTRS1989

Model-free simulations of turbulent reactive flows

The current computational methods for solving transport equations of turbulent reacting single-phase flows are critically reviewed, with primary attention given to those methods that lead to model-free simulations. In particular, consideration is given to direct numerical simulations using spectral (Galerkin) and pseudospectral (collocation) methods, spectral element methods, and Lagrangian methods. The discussion also covers large eddy simulations and turbulence modeling.

Givi, Peyman↗
NASA NTRS1990

Diagnostic studies of the HxOy-NxOy-O3 photochemical system using data from NASA GTE field expeditions

The research effort supported in part by the subject grant focused on three related topics. Our major effort was concentrated on the analysis of data gathered during GTE field expeditions. Ancillary efforts were directed at: the development and application of a Global Chemical Transport Model for the study of the atmospheric reactive nitrogen budget; the development and application of a one-dimensional, time dependent cloud model for the study of the impact of in-cloud aqueous phase chemistry on the atmospheric sulfur budget; and mechanistic studies of the chemical processes involved in dry deposition of ozone to vegetative surfaces. In the sections below, we briefly summarize the central conclusions of each of these efforts. These discussions are followed by a listing of the papers completed during the granting period and the graduate students supported by funds from the grant. Reprints and preprints of all papers completed with support from the grant are attached as appendices.

Chameides, W. L.↗
NASA NTRS2000

Evolution of Bromoform in a Global Chemistry and Transport Model

It is well known that many chlorine and bromine compounds that are inert in the troposphere are destroyed in the stratosphere and contribute to the stratospheric burden of reactive chlorine and bromine species. But the contribution from those chlorine and bromine compounds which are reactive in the troposphere is less certain because it is not known whether convection can transport these gases to the upper troposphere rapidly enough to overcome their short tropospheric lifetimes. We examine this issue using a three-dimensional chemistry and transport model to simulate the evolution of three gases which have surface sources, bromoform (CHBr3), methyl chloroform (CH3CCl3), and carbon dioxide (CO2). Our objective is to determine if CHBr3 might enhance the lower stratospheric burden of reactive bromine. The other two gases provide tests of the quality of the simulation. Both CHBr3 and CH3CCl3 are destroyed in the troposphere by reaction with hydroxyl (OH), whose latitudinal and monthly variation is provided by a two-dimensional model and upon which a diurnal variation is imposed. Comparison of the lifetime of CH3CCl3 computed from observations (5 years) with the lifetime computed from the simulation provides an integrated test of the model's transport and photochemistry. Observations also show that CO2 exhibits a strong seasonal cycle in the northern hemisphere troposphere that is not propagated directly across the tropopause into the lower stratosphere. Thus, maintenance of the observed troposphere-stratosphere distinctness of CO2 in the presence of convection is a critical benchmark for meeting our objective.

Douglass, Anne R.↗
NASA NTRS1976

Crystal growth from the vapor phase experiment MA-085

Three vapor transport experiments on multicomponent systems were performed during the Apollo Soyuz mission to determine the effects of microgravity forces on crystal morphology and mass transport rates. The mixed systems used germanium selenide, tellurium, germanium tetraiodide (transport agent), germanium monosulfide, germanium tetrachloride (transport agent), and argon (inert atmosphere). The materials were enclosed in evacuated sealed ampoules of fused silica and were transported in a temperature gradient of the multipurpose electric furnace onboard the Apollo Soyuz spacecraft. Preliminary evaluation of 2 systems shows improved quality of space grown crystals in terms of growth morphology and bulk perfection. This conclusion is based on a direct comparison of space grown and ground based crystals by means of X-ray diffraction, microscopic, and chemical etching techniques. The observation of greater mass transport rates than predicted for a microgravity environment by existing vapor transport models indicates the existence of nongravity caused transport effects in a reactive solid/gas phase system.

Wiedemeir, H.↗
NASA NTRS1972

Transport phenomena of reactive fluid flow in heterogeneous combustion processes.

A previously developed computer program was used to model two transient hybrid combustion processes involving tubes of solid Plexiglas. In the first study, representing combustion of a hybrid rocket, the oxidizing gas was oxygen, and calculations were continued sufficiently long to obtain steady-state values. Systematic variations were made in reaction rate constant, mass flow rate, and pressure, alternatively using constant and temperature dependent regression rate models for the fuel surface. Consistent results were obtained, as is evidenced by the values for the mass function of the reaction product and the flame temperature, for which plots are supplied. In the second study, fire initiation in a duct was studied, with an air mixture as the oxidizing gas. It was demonstrated that a satisfactory flame spread mechanism could be reproduced on the computer. In both of the above applications, the general, transient, two-dimensional conservation equations were represented, together with chemical reactions, solid-fuel interface conditions, and heat conduction in the solid fuel.

Hung, W. S. Y.↗
NASA NTRS1988

Calculation of the recirculating compressible flow downstream a sudden axisymmetric expansion

Significant progress has been made during the last five years to adapt conventional Navier-Stokes solver for handling nonconservative equations. A primary type of application is to use transport equation turbulence models, but the extension is also possible for describing the transport of nonpassive scalars, such as in reactive media. Among others, combustion and gas dissociation phenomena are topics needing a considerable research effort. An implicit two step scheme based on the well-known MacCormack scheme has been modified to treat compressible turbulent flows on complex geometries. Implicit treatment of nonconservative equations (in the present case a two-equation turbulence model) opens the way to the coupled solution of thermochemical transport equations.

Vandromme, D.↗
NASA NTRS2021

On the Cycling of 231 Pa and 230 Th in Benthic Nepheloid Layers

The naturally-occurring radionuclides protactinium-231 ( 231 Pa) and thorium-230 ( 230 Th) are produced at approximately uniform rates in the ocean and thought to be removed from the water column through a reversible exchange with settling particles. Recent measurements along the U.S. GEOTRACES North Atlantic transect (GA03) revealed two features which are at odds with current understanding about 231 Pa and 230 Th cycling in the ocean: (i) a sharp decrease in dissolved 231 Pa ( 231 Pa d ) and 230 Th ( 230 Th d ) activities with depth below 2000-4000 m and (ii) very high particulate 231 Pa ( 231 Pa p ) and 230 Th ( 230 Th p ) activities near the bottom, at a number of stations between the New England continental shelf and Bermuda. Concomitant measurements of light attenuation from beam transmissometry showed that both features occur in benthic nepheloid layers (BNLs), which suggests that these features may stem, at least partly, from the presence of resuspended sediment in the deep water column. Here we explore the behaviour of 231 Pa and 230 Th in BNLs by using (i) radionuclide, optical, and hydrographic data from the western segment of GA03 (west of Bermuda) and (ii) a simplified model of particle and radionuclide cycling that includes a lateral particle source. First, the BNLs observed at GA03 stations are characterized from measurements of the beam attenuation coefficient converted to particle concentrations. At all stations, particle concentrations below the clear water minimum were the highest in the bottom mixed layer, whose thickness ranged from 95 to 320 m, and decreased generally with height above the bottom. The thickness of strong BNLs varied from 482 to 1358 m and the vertical integral of particle concentration in excess to that at the clear water minimum varied from 1 x 10 4 to 2 x 10 6 mg m −2 , among different stations. Second, the particle-radionuclide model is fitted to data from stations GT11-04 (New England continental rise) and GT11-08 (Hatteras abyssal plain), where samples for radionuclide analyses were collected in the BNL. The model can reproduce simultaneously the increase of particle concentration with depth, the low 231 Pa d and 230 Th d in the BNLs, and the high 231 Pa p and 230 Th p near the bottom. According to the model, at heights less than about 300 m above the seafloor, the dissolved phase was set primarily by a balance between adsorption and desorption, with vertical turbulent mixing playing a secondary role, whilst the particulate phase behaved largely as a non-reactive constituent supplied laterally and transported vertically by particle settling and turbulent mixing. Sensitivity tests with the model suggest that lateral particle sources near continental slopes and similar reliefs can produce significant biases both in the 230 Th normalization method and in the interpretation of sediment 231 Pa/ 230 Th records. Our findings yield insights into the influence of sediment resuspension and transport on 231 Pa and 230 Th in the deep ocean and highlight the need for considering these processes in paleoceanographic applications.

Nepheloid layer↗