Search NASA⌕ Search

SEARCH · Search NASA

Results for “Solution thermodynamics”

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 181 records · Page 10

Overcoming time and complexity limitations in molecular dynamics investigations of equilibrium melting

Abstract A hybrid Monte-Carlo molecular-dynamics method for determining solidus and liquidus compositions in multicomponent systems is presented that overcomes both the time limitations in conventional molecular dynamics that prevent the evolution of distinct solid and liquid compositions via diffusion and the complexity challenge that prevents use of thermodynamic assessment in systems of many components. This hybrid method is validated in the Cu–Ni system against an independent assessment of solidus and liquidus compositions based on the regular solution model. Strategies for efficient mapping of different phase diagrams, based on the thermodynamic parameter T 0 , the temperature at which two phases of the composition X 0 have equal free energies, are presented and then demonstrated for the copper-nickel fully miscible system and the gold–silicon eutectic system. A calculation of the solidus and liquidus sampled during the equilibrium melting of equiatomic CrMnFeCoNi is performed, indicating that this method has potential to be extended to the study of many component alloys.

Au-Si↗

Non-thermal atmospheric pressure plasma–liquid synthesis of organic acids in aqueous solution from carbon monoxide

This work aims at understanding the conversion of CO to organic acids, namely oxalic acid and formic acid, using non-thermal atmospheric pressure plasma over aqueous solutions. CO exhibited significantly higher conversion to organic acids (more than 15×) compared to CO 2 under the same reaction conditions. The result bolsters a proposed two-step process for CO 2 fixation, whereby CO 2 is first converted to CO, and then CO is converted to organic acids. The organic acids produced from CO are intermediates in the water–gas shift (WGS) reaction of CO in the presence of an aqueous solution to dissolved CO 2 and hydrogen gas. Based on a simple thermodynamic analysis, the organic acid yield was increased by lowering the plasma–liquid reaction temperature using an ice bath to cool the reaction flask. The composition of the organic acids could be varied by changing the pH of the solution. Oxalate was formed in higher concentrations with increasing solution pH above the pK a of the radical species (CO 2 )˙ − . Below the pK a value, formate was the exclusive organic acid formed. The production of formate has a rather weak pH dependence but is enhanced slightly at a basic pH above 10. Furthermore, at basic pH, the effect of electrolyte concentration comes into play. Higher electrolyte concentrations, leading to shorter electrolyte Debye lengths, resulted in lowered organic acid yields. The highest yields of organic acids obtained in our system were 122 mg L −1 for oxalate and 77 mg L −1 for formate at an optimum 1 mM NaOH concentration in the starting solution. This work is a successful pioneering example of CO to organic acids conversion using non-thermal plasmas, which opens the pathway for a promising two-step conversion process of CO 2 to organic acids.

Sudagar, Alcina Johnson [Washington University in ↗

Thermodynamic assessment of the quaternary WTaCrV refractory high entropy alloy as a means to guide experimental approaches

The deployment of fusion energy poses challenges for materials in plasma facing components to withstand high temperatures and thermal gradients, particle implantation and neutron damage. The current material of choice is tungsten, although property degradation limits its consideration in future fusion reactors. Hence, materials with better resistance to harsh environments need to be developed for fusion energy to become a reality. High entropy alloys are being explored as potential candidates with some compositions showing good radiation resistance to defect cluster formation. One of these materials is the WTaCrV system, although only one composition has been tested under ion irradiation. In this work, we study the thermodynamic properties of the entire quaternary alloy composition range. Coupling first principles calculations, cluster expansion approaches, and Monte Carlo methods, we access the free energy functionals, short-range ordering as a function of temperature, and atomic configurations that can be compared to experimental observations. We use this data to inform experiments into compositions with higher propensity to form solid solutions, instead of phase separating. With this formalism we have developed thermodynamic database (TDB) files that can be used to plot quaternary phase diagrams.

Cluster Expansion↗

Energetics of oxidation and formation of uranium monocarbide

To enable better implementation of uranium monocarbide (UC) as an advanced nuclear fuel for future high-temperature reactors, it is essential to have a thorough knowledge of its thermal and thermodynamic properties under reactor operational conditions. In this work, we studied thermal bulk oxidation of UC by simultaneous thermal analysis consisting of thermogravimetric analysis – differential scanning calorimetry coupled with evolved gas mass spectrometry (TGA-DSC-MS), and we examined the thermodynamic stability of UC using high temperature oxide melt drop solution calorimetry. Further, in air, our studied UC sample (which contains ~5 mol% UO 2 ) was found to undergo a step-wise thermal oxidation process consisting of consecutive oxidations and thermal decomposition reactions: 0.95UC·0.05UO 2 → UO 3 ·0.29(C x O y ) + 0.66CO 2 → UO 3 ·0.20(C x O y ) + 0.09CO 2 → UO 3 ·0.03(C x O y ) + 0.17CO 2 → U 3 O 8 + 0.03CO 2 + 0.166O 2 . DSC was further used to determine the enthalpies of reactions associated with this series of oxidation reactions. Synchrotron X-ray diffraction (XRD) and extended X-ray absorption spectroscopy (EXAFS) were performed to characterize both the long- and short-range structures of UC. The standard enthalpy of formation (ΔH° f ) of UC was determined to be –50.7 ± 10.8 kJ/mol·atom, in good agreement with previous values measured by bomb calorimetry. Lastly, the enthalpic landscape of U-C compounds, including UC, U 2 C 3 , and α-UC 1.94 , were established based on the enthalpy normalized per mole atom, which suggests that U-C phases are thermodynamically stable at lower C/U ratios.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Lattice Boltzmann simulation of the dissolution of slag in alkaline solution using real-shape particles

Highlights: • A dissolution numerical model was proposed to capture the real dissolution kinetics of slag in alkaline solution. • The log forward dissolution rate of Si was described as a function of NBO/T and solution pH. • A threshold solid volume fraction of 0.688 was found for a voxel in 3D, 63.8% larger than that for a pixel in 2D. • The proposed dissolution numerical model provides a reliable alternative to study the dissolution kinetics of slag. A dissolution numerical model was proposed in this study to capture the real dissolution kinetics of slag in alkaline solution. It consists of three modules, i.e. (i) simulation of the initial particle parking structure of slag in alkaline solution using real-shape particles of slag, (ii) simulation of the chemical reactions between slag and solution based on the transition state theory, and (iii) simulation of the physical transport of aqueous ions using the lattice Boltzmann method. This dissolution numerical model was verified using experimental results, showing reasonable accuracy. After verification, the dissolution numerical model was implemented to study the influences of temperature and particle shape using a proper recipe of slag in alkaline solution. This recipe was designed to avoid solid phase precipitation or gel formation via thermodynamic analysis. The simulation results showed faster dissolution kinetics of slag when using higher temperatures and more irregular particle shapes.

36 MATERIALS SCIENCE↗

Atomistic Assessment of Solute-Solute Interactions during Grain Boundary Segregation

Grain boundary solute segregation is becoming increasingly common as a means of stabilizing nanocrystalline alloys. Thermodynamic models for grain boundary segregation have recently revealed the need for spectral information, i.e., the full distribution of environments available at the grain boundary during segregation, in order to capture the essential physics of the problem for complex systems like nanocrystalline materials. However, there has been only one proposed method of extending spectral segregation models beyond the dilute limit, and it is based on simple, fitted parameters that are not atomistically informed. In this work, we present a physically motived atomistic method to measure the full distribution of solute-solute interaction energies at the grain boundaries in a polycrystalline environment. We then cast the results into a simple thermodynamic model, analyze the Al(Mg) system as a case study, and demonstrate strong agreement with physically rigorous hybrid Monte Carlo/molecular statics simulations. This approach provides a means of rapidly measuring key interactions for non-dilute grain boundary segregation for any system with an interatomic potential.

36 MATERIALS SCIENCE↗

Gd-Ni-Sb-SnO 2 electrocatalysts for active and selective ozone production

Direct electrochemical production of dissolved ozone could potentially provide economic wastewater treatment and sanitation or a valuable chemical oxidant. Although Ni-Sb-SnO 2 electrocatalysts have the highest known faradaic efficiencies for electrochemical ozone production, the activity and selectivity are not yet sufficient for commercial implementation. This report finds that co-doping Ni and Gd increases the ozone selectivity by a factor of three over Ni alone. These findings are the first demonstration of an active dopant other than Ni in SnO 2 . Electrochemical and physical characterization show that trends in ozone activity are caused by chemical catalysis, not morphology effects, and that conduction band alignment is not a catalytic descriptor for the system. Selective radical quenching experiments and quantum chemistry calculations of thermodynamic energies suggest that the kinetic barriers to form solution-phase intermediates are important for understanding the role of dopants in electrochemical ozone production.

42 ENGINEERING↗

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↗

Early reactivity of sodium silicate-activated slag pastes and its impact on rheological properties

The fast loss of fluidity of sodium silicate-activated slag (SS-AAS) systems have been widely studied in the literature. In this study, the rheology of these cements has been correlated with the evolution of the reaction products formed over time. For this purpose, a combination of experimental characterization techniques, to study the pore solution composition and microstructure of SS-AAS pastes and thermodynamic calculations have been applied. The initial precipitation of ill-defined N-A-S-H and C-N-A-S-H producta leads to a fast increase of the storage modulus after 40 min of reaction. These products are formed from the reaction of the Na{sup +} ions and silicate ions of the activator with the Ca{sup 2+} and aluminate species dissolved from the slag. However, when a shear is applied, the structure breaks down and fluidity is recovered, that infers that the attractive forces between the particles are low at the early stages of hydration.

36 MATERIALS SCIENCE↗

Comparison of Heterogeneous Langmuirian Models for Mixed-Gas Adsorption Equilibria

Mixed-gas adsorption remains a very active area of research, where newly proposed adsorption equilibrium models need to be compared to standards in the field. The generalized Langmuir (gL) isotherm, proposed by Hamid et al. to address adsorption of different size molecules on heterogeneous adsorbents without requiring evaluation of the Gibbs isotherm, has thus far been tested mainly against the adsorbed solution theory (AST). The gL isotherm was shown to be thermodynamically consistent and gave either better or similar representation than AST models. Here, the gL isotherm is contrasted with the Dual-Process Langmuir isotherm popularized by Ritter and his co-workers and the Loading Ratio Correlation popularized for simple mixtures, both of which are considered the “standard” models for process simulation. Furthermore, the comparison is based on their theoretical roots, interpretations, and representation of the equilibrium loading and the isosteric enthalpy of adsorption of unary, binary, and ternary mixtures of various gaseous adsorbates on different adsorbents over a range of experimental pressures and temperatures.

42 ENGINEERING↗

Calorimetric Study of Mixed Phosphates Na 4 M 3 (PO 4 ) 2 P 2 O 7 (M = Mn 2+ , Fe 2+ , Co 2+ , Ni 2+ ) to Evaluate the Electrochemical Trends

Mixed polyanionic compounds have been studied extensively as viable cathode materials for sodium-ion batteries. Mixed phosphates, Na 4 M 3 (PO 4 ) 2 P 2 O 7 (M = Mn 2+ , Fe 2+ , Co 2+ , Ni 2+ ), provide a low barrier for Na-ion diffusion, being advantageous in comparison to phosphates and pyrophosphates. The reported order of sodium extraction is ambiguous and remains unclear. Despite being structurally similar, electrochemical performance differs for all four analogues with different degrees of (de)sodiation, according to the transition element present. High-temperature oxide melt solution calorimetry has been used to establish the relation between thermodynamic phase stability and observed capacity for this series of mixed phosphates. Thermodynamic phase stability largely depends on the kind of structure, type of bonding, and size of the cations present. So, according to our results, the thermodynamic phase stability follows the order Na 4 Mn 3 (PO 4 ) 2 P 2 O 7 > Na 4 Fe 3 (PO 4 ) 2 P 2 O 7 > Na 4 Co 3 (PO 4 ) 2 P 2 O 7 > Na 4 Ni 3 (PO 4 ) 2 P 2 O 7 . The thermodynamic studies serve as guidelines for the selection of compositions with the potential for fabricating advanced cathode materials with maximum performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Screening of metal complexes and organic solvents using the COSMOSAC-LANL model to enhance the energy density in a non-aqueous redox flow cell: an insight into the solubility

In this paper, we have proposed a first-principles methodology to screen transition metal complexes against a particular organic solvent and organic solvents against a particular transition metal complex based on their solubility information without the knowledge of heat of fusion and melting temperature. The energy density of a non-aqueous redox flow cell directly depends on the solubility of the redox active species in the non-aqueous medium. Here, we have used the “COSMOSAC-LANL” activity coefficient model which is based on first-principles COSMO calculations where the microscopic information is passed to the macroscopic world via a dielectric continuum solvation model, followed by a post-statistical thermodynamic treatment of the self-consistent properties of the solute particle to calculate the solubility. To model the activity coefficient at infinite dilution for the binary mixtures, a 3-suffix Margules (3sM) function is introduced for the quantitative estimation of the asymmetric interactions and, for the combinatorial term, the Staverman–Guggenheim (SG) form is used. The new activity coefficient model is separately called the “LANL” activity coefficient model. The metal complex and the organic solvent have been treated as a simple binary mixture. The present model has been applied to a set of 14 different organic solvents and 16 different transition metal complexes. Using the new LANL activity coefficient model in combination with the ADF-COSMOSAC-2013 model, we have shown how one can improve the solubility of a transition metal complex in an organic solvent. We applied our model to screen 84 binary mixtures to predict the compatible pair of redox active species and organic solvent to increase the energy density. The solvation mechanism of the transition metal complexes in the organic solvents was obtained using the new model. The results have been compared with the experimental and theoretical results where they are available.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Selective formation of metastable polymorphs in solid-state synthesis

Metastable polymorphs often result from the interplay between thermodynamics and kinetics. Despite advances in predictive synthesis for solution-based techniques, there remains a lack of methods to design solid-state reactions targeting metastable materials. Here, we introduce a theoretical framework to predict and control polymorph selectivity in solid-state reactions. This framework presents reaction energy as a rarely used handle for polymorph selection, which influences the role of surface energy in promoting the nucleation of metastable phases. Through in situ characterization and density functional theory calculations on two distinct synthesis pathways targeting LiTiOPO 4 , we demonstrate how precursor selection and its effect on reaction energy can effectively be used to control which polymorph is obtained from solid-state synthesis. A general approach is outlined to quantify the conditions under which metastable polymorphs are experimentally accessible. With comparison to historical data, this approach suggests that using appropriate precursors could enable targeted materials synthesis across diverse chemistries through selective polymorph nucleation.

36 MATERIALS SCIENCE↗

CO2 Capture Strategies via Mineralization with Industrial Waste Brines

Large coal-fired power plants (>500 MW) account for 30% of global CO2 emissions, and long-term management of this CO2 to is urgently needed mitigate global temperature increases. Sequestration of CO2 within stable mineral carbonates (e.g., CaCO3) represents an attractive emission reduction strategy because it offers a leakage-free alternative to geological storage of CO2 in an environmentally friendly form. We have previously described a mineralization process in which divalent cations are sourced from various waste streams (e.g., produced water and brackish water) and alkalinity is induced via regenerable ion-exchange materials (Bustillos et. al. Frontiers in Energy Research. 2020, 8, 352). In our process, aqueous carbonate-bearing streams with pH > 8 are produced by contacting fresh water and carbon dioxide with various ion-exchange materials (e.g., Na form zeolites or ion exchange resins). These streams are mixed with produced water containing varying concentrations (~0.01 – 1.0 M) of Ca2+ leading to the precipitation of solid calcium carbonate (PCC). This process has the advantages of using regenerable solids in a simple and continuous process to increase the pH of water by ion exchange instead of relying on the consumption of costly and unsustainable sources of alkalinity (e.g., sodium hydroxide). While once-through column experiments showed the above benefits, the same were yet to established in a steady-state process with recycle streams. In this work, we set up a process simulation to quantify the energy requirements and CO2 emissions associated with the process and seek optimal produced water compositions and CO2 concentrations (5 – 20 vol%). The process simulation was set up in ASPEN Plus using eRNTL as the thermodynamic property method and sequential modular strategy. Ion exchange alkaline solution was simulated using sodium hydroxide and validated against the experimental data obtained from once-through kinetic experiments. Nanofiltration and reverse osmosis membrane steps were also implemented for the separation of divalent cations and production of fresh water and a regeneration stream following mineralization. Sensitivity analysis was carried out using a range of produced water compositions (0.01 – 1.0 M Ca2+, 0.001 – 0.15 M Mg2+, 0.5 – 3.5 M Na+ and 0.0004 – 0.002 M Fe2+) according to the United States Geological Survey (USGS) database. Calcium carbonate yields increased with increasing CO2 concentrations and were maximized using produced water compositions with larger Ca2+ concentrations. Maximum calcium carbonate yields produced at 5 vol%, 12 vol% and 20 vol% CO2 were 2.3 mmol/L, 5.5 mmol/L, and 9.3 mmol/L, respectively, with the formation of brucite (a magnesium hydroxide phase, Mg(OH)2) and goethite (an iron hydroxide phase, FeOOH) as the primary contaminant phases (99% calcite, 0.6% brucite, 0.4% goethite), which agree with phases detected by XRD experimentally. These results indicate high purity calcium carbonate can be precipitated using industrial waste streams. Consequentially, energy consumption and net CO2 emissions were minimized where precipitated calcium carbonate was maximized for all produced water compositions and CO2 concentrations. Minimum energy consumptions were 0.21 kWh/ton CO2 processed, with 98% of the energy input required coming from the membrane filtration steps. Produced water compositions with large Na+ concentrations (> 0.5 M) were effective at reducing energy consumptions due to faster regeneration time of ion exchange materials. Additionally, calculated net CO2 emissions were negative for the process and ranged from -0.02 kg/ton CO2 to -0.15 kg/ton CO2 processed, indicating a low emission process. We will also present techno-economic assessment showing the economic benefits of the current process as an alternative to the addition of stoichiometric bases to induce alkalinity for the precipitation of CaCO3.

Simonetti, Dante↗

Enhancing the Performance of Plasma-facing Materials Through Solute-stabilized Nanostructured Tungsten Alloys

This report is a summary of the work completed under the auspices of the Department of Energy Early Career Award DE-SC0017899 entitled “Enhancing the Performance of Plasma-facing Materials Through Solute-stabilized Nanostructured Tungsten Alloys”. Summaries of initiatives on computational thermodynamics modeling, experimental alloy synthesis including traditional and additive manufacturing routes, several ion irradiation studies, and plasma exposures are reported.

36 MATERIALS SCIENCE↗

Potential-pH diagrams considering complex oxide solution phases for understanding aqueous corrosion of multi-principal element alloys

The potential-pH diagram, a graphical representation of the thermodynamically predominant reaction products in aqueous corrosion, is originally proposed for the corrosion of pure metals. The original approach only leads to stoichiometric oxides and hydroxides as the oxidation products. However, numerous experiments show that non-stoichiometric oxide scales are prevalent in the aqueous corrosion of alloys. In the present study, a room temperature potential-pH diagram considering oxide solid solutions, as a generalization of the traditional potential-pH diagram with stoichiometric oxides, is constructed for an FCC single-phase multi-principal element alloy (MPEA) based on the CALculation of PHAse Diagram method. The predominant reaction products, the ions in aqueous solution, and the cation distribution in oxides are predicted. The oxide solid solution is stabilized by the mixing free energy (or mixing entropy) and the stabilizing effect becomes more significant as the temperature increases. Consequently, solid solution oxides are stable in large regions of the potential-pH diagram and the mixing free energy mostly affects the equilibrium composition of the stable oxides, while the shape of stable regions for oxides is mostly determined by the structure of the stable oxides. Agreements are found for Ni 2+ , Fe 2+ , and Mn 2+ between the atomic emission spectroelectrochemistry measurements and thermodynamic calculations, while deviations exist for Cr 3+ and Co 2+ possibly due to surface complexation with species such as Cl - and the oxide dissolution. By incorporating the solution models of oxides, the current work presents a general and more accurate way to analyze the reaction products during aqueous corrosion of MPEAs.

Materials Science↗

Computational thermomechanics for crystalline rock. Part II: Chemo-damage-plasticity and healing in strongly anisotropic polycrystals

We present a thermal–mechanical–chemical-phase field model that captures the multi-physical coupling effects of precipitation creeping, crystal plasticity, anisotropic fracture, and crack healing in polycrystalline rock at various temperature and strain-rate regimes. This model is solved via a fast Fourier transfer solver with an operator-split algorithm to update displacement, temperature and phase field, and chemical concentration incrementally. In nuclear waste disposal in salt formation, brine inside the crystal salt may migrate along the grain boundary and cracks due to the gradient of interfacial energy and pressure. This migration has a significant implication on the permeability evolution, creep deformation, and crack healing within rock salt but is difficult to incorporate implicitly via effective medium theories compared with computational homogenization. As such, we introduce a thermodynamic framework and a corresponding computational implementation that explicitly captures the brine diffusion along the grain boundary and crack at the grain scale. Meanwhile, the anisotropic fracture and healing are captured via a high-order phase field that represents the regularized crack region in which a newly derived non-monotonic driving force is used to capture the fracture and healing due to the solution–precipitation. Numerical examples are presented to demonstrate the capacity of the thermodynamic framework to capture the multiphysics material behaviors of rock salt.

42 ENGINEERING↗

Development and application of ReaxFF methodology for understanding the chemical dynamics of metal carbonates in aqueous solutions

A new ReaxFF reactive force field has been developed for metal carbonate systems including Na + , Ca 2+ , and Mg 2+ cations and the CO 3 2– anion. This force field is fully transferable with previous ReaxFF water and water/electrolyte descriptions. The Me–O–C (Me = metal) three-body valence angle parameters and Me–C non-reactive parameters of the force field have been optimized against quantum mechanical calculations including equations of state, heats of formation, heats of reaction, angle distortions and vibrational frequencies. The new metal carbonate force field has been validated using molecular dynamics simulations to study the solvation and reactivity of metal and carbonate ions in water at 300 K and 700 K. The coordination radius and self-diffusion coefficient show good consistency with existing experimental and simulation results. The angular distribution analysis explains the structural preference of carbonate ions to form carbonates and bicarbonates, where Na + predominantly forms carbonates due to weaker angular strain, while Ca 2+ and Mg 2+ prefer to form bicarbonate monodentate in nature. Residence time distribution analyses on different systems reveal the role of ions in accelerating and decelerating the dynamics of water and carbonate ions under different thermodynamic conditions. The formation and dissolution of bicarbonates and carbonates in solution were explored on the basis of the protonation capability in different systems. Here, the nucleation phenomenon of metal carbonates at ambient and supercritical conditions is explained from the perspective of cluster formation over time: Ca 2+ ions can form prenucleation clusters at ambient temperature but show saturation with increasing temperature, whereas Na + and Mg 2+ ions show a rapid increase in cluster size and amount upon increasing time and temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗