首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 140 毫秒
1.
Realistic simulations of fluid flow in geologic systems have severely been hampered by the lack of a consistent formulation for fluid properties for binary salt-water fluids over the temperature-pressure-composition ranges encountered in the Earth’s crust. As the first of two companion studies, a set of correlations describing the phase stability relations in the system H2O-NaCl is developed. Pure water is described by the IAPS-84 equation of state. New correlations comprise the vapor pressure of halite and molten NaCl, the NaCl melting curve, the composition of halite-saturated liquid and vapor, the pressure of vapor + liquid + halite coexistence, the temperature-pressure and temperature-composition relations for the critical curve, and the compositions of liquid and vapor on the vapor + liquid coexistence surface. The correlations yield accurate values for temperatures from 0 to 1000 °C, pressures from 0 to 5000 bar, and compositions from 0 to 1 XNaCl (mole fraction of NaCl). To facilitate their use in fluid flow simulations, the correlations are entirely formulated as functions of temperature, pressure and composition.  相似文献   

2.
Most submarine gas hydrates are located within the two-phase equilibrium region of hydrate and interstitial water with pressures (P) ranging from 8 to 60 MPa and temperatures (T) from 275 to 293 K. However, current measurements of solubilities of methane in equilibrium with hydrate in the absence of a vapor phase are limited below 20 MPa and 283.15 K, and the differences among these data are up to 30%. When these data were extrapolated to other P-T conditions, it leads to large and poorly known uncertainties. In this study, in situ Raman spectroscopy was used to measure methane concentrations in pure water in equilibrium with sI (structure one) methane hydrate, in the absence of a vapor phase, at temperatures from 276.6 to 294.6 (±0.3) K and pressures at 10, 20, 30 and 40 (±0.4%) MPa. The relationship among concentration of methane in water in equilibrium with hydrate, in mole fraction [X(CH4)], the temperature in K, and pressure in MPa was derived as: X(CH4) = exp [11.0464 + 0.023267 P − (4886.0 + 8.0158 P)/T]. Both the standard enthalpy and entropy of hydrate dissolution at the studied T-P conditions increase slightly with increasing pressure, ranging from 41.29 to 43.29 kJ/mol and from 0.1272 to 0.1330 kJ/K · mol, respectively. When compared with traditional sampling and analytical methods, the advantages of our method include: (1) the use of in situ Raman signals for methane concentration measurements eliminates possible uncertainty caused by sampling and ex situ analysis, (2) it is simple and efficient, and (3) high-pressure data can be obtained safely.  相似文献   

3.
Soils in the McMurdo Dry Valleys, Antarctica contain ice and considerable amounts of salt. Ice often occurs at shallow depth throughout the Dry Valleys and other areas of hyperarid permafrost, notably on Mars. This common occurrence of shallow ice is enigmatic; however, since according to published sublimation models it should disappear relatively quickly (at rates of order 0.1 mm a−1) due to vapor loss to the atmosphere. This loss may be offset by recharge from snowmelt infiltrating and freezing in the soil. Herein, we present a first quantitative estimate of this recharge based on measured vertical profiles of δD and δ18O that reveal considerable detail about the sources and sinks of ice. We model these profiles, taking into account the salt content and a soil temperature record along a 1.6 m depth profile of ∼10 ka old ice-cemented soils in Victoria Valley, Antarctica. The stable isotopes of ice are enriched in heavy isotopes at the top of the ice cement (20 cm depth); both δD and δ18O values plotted against depth exhibit a concave upward curve. At depth, the isotope composition is similar to that of Lake Victoria and modern meteoric water. The concave shape of the isotope profile is suggestive of downward advection-dispersion of snowmelt water enriched in heavy isotopes into the ice cement. Our advection-dispersion model, coupled with field data, enables us to quantify the advective flux and dispersion of melt water into the ice. The advective velocity and dispersion coefficient depend on the time since advection began and the ice-to-brine ratio; they are, respectively, of the order of 10−11-10−10 m s−1 and 10−12-10−11 m2 s−1. These values suggest that over the ∼10 ka time period, a total of 190 mm water infiltrated into the ice-cemented ground. The isotope composition and deuterium excess values of the uppermost ice cement can be modeled from snowmelt water enriched in salts using open system-Rayleigh fractionation. To develop the isotopic signature of the upper ice cement requires evaporation of ∼95% of the snowmelt water. Based on 190 mm brine infiltrating into the soil requires an initial total of ∼4 m of snowmelt water. This corresponds to ∼0.4 mm a−1 suggesting that, under the current climate condition, water from snowmelt is sufficient to compensate modeled sublimation rates, and therefore conserve ground ice in Victoria Valley.  相似文献   

4.
A commercial (Setaram C80) calorimeter has been modified to measure the heat capacities of highly caustic solutions at temperatures up to 300 °C and pressures up to 20 MPa. The improvements have allowed more accurate determination of the isobaric volumetric heat capacities of chemically aggressive liquids at high temperatures. Test measurements with aqueous solutions of sodium chloride showed a reproducibility of about ±0.1%, with an accuracy of ∼0.3% or better, over the whole temperature range. Heat capacities of aqueous solutions of sodium hydroxide at concentrations from 0.5 to 8 mol/kg were measured at temperatures from 50 to 300 °C and a pressure of 10 MPa. Apparent molar isobaric heat capacities of NaOH(aq) were calculated using densities determined previously for the same solutions by vibrating-tube densimetry. Standard state (infinite dilution) partial molar isobaric heat capacities of NaOH(aq) were obtained by extrapolation using an extended Redlich-Meyer equation. Values of the standard heat capacity change for the ionization of water up to 300 °C were derived by combining the present results with the literature data for HCl(aq) and NaCl(aq).  相似文献   

5.
Preliminary measurements were carried out of the solubility of the O2-buffering assemblage bismuth + bismite (Bi2O3) in aqueous liquid–vapor and vapor-only systems at temperatures of 220, 250 and 300 °C. All experiments were carried out in Ti reaction vessels and were designed such that the Bi solids were contained in a silica tube that prevented contact with liquid water at any time during the experiment. Two blank (no Bi solids present) liquid–vapor experiments at 220 °C yielded Bi concentrations (±1σ) in the condensed liquid of 0.22 ± 0.02 mg/L, whereas the solubility measurements at this temperature yielded an average value of approximately 6 ± 9 mg/L, with replicate experiments ranging from 0.3 to 26 mg/L. Although the 6 mg/L value is associated with a considerable degree of uncertainty, the experiments do indicate transport of Bi through the vapor phase. Measured Bi concentrations in the condensed liquid at 250 °C were in the same range as those at 220 °C, whereas those at 300 °C were significantly lower (i.e., all below the blank value). Vapor-only experiments necessarily contained much smaller initial volumes of water, thereby making the results more susceptible to contamination. Single blank runs at 220 and 300 °C yielded Bi concentrations of 82 and 16 mg/L, respectively. Measured concentrations (±1σ) of Bi in the vapor-only solubility experiments at 220 °C were 235 ± 78 mg/L for an initial water volume of 0.5 mL, and at 300 °C were 56 ± 30 mg/L and 33 ± 21 for initial water volumes of 1 and 2 mL, respectively, suggesting strong preferential partitioning of Bi into the vapor. The results indicate a negative dependence of Bi solubility on temperature, but are inconclusive with respect to the dependence of Bi solubility on water density or fugacity.  相似文献   

6.
This work reports the application of thermodynamic models, including equations of state, to binary (salt-free) CH4-H2O fluid inclusions. A general method is presented to calculate the compositions of CH4-H2O inclusions using the phase volume fractions and dissolution temperatures of CH4 hydrate. To calculate the homogenization pressures and isolines of the CH4-H2O inclusions, an improved activity-fugacity model is developed to predict the vapor-liquid phase equilibrium. The phase equilibrium model can predict methane solubility in the liquid phase and water content in the vapor phase from 273 to 623 K and from 1 to 1000 bar (up to 2000 bar for the liquid phase), within or close to experimental uncertainties. Compared to reliable experimental phase equilibrium data, the average deviation of the water content in the vapor phase and methane solubility in the liquid phase is 4.29% and 3.63%, respectively. In the near-critical region, the predicted composition deviations increase to over 10%. The vapor-liquid phase equilibrium model together with the updated volumetric model of homogenous (single-phase) CH4-H2O fluid mixtures (Mao S., Duan Z., Hu J. and Zhang D. (2010) A model for single-phase PVTx properties of CO2-CH4-C2H6-N2-H2O-NaCl fluid mixtures from 273 to 1273 K and from 1 to 5000 bar. Chem. Geol.275, 148-160), is applied to calculate the isolines, homogenization pressures, homogenization volumes, and isochores at specified homogenization temperatures and compositions. Online calculation is on the website: http://www.geochem-model.org/.  相似文献   

7.
A set of correlations for the volumetric properties and enthalpies of phases in the system H2O-NaCl as a function of temperature, pressure, and composition has been developed that yields accurate values from 0 to 1000 °C, 1 to 5000 bar, and 0 to 1 XNaCl. The volumetric properties of all fluid phases from low-density vapor to hydrous salt melts and single-phase binary fluids at high pressures and temperatures, can be described by a simple equation
  相似文献   

8.
The effect of water on heat capacity has been determined for four series of hydrated synthetic aluminosilicate glasses and supercooled liquids close to albite, phonolite, trachyte, and leucogranite compositions. Heat capacities were measured at atmospheric pressure by differential scanning calorimetry for water contents between 0 and 4.9 wt % from 300 K to about 100 K above the glass transition temperature (Tg). The partial molar heat capacity of water in polymerized aluminosilicate glasses, which can be considered as independent of composition, is (J/mol K). In liquids containing at least 1 wt % H2O, the partial molar heat capacity of water is about 85 J/mol K. From speciation data, the effects of water as hydroxyl groups and as molecular water have tentatively been estimated, with partial molar heat capacities of 153 ± 18 and 41 ± 14 J/mol K, respectively. In all cases, water strongly increases the configurational heat capacity at Tg and exerts a marked depressing effect on Tg, in close agreement with the results of viscosity experiments on the same series of glasses. Consistent with the Adam and Gibbs theory of relaxation processes, the departure of the viscosity of hydrous melts from Arrhenian variations correlates with the magnitude of configurational heat capacities.  相似文献   

9.
The heat capacity (Cp) of two synthetic spessartine samples (Sps) was measured on 20-30 mg-size samples in the temperature range 2-864 K by relaxation calorimetry (RC) and differential scanning calorimetry (DSC). The polycrystalline spessartine samples were synthesized in two different laboratories at high pressures and temperatures from glass and oxide-mixture starting materials and characterized by X-ray powder diffraction and electron-microprobe analysis. The low-temperature heat capacity data show a prominent lambda transition with a peak at 6.2 K, which is interpreted to be the result of a paramagnetic-antiferromagnetic phase transition. The DSC data around ambient T agree excellently with the RC data and can be represented by the Cp polynomial for T > 250 K:
  相似文献   

10.
We perform first principles molecular dynamics simulations of Mg2SiO4 liquid and crystalline forsterite. On compression by a factor of two, we find that the Grüneisen parameter of the liquid increases linearly from 0.6 to 1.2. Comparison of liquid and forsterite equations of state reveals a temperature-dependent density crossover at pressures of ∼12-17 GPa. Along the melting curve, which we calculate by integration of the Clapeyron equation, the density crossover occurs within the forsterite stability field at P = 13 GPa and T = 2550 K. The melting curve obtained from the root mean-square atomic displacement in forsterite using the Lindemann law fails to match experimental or calculated melting curves. We attribute this failure to the liquid structure that differs significantly from that of forsterite, and which changes markedly upon compression, with increases in the degree of polymerization and coordination. The mean Si coordination increases from 4 in the uncompressed system to 6 upon twofold compression. The self-diffusion coefficients increase with temperature and decrease monotonically with pressure, and are well described by the Arrhenian relation. We compare our equation of state to the available highpressure shock wave data for forsterite and wadsleyite. Our theoretical liquid Hugoniot is consistent with partial melting along the forsterite Hugoniot at pressures 150-170 GPa, and complete melting at 170 GPa. The wadsleyite Hugoniot is likely sub-liquidus at the highest experimental pressure to date (200 GPa).  相似文献   

11.
This article presents methods for predicting the standard partial molar Gibbs energy (standard chemical potential) and related derivative properties of aqueous hydroxy and aminoderivatives of (alkyl)benzenes over a wide range of temperatures and pressures. A thorough literature overview was conducted for collecting all available experimental data resulting from phase equilibrium, calorimetric and volumetric measurements that allow calculation of the thermodynamic properties of hydration. New experimental values are presented for solubility in water of isomeric toluidines and for the partial molal volume of phenol and cresols at high temperatures. Building upon the acquired database several prediction schemes were developed and tested for calculating the standard thermodynamic properties (and namely the Gibbs energy of hydration) of aqueous alkylphenols and alkylanilines as a function of temperature and pressure. First, a simple group contribution method was proposed for estimations at 298.15 K and 0.1 MPa using the simultaneous treatment of all available data on hydration properties at near ambient conditions. Second, this group contribution method allowed re-adjustment of the parameters of the Helgeson-Kirkham-Flowers model (HKF) using a new procedure proposed recently by Plyasunov and Shock [Plyasunov, A.V., Shock, E.L., 2001b. Correlation strategy for determining the parameters of the revised Helgeson-Kirkham-Flowers model for aqueous nonelectrolytes. Geochim. Cosmochim. Acta65, 3879-3900]. Third, using the Sedlbauer-O’Connell-Wood equation of state for aqueous species (SOCW), group contributions were determined for predictions at high temperatures and pressures by simultaneous correlation of all available thermodynamic data on hydration properties. The latter method was constrained by the group contributions at 298.15 K and 0.1 MPa making both group contribution schemes consistent at near ambient conditions. The calculations from the HKF and SOCW equations of state and those from the simple thermodynamic integration of the data at 298.15 K and 0.1 MPa were compared for several alkylphenols and alkylanilines. Equilibrium constants for hydration reactions obtained from the three approaches are in very good agreement at temperatures to at least 400 K. At higher temperatures we assess the accuracy of different predictive schemes and their associated uncertainties. The reliable predictions of the standard chemical potentials to at least 573 K and 100 MPa are possible by the group contribution method using the SOCW equation of state.  相似文献   

12.
Ordering of interfacial water at the hematite and corundum (0 0 1)-water interfaces has been characterized using in situ high resolution specular X-ray reflectivity measurements. The hematite (0 0 1) surface was prepared through an annealing process to produce a surface isostructural with corundum (0 0 1), facilitating direct comparison. Interfacial water was found to display a similar structure on this pair of isostructural surfaces. A single layer of adsorbed water having a large vibrational amplitude was present on each surface and additional ordering of water extended at least 1 nm into the bulk fluid, with the degree of ordering decreasing with increasing distance from the surfaces. Consistent with prior studies of the (0 1 2) and (1 1 0) surfaces of hematite and corundum, the configuration of water above the (0 0 1) surfaces is primarily controlled by the surface structure, specifically the arrangement of surface functional groups. However, interfacial water at the (0 0 1) surfaces displayed significantly larger vibrational amplitudes throughout the interfacial region than at other isostructural sets of hematite and corundum surfaces, indicating weaker ordering. Comparison of the vibrational amplitudes of adsorbed water on a series of oxide, silicate, and phosphate mineral surfaces suggests that the presence or absence of a substantial interfacial electrostatic field is the primary control on water ordering and not the surface structure itself. On surfaces for which charge originates dominantly through protonation-deprotonation reactions the controlling factor appears to be whether conditions exist where most functional groups are uncharged as opposed to the net surface charge. The doubly coordinated functional groups on hematite and corundum (0 0 1) surfaces are largely uncharged under slightly acidic to circumneutral pH conditions, leading to weak ordering, whereas singly coordinated groups on (0 1 2) and (1 1 0) surfaces of these phases are always charged, even when the net surface charge is zero, and induce strong water ordering. Surfaces lacking structural charge can thus be divided into two distinct classes that induce either strong or weak ordering of interfacial water. Surface functional group coordination is the ultimate control on this division as it determines the charge state of such groups under different protonation configurations. Ion adsorption and electron transfer processes may differ between these classes of surfaces because of the effect of water ordering strength on interfacial capacitances and hydrogen bonding.  相似文献   

13.
Molecular dynamics simulations of water in contact with the (0 0 1) and (0 1 0) surfaces of orthoclase (KAlSi3O8) were carried out to investigate the structure and dynamics of the feldspar-water interface, contrast the intrinsic structural properties of the two surfaces, and provide a basis for future work on the diffusion of ions and molecules in microscopic mineral fractures. Electron density profiles were computed from the molecular dynamics trajectories and compared with those derived experimentally from high-resolution X-ray reflectivity measurements by Fenter and co-workers [Fenter P., Cheng L., Park C., Zhang H. and Sturchio N. C. (2003a) Structure of the orthoclase (0 0 1)- and (0 1 0)-water interfaces by high-resolution X-ray reflectivity. Geochim. Cosmochim. Acta67, 4267-4275]. For each surface, three scenarios were considered whereby the interfacial species is potassium, water, or a hydronium ion. Excellent agreement was obtained for the (0 0 1) surface when potassium is the predominant interfacial species; however, some discrepancies in the position of the interfacial peaks were obtained for the (0 1 0) surface. The two surfaces showed similarities in the extent of water ordering at the interface, the activation energies for water and potassium desorption, and the adsorption localization of interfacial species. However, there are also important differences between the two surfaces in the coordination of a given adsorbed species, adsorption site densities, and the propensity for water molecules in surface cavities and those in the first hydration layer to coordinate to surface bridging oxygen atoms. These differences may have implications for the extent of dissolution in the low-pH regime since hydrolysis of Si(Al)OSi(Al) bonds is a major dissolution mechanism.  相似文献   

14.
A thermodynamic model and equation of state (EOS) is developed from the molecular dynamics simulation experiments of Spera et al. (2009) for CaAl2Si2O8 liquid over the temperature range 3500-6000 K and pressure interval 0-125 GPa. The model is constructed utilizing the isothermal Universal EOS of Vinet et al. (1986) combined with an expression for the temperature-dependence of the internal energy derived from density functional theory (Rosenfeld and Tarazona, 1998). It is demonstrated that this model is more successful at reproducing the data than the temperature-dependent Universal EOS (Vinet et al., 1987) or the volume-explicit EOS of Ghiorso (2004a). Distinct parameterizations are required to model low (<20 GPa) and high (>20 GPa) pressure regimes. This result is ascribed to the affect of liquid structure on macroscopic thermodynamic properties, specifically the interdependence of average cation-oxygen coordination number on the bulk modulus. The thermodynamic transition between the high- and low-pressure parameterizations is modeled as second order, although the nature of the transition is open to question and may well be first order or lambda-like in character.Analysis of the thermodynamic model reveals a predicted region of liquid-liquid un-mixing at low-temperatures (<1624 K) and pressures (<1.257 GPa). These pressure-temperature conditions are above the glass transition temperature but within the metastable liquid region. They represent the highest temperatures yet suggested for liquid-liquid un-mixing in a silicate bulk composition. A shock wave Hugoniot curve is calculated for comparison with the experimental data of Rigden et al. (1989) and of Asimow and Ahrens (2008). The comparison suggests that the model developed in this paper underestimates the density of the liquid by roughly 10% at pressures greater than ∼20 GPa.  相似文献   

15.
A large number of aqueous metal complexes contribute significantly to hydrothermal, metamorphic, and magmatic processes in the crust of the Earth. Nevertheless, relatively few thermodynamic data other than dissociation constants (K) for a few dozen of these complexes have been determined experimentally at elevated temperatures and pressures. The calculations summarized below are intended to supplement these experimental data by providing interim predictions of the thermodynamic properties of supercritical aqueous metal complexes using the revised HKF (Helgeson et al., 1981) equations of state for aqueous species (Tanger and Helgeson, 1988; Shock et al., 1992) and correlations among equations of state parameters and standard partial molal properties at 25 degrees C and 1 bar (Shock and Helgeson, 1988, 1990; Shock et al., 1989). These equations and correlations permit retrieval of the conventional standard partial molal entropies (mean S0), volumes (mean V0), and heat capacities (mean C0P) of aqueous metal complexes at 25 degrees C and 1 bar from published values of log K in the supercritical region and the limited number of experimental dissociation constants available in the literature over relatively short ranges of elevated temperature at PSAT (PSAT and SAT are used in the present communication to refer to pressures corresponding to liquid-vapor equilibrium for the system H2O except at temperatures <100 degrees C, where they refer to the reference pressure of 1 bar). The standard partial molal properties computed in this way can then be used to generate corresponding values of delta mean S0, delta mean V0, and delta mean C0P of association, which for similar complexes correlate linearly with mean S0, mean V0 and mean C0P, respectively, of the constituent cations and ligands at 25 degrees C and 1 bar. Generalizing these correlations and combining them with the equations of state permits prediction of the temperature and pressure dependence of log K and other thermodynamic properties of a large number of aqueous metal complexes. As a consequence, it is possible to retrieve values of log K at 25 degrees C and 1 bar from the results of hydrothermal experiments at higher temperatures and pressures or to predict values of log K at hydrothermal conditions when no experimental data are available at temperatures and pressures above 25 degrees C and 1 bar. Such predictions can be made for temperatures and pressures from 0 degrees C and 1 bar to 1000 degrees C and 5000 bars.  相似文献   

16.
A thermodynamic model is presented to calculate methane solubility, liquid phase density and gas phase composition of the H2O-CH4 and H2O-CH4-NaCl systems from 273 to 523 K (possibly up to 573 K), from 1 to 2000 bar and from 0 to 6 mol kg−1 of NaCl with experimental accuracy. By a more strict theoretical approach and using updated experimental data, this model made substantial improvements over previous models: (1) the accuracy of methane solubility in pure water in the temperature range between 273 and 283 K is increased from about 10% to about 5%, but confirms the accuracy of the Duan model [Duan Z., Moller N., Weare J.H., 1992a. Prediction of methane solubilities in natural waters to high ionic strength from 0 to 250 °C and from 0 to 1600 bar. Geochim. Cosmochim. Acta56, 1451-1460] above 283 K up to 2000 bar; (2) the accuracy of methane solubility in the NaCl aqueous solutions is increased from >12% to about 6% on average from 273 K and 1 bar to 523 K and 2000 bar; (3) this model is able to calculate water content in the gas phase and liquid phase density, which cannot be calculated by previous models; and (4) it covers a wider range of temperature and pressure space. With a simple approach, this model is extended to predict CH4 solubility in other aqueous salt solutions containing Na+, K+, Mg2+, Ca2+, Cl− and , such as seawater and geothermal brines, with excellent accuracy. This model is also able to calculate homogenization pressure of fluid inclusions (CH4-H2O-NaCl) and CH4 solubility in water at gas-liquid-hydrate phase equilibrium. A computer code is developed for this model and can be downloaded from the website: www.geochem-model.org/programs.htm.  相似文献   

17.
A unified equation has been derived by using all available data for calculating methane vapor pressures with measured Raman shifts of C-H symmetric stretching band (υ1) in the vapor phase of sample fluids near room temperature. This equation eliminates discrepancies among the existing data sets and can be applied at any Raman laboratory. Raman shifts of C-H symmetric stretching band of methane in the vapor phase of CH4-H2O mixtures prepared in a high-pressure optical cell were also measured at temperatures between room temperature and 200 °C, and pressures up to 37 MPa. The results show that the CH4υ1 band position shifts to higher wavenumber as temperature increases. We also demonstrated that this Raman band shift is a simple function of methane vapor density, and, therefore, when combined with equation of state of methane, methane vapor pressures in the sample fluids at elevated temperatures can be calculated from measured Raman peak positions. This method can be applied to determine the pressure of CH4-bearing systems, such as methane-rich fluid inclusions from sedimentary basins or experimental fluids in hydrothermal diamond-anvil cell or other types of optical cell.  相似文献   

18.
Viscosity of silicate melts is a critical property for understanding volcanic and igneous processes in the Earth. We investigate the pressure effect on the viscosity of rhyolitic melts using two methods: indirect viscosity inference from hydrous species reaction in melts using a piston cylinder at pressures up to 2.8 GPa and direct viscosity measurement by parallel-plate creep viscometer in an internally-heated pressure vessel at pressures up to 0.4 GPa. Comparison of viscosities of a rhyolitic melt with 0.8 wt% water at 0.4 GPa shows that both methods give consistent results. In the indirect method, viscosities of hydrous rhyolitic melts were inferred based on the kinetics of hydrous species reaction in the melt upon cooling (i.e., the equivalence of rheologically defined glass transition temperature and chemically defined apparent equilibrium temperature). The cooling experiments were carried out in a piston-cylinder apparatus using hydrous rhyolitic samples with 0.8-4 wt% water. Cooling rates of the kinetic experiments varied from 0.1 K/s to 100 K/s; hence the range of viscosity inferred from this method covers 3 orders of magnitude. The data from this method show that viscosity increases with increasing pressure from 1 GPa to 3 GPa for hydrous rhyolitic melts with water content ?0.8 wt% in the high viscosity range. We also measured viscosity of rhyolitic melt with 0.13 wt% water using the parallel-plate viscometer at pressures 0.2 and 0.4 GPa in an internally-heated pressure vessel. The data show that viscosity of rhyolitic melt with 0.13 wt% water decreases with increasing pressure. Combining our new data with literature data, we develop a viscosity model of rhyolitic melts as a function of temperature, pressure and water content.  相似文献   

19.
We have performed first-principles molecular dynamics simulations of CaAl2Si2O8 (anorthite) liquid at pressures up to 120 GPa and temperatures of 3000, 4000 and 6000 K. At the lowest degrees of compression the liquid is seen to accommodate changes in density through decreasing the abundance of 3- and 4-membered rings, while increases in coordination of network forming cations take effect at somewhat higher degrees of compression. Results are fit to a fundamental thermodynamic relation with 4th order finite strain and 1st order thermal variable expansions. Upon compression by a factor of two, the Grüneisen parameter (γ) is found to increase continuously from 0.35 to 1.10. Weak temperature dependence in γ is thermodynamically consistent with a slight decrease in isochoric heat capacity (CV), for which values of between 4.4 and 5.2 NkB are obtained, depending on the temperature. Pressure and temperature dependence of self-diffusivities is found to be well represented by an Arrhenius relation, except at 3000 K and pressures lower than 5 GPa, where self-diffusivities of Si, Al, and O increase with pressure. Analysis of the lifetimes of individual coordination species reveals that this phenomenon arises due to the disproportionately high stability of 4-fold coordinated Si, and to a lesser extent 4-fold coordinated Al. Our results represent a marked improvement in accuracy and reliability in describing the physics of CaAl2Si2O8 liquid at deep mantle pressures, pointing the way to a general thermodynamic model of melts at extreme pressures and temperatures relevant to planetary-scale magma oceans and deep mantle partial melting.  相似文献   

20.
The structure of silicate melts in the system Na2O·4SiO2 saturated with reduced C-O-H volatile components and of coexisting silicate-saturated C-O-H solutions has been determined in a hydrothermal diamond anvil cell (HDAC) by using confocal microRaman and FTIR spectroscopy as structural probes. The experiments were conducted in-situ with the melt and fluid at high temperature (up to 800 °C) and pressure (up to 1435 MPa). Redox conditions in the HDAC were controlled with the reaction, Mo + H2O = MoO2 + H2, which is slightly more reducing than the Fe + H2O = FeO + H2 buffer at 800 °C and less.The dominant species in the fluid are CH4 + H2O together with minor amounts of molecular H2 and an undersaturated hydrocarbon species. In coexisting melt, CH3 - groups linked to the silicate melt structure via Si-O-CH3 bonding may dominate and possibly coexists with molecular CH4. The abundance ratio of CH3 - groups in melts relative to CH4 in fluids increases from 0.01 to 0.07 between 500 and 800 °C. Carbon-bearing species in melts were not detected at temperatures and pressures below 400 °C and 730 MPa, respectively. A schematic solution mechanism is, Si-O-Si + CH4?Si-O-CH3+H-O-Si. This mechanism causes depolymerization of silicate melts. Solution of reduced (C-O-H) components will, therefore, affect melt properties in a manner resembling dissolved H2O.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号