Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “experimental data”

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,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A Semi‐Automated, High‐Throughput Approach for the Synthesis and Identification of Highly Photo‐Cytotoxic Iridium Complexes

Abstract The discovery of new compounds with pharmacological properties is usually a lengthy, laborious and expensive process. Thus, there is increasing interest in developing workflows that allow for the rapid synthesis and evaluation of libraries of compounds with the aim of identifying leads for further drug development. Herein, we apply combinatorial synthesis to build a library of 90 iridium(III) complexes (81 of which are new) over two synthesise‐and‐test cycles, with the aim of identifying potential agents for photodynamic therapy. We demonstrate the power of this approach by identifying highly active complexes that are well‐tolerated in the dark but display very low nM phototoxicity against cancer cells. To build a detailed structure–activity relationship for this class of compounds we have used density functional theory (DFT) calculations to determine some key electronic parameters and study correlations with the experimental data. Finally, we present an optimised semi‐automated synthesise‐and‐test protocol to obtain multiplex data within 72 hours.

Kench, Timothy↗
DOE OSTI2024

Propane Dehydrogenation on Pt x Zn y Active Sites in Silicalite‐1

Abstract The improvement of Pt‐based catalysts for propane dehydrogenation (PDH) has progressed by recent investigations that have identified Zn as a promising promoter for Pt subnanometer catalysts. It is desirable to gain insights into the structure, stability, and activity of such active sites and the factors that influence them, such as Zn : Pt ratio, Pt coordination and nuclearity. Here, we employ density functional theory and microkinetic simulations to investigate the stability of Pt x Zn y ( x =1–3, y=0–3) active sites grafted on silanols of Silicalite‐1 and the PDH activity of Pt. We find that the coordination of a Pt atom to a nest of grafted Zn(II) atoms increases the stability of the Pt 1 Zn y sites, whose activity is similar for y=0–2 and drops dramatically for y>2. We further demonstrate, via linear scaling relations and microkinetic simulations, that the turnover frequency obeys a volcano law as a function of propylene binding strength. The Pt 2 Zn 1 and Pt 3 Zn 1 sites are stable and exhibit activity similar to Pt 1 Zn 2 , but only Pt 1 Zn 2 manifests reaction kinetics consistent with experimental data, strongly suggesting the active site composition in the synthesized catalyst samples. The methodology presented here suggests a general strategy for deducing active site information such as composition through simple kinetic experiments.

Liu, Yilang↗
DOE OSTI2023

Automated Generation of Microkinetics for Heterogeneously Catalyzed Reactions Considering Correlated Uncertainties**

Abstract The study presents an ab‐initio based framework for the automated construction of microkinetic mechanisms considering correlated uncertainties in all energetic parameters and estimation routines. 2000 unique microkinetic models were generated within the uncertainty space of the BEEF‐vdW functional for the oxidation reactions of representative exhaust gas emissions from stoichiometric combustion engines over Pt(111) and compared to experiments through multiscale modeling. The ensemble of simulations stresses the importance of considering uncertainties. Within this set of first‐principles‐based models, it is possible to identify a microkinetic mechanism that agrees with experimental data. This mechanism can be traced back to a single exchange‐correlation functional, and it suggests that Pt(111) could be the active site for the oxidation of light hydrocarbons. The study provides a universal framework for the automated construction of reaction mechanisms with correlated uncertainty quantification, enabling a DFT‐constrained microkinetic model optimization for other heterogeneously catalyzed systems.

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

A Semi‐Automated, High‐Throughput Approach for the Synthesis and Identification of Highly Photo‐Cytotoxic Iridium Complexes

Abstract The discovery of new compounds with pharmacological properties is usually a lengthy, laborious and expensive process. Thus, there is increasing interest in developing workflows that allow for the rapid synthesis and evaluation of libraries of compounds with the aim of identifying leads for further drug development. Herein, we apply combinatorial synthesis to build a library of 90 iridium(III) complexes (81 of which are new) over two synthesise‐and‐test cycles, with the aim of identifying potential agents for photodynamic therapy. We demonstrate the power of this approach by identifying highly active complexes that are well‐tolerated in the dark but display very low nM phototoxicity against cancer cells. To build a detailed structure–activity relationship for this class of compounds we have used density functional theory (DFT) calculations to determine some key electronic parameters and study correlations with the experimental data. Finally, we present an optimised semi‐automated synthesise‐and‐test protocol to obtain multiplex data within 72 hours.

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

Propane Dehydrogenation on Pt x Zn y Active Sites in Silicalite‐1

Abstract The improvement of Pt‐based catalysts for propane dehydrogenation (PDH) has progressed by recent investigations that have identified Zn as a promising promoter for Pt subnanometer catalysts. It is desirable to gain insights into the structure, stability, and activity of such active sites and the factors that influence them, such as Zn : Pt ratio, Pt coordination and nuclearity. Here, we employ density functional theory and microkinetic simulations to investigate the stability of Pt x Zn y ( x =1–3, y=0–3) active sites grafted on silanols of Silicalite‐1 and the PDH activity of Pt. We find that the coordination of a Pt atom to a nest of grafted Zn(II) atoms increases the stability of the Pt 1 Zn y sites, whose activity is similar for y=0–2 and drops dramatically for y>2. We further demonstrate, via linear scaling relations and microkinetic simulations, that the turnover frequency obeys a volcano law as a function of propylene binding strength. The Pt 2 Zn 1 and Pt 3 Zn 1 sites are stable and exhibit activity similar to Pt 1 Zn 2 , but only Pt 1 Zn 2 manifests reaction kinetics consistent with experimental data, strongly suggesting the active site composition in the synthesized catalyst samples. The methodology presented here suggests a general strategy for deducing active site information such as composition through simple kinetic experiments.

Liu, Yilang↗
DOE OSTI2022

Resolving Graphite‐Electrolyte Interphase of Lithium‐Ion Batteries using Air‐Tight Ambient Mass Spectrometry

Abstract The formation of a stable solid electrolyte interphase (SEI) at graphite electrode is necessary for lithium‐ion batteries (LIB) to prevent excessive degradation of electrolytes. Until now, the exact organic composition of SEI has not been decisively determined with lithium ethylene mono‐carbonate (LEMC) being one of the debated components. An air‐tight ambient mass spectrometry (MS) platform is developed to directly characterize the air‐sensitive organic components present in SEI. Ionization occurs via a solvent‐free field‐induced process at atmospheric pressure. Our experimental data, based on tandem MS analysis, revealed LEMC to be the major organic component in SEI. This result resolves a long‐standing mystery surrounding the composition of native SEI in LIB, where no trace of the other debated compound, lithium ethylene dicarbonate, was detected. The energetics of fragmentation of LEMC and the corresponding stabilities of various ionic species detected during the air‐tight ambient MS analysis were calculated using density functional theory (DFT). Our findings confirmed favorable protonation energies and suggested the molecular structures for many of the detected species. Solvent (i. e., ethylene carbonate)‐induced stabilization effects observed during the ionization process were also confirmed by DFT. The reported study encourages the routine use of ambient mass spectrometry to characterize both anodic and cathodic electrodes of all battery types as well as to study other air‐sensitive materials in their native state.

25 ENERGY STORAGE↗
DOE OSTI2020

Kinetic modeling of microbial growth, enzyme activity, and gene deletions: An integrated model of β-glucosidase function in Cellvibrio japonicus

Understanding the complex growth and metabolic dynamics in microorganisms requires advanced kinetic models containing both metabolic reactions and enzymatic regulation to predict phenotypic behaviors under different conditions and perturbations. Most current kinetic models lack gene expression dynamics and are separately calibrated to distinct media, which consequently makes them unable to account for genetic perturbations or multiple substrates. This challenge limits our ability to gain a comprehensive understanding of microbial processes towards advanced metabolic optimizations that are desired for many biotechnology applications. Here, we present an integrated computational and experimental approach for the development and optimization of mechanistic kinetic models for microbial growth and metabolic and enzymatic dynamics. Our approach integrates growth dynamics, gene expression, protein secretion, and gene-deletion phenotypes. We applied this methodology to build a dynamic model of the growth kinetics in batch culture of the bacterium Cellvibrio japonicus grown using either cellobiose or glucose media. The model parameters were inferred from an experimental data set using an evolutionary computation method. The resulting model was able to explain the growth dynamics of C. japonicus using either cellobiose or glucose media and was also able to accurately predict the metabolite concentrations in the wild-type strain as well as in β-glucosidase gene deletion mutant strains. We validated the model by correctly predicting the non-diauxic growth and metabolite consumptions of the wild-type strain in a mixed medium containing both cellobiose and glucose, made further predictions of mutant strains growth phenotypes when using cellobiose and glucose media, and demonstrated the utility of the model for designing industrially-useful strains. Importantly, the model is able to explain the role of the different β-glucosidases and their behavior under genetic perturbations. Here, this integrated approach can be extended to other metabolic pathways to produce mechanistic models for the comprehensive understanding of enzymatic functions in multiple substrates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Protective effects of non–catalytic proteins on endoglucanase activity at air and lignin interfaces

The manner in which added non-catalytic proteins during enzymatic hydrolysis of lignocellulosic substrates enhances hydrolysis mechanisms is not completely understood. Prior research has indicated that a reduction in the non-specific adsorption of enzymes on lignin, and deactivation of enzymes exposed to air–liquid interface provide rationale. This work investigated root causes including effects of the air–liquid interface on non-catalytic proteins, and effects of lignin on endoglucanase. Three different experimental designs and three variables (air–liquid interfacial area, the types of lignin (acid or enzymatic lignin), and the presence of non-enzymatic protein (bovine serum albumin [BSA] or soy proteins ) were used. The results showed that acid isolated lignin adsorbed almost all endoglucanase activity initially present in supernatant, independent of air interface conditions (25 or 250 ml flasks) with the presence of BSA preventing this effect. Endoglucanase lost 30%–50% of its activity due to an air–liquid interface in the presence of lignin while addition of non-enzymatic protein helped to preserve this enzyme's activity. Langmuir and Freundlich models applied to experimental data indicated that the adsorption increases with increasing temperature for both endoglucanase and BSA. Adsorption of the enzyme and protein were endothermic with an increase in entropy. Finally, these results, combined, show that hydrophobicity plays a strong role in the adsorption of both endoglucanase and BSA on lignin.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

CheKiPEUQ Intro 1: Bayesian Parameter Estimation Considering Uncertainty or Error from both Experiments and Theory**

A common goal is extraction of physico-chemical parameter values such as pre-exponentials and activation energies from experiment. Ever increasing knowledge from experiments and computations is enabling semi-quantitative prior predictions of such values. When prior knowledge of physically realistic ranges is available, a method named Bayesian parameter estimation (BPE) enables more physically realistic parameter estimation relative to unsophisticated fitting by seeking the most probable value when considering together the uncertainties from prior knowledge, experimental data, and approximations in the model. An impediment to widespread use of BPE is a lack of understanding, training, and user-friendly software. Along with this invited publication, a general software package for BPE is being released that is user-friendly and that does not require understanding of the math behind the methodology. Overall, two previously unpublished catalysis science examples are provided along with considerations and guidelines for successful application of BPE. Following this work, BPE can become more widespread to enable extraction of physically meaningful parameter values.

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

Fast Correction of Errors in the DFT‐Calculated Energies of Gaseous Nitrogen‐Containing Species

Abstract Modeling adsorption phenomena on surfaces by DFT calculations often involves substantial errors, resulting in inaccurate predictions of catalytic activities. Such errors partly stem from the inaccurate description of the energetics of free molecules. Herein, we use a semiempirical group‐additivity method to correct the DFT‐calculated heats of formation of 106 carbon‐ and nitrogen‐containing gaseous compounds belonging to 15 different chemical families. PBE, PW91, RPBE and BEEF‐vdW initially yield mean absolute errors (MAEs) with respect to experiments in the range of 0.32–0.75 eV. After correcting the systematic errors, the overall MAEs decrease to ∼0.05 eV. Additionally, upon applying the corrections to three types of reaction enthalpies, the resulting MAEs are below 0.10 eV. These functional‐group corrections can be used in (electro)catalysis to correct the gas‐phase references necessary to evaluate equilibrium potentials and adsorption energies, predict error cancellation, and assess conflicting experimental data.

Urrego‐Ortiz, Ricardo↗
DOE OSTI2021

Solvothermal Synthesis of [Cr 7 S 8 (en) 8 Cl 2 ]Cl 3 ⋅ 2H 2 O with Magnetically Frustrated [Cr 7 S 8 ] 5+ Double‐Cubes**

Abstract A novel transition metal chalcohalide [Cr 7 S 8 (en) 8 Cl 2 ]Cl 3 ⋅ 2H 2 O, with [Cr 7 S 8 ] 5+ dicubane cationic clusters, has been synthesized by a low temperature solvothermal method, using dimethyl sulfoxide (DMSO) and ethylenediamine ( en ) solvents. Ethylenediamine ligand exhibits bi‐ and monodentate coordination modes; in the latter case ethylenediamine coordinates to Cr atoms of adjacent clusters, giving rise to a 2D polymeric structure. Although magnetic susceptibility shows no magnetic ordering down to 1.8 K, a highly negative Weiss constant, θ =−224(2) K, obtained from Curie‐Weiss fit of inverse susceptibility, suggests strong antiferromagnetic (AFM) interactions between S =3/2 Cr(III) centers. Due to the complexity of the system with (2 S +1) 7 =16384 microstates from seven Cr 3+ centers, a simplified model with only two exchange constants was used for simulations. Density‐functional theory (DFT) calculations yielded the two exchange constants to be J 1 =−21.4 cm −1 and J 2 =−30.2 cm −1 , confirming competing AFM coupling between the shared Cr 3+ center and the peripheral Cr 3+ ions of the dicubane cluster. The best simulation of the experimental data was obtained with J 1 =−20.0 cm −1 and J 2 =−21.0 cm −1 , in agreement with the slightly stronger AFM exchange within the triangles of the peripheral Cr 3+ ions as compared to the AFM exchange between the central and peripheral Cr 3+ ions. This compound is proposed as a synthon towards magnetically frustrated systems assembled by linking dicubane transition metal‐chalcogenide clusters into polymeric networks.

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

Binding, Release and Functionalization of Intact Pnictogen Tetrahedra Coordinated to Dicopper Complexes

The bridging MeCN ligand in the dicopper(I) complexes [(DPFN)Cu 2 (μ,η 1 : η 1 -MeCN)][X] 2 (X=weakly coordinating anion, NTf 2 (1 a), FAl[OC 6 F 10 (C 6 F 5 )] 3 (1 b), Al[OC(CF 3 ) 3 ] 4 (1 c)) was replaced by white phosphorus (P 4 ) or yellow arsenic (As 4 ) to yield [(DPFN)Cu 2 (μ,η 2 : η 2 -E 4 )][X] 2 (E=P (2 a-c), As (3 a-c)). The molecular structures in the solid state reveal novel coordination modes for E 4 tetrahedra bonded to coinage metal ions. Experimental data and quantum chemical computations provide information concerning perturbations to the bonding in coordinated E 4 tetrahedra. Reactions with N-heterocyclic carbenes (NHCs) led to replacement of the E 4 tetrahedra with release of P 4 or As 4 and formation of [(DPFN)Cu 2 (μ,η 1 : η 1 - Me NHC)][X] 2 (4 a,b) or to an opening of one E-E bond leading to an unusual E 4 butterfly structural motif in [(DPFN)Cu 2 (μ,η 1 : η 1 -E 4 Dipp NHC)][X] 2 (E=P (5 a,b), E=As (6)). With a cyclic alkyl amino carbene ( Et CAAC), cleavage of two As-As bonds was observed to give two isomers of [(DPFN)Cu 2 (μ,η 2 : η 2 -As 4 Et CAAC)][X] 2 (7 a,b) with an unusual As 4 -triangle+1 unit.

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

Gamma Radiation-Induced Degradation of Acetohydroxamic Acid (AHA) in Aqueous Nitrate and Nitric Acid Solutions Evaluated by Multiscale Modelling

Acetohydroxamic acid (AHA) has been proposed for inclusion in advanced, single-cycle, used nuclear fuel reprocessing solvent systems for the reduction and complexation of plutonium and neptunium ions. For this application, a detailed description of the degradation of AHA in dilute aqueous nitric acid is required. To this end, we present a comprehensive, multiscale computer model for the coupled radiolytic and hydrolytic degradation of AHA in aqueous sodium nitrate and nitric acid solutions. Rate coefficients for the reactions of AHA and HA with the oxidizing nitrate radical were measured for the first time using pulsed electron radiolysis and used as inputs in the kinetic model. The computer model results are validated by comparison to experimental data from steady state gamma irradiations, for which the agreement is excellent. Here, the presented model accurately predicts the yields of the major degradation products of AHA: acetic acid, hydroxylamine, nitrous oxide, and molecular hydrogen.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Molten‐Salt‐Mediated Chemical Looping Oxidative Dehydrogenation of Ethane with In‐Situ Carbon Capture and Utilization

Abstract The molten‐salt‐mediated oxidative dehydrogenation (MM‐ODH) of ethane (C 2 H 6 ) via a chemical looping scheme represents an effective carbon capture and utilization (CCU) method for the valorization of ethane‐rich shale gas and concurrent mitigation of carbon dioxide (CO 2 ) emissions. Here, stepwise experimentation with Li 2 CO 3 ‐Na 2 CO 3 ‐K 2 CO 3 (LNK) ternary salts (i) assessed how each component of the LNK mixture impacted ethane MM‐ODH performance and (ii) explored physicochemical and thermodynamic mechanisms behind melt‐induced changes to ethylene (C 2 H 4 ) and carbon monoxide (CO) yields. Of fifteen screened LNK compositions, nine exhibited ethylene yields greater than 50 % at 800 °C while maintaining C 2 H 4 selectivities of 85 % or higher. LNK salts rich in Li 2 CO 3 content yielded more ethylene and CO on average than their counterparts, and net CO 2 capture per cycle reached a maximum of ~75 %. Extended MM‐ODH cycling also demonstrated long‐term stability of a high‐performing LNK medium. Density functional theory (DFT) calculations andab initiomolecular dynamics (AIMD) simulations suggested that the molten salt does not directly activate C 2 H 6 . Meanwhile, an empirical model informed by experimental data and reaction thermodynamics adequately predicted overall MM‐ODH performance from LNK composition and provided insights into the system′s primary drivers.

Chemistry↗
DOE OSTI2022

Nitrification, denitrification, and competition for soil N : Evaluation of two Earth System Models against observations

Abstract Earth System Models (ESMs) have implemented nitrogen (N) cycles to account for N limitation on terrestrial carbon uptake. However, representing inputs, losses, and recycling of N in ESMs is challenging. Here, we use global rates and ratios of key soil N fluxes, including nitrification, denitrification, mineralization, leaching, immobilization, and plant uptake (both NH 4 + and NO 3 − ), from the literature to evaluate the N cycles in the land model components of two ESMs. The two land models evaluated here, E3SM Land Model version 1 (ELMv1)‐ECA and CLM5.0, originated from a common model but have diverged in their representation of plant–microbe competition for soil N. The models predict similar global rates of gross primary productivity (GPP) but have approximately two‐fold to three‐fold differences in their underlying global mineralization, immobilization, plant N uptake, nitrification, and denitrification fluxes. Both models dramatically underestimate the immobilization of NO 3 − by soil bacteria compared with literature values and predict dominance of plant uptake by a single form of mineral nitrogen (NO 3 − for ELM, with regional exceptions, and NH 4 + for CLM5.0). CLM5.0 strongly underestimates the global ratio of gross nitrification:gross mineralization and both models are likely to substantially underestimate the ratio of nitrification:denitrification. Few experimental data exist to evaluate this last ratio, in part because nitrification and denitrification are quantified using different techniques and because denitrification fluxes are difficult to measure at all. More observational constraints on soil nitrogen fluxes such as nitrification and denitrification, as well as greater scrutiny of the functional impact of introducing separate NH 4 + and NO 3 − pools into ESMs, could help to improve confidence in present and future simulations of N limitation on the carbon cycle.

Nevison, Cynthia↗
DOE OSTI2020

A pyrolysis study of allylic hydrocarbon fuels

The pyrolysis of selected C 3 –C 5 allylic hydrocarbons has been studied using a single-pulse shock tube. A new single-pulse shock tube has been designed and constructed by recommissioning an existing conventional shock tube. This facility enables the investigation of high-temperature chemical kinetics with an emphasis on combustion chemistry. The modifications performed on the existing shock tube are described, and the details of the sampling system to analyze the species concentration using a gas chromatography-mass spectrometry-flame ionization detection (GC-MS with a flame ionization detector) system are also provided. This facility is characterized and validated by performing cyclohexene pyrolysis experiments. Furthermore, the performance of the shock tube is demonstrated by reproducing previous literature measurements on the pyrolysis of isobutene. Postvalidation, this setup is used to study the pyrolysis of trans-2-butene and 2-methyl-2-butene (2M2B). A newly developed mechanism, NUIGMech1.0, is used to simulate the experimental data of propene, isobutene, 2-butene, and 2M2B, allylic hydrocarbon fuels. A description using two different kinetic simulation approaches is provided using our isobutene experiments as a reference. We found no significant differences between the two methods. Additionally, the contribution of different reaction classes on fuel consumption is detailed and the influence of geometry on fuel consumption and first aromatic ring: benzene is discussed.

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

Elucidating structure and dynamics of crystalline α-zirconium phosphates intercalated with water and methanol by multinuclear solid-state MAS NMR: A comprehensive NMR approach

Solid-state NMR experiments on 2 H, 31 P, 13 C, and 1 H nuclei, including 31 P T 1 , 1 H T 1 , and 1 H T 1ρ measurements, as well as on the kinetics of proton-phosphorus cross-polarization have been performed to characterize the crystalline and amorphous α-zirconium phosphates, which were intercalated with D 2 O and/or CD 3 OD. The 13 C{ 1 H} CP MAS NMR experiment performed for compound 1-CD 3 OD (Zr (HPO 4 ) 2 .0.2CD 3 OD) with carbon cross-polarization via protons of phosphate groups has provided a prove that the methanol was intercalated into the interlayer spaces of this compound. The variable-temperature 2 H solid-echo MAS NMR spectra of intercalated compounds demonstrated that the methanol molecules, in contrast to the mobile water, were immobile, keeping, however, free CD 3 rotations around the C 3 -axis. It has been demonstrated that the intercalated species, D 2 O and CD 3 OD, do not affect the high-frequency motions of the phosphate groups. By utilizing local structural models that satisfy the constraints of the experimental data, it has been suggested that the immobile methanol molecules are located in the cavity between two neighboring layers of the zirconium phosphates. Here, the present work illustrates the reliable criteria in a comprehensive NMR approach to structural and dynamic studies of such systems.

1H↗
DOE OSTI2021

Evolution of anisotropy with saturation and its implications for the elastoplastic responses of clay rocks

Many clay rocks have distinct bedding planes. Experimental studies have shown that their mechanical properties evolve with the degree of saturation (DOS), often with higher stiffness and strength after drying. For transversely isotropic rocks, the effects of saturation can differ between the bed-normal (BN) and bed-parallel (BP) directions, which gives rise to saturation-dependent stiffness and strength anisotropy. Accurate prediction of the mechanical behavior of clay rocks under partially saturated conditions requires numerical models that can capture the evolving elastic and plastic anisotropy with DOS. In this study, we present an anisotropy framework for coupled solid deformation-fluid flow in unsaturated elastoplastic media. Here we incorporate saturation-dependent strength anisotropy into an anisotropic modified Cam-Clay (MCC) model and consider the evolving anisotropy in both the elastic and plastic responses. The model was calibrated using experimental data from triaxial tests to demonstrate its capability in capturing strength anisotropy at various levels of saturation. Through numerical simulations, we demonstrate the role of evolving stiffness and strength anisotropy in the mechanical behavior of clay rocks. Plane strain simulations of triaxial compression tests were also conducted to demonstrate the impacts of material anisotropy and DOS on the mechanical and fluid flow responses.

36 MATERIALS SCIENCE↗