Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gibbs”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Electronic structural and lattice thermodynamic properties of MAlO2 and M5AlO4 (M = Li, Na, K) sorbents for CO2 capture applications

Abstract The electronic properties and thermal stabilities of MAlO 2 and M 5 AlO 4 (M = Li, Na, K) are investigated by density functional theory and lattice phonon dynamics. Based on the calculated electronic and lattice thermodynamic properties, their abilities to capture CO 2 as solid sorbents are analyzed. The calculated electronic structural properties of MAlO 2 and M 5 AlO 4 indicate that all these alkali aluminates are semiconductors with a bandgap range of 2.4 ~ 6.4 eV. The 1st valence bands of these alkali aluminates are located 0 ~ − 6 eV under Fermi levels and are mainly contributed by p orbitals of O, s and p orbitals of Al and M. The phonon vibrational frequencies of M 5 AlO 4 spread at a lower frequency range compared to their MAlO 2 phases. With increasing temperature, the calculated phonon free energies of M 5 AlO 4 decrease faster than their corresponding MAlO 2 while their entropies have opposite trends. The reaction 2MAlO 2 + CO 2 = M 2 CO 3 + Al 2 O 3 has higher reaction heat and Gibbs free energy change than those of corresponding reaction 2 / 5 M 5 AlO 4 + CO 2 = M 2 CO 3 + 1 / 5 Al 2 O 3 , which shows the former reaction possesses lower turnover temperature. Among the alkali aluminates studied, the β-NaAlO 2 , lt-KAlO 2 , and γ-LiAlO 2 are better candidates that could be applied for CO 2 capture technologies. Graphical Abstract

36 MATERIALS SCIENCE↗
DOE OSTI2021

Phase evolution in two-phase alloys during severe plastic deformation

Herein, phase evolution in FCC metals with strongly interacting alloy components during severe plastic deformation is investigated using molecular dynamics simulations. Specifically, we study the alloy microstructure in steady state, nucleation and growth of precipitates in supersaturated alloys, and the decomposition of precipitates in undersaturated alloys. The results are compared to a modified effective temperature model, providing a physical understanding for the atomic processes underlying the model and a perspective on its strengths and weaknesses. Key observations in this work are nucleation and growth of precipitates during SPD at a temperature of 100 K; Gibbs-Thomson-like behavior relating solubility to precipitate size under steady-state shearing; a direct relationship between the effective temperature and the shear modulus; and the importance of cluster agglomeration during precipitate growth. The study also reveals that the mechanisms of forced chemical mixing depends on precipitate size, adding complications for effective temperature models describing inhomogeneous systems. The simulations are shown to provide good semiquantitative agreement with experimental findings reported in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Theory of strain phase equilibria and diagrams

Chemical phase equilibria and phase diagrams are well established and have served as indispensable guides for designing materials through chemical tuning. However, a general thermodynamic theory for strain phase equilibria remains elusive despite extensive studies which have revealed dramatic impacts of mechanical strain on the stability of phases and domain states in solid-state materials. Here, we establish the thermodynamic theory of strain equilibria and a general framework for efficiently constructing multidomain and multiphase diagrams of arbitrarily strained solids under incoherent conditions. As examples, we obtain temperature-strain phase diagrams of ferroelectric PbTiO 3 , strongly correlated VO 2 , and unconventional ferroelectric Hf 0.5 Zr 0.5 O 2 . We reveal the analogs of the Gibbs phase rules, multiple- and multi-critical points, common-tangent construction, and level rule in strain equilibria to those in the familiar chemical phase equilibria. Finally, our strain equilibria theory offers a powerful framework for predicting the thermodynamic stability of structural phases and domains in strained solids and for guiding the strain engineering of their functional properties.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Towards defect engineering in hexagonal MoS 2 nanosheets for tuning hydrogen evolution and nitrogen reduction reactions

Combined computational and experimental approaches were used to evaluate defective 2H-MoS 2 nanosheets for their activity and selectivity for hydrogen evolution reaction (HER) and nitrogen reduction reaction (NRR). Density functional theory calculations were used to understand the relationship between HER and NRR activity on the ideal basal MoS 2 plane, seven grain-boundaries, ten single-/few-atom vacancies and anti-sites, and zigzag and armchair edge sites. The results confirm that 2H-MoS 2 should contain several defects with high activity for HER: armchair and zigzag edges, VS vacancy, MoS 2 anti-site, and S-S and Mo-Mo grain boundaries. Considering Gibbs free energy change for all the steps in the NRR mechanism and kinetic barriers for a key NRR step, we have found that activation of perspective NRR selective sites in 2H-MoS 2 , namely V MoS6 and clusters of S-vacancies, would require large overpotential, conditions at which HER dominates. The DFT conclusions are supported by the electrochemical studies of NRR activity and selectivity under aqueous conditions, which show an increase in NRR activity but a decrease in Faradaic efficiency as applied cell potential becomes more negative. The results of this work therefore highlight the challenges in activating natural 2H-MoS 2 for NRR, which would require additional material engineering or reaction condition optimization as a way to suppress HER, decrease the NRR overpotential or preferentially both.

2H-MoS2↗
DOE OSTI2022

Probing the initial stages of iron surface corrosion: Effect of O 2 and H 2 O on surface carbonation

Iron plays a vital role in natural processes such as water, mineral, iron, and nutrient cycles. Iron undergoes reduction-oxidation and catalytic reactions to produce various corrosion films depending on its chemical environment. Near ambient pressure X-ray photoelectron spectroscopy, polarized modulated infrared reflection absorption spectroscopy, and Auger electron spectroscopy were used to study the key reactants, from O 2 (g), H 2 O vapor, Na + and Cl – on the initial stages of iron surface corrosion. Here, with increasing the ratio of O 2 and H 2 O, surface hydrocarbons were shown to oxidize into carbonates, while the Cl – was found to migrate into the interface. The effect of each individual reactant was measured separately and water was shown to have a first order rate dependence on the carbonate growth at low pressures, with little dependence for O 2 . Near ambient pressures, both H 2 O and O 2 were found to increase the carbonate growth, which was estimated using the Langmuir isotherm model, yielding Gibbs energies between –9.8 to –8.5 kJ/mol. A mechanism is suggested to explain the oxidation is catalyzed by NaCl on iron surfaces and the adventitious hydrocarbons served as the source for surface carbonation. These findings have implications for understanding other surface catalytic and redox interface chemistry in complex environments.

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

Site-specific surface reactivity on MgO for atomic layer deposition via selective hydration

Atomic layer deposition (ALD) is a powerful technique for thin film synthesis, offering atomic-scale precision and conformality. While ALD of MgO has been widely studied for applications in energy storage and microelectronics, its potential as surface on which deposition may be selective and defects repaired remains underexplored. Here, we present a combined theoretical and experimental investigation of MgO surface hydration and its implications for targeted ALD growth using water and dimethyl aluminum isopropoxide (DMAI) as reactants. We perform density functional theory (DFT) calculations to examine molecular and dissociative H 2 O adsorption on MgO (100) terraces and step-edge sites, including pristine surfaces and those with Mg/O vacancies. Reaction Gibbs free energies are calculated under various conditions to quantify surface reactivity. Our findings reveal facet- and defect-dependent hydration behaviors that align with experimental ALD growth trends on MgO (100). This study provides a molecular-level understanding of MgO surface chemistry critical for optimizing ALD processes for thin film growth and defect repair.

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

Ab initio simulations on the pure Cr lattice stability at 0K: Verification with the Fe-Cr and Ni-Cr binary systems

Significant discrepancies have been observed and discussed on the lattice stability of Cr between the predictions from the ab initio calculations and the CALPHAD approach. In the current work, we carefully examined the possible structures for pure Cr and reviewed the history back from how Kaufman originally determined the Gibbs energy of FCC-Cr in the 1970s. The reliability of Cr lattice stability derived by the CALPHAD and ab initio approaches was systematically discussed. It is concluded that the Cr lattice stability based on the CALPHAD approach has large uncertainty. Meanwhile, we cannot claim that the ab initio H FCC-Cr is error-free as FCC-Cr is an unstable phase under ambient conditions. The present work shows that the ab initio HFCC-Cr can be a viable scientific approach. As both approaches have their limitations, the present work propose to integrate the ab initio results into the CALPHAD platform for the development of the next generation CALPHAD database. As a result, the Fe-Cr and Ni-Cr binary systems were chosen as two case studies demonstrating the capability to adopt the ab initio Cr lattice stability directly into the current CALPHAD database framework.

36 MATERIALS SCIENCE↗
DOE OSTI2022

CALPHAD modeling of uranium nitride (UN) fabrication routes enabled by first-principles calculations

The thermochemical details of fabricating uranium nitride (UN) by ammonolysis of uranium tetraflouride (UF 4 ) were determined using density functional theory (DFT) and CALculation of PHAse Diagrams (CALPHAD) computational methods. The thermochemical data of all binary, ternary, and quaternary U-H-N-F phases were computed using DFT, and the data for the phases that have not been measured experimentally, including UN 2 and NH 4 F(g), were combined with existing experimentally-determined data for CALPHAD modeling. The DFT data were benchmarked using experimental Gibbs energy of reaction and experimental thermochemical data for individual species. Phase diagrams relevant to the ammonolysis reaction are depicted, showing regions of stability for solid U-N, U-F and U-N-F phases. An unidentified phase produced in a previous experiment was identified as UN 0.95 F 1.2 (UNF) by comparing its X-ray diffraction spectrum to the experimental spectrum, and its formation during the fabrication of UN from UF 4 is supported by the simulated phase diagram. Here, it is calculated that UN 2 can be produced by the ammonolysis of UF 4 , but requires elevated temperatures, high NH 3 (g) partial pressure, and large amounts of flowing NH 3 (g) to avoid solid flu oride impurities in the uranium nitride. Likewise, U 2 N 3 can be produced instead at temperatures greater than 980 K. The use of silane (SiH 4 ) gas was investigated as a potential additive in the ammonolysis fabrication route to speed removal of fluorine. The addition of SiH 4 (g) offers little advantage to the removal of fluorine, and adds the complication of Si 3 N 4 formation. The use of DFT to fill in missing data to perform CALPHAD calculations demonstrated here allows for the determination of more comprehensive and trustworthy phase diagrams than the use of existing experimental data alone.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Thermodynamic modeling of KCl-PrCl 3 and KCl-LiCl-PrCl 3 systems

Molten salt electrolysis can recover the actinides from spent nuclear fuels, and it involves a eutectic LiCl-KCl in molten form as an electrolyte. During reprocessing, the concentration of fission products such as La, Nd, Pr, etc., increases and affects the recovery efficiency of the electrolyte. In this work, thermodynamic modeling of KCl-PrCl 3 and KCl-LiCl-PrCl 3 systems was carried out using the CALPHAD (Calculation of Phase Diagrams) approach for the first time. The thermodynamic functions for the pure salts were taken from the SGTE (Scientific Group Thermodata Europe) Substances (SSUB) database. The experimental thermochemical and phase equilibria data available in the literature were used as input for the assessment of KCl-PrCl 3 and KCl-LiCl-PrCl 3 systems. The model parameters for the KCl-LiCl system were adjusted to include the new Gibbs energy descriptions for the pure salts. In addition, the sublattice model for the liquid phase in the LiCl-PrCl 3 system was modified and reassessed to ensure the model compatibility for higher-order extrapolation. There is a good agreement between the experimental and calculated thermochemical and phase diagram data for all the optimized constituent binaries and the ternary system. Furthermore, this work will be beneficial for determining the solubility limit of PrCl 3 in molten LiCl-KCl electrolytes and their thermodynamic properties for improving the efficiency of the pyrochemical process.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Thermodynamic modeling with uncertainty quantification using the modified quasichemical model in quadruplet approximation: Implementation into PyCalphad and ESPEI

The modified quasichemical model in the quadruplet approximation (MQMQA) considers the first- and the second-nearest-neighbor coordination and interactions, particularly useful in describing short-range ordering (SRO) in complex liquids such as molten salts, slag in metal processing, and electrolytic solutions. Here, the present work implements the MQMQA into the Python based open-source software PyCalphad for thermodynamic calculations. This endeavor facilitates the development of MQMQA-based thermodynamic database with uncertainty quantification (UQ) and propagation (UP) using the open-source software ESPEI. A new database structure based on Extensible Markup Language (XML) is proposed for ESPEI evaluation of MQMQA model parameters. Using the KF-NiF 2 , KCl-NaCl-MgCl 2 , and CaCl 2 -CaF 2 -LiCl-LiF salt systems as examples, we demonstrate the successful implementation of MQMQA in PyCalphad through thermodynamic calculations of Gibbs energy, equilibrium quadruplet fractions, and phase diagram, as well as database development with UQ and UP using ESPEI. Furthermore, as an application of the present implementation, both the LiF–TbF 3 and LiF-HoF 3 systems have been modeled by MQMQA for the first time, which are in good agreement with experiments. The present implementation hence offers an open-source capability for performing CALPHAD modeling for complex liquids with SRO using MQMQA plus a new XML database structure.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Enhancing perovskite electrocatalysis through synergistic functionalization of B-site cation for efficient water splitting

The family of perovskite oxides is a promising class of catalysts for diverse energy conversion processes including water splitting. In this work, a facile two-step manipulation (in-situ exsolution and post-sulfurization) strategy was proposed and applied to LaCo 0.2 Fe 0.8 O 3 (LCF) perovskite parent, through which, the electronic state, spatial immersion and intrinsic activity of B-site cobalt (Co) were stepwise tuned at nanoscale proximity accordingly (i.e., lattice Co ions segregated Co 0 → embedded CoS 2 ). Impressively, the as-prepared catalyst (S-LCF) obtains an emergent oxygen deficient microstructure seamlessly pinned with uniformly distributed CoS 2 nanoparticles (NPs), which demonstrates enhanced performance toward both oxygen evolution reaction (OER) and hydrogen evolution reaction (HER), and shows good stability in overall water splitting. The density functional theory (DFT) calculations illustrate the optimized metal-oxygen covalency and hydrogen adsorption Gibbs free energy (Δ G H *) on S-LCF, which further buttresses the prominence of our B-site cation engineering tactics.

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

Multifunctional transitional metal-based phosphide nanoparticles towards improved polysulfide confinement and redox kinetics for highly stable lithium-sulfur batteries

The shuttle effect and the sluggish redox kinetics of lithium polysulfides (LiPSs) are the major issues impeding the practical applications of lithium-sulfur batteries (LSBs). Herein, a highly-efficient Ni 2 P electrocatalyst supported on N, P co-doped graphene (Ni 2 P@NPG) is developed via a simple “recrystallization-self-assembly” method to address the above issues. Here, the ultrafine Ni 2 P nanoparticles ensure abundant adsorption-diffusion-conversion interfaces for accelerating LiPSs transformation and Li 2 S deposition, which extremely decreases the accumulation of LiPSs in the electrolyte and therefore prevents the migration of LiPSs. Their superior catalytic performance is demonstrated by reduced Gibbs free energy changes of rate-limiting step based on the systematic theoretical calculations and the reduced shuttle effect is tested by the three-dimensional reconstructions of Raman depth profiles. Benefiting from these synergistic effects, the LSBs with Ni 2 P@NPG modified separators present a superior cycling performance with an average capacity decay rate of 0.048 % per cycle at 1C around the 400 cycles and a high-rate capacity of 731 mAh/g at 2C. Even with a high-sulfur loading of 3.53 mg cm –2 , the cell can still contain a reversible capacity of 809 mAh/g at 0.2C with a remarkable columbic efficiency of 98.4 %.

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

Modeling of aqueous species interaction energies prior to nucleation in cement-based gel systems

Arguably the most ubiquitous construction material in modern civilization, concrete is enabling the development of megacities around the globe together with increasing living standards in developing nations. However, the exact formation mechanisms of the strength-giving calcium-rich gels remain a topic of debate. Using density functional modeling, we simulate the fundamental solution-based building blocks of cement hydrates (calcium ions and silicate and aluminate monomers) and their propensity to form pair-wise complexes with bonding environments characteristic of those found in calcium-silicate-hydrate, calcium-alumino-silicate-hydrate and sodium-containing calcium-alumino-silicate-hydrate gels, as assessed from Gibbs free energies of chemical reactions. By accurately simulating the high pH pore solution chemistry in Portland cements and related systems, along with discrete solvation of the species, we hypothesize potential early age formation routes of the gels and discuss limitations and future work associated with this approach.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Optimized thermodynamic properties of REE aqueous species (REE 3+ and REEOH 2+ ) and experimental database for modeling the solubility of REE phosphate minerals (monazite, xenotime, and rhabdophane) from 25 to 300 °C

Rare earth elements (REE) are critical elements found in monazite, xenotime, and hydrated REE phosphates which typically form in hydrothermal mineral deposits. Accurate predictions of the solubility of these REE phosphates and the speciation of REE in aqueous fluids are both key to understanding the controls on the transport, fractionation, and deposition of REE in natural systems. Previous monazite and xenotime solubility experiments indicate the presence of large discrepancies between experimentally derived solubility constants versus calculated solubilities by combining different data sources for the thermodynamic properties of minerals and aqueous species at hydrothermal conditions. In this study, these discrepancies were resolved by using the program GEMSFITS to optimize the standard partial molal Gibbs energy of formation (Δ f G° 298 ) of REE aqueous species (REE 3+ and REE hydroxyl complexes) at 298.15 K and 1 bar while keeping the thermodynamic properties fixed for the REE phosphates. A comprehensive experimental database was compiled using solubility data available between 25 and 300 °C. The latter permits conducting thermodynamic parameter optimization of Δ f G° 298 for REE aqueous species. Optimal matching of the rhabdophane solubility data between 25 and 100 °C requires modifying the Δ f G° 298 values of REE 3+ by 1–6 kJ/mol, whereas matching of the monazite solubility data between 100 and 300 °C requires modifying the Δ f G° 298 values of both REE 3+ and REEOH 2+ by ~15–31 kJ/mol and ~2–10 kJ/mol, respectively. For xenotime, adjustments of Δ f G° 298 values by 1–26 kJ/mol are only necessary for the REE 3+ species. The optimizations indicate that the solubility of monazite in acidic solutions is controlled by the light (L)REE 3+ species at <150 °C and the LREEOH 2+ species at >150 °C, whereas the solubility of xenotime is controlled by the heavy (H)REE 3+ species between 25 and 300 °C. Based on the optimization results, we conclude that the revised Helgeson-Kirkham-Flowers equation of state does not reliably predict the thermodynamic properties of REE 3+ , REEOH 2+ , and likely other REE hydroxyl species at hydrothermal conditions. We therefore provide an experimental database (ThermoExp_REE) as a basic framework for future updates, extensions with other ligands, and optimizations as new experimental REE data become available. As a result, the optimized thermodynamic properties of aqueous species and minerals are available open access to accurately predict the solubility of REE phosphates in fluid-rock systems.

58 GEOSCIENCES↗
DOE OSTI2024

Hydrothermal solubility of Dy hydroxide as a function of pH and stability of Dy hydroxyl aqueous complexes from 25 to 250 °C

The rare earth elements (REE) have important applications in green energy technologies. The formation of mineral deposits in geologic systems commonly involves hydrothermal fluids which can mobilize the REE. However, the REE speciation is not well known as a function of pH. The thermodynamic properties of REE hydroxyl complexes used in geochemical models are based on the Helgeson-Kirkham-Flowers (HKF) equation of state parameters which were derived by extrapolation of low temperature experimental and estimated data. In this study, Dy hydroxide solubility experiments are combined with available literature data to improve these models from 25 to 250 °C and optimize the thermodynamic properties of Dy 3+ and Dy hydroxyl complexes using GEMSFITS. Batch-type solubility experiments were conducted from 150 to 250 °C and at saturated water vapor pressure in perchloric acid solutions with initial pH values of 2 to 5 in 0.5 pH unit increments. The measured solubility of Dy hydroxide is retrograde with temperature and decreases with pH. The logarithm of total dissolved Dy molality ranges from –2.3 to –5.3 at 150 °C (pH 4.7–5.5), from –2.4 to –5.6 at 200 °C (pH 3.9–5.1), and from –3.7 to –6.9 at 250 °C (pH of 3.4 and 5.0). The optimized standard partial molal Gibbs energies of formation (Δ f G° T ) derived for Dy 3+ and DyOH 2+ display a close to linear relationship with temperature, fitting with previous optimizations based on DyPO 4 solubility data in the literature. A comparison of the optimized ΔfG°T values for aqueous Dy species with predictions from available HKF parameters indicates significant differences ranging from +11 to –26 kJ/mol between 25 and 250 °C. The experimental fits are used to derive the Dy hydroxide solubility products (K s0 ) and formation constants for the hydrolysis of Dy (β n with n = 1 to 3; Dy 3+ + nOH – = DyOH n 3-n ) as a function of temperature. The optimization method presented yields accurate thermodynamic properties for the Dy 3+ aqua ions and the DyOH 2+ species at the acidic to mildly acidic pH studied whereas more experimental work is needed at near-neutral and alkaline conditions to better constrain the other hydroxyl complexes. Furthermore, the optimized thermodynamic data have a significant impact on geochemical modeling of the mobility and solubility of REE minerals in acidic hydrothermal fluids.

58 GEOSCIENCES↗
DOE OSTI2026

Ab initio thermodynamics of Ni and Co incorporation in Mg hydroxide, carbonate, and hydroxycarbonate minerals

Ni and Co are critical elements needed for modern technologies, and a better understanding of the ability of Mg-based minerals to incorporate these elements would benefit strategy development for Ni and Co recovery from mafic and ultramafic deposits. Here, in this work, we performed density functional theory (DFT) calculations of Ni and Co incorporation in six potential products of the carbonation of mafic and ultramafic silicates: brucite (Mg(OH) 2 ), magnesite (MgCO 3 ), nesquehonite (MgCO 3 ⸱3H 2 O), lansfordite (MgCO 3 ⸱5H 2 O), artinite (Mg 2 CO 3 (OH) 2 ⸱3H 2 O), and hydromagnesite (Mg 5 (CO 3 ) 4 (OH) 2 ⸱4H 2 O). The DFT results were used in an ab initio thermodynamics framework to explore the pH 2 O–pCO 2 conditions at which the Mg-based minerals were predicted to be thermodynamically stable and to quantify the Gibbs free energy of Ni and Co substitution at Mg sites. Among the six Mg-based minerals, brucite and magnesite were predicted to have the lowest Ni and Co substitution free energy. An analysis of the effect of temperature indicated that, at low temperature (<100 K), brucite more readily accommodated Ni and Co, while, at higher temperature (>335 K), magnesite more favorably incorporated Ni and Co. Between 100 K and 335 K, Ni was predicted to preferentially substitute for Mg in brucite and Co for Mg in magnesite, thus leading to a driving force for separating Ni and Co in conditions where brucite and magnesite both form. Insights gained in this work could therefore help select experimental conditions that either promote or inhibit incorporation of these critical elements into Mg-based mineral phases.

Critical elements↗
DOE OSTI2024

Adsorption of terbium (III) on DGA and LN resins: Thermodynamics, isotherms, and kinetics

Two commercially available extraction chromatography (EXC) resins containing N,N,N’,N’-tetra-n-octyldiglycolamide (DGA Resin, Normal, 50 – 100 μm) and Bis(2-ethylhexyl) phosphate (LN Resin, 100 – 150 μm) were used as adsorbents to study fundamental adsorption properties such as thermodynamic values, equilibrium isotherms, and kinetic uptake models for terbium(III) adsorption. Weight distribution ratios (D w ) for terbium on DGA and LN resins were measured using a [ 160 Tb]Tb 3+ radiometric tracer in nitric acid as a function of acidity, temperature, initial analyte concentration, and equilibrium time. The D w values showed increasing binding affinity for DGA resin at high nitric acid concentrations and decreasing binding affinity for LN resins. Thermodynamic studies for DGA and LN resins revealed that the Gibbs free energy (ΔG) increased consistently with temperature. To model equilibrium data, increasingly higher parameter equilibrium isotherm models (Henry (1) < Langmuir, Freundlich (2) < Redlich-Peterson (3) < Fritz-Schluender (4)) were compared on their root mean squared errors (RMSE) and adjusted determination coefficients to determine the most applicable model. In all cases, the empirical four-parameter Fritz-Schluender isotherm demonstrated a superior fit. Similar comparisons for reaction-based kinetic models (Pseudo-first-order < Pseudo-second-order < Pseudo-n-order) revealed that the higher-order PNO model yielded a superior fit of kinetic data for both resins. Furthermore, in some cases, adsorption isotherms and kinetic models could also be modeled by a lower-order model with minimal change in error parameters. Weber-Morris plots revealed that two linear sections are observed for each resin, where the first linear segment is attributed to fast (film diffusion) adsorption of terbium, followed by slower intraparticle diffusion of terbium through the pores as the rate-limiting step. Based on the Weber-Morris plot, both film and intraparticle diffusion are involved in controlling the kinetic rate of adsorption for DGA and LN resins.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

FFT-based solver for higher-order and multi-phase-field fracture models applied to strongly anisotropic brittle materials

This paper presents the application of a fast Fourier transform (FFT) based method to solve two phase field models designed to simulate crack growth of strongly anisotropic materials in the brittle regime. By leveraging the ability of the FFT-based solver to generate solutions with higher-order and global continuities, we design two simple algorithms to capture the complex fracture patterns (e.g. sawtooth, and curved crack growth) common in materials with strongly anisotropic surface energy via the multi-phase-field and high-order phase-field frameworks. A staggered operator-split solver is used where both the balance of linear momentum and the phase field governing equations are formulated in the periodic domain. The unit phase field of the initial failure region is prescribed by the penalty method to alleviate the sharp material contrast between the initial failure region and the base material. The discrete frequency vectors are generalized to estimate the second and fourth order gradients such that the Gibbs effect near shape interfaces or jump conditions can be suppressed. Furthermore, a preconditioner is adopted to improve the convergence rate of the iterative linear solver. Three numerical experiments are used to systematically compare the performance of the FFT-based method in the multi-phase-field and high-order phase-field frameworks.

42 ENGINEERING↗