Contamination of a well in a uniform background flow

In this paper we describe the transport of pollution in groundwater in the neighbourhood of a well in a uniform background flow. We compute the rate at which contaminated particles reach the well as a function of the place of the source of pollution. The motion of a particle in a dispersive flow is seen as a random walk process. The Fokker-Planck equation for the random motion of a particle is transformed using the complex potential for the advective flow field. The resulting equation is solved asymptotically after a stretching transformation. Finally, the analytical solution is compared with results from Monte Carlo simulations with the random walk model. The method can be extended to arbitrary flow fields. Then by a numerical coordinate transformation the analytical results can still be employed.     

关于数字相位式测距仪的频率改正

本文在分析自动脉冲数字测相基本原理的基础之上，推出了不仅与测距频率有关，而且还与填充脉冲频率及低频测相频率有关的新的频率改正公式，并给出新的频率改正公式与传统频率改正公式之间等同的条件，即有锁相环路（或频率跟踪电路）的正常工作来保证。     

Three-dimensional structure of the crust and upper mantle beneath the south China area from long period surface waves

In the present study the long period surface wave records of 238 wave-paths from 79 earthquakes within China and its adjacent regions received by 30 seismic network stations are measured by using the improved match-filtering frequency-time analysis technique and the grid dispersion inversion method to obtain the rayleigh pure-path dispersion values for 147 slant grids of 4° × 4° in this area, then a three-dimensional shear wave velocity model of the crust and upper mantle beneath south China area to a depth of 170 km is inversed. It is found that there are obvious differences among the main structural units, and there are also certain differences among the subordinate elements even in the individual unit. The crustal thickness of this area is ranging from 30 to 43 km, and is getting thicker gradually from the east to the west. The average shear velocity of crust is ranging form 3.48 to 3.68 km/s with the lowest in the northeast part and highest in the west part. No obvious crustal low velocity layer of large scale is detected. There exist upper mantle low velocity zones in the most of south China area with the starting depth ranging from 75 to 106 km. The lowest shear velocity within the low velocity zones is about 4.28–4.38 km/s. Despite of the existing of upper mantle low velocity zones beneath the most of south China area, the interfaces between the important layers are quite clear, the variation of the bedding surfaces is very gentle, and the lateral changes measured in a larger scale of the underground structure are rather small. It may indicate that the crustal and upper mantle structure of the main part of south China area belongs to the relatively stable structure of the continental blocks except for the fringe areas such as the fold-faulted region in the west part and the fault system along the southeastern coast which may belong to the tectonically active area. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 159–167, 1993. This subject is supported by the National Science Foundation of China.     

A positivity preserving method for simulation of steam injection for NAPL site remediation

Numerical modelling of steam injection methods for cleanup of non-aqueous phase liquid (NAPL) contamination of groundwater requires consideration of multiphase, multicomponent convective and dispersive transport. Standard techniques do not ensure that the solution of the discrete equations has positive mole fractions, for finite mesh sizes. Negative mole fractions may cause the simulation to abort due to failure of the Newton iteration. A method for alleviating this problem is described. This method ensures that the mole fractions are positive, and results in an error that is the same size as the usual finite element discretization error. Example computations are presented for cartesian and axisymmetric two-dimensional geometries.     

Estimation of Scalar Source/Sink Distributions in Plant Canopies Using Lagrangian Dispersion Analysis: Corrections for Atmospheric Stability and Comparison with a Multilayer Canopy Model

Source/sink distributions of heat, water vapour andCO₂ within a rice canopy were inferred using aninverse Lagrangian dispersion analysis and measuredmean profiles of temperature, specific humidity andCO₂ mixing ratio. Monin–Obukhov similarity theorywas used to account for the effects of atmosphericstability on _w(z), the standard deviation ofvertical velocity and _L(z), the Lagrangian timescale of the turbulence. Classical surface layer scaling was applied in the inertial sublayer (z > z_ruf)using the similarity parameter = (z - d)/L, where z is height above ground, d is the zero plane displacementheight for momentum, L is the Obukhov length,and z_ruf 2.3h_c, where h_c iscanopy height. A single length scale h_c, was usedfor the stability parameter 3 = h_c/L in the height range 0.25 < z/h_c < 2.5. This choice is justified by mixing layer theory, which shows that within the roughness sublayer there is one dominant turbulence length scaledetermined by the degree of inflection in the windprofile at the canopy top. In the absence of theoretical or experimental evidence for guidance,standard Monin–Obukhov similarity functions, with = h_c/L, were used to calculate the stabilitydependence of _w(z) and _L(z) in the roughness sublayer. For z/h_c < 0.25 the turbulence length and time scales are influenced by the presence of the lowersurface, and stability effects are minimal. With theseassumptions there was excellent agreement between eddycovariance flux measurements and deductions from theinverse Lagrangian analysis. Stability correctionswere particularly necessary for night time fluxes whenthe atmosphere was stably stratified.The inverse Lagrangian analysis provides a useful toolfor testing and refining multilayer canopy models usedto predict radiation absorption, energy partitioningand CO₂ exchanges within the canopy and at thesoil surface. Comparison of model predictions withsource strengths deduced from the inverse analysisgave good results. Observed discrepancies may be dueto incorrect specification of the turbulent timescales and vertical velocity fluctuations close to theground. Further investigation of turbulencecharacteristics within plant canopies is required toresolve these issues.     

甘肃河东年降水量的周期变化

利用甘肃河东建站最早的９个代表站的年降水量资料,采用谐波分析和小波分析方法,计算了年降水量序列各种时间尺度的周期,以及１０ａ和２０ａ相空间时间序列的２～１０ａ周期试验;陇南无明显周期;陇南无明显周期。相空间序列的周期试验显示,甘肃中部的准３ａ周期主要反映在５０～７０年代,近２０～３０ａ则以８～９ａ周期较为明显;陇东近１０ａ主要表现为准２ａ振荡。进一步的分析表明,甘肃中部年降水量的准３ａ周期振荡目处     

半拉格朗日、半隐式欧拉方程组大气数值模式研究

介绍了一个具有较高时间积分效率的三维弹性大气数值模式。其中,控制方程为能够描述大气非静力及可压缩性的欧拉方程组,时间积分采用了高效率的“半拉格朗日、半隐式”方案。通过引用一个简化的“云物理过程”参数化方案,针对大气中深厚湿对流过程的若干问题进行了数值试验,结果表明,模式能够对该过程中的一些基本和复杂现象进行有效地模拟。试验结果还表明,动力学框架具有很好的稳定性,能够实现高效率的时间积分,它还具有较好的频散特性,保证了空间计算的精度,从而完好地描述了风暴的形态。     

Defect structure of the low temperature α-cristobalite phase and the cristobalite ⇆ tridymite transformation in (Si-Ge)O2

Using minimum exposure techniques, it is feasible to perform high resolution electron microscopy on the α-cristobalite phase of (Si_0.9 Ge_0.1)O₂, which is extremely radiation sensitive. Such images reveal atomic scale information of twins and tridymite-like stacking faults on (1 1 1)_β planes, as well as of domain boundaries resulting from the β→α transition. Polytype structures are formed in certain cases. Morphological features suggest that the phase transformation cristobalite → tridymite proceeds by means of a zonal dislocation mediated synchro-shear process on (1 1 1)_β planes; the geometry of this process is analyzed. Received: 13 June 1999 / Accepted: 30 October 1999     

Experimental determination of saltating glass particle dispersion in a turbulent boundary layer

A horizontal saltation layer of glass particles in air is investigated experimentally over a flat bed and also over a triangular ridge in a wind tunnel. Particle concentrations are measured by light scattering diffusion (LSD) and digital image processing, and velocities using particle image velocimetry (PIV). All the statistical moments of the particle concentration are determined such as mean concentration, root mean square concentration fluctuations, skewness and flatness coefficients. Over the flat bed, it is confirmed that the mean concentration decreases exponentially with height, the mean dispersion height being a significant length scale. It is shown that the concentration distribution follows quite well a lognormal distribution. Over the ridge, measurements were made at the top of the ridge and in the cavity region and are compared with measurements without the ridge. On the hill crest, particles are retarded, the saltation layer decreases in thickness and concentration is increased. Downwind of the ridge, particle flow behaves like a jet, in particular no particle return flow is observed. Copyright © 2006 John Wiley & Sons, Ltd.