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.

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.

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.

Seasonal Prediction of Tropical Cyclones and Storms over the Southwestern Indian Ocean Region Using the Generalized Linear Models

Tropical cyclones (TCs) and storms (TSs) are among the devastating events in the world and southwestern Indian Ocean (SWIO) in particular. The seasonal forecasting TCs and TSs for December to March (DJFM) and November to May (NM) over SWIO were conducted. Dynamic parameters including vertical wind shear, mean zonal steering wind and vorticity at 850 mb were derived from NOAA (NCEP-NCAR) reanalysis 1 wind fields. Thermodynamic parameters including monthly and daily mean Sea Surface Temperature (SST), Outgoing Longwave Radiation (OLR) and equatorial Standard Oscillation Index (SOI) were used. Three types of Poison regression models (i.e. dynamic, thermodynamic and combined models) were developed and validated using the Leave One Out Cross Validation (LOOCV). Moreover, 2 × 2 square matrix contingency tables for model verification were used. The results revealed that, the observed and cross validated DJFM and NM TCs and TSs strongly correlated with each other (p ≤ 0.02) for all model types, with correlations (r) ranging from 0.62 - 0.86 for TCs and 0.52 - 0.87 for TSs, indicating great association between these variables. Assessment of the model skill for all model types of DJFM and NM TCs and TSs frequency revealed high skill scores ranging from 38% - 70% for TCs and 26% - 72% for TSs frequency, respectively. Moreover, results indicated that the dynamic and combined models had higher skill scores than the thermodynamic models. The DJFM and NM selected predictors explained the TCs and TSs variability by the range of 0.45 - 0.65 and 0.37 - 0.66, respectively. However, verification analysis revealed that all models were adequate for predicting the seasonal TCs and TSs, with high bias values ranging from 0.85 - 0.94. Conclusively, the study calls for more studies in TCs and TSs frequency and strengths for enhancing the performance of the March to May (MAM) and December to October (OND) seasonal rainfalls in the East African (EA) and Tanzania in particular.

Calculation of Vapour Pressure of Metals by Statistical-Mechanical Method With the Debye Model

Statistical expression of vapour pressure equations of metals is derived from the Debye model.The statistical distribution of T_(-p) ensemble is presented in an in-elab- orate mode and the partition function is defined.The vapour pressure of eleven metals have been calculated with the Debye equation and compared with those given by the E- instein equation and empirical equation.Comparison of results of calculation from dif- ferent methods show their evident accordance within the same orders of magnitude.

Local structure of calcium silicate melts from classical molecular dynamics simulation and a newly constructed thermodynamic model

The distributions of local structural units of calcium silicate melts were quantified by means of classical molecular dynamics simulation and a newly constructed structural thermodynamic model. The distribution of five kinds of Si-O tetrahedra Qi from these two methods was compared with each other and also with the experimental Raman spectra, an excellent agreement was achieved. These not only displayed the panorama distribution of microstructural units in the whole composition range, but also proved that the thermodynamic model is suitable for the utilization as the subsequent application model of spectral experiments for the thermodynamic calculation. Meanwhile, the five refined regions mastered by different disproportionating reactions were obtained. Finally, the distributions of two kinds of connections between Qi were obtained, denoted as Qi-Ca-Qj and Qi-[Ob]-Qj, from the thermodynamic model, and a theoretical verification was given that the dominant connections for any composition are equivalent connections.

A universal thermodynamic model of calculating mass action concentrations for structural units or ion couples in aqueous solutions and its applications in binary and ternary aqueous solutions

A universal thermodynamic model of calculating mass action concentrations for structural units or ion couples in ternary and binary strong electrolyte aqueous solution was developed based on the ion and molecule coexistence theory and verified in four kinds of binary aqueous solutions and two kinds of ternary aqueous solutions. The calculated mass action concentrations of structural units or ion couples in four binary aqueous solutions and two ternary solutions at 298.15 K have good agreement with the reported activity data from literatures after shifting the standard state and concentration unit. Therefore, the calculated mass action concentrations of structural units or ion couples from the developed universal thermodynamic model for ternary and binary aqueous solutions can be applied to predict reaction ability of components in ternary and binary strong electrolyte aqueous solutions. It is also proved that the assumptions applied in the developed thermodynamic model are correct and reasonable, i.e., strong electrolyte aqueous solution is composed of cations and anions as simple ions, H2O as simple molecule and other hydrous salt compounds as complex molecules. The calculated mass action concentrations of structural units or ion couples in ternary and binary strong electrolyte aqueous solutions strictly follow the mass action law.

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.

Thermodynamic model of lead oxide activity in PbO-CaO-SiO_2-FeO-Fe_2O_3 slag system

According to the ion and molecule coexistence theory, a thermodynamic model of lead oxide activity in PbO-CaO-SiO2-FeO-Fe2O3 slag system was established at the temperature of 1273-1733 K. The activities of Pb O in slag were calculated, and their equal activity curves were plotted. The influences of slag basicity Q, iron oxide rate R and temperature T on activity NPb O and activity coefficient γPbO were also investigated. Results show that the calculated values of γPb O are in good agreement with the reported experimental data, showing that the model can wholly embody the slag structural characteristics. NPbO departures positively from Raoult values, and increases with increasing Pb O content in slag but changes little with T. γPbO increases with increasing Q, and goes through the maximum with increasing R for basic slag（Q0.3）. Results can be applied to the thermodynamic research and operational optimization of modern lead smelting technologies.

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.

A general thermodynamic model for the long-period stacking ordered phases in magnesium alloys

A thermodynamic model Mg x(Xs,Mg)6(Xl,Mg)8(Xs and Xl are elements smaller and larger than Mg)for long-period stacking ordered phases(LPSOs)was proposed based on two key factors:the Xs 6 Xl 8-type L12 clusters and the variation of chemical compositions.In general,all available LPSOs can be described with this model.As a representative system,Mg-Y-Zn with three LPSOs was investigated using the CALPHAD(calculation of phase diagram)approach aided with first-principles calculations.Two new three-phase equilibria were predicted and were validated by key experiments.The model-based descriptions will be the basis for the research and development of magnesium alloys.

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.

A novel hybrid thermodynamic model for pore size distribution characterisation for shale

Scholars often see the gas adsorption technique as a straight-to-interpret technique and adopt the pore size distribution(PSD)given by the gas adsorption technique directly to interpret pore-structure-related issues.The oversimplification of interpreting shale PSD based on monogeometric thermodynamic models leads to apparent bias to the realistic pore network.This work aims at establishing a novel thermodynamic model for shale PSD interpretation.We simplified the pore space into two geometric types—cylinder-shaped and slit-shaped.Firstly,Low-temperature Nitrogen Adsorption data were analyzed utilizing two monogeometric models(cylindrical and slit)to generate PSD_(cyl).and PSD_(slit);Secondly,pore geometric segmentation was carried out using Watershed by flooding on typical SEM images to obtain the ratio of slit-shaped(∅_(s))and cylinder-shaped pores(∅_(c)).Combining the results of the two,we proposed a novel hybrid model.We performed pyrolysis,XRD,FE-SEM observation,quantitative comparison with the results obtained by the DFT model,and fractal analysis to discuss the validity of the obtained PSD_(Hybrid).The results showed that:the hybrid model proposed in this work could better reflect the real geometry of pore space and provide a more realistic PSD;compared with thermodynamic monogeometric models,PSD obtained from the hybrid model are closer to that from the DFT model,with an improvement in the deviation from the DFT model from 5.06%to 68.88%.The proposed hybrid model has essential application prospects for better interpretation of shale pore space.It is also worth noting that we suggest applying the proposed hybrid model for PSD analysis in the range of 5-100 nm.

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.

Design and Numerical Simulation of an Arctic Ocean Circulation and Thermodynamic Sea－Ice Model

In this paper, the first version of a new Arctic Ocean circulation and thermodynamic sea-ice model is presentedby the authors based on the framework of a twenty-layer World Oceanic general circulation model developed byZhang et al. in 1994. The model's domain covers the Arctic Ocean and Greenland-Norwegian Seas with the horizontal resolution of 200 km×200 km on a stereographic projection plane. In vertical, the model uses the Eta-coordinate(Sigma modified to have quasi-horizontal coordinate surfaces) and has ten unevenly-spaced layers to cover the deepest water column of 3000 m. Two 150-year integrations of coupling the ocean circulation model with the sea-icemodel have been performed with seasonally cyclic surface boundary conditions. The only difference between the tWoexperiments is in the model's geography. Some preliminary analyses of the experimental results have been done focused on the following aspects: (1) surface layer temperature, salinity and current; (2) the' Atlantic Layer'; (3)sea-ice cover and its seasonal variation. In comparison with the available observational data, these results are acceptable with reasonable accuracy.

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.

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.

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.

A Molecular Thermodynamic Model for Interfacial Tension in Surfactant-Oil-Water System

An interfacial equation of state based on perturbation theory for surfactant-oil-water system has been developed. By combining the interfacial equation of state with Boudh-Hir and Mansoori's model, a molecular ther-modynamic model has been proposed. The interfacial tension of surfactant-oil-water systems can be calculated from the surface tensions of pure oil and water by this model. The interfacial tension data for sodium dodecyl sulphate-heptane-water system, polyoxyethylene n-octylphenol-heptane-water system and hexadecyl trimethyl ammonium bromide-heptane-water system have been correlated. By using the adjustable parameters obtained, the interfacial tensions of these systems at other temperatures have been predicted. Both the correlated and the predicted values are satisfactory.

Simplified Thermodynamic Model for Pro-Eutectoid Ferrite Formation in Multicomponent Structural Steel

By introducing aparameter of difference in ferrite formation temperature between binary Fe-C and multicomponent system,and referring to the thermodynamic model for Fe-C binary system,a simplified thermodynamic model for pro-eutectoid ferrite formation in Fe-ΣXiC multicomponent structural steels(Xi=Mn,Si,Mo,Cr,Ni or Ti,etc)was suggested.The comparison of the calculated Ae3 temperatures with the measured data of steels 42 shows that the relative standard deviation and root-mean-square(RMS)error between them are only 0.71% and 8.92 K,respectively.However,the deviations between the same measured data and the values calculated from the superelement model are as high as 1.86% and 23.83 K,respectively.It can be concluded that the simplified thermodynamic model for pro-eutectoid ferrite formation in multicomponent structural steels is acceptable and the calculated Ae3 temperatures are in good agreement with the experimental data.