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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Hole Trapping in Halide Perovskites Induces Phase Segregation

Metal halide perovskites have garnered a great deal of attention for their applications in photovoltaics, LEDs, and radiation detection. The ease of solution processing high-quality perovskite semiconductors with large absorption coefficients and tolerance to native defects is decidedly attractive. Additionally, the ability to precisely tune the band gap of halide perovskites through compositional alloying of the halide ion is of particular interest for a range of applications, especially for tandem solar cells. However, under steady state light irradiation, an initially homogeneous mixed halide perovskite (MHP) will form local domains that are rich in one halide ion (e.g., Br or I). This light-induced phase segregation in MHPs forms iodide-rich domains that act as charge carrier traps and lowers the efficiency of perovskite-based devices. Thus, phase segregation poses a serious challenge to the implementation of MHPs in real-world device settings. Interestingly, when a phase segregated MHP film is placed in the dark, entropic driving forces become dominant and the segregated perovskite remixes and returns to its initially homogeneous state. Several key mechanistic details of phase segregation have been elucidated over the years. However, there are still aspects of halide segregation that are not clear, and there is ongoing debate in the literature as to what are the key factors that contribute to the mechanism. This Account discusses recent results that point to the specific role of hole trapping in phase segregation. Interestingly, generation of holes through above-band-gap excitation or through electrochemical injection increases ion migration and leads to phase segregation. The thermodynamic and redox properties of halide perovskites provide a strong driving force for hole trapping and oxidation of iodide species in MHPs. However, mobile halide species within the perovskite lattice take time to migrate and generate halide-rich domains. When in contact with a nonpolar solvent, the migration of iodine species is further extended to expulsion of iodine from the perovskite film. Thus, the mobility of halides and their susceptibility to hole-induced oxidation play a crucial role in determining the long-term stability of metal halide perovskites. Strategies to gain kinetic control over ion migration to slow phase segregation are needed to overcome these hurdles and achieve stable mixed halide perovskites. Modification of the perovskite composition through introduction of different cations or halide ions, or introduction of low-dimensional perovskite phases may suppress phase segregation. Furthermore, in achieving stability and improving the efficiency of perovskite solar cells and light emitting devices with minimal impacts, suppression of segregation remains the key factor.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Thermodynamics of Tritium Trapping by Point Defects at Interfacial Cr-based Phases of the Al Coating

Density functional theory (DFT) simulations have been carried out to evaluate the potential for tritium trapping by metal vacancies in four phases Cr-containing phases (i.e., Cr3Si, Al0.3Cr0.7, Al8Cr5-HT, and Al8Cr5) 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 T2-partial pressure dependence on the thermodynamics of singly tritiated defects. Key results in this work suggest that metal vacancies in the four phases have the potential to favorably trap tritium species, especially Si vacancies in Cr3Si phase. This overall thermodynamic trend can be correlated to the energy cost of having an interstitial tritium in the lattice, which has been calculated to range from 0.17 eV in Al8Cr5-HT to 0.89 eV in Cr3Si. A comparison of the results from this study with previous theoretical works investigating tritium behavior in other Fe-Al phases identified in the aluminide coating, suggests that metal vacancies are generally able to trap tritium in various Fe-Al aluminide phases. Especially, it was found that Si and Al vacancies would be the most efficient to trap for tritium, followed by Cr vacancies, then Fe and Ni vacancies. In the four Cr-based material phases investigated in this work, it is interesting to note that, in the absences of Fe or Ni species, strong interactions between tritium and the metal vacancies are always occurring by the formation of preferential Cr—T bonds (i.e., no Si—T or Al—T bonds were formed).

Sassi, Michel J.↗
DOE OSTI2022

Conjugated Polyelectrolyte-Based Complex Fluids as Aqueous Exciton Transport Networks

The ability to assemble artificial systems that mimic aspects of natural light-harvesting functions is fascinating and attractive for materials design. Given the complexity of such a system, a simple design pathway is desirable. Here, we argue that associative phase separation of oppositely charged conjugated polyelectrolytes (CPEs) can provide such a path in an environmentally benign medium: water. We find that complexation between an exciton–donor and acceptor CPE leads to formation of a complex fluid. We interrogate exciton transfer from the donor to the acceptor CPE within the complex fluid and find that transfer is highly efficient. We also find that excess molecular ions can tune the modulus of the inter-CPE complex fluid. Even at high ion concentrations, CPEs remain complexed with significantly delocalized electronic wavefunctions. In conclusion, our work lays the rational foundation for complex, tunable aqueous light-harvesting systems via the intrinsic thermodynamics of associative phase separation.

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

Influence of processing on secondary phase formation and microstructural evolution at U-10Mo alloy and Zr interlayer interfaces

Thermo-mechanical processing of uranium-10 wt. % molybdenum (U-10Mo) fuel plates leads to microstructure changes at the U-10Mo/Zr interfaces. Secondary phases formed at this interface are particularly important to interfacial bond strength, process optimization, and maintaining structural integrity of the U-10Mo fuel plates during irradiation. In this work, we determined the phases and phase transformation products occurring at the interface of the U-10Mo fuel and Zr interlayer when the fuel plate is subjected to short and long hot isostatic pressure times. Interfacial morphology, structure and composition of phases formed, and relative hardness across the U-10Mo/Zr interfaces were studied using a multi-length scale, multi-modal characterization approach involving electron microscopy, atom probe tomography, and atomic force microscopy. Here, results highlight that the extent of phase transformations, secondary phase formation, and hardness variability across interfaces can be controlled by modifying processing parameters. Phase diagram construction and thermodynamic calculations were performed using the Thermocalc software to identify expected phases formed at interfaces during the maximum hold temperature of 560 °C experienced during HIP.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Experimental investigation and thermodynamic modeling of the binary Pt–Sn system

Here, the thermodynamics of the Pt–Sn system is revisited with attention to phase equilibria involving PtSn 4 , a topological quantum material with interesting physical properties. The phase relationship in the Sn-rich region of the Pt–Sn binary system was investigated by a combination of differential scanning calorimetry and equilibration experiments. Based on the present experiments and a critical evaluation of the available thermodynamic and phase diagram data reported in the literature, thermodynamic assessments of the entire Pt–Sn system were performed using the CALPHAD approach. Two sets of thermodynamic parameters were optimized, where the liquid phase was described using the Bragg-Williams random mixing model and the associate solution model, respectively. The modeling results are compared with the reliable experimental data and are shown to provide good representations of the phase equilibria and thermodynamic properties in this system.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Integrated low carbon H 2 conversion with in situ carbon mineralization from aqueous biomass oxygenate precursors by tuning reactive multiphase chemical interactions

Meeting our rising demand for clean energy carriers such as H 2 from renewable biomass resources is challenged by the co-emission of CO 2 and CH 4 . To address this challenge, we design novel reactive separation pathways that integrate multiphase chemical reactions by harnessing Ca and Mg bearing minerals as a sorbent to capture CO 2 released during the hydrothermal deconstruction of aqueous biomass oxygenates to produce H 2 and solid carbonates via low temperature aqueous phase reforming and thermodynamically downhill carbon mineralization. Earth abundant catalysts such as Ni/Al 2 O 3 are effective in producing H 2 yields as high as 79% and 74% using ethylene glycol and methanol in the presence of Ca(OH) 2 as an alkaline sorbent, without contaminating or deactivating the catalyst. H 2 yields with in situ carbon mineralization using a Ni or Pt/Al 2 O 3 catalyst are enhanced based on the following order of reactivity: acetate < glycerol < methanol < formate < ethylene glycol. These studies demonstrate that the multiphase chemical interactions can be successfully tuned to enhance H 2 yields through the selective cleavage of C–C bonds using Ni/Al 2 O 3 catalysts to deconstruct biomass oxygenates for producing H 2 and CO 2 , and in situ carbon mineralization by harnessing abundant alkaline materials, as demonstrated using ladle slag. This approach unlocks new scientific possibilities for harnessing multiple emissions including abundant organic-rich wastewater streams and alkaline industrial residues to co-produce low carbon H 2 and carbonate-bearing materials for use in construction by using renewable solar thermal energy resources.

09 BIOMASS FUELS↗
DOE OSTI2022

Density functional thermodynamic description of spin, phonon and displacement degrees of freedom in antiferromagnetic-to-paramagnetic phase transition in YNiO 3

This work herein demonstrates a direct density functional description of the finite-temperature thermodynamic properties of solids exhibiting phase transitions through positional and spin symmetry breaking degrees of freedom. A classic example addressed here is the rare-earth (R) nickelates RNiO 3 where the ground state is characterized by crystallographic and magnetic (e.g., antiferromagnetic) long-range order (LRO), whereas the higher temperature paramagnetic phase manifests a range of local spin and positional symmetry breaking motifs with short-range order (SRO). Unlike time-dependent simulations of spin and positional degrees of freedom, in the present work, phases are described via a superposition of static configurations constructed by populating a periodic base lattice supercell allowing for the formation of energy lowing distribution of positional and spin local motifs. The thermal populations of the configurations in such a superposition phase are obtained from the energy-minimized Density Functional Theory (DFT)-calculated partition functions at different temperatures. This approach offers flexible inclusion of different physical contributions to the free energy, such as elastic, electronic and phonon free energies, all obtained from the same underlying DFT total energy calculations of periodic structures. The thermodynamic and magnetic properties of both LRO and SRO crystallographic and spin phases, including antiferromagnetic (AFM) to paramagnetic (PM) Néel phase transition in YNiO3 are studied. Including spin and phonon contributions, we find a DFT-calculated Néel temperature to be 144 K in satisfactory agreement with the experimental value of 145 K; whereas omitting the phonon contribution, one obtains a Néel temperature of 81 K. We present phonon contributions to the DFT-calculated temperature-dependent SRO, heat capacities, and the polymorphous distribution of nonzero local magnetic moments in the PM phase. This approach thus extends to finite temperatures the symmetry-broken DFT description of both the AFM and PM phases, demonstrating that a thermodynamic superposition approach based on symmetry broken configurations evaluated by a mean-field like DFT is sufficient to obtain a consistent description of the thermal physics of the AFM, PM phases and their interconversion in 3d oxides illustrated by YNiO 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2020

Competing Interactions between Mesoscale Length-Scales, Order-Disorder, and Martensitic Transformation in Ferromagnetic Shape Memory Alloys

In the present study, the effects of composition and heat treatments and the resulting microstructural changes on the martensitic transformation and ferromagnetic transition have been investigated in NiCoMnIn metamagnetic shape memory alloys. In this shape memory alloy system, it is observed that upon heat treatments at a wide temperature range, the onset temperature of the martensitic transformation follows a non-monotonic behavior with respect to heat treatment time, and the nature of the non-monotonic behavior is also a function of composition. This behavior cannot be attributed to well-known factors such as precipitation, change in local composition due to the precipitation and/or global degree of order. In this work, a systematic investigation through synthesis, thermal processing and characterization via thermo-physical measurements, transmission electron microscopy and in-situ synchrotron x-ray diffraction experiments has been used to correlate the non-monotonic dependence of the martensitic transformation and ferromagnetic transition temperatures to the evolution of L2 1 domains arising from a order/disorder phase transition. A thermodynamic model for the magneto-structural transition is combined with classical nucleation theory to further ascertain the role of microstructural length-scales on the onset of the martensitic transformation. This work thus provides understanding of the thermodynamic and kinetic factors that can be controlled to tune the coupled magneto-structural transformations in NiCoMnIn metamagnetic shape memory alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Flow visualisation in real-size optical injectors of conventional, additised, and renewable gasoline blends

Research on renewable and alternative fuels is crucial for improving the energy and environmental efficiency of modern gasoline internal combustion engines. To highlight the influence of fuel rheological and thermodynamic properties on phase change and atomisation processes, three types of gasoline blends were tested. More specifically, the campaign comprised a reference gasoline, an ethanol/gasoline blend (10% v/v) representative of renewable fuels, and an additised gasoline sample treated with viscoelasticity-inducing agents. High-speed imaging of the transient two-phase flow field arising in the internal geometry and the near-nozzle spray region of gasoline injectors was performed employing Diffuse Backlight Illumination. The metallic body of a commercial injector was modified to fit transparent tips realising two nozzle layouts, namely a two-hole real size model resembling the Engine Combustion Network spray G injector and an enraged replica with an offset hole. Experiments were conducted at realistic operating conditions comprising an injection pressure of 100 bar and ambient pressures in the range of 0.1–6.0 bar to cover the entire range of chamber pressures prevailing in Gasoline Direct Injection engines. The action of viscoelastic additives was verified to have a suppressive effect on in-nozzle cavitation (6% reduction in cavitation extent) , while also enhancing spray atomisation at flash-boing conditions, in a manner resembling the more volatile gasoline/ethanol blends. Finally, persisting liquid ligaments were found to form after the end of injection for the additised sample, owing to the surfactant nature of the additives.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

Parton distribution functions from scalar light-front parton gas model

Here, we propose an application of a microcanonical ensemble with light-front kinematics to model the phase-space distribution of relativistic constituents of a bound state. These constituents denoted by partons are treated as classical spin-zero particles confined inside the bound state with inter-parton collisions as their only interaction. The microcanonical molecular dynamics ensemble is applied to obtain the phase-space distribution of such a thermodynamic system. We sample this phase-space distribution using Monte Carlo algorithms to obtain the parton distribution functions (PDFs) in scenarios with 3, 4, and 5 identical partons. In addition PDFs when a selected number of massless partons are mixed with 3 massive partons are also presented.

Microcanonical ensemble↗
DOE OSTI2020

Kinetic Pathways Templated by Low-Temperature Intermediates during Solid-State Synthesis of Layered Oxides

Layered oxides have been the dominant cathodes in lithium-ion batteries, and among them, high-nickel (Ni) systems are attractive because of their high capacity. For practical use, synthetic control of stoichiometry and structural ordering is crucial but has been nontrivial due to the complexity inherent to synthesis reactions, which often proceed via nonequilibrium pathways. In this paper, we report here a combined in situ synchrotron X-ray diffraction and ab initio study of solid-state synthesis of layered oxides starting from acetate precursors for LiCoO 2 and LiNiO 2 and their solid solution LiNi 0.8 Co 0.2 O 2 . While all three systems ultimately evolve into the same thermodynamically stable layered phase ($R\overline{3}m$), each chemistry involves distinct metastable intermediates. We explain the phase progressions using a structural template model, demonstrating that during the synthesis of LiCoO 2 , the formed metastable spinel polymorph (Li 2 Co 2 O 4 ; $Fd\overline{3}m$) is a kinetically facile lithiation product of spinel Co 3 O 4 —the low-temperature (LT) intermediate from the decomposed Co-acetate. Similarly, in the Ni-based systems, the acetate decomposition products, rocksalts (Ni,Co)O, topotactically template the kinetic pathways of forming disordered rocksalts (Li x (Ni,Co) 2–x O 2 ; $Fm\overline{3}m$), consequently leading to off-stoichiometric Li x (Ni,Co)O 2 with undesired high Li/Ni mixing. These findings highlight new opportunities for engineering precursors to form LT intermediates that template the synthesis of target phases and structural properties.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Electrochemical Oxidation in Garnet-Type Solid Electrolyte by Formation of Point Defects

All-solid-state batteries hold greater promise for improving safety and energy density over conventional battery technology employing organic liquid electrolytes. One of the required features of a Li + conducting solid electrolyte is electrochemical stability, attained thermodynamically or kinetically, within the targeted operating voltage and temperature ranges. Therefore, understanding of the oxidative or reductive degradation mechanism is important to allow the design of stable solid electrolyte materials. This work contributes to building an understanding of the oxidative degradation mechanism in lithium solid electrolytes at cell operating conditions. Here, we have focused on resolving the oxidative decomposition mechanism of Al-doped lithium garnet Li 6.28 Al 0.24 La 3 Zr 2 O 12 (LLZO) as a state-of-the-art inorganic ceramic electrolyte. By combining experimental and computational analyses, we show that oxidation of LLZO occurs by simultaneous loss of oxygen and lithium from the structure, resulting in substoichiometric LLZO, at a moderate temperature (80 °C) and a high electrode potential (4.3 V vs Li/Li + ). Based on X-ray absorption and diffraction analyses, we find that the zirconium coordination shells in LLZO contract while the crystal structure experiences positive chemical strain upon electrochemical oxidation. The results from ex situ structural characterization of both the local structure and crystal symmetry are supported by a substoichiometric LLZO with lithium and oxygen vacancies, modeled by density functional theory (DFT) calculations. These chemical and structural changes in LLZO suppress effective lithium-ion conductivity by an order of magnitude. Formation of lithium and oxygen vacancies in LLZO upon electrochemical oxidation is different from prior thermodynamic predictions of phase decomposition of LLZO. The difference here is that the experiments were conducted at near-room temperature, which can hinder the kinetics of phase separation, and thus, the resultant LLZO solid electrolyte is still single-phase but substoichiometric in Li and O. In conclusion, these findings contribute an important degradation mechanism of the electrolyte, relevant for practical operational conditions of solid-state batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Polymorphic Control of CsGaE 2 (E = S, Se) Using Diorganyl Dichalcogenide Precursors

Diorganyl dichalcogenides (R–E–E–R; R = organyl, E = S, Se, Te) are widely employed for the polymorphic control of diamond-like crystal structures owing to their tunable C–E bond reactivity. In this work, we extend their application to non-close-packed alkali metal-based chalcogenides by demonstrating the selective synthesis of CsGaE 2 (E = S, Se) polymorphs. Weaker C–E bonds or higher temperatures promote the formation of the thermodynamically stable mC64 phases, while stronger C–E bonds or lower temperatures allow for the isolation of the metastable mC16 phases. Furthermore, heat-treatment studies reveal the pronounced metastability of CsGaSe 2 -mC16 particles and the irreversible phase transition of CsGaS 2 analogous to its bulk counterpart.

Cancer↗
DOE OSTI2022

Ion-Assisted Ligand Exchange for Efficient and Stable Inverted FAPbI 3 Quantum Dot Solar Cells

Perovskite quantum dot (QD) has emerged as a promising material for photovoltaics with its superior stability compared to their three-dimensional bulk counterparts, owing to its thermodynamically stabilized photoactive phase. However, ligand management on perovskite QD surfaces is extremely difficult due to the ionic nature of the perovskite lattice, which often leads to either incomplete removal of the native insulating ligands or formation of trap states during the ligand-exchange process, greatly hampering the photovoltaic performances. In this work, we report an ion-assisted ligand-exchange method for FAPbI 3 QDs using AOAc (A = formamidinium (FA + ), guanidinium (GA + ), and phenethylammonium (PEA + ), OAc = acetate), with the A + and OAC - ions promoting the removal of native long-chain insulating ligands. Further, the more complete ligand exchange results in dense and well-oriented packing of QDs, together with the enhancement of electronic coupling and charge transport across QDs. In addition, the A + and OAC - ions can fill the surface A-/X-site vacancies, respectively, reducing the QD surface trap state density and hence suppressing charge recombination and iodide migration. The p-i-n inverted QD solar cells fabricated with this ligand-exchange method exhibit significant enhancement in shortcircuit current density (J SC ), reaching a PCE of 10.13%. Moreover, unencapsulated devices show impressive stabilities of more than 7300 h (10 months) storage time in a N 2 -filled glovebox.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Liquid–Solid Mixtures of Ga Metal Infused with Cu Microparticles and Nanoparticles for Microscale and Nanoscale Patterning of Solid Metals at Room Temperature

This paper demonstrates and characterizes a simple ink for nanopatterning of solid metallic structures under ambient conditions by taking advantage of the low melting point of gallium and its affinity to form intermetallics with other metals. Bare copper micro- and nanoparticles readily mix with liquid gallium near room temperature to form a paste that ultimately amalgamates into a dense solid when mixed at the appropriate concentration ratios. The paste has sufficient shelf life that can be extruded from a nozzle at modest temperatures to form solid three-dimensional (3D) shapes. Additionally, the paste can be molded at room temperature to replicate feature sizes ranging from a few millimeters down to hundreds of nanometers. In situ X-ray diffraction (XRD) and thermo-mechanical analysis (TMA) data show gallium and copper readily interdiffuse to form the thermodynamically expected intermetallic phase. We describe the capabilities and limitations of a simple way to pattern solid metals in an additive fashion (syringe-based extrusion) and with high resolution (molding) at or near room temperature. Here, the use of a paste that solidifies provides a novel route for 3D printing of solid metals at ambient temperatures as well as the creation of micro- and nanostructured metallic surfaces that may be useful for optics, non-wetting surfaces, or electronic microcomponents.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Substitution-Mediated Calcination of Nickel-Based Cathodes: Decoupling Lithiation and Crystallization

Nickel-based layered cathodes such as LiNiO 2 offer high energy density for lithium-ion batteries, yet improvements in cycling performance and safety are required for practical use–often achieved through manganese and cobalt substitution as in LiNi 0.80 Mn 0.10 Co 0.10 O 2 (NMC811). However, how such substitution impacts calcination, the key process that governs lithiation, structural ordering, crystallization, and ultimately the resulting material properties, remains unclear. Here, in this study, we investigate substitution-mediated calcination dynamics in NMC811 compared to LiNiO 2 using multiscale-correlated in situ spectroscopy and atomistic-to-mesoscale modeling. While both systems progress through the same sequence of intermediates toward the thermodynamically favored layered phase, NMC811 exhibits an earlier onset of layering, concurrent with hydroxide decomposition followed by sluggish crystallization. Modeling reveals that Mn and Co lower the energy barrier for lithium incorporation and ordering but increase the penalty for interlayer gliding, thereby slowing crystal growth at elevated temperatures. This substitution-mediated decoupling of lithiation and crystallization explains the fine-grained microstructure observed in NMC811 versus coarsened particles in LiNiO 2 and establishes a mechanistic framework for predictive microstructure engineering of Ni-based cathodes.

Calcination↗
DOE OSTI2021

Enhanced Visibility of Subduction Slabs by the Formation of Dense Hydrous Phase A

Abstract Phase A (Mg 7 Si 2 O 8 (OH) 6 ) is one of the important dense hydrous magnesium silicates in subducting slab, because it forms after the breakdown of antigorite serpentine and could be the dominant hydrous phase in the upper‐mantle deep slab for the water transportation into deep Earth. In this study, the compressional ( P ) and shear ( S ) wave velocities of phase A were measured at simultaneous pressure and temperature conditions up to 11 GPa and 1073 K. Combined with elastic properties of olivine and pyroxenes, we calculate the hydration effect on the velocities of harzburgite lithology in cold subduction zones throughout the depth ranges where phase A is thermodynamically stable. Our calculations suggest that the hydration increases both P ‐ and S ‐wave velocities of harzburgite; at ∼5 wt% hydration, its seismic detectability is enhanced by 1%–1.5% in velocity contrasts relative to its anhydrous counterpart.

Cai, Nao↗
DOE OSTI2020

A higher-order finite element reactive transport model for unstructured and fractured grids

Abstract This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.

58 GEOSCIENCES↗