Statistical characteristics of flow as indicators of channeling in heterogeneous porous and fractured media

We introduce two new channeling indicators D_ic and D_cc based on the Lagrangian distribution of flow rates. On the basis of the participation ratio, these indicators characterize the extremes of both the flow-tube width distribution and the flow rate variation along flow lines. The participation ratio is an indicator biased toward the larger values of a distribution and is equal to the normalized ratio of the square of the first-order moment to the second-order moment. Compared with other existing indicators, they advantageously provide additional information on the flow channel geometry, are consistently applicable to both porous and fractured media, and are generally less variable for media generated using the same parameters than other indicators. Based on their computation for a broad range of porous and fracture permeability fields, we show that they consistently characterize two different geometric properties of channels. D_ic gives a characteristic scale of low-flow zones in porous media and a characteristic distance between effectively flowing structures in fractured cases. D_cc gives a characteristic scale of the extension of high-flow zones in porous media and a characteristic channel length in fractured media. D_ic is mostly determined by channel density and permeability variability. D_cc is, however, more affected by the nature of the correlation structure like the presence of permeability channels or fractures in porous media and the length distribution in fracture networks.     

An adaptive approach to selecting a flow‐partition exponent for a multiple‐flow‐direction algorithm

Most multiple‐flow‐direction algorithms (MFDs) use a flow‐partition coefficient (exponent) to determine the fractions draining to all downslope neighbours. The commonly used MFD often employs a fixed exponent over an entire watershed. The fixed coefficient strategy cannot effectively model the impact of local terrain conditions on the dispersion of local flow. This paper addresses this problem based on the idea that dispersion of local flow varies over space due to the spatial variation of local terrain conditions. Thus, the flow‐partition exponent of an MFD should also vary over space. We present an adaptive approach for determining the flow‐partition exponent based on local topographic attribute which controls local flow partitioning. In our approach, the influence of local terrain on flow partition is modelled by a flow‐partition function which is based on local maximum downslope gradient (we refer to this approach as MFD based on maximum downslope gradient, MFD‐md for short). With this new approach, a steep terrain which induces a convergent flow condition can be modelled using a large value for the flow‐partition exponent. Similarly, a gentle terrain can be modelled using a small value for the flow‐partition exponent. MFD‐md is quantitatively evaluated using four types of mathematical surfaces and their theoretical ‘true’ value of Specific Catchment Area (SCA). The Root Mean Square Error (RMSE) shows that the error of SCA computed by MFD‐md is lower than that of SCA computed by the widely used SFD and MFD algorithms. Application of the new approach using a real DEM of a watershed in Northeast China shows that the flow accumulation computed by MFD‐md is better adapted to terrain conditions based on visual judgement.     

Small-scale patterns of meiofauna in a bryophyte covered tufa barrier (Plitvice Lakes,Croatia)

Meiobenthic fauna is complex and essential part of the stream benthos and it is known as intermediaries from bacteria to higher trophic levels. However, the dynamics and mechanisms regulating this community at small-scale largely have been neglected. This study was carried out to determine meiofaunal dynamics in a small scale-patterns influenced by flow velocity and other abiotic and biotic factors. We examined meiofauna within bryophytes on a tufa barrier in Plitvice Lakes National Park using a detailed taxonomic approach of various meiofaunal groups and their feeding guilds. We choose three microhabitats (slow, medium, fast) differing significantly in flow velocity above bryophytes. Bdelloid rotifers were the most abundant group in microhabitat with highest flow velocity, while in other two microhabitats nematodes and monogonont rotifers prevailed in abundance. Data on environmental variables and main meiofaunal taxa and feeding guilds were analyzed using redundancy analysis. This analysis indicated that microfilter feeding guild (e.g. bdelloid rotifers) was strongly affected by interaction of flow velocity and POM fractions. Other feeding guilds were influenced by temperature, oxygen and/or pH and did not prefer high flow velocity. Suction-feeder nematodes and microfilter-feeder rotifers were dominant on temporal and spatial scale, indicating their good adaptations on frequently disturbed conditions that prevailed on bryophyte covered tufa barrier. Our results provide comprehensive survey of diversity, density as well as trophic structure of meiofauna in aquatic bryophytes. Differences in meiofaunal composition and density between three microhabitats suggest that the meiofauna is relevant indicator of environmental changes even at small-scale pattern.     

Subharmonic Dynamo Action in the Roberts Flow

The paper deals with the dynamo action of the Roberts flow, that is, a flow depending periodically on two cartesian coordinates, X and Y , but being independent of the third one, Z . In particular the case is considered in which the magnetic fields, which are periodic in X, Y and Z , have period lengths in the XY -plane being integer multiples of that of the flow. Two approaches are used. Firstly, the equations governing the magnetic field are reduced to a matrix eigenvalue problem, which is solved numerically. Secondly, a mean magnetic field is defined by averaging over proper areas in the XY -plane, corresponding equations are derived, in which the induction effect of the flow occurs as an anisotropic f -effect, and analytic solutions are given. The results are of particular interest for the Karlsruhe dynamo experiment, which works with a Roberts type flow consisting of 52 cells inside a cylindrical volume. In order to check the reliability of predictions concerning self-excitation based on the mean-field approach, analogous predictions are derived for a rectangular box containing 50 cells, and are compared with results obtained with the help of direct solutions of the eigenvalue problem mentioned. It turns out that the simple mean-field approach in general underestimates the requirements for self-excitation. The corresponding results agree with those obtained in the subharmonic approach only if the side length L of the box, its height H and the edge length l of a spin generator satisfy $ L \gg H \gg l $ . In Appendix B, some comments on previous results concerning $\cal {ABC}$ dynamos are made in the light of the subharmonic formalism used in the paper.     

Dynamic p-enrichment schemes for multicomponent reactive flows

We present a family of p-enrichment schemes. These schemes may be separated into two basic classes: the first, called fixed tolerance schemes, rely on setting global scalar tolerances on the local regularity of the solution, and the second, called dioristic schemes, rely on time-evolving bounds on the local variation in the solution. Each class of p-enrichment scheme is further divided into two basic types. The first type (the Type I schemes) enrich along lines of maximal variation, striving to enhance stable solutions in “areas of highest interest.” The second type (the Type II schemes) enrich along lines of maximal regularity in order to maximize the stability of the enrichment process. Each of these schemes are tested on three model systems. The first is an academic exact system where basic analysis is easily performed. Then we discuss a pair of application model problems arising in coastal hydrology. The first being a contaminant transport model, which addresses a declinature problem for a contaminant plume with respect to a bay inlet setting. And the second, a multicomponent chemically reactive flow model of estuary eutrophication arising in the Gulf of Mexico.     

基于物质流分析的区域循环经济评价指标体系构建——以陕西省榆林市为例

以物流分析为基础,兼顾循环经济的3R原则,利用循环经济发展中涉及的物质流分析指标,以及物质流指标与人、经济、单位投入、产出等外延领域的结合衍生指标构建循环经济评价指标体系,并对各指标作了解释,对数据来源途径做了分析.     

基于水电比拟方法的山东南四湖上级湖流场研究

利用水电比拟方法模拟了南四湖流场结构。同时,基于HAMSOM模式,对其流场进行了数值模拟,所得结果与水电比拟结果一致,两者平均相关系数在0.94以上,表明水电比拟实验技术路线可行,流场基本结构准确。该方法适用于低流速大流场的模拟。流场结构分析表明,湖区流速在汛期为cm/s量级,非汛期仅为mm/s量级。其中,在南阳湖北部...     

VIV and galloping of single circular cylinder with surface roughness at 3.0×10≤Re≤1.2×10

Passive Turbulence Control (PTC) in the form of selectively distributed surface roughness is used to alter Flow Induced Motion (FIM) of a circular cylinder in a steady flow. The objective is to enhance FIM's synchronization range and amplitude, thus maximizing conversion of hydrokinetic energy into mechanical energy by oscillator in vortex-induced vibration and/or galloping. Through additional viscous damping, mechanical energy is converted to electrical harnessing clean and renewable energy from ocean/river currents. High Reynolds numbers (Re) are required to reach the high-lift TrSL3 (Transition-Shear-Layer-3) flow regime. PTC trips flow separation and energizes the boundary layer, thus inducing higher vorticity and consequently lift. Roughness location, surface coverage, and size are studied using systematic model tests with broad-field laser visualization at 3.0×10⁴<Re<1.2×10⁵ in the low-turbulence free-surface water-channel of the Marine Renewable Energy Laboratory of the University of Michigan. Test results show that 16° roughness coverage is effective in the range (10°-80°) inducing reduced vortex-induced vibration (VIV), enhanced VIV, or galloping. Range of synchronization may increase or decrease, galloping amplitude of oscillation reaches three diameters; wake structures change dramatically reaching up to ten vortices per cycle. Conversion of hydrokinetic energy to mechanical is enhanced strongly with proper PTC.     

New model to determine forces at on-bottom slender pipelines

The present paper proposes a numerical model to determine horizontal and vertical components of the hydrodynamic forces on a slender submarine pipeline lying at the sea bed and exposed to non-linear waves plus a current. The new model is an extension of the Wake II type model, originally proposed for sinusoidal waves (Soedigdo et al., 1999) and for combined sinusoidal waves and currents (Sabag et al., 2000), to the case of periodic or random waves, even with a superimposed current. The Wake II type model takes into account the wake effects on the kinematic field and the time variation of drag and lift hydrodynamic coefficients. The proposed extension is based on an evolutional analysis carried out for each half period of the free stream horizontal velocity at the pipeline. An analytical expression of the wake velocity is developed starting from the Navier–Stokes and the boundary layer equations. The time variation of the drag and lift hydrodynamic coefficients is obtained using a Gaussian integration of the start-up function. A reduced scale laboratory investigation in a large wave flume has been conducted in order to calibrate the empirical parameters involved in the proposed model. Different wave and current conditions have been considered and measurements of free stream horizontal velocities and dynamic pressures on a bottom-mounted pipeline have been conducted. The comparison between experimental and numerical hydrodynamic forces shows the accuracy of the new model in evaluating the time variation of peaks and phase shifts of the horizontal and vertical wave and current induced forces.