Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CALPHAD modeling”

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

High-Temperature Thermodynamics Modeling of Graphite

We present high-temperature thermodynamic properties for graphite from first-principles anharmonic theory. The ab initio electronic structure is obtained from density-functional theory coupled to a lattice dynamics method that is used to model anharmonic lattice vibrations. This combined approach produces free energies and specific heats for graphite that compare well with available experiments and results from models that empirically represent experimental data, such as CALPHAD. We show that anharmonic theory for the phonons is essential for accurate thermodynamic quantities above about 1000 K.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Development of Innovative Non-Concentrated Alloy (NCA) Fuel Cladding for Advanced Nuclear Reactors

As the nuclear energy sector advances toward next-generation reactors, the need for high-performance fuel cladding materials has become increasingly urgent. Traditional alloys like zirconium and stainless steel are reaching their performance limits under higher temperatures, more corrosive coolants, and extended irradiation. This report presents the development of a new class of fuel cladding materials based on Non-Concentrated Alloys (NCAs) ? multi-element systems designed to deliver enhanced mechanical strength, corrosion resistance, and radiation tolerance. Through a combination of computational modeling (Computer Coupling of Phase Diagrams and Thermochemistry: CALPHAD), simulation-guided alloy selection, and experimental fabrication via arc melting and spark plasma sintering (SPS), three strategic alloy design paths were explored: (1) FeCrAl-based NCAs, (2) refractory-lean neutron-efficient alloys, and (3) equimolar high-entropy compositions. Microstructural analysis confirmed the formation of stable body-centered cubic BCC_A2 phases, while mechanical testing demonstrated hardness values significantly exceeding those of conventional cladding materials. The results highlight the tunability of NCA systems and their potential for balancing strength and ductility ? a critical consideration for in-reactor performance. Looking forward, future work will focus on thermomechanical optimization, CALPHAD refinement, and benchmarking against industry standards to enable scalable deployment. This work not only advances the science of nuclear materials but also supports broader goals in nuclear safety, performance, and nuclear energy innovation.

36 - MATERIALS SCIENCE↗
DOE OSTI2023

Predicting Sintering Window of Binder Jet Additively Manufactured Parts Using a Coupled Data Analytics and CALPHAD Approach

Batch-to-batch variation in powder compositions for binder jet additive manufacturing (BJAM) can significantly deter defining an “ideal” sintering window for a given alloy. One way to overcome the problem is by running sintering experiments at various temperatures for each batch of the powder. However, such an approach increases the time required to achieve large-scale production of parts. The predictive capabilities of computational thermodynamic tools like CALPHAD can be leveraged to overcome the challenge, especially for binder jet additive manufacturing, since the process occurs under near-equilibrium conditions. However, calculating the sintering window using CALPHAD can be computationally expensive, considering many possible feedstock compositions within “specification”. Here, we generate high throughput CALPHAD data for nickel-based superalloys to develop machine learning models to predict the sintering window rapidly. The predictive capability of the models has been validated using published results on BJAM of Inconel 718 and 625. Further, validated models are lightweight and can be deployed in an industrial setting to get sintering window in an accelerated manner.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Computational and Physics-Based Modeling for the Development of in-Space Welding Technology

Welding will be an enabling technology for in-space assembly and manufacturing (ISAM) and in-space repair operations in support of the growing Space Economy. Prospective use cases of welding include metal repair on long-duration lunar surface missions and assembly of large structures that exceed dimensions of the currently available launch vehicle envelopes such as space arrays, tall lunar transmission towers, or even large transit vehicles powered by nuclear propulsion. There is a brief history of beam welding for in-space use. Soviet and NASA astronauts performed electron beam welds (EBW) in space in the late-1960s through mid-1980s and several parabolic flights of laser beam welding (LBW) were undertaken in the USA and Canada in the early-1990s. Work has since stagnated until the recent boom of the in-space economy has demonstrated need for continued development of EBW and LBW processes, which are desirable for their compactness, low mass, low power consumption, and high wall plug efficiency. No space-based manufacturing laboratory capable of demonstrating, developing, and applying welding processes exists, and costs are high to develop and fly individual space-based trials. This drives the current approach to develop welding technology with computational and physics-based models in tandem with earth-based experiments. Several approaches are taken to complement and drive the selection of instrumentation for planned parabolic and suborbital flights. First, a thermal model of Skylab EBW experiments is established and compared to empirical data. Next, the model is used to show variation between Earth-based and in-space weld structure and properties through exploration of key variables that define the Low earth orbit and lunar surface space environments, namely temperature (40 – 400 K), pressure (10-10 to 104 Pa), and gravity (10-6 g to 1 g). Finally, physical processes associated with heat and mass transport during welding are evaluated through several case studies across the above variable ranges demonstrating need for further research.

in-space welding↗
DOE OSTI2024

An integrated modeling framework with open architecture for phase field simulation of multi-component alloys

An integrated modeling framework (PanPhaseField) has been developed, which enables a direct and fast coupling between CALPHAD calculations and large-scale phase field simulations for multi-component alloys. Further, it adopts an open architecture allowing for integration of user-defined phase field models in a plug-and-play manner by taking full advantage of the user-friendly graphical interface of Pandat software. The developed modeling platform becomes an enabling tool that can be used to simulate the evolution of spatially varying microstructures of industrial complex alloys for various engineering applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Thermodynamic description of molten salt systems: KCl-LiCl-NaCl and KCl-LiCl-NdCl 3

The mixture of KCl and LiCl has been used as electrolyte in the electrorefining process to recover uranium from used nuclear fuels due to the low melting point. However, lanthanides and sodium in the reactor waste continuously dissolve into it and thus alter its thermodynamic properties. To understand the thermodynamic behavior evolution of the electrolyte with the accumulation of impurities, thermodynamic modeling for KCl-LiCl-NaCl and KCl-LiCl-NdCl 3 and four constituent binary systems in the entire composition space was performed using the CALPHAD (CALculation of PHAse Diagrams) approach. The ionic liquid was described by the two-sublattice model, where neutral species were introduced to consider short-range ordering (SRO) within the melt, whereas the solid solution was modeled based on the Compound Energy Formalism. Literature data on phase equilibria and thermochemical properties were critically evaluated and used during the optimization of thermodynamic parameters for KCl-LiCl-NaCl and KCl-LiCl-NdCl 3 and their subsystems. The calculated phase diagrams and mixing enthalpies are in good agreement with the experimental data. The thermodynamic modeling for the KCl-LiCl-NdCl 3 system was carried out for the first time. To fill the gap in experimental measurement, enthalpy of mixing for the KCl-LiCl-NdCl 3 melt was estimated using the surrounded-ion model. These data then served as critical inputs for thermodynamic optimization. Furthermore, the present study can provide insights into thermodynamic property evolution of the electrolyte and solubility limit of various impurities during the electrorefining process.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

CALPHAD Uncertainty Quantification and TDBX

CALPHAD uncertainty quantification (UQ) is the foundation of materials design with quantified confidence. We report a framework and software packages to enable CALPHAD UQ assessment and calculation using commercial CALPHAD software (Thermo-Calc). This Bayesian inference framework is coupled with a Markov chain Monte Carlo algorithm to establish uncertainty traces with a given thermodynamic database file (TDB) and corresponding experimental data points. This general framework is demonstrated with the Ni–Cr binary system. The algorithm is firstly validated on synthetic data with known ground truth. Then it is applied to real experimental data to generate posterior traces. We develop a file format named TDBX, which provides a single source of truth by combining the original TDB content and the traces for each assessed Gibbs energy parameter. CALPHAD UQ calculations are performed based on the TDBX file, from which uncertainties for phase boundaries, enthalpy curves, and solidification range are collected as examples of basic design parameters. This TDBX file with corresponding scripts are made open-source. Finally, the combination of CALPHAD UQ assessments and calculations connected by TDBX supports uncertainty-assisted modeling, enabling the integrated application of modern design with uncertainty methodologies to computational materials design.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thermodynamics of Plutonium Monocarbide from Anharmonic and Relativistic Theory

Thermodynamics of plutonium monocarbide is studied from first-principles theory that includes relativistic electronic structure and anharmonic lattice vibrations. Density-functional theory (DFT) is expanded to include orbital-orbital coupling in addition to the relativistic spin-orbit interaction for the electronic structure and it is combined with anharmonic, temperature dependent, lattice dynamics derived from the self-consistent ab initio lattice dynamics (SCAILD) method. The obtained thermodynamics are compared to results from simpler quasi-harmonic theory and experimental data. Formation enthalpy, specific heat, and Gibbs energy calculated from the anharmonic model are validated by direct comparison with a calculation of phase diagram (CALPHAD) assessment of PuC and sub-stochiometric PuC0.896. Overall, the theory reproduces CALPHAD results and measured data for PuC rather well, but the comparison is hampered by the sub-stoichiometric nature of plutonium monocarbide. It was found that a bare theoretical approach that ignores spin-orbit and orbital-orbital coupling (orbital polarization) of the plutonium 5f electrons promotes too soft phonons and Gibbs energies that are incompatible with that of the CALPHAD assessment of the experimental data. The investigation of PuC suggests that the electronic structure is accurately described by plutonium 5f electrons as “band like” and delocalized, but correlate through spin polarization, orbital polarization, and spin-orbit coupling, in analogy to previous findings for plutonium metal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRSDate not supplied

Open-source Numerical Modeling of Solidification Cracking Susceptibility: Application to Refractory Alloy Systems

Introduction. Alloys such as aluminum, nickel-base, and austenitic stainless steels are susceptible to solidification cracking during welding and 3D printing. Compositional optimization is one method used to effectively mitigate solidification cracking of those alloy systems. With the surge in hypersonic and in-space propulsion activities, refractory metals (Nb, Mo, Ta, W, and Re) and their alloy derivatives are increasing in importance due to their extreme high melting point and retention of high-temperature strength; however, their chemistry was most typically optimized to promote ductility during mechanical operations such as drawing and forming. Welding of such alloys has been a challenge due to a number of issues including solidification cracking, atmospheric contamination (O, C, and N), as well as a shift in ductile-to-brittle transition to higher temperature following grain growth induced by welding. Compositional optimization of refractory alloys for solidification cracking resistance in particular is desirable as their usage increases with the advent of advanced manufacturing methods such as 3D printing. This work evaluates the effect of compositional variation in refractory metal systems on the solidification cracking susceptibility with the goals of optimizing existing alloys and joining process techniques, and formulating new alloys with increased solidification cracking resistance. Experimental Procedures. A python code was developed in a Jupyter notebook environment (Michael and Sowards, 2023) to facilitate the calculation of crack susceptibility index proposed by Kou (2015). Composition is entered as a single point, or as a 1-D or 2-D array. The notebook calls pycalphad (Otis and Liu, 2017 and Bocklund et al, 2020) to calculate the evolution of fraction solid as a function of temperature (under either Scheil or equilibrium assumptions) and then evaluates steepness of the fraction solid curve near the terminal stage of solidification to predict solidification cracking resistance. Open source thermodynamic databases available at online repositories are used (van de Walle). The process is setup in an automated fashion to generate plots that show variation in solidification cracking susceptibility according to composition on 1-D line plots or 2-D contour plots. The Jupyter notebook and crack susceptibility algorithm was also integrated with a widely used commercial CALPHAD code for validation and alloy exploration. Results and Discussion. The crack susceptibility model was first validated against a series of refractory alloy compositions evaluated in past work which utilized a specialized Varestraint test built inside a vacuum chamber environment (Lessman and Gold, 1971). The alloys tested in the Varestraint apparatus included T-111 (Ta-8W-2Hf), ASTAR-811C (Ta-8W-1Re-0.7Hf-0.025C), FS-85 (Nb-27Ta-10W-1Zr), T-222 (Ta-9.6W-2.4Hf-0.01C), Ta-10W, B-66 (Nb-5Mo-5V-1Zr), and SCb-291 (Nb-10W-10Ta). The initial test of the model showed a strong correlation with empirical Varestraint data, i.e., a Spearman rank correlation between model predictions and hot cracking measurements was observed to be greater than 0.8. Following the validation, a set of refractory metal binary mixtures was investigated to evaluate sensitivity of Nb, Mo, W, and Ta to C, N, and O content. A series of plots were produced that suggest ppmw ranges of C, N, and O where solidification cracking increases significantly and reaches a maximum. Also comparative ranking of each primary refractory metal to each interstitial was produced. For example C produces greater cracking response in Mo whereas O produces greater cracking response in Ta and Nb. Such compositional values have utility in setting limits on pickup of these interstitial elements during welding and printing rather than using a one-size-fits-all approach. Furthermore, the results have use in determining additive powder recycling requirements, which is especially pertinent for refractory metal powders due to their high cost compared to conventional alloys. Another application created thousands of hypothetical alloys within the nominal specified composition range of two widely used refractory alloys C103 (Nb-10Hf-1Ti) and TZM (Mo-0.5Ti-0.1Zr). The cracking index was calculated for the alloys and results were fed into machine learning regression techniques including Multiple Linear Regression, Ridge Regression, and Lasso Regression to determine relative potency each alloying element had on computed solidification cracking index. A series of linear equations were produced that relate composition of C103 and TZM to solidification cracking index. The crack susceptibility of C103 for example is described by an equation of the form: cracking index ~ O + 0.667*C + 0.635*N + 0.00037*Ta – 0.0008*Hf (in wt.%) From that equation, it is clear that O has strong propensity to induce solidification cracking. Interestingly, Hf is shown to reduce calculated cracking response. Finally, realizing the potential of this method to discover new refractory alloy formulations across the period table that have low solidification cracking sensitivity, the code was applied to new untested alloy systems including W-Zr-C, W-Ta-C, and others. Conclusions. In summary, an open source numerical method has been developed using Python code to calculate Kou’s crack susceptibility index. The method was applied to refractory metals which are inherently difficult to study from a weldability testing standpoint since inert shielding gas is not sufficient and welding is typically done in vacuum, especially in light of findings presented here where oxygen has profound influence on solidification cracking. This work revealed the effect of compositional variations on a series of refractory metals and showed the framework defined here will be useful in 1) the development of new alloys that have improved weldability and 3D printability, 2) placing compositional limits on existing alloys, and 3) ensuring adequate controls of manufacturing processes such as 3D printing where powder reuse is critical. Keywords. pycalphad; Python; refractory metals; solidification cracking. References. B. Bocklund et. al. (2020) http://doi.org/10.5281/zenodo.3630657. S. Kou. (2015) https://doi.org/10.1016/j.actamat.2015.01.034. G.G. Lessmann and R.E. Gold. Welding Journal, issue 1, pp. 1-s – 8-s (1971). F.N. Michael and J.W. Sowards. NASA/TM-20230002218 (2023). R. Otis and Z.-K. Liu. (2017) http://doi.org/10.5334/jors.140. A. Van de Wallle et. al. (2018) https://doi.org/10.1016/j.calphad.2018.04.003.

pycalphad↗
DOE OSTI2022

Oxidation Lifetime Modeling of 625 and 120 Foils After Long-term Exposure in Flowing Air + 10% ${\hbox {H}_2\hbox {O}}$ at 700 and 800 $^\circ$C

A model describing oxidation kinetics considering simultaneous oxidation and volatilization of external oxide scales was integrated with a CALPHAD-based coupled thermodynamic and kinetic approach to predict the oxidation-induced lifetime of foils of alloys 625 and 120. Long-term exposures in flowing air +10% H 2 O at 700 and 800 °C under two flow velocities were used to validate the modeling results. The time to a critical Cr concentration of 10 wt% at the oxide/alloy interface, t 10 , was shown to correlate with the breakdown of the Cr 2 O 3 scale for the studied experimental conditions. Finally, the modeling approach was extended to consider time-dependent oxidation and volatilization rate constants to develop a more appropriate lifetime criterion for microturbine recuperator foils.

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

Thermodynamics Modeling for Actinide Monocarbides and Mononitrides from First Principles

The high-temperature thermodynamical properties for the actinide monocarbides and mononitrides ThC, ThN, UC, UN, PuC, and PuN are calculated from first-principles electronic-structure theory. The electronic structure is modeled with density-functional theory (DFT) and is fully relativistic, including the spin-orbit interaction. Furthermore, the DFT is extended to account for orbital–orbital interactions, by means of a parameter-free orbital-polarization (OP) technique, that has proven to be essential for the 5f electrons in plutonium. Strong anharmonicity and the temperature dependence of the lattice vibrations are captured with the self-consistent ab initio lattice dynamics (SCAILD) method. The calculated free energies and heat capacities are compared to published results from quasi-harmonic (QH) theory, and experiments, where available. For the uranium and plutonium compounds, we make use of CALPHAD assessments to help evaluate the theory. Generally, our anharmonic relativistic approach compares well with both CALPHAD and experiments. For the thorium compounds, our theory is in good accord with QH modeling of the free energy at lower temperatures but for the heat capacity the comparison is less favorable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Additively manufactured Ni-20Cr to V functionally graded material: Computational predictions and experimental verification of phase formations

Throughout this work, a database for the Cr-Ni-V system was constructed by modeling the binary Cr-V and ternary Cr-Ni-V systems using the CALPHAD approach aided by density functional theory (DFT)-based first-principles calculations and ab initio molecular dynamics (AIMD) simulations. To validate this new database, a functionally graded material (FGM) using Ni-20Cr and elemental V was fabricated using directed energy deposition additive manufacturing (DED AM) and experimentally characterized. The deposited Ni-20Cr was pure fcc phase, while increasing the amount of V across the gradient resulted in the formation of sigma phase, followed by the bcc phase. The experimentally measured phase data was compared with computational predictions made using a Cr-Ni-V thermodynamic database from the literature as well as the database developed in the present work. The newly developed database was shown to better predict the experimentally observed phases due to its accurate modeling of binary systems within the database and the ternary liquid phase, which is critical for accurate Scheil calculations.

36 MATERIALS SCIENCE↗
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

Simulations of self- and Xe diffusivity in uranium mononitride including chemistry and irradiation effects

A combination of density functional theory and empirical potential atomic scale simulations have been used to determine a model for defect stability and mobility in uranium mononitride (UN), as a function of temperature (T) and N 2 partial pressure (p N 2 ). Using the model, predictions of hypo-stoichiometry under U-rich conditions compare favorably to CALPHAD calculations using the TAF-ID database. Furthermore, our predictions of U and N self-diffusivity are in good agreement with experiments carried out as a function of T at specific partial pressures under thermal equilibrium. The validated atomic scale data have then been implemented within a cluster dynamics method to simulate irradiation-enhanced defect concentrations. All defects and clusters studied have significantly enhanced concentrations, with respect to thermal equilibrium, as T is lowered. The irradiation-enhanced Xe diffusivity is compared to post-irradiation annealing and in-pile experiments. In conclusion, the contributions of various defects and clusters to non-stoichiometry, self-diffusivity, and Xe diffusivity are discussed.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

CALPHAD-based ICME design of single-step aging to enhance mechanical strength of WAAM Haynes 282

To match the strength of wire-arc additive manufactured Haynes 282 to its wrought counterpart via a single-step aging heat treatment, the CALPHAD (Calculation of Phase Diagrams) method is integrated with physics-based process-structure-property models and experimental validation. The integrated computational materials engineering (ICME) framework simulates the effects of aging on γ′ and M 23 C 6 precipitation and the resulting yield strength. To improve simulation reliability, the interfacial energies between γ/γ′ and γ/M 23 C 6 carbides were estimated by comparison with precipitation kinetic modeling and measured precipitate sizes. γ′ and M23C6 were found to precipitate simultaneously between 640 and 860 °C, producing microstructures similar to those produced by two-step aging. The optimal γ′ size for peak yield stress was calculated to be 20–23 nm. WAAM Haynes 282 aged at 780 °C for 50 h exceeded the mechanical performance of its wrought counterpart subjected to two-step aging, though desired properties can also be achieved at 800 °C for 16 h or less. The error in yield strength is less than 20 MPa, demonstrating good agreement between the modeling framework and experiments. Creep studies showed that WAAM Haynes 282 exceeded the calculated rupture time, reaching 481 h. This proposed methodology can accelerate the design of aging heat treatments for any γ′-strengthened nickel-base alloy, minimizing the resources required for trial-and-error experiments.

CALPHAD↗
DOE OSTI2026

An Integrated Computational Materials Engineering (ICME) Approach to Design Nonlinear Transition Zones Between Dissimilar Metals

Current approaches to designing graded transition joints (GTJs) between dissimilar metals often rely on linear changes in both composition profiles and thickness of each sublayer. This increases fabrication cost and may not be optimal with respect to residual stress or the formation of undesirable phases. Here, in this study, GTJs between P91 ferritic/martensitic steel and 347H austenitic stainless steel were designed using Integrated Computational Materials Engineering (ICME) principles with nonlinear composition and length profiles. Guided by inputs from classical mechanics and CALPHAD predictions of carbon chemical potential, a novel transition zone consisting of five discrete compositions was proposed, with the thickness of each sublayer varying according to a brachistochrone-inspired distribution. In addition to carbon potential gradients, CALPHAD was used to predict coefficients of thermal expansion, which were incorporated into finite element models to evaluate stress evolution. The proposed nonlinear design resulted in a smoother carbon potential gradient, lower carbon depletion at the P91 interface, and a comparable residual stress under long-term thermal exposure, compared to a conventional linear design using ten sublayers with equal thickness. This work introduces a brachistochrone-inspired distribution for GTJ design, offering a general framework for optimizing graded interfaces between dissimilar metals.

Directed Energy Deposition↗
DOE OSTI2024

Boron and lithium aqueous thermochemistry to model crud deposition in pressurized water reactors

Here, a comprehensive database of boron and lithium aqueous thermochemistry has been developed for use at elevated temperatures. The Helgeson-Kirkham-Flowers (HKF) formalism provides a framework to describe thermodynamic properties over a broad range of temperatures and pressures. Accuracy at high temperatures is vital to modeling nickel oxide and ferrite fuel deposits that occur in pressurized water nuclear reactors (PWRs). CALPHAD (CALculation of PHAse Diagrams) calculations are performed at PWR crud conditions to predict the stability regions of the solid lithium metaborate (LiBO 2 ) and lithium tetraborate (Li 2 B 4 O 7 ) precipitates. In addition, similar calculations are performed using sodium and potassium instead of lithium in order to assess the thermodynamics associated with preventing the formation of such precipitates. This database contributes to understanding of crud formation and composition and will aid in the prediction of phenomena such as Crud-Induced Power Shifts (CIPS).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗