Temperature-dependent solubility of Rebaudioside A in methanol/ethanol and ethyl acetate mixtures:Experimental measurements and thermodynamic modeling

The equilibrium solubility of Rebaudioside A(Reb A)FormⅡin binary mixtures of methanol/ethanol and ethyl acetate was quantitatively determined within the temperature range of 283.15—328.15 K at ambient pressure.The experimental findings indicate a positive correlation between the solubility of Reb A(FormⅡ)and both the temperature and the methanol/ethanol content in the solvent system.To describe the solubility data,six distinct models were employed:the modified Apelblat equation,theλh model,the combined nearly ideal binary solvent/Redlich—Kister(CNIBS/R—K)model,the van't HoffJouyban-Acree(VJA)model,the Apelblat-Jouyban-Acree(AJA)model,and the non-random two-liquid(NRTL)model.The combined nearly ideal binary solvent/Redlich—Kister model exhibited the most precise fit for solubility in methanol+ethyl acetate mixtures,reflected by an average relative deviation(ARD)of 0.0011 and a root mean square deviation(RMSD)of 12×10^(-7).Conversely,for ethanol+ethyl acetate mixtures,the modified Apelblat equation provided a superior correlation(ARD=0.0014,RMSD=4×10^(-7)).Furthermore,thermodynamic parameters associated with the dissolution of Reb A(FormⅡ),including enthalpy,entropy,and the Gibbs energy change,were inferred from the data.The findings underscore that the dissolution process is predominantly endothermic across the solvent systems examined.Notably,the entropy changes appear to have a significant influence on the Gibbs free energy associated with the dissolution of Reb A(FormⅡ),suggesting that entropic factors may play a pivotal role in the studied systems.

Machine learning prediction of methane,ethane,and propane solubility in pure water and electrolyte solutions:Implications for stray gas migration modeling

Hydraulic fracturing is an effective technology for hydrocarbon extraction from unconventional shale and tight gas reservoirs.A potential risk of hydraulic fracturing is the upward migration of stray gas from the deep subsurface to shallow aquifers.The stray gas can dissolve in groundwater leading to chemical and biological reactions,which could negatively affect groundwater quality and contribute to atmospheric emissions.The knowledge oflight hydrocarbon solubility in the aqueous environment is essential for the numerical modelling offlow and transport in the subsurface.Herein,we compiled a database containing 2129experimental data of methane,ethane,and propane solubility in pure water and various electrolyte solutions over wide ranges of operating temperature and pressure.Two machine learning algorithms,namely regression tree(RT)and boosted regression tree(BRT)tuned with a Bayesian optimization algorithm(BO)were employed to determine the solubility of gases.The predictions were compared with the experimental data as well as four well-established thermodynamic models.Our analysis shows that the BRT-BO is sufficiently accurate,and the predicted values agree well with those obtained from the thermodynamic models.The coefficient of determination(R2)between experimental and predicted values is 0.99 and the mean squared error(MSE)is 9.97×10^(-8).The leverage statistical approach further confirmed the validity of the model developed.

Phase equilibrium data prediction and process optimizationin butadiene extraction process

In response to the lack of reliable physical parameters in the process simulation of the butadiene extraction,a large amount of phase equilibrium data were collected in the context of the actual process of butadiene production by acetonitrile.The accuracy of five prediction methods,UNIFAC(UNIQUAC Functional-group Activity Coefficients),UNIFAC-LL,UNIFAC-LBY,UNIFAC-DMD and COSMO-RS,applied to the butadiene extraction process was verified using partial phase equilibrium data.The results showed that the UNIFAC-DMD method had the highest accuracy in predicting phase equilibrium data for the missing system.COSMO-RS-predicted multiple systems showed good accuracy,and a large number of missing phase equilibrium data were estimated using the UNIFAC-DMD method and COSMO-RS method.The predicted phase equilibrium data were checked for consistency.The NRTL-RK(non-Random Two Liquid-Redlich-Kwong Equation of State)and UNIQUAC thermodynamic models were used to correlate the phase equilibrium data.Industrial device simulations were used to verify the accuracy of the thermodynamic model applied to the butadiene extraction process.The simulation results showed that the average deviations of the simulated results using the correlated thermodynamic model from the actual values were less than 2%compared to that using the commercial simulation software,Aspen Plus and its database.The average deviation was much smaller than that of the simulations using the Aspen Plus database(>10%),indicating that the obtained phase equilibrium data are highly accurate and reliable.The best phase equilibrium data and thermodynamic model parameters for butadiene extraction are provided.This improves the accuracy and reliability of the design,optimization and control of the process,and provides a basis and guarantee for developing a more environmentally friendly and economical butadiene extraction process.

FLOCKING OF A THERMODYNAMIC CUCKER-SMALE MODEL WITH LOCAL VELOCITY INTERACTIONS

In this paper, we study the flocking behavior of a thermodynamic Cucker–Smale model with local velocity interactions. Using the spectral gap of a connected stochastic matrix, together with an elaborate estimate on perturbations of a linearized system, we provide a sufficient framework in terms of initial data and model parameters to guarantee flocking. Moreover, it is shown that the system achieves a consensus at an exponential rate.

Parameterization, sensitivity, and uncertainty of 1-D thermodynamic thin-ice thickness retrieval

Retrieval of Thin-Ice Thickness(TIT)using thermodynamic modeling is sensitive to the parameterization of the independent variables(coded in the model)and the uncertainty of the measured input variables.This article examines the deviation of the classical model’s TIT output when using different parameterization schemes and the sensitivity of the output to the ice thickness.Moreover,it estimates the uncertainty of the output in response to the uncertainties of the input variables.The parameterized independent variables include atmospheric longwave emissivity,air density,specific heat of air,latent heat of ice,conductivity of ice,snow depth,and snow conductivity.Measured input parameters include air temperature,ice surface temperature,and wind speed.Among the independent variables,the results show that the highest deviation is caused by adjusting the parameterization of snow conductivity and depth,followed ice conductivity.The sensitivity of the output TIT to ice thickness is highest when using parameterization of ice conductivity,atmospheric emissivity,and snow conductivity and depth.The retrieved TIT obtained using each parameterization scheme is validated using in situ measurements and satellite-retrieved data.From in situ measurements,the uncertainties of the measured air temperature and surface temperature are found to be high.The resulting uncertainties of TIT are evaluated using perturbations of the input data selected based on the probability distribution of the measurement error.The results show that the overall uncertainty of TIT to air temperature,surface temperature,and wind speed uncertainty is around 0.09 m,0.049 m,and−0.005 m,respectively.

Thermodynamic modeling and phase diagram prediction of salt lake brine systems.Ⅰ. Aqueous Mg^2+–Ca^2+–Cl^- binary and ternary systems

Salt lake brine is a complex salt-water system under natural environment.Although many models can express the thermodynamic properties and phase equilibrium of electrolyte aqueous solution,the multi-temperature characteristics and predictability are still the goals of model development.In this study,a comprehensive thermodynamic model system is re-established based on the eNRTL model and some improvements:(1) new expression of long-range electrostatic term with symmetrical reference state is proposed to handle the electrolyte solution covering entire concentration range;(2) the temperature dependence of the binary interaction parameters is formulated with a Gibbs Helmholtz expression containing three temperature coefficients,the liquid parameters,which associated with Gibbs energy,enthalpy,and heat capacity contribution;and(3) liquid parameters and solid species data are regressed from properties and solubility data at full temperature range.Together the activity coefficient model,property models and parameters of liquid and solid offer a comprehensive thermodynamic model system for the typical bittern of MgCl2-CaCl2-H2 O binary and ternary systems,and it shows excellent agreement with the literature data for the ternary and binary systems.The successful prediction of complete phase diagram of ternary system shows that the model has the ability to deal with high concentration and high non-idealitv system,and the ability to extrapolate the temperature.

Thermodynamic modeling and phase diagram prediction of salt lake brine systemsⅡ.Aqueous Li^(+)-Na^(+)-K^(+)-SO_(4)^(2-) and its subsystems

It is still a challenging task to accurately and temperature-continuously express the thermodynamic properties and phase equilibrium behaviors of the salt-lake brine with multi-component,multitemperature and high concentration.The essential subsystem of sulfate type brine,aqueous Li^(+)-Na^(+)-K^(+)-SO_(4)^(2-) and its subsystems across a temperature range from 250 K to 643 K are investigated with the improved comprehensive thermodynamic model.Liquid parameters(Δg_(IJ),Δh_(IJ),and ΔC_(p,IJ))associated with the contributions of Gibbs energy,enthalpy,and heat capacity to the binary interaction parameters,i.e.the temperature coefficients of eNRTL parameters formulated with a Gibbs Helmholtz expression,are determined via multi-objective optimization method.The solid constantsΔ_(f)G_(k)°^((298.15))andΔ_(f)H_(k)°^((298.15))of11 solid species occurred in the quaternary system are rebuilt from multi-temperature solubilities.The modeling results show the accurate representation of(1)solution properties and binary phase diagram at temperature ranges from eutectic points to 643 K;(2)isothermal phase diagrams for Li_(2)SO_(4)-Na_(2)SO_(4)-H_(2)O,Li_(2)SO_(4)-K_(2)SO_(4)-H_(2)O and Na_(2)SO_(4)-K_(2)SO_(4)-H_(2)O ternary systems.The predicted results of complete structure and polythermal phase diagram of ternary systems and the isothermal phase diagrams of quaternary system excellently match with the experimental data.

Thermodynamic Modeling of Wax Precipitation in Crude Oils

Most of the crude oils contain waxes which precipitate when temperature drops, resulting in deposition in pipelines and production equipment. It is necessary to set up a model which can predict the wax appearance tem-perature and the amount of solid precipitated in the different conditions. A modified thermodynamic solid-liquid equilibrium model to calculate wax precipitation in crude oil systems has been developed recently. The assumption that precipitated waxes consist of several solid phases is adopted in this research, and the solid-solid transition is also considered in the modified model. The properties of the pseudo-components are determined by using empirical correlations. New correlations for properties of solid-solid and solid-liquid transitions are also established in this work on the basis of the data from the literature. The results predicted by the proposed model for three crude oil systems are compared with the experimental data and the calculated results from the literature, and good agreement is observed.

Thermodynamic modeling of the La-Mg-Y system and Mg-based alloys database

As an example of the La-Mg-Y system, the method how to set up the themaodynamic model of individual phases was introduced in the process of thermodynamic optimization. The solution phases （liquid, body-centered cubic, face-centered cubic, hexagonal close-packed and double hexagonal close-packed） were modeled with the Redlich-Kister equation. The compound energy model has been used to describe the thermodynamic functions of the intermetallic compounds in the La-Mg-Y systems. The compounds Mg2Y, Mg24Y5, Mg12La, Mg17La2, Mg41Las, Mg3La and Mg2La in the La-Mg-Y system were treated as the formulae （Mg,Y）2（La,Mg,Y）, Mg24（La,Mg,Y）4Y, Mg12（La, Y）, Mg17（La,Y）2, Mg41（La,Y）5, Mg3（La,Mg,Y） and Mg2（La, Y）, respectively. A model （La, Mg,Y）0.5（La,Mg,Y）0.5 was applied to describe the compound MgM formed by MgLa and MgY in order to cope with the order-disorder transition between body-centered cubic solution （A2） and MgM with CsCl-type structure （B2） in the La-Mg-Y system. The Gibbs energies of individual phases were optimized in the La-Mg, La-Y and La-Mg-Y systems by CALPHAD technique. The projection of the liquidus surfaces for the La-Mg-Y system was predicted. The Mg-based alloys database including 36 binary and 15 ternary systems formed by Mg, Al, Cu, Ni, Mn, Zn and rare earth elements was set up in SGTE standard.

Inter-comparisons of thermodynamic sea-ice modeling results using various parameterizations of radiative flux

Radiative fluxes are of primary importance in the energy and mass balance of the sea-ice cover. Various parameterizations of the radiative fluxes are studied in a thermodynamic sea-ice model. Model outputs of the surface radiative and heat fluxes and mass balance are compared with observations. The contribution of short-wave radiation is limited to a long part of winter. Therefore, simple schemes are often sufficient. Errors in estimations of the short-wave radiation are due mainly to cloud effects and occasionally to multi-reflection between surface and ice crystals in the air. The long-wave radiation plays an important role in the ice surface heat and mass balance during most part of a winter. The effect of clouds on the accuracy of the simple radiative schemes is critical, which needs further attention. In general, the accuracy of an ice model depends on that of the radiative fluxes.

Thermodynamic modeling and elemental migration for the early stage of rodingitization:An example from the Xialu massif of the Xigaze ophiolite,southern Tibet

The analysis of early stage rodingite from the ultramafic rocks of the Xialu Massif in the Xigaze Ophiolite,Tibet,in China shows that the rodingitization involved continuous changes in fluid composition during different stages of subduction.The early stage prehnite-bearing rodingite was produced at low pressures and temperatures along extensional fractures.Samples of rodingite were collected along a profile from the center to the margin of a rodingitized intrusive igneous rock(^10 m×30 m),and they record wide variations in bulk composition,mineralogy,and texture.The mineral assemblages,from center to margin,vary from(1)relics of primary clinopyroxene(Cpx_(r))and primary amphibole(Amp_(r))+newly formed late amphibole(Act)+primary plagioclase(Pl_(r))+clinozoisite+prehnite+albite+chlorite+titanite+ilmenite(R1 rodingite),through(2)relics of primary clinopyroxene(Cpx_(r))+newly formed late clinopyroxene(Cpx_(n))+primary and late amphiboles(Amp_(r)+Act)+clinozoisite+prehnite+albite+chlorite+titanite(R2 rodingite),to(3)newly formed late clinopyroxene(Cpx_(n))and amphibole(Act)+clinozoisite+prehnite+albite+chlorite+titanite(R3 rodingite).As a result of the metasomatic process of rodingitization,the content of CaO in the whole rock chemical composition from R1 to R3 increases,SiO_(2) decreases,and Na_(2)O+K_(2)O is almost completely removed.Massbalance diagrams show enrichments in large ion lithophile elements such as Rb,Cs,Ba,and Pb as well as Ni during rodingitization.The central part of the rodingitized intrusion(R1 rodingite)was only slightly affected by metasomatism.On the other hand,the contents of the rare earth elements(REEs),high field strength elements(HFSEs;e.g.Zr,Nb,Ta,Hf,and Y),and some highly compatible elements such as Cr and Sc decreased slightly during rodingitization.Thermodynamic modeling based on equilibrium mineral assemblages indicates that the rodingite of the Xialu Massif formed in an H_(2)O-saturated,CO_(2)-rich environment.The estimated conditions of metamorphism were-281-323℃and 0.4-3.9 kbar,representing the subgreenschist facies.In this environment,REEs and HFSEs were soluble in the fluids and partly removed.Moreover,these prehnite rodingites formed in a progressively reducing and less alkaline environment,as indicated by decreases in f(O_(2))and bulk-rock Fe^(3+)/Fe^(2+) ratios,and the records of fluidΔpH from the center to the margin of the studied rodingitized intrusion.

Dynamic Modeling of Tubular SOFC for Marine Power System

Solid oxide fuel cell (SOFC) has been identified as an effective and clean alternative choice for marine power system.This paper emphasizes on the dynamic modeling of SOFC power system and its performance based upon marine operating circumstance.A SOFC power system model has been provided considering thermodynamic and electrochemical reaction mechanism.Subcomponents of lithium ion battery, power conditioning unit, stack structure and controller are integrated in the model.The dynamic response of the system is identified according to the inertia of its subcomponent and controller.Validation of the whole system simulation at steady state and transit period are presented, concerning the effects of thermo inertia, control strategy and seagoing environment.The simulation results show reasonable accuracy compare with lab test.The models can be used to predict performance of a SOFC power system and identify the system response when part of the component parameter is adjusted.

Modeling and optimization of rotary kiln treating EAF dust

Electric arc furnace （EAF） dust from steel industries is listed by the United Sates EPA as a hazardous waste under the regulations of the Resource Conservation and Recovery Act due to the presence of lead, cadmium and chlorine. The disposal of the approximately 650000 t of EAF dust per year in the U.S. and Canada is an expensive and unresolved problem for the majority of steel companies. The Waelz process has been considered as the best process for treating the EAF dust. A process model, combined thermodynamic modeling with heat transfer calculations, has been developed to simulate the chemical reactions, mass and heat trans- fer and heat balance in the kiln. The injection of air into the slag and the temperature profile along the kiln have been modeled. The effect of （CaO＋MgO）/SiO2 on the solidus temperature of slag has also been predicted and discussed. Some optimized results have been presented.

Experimental measurement and thermodynamic modeling of binary and ternary solid–liquid phase equilibrium for the systems formed by L-arabinose,D-xylose and water

Solid–liquid phase equilibrium data for binary(L-arabinose–water) and(D-xylose–water) systems at temperatures from(269.85–298.05) K and ternary(L-arabinose–D-xylose–water) system at temperatures of 273.85 K,278.85 K and 284.45 K were measured at atmospheric pressure.The ternary phase diagrams of the systems were constructed on the base of the measured solubility.Two pure solid phases were formed at given temperatures,including pure L-arabinose and pure D-xylose,which were con firmed and determined by the method of Schreinemakers' wet residue.At the same temperature,the crystallization region of L-arabinose was larger than D-xylose's.The acquired solubility data were then correlated using the NRTL model,Wilson model and Xu model.The calculated solubility with the three models agreed well with the experimental values.

Thermal-hydraulic modeling and analysis of hydraulic system by pseudo-bond graph

To increase the efficiency and reliability of the thermodynamics analysis of the hydraulic system, the method based on pseudo-bond graph is introduced. According to the working mechanism of hydraulic components, they can be separated into two categories: capacitive components and resistive components. Then, the thermal-hydraulic pseudo-bond graphs of capacitive C element and resistance R element were developed, based on the conservation of mass and energy. Subsequently, the connection rule for the pseudo-bond graph elements and the method to construct the complete thermal-hydraulic system model were proposed. On the basis of heat transfer analysis of a typical hydraulic circuit containing a piston pump, the lumped parameter mathematical model of the system was given. The good agreement between the simulation results and experimental data demonstrates the validity of the modeling method.

Mathematical modeling and analysis of thermodynamic processes in a twin-rotor piston engine

In order to study the major performance indicators of the twin-rotor piston engine(TRPE), Matlab/simulink was used to simulate the mathematical models of its thermodynamic processes. With consideration of the characteristics of the working processes in the TRPE, corresponding differential equations were established and then simplified by period features of the TRPE. Finally, the major boundary conditions were figured out. The changing trends of mass, pressure and temperature of working fuel in the working chamber during a complete engine cycle were presented. The simulation results are consistent with the trends of an actual working cycle in the TRPE, which indicates that the method of simulation is feasible. As the pressure in the working chamber is calculated, all the performance parameters of the TRPE can be obtained. The major performance indicators, such as the indicated mean effective pressure, power to weight ratio and the volume power, are also acquired. Compared with three different types of conventional engines, the TRPE has a bigger utilization ratio of cylinder volume, a higher power to weight ratio and a more compact structure. This indicates that TRPE is superior to conventional engines.

Measurement and Modeling for the Solubility of Hydrogen Sulfide in Primene JM-T

The primary aliphatic amine Primene JM-T was investigated as a potential absorbent for H2S removal. The solubility of HzS in JM-T was measured at 298, 313,333,353, and 368 K with H2S partial pressures from 20 to 760 kPa and HzS loading from 0.02 to 0.8 mol H2S per mol JM-T. Relevant physical properties, such as density, viscoslty and dielectric constant, ot the matenal were measured. ＇The thermodynamlc model with two-suttlX Margules equation was used to correlate the experimental vapor-liquid equilibrium data. Furthermore, the absorption mechanism in non-aqueous system is suggested and the difference between non-aqueous and aqueous absorption system is pointed out.

Thermodynamic modeling of gas solubility in aqueous sodium chloride solution

An electrolyte Equation of State is presented by combining the Cubic Plus Association Equation of State,Mean Spherical Approximation and the Born equation.This new model uses experimental relative static permittivity,intend to predict well the activity coefficients of individual ions(ACI)and liquid densities of aqueous solutions.This new model is applied to model water+Na Cl binary system and water+gas+Na Cl ternary systems.The cation/anion-water interaction parameters of are obtained by fitting the experimental data of ACI,mean ionic activity coefficients(MIAC)and liquid densities of water+Na Cl binary system.The cation/anion-gas interaction parameters are obtained by fitting the experimental data of gas solubilities in aqueous Na Cl solutions.The modeling results show that this new model can correlate well with the phase equilibrium and volumetric properties.Without gas,predictions for ACI,MIAC,and liquid densities present relative average deviations of 1.3%,3.6%and 1.4%compared to experimental reference values.For most gas-containing systems,predictions for gas solubilities present relative average deviations lower than 7.0%.Further,the contributions of ACI,and salting effects of Na Cl on gases are analyzed and discussed.

A Maxwell Demon Model Connecting Information and Thermodynamics

In the past decades several theoretical Maxwell＇s demon models have been proposed to exhibit effects such as refrigerating, doing work at the cost of information, and some experiments have been carried out to realize these effects. We propose a model with a two-level demon, information represented by a sequence of bits, and two heat reservoirs. The reservoir that the demon is interacting with depends on the bit. When the temperature difference between the two heat reservoirs is large enough, the information can be erased. On the other hand, when the information is pure enough, heat transfer from one reservoir to the other can happen, resulting in the effect of refrigeration. Genuine examples of such a system are discussed.

The extraction of aromatics using N-methylpyrrolidone: Liquid-liquidequilibrium determination and mechanism exploration

LLE data of cyclooctane/3-methylpentane + benzene/toluene + N-methylpyrrolidone (NMP) at 298.15 Kand 313.15 K under a pressure of 101.3 kPa were measured in this work. The Othmer-Tobias and Handequations were adopted to validate the reliability of LLE data, where the correlation coefficients (R2) wereclose to unity, indicating the high reliability of the experimental data. The experimental data were analyzed using the distribution coefficient (D) and separation factor (S), and the effect of NMP extracting benzene and toluene from aromatics was explored. Meanwhile, the reason for the different extractionefficiencies of benzene and toluene using NMP was analyzed by quantum chemical calculations. TheNRTL and UNIQUAC thermodynamic models were used to correlate the liquid–liquid equilibrium data,and the relevant binary interaction parameters were obtained. The calculated root mean square deviation(RMSD) were all less than 0.0063, indicating that the obtained binary interaction parameters can be usedto simulate and calculate the extraction of aromatics using NMP.