On the stability of stratified plane couette flow

Abstract

Unbounded stratified plane Couette flow is shown to be stable against small amplitude disturbances. The Brunt-Väisälä frequency is assumed to be constant. Both viscosity and thermal diffusion are included, and shown to be stabilizing.     

Coupled finite‐element simulation of injection well testing in unconsolidated oil sands reservoir

This paper presents a finite‐element (FE) model for simulating injection well testing in unconsolidated oil sands reservoir. In injection well testing, the bottom‐hole pressure (BHP) is monitored during the injection and shut‐in period. The flow characteristics of a reservoir can be determined from transient BHP data using conventional reservoir or well‐testing analysis. However, conventional reservoir or well‐testing analysis does not consider geomechanics coupling effects. This simplified assumption has limitations when applied to unconsolidated (uncemented) oil sands reservoirs because oil sands deform and dilate subjected to pressure variation. In addition, hydraulic fracturing may occur in unconsolidated oil sands when high water injection rate is used. This research is motivated in numerical modeling of injection well testing in unconsolidated oil sands reservoir considering the geomechanics coupling effects including hydraulic fracturing. To simulate the strong anisotropy in mechanical and hydraulic behaviour of unconsolidated oil sands induced by fluid injection in injection well testing, a nonlinear stress‐dependent poro‐elasto‐plastic constitutive model together with a strain‐induced anisotropic permeability model are formulated and implemented into a 3D FE simulator. The 3D FE model is used to history match the BHP response measured from an injection well in an oil sands reservoir. Copyright © 2013 John Wiley & Sons, Ltd.     

Nonlinear analysis of laterally loaded rigid piles in cohesive soil

This article presents a method for the nonlinear analysis of laterally loaded rigid piles in cohesive soil. The method considers the force and the moment equilibrium to derive the system equations for a rigid pile under a lateral eccentric load. The system equations are then solved using an iteration scheme to obtain the response of the pile. The method considers the nonlinear variation of the ultimate lateral soil resistance with depth and uses a new closed‐form expression proposed in this article to determine the lateral bearing factor. The method also considers the horizontal shear resistance at the pile base, and a bilinear relationship between the shear resistance and the displacement is used. For simplicity, the modulus of horizontal subgrade reaction is assumed to be constant with depth, which is applicable to piles in overconsolidated clay. The nonlinearity of the modulus of horizontal subgrade reaction with pile displacement at ground surface is also considered. The validity of the developed method is demonstrated by comparing its results with those of 3D finite element analysis. The applications of the developed method to analyze five field test piles also show good agreement between the predictions and the experimental results. The developed method offers an alternative approach for simple and effective analysis of laterally loaded rigid piles in cohesive soil. Copyright © 2011 John Wiley & Sons, Ltd.     

Garnet crystal plasticity in the continental crust,new example from south Madagascar

Garnet (10 vol.%; pyrope contents 34–44 mol.%) hosted in quartzofeldspathic rocks within a large vertical shear zone of south Madagascar shows a strong grain‐size reduction (from a few cm to ～300 μm). Electron back‐scattered diffraction, transmission electron microscopy and scanning electron microscope imaging coupled with quantitative analysis of digitized images (PolyLX software) have been used in order to understand the deformation mechanisms associated with this grain‐size evolution. The garnet grain‐size reduction trend has been summarized in a typological evolution (from Type I to Type IV). Type I, the original porphyroblasts, form cm‐sized elongated grains that crystallized upon multiple nucleation and coalescence following biotite breakdown: biotite + sillimanite + quartz = garnet + alkali feldspar + rutile + melt. These large garnet grains contain quartz ribbons and sillimanite inclusions. Type I garnet is sheared along preferential planes (sillimanite layers, quartz ribbons and/or suitably oriented garnet crystallographic planes) producing highly elongated Type II garnet grains marked by a single crystallographic orientation. Further deformation leads to the development of a crystallographic misorientation, subgrains and new grains resulting in Type III garnet. Associated grain‐size reduction occurs via subgrain rotation recrystallization accompanied by fast diffusion‐assisted dislocation glide. This plastic deformation of garnet is associated with efficient recovery as shown by the very low dislocation densities (10¹⁰ m^?3 or lower). The rounded Type III garnet experiences rigid body rotation in fine‐grained matrix. In the highly deformed samples, the deformation mechanisms in garnet are grain‐size‐ and shape‐dependent: dislocation creep is dominant for the few large grains left (>1 mm; Type II garnet), rigid body rotation is typical for the smaller rounded grains (300 μm or less; Type III garnet) whereas diffusion creep may affect more elliptic garnet (Type IV garnet). The P–T conditions of garnet plasticity in the continental crust (≥950 °C; 11 kbar) have been identified using two‐feldspar thermometry and GASP conventional barometry. The garnet microstructural and deformation mechanisms evolution, coupled with grain‐size decrease in a fine‐grained steady‐state microstructure of quartz, alkali feldspar and plagioclase, suggests a separate mechanical evolution of garnet with respect to felsic minerals within the shear zone.     

Rate and historical change of direct recharge from precipitation constrained by unsaturated zone profiles of chloride and oxygen‐18 in dry river bed of North China Plain

Quantifying of direct recharge derived from precipitation is crucial for assessing sustainability of well‐irrigated agriculture. In the North China Plain, the land use is dominated by groundwater‐irrigated farmland where the direct recharge derived from precipitation and irrigation. To characterize the mean rate and historical variance of direct recharge derived from precipitation, unsaturated zone profiles of chloride and δ¹⁸O in the dry river bed of the Beiyishui River were employed. The results show that archival time scale of the profile covers the duration from 1980 to 2002 (corresponding to depths from 5 to 2 m) which is indicated by matching the δ¹⁸O peaks in the isotope profile with the aridity indexes gained by instrumental records of annual precipitation and annual potential evaporation. Using the chloride mass balance method, the mean rate of the direct recharge corresponding to the archival time scale is estimated to be 3·8 ± 0·8 mm year^?1, which accounts for about 0·7% of the long‐term average annual precipitation. Further, the direct recharge rates vary from 2·1 to 6·8 mm year^?1 since 1980. Despite the subhumid climate, the estimate of recharge rates is in line with other findings in semiarid regions. The low rate of direct recharge is considered as a result of the relative dry climate in recent decades. In dry river bed, unsaturated zone profiles of chloride and δ¹⁸O combined with instrumental records could offer valuable information about the direct recharge derived from precipitation during droughts. Copyright © 2011 John Wiley & Sons, Ltd.     

The role of eddies in solute transport and recovery in rock fractures: Implication for groundwater remediation

A better understanding of solute transport and retention mechanism in rock fractures has been challenging due to difficulty in their direct observations in microscale rough‐walled fractures. Six representative troughs in a rough‐walled fracture were selected for microscale observations of eddy formation with increasing flow velocity and its effect on spatiotemporal changes of solute concentration. This experimental study was enabled by a microscale visualization technique of micro particle image velocimetry. With increasing flow velocity (Re ≤ 2.86), no eddies were generated, and solutes along the main streamlines transported rapidly, whereas those near the wall moved slowly. A larger amount of solutes remained trapped at all troughs at Re = 2.86 than Re < 1. For Re = 8.57, weak eddies started to be developed at the troughs on the lee side, which little contributed to overall solute flushing in the fracture. Accordingly, a large of amount of water was needed for solute flushing. The flow condition of 1 < Re < 10, before a full development of eddies, was least favourable in terms of time and amount of remediation fluid required to reach a target concentration. After large eddies were fully developed at troughs on the lee side for Re = 17.13, solutes were substantially reduced by eddies with less amount of water. Fully developed eddies were found to enhance solute transport and recovery, as opposed to a general consensus that eddies trap and delay solutes. Direct inflow into troughs on the stoss side also made a great contribution to solute flushing out of the troughs. This study indicates that fully developed eddies or strong inflows at troughs are highly possible to form for Re > 10 and this flow range could be favourable for efficient remediation.     

A mixing-length model for predicting vertical velocity distribution in flows through emergent vegetation

Abstract

The vertical profiles of streamwise velocities are computed on flood plains vegetated with trees. The calculations were made based on a newly developed one-dimensional model, taking into account the relevant forces acting on the volumetric element surrounding the considered vegetation elements. A modified mixing length concept was used in the model. An important by-product of the model is the method for evaluating the friction velocities, and consequently bed shear stresses, in a vegetated channel. The model results were compared with the relevant experimental results obtained in a laboratory flume in which flood plains were covered by simulated vegetation.     

Late-Hercynian shearing during crystallization of granitoid magmas (Sila massif,southern Italy): regional implications

Shearing of regional extent, involving granitoids and underlying mid-crustal rocks of the Sila massif (Calabria, Italy), is analysed in this paper. The deformed granitoids are affected by a wide NNW-SSE oriented deformation zone, stretching for about 60 km, from the neighbourhood of Cecita Lake to Cropani village. Meso- and micro-structures in granitoids, close to the boundary with underlying migmatitic paragneiss, indicate that deformation developed from melt-present to solid-state conditions. Simultaneous tectonics and magmatism activated a plutonic accretionary process at mid-crustal levels. This took place at about 300 Ma and involved hybrid magmas with a dominat contribution from a mantle source. The deformation regime remained steady for a long time during magma crystallization and cooling in subsolidus conditions. The regional top-to-the-W sense of shear in the present geographic coordinates, recorded in the deformed granitoids, seems geometrically consistent with the coeval direction of maximum extension found in another sector of the southern Hercynian belt, suggesting the original position of the Sila basement in this context. Magmatic ativity ended with the intrusion of mafic and felsic magams affected by a very weak deformation, ongoing during the final strain increments of the late-Hercynian stage.     

Characterising deep vadose zone water movement and solute transport under typical irrigated cropland in the North China Plain

Understanding the dynamics and mechanisms of soil water movement and solute transport is essential for accurately estimating recharge rates and evaluating the impacts of agricultural activities on groundwater resources. In a thick vadose zone (0–15 m) under irrigated cropland in the piedmont region of the North China Plain, soil water content, matric potential, and solute concentrations were measured. Based on these data, the dynamics of soil water and solutes were analysed to investigate the mechanisms of soil water and solute transport. The study showed that the 0–15‐m vadose zone can be divided into three layers: an infiltration and evaporation layer (0–2 m), an unsteady infiltration layer (2–6 m), and a quasi‐steady infiltration layer (6–15 m). The chloride, nitrate, and sulphate concentrations all showed greater variations in the upper soil layer (0–1 m) compared to values in the deep vadose zone (below 2 m). The average concentrations of these three anions in the deep vadose zone varied insignificantly with depth and approached values of 125, 242, and 116 mg/L. The accumulated chloride, sulphate, and nitrate were 2,179 ± 113, 1,760 ± 383, and 4,074 ± 421 kg/ha, respectively. The soil water potential and solute concentrations indicated that uniform flow and preferential flow both occurred in the deep vadose zone, and uniform flow was the dominant mechanism of soil water movement in this study. The piston‐like flow velocity of solute transport was 1.14 m per year, and the average value of calculated leached nitrate nitrogen was 107 kg/ha?year below the root zone. The results can be used to better understand recharge processes and improve groundwater resources management.     

宽窄相间河道水流紊动特性试验研究

为了探索宽窄相间河道的水流紊动特性，以西南地区宝兴河上游宽窄相间河段为研究对象，基于室内概化模型试验，采用多普勒声学流速仪（ADV）测量了室内模型典型断面上的三维瞬时流速，分析典型断面上的纵向时均流速、紊动强度、雷诺切应力和紊动能的分布规律。试验结果显示：宽窄相间水槽中，扩散段边壁的紊动强度大于中心区域的紊动强度，最大值位于0.2倍水深处；扩散段两侧坡脚处紊动能最大；侧壁区的平面和立面雷诺切应力最大值出现在扩散段内，中心区域最大雷诺切应力位于两槽间的中间断面处；扩散段内水流紊乱，两侧出现旋涡和涡脱，易造成侧壁侵蚀加强，引起河道拓宽。深入分析了宽窄相间河道水流的紊动特性，可为山区河流治理和自然灾害防治提供参考。