Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “experimental data”

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

Parked aeroelastic field rotor response for a 20% scaled demonstrator of a 13‐MW downwind turbine

Abstract Aeroelastic parked testing of a unique downwind two‐bladed subscale rotor was completed to characterize the response of an extreme‐scale 13‐MW turbine in high‐wind parked conditions. A 20% geometric scaling was used resulting in scaled 20‐m‐long blades, whose structural and stiffness properties were designed using aeroelastic scaling to replicate the nondimensional structural aeroelastic deflections and dynamics that would occur for a lightweight, downwind 13‐MW rotor. The subscale rotor was mounted and field tested on the two‐bladed Controls Advanced Research Turbine (CART2) at the National Renewable Energy Laboratory's Flatiron Campus (NREL FC). The parked testing of these highly flexible blades included both pitch‐to‐run and pitch‐to‐feather configurations with the blades in the horizontal braked orientation. The collected experimental data includes the unsteady flapwise root bending moments and tip deflections as a function of inflow wind conditions. The bending moments are based on strain gauges located in the root section, whereas the tip deflections are captured by a video camera on the hub of the turbine pointed toward the tip of the blade. The experimental results are compared against computational predictions generated by FAST, a wind turbine simulation software, for the subscale and full‐scale models with consistent unsteady wind fields. FAST reasonably predicted the bending moments and deflections of the experimental data in terms of both the mean and standard deviations. These results demonstrate the efficacy of the first such aeroelastically scaled turbine test and demonstrate that a highly flexible lightweight downwind coned rotor can be designed to withstand extreme loads in parked conditions.

17 WIND ENERGY↗
DOE OSTI2020

Force Balance Model Assessment for Mechanistic Prediction of Sliding Bubble Velocity in Vertical Subcooled Boiling Flow

The bubble sliding after departing from nucleation site is frequently observed in flow boiling systems and, the crucial impact on wall heat transfer has been evidenced through many experiments. As a result, the heat transfer modeling associated with sliding bubble has become one of the subjects of great attention in CFD boiling heat transfer community. The modeling efforts are primarily aimed at improving the existing Heat Flux Partitioning (HFP) model via the implementation of sliding bubble-induced heat transfer. The performance of HFP model depends inherently on the fidelity of sub-models used to predict the fundamental bubble parameters (e.g., bubble departure/lift-off diameter). In the same context, the accurate prediction of sliding bubble parameters (e.g., sliding bubble growth, sliding bubble velocity) is essential to achieving the successful heat transfer modeling associated with sliding bubble. Of many sliding bubble parameters, this paper deals with the sliding bubble velocity. Specifically, the force balance model was assessed in view of the sliding bubble velocity prediction. The parametric effect of key sub-models (e.g., drag force, bubble growth models) used in the force balance equation was investigated. The experimental data from Maity (2000) and Yoo et al. (2016) were used for demonstrating the model predictive performance. It was found that the force balance model proposed in this study was able to predict well the bubble sliding velocity based on the accurate prediction of bubble growth during sliding. The predictive performance was proven with the experimental data measured under various subcooled boiling conditions of both water and refrigerant (NOVEC-7000) flowing upward in vertical channels.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Application of an equation‐oriented framework to formulate and estimate parameters of chemical looping reaction models

Abstract Accurate, predictive reaction models are critical for the design and optimization of chemical looping combustion (CLC) reactors. The formulation and estimation of kinetic parameters for these reaction models using a first‐principles equation‐oriented (EO) approach is particularly beneficial as large amounts of experimental data spanning process‐relevant conditions can be used to estimate parameters in a computationally tractable way. This work demonstrates the application of a novel EO framework to develop reduction reaction kinetic models of an iron‐based CLC oxygen carrier (OC). An optimization problem is formulated to estimate kinetic parameters that provide the best fit to the experimental data. The model predicts the state of the OC with mean square error values of 2.5%–4.4% across the full range of validation data, including multiple reduction cycles.

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

Temperature-dependent phonon lifetimes and thermal conductivity of silicon by inelastic neutron scattering and ab initio calculations

In this paper, inelastic neutron scattering on a single crystal of silicon was performed at temperatures from 100 to 1500 K. These experimental data were reduced to obtain phonon spectral intensity at all wave vectors $\overset{→}{Q}$ and frequencies ω in the first Brillouin zone. Thermal broadenings of the phonon peaks were obtained by fitting and by calculating with an iterative ab initio method that uses thermal atom displacements on an ensemble of superlattices. Agreement between the calculated and experimental broadenings was good, with possible discrepancies at the highest temperatures. Distributions of phonon widths versus phonon energy had similar shapes for computation and experiment. These distributions grew with temperature but maintained similar shapes. Parameters from the ab initio calculations were used to obtain the thermal conductivity from the Boltzmann transport equation, which was in good agreement with experimental data. Despite the high group velocities of longitudinal acoustic phonons, their shorter lifetimes reduced their contribution to the thermal conductivity, which was dominated by transverse acoustic modes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Global tuning of hadronic interaction models with accelerator-based and astroparticle data

In high-energy and astroparticle physics, event generators play an essential role, even in the simplest data analyses. As analysis techniques become more sophisticated, e.g. based on deep neural networks, their correct description of the observed event characteristics becomes even more important. Physical processes occurring in hadronic collisions are simulated within a Monte Carlo framework. A major challenge is the modeling of hadron dynamics at low momentum transfer, which includes the initial and final phases of every hadronic collision. QCD-inspired phenomenological models used for these phases cannot guarantee completeness or correctness over the full phase space. These models usually include parameters which must be tuned to suitable experimental data. Until now, event generators have been developed and tuned mainly on the basis of data from high-energy physics experiments at accelerators. The wealth of data available from the latest generation of astroparticle experiments has not yet been fully exploited, and in many cases is not satisfactorily described. Both kinds of data sets are complementary as astroparticle experiments provide sensitivity especially to hadrons produced nearly parallel to the collision axis and cover center-of-mass energies up to several hundred TeV, well beyond those reached at colliders so far. In this report, we provide an overview of state-of-the-art event generators and their tuning, including the most relevant inputs from high-energy accelerator and astroparticle experiments. We present a road map that shows, for the first time, how the unified tuning of event generators with accelerator-based and astroparticle data can be performed.

Albrecht, J. [Ruhr U., Bochum, RAPP Ctr.; Ruhr U.,↗
DOE OSTI2020

Co-ion specific effect on sodium halides sorption and transport in a cross-linked poly(p-styrene sulfonate-co-divinylbenzene) for membrane applications

Here, sodium halide sorption and transport in a cation exchange membrane based on cross-linked poly(p-styrene sulfonate-co-divinylbenzene) (CR61) were investigated experimentally and theoretically. At fixed co-ion valence, sodium halide partitioning between aqueous electrolyte solutions and ion exchange polymers is only minimally sensitive to co-ion chemical identity. NaCl, NaF and NaI exhibit similar solubility coefficients in CR61, which agrees reasonably well with predictions of the Donnan-Manning model, especially at high external salt concentrations. Ion activity coefficients in the membrane were evaluated and ion sorption thermodynamics were investigated. Ion diffusion coefficients in CR61 were predicted using a recently developed framework based on the Manning diffusion model, which accounts for electrostatic effects on ion diffusion, and the Mackie-Meares model, which accounts for tortuosity/area effects on ion diffusion in swollen polymers. Model predictions agree reasonably well with the experimental data, especially at high external salt concentrations. The observed decrease in salt diffusion coefficients in CR61 relative to values in aqueous solutions is due primarily to tortuosity/area effects. Salt diffusion coefficients correlate with co-ion hydrated radius and membrane water content. Finally, salt permeability coefficients predicted using the solution-diffusion model exhibit good agreement with experimental data, especially at external salt concentration greater than 0.1 mol/L.

42 ENGINEERING↗
DOE OSTI2025

Direct local parametrization of nuclear state densities using the back-shifted Bethe formula

Level densities are often parametrized using the back-shifted Bethe formula (BBF) for nuclei that possess experimental data for s-wave neutron resonance average spacings and a complete discrete level sequence at low excitation energies. However, these parametrizations require the additional modeling of the dependence of the spin-cutoff parameter on excitation energy. Here, in this work, we avoid the need to model the spin distribution of level densities by using the experimental data to parametrize directly the state densities, for which the BBF does not depend on the spin-cutoff parameter. This approach allows for a local parameterization of state densities that is independent of the spin-cutoff parameter. We provide these parameters in a tabulated form for applications in nuclear reaction calculations and for testing microscopic approaches to state densities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Extended experimental inferential structure determination method in determining the structural ensembles of disordered protein states

Proteins with intrinsic or unfolded state disorder comprise a new frontier in structural biology, requiring the characterization of diverse and dynamic structural ensembles. Here we introduce a comprehensive Bayesian framework, the Extended Experimental Inferential Structure Determination (X-EISD) method, which calculates the maximum log-likelihood of a disordered protein ensemble. X-EISD accounts for the uncertainties of a range of experimental data and back-calculation models from structures, including NMR chemical shifts, J-couplings, Nuclear Overhauser Effects (NOEs), paramagnetic relaxation enhancements (PREs), residual dipolar couplings (RDCs), hydrodynamic radii (R h ), single molecule fluorescence Förster resonance energy transfer (smFRET) and small angle X-ray scattering (SAXS). We apply X-EISD to the joint optimization against experimental data for the unfolded drkN SH3 domain and find that combining a local data type, such as chemical shifts or J-couplings, paired with long-ranged restraints such as NOEs, PREs or smFRET, yields structural ensembles in good agreement with all other data types if combined with representative IDP conformers.

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

Fast permeability measurement for tight reservoir cores using only initial data of the one chamber pressure pulse decay test

Here, in this study, a mathematical model for fast determination of the permeabilities of tight rocks using measurements taken from the initial period of the One Chamber Pressure Pulse Decay (OC-PPD) test is presented. The model applies to measurements taken both before and after the pressure pulse front has reached the downstream end of the specimen. The analytical solutions for the pressure decay in the upstream chamber are derived based on a parabolic arc approximation of pore pressure distribution along the test specimen. This approximation allows converting the initial–boundary value problem of fluid diffusion in the specimen, governed by partial differential equations, to a system of ordinary differential equations that can be easily solved by explicit formulae. Thus, an explicit formula for the pressure decay rate is obtained, which enables inverse analysis of the initial experimental data to estimate the rock permeability. The proposed method expedites the pulse decay test as it does not require the system to reach equilibrium. The method is validated with three sets of experimental data of the OC-PPD test using helium as the diffusing fluid, for which the relative error of the permeability is found to be less than 6%. This method is particularly useful if the equilibrium time of the pulse decay test for rock specimens with permeabilities in the range of nano-Darcy takes hours or days.

early-time solution↗
DOE OSTI2020

Superallowed 0 + → 0 + nuclear β decays: 2020 critical survey, with implications for V ud and CKM unitarity

A new critical survey of all half-life, decay-energy and branching-ratio measurements related to 23 superallowed 0 + → 0 + β decays is presented. Included are 222 individual measurements of comparable precision obtained from 174 published references. Compared with our last survey in 2015, we have added results from 28 new publications and eliminated an approximately equal number whose results have been superseded by much more precise modern data. We obtain world-average ft values for each of the 21 transitions that have a complete set of data, then apply radiative and isospin-symmetry-breaking corrections to extract “corrected” Ft values. Fifteen of these Ft values now have a precision of 0.3% or better and all take the same value within statistics, as expected from conservation of the vector current. Their average, Ft¯, when combined with the muon lifetime, yields the up-down quark-mixing element of the Cabibbo-Kobayashi-Maskawa matrix, V ud = 0.97373 ± 0.00031. This is lower than our 2015 result by one standard deviation and its uncertainty is increased by 50%. This is a consequence, not of any shifts in the experimental data, but of new calculations for the radiative corrections. The lower V ud value now leads to greater tension in the top-row test of unitarity in the CKM matrix. Updates in experimental data have independently led to a factor-of-two tighter limit being set on the possible existence of a scalar interaction. In conclusion, the new limit on Fierz interference is b F ≤ 0.0033 at the 90% confidence level.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Geometrical defect detection for additive manufacturing with machine learning models

This study proposed a scheme based on Machine Learning (ML) models to detect geometric defects of additively manufactured objects. The ML models are trained with synthetic 3D point clouds with defects and then applied to detect defects in actual production. Using synthetic 3D point clouds rather than experimental data could save a huge amount of training time and costs associated with many prints for each design. Besides distance differences of individual points between source and target point clouds, this scheme uses a new concept called “patch” to capture macro-level information about nearby points for ML training and implementation. Numerical comparisons of prediction results on experimental data with different shapes showed that the proposed scheme outperformed the existing Z-difference method in the literature. Five ML methods (Bagging of Trees, Gradient Boosting, Random Forest, K-nearest Neighbors and Linear Supported Vector Machine) were compared under various conditions, such as different point cloud densities and defect sizes. Bagging and Random Forest were found the two best models regarding predictability; and the right patch size was found to be at 20. The proposed ML-based scheme is applicable to in-situ defect detection during additive manufacturing with the aid of a proper 3D data acquisition system.

Additive manufacturing↗
DOE OSTI2020

Isoscalar and Isovector Giant Resonances in Closed Shells Nuclei and Bulk Properties of Nuclear Matter

Centroid energies, E CEN , of the isoscalar (T = 0) and isovector (T = 1) giant resonances of multipolarities L = 0 – 3 in 40,48 Ca, 68 Ni, 90 Zr, 116 Sn, 144 Sm and 208 Pb, were calculated within the fully self-consistent spherical Hartree-Fock (HF)-based random phase approximation (RPA) theory, using 33 different energy density functionals associated with Skyrme-type effective nucleon-nucleon interactions of the standard form commonly employed in the literature. Here, we also calculate the Pearson linear correlation coefficients between each E CEN and each bulk property of nuclear matter (NM), associated with the Skyrme interactions used in the calculations, and determine the sensitivity of E CEN to bulk properties of NM. By comparing the calculated values of E CEN to the experimental data, we constrain the values of the bulk NM properties. We find that interactions associated with the values of the NM effective mass, m* / m = 0.70 to 0.90, incompressibility coefficient, K NM = 210 to 240 MeV, and the enhancement coefficient of the energy weighted sum rule of the isovector giant dipole resonance, κ = 0.25 to 0.70, best reproduce the experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Hyperfine interactions for small systems including transition-metal elements using self-interaction corrected density-functional theory

The interactions between the electronic magnetic moment and the nuclear spin moment, i.e., magnetic hyperfine (HF) interactions, play an important role in understanding electronic properties of magnetic systems and in realizing platforms for quantum information science applications. We investigate the HF interactions for atomic systems and small molecules, including Ti or Mn, by using Fermi–Löwdin orbital (FLO) based self-interaction corrected (SIC) density-functional theory. We calculate the Fermi contact (FC) and spin-dipole terms for the systems within the local density approximation (LDA) in the FLO-SIC method and compare them with the corresponding values without SIC within the LDA and generalized-gradient approximation (GGA), as well as experimental data. For the moderately heavy atomic systems (atomic number Z ≤ 25), we find that the mean absolute error of the FLO-SIC FC term is about 27 MHz (percentage error is 6.4%), while that of the LDA and GGA results is almost double that. Therefore, in this case, the FLO-SIC results are in better agreement with the experimental data. For the non-transition-metal molecules, the FLO-SIC FC term has the mean absolute error of 68 MHz, which is comparable to both the LDA and GGA results without SIC. For the seven transition-metal-based molecules, the FLO-SIC mean absolute error is 59 MHz, whereas the corresponding LDA and GGA errors are 101 and 82 MHz, respectively. Therefore, for the transition-metal-based molecules, the FLO-SIC FC term agrees better with experiment than the LDA and GGA results. We observe that the FC term from the FLO-SIC calculation is not necessarily larger than that from the LDA or GGA for all the considered systems due to the core spin polarization, in contrast to the expectation that SIC would increase the spin density near atomic nuclei, leading to larger FC terms.

Chemistry↗
DOE OSTI2025

Advancing Neutrino Simulation Modeling with MARLEY: Insights from the NNSA-MSIIP Internship

Core-collapse supernovae are intense sources of tens-of-MeV neutrinos. However, there is no experimental data to validate the current event generator MARLEY (Model of Argon Reaction Low Energy Yields) that can model supernova neutrinos. To validate the model, I developed a new muon capture feature within the MARLEY simulation framework by coding key functions in C++, Python, and ROOT. I generated and analyzed one million simulated muon capture events and comparing the results to experimental data. To improve the model, I worked on optimizing the model’s parameters to improve its precision using statistical methods.

Wong, Baker [Fermilab]↗
DOE OSTI2025

Describing hadronization via histories and observables for Monte-Carlo event reweighting

We introduce a novel method for extracting a fragmentation model directly from experimental data without requiring an explicit parametric form, called Histories and Observables for Monte-Carlo Event Reweighting (HOMER), consisting of three steps: the training of a classifier between simulation and data, the inference of single fragmentation weights, and the calculation of the weight for the full hadronization chain. We illustrate the use of HOMER on a simplified hadronization problem, a q\bar{q} q q ‾ string fragmenting into pions, and extract a modified Lund string fragmentation function f(z) f ( z ) . We then demonstrate the use of HOMER on three types of experimental data: (i) binned distributions of high-level observables, (ii) unbinned event-by-event distributions of these observables, and (iii) full particle cloud information. After demonstrating that f(z) f ( z ) can be extracted from data (the inverse of hadronization), we also show that, at least in this limited setup, the fidelity of the extracted f(z) f ( z ) suffers only limited loss when moving from (i) to (ii) to (iii). Public code is available at https://gitlab.com/uchep/mlhad.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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↗
DOE OSTI2021

Experimental and modeling studies of PR and ND oxalate solubility to high ionic strengths: Insight into actinide(III) oxalates

Actinide oxalates are chemical compounds important to nuclear industry, ranging from actinide separation in waste reprocessing, to production of specialty actinides, and to disposal of high level nuclear waste (HLW) and spent nuclear fuel (SNF). In this study, the solubility constants for Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O by performing solubility experiments in HNO 3 and mixtures of HNO 3 and H 2 C 2 O 4 at 23.0 ± 0.2 °C have been determined. The targeted starting materials, Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O, were successfully synthesized at room temperature using PrCl 3 , NdCl 3 and oxalic acid as the source metrials. Then, we utilized the targeted solubility-controlling phases to conduct solubility measurements. There was no phase change over the entire periods of experiments, demonstrating that Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O were the solubility-controlling phases in our respective experiments. Based on our experimental data, we have developed a thermodynamic model for Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O in the mixtures of HNO 3 and H 2 C 2 O 4 to high ionic strengths. The model for Pr 2 (C 2 O 4 ) 3 •10H 2 O reproduces well the reported experimental data for Pu 2 (C 2 O 4 ) 3 •10H 2 O, which are not utilized for the model development, demonstrating that Pr(III) is an excellent analog for Pu(III). Similarly, the model for Nd 2 (C 2 O 4 ) 3 •10H 2 O reproduces the solubility of Am 2 (C 2 O 4 ) 3 •10H 2 O and Cm 2 (C 2 O 4 ) 3 •10H 2 O. The Pitzer model was used for the calculation of activity coefficients. Based on the published, well established model for dissociation constants for oxalic acid and stability constants for actinide-oxalate complexes [i.e., AmC 2 O 4 + , and Am(C 2 O 4 ) 2 - ] to high ionic strengths, we have obtained the solubility constants (log 10 K 0 ) for the following reactions at 25 °C, Pr 2 (C 2 O 4 ) 3 •10H 2 O ⇌ 2Pr 3+ + 3C 2 O 4 2- + 10H 2 O(l). Nd 2 (C 2 O 4 ) 3 •10H 2 O ⇌ 2Nd 3+ + 3C 2 O 4 2- + 10H 2 O(l). to be -30.82 ± 0.30 (2σ), and - 31.14 ± 0.35 (2σ), respectively. These values for can be directly applied to Pu 2 (C 2 O 4 ) 3 •10H 2 O, Am 2 (C 2 O 4 ) 3 •10H 2 O and Cm 2 (C 2 O 4 ) 3 •10H 2 O. The model established for actinide oxalates by this study provides the needed knowledge with regard to solubilities of actinide/REE oxalates at various ionic strengths, and is expected to find applications in many fields, including the geological disposal of nuclear waste and the mobility of REE under the surface conditions, as Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O can be regarded as the pure Pr and Nd end-members of deveroite, a recently discovered natural REE oxalate with the following stoichiometry, (Ce 1.01 Nd 0.33 La 0.32 Pr 0.11 Y 0.11 Sm 0.01 Pb 0.04 U 0.03 Th 0.01 Ca 0.04 ) 2.01 (C 2 O 4 ) 2.99 •9.99H 2 O. Regarding its importance in the geological disposal of nuclear waste, Am 2 (C 2 O 4 ) 3 •10H 2 O/Pu 2 (C 2 O 4 ) 3 •10H 2 O/Cm 2 (C 2 O 4 ) 3 •10H 2 O can be the source-term phase for actinides, as demonstrated by the instance in the disposal in clay/shale formations. This is exemplified by the stability of Am 2 (C 2 O 4 ) 3 •10H 2 O in comparison with Am(OH) 3 (am), Am(OH) 3 (s) and AmCO 3 (OH)(s) under the relevant geological repository conditions.

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

Thermodynamic Modeling of Aqueous LiCl, LiBr, LiI, and LiNO 3 Solutions

Thermodynamic models are essential to facilitate the advancement of process design, optimization, and operation of electrolyte systems. In this work, a comprehensive thermodynamic framework based on the Electrolyte Nonrandom Two-Liquid model is developed to calculate phase equilibria behavior and salt solubility of aqueous LiCl, LiBr, LiI, and LiNO 3 solutions. The model describes the non-ideality of the electrolyte solutions by using two binary interaction parameters for each electrolyte-molecule pair in the system. Each binary interaction parameter is further expressed with up to three temperature coefficients which are regressed from experimental data. To take into account the hydration of lithium ion, two separate chemistries for the dissociation of lithium salts are investigated. In the first case, the lithium ion is considered as a bare ion, Li + , while in the second case hydration of the lithium ion from Li + to Li(H 2 O) + is considered. Here, the calculated thermodynamic properties compare adequately with the experimental data for both sets of chemistries for concentrations up to saturation and temperatures from 273.15 K up to 623.15 K. Moderate to significant improvements are observed with the incorporation of the hydration chemistry for aqueous LiCl, LiBr, and LiI solutions when compared to the non-hydrated lithium ion model results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗