SEARCH · Engineering Papers
Results for “pathway analysis”
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.
Life Cycle Analysis of Electrofuels: Fischer–Tropsch Fuel Production from Hydrogen and Corn Ethanol Byproduct CO 2
Electrofuels from renewable H 2 and waste CO 2 streams are of increasing interest because of their CO 2 emissions reduction potentials compared to fossil counterparts. This study evaluated the well-to-wheel (WTW) greenhouse gas (GHG) emissions of Fischer-Tropsch (FT) fuels from various electrolytic H 2 pathways and CO 2 sources, using various process designs (i.e., with and without H 2 recycle) and system boundaries. Two systems with different boundaries were considered: a stand-alone plant (with CO 2 from any source) and an integrated plant with corn ethanol production (supplying CO 2 ). The FT fuel synthesis process was modeled using Aspen Plus, which showed that 45% of the carbon in CO 2 can be fixed in the FT fuel, with a fuel production energy efficiency of 58%. Using nuclear or solar/wind electricity, the stand-alone FT fuel production from various plant designs can reduce WTW GHG emissions by 90-108%, relative to petroleum fuels. When integrating the FT fuel production process with corn ethanol production, the WTW GHG emissions of FT fuels are 57-65% lower compared to petroleum counterparts. This study highlights the sensitivity of the carbon intensity of FT fuels to the system boundary selection (i.e., stand-alone vs integrated), which has different implications under various GHG emission credit frameworks.
Modeling Non-homologous End Joining
Non-homologous end joining (NHEJ) is the dominant DNA double strand break (DSB) repair pathway and involves several NHEJ proteins such as Ku, DNA-PKcs, XRCC4, Ligase IV and so on. Once DSBs are generated, Ku is first recruited to the DNA end, followed by other NHEJ proteins for DNA end processing and ligation. Because of the direct ligation of break ends without the need for a homologous template, NHEJ turns out to be an error-prone but efficient repair pathway. Some mechanisms have been proposed of how the efficiency of NHEJ repair is affected. The type of DNA damage is an important factor of NHEJ repair. For instance, the length of DNA fragment may determine the recruitment efficiency of NHEJ protein such as Ku [1], or the complexity of the DNA breaks [2] is accounted for the choice of NHEJ proteins and subpathway of NHEJ repair. On the other hand, the chromatin structure also plays a role of the accessibility of NHEJ protein to the DNA damage site. In this talk, some mathematical models of NHEJ, that consist of series of biochemical reactions complying with the laws of chemical reaction (e.g. mass action, etc.), will be introduced. By mathematical and numerical analysis and parameter estimation, the models are able to capture the qualitative biological features and show good agreement with experimental data. As conclusions, from the viewpoint of modeling, how the NHEJ proteins are recruited will be first discussed for connection between the classical sequential model [4] and recently proposed two-phase model [5]. Then how the NHEJ repair pathway is affected, by the length of DNA fragment [6], the complexity of DNA damage [7] and the chromatin structure [8], will be addressed
The role of AdhE mutations in Thermoanaerobacterium saccharolyticum
ABSTRACT Thermoanaerobacterium saccharolyticum is a thermophilic anaerobic bacterium that natively ferments a variety of hemicellulose substrates to organic acids and alcohols. It has recently been engineered to produce ethanol at high yield and titer; however, it uses a unique metabolic pathway for ethanol production that is poorly characterized. One of the distinctive aspects of this pathway is the presence of acetyl-CoA as an intermediate metabolite. In this organism, acetyl-CoA is converted to ethanol by a bifunctional AdhE enzyme. This enzyme has been a frequent target for mutations, and in many cases, the function of these mutations was unknown. Using a combination of genetic modifications, enzyme assays, and computational analysis, we have developed a better understanding of how mutations in AdhE affect ethanol production in the engineered homoethanologen strain. We identify a set of approximately interchangeable AdhE mutations (G544D, T597K, T597I, and T605I), whose function is to disrupt the activity of the alcohol dehydrogenase (ADH) domain of AdhE. This reduces NADH-linked ADH activity, which dramatically increases ethanol tolerance and changes the overall stoichiometry of acetaldehyde to ethanol conversion. Furthermore, our improved understanding of the function of these AdhE mutations calls into question a proposed feature of AdhE enzymes known as substrate channeling—direct transfer of acetaldehyde between the two domains of the AdhE enzyme. This improved the understanding of the role of AdhE mutations in T. saccharolyticum and provides deeper insights into the function of the unique ethanol production pathway in this organism. IMPORTANCE Many anaerobic bacteria maintain redox equilibrium by producing reduced organic compounds such as ethanol. The final two steps of ethanol production are mediated by a bifunctional enzyme, AdhE, and this enzyme is a frequent target of mutations in strains engineered for increased ethanol production. Paradoxically, these mutations increase ethanol production by eliminating the activity of one domain of the AdhE enzyme (the ADH domain). This provides additional support for a redox-imbalance theory of alcohol tolerance, which challenges the prevailing hypothesis that alcohol tolerance is associated with cell membrane effects.
Probing the Behavior of Composition‐Tunable Ultrathin PtNi Nanowires for CO Oxidation and Small‐Molecule Electrocatalytic Reactions
We have successfully synthesized ultrathin nanowires of pure Pt, Pt 99 Ni 1 , Pt 9 Ni 1 , and Pt 7 Ni 3 using a modified room-temperature soft-template method. Analysis of both methanol oxidation reaction (MOR) and ethanol oxidation reaction (EOR) results found that the Pt 7 Ni 3 samples yielded the best performance with specific activities of 0.36 and 0.34 mA/cm 2 respectively. Additionally, formic acid oxidation reaction (FAOR) tests noted that both Pt and PtNi nanowires oxidize small organic molecules (SOMs) via an indirect pathway. CO oxidation data suggests little measurable performance without any pre-reduction treatment; however, after annealing in H 2 , we detected significantly improved CO 2 formation for both Pt 9 Ni 1 and Pt 7 Ni 3 motifs. These observations highlight the importance of pre-treating these nanowires under a reducing atmosphere to enhance their performance for CO oxidation. To explain these findings, we collected extended x-ray adsorption fine structure (EXAFS) spectroscopy data, consistent with the presence of partial alloying with a tendency for Pt and Ni to segregate, thereby implying the formation of a Pt-rich shell coupled with a Ni-rich core. Here, we also observed that the degree of alloying within the nanowires increased after annealing in a reducing atmosphere, a finding deduced through analysis of the coordination numbers and calculations of Cowley's short range order parameters.
Multimodal quantification of degradation pathways during extreme fast charging of lithium-ion batteries
Enabling fast charging of Li-ion batteries will be a key step towards realizing the technology's full potential in electric vehicles. Currently, fast charging is limited by a variety of processes that reduce cell capacity upon extended cycling. Using a multimodal approach combining incremental capacity analysis (dQ/dV), high energy X-ray diffraction (HEXRD), and mass spectrometry titration (MST), we identify specific degradation mechanisms—including Li plating, dead Li x C 6 formation, Li 2 C 2 formation, solid carbonate solid-electrolyte interphase (SEI) deposition, and loss of positive electrode active material (LAM PE )—that occur during extended fast charge cycling. We find that Li plating is the major source of capacity loss in cells cycled at 6C, while non-carbonate SEI species deposition on the graphite anode is the main source of capacity loss when cycled at 4C. We also study local degradative phenomena by examining specific ~1–5 cm 2 regions of the cells using HEXRD and MST. Here, we find that plated Li is often collocated with dead Li x C 6 , Li 2 C 2 , and solid carbonate SEI species, and these additional species cumulatively account for ~20% of the capacity lost during 6C cycling. Finally, in a cell with an anomalously high amount of LAM PE (quantified via dQ/dV), we find that regions of cathode degradation were accompanied by non-carbonate SEI products on the adjacent region of the anode. We postulate that this phenomenon arises due to crosstalk between the electrodes, wherein soluble electrolyte oxidation products formed at the delithiated cathode migrate to the graphite anode and are ultimately deposited on the graphite surface. This work demonstrates the utility of combining multiple characterization techniques to reveal a more holistic understanding of degradative phenomena that occur across multiple length scales during fast charge.
Coordination-Induced N–H Bond Weakening in a Molybdenum Pyrrolidine Complex: Isotopic Labeling Provides Insight into the Pathway for H 2 Evolution
The synthesis and characterization of a cationic molybdenum pyrrolidine complex is described that exhibits significant coordination-induced N–H bond weakening. Here, the N–H bond dissociation free energy (BDFE) of the coordinated pyrrolidine in [( Ph Tpy)(PPh 2 Me) 2 Mo(NH(pyrr))][BArF 24 ] ([1-NH(pyrr)] + ; PhTpy = 4'-Ph-2,2',6',2''-terpyridine, NH(pyrr) = pyrrolidine, ArF 24 = [C 6 H 3 -3,5-(CF 3 ) 2 ] 4 ) was determined to be between 41–51 kcal mol -1 by thermochemical analysis and supported by a density functional theory (DFT) computed value of 48 kcal mol -1 . The complex [1-NH(pyrr)] + underwent proton-coupled electron transfer (PCET) to 2,4,6-tri-tert-butylphenoxyl radical, as well as spontaneous H 2 evolution upon gentle heating to furnish the corresponding molybdenum pyrrolidide complex, [( Ph Tpy)(PPh 2 Me) 2 Mo(N(pyrr))][BArF 24 ] ([1-N(pyrr)] + ). Thermolysis of the deuterated isotopolog, [1-ND(pyrr)] + still produced H 2 with concomitant incorporation of the isotopic label into the pyrrolidide ligand in the product [(1-N(pyrr-d n )] + (n = 0–2), consistent with an H 2 evolution pathway involving intramolecular H–H bond formation followed by an intermolecular product-forming PCET step. These observations provide context for understanding H 2 evolution in the nonclassical ammine complex [( Ph Tpy)(PPh 2 Me) 2 Mo(NH 3 )][BArF 24 ] ([1-NH 3 ] + ) and are supported by DFT-computed reaction thermochemistry. Overall, these studies offer rare insight into the H 2 formation pathway in nonclassical amine complexes with N–H BDFEs below the thermodynamic threshold for H 2 evolution and inform the development of well-defined, thermody-namically potent PCET reagents.
Systematic engineering for production of anti-aging sunscreen compound in Pseudomonas putida
Sunscreen has been used for thousands of years to protect skin from ultraviolet radiation. However, the use of modern commercial sunscreen containing oxybenzone, ZnO, and TiO 2 has raised concerns due to their negative effects on human health and the environment. In this study, we aim to establish an efficient microbial platform for production of shinorine, a UV light absorbing compound with anti-aging properties. First, we methodically selected an appropriate host for shinorine production by analyzing central carbon flux distribution data from prior studies alongside predictions from genome-scale metabolic models (GEMs). We enhanced shinorine productivity through CRISPRi-mediated downregulation and utilized shotgun proteomics to pinpoint potential competing pathways. Simultaneously, we improved the shinorine biosynthetic pathway by refining its design, optimizing promoter usage, and altering the strength of ribosome binding sites. Finally, we conducted amino acid feeding experiments under various conditions to identify the key limiting factors in shinorine production. The study combines meta-analysis of 13 C-metabolic flux analysis, GEMs, synthetic biology, CRISPRi-mediated gene downregulation, and omics analysis to improve shinorine production, demonstrating the potential of Pseudomonas putida KT2440 as platform for shinorine production.
Eliminating chemo-mechanical degradation of lithium solid-state battery cathodes during >4.5 V cycling using amorphous Nb2O5 coatings
Abstract Lithium solid-state batteries offer improved safety and energy density. However, the limited stability of solid electrolytes (SEs), as well as irreversible structural and chemical changes in the cathode active material, can result in inferior electrochemical performance, particularly during high-voltage cycling (>4.3 V vs Li/Li + ). Therefore, new materials and strategies are needed to stabilize the cathode/SE interface and preserve the cathode material structure during high-voltage cycling. Here, we introduce a thin (~5 nm) conformal coating of amorphous Nb 2 O 5 on single-crystal LiNi 0.5 Mn 0.3 Co 0.2 O 2 cathode particles using rotary-bed atomic layer deposition (ALD). Full cells with Li 4 Ti 5 O 12 anodes and Nb 2 O 5 -coated cathodes demonstrate a higher initial Coulombic efficiency of 91.6% ± 0.5% compared to 82.2% ± 0.3% for the uncoated samples, along with improved rate capability (10x higher accessible capacity at 2C rate) and remarkable capacity retention during extended cycling (99.4% after 500 cycles at 4.7 V vs Li/Li + ). These improvements are associated with reduced cell polarization and interfacial impedance for the coated samples. Post-cycling electron microscopy analysis reveals that the Nb 2 O 5 coating remains intact and prevents the formation of spinel and rock-salt phases, which eliminates intra-particle cracking of the single-crystal cathode material. These findings demonstrate a potential pathway towards stable and high-performance solid-state batteries during high-voltage operation.
Kinetics and mechanisms of polycondensation reactions between aryl halides and bisphenol A
Aryl chlorides (ArCl) or aryl fluorides (ArF) were used in polycondensation reactions to form poly(arylene ether sulfone)s (PAES). Interestingly, the kinetics of the ArF reaction fit a third-order rate law, which is attributed to the activation of the carbon–fluorine bond by two potassium cations (at least one bound to phenolate), which form a three-body complex. Here, the ArCl monomer follows a second-order rate law, where a two-body complex forms at the initial state of the aromatic nucleophilic substitution (S N Ar) pathway. These metal cation-activated complexes act as intermediates during the attack by the nucleophile. This finding was reproduced with either the potassium or the sodium counterion (introduced via potassium carbonate or sodium carbonate). Through a combination of experimental analysis of reaction kinetics and computational calculations with density functional theory (DFT) methods, the present work extends the fundamental understanding of polycondensation mechanisms for two aryl halides and highlights the importance of the CX–metal interaction(s) in the S N Ar reaction, which is translational to other ion-activated substitution reactions.
An unwanted guest in the electrochemical oxidation of high-voltage Li-ion battery electrolytes: the life of highly reactive protons
Lithium-ion batteries (LIBs) are central to the urgent societal need to decarbonize both transportation and energy storage on the grid. Unfortunately, despite their attractive energy/power density, as well as high coulombic and energy efficiencies, further improvement of this technology – especially their durability – is desperately needed. To support these efforts, our study focuses on fundamental understanding of the decomposition pathways for LIB electrolytes at the cathode–electrolyte interface (CEI), as the nature of these reactions directly controls the extent to which cell capacity and voltage decays in these systems. In this study, we employ electrochemical methods, coupled with product analysis using NMR spectroscopy and mass spectrometry, to determine the decomposition mechanisms in both model and technologically relevant electrolytes. Remarkably, we discovered the electrochemical formation of protons with high chemical activity, comparable to known superacids, at potentials relevant to practical Li-ion batteries. Their reactivity toward every individual component of the CEI provides a unified thermochemical origin for a myriad of side reactions that are commonly associated with the electrochemical reaction. In particular, electrochemically generated protons react with intact EC molecules to form CO 2 and other short and long chain ethers. They also undergo an acid–base reaction with LiPF 6 , to form the weaker acid HF, and with the cathode active material, leaching transition metals into the electrolyte. Collectively, the results of this study all point to the urgent need to either mitigate this proton formation or introduce benign harvesting additives via new electrolyte design strategies.
INDETERMINATE1 -mediated expression of FT family genes is required for proper timing of flowering in Brachypodium distachyon
The transition to flowering is a major developmental switch in plants. In many temperate grasses, perception of indicators of seasonal change, such as changing day-length and temperature, leads to expression of FLOWERING LOCUS T1 ( FT1 ) and FT-Like ( FTL ) genes that are essential for promoting the transition to flowering. However, little is known about the upstream regulators of FT1 and FTL genes in temperate grasses. Here, we characterize the monocot-specific gene INDETERMINATE1 ( BdID1 ) in Brachypodium distachyon and demonstrate that BdID1 is a regulator of FT family genes. Mutations in ID1 impact the ability of the short-day (SD) vernalization, cold vernalization, and long-day (LD) photoperiod pathways to induce certain FTL genes. BdID1 is required for upregulation of FTL9 ( FT-LIKE9 ) expression by the SD vernalization pathway, and overexpression of FTL9 in an id1 background can partially restore the delayed flowering phenotype of id1 . We show that BdID1 binds in vitro to the promoter region of FTL genes suggesting that ID1 directly activates FTL expression. Transcriptome analysis shows that BdID1 is required for FT1 , FT2 , FTL12 , and FTL13 expression under inductive LD photoperiods, indicating that BdID1 is a regulator of the FT gene family. Moreover, overexpression of FT1 in the id1 background results in rapid flowering similar to overexpressing FT1 in the wild type, demonstrating that BdID1 is upstream of FT family genes. Interestingly, ID1 negatively regulates a previously uncharacterized FTL gene, FTL4, and we show that FTL4 is a repressor of flowering. Thus, BdID1 is critical for proper timing of flowering in temperate grasses.
Advancing the Predictability of Water Cycle Phenomena via the Application of AI to Model Ensemble Simulations and Observations
AI has the power to identify new pathways to extended predictability of the water cycle via its application to model ensembles together with accompanying observations. We discuss both scientific and technological aspects of this challenge, and address the parts of the MODEX approach involving “Model simulations, evaluation, analysis, and benchmarking” and “Identification of key knowledge gaps”. The widespread incorporation of AI into model analysis will significantly advance Earth system predictability and our predictive capabilities.
Model form and sensitivity analysis of CALPHAD-based nucleation models in b-stabilized Ti alloys
Accurate prediction of α-phase nucleation and growth in β-stabilized titanium alloys is crucial for designing heat treatments to optimize mechanical properties in additively manufactured lightweight components. Ideally, predictions of nucleation and growth would incorporate both top-down observations of past experimental heat treatments and bottom-up modeling of phase transformations; however, the appropriate method of combining these information sources is not self-evident. Combining top-down and bottom-up information requires a unified form of model that can connect between spatiotemporal scales, as well as sets of fitting parameters that can be identified by each data source. The selection of which parameters to fit to which data source can be made based on expert opinion, or by performing a sensitivity analysis. In solid-solid nucleation, direct observation of the nucleation and growth process is challenging. Most data on the heat treatment-controlled phase transformations are not in-situ. To predict the process and outcome of the nucleation, growth and coarsening of precipitates, theoretical models of the nucleation pathway are used to bridge the gap. Many sources of uncertainty affect the modeling of this nucleation process. It can be influenced by small variations in the thermomechanical processing history, chemical composition, and initial microstructure. If molecular dynamics (MD) simulations are used to determine thermodynamic quantities and inform CALPHAD modeling, additional uncertainty can be introduced and accounted for using Bayesian methods. Top-down uncertainties require additional steps to quantify. The influence of nucleation model form on the sensitivity of predictions to input parameters and physical conditions is the focus of this study. Classical nucleation theory (CNT) allows modeling to formulate the nucleation as homogeneous or, more commonly, heterogeneous. Non-classical nucleation models are also increasingly explored as a means of reconciling top-down and bottom-up data. In this study, the sensitivity of the intragranular nucleation of α in a β-annealed, slow-cooled aging (BASCA) heat treatment of β-stabilized Ti5553 alloy is explored using CNT and both heterogeneous and homogeneous assumptions. The Kampmann-Wagner Numerical model of precipitate nucleation and growth is employed. Using open-source tools (pyCalphad and thermodynamic modeling of TiMo as a surrogate system, a sensitivity analysis is performed to measure variations in key parameters, including chemical driving force, interfacial energy, and diffusivity, as they relate to predictions of precipitate number density. The inclusion of top-down and bottom-up data in selection of nucleation model form is discussed.
A system and methodology for measuring volatile organic compounds produced by hydroponic lettuce in a controlled environment
A system and methodology were developed for the nondestructive qualitative and quantitative analysis of volatile emissions from hydroponically grown 'Waldmann's Green' leaf lettuce (Lactuca sativa L.). Photosynthetic photon flux (PPF), photoperiod, and temperature were automatically controlled and monitored in a growth chamber modified for the collection of plant volatiles. The lipoxygenase pathway products (Z)-3-hexenal, (Z)-3-hexenol, and (Z)-3-hexenyl acetate were emitted by lettuce plants after the transition from the light period to the dark period. The volatile collection system developed in this study enabled measurements of volatiles emitted by intact plants, from planting to harvest, under controlled environmental conditions.
Overview and Initial Results of SAND-E: Semi-Autonomous Navigation for Detrital Environments
Unmanned aerial systems (UAS) and automated terrain analysis for science and navigation are new technologies for planetary exploration. The Mars Helicopter will fly with the Mars2020 rover, the Dragonfly quadcopter will explore Titan, and Soil Properties and Object Classification (SPOC) software will be used for path planning and navigation on the Mars2020 rover. Using an Argo J5 rover instrumented with stereo cameras and Autonomous Soil Assessment System (ASAS) software, and an off the shelf quadcopter, SAND-E tested the use of automated terrain analysis and UAS data for science operations in a Mars-analog environment in Iceland during July of 2019. Scientifically, we sought to determine changes in the physical and chemical properties of sediments along a glacial-fluvial-aeolian transport pathway. Operationally, we tested rover mission-like scenarios that included UAS images and classified terrain images. Here, we present the initial results for both the operations and science elements of the study. Site Selection: A goal of SAND-E is examine sorting and alteration of sediments in fluvial and aeolian environments in both mineral-dominated and glass-dominated basaltic settings. During the first year of the project we focused on a mineral-dominated environment. Selection of the location was based on prior publications that indicated our selected region had a greater abundance of crystalline sediments than other areas fluvial-aeolian settings in Iceland. Other criteria included the presence of both fluvial and aeolian landforms along a transport pathway such that the sediments in transport could be linked to their source rocks. We chose the Skjaldbreidauhraun glacial outwash plain, which sits at the base of Thórisjökull glacier. The site is 30 km north of Thingvellir National Park and ~2 hours from Reykjavik. The outwash plain is fed by two small catchments that drain from the base of the glacier and cut through hyaloclastite and shield volcano bedrock. The drainage progresses from steep alluvial fans near the glacier into a low-sloping fluvial braidplain that becomes confined by the Skjaldbreidur shield volcano and creates a shallow canyon cut into lava bedrock. The fluvial system was a typical braided alluvial environment composed pebble- and cobble-bedded longitudinal bars and sandy channel beds. The river remained active and fluctuated in response to diurnal runoff cycles near the glacier before disappearing into the sandy substrate downstream. The high concentration of suspended sediment in the river was evident by the cloudy water and the silt and clay-sized sediments that draped the channel beds after abandonment and created playas in the lowest sloping areas of the catchment. The entire fluvial system was affected by the winds generated by frontal systems and katabatic flows descending the glacier. This resulted in the formation of aeolian lag deposits and a wind-deflation plain where the fluvial system was not active. Wind ripples and drifts formed in abandoned fluvial channels from aeolian reworking of the sand-sized fluvial sediments. The silt- and clay-sized sediments found in fluvial channels, bar tops, and playas generated dust plumes during high wind events. Our operation sought to capture the variability in this system by sampling from the range of fluvial and aeolian features 6.3 km (proximal), 11.3 km (medial), and 14.4 km (distal) along the river from its origin at the base of glacier.
Engineering isoprenoids production in metabolically versatile microbial host Pseudomonas putida
Abstract With the increasing need for microbial bioproduction to replace petrochemicals, it is critical to develop a new industrial microbial workhorse that improves the conversion of lignocellulosic carbon to biofuels and bioproducts in an economically feasible manner. Pseudomonas putida KT2440 is a promising microbial host due to its capability to grow on a broad range of carbon sources and its high tolerance to xenobiotics. In this study, we engineered P. putida KT2440 to produce isoprenoids, a vast category of compounds that provide routes to many petrochemical replacements. A heterologous mevalonate (MVA) pathway was engineered to produce potential biofuels isoprenol (C 5 ) and epi-isozizaene (C 15 ) for the first time in P. putida . We compared the difference between three different isoprenoid pathways in P. putida on isoprenol production and achieved 104 mg/L of isoprenol production in a batch flask experiment through optimization of the strain. As P. putida can natively consume isoprenol, we investigated how to prevent this self-consumption. We discovered that supplementing l -glutamate in the medium can effectively prevent isoprenol consumption in P. putida and metabolomics analysis showed an insufficient energy availability and an imbalanced redox status during isoprenol degradation. We also showed that the engineered P. putida strain can produce isoprenol using aromatic substrates such as p -coumarate as the sole carbon source, and this result demonstrates that P. putida is a valuable microbial chassis for isoprenoids to achieve sustainable biofuel production from lignocellulosic biomass. Graphical Abstract
Phenomenological Model of Nonlinear Dynamics and Deterministic Chaotic Gas Migration in Bentonite: Experimental Evidence and Diagnostic Parameters
Understanding gas migration in compacted clay materials, e.g., bentonite and claystone, is important for the design and performance assessment of an engineered barrier system of a radioactive waste repository system, as well as many practical applications. Existing field and laboratory data on gas migration processes in low-permeability clay materials demonstrate the complexity of flow and transport processes, including various types of instabilities, caused by nonlinear dynamics of coupled processes of liquid–gas exchange, dilation, fracturing, fracture healing, etc., which cannot be described by classical models of fluid dynamics in porous media. We here show that the complexity of gas migration processes can be explained using a phenomenological concept of nonlinear dynamics and deterministic chaos theory. To do so, we analyzed gas pressure and gas influx (i.e., input) and outflux (i.e., output), recorded during the gas injection experiment in the compact Mx80-D bentonite sample, and calculated a set of the diagnostic parameters of nonlinear dynamics and chaos, such a global embedding dimension, a correlation dimension, an information dimension, and a spectrum of Lyapunov exponents, as well as plotted 2D and 3D pseudo-phase-space strange attractors, based on the univariate influx and outflux time series data. These results indicate the presence of phenomena of low-dimensional deterministic chaotic behavior of gas migration in bentonite. In particular, during the onset of gas influx in the bentonite core, before the breakthrough, the development of gas flow pathways is characterized by the process of chaotic gas diffusion. After the breakthrough, with inlet-to-outlet movement of gas, the prevailing process is chaotic advection. During the final phase of the experiment, with no influx to the sample, the relaxation pattern of gas outflux is resumed back to a process of chaotic diffusion. The types of data analysis and a proposed phenomenological model can be used to establish the basic principles of experimental data-gathering, modeling predictions, and a research design.