A Microtube Viscometer with a Thermostat

The viscometer presented in this paper is suit-able for measuring the viscosity of liquids in micro-litre quantities. It consists of a micro-flow experimental system with a thermostat. Using the measurements of the flow rates and pressure drops of a liquid passing through a microtube, the liquid's viscosity can be calculated from on Hagen-Poiseuille theory. After calibration, the viscometer was used to measure viscosities of deionized water and ethyl alcohol at temperatures ranging from 0 to 40 "C. For both test liquids, the relative deviation of the measured values from those quoted in the literature (obtained using other viscometers) was less than 2.6o/o. The relative uncertainty of the experimental system was reduced to +-l.8% using the relative measuring method. Due to the micro-scale of the test section, only a micro-litre quantity of liquid is needed for a test, this is a potential advantage for measurement of bio-liquid viscosities.

Quantum Confinement Effects and Electronic Properties of SnO2 Quantum Wires and Dots

On the basis of the density functional theory (DFT) within local density approximations (LDA) approach, we calculate the band gaps for different size SnO2 quantum wire (QWs) and quantum dots (QDs). A model is proposed to passivate the surface atoms of SnO2 QWs and QDs. We find that the band gap increases between QWs and bulk evolve as Delta E-g(wire) = 1.74/d(1.20) as the effective diameter d decreases, while being Delta E-g(dot) = 2.84/d(1.26) for the QDs. Though the similar to d(1.2) scale is significantly different from similar to d(2) of the effective mass result, the ratio of band gap increases between SnO2 QWs and QDs is 0.609, very close to the effective mass prediction. We also confirm, although the LDS calculations underestimate the band gap, that they give the trend of band gap shift as much as that obtained by the hybrid functional (PBE0) with a rational mixing of 25% Fock exchange and 75% of the conventional Perdew-Burke-Ernzerhof (PBE) exchange functional for the SnO2 QWs and QDs. The relative deviation of the LDA calculated band gap difference Lambda E-g compared with the corresponding PBE0 results is only within 5%. Additionally, it is found the states of valence band maximum (VBM) and conduction band minimum (CBM) of SnO2 QWs or QDs have a mostly p- and s-like envelope function symmetry, respectively, from both LDA and PBE0 calculations.

Analysis of sulfur-containing compounds in crude oils by comprehensive two-dimensional gas chromatography with sulfur chemiluminescence detection

This paper reports an analytical method for separating, identifying, and quantifying sulfur-containing compounds in crude oil fraction (IBP-360degreesC) samples based on comprehensive two-dimensional gas chromatography coupled with a sulfur chemiluminescence detector. Various sulfur-containing compounds and their groups were analyzed with one direct injection. 3620 peaks were detected including 1722 thiols/thioethers/ disulfides/1-ring thiophenes, 953 benzothiophenes, 704 dibenzothiophenes, and 241 benzonaphthothiophenes. The target sulfur compounds and their groups were identified based on the group separation feature and structured retention of comprehensive two-dimensional gas chromatography as well as standard substances. The quantitative analysis of major sulfur-containing compounds and total sulfur was based on the linear response of the sulfur chemiluminescence detector using the internal standard method. The sulfur contents of target sulfur compounds and their groups in 4 crude oil fractions were also determined. The recoveries for standard sulfur-containing compounds were in the range of 90-102%. The quantitative result of total sulfur in the Oman crude oil fraction sample was compared with those from ASTM D 4294 standard method (total S by X-ray fluorescence spectrometry), the relative deviation (RD%) was 4.2% and the precision of the method satisfactory.

潮波模型的优化开边界方法研究

数值模式是潮波研究的一种有利手段，但在研究中会面临各种具体问题，包括开边界条件的确定、底摩擦系数和耗散系数的选取等。数据同化是解决这些问题的一种途径，即利用有限数量的潮汐观测资料对潮波进行最优估计，其根本目的是迫使模型预报值逼近观测值，使模式不要偏离实际情况太远。本文采用了一种优化开边界方法，沿着数值模型的开边界优化潮汐水位信息，目的是设法使数值解在动力约束的意义下接近观测值，获得研究区域的潮汐结果。边界值由指定优化问题的解来定，以提高模拟区域的潮汐精度，最优问题的解是基于通过开边界的能量通量的变化，处理开边界处的观测值与计算值之差的最小化。这里提供了辐射型边界条件，由Reid 和Bodine（本文简称为RB）推导，我们将采用的优化后的RB方法（称为ORB）是优化开边界的特殊情况。 本文对理想矩形海域（ E- E， N- N， 分辨率 ）进行了潮波模拟，有东部开边界，模式采用ECOM3D模式。对数据结果的误差分析采用，振幅平均偏差，平均绝对偏差，平均相对误差和均方根偏差四个值来衡量模拟结果的好坏程度。 需要优化入开边界的解析潮汐值本文采用的解析解由方国洪《海湾的潮汐与潮流》（1966年）方法提供，为验证本文所做的解析解和方文的一致，本文做了其第一个例子的关键值a,b,z,结果与其结果吻合的相当好。但略有差别，分析的可能原因是两法在具体迭代方案和计算机保留小数上有区别造成微小误差。另外，我们取m=20，得到更精确的数值，我们发现对前十项的各项参数值，取m=10,m=20各项参数略有改进。当然我们可以获得m更大的各项参数值。 同时为了检验解析解的正确性讨论m和l变化对边界值的影响，结果指出，增大m，m=20时，u的模最大在本身u1或u2的模的6%；m=100时，u的模最大在本身u1或u2的模的4%；m再增大，m=1000时，u的模最大在本身u1或u2的模的4%，改变不大。当l<1时， =0处u的模最大为2。当l=1时， =0处u的模最大为0.1，当l>1时，l越大，u的模越小，当l=10时，u的模最大为0.001，可以认为为0。 为检验该优化方法的应用情况，我们对理想矩形区域进行模拟，首先将本文所采用的优化开边界方法应用于30m的情况，在开边界优化入开边界得出模式解，所得模拟结果与解析解吻合得相当好，该模式解和解析解在整个区域上，振幅平均绝对偏差为9.9cm,相位平均绝对偏差只有4.0 ，均方根偏差只有13.3cm，说明该优化方法在潮波模型中有效。 为验证该优化方法在各种条件下的模拟结果情况，在下面我们做了三类敏感性试验： 第一类试验：为证明在开边界上使用优化方法相比于没有采用优化方法的模拟解更接近于解析解，我们来比较ORB条件与RB条件的优劣，我们模拟用了两个不同的摩擦系数，k分别为：0,0.00006。 结果显示，针对不同摩擦系数，显示在开边界上使用ORB条件的解比使用RB条件的解无论是振幅还是相位都有显著改善，两个试验均方根偏差优化程度分别为84.3%，83.7%。说明在开边界上使用优化方法相比于没有采用优化方法的模拟解更接近于解析解，大大提高了模拟水平。上述的两个试验得出， k=0.00006优化结果比k=0的好。 第二类试验，使用ORB条件确定优化开边界情况下，在东西边界加入出入流的情况，流考虑线性和非线性情况，结果显示，加入流的情况，潮汐模拟的效果降低不少，流为1Sv的情况要比5Sv的情况均方根偏差相差20cm,而不加流的情况只有0.2cm。线性流和非线性流情况两者模式解相差不大，振幅，相位各项指数都相近, 说明流的线性与否对结果影响不大。 第三类试验，不仅在开边界使用ORB条件，在模式内部也使用ORB条件，比较了内部优化和不优化情况与解析解的偏差。结果显示，选用不同的k,振幅都能得到很好的模拟，而相位相对较差。另外，在内部优化的情况下，考虑不同的k的模式解， 我们选用了与解析解相近的6个模式解的k,结果显示，不同的k,振幅都能得到很好的模拟，而相位较差。 总之，在开边界使用ORB条件比使用RB条件好，振幅相位都有大幅度改进，在加入出入流情况下，流的大小对模拟结果有影响，但线形流和非线性流差别不大。内部优化的结果显示，模式采用不同的k都能很好模拟解析解的振幅。

利用MC-ICP-MS测定硼同位素及其在西藏地热和盐湖中的初步应用

Based on previous studies, boron can be separated from aqueous samples with Amberlite IRA-743 resin. Experiments on the elute temperature, elute volume and the dynamic resin exchange capacity have been performed in this study. Results show that the dynamic exchange capacity of the resin is 4.2mg B/g and at room temperature, boron fixed on the resin within this capacity level can be extracted quantitatively by using 5ml 2%HNO3. A new procedure has been developed for the measurement of boron isotope ratios in water samples using a Neptune MC-ICP-MS, after resolving the memory effect, which is a key problem, and investigating the impacts of mass bias and Si matrix effect. Using this method, it usually takes 20 min to perform one measurement on 0.1ppm boron solution with a precision of 0.23‰ (SD). If the relative deviation between a sample and the standard is large, the washout time needs to be doubled to achieve a higher precision. δ11B values of water samples from Yangbajing geothermal field vary from -10.53 to -9.13‰. Owing to the large difference B concentration and the small B isotope difference between deep geothermal water and surface water, B isotope ratios of the shallow geothermal fluids are dominated by the deep end member rather than the shallower one in the mixing process. As a consequence, δ11B-B relation is indicative basically of a dilution process. Vapor-liquid separation and calcite scaling also greatly influence B isotope fractionation. δ11B values of water samples from Dagejia geothermal field are from -15.98‰ to -11.67‰. Boron in Changma River near the field has two sources, freshwater lakes (Dajiamang Lake and Canke Lake) and geothermal waters. Finally, a preliminary discussion is included on boron geochemical characteristics of the salt lakes in Shuanghu area and other geothermal fields, to provide information for future studies on boron isotope geochemistry of geothermal systems and salt lakes in Tibet.

磁层-电离层耦合现象的分析研究

The magnetosphere-ionosphere coupling is mainly manifested by the trans- porting processes of energy into the ionosphere , the energy is carried by solar wind and firstly accumulate at the magnetosphere, and the coupling processes also significantly include the interaction between the magnetosphere and ionosphere for mass and energy. At the quiet condition, energy is delivered by the large-scale convection of the geomagnetic field; the huge energy from solar wind bulk will be injected into and consumed at the near magnetosphere and ionosphere by the geomagnetic storm and substorm activities. Aurorae and FACs (Field-aligned currents) are the important phenomena in the coupling processes. In the present work, firstly, we analyze the activity characteristics of auroral precipitating particle, secondly, we study the distribution characters of large-scale field aligned currents (LS FACs) at storm-time using the observations from different satellites at different altitudes. Finally, we investigate the evolution of the geomagnetic field configuration at the nightside sector on the onset of the expansion phase in a substorm event, the substorm event happened at 0430UT to 0630UT on 8th Nov. 2004. The main results as follows： At the first, the data of the estimated power input (EPI) of auroral particles from NOAA/POES (Polar orbiting environmental satellite) for some 30 years have been analyzed. The variation tendencies of the EPI generally coincide with aa, AE and Dst indices. The annual variation of EPI shows equinox peaks and an asymmetric-activity with a higher peak in the winter-hemisphere than in the summer-hemisphere. The diurnal UT variations are different from north and south hemisphere: for north hemisphere, the peak appears at 1200UT, and the relative deviation is 22% to the daily average of the north hemisphere. For south hemisphere, the maximal deviation is 22% at 2000UT. So the diurnal variation of EPI is more dominant than the annual variation which maximal deviation is 3% to 12% for different seasons. Studies on correlations of the hourly average of EPI, Pa, with AE and Dst indices show a correlation coefficient r=0.74 of Pa and AE, and r=-0.55 of Pa and Dst. The hourly EPIs for north and south polar regions, NPa and SPa, show a north-south asymmetry with a higher correlation of SPa and AE (or Dst). Time delays of EPI with respect to magnetic indices are examined, the maximum correlation coefficient of Pa with AE (r=0.78) occurs when the time delay =0, suggesting a synchronous activity of auroral electrojet and auroral precipitating particles, while =1-2h, the correlation coefficient of Pa with Dst is maximum (r=0.57), suggesting that the activity of auroral particle precipitating may influence the ring current on some extent. Sencondly, we use the high-resolution magnetic field vector data of the CHAMP satellite to investigate the distribution of large-scale FACs during the great magnetic storm on 7th to 8th Nov. 2004. The results show that, whether in the northern or southern hemisphere, the number and density of large-scale FACs during the main-phase are more and bigger than these during the recover-phase, and the number of large-scale FACs in morning sector obviously is more than that in afternoon sector. In terms of the magnetic indices, we find that large-scale FACs in morning sector significantly affected by the substorm activities, while in afternoon sector the large-scale FACs mainly indicate the fluctuations of the ring-current in storm time. Accordingly to the former studies, similarly, we find that in the morning sector, the scale of the large-scale FACs move to the high-latitude region, and in the afternoon sector, large-scale FACs distinctly expand to the low-latitude region. During the time periods that the NOAA/POES auroral precipitating particle power data temporally correspond to the large-scale FACs, the more the power of auroral particle is, the more and bigger the number and density of FACs are. At the same time, we use the magnetic field vector data of POLAR obtain a good form of region 1, region 2, and three pieces of cusp FACs during a single transit at 1930UT-2006UT on 07th. And the characteristics of simultaneous electric field and energy particles observations on Polar are coincide with the five FACs pieces. Finally, by means of the observation of Cluster 4 and Goes 10、 Goes 12, we analyze the evolution process of the change of the magnetic field configuration at night sector at the expansion phase of a substorm event which happened during 0430UT to 0630UT on 8th Nov. 2004, we find that the times of the beginning of the polarizations of magnetic field are observed from Goes 10 to Goes 12 then to Cluster 4. So, at the synchronous orbit ( 6.6 RE) to 10RE distance scale of the neutral sheet, the current disruption spread tailward. Simultaneously, the strengthen of the FACs deduced from these satellites’ magnetic field observations are almost consistent with the times of polarizations, as well as the high energy particles injection and the electric field dominant variation. The onset times determined by the magnetic field polarizations from these satellites are all ahead of the onset time that confirmed from the auroral electrojet indices. So, these characters of different observations can be used as the criterions to determine the onset time for the substorms of such type as we studied.

Optimized Group Velocity Control Scheme and DNS of Decaying Compressible Turbulence of Relative High Turbulent Mach Number

To overcome the difficulty in the DNS of compressible turbulence at high turbulent Mach number, a new difference scheme called GVC8 is developed. We have succeeded in the direct numerical simulation of decaying compressible turbulence up to turbulent Mach number 0.95. The statistical quantities thus obtained at lower turbulent Mach number agree well with those from previous authors with the same initial conditions, but they are limited to simulate at lower turbulent Mach numbers due to the so-called start-up problem. The energy spectrum and coherent structure of compressible turbulent flow are analysed. The scaling law of compressible turbulence is studied. The computed results indicate that the extended self-similarity holds in decaying compressible turbulence despite the occurrence of shocklets, and compressibility has little effects on relative scaling exponents when turbulent Mach number is not very high.

3D simulation of the micro indentation process and relative problems research

In this paper the finite element method was used to simulate micro-scale indentation process. The several standard indenters were simulated with 3D finite element model. The emphasis of this paper was the differences between 2D axisymmetric cone model and

Analysis of Spacecraft Relative Motion and Applicability of Hill's Solution

研究了空间飞行器编队中最具基础性的问题之一,即相对运动的解析表达及Hill方程的适用条件。通过建立相对运动的通解公式,针对不同性质的初值深入地分析了其相对运动轨迹的本质特征,并给出了Hill方程的适用条件。此外,文中还给出了一个新的编队设计简化公式。

The effect of the incident relative phase on the four-wave mixing field and the electromagnetically induced transparency

In a A-type system employing a two-photon pump field, a four-wave mixing field can be generated simultaneously and, hence, a closed-loop system forms. We study theoretically the effect of the relative phase between the two incident fields on the generated four-wave mixing field and the electromagnetically induced transparency. It is found that the phase of the generated four-wave mixing field is the sum of the incident relative phase and a fixed phase that is irrelative to the incident relative phase. Hence, the total phase of the closed-loop system is independent of the incident relative phase. As a result, the incident relative phase has no effect on the electromagnetically induced transparency, which is different from the case of a A-type loop system closed by a third incident field. (c) 2005 Pleiades Publishing, Inc.

Enhancement and broadening of extreme-ultraviolet supercontinuum in a relative phase controlled two-color laser field

We experimentally demonstrate the generation of an extreme-ultraviolet (XUV) supercontinuum in argon with a two-color laser field consisting of an intense 7 fs pulse at 800 nm and a relatively weak 37 fs pulse at 400 nm. By controlling the relative time delay between the two laser pulses, we observe enhanced high-order harmonic generation as well as spectral broadening of the supercontinuum. A method to produce isolated attosecond pulses with variable width and intensity is proposed. (C) 2008 Optical Society of America.

Effects of relative phase in an open ladder system without incoherent pumping

We studied effects of the relative phase between the probe and driving fields on the absorption and dispersion properties in an open three-level ladder system with spontaneously generated coherence but without incoherent pumping. It is shown that by the phase controlling, switching from absorption to lasing without inversion (LWI) and enhancing remarkablely LWI gain can be realized; large index of refraction with zero absorption and the electromagnetically induced transparency can be obtained. We also find that varying the atomic injection and exit rates has a considerable influence on the phase dependent-absorption property of the probe field, existent of the atomic injection and exit rates gives the necessary condition of the realization of LWI, getting LWI is impossible in the corresponding closed system without incoherent pumping. We studied effects of the relative phase between the probe and driving fields on the absorption and dispersion properties in an open three-level ladder system with spontaneously generated coherence but without incoherent pumping. It is shown that by the phase controlling, switching from absorption to lasing without inversion (LWI) and enhancing remarkablely LWI gain can be realized; large index of refraction with zero absorption and the electromagnetically induced transparency can be obtained. We also find that varying the atomic injection and exit rates has a considerable influence on the phase dependent-absorption property of the probe field, existent of the atomic injection and exit rates gives the necessary condition of the realization of LWI, getting LWI is impossible in the corresponding closed system without incoherent pumping.

Coma measurement of projection optics in lithographic tools based on relative image displacements at multiple illumination settings

In this paper, we propose a novel method for measuring the coma aberrations of lithographic projection optics based on relative image displacements at multiple illumination settings. The measurement accuracy of coma can be improved because the phase-shifting gratings are more sensitive to the aberrations than the binary gratings used in the TAMIS technique, and the impact of distortion on displacements of aerial image can be eliminated when the relative image displacements are measured. The PROLITH simulation results show that, the measurement accuracy of coma increases by more than 25% under conventional illumination, and the measurement accuracy of primary coma increases by more than 20% under annular illumination, compared with the TAMIS technique. (c) 2007 Optical Society of America.