Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reactive Flow”

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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

High rate reactive sputtering of MoN(x) coatings

High rate reactive sputtering of MoN(x) films was performed using feedback control of the nitorgen partial pressure. Coatings were made at four different target powers: 2.5, 5.0, 7.5 and 10 kW. No hysteresis was observed in the nitrogen partial pressure vs. flow plot, as is typically seen for the Ti-N system. Four phases were determined by X-ray diffraction: molybdenum, Mo-N solid solution, Beta-Mo2N and gamma-Mo2N. The hardness of the coatings depended upon composition, substrate bias, and target power. The phases present in the hardest films differed depending upon deposition parameters. For example, the Beta-Mo2N phase was hardest (load 25 gf) at 5.0 kW with a value of 3200 kgf/sq mm, whereas the hardest coatings at 10 kW were the gamma-Mo2N phase (3000 kgf/sq mm). The deposition rate generally decreased with increasing nitrogen partial pressure, but there was a range of partial pressures where the rate was relatively constant. At a target power of 5.0 kW, for example, the deposition rates were 3300 A/min for a N2 partial pressure of 0.05 - 1.0 mTorr.

Rudnik, Paul J.↗
NASA NTRS2000

Resin Film Infusion (RFI) Process Modeling for Large Transport Aircraft Wing Structures

This investigation completed the verification of a three-dimensional resin transfer molding/resin film infusion (RTM/RFI) process simulation model. The model incorporates resin flow through an anisotropic carbon fiber preform, cure kinetics of the resin, and heat transfer within the preform/tool assembly. The computer model can predict the flow front location, resin pressure distribution, and thermal profiles in the modeled part. The formulation for the flow model is given using the finite element/control volume (FE/CV) technique based on Darcy's Law of creeping flow through a porous media. The FE/CV technique is a numerically efficient method for finding the flow front location and the fluid pressure. The heat transfer model is based on the three-dimensional, transient heat conduction equation, including heat generation. Boundary conditions include specified temperature and convection. The code was designed with a modular approach so the flow and/or the thermal module may be turned on or off as desired. Both models are solved sequentially in a quasi-steady state fashion. A mesh refinement study was completed on a one-element thick model to determine the recommended size of elements that would result in a converged model for a typical RFI analysis. Guidelines are established for checking the convergence of a model, and the recommended element sizes are listed. Several experiments were conducted and computer simulations of the experiments were run to verify the simulation model. Isothermal, non-reacting flow in a T-stiffened section was simulated to verify the flow module. Predicted infiltration times were within 12-20% of measured times. The predicted pressures were approximately 50% of the measured pressures. A study was performed to attempt to explain the difference in pressures. Non-isothermal experiments with a reactive resin were modeled to verify the thermal module and the resin model. Two panels were manufactured using the RFI process. One was a stepped panel and the other was a panel with two 'T' stiffeners. The difference between the predicted infiltration times and the experimental times was 4% to 23%.

Loos, Alfred C.↗
DOE OSTI2019

Catalytic Hot-Gas Filtration with a Supported Heteropolyacid Catalyst for Preconditioning Biomass Pyrolysis Vapors

During ex situ catalytic fast pyrolysis (CFP) of biomass, the separation of reactive char and alkali/alkaline particulates from biomass pyrolysis vapors by hot-gas filtration (HGF) leads to improved vapor stability and quality. HGF in tandem with chemical tailoring (e.g., partial deoxygenation) of the clean pyrolysis vapors, denoted as catalytic hot-gas filtration (CHGF), has the potential to further improve vapor composition by removing reactive oxygen moieties and protect downstream upgrading catalysts from fouling. Downstream upgrading refers to both vapor phase upgrading (e.g., ex situ CFP) and condensed phase upgrading (e.g., hydrotreating). Consequently, CHGF (as a single unit operation) was evaluated for preconditioning pyrolysis vapors for downstream upgrading processes. In order to understand the effective operating conditions that successfully filter and partially deoxygenate pyrolysis vapors, a titania-supported molybdenum heteropolyacid (Mo-HPA/TiO 2 ) catalyst was studied for use in CHGF. Here, pine pyrolysis vapors were generated in a small pilot-scale pyrolyzer and transferred to a CHGF unit via a continuous-flow slipstream. In the CHGF unit, the pyrolysis vapors were filtered and upgraded over a packed Mo-HPA/TiO 2 catalyst bed. Real-time monitoring and identification of the products formed were achieved by molecular beam mass spectrometry. The results showed that under a hydrogen-rich environment, the pine vapors were partially deoxygenated and alkylated over the Mo-HPA/TiO 2 catalyst. Reactivity studies revealed that an increase in hydrogen concentration and a reduction in weight-hourly space velocity enhanced deoxygenation and alkylation. Time-on-stream (TOS) studies showed stable product formation up to 1 h with little change in catalyst activity. Additionally, the liquid product was collected using a custom fractional condensation unit (built in-house) and analyzed by gas chromatography mass spectrometry to confirm that the product was partially deoxygenated and alkylated. The combination of CHGF and fractional condensation allowed for chemical and physical removal of both foulant and value-added compounds (e.g., phenols, alkylphenols, methoxyphenols, cyclopentenones) for additional enhancement of downstream upgrading processes. The pre- and postreaction catalysts were characterized using temperature-programmed desorption, N 2 physisorption, and elemental analysis with results indicating some catalyst coking. A hydrogen-based catalyst regeneration procedure restored the reacted catalyst activity to that of fresh Mo-HPA/TiO 2 .

09 BIOMASS FUELS↗
DOE OSTI2014

Simulation of Coupled Processes of Flow, Transport, and Storage of CO 2 in Saline Aquifers

This report is the final scientific one for the award DE- FE0000988 entitled “Simulation of Coupled Processes of Flow, Transport, and Storage of CO 2 in Saline Aquifers.” The work has been divided into six tasks. In task, “Development of a Three-Phase Non-Isothermal CO 2 Flow Module,” we developed a fluid property module for brine-CO 2 mixtures designed to handle all possible phase combinations of aqueous phase, sub-critical liquid and gaseous CO 2 , supercritical CO 2 , and solid salt. The thermodynamic and thermophysical properties of brine-CO 2 mixtures (density, viscosity, and specific enthalpy of fluid phases; partitioning of mass components among the different phases) use the same correlations as an earlier fluid property module that does not distinguish between gaseous and liquid CO 2 -rich phases. We verified the fluid property module using two leakage scenarios, one that involves CO 2 migration up a blind fault and subsequent accumulation in a secondary “parasitic” reservoir at shallower depth, and another investigating leakage of CO 2 from a deep storage reservoir along a vertical fault zone. In task, “Development of a Rock Mechanical Module,” we developed a massively parallel reservoir simulator for modeling THM processes in porous media brine aquifers. We derived, from the fundamental equations describing deformation of porous elastic media, a momentum conservation equation relating mean stress, pressure, and temperature, and incorporated it alongside the mass and energy conservation equations from the TOUGH2 formulation, the starting point for the simulator. In addition, rock properties, namely permeability and porosity, are functions of effective stress and other variables that are obtained from the literature. We verified the simulator formulation and numerical implementation using analytical solutions and example problems from the literature. For the former, we matched a one-dimensional consolidation problem and a two-dimensional simulation of the Mandel-Cryer effect. For the latter, we obtained a good match of temperature and gas saturation profiles, and surface uplift, after injection of hot fluid into a model of a caldera structure. In task, “Incorporation of Geochemical Reactions of Selected Important Species,” we developed a novel mathematical model of THMC processes in porous and fractured saline aquifers, simulating geo-chemical reactions associated with CO 2 sequestration in saline aquifers. Two computational frameworks, sequentially coupled and fully coupled, were used to simulate the reactions and transport. We verified capabilities of the THMC model to treat complex THMC processes during CO 2 sequestration by analytical solutions and we constructed reactive transport models to analyze the THMC process quantitatively. Three of these are 1D reactive transport under chemical equilibrium, a batch reaction model with equilibrium chemical reactions, and a THMC model with CO 2 dissolution. In task “Study of Instability in CO 2 Dissolution-Diffusion-Convection Processes,” We reviewed literature related to the study of density driven convective flows and on the instability of CO 2 dissolution-diffusion-convection processes. We ran simulations that model the density-driven flow instability that would occur during CO 2 sequestration. CO 2 diffused through the top of the system and dissolved in the aqueous phase there, increasing its density. Density fingers formed along the top boundary, and coalesced into a few prominent ones, causing convective flow that forced the fluid to the system bottom. These simulations were in two and three dimensions. We ran additional simulations of convective mixing with density contrast caused by variable dissolved CO 2 concentration in saline water, modeled after laboratory experiments in which supercritical CO2 was circulated in the headspace above a brine saturated packed sand in a pressure vessel. As CO 2 dissolved into the upper part of the saturated sand, liquid phase density increases causing instability and setting off convective mixing. We obtained good agreement with the laboratory experiments, which were characterized by finger development and associated mixing of dissolved CO 2 into the system. We then varied a wide range of parameters and conceptual models in order to analyze the possibility of convective mixing under different conditions, such as various boundary conditions, and chemical reaction conditions. The CO 2 fingers from different simulations showed great differences as time progressed, caused by permeability heterogeneity. The early time diffusive phenomenon was captured by fine grid resolution, and the permeability heterogeneity affected the pattern of the CO 2 fingers. In addition, the fingers from three-dimensional simulations tended to be larger and flatter than the two-dimensional ones. In task “Implementation of Efficient Parallel Computing Technologies,” we made enhancements and modifications to our code in order to substantially increase the grid size that could be run. We installed and ran it on various platforms, including a multi-core PC and a cluster, and verified the numerical implementation and parallel code using an example problem from the literature. This problem, with a grid size of sixty million, utilized the cluster’s entire memory and all of its processors. In task “Implementation of General Fracture Conceptual Models,” we used the MINC approach, a generalization of the double-porosity concept, to model flow through porous and fractured media. In this approach, flow within the matrix is described by subdividing the matrix into nested volumes, with flow occurring between adjacent nested matrix volumes as well as between the fractures and the outer matrix volume. We generalized Hooke’s law to a thermo-multi- poroelastic medium, and derived from the fundamental equations describing deformation of porous and fractured elastic media a momentum conservation equation for thermo-multi- poroelastic media. This equation is a generalization to multi-poroelastic media of the one derived in Task 3.0 for single porosity media. We describe two simulations to provide model verification and application examples. The first, one-dimensional consolidation of a double-porosity medium, is compared to an analytical solution. The second is a match of published results from the literature, a simulation of CO 2 injection into hypothetical aquifer-caprock systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Protection against Chemical Warfare Agents and Biological Threats Using Metal–Organic Frameworks as Active Layers

The SARS-CoV-2 pandemic outbreak and the unfortunate misuse of toxic chemical warfare agents (CWAs) highlight the importance of developing functional materials to protect against these chemical and pathogen threats. Metal–organic frameworks (MOFs), which comprise a tunable class of crystalline porous materials built from inorganic nodes and organic linkers, have emerged as a class of heterogeneous catalysts capable of rapid detoxification of multiple classes of these harmful chemical or biological hazards. In particular, zirconium-based MOFs (Zr-MOFs) feature Lewis acidic nodes that serve as active sites for a wide range of catalytic reactions, including the hydrolysis of organophosphorus nerve agents within seconds in basic aqueous solutions. In addition, postsynthetic modification of Zr-MOFs enables the release of active species capable of reacting with and deactivating harmful pathogens. Despite this impressive performance, utilizing Zr- MOFs in powder form is not practical for application in masks or protective uniforms. To address this challenge, our team sought to develop MOF/fiber composite systems that could be adapted for use under realistic operating conditions to protect civilians, military personnel, and first responders from harmful pathogens and chemical warfare agents. Over the last several years, our group has designed and fabricated reactive and biocidal MOF/fiber composites that effectively capture and deactivate these toxic species. In this Account, we describe the evolution of these porous and reactive MOF/fiber composites and focus on key design challenges and considerations. First, we devised a scalable method for the integration of Zr-MOFs onto textile substrates using aqueous precursor solutions and without using pretreated textiles, highlighting the potential scalability of this method. Moving beyond standard textiles, we also developed a microbial synthesis strategy to prepare hierarchically porous MOF/bacterial cellulose nanofiber composite sponges that can both capture and detoxify nerve agents when exposed to contaminated gas flows. The mass loading of the MOF in the nanofibrous composite sponge is up to 90%, affording higher work capacities compared to those of textile-fiber-based composites with relatively lower MOF loadings. Next, we demonstrated that heterogeneous polymeric bases are suitable replacements for volatile liquid bases typically used in solution-phase reactions, and we showed that these composite systems are capable of effectively hydrolyzing nerve agents in the solid state by using only water that is present as humidity. Moreover, incorporating a reactive dye precursor into the composite affords a dual function sensing and detoxifying material that changes color from white to orange upon reaction with the byproduct following nerve agent hydrolysis, demonstrating the versatility of this platform for use in decontamination applications. We then created chlorine-loaded MOF/fiber composites that act as biocidal and reactive textiles that are capable of not only detoxifying sulfur-mustard-based chemical warfare agents and simulants but also deactivating both bacteria and the SARS-CoV-2 virus within minutes of exposure. Lastly, we synthesized a mixed-metal Ti/Zr-MOF coating on cotton fibers to afford a photoactive biocidal cloth that shows fast and broad-spectrum biocidal performance against viruses and Gram-positive and Gram-negative bacteria under visible light irradiation. Given the tunable, multifunctional nature of these MOF/fiber composites, we believe that this Account will offer new insights for the rational design and preparation of functional MOF/fiber composites and pave the way toward the development of next-generation reactive and protective textiles.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Efficiencies of Eu{sup 3+} ions and hydrogen atoms as donors in ZnO thin films

The donor efficiencies of Eu{sup 3+} ions and hydrogen atoms in ZnO crystalline films were investigated with reference to that of Ga{sup 3+} ions. It was found that Eu{sup 3+} ions acted as extrinsic donors in ZnO:Eu films, yielding a resistivity of 1.8 × 10{sup −3} Ω cm at a doping level of 1 at. %. This value is comparable to one for intrinsic donors in undoped ZnO films. The conductivity was maintained as the deposition temperature was increased to 200 °C, and this is evidence for the contribution of extrinsic donors. Deposition of Ga-doped and Eu-doped ZnO films in an H{sub 2}O gas flow produced oxyhydrogenated ZnO:(Ga, H) and ZnO:(Eu, H) films in which the Ga{sup 3+} and Eu{sup 3+} donors were deactivated by oxidization. Nevertheless, hydrogen donors contributed to electrical conduction yielding a resistivity of 1 × 10{sup −2} Ω cm. Postannealing in an H{sub 2} gas ambient alleviated the excessive oxidization of the films and thereby reactivated the donor action of Ga{sup 3+} and Eu{sup 3+} ions, causing the resistivity to recover to 10{sup −3} Ω cm for ZnO:(Ga, H) and 10{sup −2} Ω cm for ZnO:(Eu, H). In contrast, vacuum annealing of ZnO:(Ga, H) and ZnO:(Eu, H) films increased resistivity through removal of hydrogen donors while not affecting the oxidized condition of the samples.

42 ENGINEERING↗
DOE OSTI2017

Coupling a three-dimensional subsurface flow and transport model with a land surface model to simulate stream–aquifer–land interactions (CP v1.0)

A fully coupled three-dimensional surface and subsurface land model is developed and applied to a site along the Columbia River to simulate three-way interactions among river water, groundwater, and land surface processes. The model features the coupling of the Community Land Model version 4.5 (CLM4.5) and a massively parallel multiphysics reactive transport model (PFLOTRAN). The coupled model, named CP v1.0, is applied to a 400 m × 400 m study domain instrumented with groundwater monitoring wells along the Columbia River shoreline. CP v1.0 simulations are performed at three spatial resolutions (i.e., 2, 10, and 20 m) over a 5-year period to evaluate the impact of hydroclimatic conditions and spatial resolution on simulated variables. Results show that the coupled model is capable of simulating groundwater–river-water interactions driven by river stage variability along managed river reaches, which are of global significance as a result of over 30 000 dams constructed worldwide during the past half-century. Our numerical experiments suggest that the land-surface energy partitioning is strongly modulated by groundwater–river-water interactions through expanding the periodically inundated fraction of the riparian zone, and enhancing moisture availability in the vadose zone via capillary rise in response to the river stage change. Meanwhile, CLM4.5 fails to capture the key hydrologic process (i.e., groundwater–river-water exchange) at the site, and consequently simulates drastically different water and energy budgets. Furthermore, spatial resolution is found to significantly impact the accuracy of estimated the mass exchange rates at the boundaries of the aquifer, and it becomes critical when surface and subsurface become more tightly coupled with groundwater table within 6 to 7 meters below the surface. Inclusion of lateral subsurface flow influenced both the surface energy budget and subsurface transport processes as a result of river-water intrusion into the subsurface in response to an elevated river stage that increased soil moisture for evapotranspiration and suppressed available energy for sensible heat in the warm season. The coupled model developed in this study can be used for improving mechanistic understanding of ecosystem functioning and biogeochemical cycling along river corridors under historical and future hydroclimatic changes. The dataset presented in this study can also serve as a good benchmarking case for testing other integrated models.

58 GEOSCIENCES↗
DOE OSTI2021

A second order yield-temperature relation for accurate inference of burn-averaged quantities in multi-species plasmas

Measured yields and ion temperatures inferred from the fusion product energy spectra can be used as metrics for the performance of an ICF implosion. This can be to infer species separation, thermal decoupling, flows or other effects that can cause the inferred ion temperatures to deviate from the true underlying thermal temperature and the yield ratio to deviate from the expected value. Direct inference of the impact of these effects on observed temperatures and yields can be difficult to uncover due to underlying dependence on the shape and time evolution of the temperature and density proles of the fusing plasma. Due to differences in the temperature dependence of the reactivities, different fusion products are emitted from different regions and times within the implosion. In order to properly account for this, a second order analytic expression relating the apparent temperatures and yield ratios is developed. This expression can be coupled to models of yield and/or temperature altering effects to infer their burn-averaged impact on an implosion. The second order expression shows significant improvement over lower order expressions in synthetic data studies. Demonstrations of its applications to synthetic data coupled with models of ion thermal decoupling and radial flows are presented. In the case of thermal decoupling both first and second order expressions show reasonable levels of accuracy. To consistently infer the amplitude of radial flow with <10% error the second order equation is required.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Representative Elementary Volume (REV) in spatio-temporal domain: A method to find REV for dynamic pores

One of the potential risks associated with subsurface storage of CO 2 is the seepage of CO 2 through existing faults and fractures. There have been a number of studies devoted to this topic. Some of these studies show that geochemistry, especially mineralization, plays an important role in rendering the faults as conduits for CO 2 movement while others show that mineralization due to CO 2 injection can result in seep migration and flow diversion. Therefore, understanding the changes in reservoir properties due to pore alterations is important to ensure safe long term CO 2 storage in the subsurface. We study the changes in the Representative Elementary Volume (REV) of a rock due to reactive kinetics over a time, using a statistical approach and pore-scale CO 2 -rock interactiondata. The goal of this study is to obtain the REV of a rock property that accounts for pore-scale changes over time due to reactive kinetics, and we call this as spatiotemporal REV. Scale-up results suggest that the REV changes with time when CO 2 -rock interaction is considered. It is hypothesized that the alteration in pore structure introduces more heterogeneity in the rock, and because of this the magnitude of REV increases. Furthermore, it is possible that these noticeable changes in REV at pore-scale may have an impact when analyzed at the reservoir scale.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Numerical Simulation of Carbon Dioxide Injection in the Western Section of the Farnsworth Unit

Numerical simulation is an invaluable analytical tool for scientists and engineers in making predictions about of the fate of carbon dioxide injected into deep geologic formations for long-term storage. Current numerical simulators for assessing storage in deep saline formations have capabilities for modeling strongly coupled processes involving multifluid flow, heat transfer, chemistry, and rock mechanics in geologic media. Except for moderate pressure conditions, numerical simulators for deep saline formations only require the tracking of two immiscible phases and a limited number of phase components, beyond those comprising the geochemical reactive system. The requirements for numerically simulating the utilization and storage of carbon dioxide in partially depleted petroleum reservoirs are more numerous than those for deep saline formations. The minimum number of immiscible phases increases to three, the number of phase components may easily increase fourfold, and the coupled processes of heat transfer, geochemistry, and geomechanics remain. Public and scientific confidence in the ability of numerical simulators used for carbon dioxide sequestration in deep saline formations has advanced via a natural progression of the simulators being proven against benchmark problems, code comparisons, laboratory-scale experiments, pilot-scale injections, and commercial-scale injections. This paper describes a new numerical simulator for the scientific investigation of carbon dioxide utilization and storage in partially depleted petroleum reservoirs, with an emphasis on its unique features for scientific investigations; and documents the numerical simulation of the utilization of carbon dioxide for enhanced oil recovery in the western section of the Farnsworth Unit and represents an early stage in the progression of numerical simulators for carbon utilization and storage in depleted oil reservoirs.

carbon sequestration↗
DOE OSTI2016

A 2.5D Reactive Transport Model for Fracture Alteration Simulation

Understanding fracture alteration resulting from geochemical reactions is critical in predicting fluid migration in the subsurface and is relevant to multiple environmental challenges. Here in this paper, we present a novel 2.5D continuum reactive transport model that captures and predicts the spatial pattern of fracture aperture change and the development of an altered layer in the near-fracture region. The model considers permeability heterogeneity in the fracture plane and updates fracture apertures and flow fields based on local reactions. It tracks the reaction front of each mineral phase and calculates the thickness of the altered layer. Given this treatment, the model is able to account for the diffusion limitation on reaction rates associated with the altered layer. The model results are in good agreement with an experimental study in which a CO 2 -acidified brine was injected into a fracture in the Duperow Dolomite, causing dissolution of calcite and dolomite that result in the formation of a preferential flow channel and an altered layer. Finally, with an effective diffusion coefficient consistent with the experimentally observed porosity of the altered layer, the model captures the progressive decrease in the dissolution rate of the fast-reacting mineral in the altered layer.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1999

Approaches to New Endcaps for Improved Oxidation Resistance

Norbornenyl-end capped PMR polyimide resins are widely used as polymer matrix composite materials for aircraft engine applications, since they combine ease of processing with good oxidative stability up to 300 C. PMR resins are prepared by a two-step approach involving the initial formation of oligomeric pre-polymers capped at both ends by a latent reactive end cap. The end cap undergoes cross-linking during higher temperature processing, producing the desired low density, high specific strength materials, for PMR- 15. The end cap facilitates processing by controlling the molecular weight of the oligomer and allowing flow before it cross-links. However, after cross-linking, this very end cap accounts for much of the weight loss in the polymer on aging in air at elevated temperatures. Understanding this degradation provides clues for designing new end caps to slow down degradation, and prolong the lifetime of the material.

Meador, Mary Ann B.↗
DOE OSTI2024

Catalytic Resonance Theory: Forecasting the Flow of Programmable Catalytic Loops

Chemical transformations on catalyst surfaces occur through series and parallel reaction pathways. These complex networks and their behavior can be most simply evaluated through a three-species surface reaction loop (A* to B* to C* to A*) that is internal to the overall chemical reaction. Application of an oscillating dynamic catalyst to this reactive loop has been shown to exhibit one of three types of behavior: (1) a positive net flux of molecules about the loop in the clockwise direction, (2) a negative net flux of molecules about the loop in the counterclockwise direction, or (3) negligible flux of molecules about the loop at the limit cycle of reaction. Three-species surface loops were simulated with microkinetic modeling to assess the reaction loop behavior resulting from a catalytic surface oscillating between two or more catalyst surface energy states. Selected input parameters for the simulations spanned an 11-dimensional parameter space using 127 688 different parameter combinations. Their converged limit cycle solutions were analyzed for their loop turnover frequencies, the majority of which were found to be approximately zero. Classification and regression machine learning models were trained to predict the sign and magnitude of the loop turnover frequency and successfully performed above accessible baselines. Notably, the classification models exhibited a baseline weighted F1 score of 0.49, whereas trained models achieved weighted F1 scores of 0.94 and 0.96 when trained on the parameters used to define the simulations and derived rate constants, respectively. The trained models successfully predicted catalytic loop behavior, and interpretation of these models revealed all input parameters to be important for the prediction and performance of each model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Coupling a three-dimensional subsurface flow and transport model with a land surface model to simulate stream–aquifer–land interactions (CP v1.0)

Abstract. A fully coupled three-dimensional surface and subsurface land model is developed and applied to a site along the Columbia River to simulate three-way interactions among river water, groundwater, and land surface processes. The model features the coupling of the Community Land Model version 4.5 (CLM4.5) and a massively parallel multiphysics reactive transport model (PFLOTRAN). The coupled model, named CP v1.0, is applied to a 400 m × 400 m study domain instrumented with groundwater monitoring wells along the Columbia River shoreline. CP v1.0 simulations are performed at three spatial resolutions (i.e., 2, 10, and 20 m) over a 5-year period to evaluate the impact of hydroclimatic conditions and spatial resolution on simulated variables. Results show that the coupled model is capable of simulating groundwater–river-water interactions driven by river stage variability along managed river reaches, which are of global significance as a result of over 30 000 dams constructed worldwide during the past half-century. Our numerical experiments suggest that the land-surface energy partitioning is strongly modulated by groundwater–river-water interactions through expanding the periodically inundated fraction of the riparian zone, and enhancing moisture availability in the vadose zone via capillary rise in response to the river stage change. Meanwhile, CLM4.5 fails to capture the key hydrologic process (i.e., groundwater–river-water exchange) at the site, and consequently simulates drastically different water and energy budgets. Furthermore, spatial resolution is found to significantly impact the accuracy of estimated the mass exchange rates at the boundaries of the aquifer, and it becomes critical when surface and subsurface become more tightly coupled with groundwater table within 6 to 7 meters below the surface. Inclusion of lateral subsurface flow influenced both the surface energy budget and subsurface transport processes as a result of river-water intrusion into the subsurface in response to an elevated river stage that increased soil moisture for evapotranspiration and suppressed available energy for sensible heat in the warm season. The coupled model developed in this study can be used for improving mechanistic understanding of ecosystem functioning and biogeochemical cycling along river corridors under historical and future hydroclimatic changes. The dataset presented in this study can also serve as a good benchmarking case for testing other integrated models.

58 GEOSCIENCES↗
DOE OSTI2021

A second order yield-temperature relation for accurate inference of burn-averaged quantities in multi-species plasmas

Measured yields and ion temperatures inferred from the fusion product energy spectra can be used as metrics for the performance of an ICF implosion. This can be used to infer species separation, thermal decoupling, flows, or other effects that can cause the inferred ion temperatures to deviate from the true underlying thermal temperature and the yield ratio to deviate from the expected value. Direct inference of the impact of these effects on observed temperatures and yields can be difficult to uncover due to the underlying dependence on the shape and time evolution of the temperature and density profiles of the fusing plasma. Due to differences in the temperature dependence of the reactivities, different fusion products are emitted from different regions and times within the implosion. In order to properly account for this, a second-order analytical expression relating the apparent temperatures and yield ratios is developed. This expression can be coupled to models of yield and/or temperature altering effects to infer their burn-averaged impact on an implosion. The second-order expression shows significant improvement over lower-order expressions in synthetic data studies. Demonstrations of its applications to synthetic data coupled with models of ion thermal decoupling and radial flows are presented. In the case of thermal decoupling, both first and second-order expressions show reasonable levels of accuracy. To consistently infer the amplitude of radial flow with a <10% error, the second-order equation is required.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Nitrogen‐Doped Nanoporous Graphenic Carbon: An Efficient Conducting Support for O 2 Cathode

Abstract Herein, we synthesize a nanoporous N‐doped graphenic carbon with high surface area and a higher graphenization/graphitic degree by further developing the metallothermic reduction of gaseous CO 2 by adding N 2 into the gas flow. The N‐doped nanoporous carbon is composed of both a highly porous graphenic and non‐graphitic matrix and homogeneously dispersed ordered graphitic nanodomains, which constitute a very unique composite carbon structure. The resulting N‐doped graphenic carbon exhibits much more favorable reactivity as a carbon conducting support in the O 2 cathode of Li‐O 2 batteries, increasing the specific capacity of the GC electrode from 5300 to 9600 mAh g −1 . The N‐doped carbon also exhibits lower overpotentials during initial cycling for the charging process as well as an enhanced cycling performance compared to the undoped carbon. These results demonstrate that N‐doping has a strong correlation with the enhanced performance of O 2 cathode of Li‐O 2 batteries.

Xing, Zhenyu↗
DOE OSTI2022

TOUGH3-FLAC3D: a modeling approach for parallel computing of fluid flow and geomechanics

The recent development of the TOUGH3 code allows for a faster and more reliable fluid flow simulator. At the same time, new versions of FLAC3D are released periodically, allowing for new features and faster execution. In this paper, we present the first implementation of the coupling between TOUGH3 and FLAC3Dv6/7, maintaining parallel computing capabilities for the coupled fluid flow and geomechanical codes. We compare the newly developed version with analytical solutions and with the previous approach, and provide some performance analysis on different meshes and varying the number of running processors. Finally, we present two case studies related to fault reactivation during CO 2 sequestration and nuclear waste disposal. The use of parallel computing allows for meshes with a larger number of elements, and hence more detailed understanding of thermo-hydro-mechanical processes occurring at depth.

58 GEOSCIENCES↗
DOE OSTI2016

Effects of Asperity Distribution on Fluid Flow and Induced Seismicity During Deep Geothermal Exploitation

This work investigates the injection-induced seismic response of a heterogeneous fault plane, featuring low-permeability asperities embedded into a high-permeability damage zone. We simulate the pressure evolution with a hydrogeological simulator, accounting for the heterogeneous fault plane. Seismicity occurs then on the asperities, represented as unstable patches reactivating by means of the Mohr-Coulomb criterion. The hydrological and seismic modules are implicitly coupled to account for effects of asperity reactivation on the permeability. Results show that permeability changes may cause at a later time a change in seismicity propagation. Finally, we also investigated such effects by varying the density of asperities.

58 GEOSCIENCES↗