Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermodynamic phase”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Neural network representations of multiphase Equations of State

Abstract Equations of State model relations between thermodynamic variables and are ubiquitous in scientific modelling, appearing in modern day applications ranging from Astrophysics to Climate Science. The three desired properties of a general Equation of State model are adherence to the Laws of Thermodynamics, incorporation of phase transitions, and multiscale accuracy. Analytic models that adhere to all three are hard to develop and cumbersome to work with, often resulting in sacrificing one of these elements for the sake of efficiency. In this work, two deep-learning methods are proposed that provably satisfy the first and second conditions on a large-enough region of thermodynamic variable space. The first is based on learning the generating function (thermodynamic potential) while the second is based on structure-preserving, symplectic neural networks, respectively allowing modifications near or on phase transition regions. They can be used either “from scratch” to learn a full Equation of State, or in conjunction with a pre-existing consistent model, functioning as a modification that better adheres to experimental data. We formulate the theory and provide several computational examples to justify both approaches, highlighting their advantages and shortcomings.

Science & Technology - Other Topics↗
DOE OSTI2020

Synthesis of Honeycomb-Structured Beryllium Oxide via Graphene Liquid Cells

Using high-resolution transmission electron microscopy and electron energy-loss spectroscopy, we demonstrate that beryllium oxide crystallizes in the planar hexagonal structure in a graphene liquid cell by a wet-chemistry approach. These liquid cells can feature van-der-Waals pressures up to 1 GPa, producing a miniaturized high-pressure container for the crystallization in solution. The thickness of as-received crystals is beyond the thermodynamic ultra-thin limit above which the wurtzite phase is energetically more favorable according to the theoretical prediction. The crystallization of the planar phase is ascribed to the near-free-standing condition afforded by the graphene surface. Our calculations show that the energy barrier of the phase transition is responsible for the observed thickness beyond the previously predicted limit. These findings open a new door for exploring aqueous-solution approaches of more metal-oxide semiconductors with exotic phase structures and properties in graphene-encapsulated confined cells.

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

C-FOG: Life of Coastal Fog

C-FOG is a comprehensive bi-national project dealing with the formation, persistence, and dissipation (life cycle) of fog in coastal areas (coastal fog) controlled by land, marine, and atmospheric processes. Given its inherent complexity, coastal-fog literature has mainly focused on case studies, and there is a continuing need for research that integrates across processes (e.g., air–sea–land interactions, environmental flow, aerosol transport, and chemistry), dynamics (two-phase flow and turbulence), microphysics (nucleation, droplet characterization), and thermodynamics (heat transfer and phase changes) through field observations and modeling. Central to C-FOG was a field campaign in eastern Canada from 1 September to 8 October 2018, covering four land sites in Newfoundland and Nova Scotia and an adjacent coastal strip transected by the Research Vessel Hugh R. Sharp . An array of in situ, path-integrating, and remote sensing instruments gathered data across a swath of space–time scales relevant to fog life cycle. Satellite and reanalysis products, routine meteorological observations, numerical weather prediction model (WRF and COAMPS) outputs, large-eddy simulations, and phenomenological modeling underpin the interpretation of field observations in a multiscale and multiplatform framework that helps identify and remedy numerical model deficiencies. An overview of the C-FOG field campaign and some preliminary analysis/findings are presented in this paper.

58 GEOSCIENCES↗
DOE OSTI2021

Ice and Supercooled Liquid Water Distributions over the Southern Ocean based on In Situ Observations and Climate Model Simulations

Three climate models are evaluated using in situ airborne observations from the Southern Ocean Clouds, Radiation, Aerosol Transport Experimental Study (SOCRATES) campaign. The evaluation targets cloud phases, microphysical properties, thermodynamic conditions, and aerosol indirect effects from –40°C to 0°C. Compared with 580-s averaged observations (i.e., 100 km horizontal scale), the Community Atmosphere Model version 6 (CAM6) shows the most similar result for cloud phase frequency distribution and allows more liquid-containing clouds below –10°C compared with its predecessor—CAM5. The Energy Exascale Earth System Model (E3SM) underestimates (overestimates) ice phase frequencies below (above) –20°C. CAM6 and E3SM show liquid and ice water contents (i.e., LWC and IWC) similar to observations from –25°C to 0°C, but higher LWC and lower IWC than observations at lower temperatures. Simulated in-cloud RH shows higher minimum values than observations, possibly restricting ice growth during sedimentation. As number concentrations of aerosols larger than 500 nm (Na 500 ) increase, observations show increases of LWC, IWC, liquid, and ice number concentrations (N liq , N ice ). Number concentrations of aerosols larger than 100 nm (Na 100 ) only show positive correlations with LWC and N liq . From –20°C to 0°C, higher aerosol number concentrations are correlated with lower glaciation ratio and higher cloud fraction. From –40°C to –20°C, large aerosols show positive correlations with glaciation ratio. CAM6 shows small increases of LWC and N liq with Na 500 and Na 100 . E3SM shows small increases of N ice with Na 500 . Altogether, CAM6 and E3SM underestimate aerosol indirect effects on ice crystals and supercooled liquid droplets over the Southern Ocean.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Peptide programming of supramolecular vinylidene fluoride ferroelectric phases

Ferroelectric structures have spontaneous macroscopic polarization that can be inverted using external electric fields and have potential applications including information storage, energy transduction, ultralow-power nanoelectronics and biomedical devices. These functions would benefit from nanoscale control of ferroelectric structure, the ability to switch polarization with lower applied fields (low coercive field) and biocompatibility. Soft ferroelectrics based on poly(vinylidene fluoride) (PVDF) have a thermodynamically unstable ferroelectric phase in the homopolymer, complex semi-crystalline structures, and high coercive fields. Here, in this work, we report on ferroelectric materials formed by water-soluble molecules containing only six VDF repeating units covalently conjugated to a tetrapeptide, with the propensity to assemble into the β-sheet structures that are ubiquitous in proteins. This led to the discovery of ribbon-shaped ferroelectric supramolecular assemblies that are thermodynamically stable with their long axes parallel to both the preferred hydrogen-bonding direction of β-sheets and the bistable polar axes of VDF hexamers. Relative to a commonly used ferroelectric copolymer, the biomolecular assemblies exhibit a coercive field that is two orders of magnitude lower, as the result of supramolecular dynamics, and a similar level of remnant polarization, despite having a peptide content of 49 wt%. Furthermore, the Curie temperature of the assemblies is about 40 °C higher than that of a copolymer containing a similar amount of VDF. This supramolecular system was created using a biologically inspired strategy that is attractive in terms of sustainability and that could lead to new functions for soft ferroelectrics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Development of a CPU/GPU portable software library for Lagrangian–Eulerian simulations of liquid sprays

The Lagrangian–Eulerian method is widely used for simulations of fuel sprays in turbulent combustion because of the advantage of treating the spray droplets as discrete points. One challenge of the Lagrangian–Eulerian method is the intense computational requirement when tracking the large number of Lagrangian particles needed for high fidelity. We have developed a performance-portable library, Grit, to track the Lagrangian particles in parallel on central processing unit (CPU) and graphics processing unit (GPU) accelerated high performance computing (HPC) architectures. Grit is a C++ library which employs Message Passing Interface (MPI) for distributed memory parallelism and Kokkos programming model for on-node shared memory parallelism with performance portability across different architectures of GPUs and multi-core/manycore CPUs. The parallel algorithms, key parallel kernels, and their performances on the pre-exascale supercomputer, Summit, are presented. Grit is coupled with a direct numerical simulation (DNS) solver, S3D, for multiphase simulations. A conservative formulation has been developed and implemented in Grit for phase coupling with thermodynamic consistency. The formulation separates the conservation of mass, momentum and energy from the physical models to prevent accidental violation of conservation laws due to inconsistent models. The formulation also enforces consistent definitions of enthalpies of the fuel for both phases and the latent heat of evaporation. Finally, simulations of turbulent particle-laden flow and the evaporation of dilute turbulent spray jet are performed to verify the software implementation, and to demonstrate the scalability of Grit for large-scale multiphase simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Using a Chromatographic Pseudophase Model To Elucidate the Mechanism of Olefin Separation by Silver(I) Ions in Ionic Liquids

Silver(I) ions undergo selective olefin complexation and have been utilized in various olefin/paraffin separation techniques such as argentation chromatography and facilitated transport membranes. Ionic liquids (ILs) are solvents known for their low vapor pressure, high thermal stability, low melting points, and ability to promote a favorable solvation environment for silver(I) ion–olefin interactions. To develop highly selective separation systems, a fundamental understanding of analyte partitioning to the stationary phase and the thermodynamic driving forces behind solvation is required. In this study, a chromatographic model treating silver(I) ions as a pseudophase is constructed and employed for the first time to investigate the olefin separation mechanism in silver(I) salt/IL mixtures. Stationary phases containing varying amounts of noncoordinated silver(I) salt ([Ag + ][NTf 2 – ]) dissolved in the 1-decyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl]imide ([C 10 MIM + ][NTf 2 – ]) IL are utilized to determine the partition coefficients of various analytes including alkanes, alkenes, alkynes, aromatics, aldehyde, esters, and ketones. As ligand coordination to silver(I) ions is known to lower its olefin complexation capability, this study also examines two different types of coordinated silver(I) ion pseudophases, namely, monocoordinated silver(I) salt ([Ag + (1-decyl-2-methylimidazole, DMIM)][NTf 2 – ]) and dicoordinated silver(I) salt ([Ag + (1-methylimidazole, MIM)(DMIM)][NTf 2 – ]). The extent of olefin partitioning to the coordinated silver(I) ion pseudophases over the carrier gas and IL decreased by up to two orders of magnitude. Values for enthalpy, entropy, and free energy of solvation were determined for the three silver(I) ion-containing systems. Olefin retention was observed to be enthalpically dominated, while ligand coordination to the silver(I) ion pseudophase resulted in variations for both enthalpic and entropic contributions to the free energy of solvation. Furthermore, the developed model can be used to study chemical changes that occur in silver(I) ions over time as well as identify optimal silver(I) salt/IL mixtures that yield high olefin selectivity.

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

Disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y alloys

This dataset contains input and output files from density functional theory (DFT) simulations used to study the disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y (0 ≤ x + y ≤ 1) alloys and their binary end members MnTe, MnSb, and MnBi. The alloys adopt the hexagonal NiAs-type (nickeline) structure and span ternary (MnTexSb1-x, MnTexBi1-x, MnBixSb1-x), and quaternary compositions across the full MnTe–MnSb–MnBi composition triangle. For each alloy composition, the dataset provides DFT calculations in three magnetic configurations: A-type antiferromagnetic (AFM), C-type AFM, and ferromagnetic (FM). Every magnetic configuration folder contains the fully relaxed crystal structure (CONTCAR), VASP input parameters (INCAR), and the main VASP output file (OUTCAR), from which total electronic energies, Mn magnetic moments, lattice parameters, and percent volume changes between magnetic states are extracted. These data are used to construct compositional phase diagrams, evaluate thermodynamic stability (formability), and map magnetoelastic responses across the alloy space. For A-type AFM and FM configurations, additional data are provided as follows: (i) FORCE_CONSTANTS and thermal_properties.yaml files at the top level of A-type_AFM/ and FM/ folders — present only for compositions marked with an asterisk (*) in Table I of the main text. These are derived from Phonopy finite-displacement calculations on full disordered 128-atom supercells and provide vibrational free energy, entropy (Svib)contribution from explicit disorder calculations. (Table I of the associated main manuscript) (ii) A VCA/ subfolder within A-type_AFM/ and FM/, containing FORCE_CONSTANTS and thermal_properties.yaml from Virtual Crystal Approximation phonon calculations (without spin-orbit coupling). VCA data are available for all compositions and are used to estimate vibrational contributions to the Gibbs free energy across the full composition space. (iii) A SOC/ subfolder containing CONTCAR, INCAR, and OUTCAR from spin-orbit coupling calculations, providing relativistic corrections to electronic energies and lattice parameters (Tables S2–S3 of the SM, and Table I of the main manuscript). (iv) A SOC/VCA/ subfolder containing FORCE_CONSTANTS and thermal_properties.yaml from VCA phonon calculations performed within the SOC framework, combining relativistic and vibrational thermodynamic corrections. The computed properties are used to map the AFM–FM magnetic crossover near MnTe0.75Sb0.25, demonstrate disorder- and spin-induced phonon broadening, identify a semiconductor-to-metal crossover, and quantify the pronounced magnetoelastic volume response near the magnetic phase boundary.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Mo Atom Rearrangement Drives Layer-Dependent Reactivity in Two-Dimensional MoS 2

Two-dimensional (2D) materials offer a valuable platform for manipulating and studying chemical reactions at the atomic level, owing to the ease of controlling their microscopic structure at the nanometer scale. While extensive research has been conducted on the structure-dependent chemical activity of 2D materials, the influence of structural transformation during the reaction has remained largely unexplored. In this work, we report the layer-dependent chemical reactivity of MoS 2 during a nitridation atomic substitution reaction and attribute it to the rearrangement of Mo atoms. Our results show that the chemical reactivity of MoS 2 decreases as the number of layers is reduced in the few-layer regime. In particular, monolayer MoS 2 exhibits significantly lower reactivity compared with its few-layer and multilayer counterparts. Atomic-resolution transmission electron microscopy (TEM) reveals that MoN nanonetworks form as reaction products from monolayer and bilayer MoS 2 , with the continuity of the MoN crystals increasing with layer number, consistent with the local conductivity mapping data. The layer-dependent reactivity is attributed to the relative stability of the hypothetically formed MoN phase, which retains the number of Mo atomic layers present in the precursor. Specifically, the low chemical reactivity of monolayer MoS 2 is attributed to the high energy cost associated with Mo atom diffusion and migration necessary to form multilayer Mo lattices in the thermodynamically stable MoN phase. In conclusion, this study underscores the critical role of lattice rearrangement in governing chemical reactivity and highlights the potential of 2D materials as versatile platforms for advancing the understanding of materials chemistry at the atomic scale.

Chemical reactivity↗
DOE OSTI2025

Tritium Trapping Thermodynamics by Point Defects at Interfacial Fe-Al Phases of the Aluminide Coating

Density functional theory (DFT) simulations have been carried out to evaluate the potential for tritium trapping by metal vacancies in four phases (i.e., (Fe, Cr, Ni), AlFe 3 , NiAl, and AlFe) identified near the interface between the Al coating and 316 stainless steel (316 SS) cladding. In addition, an ab initio thermodynamics approach has been employed to predict the temperature and T 2 -partial pressure dependence on the thermodynamics of tritiated defects. Key results in this work suggest that metal vacancies in the four phases have the potential to favorably trap tritium species. This thermodynamic trend can be correlated to the high energy cost of having an interstitial tritium in the lattice, which has been calculated to be at least 0.48 eV. By combining the results of this study with previous theoretical works investigating tritium behavior in other Fe-Al phases identified in the aluminide coating, suggest that metal vacancies are generally able to trap tritium in various Fe-Al aluminide phases. Especially, it was found that Al vacancies are the most efficient to trap tritium, followed by Fe vacancies, then Ni and Cr vacancies. Strong interactions between tritium and the metal vacancies are occurring by the formation of Fe—T bonds. The formation of Cr—T, Ni—T, or Al—T bonds are found less energetically favorable than Fe—T bonds.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Thermodynamic modeling of calcium carbonate scale precipitation: aqueous Na + -Ca 2+ -Cl – -HCO 3 – -CO 3 2– -CO 2 system

To allow for accurate calculations of calcium carbonate scaling in highly saline produced waters, we present a comprehensive thermodynamic model based on the electrolyte nonrandom two-liquid (eNRTL) activity coefficient equation for the aqueous Na ⁺ -Ca ²⁺ -Cl – -HCO 3 – -CO 3 ²– -CO 2 system. The eNRTL binary interaction parameters for the H 2 O:(Na ⁺ -CO 3 2– ) pair, the H 2 O:(Na ⁺ -HCO 3 – ) pair, the (Na⁺-Cl–):(Na ⁺ -CO 3 2– ) pair, and the (Na ⁺ -Cl – ):(Na ⁺ -HCO 3 – ) pair are identified in this work via the regression of thermodynamic, calorimetric, and phase equilibria experimental data. The binary interaction parameters associated with the H 2 O:(Na ⁺ -Cl – ) pair, the CO 2 :(Na ⁺ -Cl – ) pair, the H 2 O:(Ca ²⁺ -Cl – ) pair, and the (Na ⁺ -Cl – ):(Ca ²⁺ -Cl – ) pair are retrieved from the literature. The remaining binary interaction parameters are retrieved from Aspen Plus or set to zero. In addition, the solubility product constants are identified for Na 2 CO 3 ·10H 2 O (s) , Na 2 CO 3 · 7H 2 O (s) , Na 2 CO 3 ·H 2 O (s) , Na 2 CO 3 ·NaHCO 3 2H 2 O (s) , Na 2 CO 3 ·3NaHCO 3 (s) , and CaCO 3(s) via regression of solubility data. Here, the model is capable of accurately calculating all phase equilibria and calorimetric properties at temperatures up to 473.15 K and salt concentrations up to saturation.

42 ENGINEERING↗
DOE OSTI2021

Various Low-Symmetry Phases in High- χ and Conformationally Asymmetric PDMS- b -PTFEA Copolymers

Since the discovery of the Frank–Kasper (FK) phase a decade ago, single-component block copolymer systems have offered unique opportunities to explore the structural and phase-transition properties of topologically close-packed FK phases composed of point-particle motifs. Here, we present the formation of complex low-symmetry phases and their thermodynamic stability using high-χ and high-ε linear diblock copolymers, where χ and ε are the Flory–Huggins interaction parameter and conformational asymmetry between the two blocks, respectively. In this work, a series of polydimethylsiloxane-b-poly(2,2,2-trifluoroethyl acrylate) (PDMS-b-PTFEA) copolymers are prepared by tuning the volume fraction ($f_{\text{PTFEA}}$) of a minor PTFEA block. The packing structures of PDMS-b-PTFEAs are mapped in the space of temperature versus $f_{\text{PTFEA}}$ using X-ray scattering measurements. In addition to the packing structures reported earlier, we identify another A15 phase at higher $f_{\text{PTFEA}}$ close to a hexagonally packed cylinder (HEX) phase and establish the phase sequence of HEX–A15–σ–C14–BCC with decreasing $f_{\text{PTFEA}}$. Theoretical calculation of relative free-energy densities confirms the experimental phase sequence of the composition-dependent sphere-packing structures. Interestingly, the domain sizes of three-dimensional (3D) sphere-packing structures are nearly temperature-invariant compared to that of a two-dimensional HEX structure, and we attribute this dimensional stability to the strong enthalpic effects in the formation of 3D packing structures and the geometrical characteristics of spherical domains. Our results demonstrate that the interaction parameter χ, together with the known role of conformational asymmetry ε, serves as the primary thermodynamic parameter for selecting and stabilizing various low-symmetry packing structures in sphere-forming linear diblock copolymer systems.

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

Mixed-Chalcogen 2D Silver Phenylchalcogenides (AgE 1– x E x Ph; E = S, Se, Te)

Alloying is a powerful strategy for tuning the electronic band structure and optical properties of semiconductors. Here, we investigate the thermodynamic stability and excitonic properties of mixed-chalcogen alloys of two-dimensional (2D) hybrid organic-inorganic silver phenylchalcogenides (AgEPh; E = S, Se, Te). Using a variety of structural and optical characterization techniques, we demonstrate that the AgSePh-AgTePh system forms homogeneous alloys (AgSe 1-x Te x Ph, 0 ≤ x ≤ 1) across all compositions, whereas the AgSPh-AgSePh and AgSPh-AgTePh systems exhibit distinct miscibility gaps. Density functional theory calculations reveal that chalcogen mixing is energetically unfavorable in all cases, but comparable in magnitude to the ideal entropy of mixing at room temperature. Because AgSePh and AgTePh have the same crystal structure (which is different from AgSPh), alloying is predicted to be thermodynamically preferred over phase separation in the case of AgSePh-AgTePh, whereas phase separation is predicted to be more favorable than alloying for both the AgSPh-AgSePh and AgSPh-AgTePh systems, in agreement with experimental observations. Homogeneous AgSe1-xTexPh alloys exhibit continuously tunable excitonic absorption resonances in the ultraviolet-visible range, while the emission spectrum reveals competition between exciton delocalization and self-trapping behavior. Altogether, these observations provide new insight into the thermodynamics of 2D silver phenylchalcogenides and the effect of lattice composition on electron-phonon interactions in 2D hybrid organic-inorganic semiconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Theoretical predictive screening of noble-metal-containing M 3 AuC 2 (M = Ti, V, and Cr) MAX phases

Using first-principles density functional theory (DFT), the ground state physical properties of the newly synthesized noble-metal-containing Ti 3 AuC 2 MAX phase have been investigated. The effect of transition-element Ti replacement (V and Cr) on physical properties, including structural, elastic, electronic, thermal, and optical are presented. The optimized lattice parameters of Ti 3 AuC 2 MAX phase are in good agreement with the experimental values and decrease when we replace Ti with V and then Cr. The magnetic properties of M 3 AuC 2 are predicted by using generalized gradient approximation (GGA) incorporating onsite electron correlation parameter (Hubbard U parameter) within DFT. Out of three studied MAX phases, V 3 AuC 2 and Cr 3 AuC 2 phases exhibit magnetism. The thermodynamical stability of M 3 AuC 2 (M = Ti, V, and Cr) phases is discussed by the formation enthalpy with respect to their most competing phases. The electronic structure reveals that these phases have metallic nature and are electrically anisotropic. The M-element replacement has an effect on the bonding properties of M 3 AuC 2 . The bonding between M-C is in the order of Cr 3 AuC 2 > V 3 AuC 2 > Ti 3 AuC 2 according to their peak positions and heights of density of states in the occupied site, which is also confirmed from the charge density distribution. The elastic properties indicate that M 3 AuC 2 phases are ductile, machinable, less stiff, and better resistant to thermal shock and are elastically anisotropic. To understand the properties of M 3 AuC 2 for the extreme environment, we employed a quasi-harmonic Debye model at the pressure and temperature range of 0–50 GPa and 0–1600 K, respectively. Additionally, we evaluated thermal conductivity and melting temperature to further understand the potential of these MAX phases in coating applications at elevated temperature. In relation to optical properties, these MAX phases have a reasonable absorption coefficient in the visible as well as in the ultra-violet regions. The reflectivity of M 3 AuC 2 phases is polarization dependent and of 45 % against the visible region, which reveals its potential as coating material to minimize solar heating.

312 MAX phases↗
DOE OSTI2023

Stability study of aqueous foams under high-temperature and high-pressure conditions relevant to Enhanced Geothermal Systems (EGS)

In Enhanced Geothermal System (EGS) hydraulic fracturing is carried out by injecting cold water into deep Hot Dry Rocks (HDR) under carefully controlled conditions to create new or reopen existing fractures. Water-based fracturing fluids demonstrate some challenges including immense quantity of water usage, water sensitivity of the formations, water blocking, and lack of proppant carrying capacity and transportation. Thus, an alternative is to use foam-based fracturing fluid which offers potential advantage over conventional water-based fracturing fluid such as minimum water usage, reduced wellbore damage, high proppant carrying capacity, and less environmental damage. However, foams are complex mixture of gaseous phase and liquid phase which are thermodynamically unstable at downhole conditions, and their stability can decrease over time due to liquid drainage, bubble coarsening, and coalescence. This paper shows laboratory experiments executed to study foam stability at high temperature (200°C) and high pressure (6.9MPa) conditions which simulates the geothermal environment. Foam stability was characterized by half-life of foam, which is defined as the time taken by the foam to decreases by 50% of its original height due to drainage. In this paper, two types of gaseous phases, nitrogen (N 2 ) and carbon dioxide (CO 2 ) were investigated. Also, based on successful practice of foam-based fracturing fluid in oil and gas industries, four surfactants, including Alpha olefin sulfonate (AOS), Sodium dodecyl sulfonate (SDS), Tergitol™ (NP – 40), and Cetyltrimethylammonium chloride (CTAC) at optimum concentration of 1wt.% were tested for best stability performance. In addition, different stabilizing agents including guar gum, bentonite clay, crosslinker, silicon dioxide nanoparticles (SiO 2 ), graphene oxide (GO) were also studied. Experimental results showed that N 2 foams were more stable than CO 2 foams. It was observed that foam half-life decreased with the increase in temperature. Among all the surfactants, AOS foams showed the most promising thermal stability at high temperatures. Moreover, with the addition of stabilizing agents, foam's half-life was enhanced. Stabilizing agents such as crosslinker and GO dispersion showed the most stable foams with half-life recorded at 20min and 17min, respectively, at 200°C and 6.9MPa. Finally, pressure also showed a positive effect on foam stability; with increased pressure, foam half-life was increased. Based on the experimental data, analytical models for the effect of temperature and pressure were developed, considering foam degradation is a first-order kinetic reaction that linearly depends on the foam drainage mechanism. The effect of temperature on foam half-life was studied as an exponential decay model. In this model, foam half-life is a function of drainage rate constant (D A ) and activation energy (E a ) of the foam system. The effect of pressure on foam half-life was found to obey a power-law model where an increase in pressure showed an increase in foam half-life. Furthermore, a linear relation was studied for the effect of pressure on foam activation energy and drainage rate. Then the combined effects of temperature and pressure were studied, which yielded an analytical model to predict the foam stabilities in terms of half-life for different foam compositions. In conclusion, this research indicates that with an appropriate selection of surfactants and stabilizing agents, it is possible to obtain stable foams, which could replace conventional water fracturing fluid under EGS conditions.

15 GEOTHERMAL ENERGY↗