基于辅助变量的森林半腐层厚度空间插值精度

基于地统计学方法,利用3种以海拔作为辅助变量的空间插值算法[局部平均的简单kriging法(simple kriging with locally varying mean,SKlm)、带有外部漂移的kriging法(krig-ing with an external drift,KED)和协kriging法(cokrging,COK)]计算了森林半腐层厚度的空间插值精度,并进行了交叉验证.结果表明:KED法既考虑了变量之间的空间变异,又考虑到影响局部空间变化的因素,与其他插值方法相比,其精度有很大提高;由于海拔与半腐层厚度之间的相关关系较弱,导致SKlm法的插值精度没有达到预期效果;COK法直接将海拔用于估计半腐层厚度,由于在边界地区缺乏采样点数据,因此边界地区的插值出现了多处突变区域.对比地统计学方法与距离反比权重法(inverse distance weighting,IDW)在本研究中的插值精度,除了KED方法的插值精度较高外,其余方法的插值精度均不及IDW,原因可能是利用辅助变量辅助地统计学插值时,主、辅助变量之间的相关关系在插值中起着重要作用.

Simple lattice Boltzmann model for traffic flows

A lattice Boltzmann model with 5-bit lattice for traffic flows is proposed. Using the Chapman-Enskog expansion and multi-scale technique, we obtain the higher-order moments of equilibrium distribution function. A simple traffic light problem is simulated by using the present lattice Boltzmann model, and the result agrees well with analytical solution.

A simple approach to the evaluation of fiber/matrix interfacial shear strength and fracture toughness

On the basis of microscopical analyses of the fiber distribution and longitudinal shear deformation in unidirectional fiber composites, a simple approach is presented for characterizing the interfacial sheer strength and fracture toughness.

On the Evolution of Simple Shear in Saturated Soils

Development of shear bands in saturated soils is a multi-stage process based on the theoretical and numerical investigations in this paper. The soil is initially in homogenous shear strain state, and the instability can be characterized by a dimensionless number D. The inhomogenous distribution of shear strains appears when D>1, and the shear band will initiate and develop gradually. Numerical solutions show that only single shear band that is finally formed in the central region of the specimen even several disturbances (distributed along the specimen) appear in the beginning.

Simple phenomenological determination of contact stiffness and elastic modulus of Ce-based bulk metallic glasses through nanoindentation

The viscoelastic deformation of Ce-based bulk metallic glasses (BMGs) with low glass transition temperature is investigated at room temperature. Contact stiffness and elastic modulus of Ce-based BMGs cannot be derived using the conventional Oliver-Pharr method [W. C. Oliver and G. M. Pharr, J. Mater. Res. 7, 1564 (1992)]. The present work shows that the time dependent displacement of unloading segments can be described well by a generalized Kelvin model. Thus, a modified Oliver-Pharr method is proposed to evaluate the contact stiffness and elastic modulus, which does, in fact, reproduce the values obtained via uniaxial compression tests. (c) 2007 American Institute of Physics.

A simple approach to solve boundary-value problems in gradient elasticity

We outline a procedure for obtaining solutions of certain boundary value problems of a recently proposed theory of gradient elasticity in terms of solutions of classical elasticity. The method is applied to illustrate, among other things, how the gradient theory can remove the strain singularity from some typical examples of the classical theory.

Thermo-plastic instability in simple shear

A theoretical description of thermo-plastic instability in simple shear is presented in a system of equations describing plastic deformation, the first law of thermodynamics and Fourier's heat transfer rule. Both mechanical and thermodynamical parameters influence instability and it is shown that two different modes of instability may exist. One of them is dominated by thermal softening and has a characteristic time and length, connected to each other by thermal diffusion.A criterion combining thermal softening, current stress, density, specific heat, work-hardening, thermal conductivity and current strain rate is obtained and practical implications are discussed.

非结构同位网格SIMPLE类算法收敛性能比较

以非结构网格的SIMPLE算法为基础,将算法扩展为SIMPLEC算法.利用30°角的斜方腔流动计算成果,分析了非结构同位网格的SIMPLE/SIMPLEC算法的收敛性能;比较了因网格的非正交而引入的非正交项的取舍对该算法收敛性能的影响;并采用显式校正步法对SIMPLEC算法进行了显式校正.比较表明,在非结构同位网格SIMPLEC算法中可忽略非正交项,但有必要对压力作亚松弛.显式校正步法可显著地加速在非结构网格上求解N-S方程的收敛性能,而且在不同的松弛因子组合下,均有较好的收敛速率.

基于非结构化同位网格的SIMPLE算法

通过基于非结构化网络的有限体积法对二维稳态Navier-Stokes方程进行了数值求解。其中对流项采用延迟修正的二阶格式进行离散；扩散项的离散采用二阶中心差分格式；对于压力-速度耦合利用SIMPLE算法进行处理；计算节点的布置采用同位网格技术，界面流速通过动量插值确定。本文对方腔驱动流、倾斜腔驱动流和圆柱外部绕流问题进行了计算，讨论了非结构化同位网格有限体积法在实现SIMPLE算法时，迭代次数与欠松弛系数的关系、不同网格情况的收敛性、同结构化网格的对比以及流场尾迹结构。通过和以往结果比较可知，本文的方法是准确和可信的。

Experiments about Simple Liquid Flows in Microtubes

The behavior of micro-scale flow is significant for the performance of Micro-Electro-Mechanical- Systems (MEMS) devices. Some experiments about liquid flow through microtubes with diameters about 3similar to20mum are presented here. The liquids used in our experiments include some simple liquids with small molecules, such as non-ion water and several kinds of organic liquids (CCL4, C6H5C2H5 and Isopropanol etc.). The flow rate and the normalized friction cocfficients were measured in micro-flow experimental apparatus. The results show that when the driven pressure varies from 0 to 1Mpa, the flow behaviors in 20mum microtube for both polar and non-polar liquids are in agreement with Hagen-Poiseuille law of the classical theory. It means that N-S equation based on continuous medium still acts well in this case. For higher pressure drop from 1 to 30Mpa, in the microtubes with diameter of 3similar to10mum, the normalized friction coefficients of organic liquids can't keep constant with pressure increases. However the non-ion water reveals different trends.

A simple method to simulate shrinkage-induced cracking in cement-based composites by lattice-type modeling

A new numerical procedure is proposed to investigate cracking behaviors induced by mismatch between the matrix phase and aggregates due to matrix shrinkage in cement-based composites. This kind of failure processes is simplified in this investigation as a purely spontaneous mechanical problem, therefore, one main difficulty during simulating the phenomenon lies that no explicit external load serves as the drive to propel development of this physical process. As a result, it is different from classical mechanical problems and seems hard to be solved by using directly the classical finite element method (FEM), a typical kind of "load -> medium -> response" procedures. As a solution, the actual mismatch deformation field is decomposed into two virtual fields, both of which can be obtained by the classical FEM. Then the actual response is obtained by adding together the two virtual displacement fields based on the principle of superposition. Then, critical elements are detected successively by the event-by-event technique. The micro-structure of composites is implemented by employing the generalized beam (GB) lattice model. Numerical examples are given to show the effectiveness of the method, and detailed discussions are conducted on influences of material properties.

A Simple Model For Predicting The Void Fraction Of Gas/Non-Newtonian Fluid Intermittent Flows In Upward Inclined Pipes

In this work, a simple correlation, which incorporates the mixture velocity, drift velocity, and the correction factor of Farooqi and Richardson, was proposed to predict the void fraction of gas/non-Newtonian intermittent flow in upward inclined pipes. The correlation was based on 352 data points covering a wide range of flow rates for different CMC solutions at diverse angles. A good agreement was obtained between the predicted and experimental results. These results substantiated the general validity of the model presented for gas/non-Newtonian two-phase intermittent flows.

Pure Coulomb explosions of highly charged methane clusters investigated by a simple electrostatic model

The pure Coulomb explosions of the methane clusters (CA(4))(n), (light atom A = H or D) have been investigated by a simplified electrostatic model for both a single cluster and an ensemble of clusters with a given cluster size distribution. The dependence of the energy of ions produced from the explosions on cluster size and the charge state of the carbon ions has been analysed. It is found that, unlike the average proton energy which increases with the charge q of the carbon ions, the average deuteron energy tends to saturate as q becomes larger than 4. This implies that when the laser intensity is sufficiently high for the (CD4)(n) to be ionized to a charge state of (C4+D4+)(n), the neutron yield from a table-top laser-driven Coulomb explosion of deuterated methane clusters (CD4)(n) could be increased significantly by increasing the interaction volume rather than by increasing the laser intensity to produce the higher charge state (C6+D4+)(n). The flight-time spectra of the carbon ions and the light ions have also been studied.

Simple ABCD matrix method for evaluating optical coupling system of laser diode to single-mode fiber with a lensed-tip

In this paper, we present a simple technique to determine the coupling efficiency between a laser diode and a lensed-tip based on the ABCD transformation matrix method. We have compared our analysis technique to that of previous work and have found that the presented method is reliable in predicting the coupling efficiency of lensed-tip and has the advantage of simplicity of coupling efficiency calculation even by a pocket calculator. The results can be useful for designing coupling optics. (c) 2005 Elsevier GmbH. All rights reserved.