Search NASA⌕ Search

SEARCH · Search NASA

Results for “Thermodynamic 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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Performance of sweeping gas membrane distillation for treating produced water: Modeling and experiments

Membrane distillation (MD) has shown promise for desalination of produced water (PW), but there are very few studies reported on the performance of MD with real PW and the characteristics of the fouling and scaling in MD. Among different MD configurations, sweeping gas MD (SGMD) is the least studied configuration. Here we report the application of SGMD for treating high-salinity PW (123 g L -1 ). The SGMD operating conditions were first carefully optimized through pure water MD experiments by varying operating conditions. The flux evaluation was accompanied by a detailed Aspen Custom Modeler simulation based on the electrolyte-NRTL thermodynamic model to better understand the heat and mass transfer in the feed and permeate channel. Supported by experiments, modeling results confirmed the partial condensation of the vapor in the permeate channel. Furthermore, it was concluded that the saturation of the sweeping gas is controlling the overall mass transfer through the membrane. Selected membranes and optimized conditions were then used to conduct short-term and long-term performance tests with PW. In long-term MD test, the concentration of the PW increased by a factor of two (50% water recovery), while the membrane exhibited small flux decline. Elemental analysis showed that strontium and sodium were the major scalant compounds.

42 ENGINEERING↗
DOE OSTI2021

Modeling Phase Equilibrium of Common Sugars Glucose, Fructose, and Sucrose in Mixed Solvents

The industrial processing of sugars and sugar-containing mixtures is gaining widespread use as a form of renewable manufacturing from biomass. To aid in process modeling for design and optimization, a commonly available thermodynamic model is needed that describes the phase equilibrium of these compounds. This work compiles and compares models for solid–liquid and vapor–liquid phase equilibrium from the available data for the representative sugars glucose, fructose, and sucrose in the representative solvents water, methanol, and ethanol, including data for multisugar, multisolvent systems. Additionally, the nonrandom two-liquid (NRTL) model was chosen for these systems because of its widespread use in industry and the availability of parameters for many solvents and cosolutes. The association-NRTL (aNRTL) model was investigated as an improvement for modeling sugars, which may experience a high degree of association because of their many hydroxy groups. Both models accurately capture the data and are able to predict the behavior of multisugar systems with only solute–solvent interaction parameters. The aNRTL model shows an improvement over the baseline NRTL model that is most significant for sucrose, the component with the highest association strength, and least significant for fructose, which has the lowest association strength.

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

First-Principles Thermodynamic Assessments of Sr-Containing Secondary Phase Formation in La1-xSrxMnO3±δ Perovskites for Solid Oxide Cell Applications

Sr-secondary phase formation is a potentially significant degradation mode threatening solid-oxide cell (SOC) commercial viability. A first-principles thermodynamic study was performed for rhombohedral perovskite (La1-xSrx) MnO3±δ (LSM) to assess its stability against Sr secondary phase formation in SOC applications. In this work, the Sr secondary phase formation reaction free energies were determined by combining ab initio lattice dynamics calculations for the solid phases and an ab initio thermodynamics approach for the gas phases. Furthermore, this approach goes beyond previous thermodynamic modeling studies by integrating first-principles based point-defect equilibria into the analyses. The modeling results indicate an increased tendency to form SrO oxide from LSM upon decreasing the oxygen partial pressure. Additionally, enhancing factors to form the Sr-related secondary phase from the associated SrO activity in LSM are further quantified by considering the equilibrium of SrO reacting with contaminant gas species as a function of temperature and gas pressure.

Defect and phase stability↗
DOE OSTI2022

Generalization of thermodynamic Langmuir isotherm for mixed‐gas adsorption equilibria

Abstract This work presents a comprehensive thermodynamic model for both pure component isotherms and mixed‐gas adsorption equilibria. A generalization of thermodynamic Langmuir isotherm, the proposed model assumes competitive adsorption of multiple adsorbates on adsorbent surface for mixed‐gas adsorption equilibria, and it applies an area‐based adsorption nonrandom two‐liquid activity coefficient model in the activity coefficient calculations for the adsorbate phase. The resulting generalized Langmuir isotherm properly captures both surface loading dependence and adsorbate phase composition dependence for mixed‐gas adsorption equilibria. The model is validated with accurate representations of gas adsorption equilibrium data for varieties of unary, binary, and ternary gas systems. The model results are further compared with those calculated from extended Langmuir isotherm and Ideal Adsorbed Solution Theory.

Hamid, Usman↗
DOE OSTI2024

Lambda-PFLOTRAN 1.0: a workflow for incorporating organic matter chemistry informed by ultra high resolution mass spectrometry into biogeochemical modeling

Abstract. Organic matter (OM) composition plays a central role in microbial respiration of dissolved organic matter and subsequent biogeochemical reactions. Here, a direct connection of organic matter chemistry and thermodynamics to reactive transport simulators has been achieved through the newly developed Lambda-PFLOTRAN workflow tool that succinctly incorporates carbon chemistry data generated from Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) into reaction networks to simulate organic matter degradation and the resulting biogeochemistry. Lambda-PFLOTRAN is a Python-based workflow, executed through a Jupyter notebook interface, that digests raw FTICR-MS data, develops a representative reaction network based on substrate-explicit thermodynamic modeling (also termed lambda modeling due to its key thermodynamic parameter λ used therein), and completes a biogeochemical simulation with the open source, reactive flow and transport code PFLOTRAN. The workflow consists of the following five steps: configuration, thermodynamic (lambda) analysis, sensitivity analysis, parameter estimation, and simulation output and visualization. Two test cases are provided to demonstrate the functionality of the Lambda-PFLOTRAN workflow. The first test case uses laboratory incubation data of temporal oxygen depletion to fit lambda parameters (i.e., maximum utilization rate and microbial carrying capacity). A slightly more complex second test case fits multiple lambda formulation and soil organic matter release parameters to temporal greenhouse gas generation measured during a soil incubation. Overall, the Lambda-PFLOTRAN workflow facilitates upscaling by using molecular-scale characterization to inform biogeochemical processes occurring at larger scales.

58 GEOSCIENCES↗
DOE OSTI2020

Nonrandom two‐liquid activity coefficient model with association theory

Abstract Considered one of the most versatile and widely used classical thermodynamic models to correlate phase equilibria behavior of nonideal systems, the nonrandom two‐liquid (NRTL) theory does not explicitly account for specific chemical associations such as hydrogen bonding. This deficiency has been recognized as the cause for unsatisfactory representation of association systems like methanol–alkanes binaries. This work presents a practical approach to integrate Wertheim's perturbation theory for association contribution with the classical NRTL model. Specifically, the association contribution is calculated from pre‐determined molecule‐specific association strengths while the physical interaction contribution is captured with NRTL binary interaction parameters. The resulting association NRTL model correlates fluid phase equilibria of association systems with few adjustable parameters and offers improved predictive capability for higher order systems.

Hao, Yifan↗
DOE OSTI2025

An eigenvalue-based method for computing the relaxed pressure in compressible multiphase flow with N phases

The modeling of compressible multiphase flows is a decades-old area of study with many applications across various fields. Many of these application areas use stiff pressure relaxation. This process involves the solution of a nonlinear system with N + 1 equations and N + 1 unknowns, where N is the number of phases. The resolution of this system with general equations of state (EOSs) is difficult. Furthermore, nonlinear systems can admit multiple solutions, and current solution methods do not address this possibility. Very recently, a thermodynamic relaxation method was introduced, which effectively maps a relatively simple predictor equation of state onto a more complex target equation of state. In this context, the target EOSs are the chosen EOSs for the thermodynamic model. Furthermore, this thermodynamic relaxation has the benefit of simplifying the stiff pressure relaxation system of equations. In this article, we show this system reduces to a polynomial of degree N, which can be recast as an eigenvalue problem through the use of the associated companion matrix. We show that although this eigenvalue method is generally less efficient than Newton–Raphson iteration, it does not suffer from convergence issues and finds all N roots of the polynomial. Hence, the method provides a fail-safe for root-finding iterative methods and a way to address the issue of multiple solutions to the nonlinear system of equations in stiff pressure relaxation.

Eigenvalue algorithm↗
DOE OSTI2021

Low global warming potential (GWP) refrigerant supermarket refrigeration system modeling and its application

As an environmentally friendly low global warming potential (GWP) refrigerant, Carbon dioxide (CO 2 ) has continuously gained popularity and research attention as alternative refrigerant for supermarket refrigeration system. In this paper, to fulfill the increasing need of accurate Low-GWP supermarket refrigeration models for development of supervisory level control and optimization strategies, a high fidelity model is developed for CO 2 transcritical supermarket refrigeration system which includes compressor rack of low temperature (LT) compressors and medium temperature (MT) compressors, air-cooled gas cooler, evaporator, expansion valves and other auxiliary equipment. A resistance-capacity model structure is proposed to simulate the display cases. Semi-thermodynamic models are proposed to estimate reciprocating compressors volumetric efficiency and power consumption. Furthermore the zone modeling approach is used for evaporator simulation, and air-cooled gas cooler is modeled with distributed modeling method. The expansion valve simulation is based on orifice flow model. To calibrate these models, both manufacture performance data and experimental data are used. The experiments are conducted with a full instrumental CO 2 supermarket refrigeration system installed in Oak Ridge National Lab Flexible Research Platform (FRP). The simulation model can predict the system performance, including power consumption, cooling capacity, mass flow rate, temperature, and pressure, with high accuracy (within ±4%) compared to experimental data. In addition, this developed model has been used to create the system optimum high side pressure for high side expansion valve control, and to generate wide operating range simulation data for developing the batter-equivalent commercial refrigeration system model which can be used in grid interactive control development.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Developing Practical Models of Complex Salts for Molten Salt Reactors

Molten salt reactors (MSRs) utilize salts as coolant or as the fuel and coolant together with fissile isotopes dissolved in the salt. It is necessary to therefore understand the behavior of the salts to effectively design, operate, and regulate such reactors, and thus there is a need for thermodynamic models for the salt systems. Molten salts, however, are difficult to represent as they exhibit short-range order that is dependent on both composition and temperature. A widely useful approach is the modified quasichemical model in the quadruplet approximation that provides for consideration of first- and second-nearest-neighbor coordination and interactions. Its use in the CALPHAD approach to system modeling requires fitting parameters using standard thermodynamic data such as phase equilibria, heat capacity, and others. A shortcoming of the model is its inability to directly vary coordination numbers with composition or temperature. Another issue is the difficulty in fitting model parameters using regression methods without already having very good initial values. The proposed paper will discuss these issues and note some practical methods for the effective generation of useful models.

Besmann, Theodore M. (ORCID:0000000155980550)↗
DOE OSTI2021

A method to apply Zavitsas' aqueous electrolyte model to multicomponent solutions and its equivalence to Zdanovskii's rule

Abstract Zavitsas developed an aqueous electrolyte thermodynamic model and applied it to single electrolyte solutions previously. Here, a method to apply the model to multicomponent solutions is derived by setting the water activity to the mole fraction of free water in the mixture. The method is shown to work effectively for the CsBr‐KNO 3 ‐H 2 O and CsCl‐KNO 3 ‐H 2 O systems at 100.3°C even at concentrations greater than 6 molal. The largest error in the prediction of water activities was just −0.00365 for the CsBr‐KNO 3 ‐H 2 O system and −0.00631 for the CsCl‐KNO 3 ‐H 2 O system. Zavitsas' model naturally obeys Zdanovskii's rule because when two electrolyte solutions with the same fraction of free water are mixed together, water is obviously going to have that same fraction in the mixture. This means that Zavitsas' model will likely work well for the many electrolytes that obey Zdanovskii's rule.

Reynolds, Jacob G.↗
DOE OSTI2025

Experiments to validate Thermodynamics and Transport models of Strongly Coupled Dusty Plasma Matter (Final Technical Report for DE-SC0023416)

The goal of this two-year grant is to provide access to the PI to dusty plasma experimental facilities at the DOE-funded Collaborative Research Facility Magnetized Plasma Research Laboratory, Auburn University to become a user of that facility, to obtain experimental data to support another ongoing grant DE-SC0021146 (an Early Career Award to the PI that is focused on modeling of dusty plasma thermodynamics and transport processes), to generate experimental data for funding proposals, and to provide exposure to University of Memphis students to advanced experimental techniques. The following technical accomplishments were made: 1. Development of a novel Bidirectional Electrode Control Arms Assembly (BECAA) for producing perfect 2D grain layers for complex plasma experimentation. BECAA uses movable electrode arms to tilt or move the electrode in a RF discharge from outside the chamber, allowing for the manipulation of grain clouds without needing to change the plasma parameters or gas pressure. This work addresses a longstanding gap in the literature for a method to produce clusters of selectable number of grains and that are perfectly two dimensional as opposed to being only quasi-2D. 2. Experimental investigation of the structural properties of finite-N clusters with N=2 to 50. Individual particle behavior in clusters could vary from grain to grain and this study measured systematically produced clusters for two different grain sizes. Analysis (funded by another grant DE-SC0021146) is currently underway to quantify the differences between grains that are found on the surface vs. the interior of clusters, the shell structure, and the decay of correlations in position, velocity, and kinetic energy. 3. An experimental method to measure the structural entropy of clusters was developed by observing the self-induced structural transitions between various possible arrangements. In a series of heating and cooling cycles, the number of times each possible arrangement was attained was experimentally observed and used to compute the probability of existence of that arrangement, and subsequently the configurational entropy of the cluster. Analysis (funded by another grant DE-SC0021146) is currently underway to produce the entropy of clusters as a function of the number of grains and use the same to compute thermodynamic state variables for 2D complex plasma/grain clusters. 4. A preliminary experimental study of multibody collisions between N grains (N=2 – 10) was conducted. The clusters were produced using the BECAA technique and velocities were imparted to the grains using manipulation laser pulses. Analysis (funded by another grant DE-SC0021146) is currently underway to develop a theoretical framework to describe multibody collisions analogous to classical two-body interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Topical Report on Water Balance Model Development and Initial Application

This topical report summarizes results from developing a generic dynamic water balance tool which can be used by coal-fired power plant operators to evaluate FGD water management strategies. EPRI is using a statistical software model created on the GoldSim platform and data from coal-fired energy plants to determine the fraction of power plant water needs that can be satisfied by electrodialysis reversal treatment and, correspondingly, the discharge volume and/or byproducts that still need to be managed. The generic model is designed to accept a number of easily obtainable process variables and generate probabilities of resultant water characteristics. The model is also constructed as a universal process flow; one that can be adjusted to accommodate most major components of an FGD blowdown handling system. The generic model is dynamic and utilizes user-inputted operational variation in water quality and quantity to best provide a probabilistic representation of potential performance outcomes. In order to account for complexities of water composition and treatment, a separate chemical modeling software is used to evaluate water chemistry and unit operation performance. The thermodynamic modeling software, OLI Studio, is employed to convert ion data from laboratory analyses into balanced molecular concentrations that are classified as dissolved and suspended solids. The dynamic water model is designed to import a total of 11 species commonly found in FGD blowdown water as molecular concentrations.Water management unit operations are used in the model and categorized by either physical unit operation blocks or chemical unit operations blocks. As the name implies, the physical water treatment blocks are unit operations that are exclusively physical changes to the system (e.g., separations and volume changes). The chemical water treatment blocks are unit operations where chemical changes occur and require the use of OLI Studio to characterize subsequent streams. Water and material balances from three unique coal-fired power plants were obtained and preliminary modeling on the dynamic water balance tool was performed using the extracted data from these power plants. These three power plants were chosen due to their unique configurations and water profiles that represented a variety of operational changes. Results from using the dynamic tool with the three plant configurations are presented and discussed in the context of the existing tool and planned future capabilities.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Modeling dissociation of ionic liquids with electrolyte NRTL model

Dissociation of ionic liquids (ILs) is critical to electrochemical process performance by affecting ionic conductivity. Prior studies of thermodynamic modeling either assumed no dissociation or complete dissociation that were not consistent with experimental observations. This work presented an eNRTL model for three aqueous imidazolium-based IL systems: [Emim][EtSO 4 ], [Emim][TFA], and [Emim][TfO]. The model accurately correlates the IL dissociation and VLE data in the entire concentration range from pure IL to infinite dilution aqueous solution. The system non-ideality is dominated by the dissociation chemistry. Dissociation behavior is driven by strong short-range interaction of the molecular IL-water pair relative to the dissociated IL-water pair, consistent with prior studies indicating that molecular IL and water can form an energetically favorable complex. The interaction energy parameters calculated in this work qualitatively agrees with the hydrogen bonding ability of the anions previously reported. This research not only provides a model that well describes the IL dissociation but also demonstrates the physical significance of eNRTL model.

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

Optimization of Desalination Systems with Detailed Water Chemistry through Integration of Reaktoro in WaterTAP

Chemistry predictions are critical for an accurate estimation of performance and costs in desalination process models, which allows for the estimation of the value of new technologies and the viability of treating new water sources. Herein, we present how an implicit function formulation can be used to integrate the chemical modeling package, Reaktoro, into the techno-economic assessment and modeling platform, WaterTAP. This approach resolves the critical issues of integrating large-scale thermodynamic models and databases into equation-oriented process models while allowing more flexibility relative to previously presented surrogate-based methods. We describe how this integration into Pyomo and WaterTAP models is implemented and used through the open-source package Reaktoro-PSE . We first validate this integration approach by performing optimization on a previously presented desalination treatment train with softening and acid addition as the pretreatment steps. Then, to demonstrate the value of this approach, we extend the cost-optimization problem to include the simultaneous addition of lime and soda ash for softening, and HCl and H 2 SO 4 in the acidification steps. Finally, we were able to confirm the previously established results that were obtained by using surrogate models and demonstrate that the implicit function approach enables exploration of different feedwater compositions and a larger number of chemicals and their combinations.

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

Uncertainty propagation in pore water chemical composition calculation using surrogate models

Performance assessment in deep geological nuclear waste repository systems necessitates an extended knowledge of the pore water chemical conditions prevailing in host-rock formations. In the last two decades, important progress has been made in the experimental characterization and thermodynamic modeling of pore water speciation, but the influence of experimental artifacts and uncertainties of thermodynamic input parameters are seldom evaluated. In this respect, we conducted an uncertainty propagation study in a reference geochemical model describing the pore water chemistry of the Callovian-Oxfordian clay formation. Nineteen model input parameters were perturbed, including those associated to experimental characterization (leached anions, exchanged cations, cation exchange selectivity coefficients) and those associated to generic thermodynamic databases (solubilities). A set of 13 quantities of interest were studied by the use of polynomial chaos expansions built non-intrusively with a least-squares forward stepwise regression approach. Training and validation sets of simulations were carried out using the geochemical speciation code PHREEQC. The statistical results explored the marginal distribution of each quantity of interest, their bivariate correlations as well as their global sensitivity indices. The influence of the assumed distributions for input parameters uncertainties was evaluated by considering two parametric domain sizes.

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

Thermodynamic re-modeling of the Yb-Sb system aided by first-principles calculations

Here, the thermodynamic description of the Yb-Sb binary system is developed by means of the CALculations of PHAse Diagrams (CALPHAD) method by combining experimental data in the literature and predictions from first-principles calculations based on density functional theory (DFT) in the literature and the present work. Two pseudopotentials of Yb are compared in the present DFT-based calculations with 14 and 13 f-electrons frozen in the core, i.e., 5p 6 6s 2 and 5p 6 6s 2 5 d 1 electrons as valence electrons, termed Yb_2 and Yb_3, respectively. It is shown that the phonon spectrum of the YbSb phase calculated using the Yb_3 pseudopotential does not have imaginary phonon modes and is subsequently used to predict its temperature dependent thermodynamic properties by the DFT-based quasiharmonic phonon calculations. The present thermodynamic database includes the Yb 16 Sb 11 phase in addition to five intermetallic phases that were considered in previous modeling studies, i.e., YbSb 2 , YbSb, Yb 11 Sb 10 , Yb 4 Sb 3 , and Yb 5 Sb 3 . The high temperature orthorhombic structure of the Yb 5 Sb 3 phase is not considered in the present work as it was stabilized by hydrogen. The associate solution model is used to describe the short-range ordering behavior in the liquid phase. The calculations from the present thermodynamic model show good agreement with thermochemical and phase equilibrium data from both the present work and the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A purely kinetic description of the evaporation of water droplets

The process of water evaporation, although deeply studied, does not enjoy a kinetic description that captures known physics and can be integrated with other detailed processes such as drying of catalytic membranes embedded in vapor-fed devices and chemical reactions in aerosol whose volumes are changing dynamically. In this work, we present a simple, three-step kinetic model for water evaporation that is based on theory and validated by using well-established thermodynamic models of droplet size as a function of time, temperature, and relative humidity as well as data from time-resolved measurements of evaporating droplet size. Additionally, the kinetic mechanism for evaporation is a combination of two limiting processes occurring in the highly dynamic liquid-vapor interfacial region: direct first order desorption of a single water molecule and desorption resulting from a local fluctuation, described using third order kinetics. The model reproduces data over a range of relative humidities and temperatures only if the interface that separates bulk water from gas phase water has a finite width, consistent with previous experimental and theoretical studies. The influence of droplet cooling during rapid evaporation on the kinetics is discussed; discrepancies between the various models point to the need for additional experimental data to identify their origin.

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

A Landau-Devonshire Approach of Elastic Strain Effects on Ferroelectric Al1-xScxN

Landau-Devonshire thermodynamic modeling enables the phenomenological description of ferroelectrics. 1 The approach has been validated across a wide range of ferroelectric materials from the first ferroelectric discovered, Rochelle Salt, to recent HfO2-based materials. 2-4 Extension to what is now referred to as the Landau-Ginsburg-Devonshire model enables rigorous investigation of polarization boundaries and reorientation, phase transformations, and domain formation/evolution. We present a thermodynamic analysis of the recently discovered nitride ferroelectric materials 5-7 using the classic Landau-Devonshire approach. The electrostrictive and dielectric stiffness coefficients of Al1-xScxN with wurtzite structure (6mm) are determined using a free energy density function assuming a hexagonal parent phase (6/mmm), with the first order phase transition based on the dielectric stiffness relationships (Fig. 1). The results of this analysis show that the strain sensitivity of the energy barrier is one order of magnitude larger than that of the spontaneous polarization in these wurtzite ferroelectrics, yet both are less sensitive to strain compared to classic perovskite ferroelectrics (Fig. 2). These analysis results reported here explain experimentally reported sensitivity of coercive field to elastic strain/stress in Al1-xScxN films, and would enable further thermodynamic analysis via phase field simulation and related methods.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗