The effect of expansion ratio on the critical Reynolds number in single fracture flow with sudden expansion

Flow in a single fracture (SF) is an important research subject in groundwater hydrology, hydraulic engineering, radioactive nuclear waste repository and geotechnical engineering. An abruptly changing aperture is a unique type of SF. This study discusses the relation between the values of the critical Reynolds number (Re_c) for the onset of symmetry breaking of flow and the expansion ratio (E) of SF, which is defined as the ratio between the outlet (D) and inlet (d) apertures. This study also investigates the effect of inlet aperture d on Re_c for flow in an SF with abruptly changing apertures (SF‐ACA) using the finite volume method. Earlier numerical and experimental results showed that flow is symmetric in respect to the central plane of the SF‐ACA at small Reynolds number (Re) but becomes asymmetric when Re is sufficiently large. Our simulations show that the value of Re_c decreases with the increasing E, and the relationship between the logarithm of Re_c and E can be described accurately using either a quadratic polynomial function or a logarithmic function. However, the relationship of Re_c and d for a given E value is vague, and Re_c becomes even less sensitive to d when E increases. This study also reveals that the hydraulic gradient (J) and flow velocity (v) follow a super‐linear relationship that can be fitted almost perfectly by the Forchheimer equation. The inertial component (J_i) of J increases monotonically with Re, whereas the viscous component (J_v) of J decreases monotonically with Re. The Re value corresponding to equal inertial and viscous components of J (named as the transitional point Re) decreases when E increases, and such a transitional point Re should be closely related to the critical Reynolds number Re_c, although a rigorous theoretical proof is not yet available. Copyright © 2015 John Wiley & Sons, Ltd.     

A new approach to estimate canopy evaporation and canopy interception capacity from evapotranspiration and sap flow measurements during and following wetting

Canopy interception and its evaporation into the atmosphere during irrigation or a rainfall event are important in irrigation scheduling, but are challenging to estimate using conventional methods. This study introduces a new approach to estimate the canopy interception from measurements of actual total evapotranspiration (ET) using eddy covariance and estimation of the transpiration from measurements of sap flow. The measurements were conducted over a small‐scale sprinkler‐irrigated cotton field before, during and after sprinkler irrigation. Evaporation and sap flow dynamics during irrigation show that the total ET during irrigation increased significantly because of the evaporation of free intercepted water while transpiration was suppressed almost completely. The difference between actual ET and transpiration (sap flow) during and immediately following irrigation (post irrigation) represents the total canopy evaporation while the canopy interception capacity was calculated as the difference between actual ET and transpiration (sap flow) during drying (post irrigation) following cessation of the irrigation. The canopy evaporation of cotton canopy was calculated as 0.8 mm, and the interception capacity was estimated to be 0.31 mm of water. The measurement uncertainty in both the non‐dimensional ET and non‐dimensional sap flow was shown to be very low. Copyright © 2015 John Wiley & Sons, Ltd.     

Parallelized combined finite-discrete element (FDEM) procedure using multi-GPU with CUDA

This paper focuses on the efficiency of finite discrete element method (FDEM) algorithmic procedures in massive computers and analyzes the time-consuming part of contact detection and interaction computations in the numerical solution. A detailed operable GPU parallel procedure was designed for the element node force calculation, contact detection, and contact interaction with thread allocation and data access based on the CUDA computing. The emphasis is on the parallel optimization of time-consuming contact detection based on load balance and GPU architecture. A CUDA FDEM parallel program was developed with the overall speedup ratio over 53 times after the fracture from the efficiency and fidelity performance test of models of in situ stress, UCS, and BD simulations in Intel i7-7700K CPU and the NVIDIA TITAN Z GPU. The CUDA FDEM parallel computing improves the computational efficiency significantly compared with the CPU-based ones with the same reliability, providing conditions for achieving larger-scale simulations of fracture.     

运用ENVI实现SPOT-5卫星影像的RVI提取

运用ENVI软件处理长沙地区的SPOT-5影像,将全色影像和多波段影像进行融合,分析融合后的影像,再选择感兴趣区域进行影像裁剪;采用ENVI的波段运算(Band3/Band2)提取该地区比值植被指数并进行分析。结果表明:比值植被指数对植被覆盖率较高的区域有非常高的敏感度;在植被覆盖率小于50%的时候,其敏感度明显降低,在绿色植被覆盖区域的比值植被指数远大于1,在裸露地表、建筑物、水体等没有植被覆盖区域的比值植被指数则在1附近,而比值植被指数大于2的则出现在高覆盖健康绿色植被上。     

The application of remote sensing data to diagnose soil degradation in the Dakhla depression – Western Desert,Egypt

Dakhla depression in Egypt’s Western Desert is experiencing two soil degradation processes, notably: soil salinization and sand encroachment. The present study aimed to diagnose the severity of these processes using remote sensing. Soil salinity was determined by spectral regression analysis between tasselled cap spectral transform extracted from a Landsat-8 image acquired in September 2013 along with synchronized soil salinity measurements. Assessment of sand advance rate was conducted by temporal change detection of brilliant crescentic sand dune visualized by Google Earth in old (2002) and recent (2013) images. Results showed that salinized soils (dS/m4<) represent 91% of bare lands and salinization is attributed to aridity, topography and poor drainage. Barchan dunes north and south of Abu Tartur escarpment moved at rates of 5.9 and 3.6 m/year, respectively. The escarpment protected the majority of the depression from massive dune invasion. However, sand encroachment is clearly observed west of the depression.     

A note on “Origin of components in Chilean thermal waters” by

Risacher et al. (2011) have presented voluminous data on thermal waters of Chile that is quite appreciable; however, their work still falls short on several counts. The most notable shortcoming of the work is the presentation and treatment of data. The interpretations are based on questionable premises (viz., extent of seawater intrusion) and considerations (viz., using average chemical composition of rock types for geochemical modeling, considering volcanic rocks as the only rock type in contact with the geothermal fluids) and assumptions not always substantiated by facts (with proper references) could have been corroborated. Use of Cl/Br ratio for discrimination purpose is unconvincing, considering the uncertainty in the measurement of low (<1 mg/L) Br concentration.     

直线型海洋重力仪负漂移原因分析及处理

分析了直线型KSS31M(SN035)海洋重力仪短期负漂移现象的原因、采取的处理措施和经验总结,认为在减震措施不足的情况下,重力仪长期经受震动,而频繁的搬动会使重力仪弹性系统中的金属弹簧产生较大的弹性疲劳,这种影响在短时间内无法完全消除,由此可能导致月漂移超标;陀螺稳定平台轴因承受震动磨损较严重,轴承间隙变大,导致安装在平台之上的重力传感器经过抛物线测试之后不能完全垂直指向地心,重力在水平方向存在一个微小的分量,最终导致重力读数偏小,即短期负漂移;重力基点环境的变化等对漂移也有一定的影响。建议在搬动重力仪后将仪器保持足够长时间的恒温后再开机,密切监测一段时间,以确保仪器处于良好的工作状态。     

延边开山屯花岗闪长岩-石英闪长岩体年代学、地球化学特征及对古太平洋板块俯冲作用时限的制约

近年来古太平洋构造域的构造演化备受学者关注。本文报道的延边开山屯地区花岗闪长岩-石英闪长岩体LAICP-MS U-Pb年龄表明其形成时间为早侏罗世早期(198±1Ma),所采样品可根据Zr/Hf值分为高Zr/Hf值组花岗闪长岩和低Zr/Hf值组石英闪长岩。高Zr/Hf值组花岗闪长岩起源深度浅,富集Rb、Th、U、K等大离子亲石元素(LILEs),贫Nb、Ta、Ti等高场强元素(HFSEs),具壳源岩浆的特点。低Zr/Hf值组为壳源岩浆与来自深部的亏损地幔岩浆混合而成,岩石亏损Nb、Ta、Zr、Hf、Ti等高场强元素,具有典型的弧型岩浆地球化学特征。岩体中存在细粒闪长质包体,镜下可见针柱状磷灰石。开山屯岩体属钙碱性系列岩石,结合前人资料,认为其与该地区同时代火成岩组成北-东向分布的早侏罗世活动大陆边缘型火成岩带,而位于该带西侧的小兴安岭-张广才岭地区存在同时代弧后拉张带,两者构成典型的大陆弧与弧后拉张带模型,共同揭示了早侏罗世早期古太平洋板块对东北地区的俯冲作用。     

黄海南部小黄鱼洄游分布及种群结构特征

Monthly fishery survey data of the small yellow croaker Larimichthys polyactis in the southern Yellow Sea from2003 and 2013 were employed to evaluate the variation in the resource distribution and biological characteristics(especially body length and sex ratio) in the population on a decadal scale.The results indicated that the small yellow croaker migrated from the Shawai fishing ground to the Dasha fishing ground in spring and was mainly distributed in the central and western parts of the Dasha fishing ground in April and May.Larimichthys polyactis in the Dasha fishing ground migrated eastward to offshore wintering grounds in autumn and reached the central Dasha fishing ground in October and November.The small yellow croaker entered the western waters of the Shawai fishing ground in winter.A large number of age 0+ fish occurred in the Shawai and Jiangwai fishing grounds in October of 2003 and 2013.The body lengths of the spawning stock and wintering stock in 2013 were larger than those in 2003,and the monthly sex ratios(female to male) were significantly less than 1 in both years.The monthly distribution of this fish in the southern Yellow Sea was consistent with a previous finding that "the stock migrated between the wintering grounds in the west of Jeju Island and the Lüsi spawning grounds" but tended to move more northward,with the spawning grounds extending outward.In the past decade,body length variation experienced a decline after an increase,rather than a steady decrease.The sex ratio in the single-stick stow net showed a tendency to increase over the decade,but was either less than or more than 1 depending on the fishing gear;therefore,further studies should be conducted to determine the sex structure.     

应用两次展开法求解系泊柱体慢漂运动方程

在深海中系泊的海洋平台,如Spar平台,水下部分为带有系泊的圆柱结构,其水平方向运动响应往往具有较低的自振频率,容易在低频波浪力(源于非线性的差频效应)作用下发生共振响应,使结构发生大幅水平慢漂。当浮体的瞬时位置大幅偏离初始位置时,基于初始平衡位置的摄动展开法会存在较大误差。针对这一问题,采用两次展开方法,对大幅慢漂运动开展时域模拟研究。对双色波作用下自由漂浮圆柱的大幅运动响应问题进行数值分析,并与采用基于初始平衡位置的摄动展开法的计算结果进行了对比。结果表明,采用新的两次展开法可以计算出波浪遭遇频率的变化和波浪漂移阻尼,而这无法从基于初始平衡位置的常规摄动展开法中得到,体现了两次展开法在分析大幅慢漂问题上的优势。