锡在流体和花岗质熔体间分配行为的实验研究

锡矿床是与花岗岩在时间、空间、成因上有着密切联系的典型矿种之一。与锡矿有关的花岗岩多具有过铝、富钾、硅含量高的特征。传统观点认为与锡矿有关的花岗岩主要是S型花岗岩，可是近年在国内外相继发现了许多具有重要经济价值的锡矿床与富碱侵入岩有着密切成因联系。与富碱侵入岩有关的锡成矿作用日益受到地质学家的重视，锡矿床与富碱侵入岩的关系已成为研究热点之一，相关的研究工作虽然取得了很大的进展，但是富碱侵入岩体能否分异出富锡成矿流体还存在争议。研究表明，与岩浆岩有关的成矿与岩浆演化过程中成矿元素在流体-熔体相间的分配行为有着密切的关系。成矿元素在流体-熔体相间的分配行为除受到温度、压力及氧逸度等物理化学条件的制约外，还受到岩浆熔体成份及岩浆分异出来的流体化学组成的影响。以往有关锡在流体-熔体相间分配行为的实验研究主要侧重于改变流体相来观测锡的分配系数，且多为单一的含氯或含氟岩浆体系，这制约了对岩浆演化过程中元素在流体-熔体相间分配行为的深入认识。本文通过改变流体相、熔体相的化学组成，开展了一系列锡在流体和花岗质熔体相间分配行为的实验研究。综合分析了锡在晶体-熔体-流体间的分配行为，并结合地质实际探讨与富碱侵入岩有关的锡成矿的物理化学条件和成矿机理。研究成果对深入认识与花岗岩有关锡矿的成矿机理、丰富和完善与花岗岩有关的锡成矿理论、为进一步探索与花岗岩有关的锡成矿规律提供重要的实验依据。此外，实验对进一步推动实验地球化学学科发展具有重要意义。 实验在中国科学院地球化学研究所矿床地球化学国家重点实验室的成矿实验室完成，主要实验设备为快速内冷（RQV）高压釜。实验的温度为850℃，压力100MPa、氧逸度接近NNO。实验首先采用人工合成硅酸盐凝胶的方法制成具有不同化学组成的花岗质熔体，使用分析纯化学试剂配制不同成分和不同浓度的溶液，分别作为实验初始固液相。主要开展了三方面的实验研究：1．熔体相组成不变，以改变流体相组成来观察锡的分配行为。这组实验固相初始物为过碱质富钾的硅酸盐，初始液相分别为NaCl、KCl、HCl、HF、去离子水溶液；2．流体相组成不变，改变熔体化学组成观察锡分配行为。初始液相选用低浓度的0.1mol/L HCl溶液，熔体相为具不同化学组成的凝胶（其中一组改变熔体碱质含量和铝饱和指数ASI、另一组改变熔体钠钾摩尔比值）；3．氟氯共存含水的花岗质岩浆体系中氟氯含量相对变化时锡分配行为。实验通过改变熔体相中氟含量和液相盐酸溶液的浓度来观察锡在含氟硅酸盐熔体和不同浓度盐酸溶液间的分配行为。氟主要以(NaF+KF)混合物的形式加入初始固相中。实验研究结果表明： 1．流体相络阴离子种类及含量对锡在流体-熔体相间的分配行为有着明显的影响。当流体相中络阴离子Cl－、F－含量增大时，有利于增大锡在流体-熔体相间的分配系数；尤其当流体为富氯的酸性流体时，锡在流体-熔体相间的分配系数随液相中HCl浓度的增大而增大并存在关系式logD Sn=2.0247×log[HCl]＋0.6717（[HCl]的单位为mol/L），锡在流体相中主要以二价锡氯配合物的形式迁移，锡倾向于分配进入富氯的酸性流体中。此外，富氯酸性流体与共存的熔体反应后，熔体中的碱质含量降低，铝饱和指数增大。 2．熔体化学组成对锡在熔体相/流体相的分配行为有着明显的影响。D Sn随着熔体中碱质含量增大而减小：D Sn=－0.0489×MAlk＋0.4516, R2=0.98（MAlk为熔体中Na2O+K2O摩尔含量），表明富碱质熔体有利于锡在熔体相中富集，从而可能为锡矿形成提供矿质来源。D Sn随熔体ASI值的增大而增大：D Sn=0.1886×ASI－0.1256, R2=0.99，即过铝质熔体相对有利于锡分配进入流体相中。过铝质熔体中碱质总量及其它组分相对不变的前提下，熔体钠钾摩尔比值越高D Sn越小：D Sn=－0.0314×RNa/K＋0.0483, R2=0.82（RNa/K为Na/K摩尔比值），富钠的熔体有利于锡分配进入熔体相，而富钾的熔体却相对有利于锡分配进入流体相中。 3．在氟氯共存花岗质岩浆体系中：①熔体相中氟含量对氯在流体-熔体相间的分配有着明显影响，熔体中氟含量降低有利于氯分配进入流体相。②熔体中氟含量大于约1 wt%后，D Sn小于0.1且变化不大，当液相富含HCl且熔体中氟含量从约1 wt%降低后，D Sn 迅速增大，即熔体中氟含量小于约1 wt%后锡倾向于分配进入富氯的酸性流体中。而富氟（F含量大于约1 wt%）的熔体有利于萃取锡并使锡在熔体相中富集。③熔体铝饱和指数ASI值越大，相应锡的分配系数越大；流体相中HCl浓度增大时，锡分配系数随之增大；当熔体为过铝质的花岗质熔体、流体富含HCl时有利于锡分配进入流体相。 分析总结与花岗岩有关的锡成矿特征和锡在不同晶体相和熔体相间的分配行为得出：壳源铝质、富碱、富挥发份、贫钙铁镁的岩浆在结晶分异演化过程中相对有利于锡在残余熔流体相中富集。因此，具有这些特征的岩浆结晶分异演化产生的晚期岩浆可富含锡，能为后期锡矿床的形成提供矿质来源。这种富锡富挥发份的岩浆在上侵过程中，当温度压力降低、岩浆水饱和度增大、硅含量增大、熔体相氟含量降低时，可分异出含氯富锡的成矿流体。 根据上述结论，分析了与湖南芙蓉锡矿床有着密切成因联系的骑田岭花岗岩的岩石化学特征、成岩成矿物理化学条件，得出芙蓉锡矿床成矿流体可由骑田岭晚期岩浆分异产生。 最后得出如下认识：1)当花岗质岩浆体系水不饱和、流体相络阴离子浓度低的情况下，锡倾向于分配进入熔体相中；2)水饱和富含挥发份的过铝、富钾的岩浆体系有利于锡分配进入流体相；3) 铝质、富钾、富挥发份的富碱侵入岩岩浆演化过程中可在有利的物理化学条件下分异出富锡的流体相，与芙蓉超大型锡矿床有成因联系的骑田岭富碱侵入岩体成岩过程中可分异出富锡的成矿流体。

Investigation of Interfacial Tensions for Carbon Dioxide Aqueous Solutions by Perturbed-Chain Statistical Associating Fluid Theory Combined with Density-Gradient Theory

The perturbed-chain statistical associating fluid theory and density-gradient theory are used to construct an equation of state (EOS) applicable for the phase behaviors of carbon dioxide aqueous solutions. With the molecular parameters and influence parameters respectively regressed from bulk properties and surface tensions of pure fluids as input, both the bulk and interfacial properties of carbon dioxide aqueous solutions are satisfactorily correlated by adjusting the binary interaction parameter (k(ij)). Our results show that the constructed EOS is able to describe the interfacial properties of carbon dioxide aqueous solutions in a wide temperature range, and illustrate the influences of temperature, pressure, and densities in each phase on the interfacial properties.

Partitioning behavior of amino acids in aqueous two-phase systems containing polyethylene glycol and phosphate buffer

The partitioning behavior of four amino acids, cysteine, phenylalanine, methionine, and lysine in 15 aqueous two-phase systems (ATPSs) with different polyethylene glycol (PEG) molecular weights and phosphate buffers has been studied in the present paper. The phase diagrams of the systems are investigated together with the effect of the PEG molecular weight and pH of the phosphate solutions. The composition of these systems and some parameters such as density and refractive index are determined. The influences of salts in ATPSs, side chain structure of the amino acids, pH of ATPSs, and the PEG molecular weight on the distribution ratios of the amino acids have been studied. This work is useful for the purification of amino acids and the separation of some proteins whose main surface exposed amino acid residues are these four amino acids, respectively.

Spectroscopic investigation of the function of aqueous 2-hydroxyethylmethacrylate/glutaraldehyde solution as a dentin desensitizer

Fourier-transform (FT)-Raman and -infrared (IR) spectroscopy were employed to investigate the function of the aqueous 2-hydroxyethylmethacrylate/glutaraldehyde solution (Gluma) as a desensitizer. 2-Hydroxyethylmethacrylate (HEMA), glutaraldehyde (GA), and the mixture of HEMA/GA (i.e. Gluma) were used to interact with dentin, collagen, hydroxyapatite (HAP), and bovine serum albumin (BSA) individually. All the interactions were monitored by an FT-Raman spectrometer. FT-IR spectroscopy was also used in this study. The results show that HEMA could be absorbed by dentin and collagen; GA could cross-link collagen and BSA; and when BSA was added to Gluma, polymerization of HEMA occurred. The results suggest that Gluma acts as a desensitizer whereby, first, GA reacts with part of the serum albumin in dentinal fluid, which induces a precipitation of serum albumin, then, second, a reaction of GA with serum albumin induces polymerization of HEMA. The function of Gluma as a desensitizer to block dentinal tubules occurs via these two reactions.

Micro-shock waves generated inside a fluid jet impinging on plane wall

An Nd:glass laser pulse (18 ns, 1.38 J) is focused in a tiny area of about 100-mum diam under ambient conditions to produce micro-shock waves. The laser is focused above a planar surface with a typical standoff distance of about 4 mm, The laser energy is focused inside a supersonic circular jet of carbon dioxide gas produced by a nozzle with internal diameter of 2.9 mm and external diameter of 8 mm, Nominal value of the Mach number of the jet is around 2 with the corresponding pressure ratio of 7.5 (stagnation pressure/static pressure at the exit of the nozzle), The interaction process of the micro-shock wave generated inside the supersonic jet with the plane wall is investigated using double-pulse holographic interferometry. A strong surface vortex field with subsequent generation of a side jet propagating outward along the plane wail is observed. The interaction of the micro-shock wave with the cellular structure of the supersonic jet does not seem to influence the near surface features of the flowfield. The development of the coherent structures near the nozzle exit due to the upstream propagation of pressure waves seems to be affected by the outward propagating micro-shock wave. Mach reflection is observed when the micro-shock wave interacts with the plane wall at a standoff distance of 4 mm, The Mach stem is slightly deflected, indicating strong boundary-layer and viscous effects near the wall. The interaction process is also simulated numerically using an axisymmetric transient laminar Navier-Stokes solver. Qualitative agreement between experimental and numerical results is good.

Parameters Inversion of Fluid-Saturated Porous Media Based on Neural Networks

The multi-layers feedforward neural network is used for inversion of material constants of fluid-saturated porous media. The direct analysis of fluid-saturated porous media is carried out with the boundary element method. The dynamic displacement responses obtained from direct analysis for prescribed material parameters constitute the sample sets training neural network. By virtue of the effective L-M training algorithm and the Tikhonov regularization method as well as the GCV method for an appropriate selection of regularization parameter, the inverse mapping from dynamic displacement responses to material constants is performed. Numerical examples demonstrate the validity of the neural network method.

Three-Dimensional Modeling of Heat Transfer and Fluid Flow in Laminar-Plasma Material Re-Melting Processing

Modeling study is performed concerning the heat transfer and fluid flow for a laminar argon plasma jet impinging normally upon a flat workpiece exposed to the ambient air. The diffusion of the air into the plasma jet is handled by using the combined-diffusion-coefficient approach. The heat flux density and jet shear stress distributions at the workpiece surface obtained from the plasma jet modeling are then used to study the re-melting process of a carbon steel workpiece. Besides the heat conduction within the workpiece, the effects of the plasma-jet inlet parameters (temperature and velocity), workpiece moving speed, Marangoni convection, natural convection etc. on the re-melting process are considered. The modeling results demonstrate that the shapes and sizes of the molten pool in the workpiece are influenced appreciably by the plasma-jet inlet parameters, workpiece moving speed and Marangoni convection. The jet shear stress manifests its effect at higher plasma-jet inlet velocities, while the natural convection effect can be ignored. The modeling results of the molten pool sizes agree reasonably with available experimental data.

Damping of Vertically Excited Surface Wave in Weakly Viscous Fluid

In a vertically oscillating circular cylindrical container, singular perturbation theory of two-time scale expansions is developed in weakly viscous fluids to investigate the motion of single free surface standing wave by linearizing the Navier-Stokes equation. The fluid field is divided into an outer potential flow region and an inner boundary layer region. The solutions of both two regions are obtained and a linear amplitude equation incorporating damping term and external excitation is derived. The condition to appear stable surface wave is obtained and the critical curve is determined. In addition, an analytical expression of damping coefficient is determined. Finally, the dispersion relation, which has been derived from the inviscid fluid approximation, is modified by adding linear damping. It is found that the modified results are reasonably closer to experimental results than former theory. Result shows that when forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when forcing frequency is high, the surface tension of the fluid is prominent.

Numerical Analysis on Intermittent Flow in a Longitudinal Section of Bingham Fluid along a Slope

The transition process of intermittent flow in a longitudinal section of Bingham fluid from initial distribution to fully developed state was numerically investigated in this paper. The influences of slope dimensionless runoff Q* and viscosity μ0* on the dimensionless surge speed U* were also presented in a wide range of parameters. By one typical example, the intermittent flow possessed wave characteristics and showed a supercritical flow conformation for a fully developed flow. The distributions of gravity and bed drag along the flow path and the velocity distribution of flow field were also analyzed.

Fluid-Solid Interaction in Particle-Laden Flows

Liu Qingquan, Singh V.P

Effects of Baffle Design on Fluid Flow and Heat Transfer in Ammonothermal Growth System of Nitrides

Efforts have been made in growing bulk single crystals of GaN front supercritical fluids using the ammonothermal method, which utilizes ammonia as fluid rather than water as in the hydrothermal process. Different mineralizers such as amide or azide and temperatures in the range of 200-600degreesC have been used to increase the solubility. The pressure is from 1 to 4 kbar. Modeling of the ammonothermal growth process has been used to identify factors which may affect the temperature distribution, fluid flow and nutrient transport. The GaN charge is considered as a porous media bed and the flow in the charge is simulated using the Darcy-Brinkman-Forchheimer model. The resulting governing equations are solved using the finite volume method. The effects of baffle design and opening on flow pattern and temperature distribution in an autoclave are analyzed. Two cases are considered with baffle openings of 15% and 20% in cross-sectional area, respectively.

Two Bifurcation Transition Processes in Floating Half Zone Convection of Larger Prandtl Number Fluid

Processes of the onset oscillation in the thermocapillaxy convection under the Earth's gravity are investigated by the numerical simulation and experiments in a floating half zone of large Prandtl number with different volume ratio. Both computational and experimental results show that the steady and axisymmetric convection turns to the oscillatory convection of m=1 for the slender liquid bridge, and to the oscillatory convection before a steady and 3D asymmetric state for the case of a fat liquid bridge. It implies that, there are two critical Marangoni numbers related, respectively, to these two bifurcation transitions for the fat liquid bridge. The computational results agree with the results of ground-based experiments.