Filtered density function for large eddy simulation of turbulent reacting flows

A methodology termed the “filtered density function” (FDF) is developed and implemented for large eddy simulation (LES) of chemically reacting turbulent flows. In this methodology, the effects of the unresolved scalar fluctuations are taken into account by considering the probability density function (PDF) of subgrid scale (SGS) scalar quantities. A transport equation is derived for the FDF in which the effect of chemical reactions appears in a closed form. The influences of scalar mixing and convection within the subgrid are modeled. The FDF transport equation is solved numerically via a Lagrangian Monte Carlo scheme in which the solutions of the equivalent stochastic differential equations (SDEs) are obtained. These solutions preserve the Itô-Gikhman nature of the SDEs. The consistency of the FDF approach, the convergence of its Monte Carlo solution and the performance of the closures employed in the FDF transport equation are assessed by comparisons with results obtained by direct numerical simulation (DNS) and by conventional LES procedures in which the first two SGS scalar moments are obtained by a finite difference method (LES-FD). These comparative assessments are conducted by implementations of all three schemes (FDF, DNS and LES-FD) in a temporally developing mixing layer and a spatially developing planar jet under both non-reacting and reacting conditions. In non-reacting flows, the Monte Carlo solution of the FDF yields results similar to those via LES-FD. The advantage of the FDF is demonstrated by its use in reacting flows. In the absence of a closure for the SGS scalar fluctuations, the LES-FD results are significantly different from those based on DNS. The FDF results show a much closer agreement with filtered DNS results. © 1998 American Institute of Physics.

An improved method for calculating activities from distribution equilibria

The method of Gibbs-Duhem integration suggested by Speiser et al. has been modified to derive activities from distribution equilibria. It is shown that, in general, the activities of components in melts with a common anion can be calculated, without using their standard Gibbs energies of formation, from eqUilibrium ratios and the knowledge of activities in the metal phase. Moreover, if systems are so chosen that the concentration of one element in the metal phase lies in the Henry's law region (less than 1 %), information on activities in the metal phase is not required. Conversely, activities of elements in an alloy can be readily calculated from equilibrium distribution ratios alone, if the salt phase in equilibrium contains very small amounts of one element. Application of the method is illustrated using distribution ratios from the literature on AgCI-CuCI, AgBr-CuBr, and CuDo.5 -PbD systems. The results indicate that covalent bonding and van der Waals repulsive interactions in certain types of fused salt melts can significantly affect the thermodynamic properties of mixing.

Solubilities of Pb+2, Bi+3, Sb+3 and Cu+1 sulphides in (NaCl KCl Na2S) melts

Solubilities of common metal sulfides have been determined in the (NaCl+KCl) eutectic melt with and without Na2S. A novel gas-phase equilibrium technique has been used for PbS, Bi2S3, and So2S3, and an improved liquid phase equilibrium technique for Cu2S, which eliminates the errors due to physical entrapment of the sulfide phase and segregation on quenching, enabling precise measurements to be made. Solubilities in the (NaCl+KCl) eutectic melt were determined as a function of temperature in the rante 700° to 950°C, and were found to be small. The partial molar heats of mixing of the sulfides in the eutectic melt have been calculated from the solubility measurements, to be 13.3, 31.4, 37.1, and 49.0 kcal for PbSs), Sb2S2(l), and Cu2S(s), respectively. Sodium sulfide addition was observed to enhance these solubilities, the effect being largest for Cu2S followed by Sb2S3, Bi2S3, and PbS. This effect is explained qualitatively. It was observed that PbS and Sb2S3 obey Henry's law up to saturation in (NaCl+KCl+Na2S) melts.

Thermodynamics of solid solutions of silver with indium and tin

A solid oxide galvanic cell and a gas-solid (View the MathML source) equilibration technique have been used to measure the activities of the solutes in the α-solid solutions of silver with indium and tin. The results are consistent with the information now available for the corresponding liquid alloys, the phase diagram and the heats of mixing of the solid alloy. When the results of this study are taken together with published data for the α-solid solutions in Ag + Cd system, it is found that the variation of the excess partial free energy of the solute with mole fraction can be correlated to the electron/atom ratio. The significant thennodynamic parameter that explains the Hume-Rothery findings in these alloys appears to be the rate of change of the excess partial free energy with composition near the phase boundary, and this in turn reflects the value of the solute-solute interaction energy.

Evidence of residual entropy in the cubic spinel Zn2TiO4

The standard free energies of formation of Zn2Ti04 and ZnTi03 have been determined in the temperature range 930° to i ioo'x from electromotive force measurements on reversible solid oxide galvanic cells;Ag-5at%znll I Pt, + CaO-Zr02 ZnO I II Ag-5at%Zn Y20r Th02 CaO-Zr02 + ,Pt Zn2Ti04+ ZnTi03 and II Ag-5at%Zn CaO-Zr02 + ,Pt ZnTi03+ Ti02 The values may be expressed by the equations,2ZnO (wurtz) + Ti02(rut) -> Zn2Ti04(sp), f:!:.Go = -750-2-46T (±75)cal;ZnO(wurtz) +Ti02(rut) -> ZnTi03(ilmen) ,f:!:.Co = -]600-0·]99T(±50)cal.Combination of the free energy values with the calorimetric heat of formation, and low-temperature and high-temperature heat capacity of Zn2Ti04 reported in literature, suggests a residual entropy of ],9 (±0·6) cal K-1 mol ? for the cubic spinel. Ideal mixing of Zn2+ and Ti4+ ions on the octahedral sites would result in a configurational contribution to the entropy of 2· 75 cal K-1 rnol ".The difference is indicative of short-range ordering of cations on octahedral sites.

Activities and their relation to cation distribution in NiAl2O4 MgAl2O4 spinel solid solutions

The activity of NiAl2O4 in NiAl2O4MgAl2O4 solid solutions has been measured by using a solid oxide galvanic cell of the type, Pt, Ni + NiAl2O4 + Al2O3(α)/CaOZrO2/Ni + NixMg1−xAl2O4 + Al2O3(α). Pt, in the temperature range 750–1150°C. The activities in the spinel solid solutions show negative deviations from Raoult's law. The cation distribution in the solid solutions has been calculated using site preference energies independent of composition for Ni2+, Mg2+, and Al3+ ions obtained from crystal field theory and measured cation disorder in pure NiAl2O4 and MgAl2O4, and assumi g ideal mixing of cations on the tetrahedral and octahedral positions. The calculated values correctly predict the decrease in the fraction, α, of Ni2+ ions on tetrahedral sites for 1>x>0.25, observed by Porta et al. [J. Solid State Chem.11, 135 (1974)] but do not support their tentative evidence for an increase in α for x < 0.25. The measured excess free energy of mixing can be completely accounted for by using either the calculated or the measured cation distributions. This suggests that the Madelung energy is approximately a linear function of composition in the solid solutions. The composition of NiOMgO solid solutions in equilibrium with NiAl2O4MgAl2O4 solid solutions has been calculated from the results and information available in literature.

Activities in the spinel solid solution, phase equilibria and thermodynamic properties of ternary phases in the system Cu Fe O

A review of the structural and thermodynamic information and phase equilibria in the Cu-Fe-O system suggested that a consistent, quantitative description of the system is hampered by lack of data on activities in the spinel solid solution CuFe2O4-Fe3O4. Therefore the activity of Fe3O4 in this solid solution is derived from measurements of the oxygen potentials established at 1000°C by mixtures containing Fe2O3 and spinel solid solutions of known composition. The oxygen pressures were measured manometrically for solid solutions rich in CuFe2O4, while for Fe3O4-rich compositions the oxygen potentials were obtained by an emf technique. The activities show significant negative deviations from Raoult’s law. The compositions of the spinel solid solutions in equilibrium with CuO + CuFeO2 and Cu + CuFeO2 were obtained from chemical analysis of the solid solution after magnetic separation. The oxygen potential of the three-phase mixture Cu + CuFeO2 + Fe3O4(spinel s.s.) was determined by a solid oxide galvanic cell. From these measurements a complete phase diagram and consistent thermodynamic data on the ternary condensed phases, CuFeO2 and CuFeO2O4, were obtained. An analysis of the free energy of mixing of the spinel solid solution furnished information on the distribution of cations and their valencies between the tetrahedral and octahedral sites of the spinel lattice, which is consistent with X-ray diffraction, magnetic and Seebeck coefficient measurements.

Electrochemical determination of activities in Cr2O3 Al2O3 solid solutions

The activity of Cr20~ in Cr20~-A12Oa solid solution has been determined in the temperature range 800~176 from electromotive force measurements on the solid oxide galvanic cell Pt,Cr + Cr2OJY~O~-ThO2/Cr + Cr~A12-xO~,Pt The activities of Cr203 and A120~ in the solid solution show both positive and negative deviations from Raoult's law. The heat and entropy of mixing of the solid Solution obtained from the temperature dependence of the emf can be expressed as AH = XCr203XA1203 [31,700Xcrzo3 -}- 37,470XA1203] J mole -I hS = -- 1.8R [Xcr2o3 In Xcr2o3 + XA12o3 In XAaos]The entropy of mixing is 10% lower than that predicted by the Temkin model.The large positive heat of mixing in the Cr2Os-A12Oa solid solution, however, suggests that this apparent: entropy discrepancy originates with the clustering of positive ions on the cation sublattice. The asymmetric miscibility gap exhibited in the CrzOa-A12Oa system below 900~ is consistent with the thermodynamic data trends recorded at the more elevated temperatures.

Thermodynamics of oxide sulfate melts: The system PbO PbSO4

Activities in the PbO-PbSO4 melts at 1253 K have been measured by emf and gas-equilibration techniques. The activity of PbO was directly obtained from the emf of the solid oxide cell, Pt, Ni-NiO/CaO-ZrO2/Auo.92PbO.08, PbOx-PbSO4(1-x), Ir, Pt for 1.0 >XPbO > 0.6. The melt and the alloy were contained in closed zirconia crucibles. Since the partial pressure of SO2 gas in equilibrium with the melt and alloy was appreciable (>0.08 atm) atXPbO < 0.6, activities at lower PbO concentrations were derived from measurements of the weight gain of pure PbO under controlled gas streans of Ar + SO2 + O2. The partial and integral free energies of mixing at 1253 K were calculated and found to fit a subregular model: ΔGEPbO =X2PbSO4 {-42,450 + 20,000X2PbSO4} J mol-1 ΔGEPbO =X2pbSO {-12,450 - 20,000XPbS} J mol-1 ΔGEpbSOXPbSO4 {-32,450XPbS - 22,450XPbSO4 } J mol-1. The standard free energy of formation of liquid PbSO4 from pure liquid PbO and gaseous SO3 at 1 atm at 1253 K was evaluated as -88.02 (±0.72) kJ mol-1.

A comparison of the structure of liquid alloys derived from thermodynamic and diffraction studies

A significant amount of research on the thermodynamic properties of molten alloys is undertaken for obtaining insights into their structure . The partial and integral molar enthalpies, entropies and volumes of mixing provide some general information on the nature and strength of atomic bonds and the distribution of atoms. However, until recently it has been difficult to derive specific quantitative information because the excess entropy of mixing contains configurational , vibrational , electronic , and sometimes magnetic contributions which cannot be easily separated.

Correlations between structure and thermodynamic properties of spinels and their applications in technology

Measurements a/the Gibbs' energy enthalpy and entrupy vffarmation oj chromites, vanadites and alumlnat.:s 0/ F", Ni. Co'. Mn, Zn Mg and Cd, using solid oxide galvanic cells over a ternperature range extending approximately lOOO°C, have shown that the '~'Ilir"!,,, J'JrIl/iJ~ tion 0/ cubic 2-3 oxide spinel phases (MX!O,), from component oxide (MO) with rock-salt and X.Os whir c(1f'l/!ldwn st!'llt'lw,·. call b,' represented by a semi-empirical correlalion, ~S~ = --LiS + L'i,SM +~S~:"d(±O.3) cal.deg-1 mol-1 where /',.SM Is the entropy 0/calian mixing oillhe tetrahedral alld octahedral sites o/the spinel and Sr:~ is tlie enfropy associaf,'d Wifh Ih,' randomization a/the lahn-Telier distortions. A review a/the methods/or evaluating the cation distriblltion lfl spille!s suggeJ{j' l/r,l! Ihe most promising scheme is based Oil octahedral site preference energies from the crystal field theory for the Iral1silioll IIIl'f"! IlIIL';. For I/""-Irallsifioll melal cal ions site preference energies are derived relative /0 thol'lt fLI, [ransilion metal ions from measured high tClllP('ftJi ure Cal iUlI disll iiJuriol1 in spine! phases thar contail! one IransilioJl metal and another non-transition metal carion. For 2-3 srinds compulatiorrs b,IS"J Oil i.!c[J;' Temkin mixing on each catioll subialtice predici JistributionJ that are In fair agreement with X-ray and 1I1'IIIrOll ditTraction, /IIdg""!ic dll.! electrical propcrries, and spectroscopic measurements. In 2-4 spineis mixing vI ions do not foliow strictly ideal slllIistli:al Jaws, Th,' OIl/up) associated with the randomizalion 0/the Jllhn-Teller dislOriioll" appear to be significant, only ill spinels witll 3d'. 3d', 3d' (ifld~UI' iOtls in tetrahedral and 3d' and 3d9 ions in octahedral positions. Application 0/this structural model for predicting the thermodynamic proputies ofspinel solid .,olutiofl5 or,' illustrated. F,lr complex systems additional contributions arising from strain fields, redox equilibria and off-center ions have to be qllalllififti. The entropy correlation for spinels provides a method for evaluating structure tran:.jormafiofl entropies in silllple o.\id.-s, ["founlllion on the relative stabilities ofoxides in different crystallCtructures is USe/III for computer ea/culaliof! a/phase dfugrullls ofIlIrer,',,1 III (N.lll1ie5 by method, similar to thost: used by Kaufman and Bernstein for refractory alloy systems. Examples oftechnoiogical appliCation tnclude the predictioll ofdeoxidation equilibria in Fe-Mn-AI-O s),slelll at 1600°C duj ,'Ulllpltfalion 0/phase relutions in Fe-Ni-Cr-S system,

A derivation of thermodynamic quantities of liquid alloys from their structure factors

A method of deriving the thermodynamic properties of mixing in liquid alloys Delta G, Delta S and Delta H, from low-Q scattering data has been presented. As an example, the method has been demonstrated with liquid Na-Ga alloys for which both thermodynamic and diffraction data have recently been obtained by the authors.

Computation of thermodynamic properties of liquid ferrous alloys from structural measurements

Relation between X-ray scattering intensities, mean square thermal fluctuations and thermodynamic properties. High temperature X-ray diffraction study of liquid Fe-Ni and Fe-Si alloys using reflection and transmission geometries. Calculation of the structure factor as a function of wave vector. Extrapolation to zero wave vector. Calculation of the concentration-concentration correlation function defined by A. B. Bhatia and D. E. Thorton. Computation of thermodynamic quantities of mixing A G, LlH and LlS for the binary alloys. Comparison with direct thermodynamic measurements reported in the literature.

Revision of thermodynamic data on MnO Al2O3 melts

The high temperature region of the MnO-A1203 phase diagram has been redetermined to resolve some discrepancies reported in the literature regarding the melting behaviour of MnA1,04. This spinel was found to melt congruently at 2108 (+ 15) K. Theactivity of MnOin MnO-Al,03 meltsand in the two phase regions, melt + MnAI,04 and MnAI2O4 + A1203, has been determined by measuring the manganese concentration in platinum foils in equilibrium under controlled oxygen potentials. The activity of MnO obtained in this study for M ~ O - A I ,m~el~ts is in fair agreement with the results of Sharma and Richardson.However. the alumina-rich melt is found to be in equilibrium with MnAl,04 rather than AI2O3. as suggested by ~ha rmaan d Richardson. The value for the acthity of MnO in the M~AI ,O,+ A1,03 two phaseregion permits a rigorous application of the Gibbs-Duhem equation for calculating the activity of A1203 and the integral Gibbs' energy of mixing of MnO-A1203 melts, which are significantly different from those reported in the literature.

Thermodynamic properties of Fe3O4 FeAl2O4 spinel solid solution

The Gibbs energy of mixing for the system Fe3O4-FeAl2O4 was determined at 1573 K using a gas-metal-oxide equilibration technique. Oxide solid solution samples were equilibrated with Pt foils under controlled CO+CO2 gas streams. The equilibrium iron concentration in the foil was determined by chemical analysis. The cation distribution between tetrahedral and octahedral sites in the spinel crystal can be calculated from site-preference energies and used as an alternate method of determining some thermodynamic properties, including the Gibbs energy of mixing. The solvus occurring at low temperatures in the system Fe3C4-FeAl2C4 was used to derive the effect of lattice distortion due to cation size difference on the enthalpy of mixing and to obtain a better approximation to the measured thermodynamic quantities.