Simulation of Ekman Boundary Layers by Large Eddy Model with Dynamic Mixed Subfilter Closure

Theoretical analysis of boundary layer turbulence has suggested a feasibility of sufficiently accurate turbulence resolving simulations at relatively coarse meshes. However, large eddy simulation (LES) codes, which employ traditional eddy-viscosity turbulence closures, fail to provide adequate turbulence statistics at coarse meshes especially within a surface layer. Manual tuning of parameters in these turbulence closures may correct low order turbulence statistics but severely harms spectra of turbulence kinetic energy (TKE). For more than decade, engineering LES codes successfully employ dynamic turbulence closures. A dynamic Smagorinsky turbulence closure (DSM) has been already tried in environmental LES. The DSM is able to provide adequate turbulence statistics at coarse meshes but it is not completely consistent with the LES equations. This paper investigates applicability of an advanced dynamic mixed turbulence closure (DMM) to simulations of Ekman boundary layers of high Reynolds number flows. The DMM differs from the DSM by explicit calculation of the Leonard term in the turbulence stress tensor. The Horizontal Array Turbulence Study (HATS) field program has revealed that the Leonard term is indeed an important component of the real turbulence stress tensor. This paper presents validation of a new LES code LESNIC. The study shows that the LES code with the DMM provides rather accurate low order turbulence statistics and the TKE spectra at very coarse meshes. These coarse LES maintain more energetic small scale fluctuations of velocity especially within the surface layer. This is critically important for success of simulations. Accurate representation of higher order turbulence statistics, however, requires essentially better LES resolution. The study also shows that LES of the Ekman boundary layer cannot be directly compared with conventionally neutral atmospheric boundary layers. The depth of the boundary layer is an important scaling parameter for turbulence statistics.     

Measurements and modeling of open-channel flows with finite semi-rigid vegetation patches

The hydrodynamics of flows through a finite length semi-rigid vegetation patch (VP) were investigated experimentally and numerically. Detailed measurements have been carried out to determine the spatial variation of velocity and turbulence profiles within the VP. The measurement results show that an intrusion region exists in which the peak Reynolds stress remains near the bed. The velocity profile is invariant within the downstream part of the VP while the Reynolds stress profile requires a longer distance to attain the spatially invariant state. Higher vegetation density leads to a shorter adjustment length of the transition region, and a higher turbulence level within the VP. The vegetation density used in the present study permits the passing through of water and causes the peak Reynolds stress and turbulence kinetic energy each the maximum at the downstream end of the patch. A 3D Reynolds-averaged Navier–Stokes model incorporating the Spalart–Allmaras turbulence closure was employed subsequently to replicate the flow development within the VP. The model reproduced transitional flow characteristics well and the results are in good agreement with the experimental data. Additional numerical experiments show that the adjustment length can be scaled by the water depth, mean velocity and maximum shear stress. Empirical equations of the adjustment lengths for mean velocity and Reynolds stress were derived with coefficients quantified from the numerical simulation results.     

Stratification and Circulation in a Shallow Turbid Waterbody

Shallow waterbodies are often assumed to be well mixed in the vertical. However, when they are characterised by high turbidity levels, absoption of solar heating within a relatively thin surface layer can produce thermal stratification. Results from an intensive monitoring program have been combined with three-dimensional circulation modelling to examine the diurnal stratification cycles in a small turbid waterbody. The waterbody, known as Rushy Billabong, is located in southeastern Australia and its high turbidity coupled with forcing by wind and solar radiation resulted in regular diurnal cycles of stratification and overturning. Under conditions of light wind and high solar radiation, the model results were generally consistent with the observed temperature field. However, under stronger winds, preferential cooling and sinking of shallow water around the edge of the lake began to contribute significantly to the interior stratification. Model estimates then became more sensitive to the detailed bathymetry and the choice of turbulence parameterisation. The level of stratification is also shown to influence the circulation in the billabong by trapping the wind-driven flow near the surface. Insights provided by the observations and modelling may have broader implications for the management of small turbid systems such as settling ponds, aquaculture ponds, and some natural wetlands.     

Analytical model for vertical velocity profiles in flows with submerged shrub-like vegetation

An analytical solution for the vertical profiles of the horizontal velocity of channel flow with submerged shrub-like vegetation is investigated in this paper. At first, a shape function is proposed to fit the diameter change of different types of shrub-like vegetation. Using the momentum theorem and the mixing-length turbulence model, an analytical solution for the vertical profile of the horizontal velocity within the vegetation is obtained. The velocity distribution of the whole column is determined in tandem with the logarithmic velocity profile above the vegetation. The solution is compared with experimental data in excellent agreement. The results show that the flow above the vegetation has a logarithmic velocity profile while the flow within the vegetation is impacted greatly by the shape and density of vegetation. The flows within shrub-like vegetations are non-uniform and vary inversely with the shrub diameter.     

不同底质和透明度下马来眼子菜的表型可塑性研究

马来眼子菜(Potamogeton malaianusMiq.)是目前太湖沉水植物优势种之一。文章比较分析了3种不同底质和透明度情况下马来眼子菜的生长特征,探索其在太湖不断扩展的原因及其适应性。分别选取粘土质粉砂、粉砂和下蜀黄土底质上生长的马来眼子菜进行观察试验。结果表明,这3种底质在粒度组成、营养盐和分布特征上具有显著差异。在太湖粘土质粉砂和粉砂底质上马来眼子菜的生物量比下蜀黄土底质上的马来眼子菜高,其生物量随着底质营养盐的增加而增加。马来眼子菜与觅光相关的形态指标,如高度、节间距、节数、叶数、叶长及叶面积均随着底质营养盐的增加而显著增加。马来眼子菜的冠层高度与水体的透明度呈显著负相关,表明其对剧烈变化的湖泊光环境具有很强的适应能力。马来眼子菜在不同底质上的形态可塑性是其优先占据湖泊资源成为优势种的重要原因。     

Effects of land use on water quality of two high-altitude lakes and catchments in Yunnan Province,China

SUMMARY

Reclaiming farmland from lakes in China in the 1950s damaged the water quality of many lakes. Tremendous efforts have been made since the late 1990s to restore vegetation around the damaged lakes. This paper examines water quality of Fuxian and Qilu Lakes and land-use characteristics within the two catchments in the high-altitude area of Yunnan Province, China. Landsat TM data acquired in 1989, 1994, 1999 and 2005 were used to extract land use and land cover (LULC) information. Measurements of five water quality indices (WQIs), BOD, COD, pH, TN and TP, for the same period of time were examined. The results showed that the area of residential and forest/shrub increased, whereas that of cropland and barren land decreased from 1989 to 2005 in both catchments. Qilu Lake was much more polluted than Fuxian Lake, and pollution worsened over time for both lakes. The differences in water quality between the two lakes were caused by differences in LULC composition and continued degradation in water quality was caused by intensive farming and urban sprawl. Unless the landscape is converted back to its pre-1950 composition and structure, water quality in both lakes cannot be significantly improved and will continue to threaten sustainable development in the region.     

Hydrodynamics and bed stability of open channel flows with submerged foliaged plants

In this work we investigate experimentally and numerically the flow structure around foliaged plants deployed in a channel with gravels on the bed under submerged and barely submerged conditions. Velocity and Reynolds stress were measured by using a NORTEK Vectrino profiler. Visual observation shows that the initial motion of gravels is easier to be triggered under the condition of flow with barely submerged vegetation. This is confirmed by the measured velocity, Reynolds stress and total kinetic energy (TKE) profiles. The velocity exhibits a speed up in the near-bed region, and the associated Reynolds stress and TKE increase there. A 3D numerical model is then verified against the experiments and used to investigate systematically the effect of degree of submergence of foliaged plants on the channel bed shear stress. The results show that the maximum bed shear stress occurs when the foliage is situated slightly below the water surface, which can enhance channel bed instability.     

A dispersion model for traffic produced turbulence in a two-way traffic scenario

Low wind speed conditions are often associated with poor air quality in urban areas, especially near roadways. Predictions of pollutant concentration under such conditions, i.e. low wind-speeds and near road locations, are, however, complicated by the role of traffic produced turbulence (TPT) on pollutant mixing and dilution. Existing dispersion models consider the effect of TPT on pollutant concentrations near roadways, accounting for parameters such as vehicle intensity, vehicle speeds, etc, but do not explicitly account for the contribution of two-way traffic interaction on the pollutant dispersion parameter. The turbulent kinetic energy (TKE) resulting from a two-way traffic condition will be higher than that with a one-way traffic pattern. Here, we obtain a simple formulation for TKE under a two-way traffic pattern from the balance of production and dissipation of turbulence. Considering the vorticity generated by the two-way traffic and determining the equivalent drag coefficient, an expression for TKE due to the two-way traffic interaction was obtained for three different traffic density regimes: light, intermediate, and heavy. The model predictions are validated by comparison with published data from a field study. An improved parameterization of the TPT considering the two-way traffic interaction effect is seen to significantly improve predictions of near roadway pollutant concentrations.     

几种植物对潜流型湖滨湿地中氮磷的处理效果比较

湖滨带是连接湖泊水域生态系统与陆地生态系统的一个功能过渡区,是湖泊的最后一道保护屏障。在河流的入湖口建造湖滨湿地,可有效净化入湖水中携带的部分有机污染物、营养盐等。湿地净化污染物的机理极为复杂,其中植被起着重要的作用。选择适当的植物是构建湿地系统的关键。以云南省抚仙湖北岸的湖滨湿地——马料河潜流型湿地为研究对象,探讨了美人蕉(Cannaindica),菖蒲(AcoruscalamusLinn),伞竹(Cyperusalternifolius),香蒲(Typhaelatifolia)四种不同植被系统的潜流型湖滨湿地去除污水中营养盐的效果,为抚仙湖湖滨带人工湿地的合理优化选择植物提供一定的依据。结果表明,四种植被系统的潜流湿地对总氮和总磷均有一定的去除效果。但是由于湿地长期运行而疏于管理,植物没有定期收割翻新,四种植物的氮磷去除率差别不太大,去除效果均不太佳。其中,美人蕉对总氮和总磷的去除效果较好,其去除率分别是30.4%,19.13%,其次为香蒲和伞竹,效果较差的是菖蒲。并且实验表明,暴雨对四种植物湿地氮的去除效果有很大的影响。连续几天的暴雨情况下,四种植物系统氮的去除率都很低。     

Exchange flow between open water and floating vegetation

This study describes the exchange flow between a region with open water and a region with a partial-depth porous obstruction, which represents the thermally-driven exchange that occurs between open water and floating vegetation. The partial-depth porous obstruction represents the root layer, which does not penetrate to the bed. Initially, a vertical wall separates the two regions, with fluid of higher density in the obstructed region and fluid of lower density in the open region. This density difference represents the influence of differential solar heating due to shading by the vegetation. For a range of root density and root depths, the velocity distribution is measured in the lab using PIV. When the vertical wall is removed, the less dense water flows into the obstructed region at the surface. This surface flow bifurcates into two layers, one flowing directly through the root layer and one flowing beneath the root layer. A flow directed out of the vegetated region occurs at the bed. A model is developed that predicts the flow rates within each layer based on energy considerations. The experiments and model together suggest that at time- and length-scales relevant to the field, the flow structure for any root layer porosity approaches that of a fully blocked layer, for which the exchange flow occurs only beneath the root layer.     

Simulation of exchange flow between open water and floating vegetation using a modified RNG <Emphasis Type="Italic">k</Emphasis>-<Emphasis Type="Italic">ε</Emphasis> turbulence model

The paper focuses on the numerical simulation of the exchange flow between open water and floating vegetation, which plays an important role in maintaining the ecological balance by transporting nutrient matter. The simulation was conducted using a new solver developed upon OpenFOAM. A modified RNG k-ε turbulence model, which is expected to model both the high- and low-Reynolds number flows correctly, was used to determine the eddy viscosity. Several particular terms were added into the momentum equations and turbulence model equations to model the effects of vegetation and buoyancy. Among these terms, the term for the effect of vegetation in the ε-equation was re-modelled. The model was validated by properly predicting the profiles of mean velocity and turbulent kinetic energy for flows through suspended canopies. The density flow between open and vegetated water was simulated with the same conditions as those of the experiment conducted by Zhang and Nepf. The predicted results agreed well with the experimental data and provided more detailed information of such exchange flow. The convection between the root layer and the layer beneath the roots, which was not observed in the experiment, was observed in the numerical simulation.     

罗非鱼-附着藻-沉水植物相互关系研究进展

在浅水富营养化湖泊中,沉水植被是决定湖泊清水态或混水态的关键因子。而附着藻对沉水植物强烈的遮阴作用以及对碳源、营养盐等资源强烈的竞争,成为限制沉水植物群落生长和发展的关键因子。罗非鱼作为一种杂食性鱼类,具有牧食附着藻的能力,其下行效应（top-down effect）可以在一定程度上减轻附着藻对沉水植物生长的这种不利影响。因此,作为我国南方水体中的优势种类,适当种群密度的罗非鱼在富营养化浅水湖泊生态修复过程中是可加以利用,并在一定程度上抑制了附着藻的生长和发展,有利于浅水湖泊的生态修复和管理。同时,罗非鱼也具有通过摄食、排泄等活动加速水体氮、磷营养盐再生,牧食浮游动物、沉水植物等不利的一面。因此,在综合考虑多种因素条件下,需要对罗非鱼-附着藻-沉水植物三者之间的相互关系进行深入研究,探讨生态系统对罗非鱼的响应,这对我国南方浅水富营养化湖泊的生态恢复与管理,尤其是沉水植被的重建与保护具有重要的理论意义和实践价值。     

草、藻型湖泊水体生态及理化特性的实验对比

2006年9月,根据营养水平和种植水草的差异设计了6个浅水湖泊模拟系统,实验用水草为菹草(Potamogeton crispusLinn.)和马来眼子菜(Potamogeton malaianus-Miq).在15个月实验期间,通过多次监测各系统的景观外貌和水质,对草、藻型湖泊生态及理化特性的差异进行研究,得出以下结论:(1)草、藻型系统分别对应清水和浊水2种状态,景观外貌差异很大.(2)水草可使湖泊系统维持在清水状态,在一定条件下,甚至可使富营养化湖泊维持在清水状态;但是水草腐烂分解等也可使水质迅速恶化,甚至引起湖泊草、藻状态的转变;关键在于,对于不断变化的环境条件,系统内水草能否健康生长.(3)由于营养和生产力水平低,贫营养系统的水质指标随时间变化较小,草、藻型系统间的差异不明显,DO变化范围分别为8.1～14.4 mg·L~(-1)、7.5～11.6 mg·L~(-1),pH 8.71～9.89、8.25～9.22,TP 0.006～0.012 mg·L~(-1)、0.006～0.053 mg·L~(-1),TN 0.11～0.71 mg·L~(-1)、0.10～0.83 mg·L~(-1),NH_4~+-N 0.01～0.17 mg·L~(-1)、0.01～0.26 mg·L~(-1),PO_4~(3-)-P 0.002～0.012 mg·L~(-1)、0.000～0.008mg·L~(-1).(4)由于水草和藻类的大量生长等,中营养与富营养系统湖水的DO、pH、水温和NH_4~+-N的日变化明显,日变化曲线呈“⌒”形,且具有季节性变化规律;由于水草向底泥中输氧气等原因,与藻型湖泊相比,草型湖泊水中TP、TN和NH_4~+-N的浓度较低,PO_4~(3-)-P浓度较高,草、藻型系统的TP均值分别为0.16、0.51 mg·L~(-1),TN 1.30、8.32 mg·L~(-1),NH_4~+-N 0.19、0.43mg·L~(-1),PO_4~(3-)-P 0.07、0.01 mg·L~(-1).     

Turbulent kinetic energy budgets in a model canopy: comparisons between LES and wind-tunnel experiments

A comparative study of turbulence in a wind-tunnel model canopy is performed, using Large eddy simulation (LES) and experimental data from PIV and hot-wire anemometry measurements. The model canopy is composed of thin cylindrical stalks. In the LES, these are represented using a plant-scale approach, while the scale-dependent Lagrangian dynamic model is used as subgrid-scale model. LES predictions of turbulence statistics and energy spectra are found to be in good agreement with the experimental data. Turbulent kinetic energy (TKE) budgets from the LES simulation are analyzed to provide more information absent in the measurements. Results confirm that sloshing motions at the low levels of the canopy are mainly driven by pressure fluctuations. A difference between the energy flux obtained from the energy spectrum and the SGS dissipation rate is observed, consistent with a spectral bypass mechanism in which the real spectral flux due to cascade is smaller than that implied by the energy-spectrum level, due to direct drain by the canopy.     

中国湖泊水生植被退化现状与对策

以中国东部湖区与云贵湖区的若干湖泊为例,论述了中国湖泊水生植被退化的现状与问题。这些问题主要是水生植被分布面积缩小,水生植物分布深度下降,总生物量下降,物种多样性下降,耐污性较强的植物发展迅速。区别于国内外其它研究主要是从水生植物自身的生理学角度阐述水生植被退化的原因,本文是从人口密度、人口压力、经济发展、城市化进程等社会因素方面深层次分析了中国水生植被退化的因素,并且提出了可能的恢复对策与恢复过程中可能存在的问题。认为水生植被恢复,在截污的基础上,营造水生植被适宜生长的小生境至关重要。由于水生植被恢复与自然演替是一个同时存在的动态过程,以及水生植被恢复指标的不可确定等,所以水生植被生态恢复无论在基础理论研究方面还是生态工程实施技术上都有许多值得进一步探讨和研究之处。     

Predictive Model of the Effects on Lake Metabolism of Decreased Airborne Litterfall through Riparian Deforestation

The importance of airborne allochthonous litter to the carbon and nutrient budgets of lakes has been seldom studied. We complied data on the input of terrestrial litter to develop a simple and speculative model to predict the potential consequences of riparian deforestation on one aspect of lake metabolism, specifically the balance between phytoplankton production and plankton respiration. During the autumn of 1992, 56 litter traps were deployed around the littoral zones of four oligotrophic lakes in a densely forested region of northwestern Ontario, Canada. The airborne litter input was estimated to be 32 g dry weight per meter of forested shoreline per year. Allochthonous litter input per unit offshore distance was related to the size of riparian trees, their proximity to the shoreline, and the elevation of their canopy. Combining our data with those from other studies suggests that terrestrial litter can contribute up to 15% of the total carbon supply to oligotrophic lakes and up to 10% of the total phosphorus supply to lakes with a large surface area relative to that of their drainage basin. These results were incorporated into a simple model that predicts that removal of shoreline trees could increase the ratio of plankton production to respiration in oligotrophic lakes situated within small drainage basins. Such lakes may therefore shift from allotrophy to increasing autotropy (energy self-sustenance) following riparian deforestation.     

Modeling approach to regime shifts of primary production in shallow coastal ecosystems

Pristine coastal shallow systems are usually dominated by extensive meadows of seagrass species, which are assumed to take advantage of nutrient supply from sediment. An increasing nutrient input is thought to favour phytoplankton, epiphytic microalgae, as well as opportunistic ephemeral macroalgae that coexist with seagrasses. The primary cause of shifts and succession in the macrophyte community is the increase of nutrient load to water; however, temperature plays also an important role. A competition model between rooted seagrass (Zostera marina), macroalgae (Ulva sp.), and phytoplankton has been developed to analyse the succession of primary producer communities in these systems. Successions of dominance states, with different resilience characteristics, are found when modifying the input of nutrients and the seasonal temperature and light intensity forcing.     

Turbulence suppression by suspended sediment within a geophysical flow

Experiments are performed in a mixing box to evaluate the effect of suspended sediment on turbulence generated by an oscillating grid. Quartz-density sand of varying sizes and concentrations is used, and particle image velocimetry is employed to quantify only the fluid phase. Results show that (1) while a relatively large secondary flow field is present in the box, turbulence is a maximum near the grid and it decreases systematically toward the water surface; (2) relatively high concentrations of fine sediment can markedly alter this secondary flow field and significantly decrease both the time-mean and turbulent kinetic energy within the flow, yet these same sediment concentrations have little effect on the integral time and length scales derived for each velocity component; and (3) the overall turbulence suppression observed can be related to the transfer of energy from the fluid to the sediment and the maintenance of a suspended sediment load rather than commonly employed turbulence modulation criteria. These experimental data demonstrate unequivocally that the presence of a suspended sediment load can significantly reduce overall turbulent kinetic energy, and these results should be applicable to a range of sediment-laden geophysical flows.     

Characteristics of flow structure of free-surface flow in a partly obstructed open channel with vegetation patch

Free-surface flows over patchy vegetation are common in aquatic environments. In this study, the hydrodynamics of free-surface flow in a rectangular channel with a bed of rigid vegetation-like cylinders occupying half of the channel bed was investigated and interpreted by means of laboratory experiments and numerical simulations. The channel configurations have low width-to-depth aspect ratio (1.235 and 2.153). Experimental results show that the adjustment length for the flow to be fully developed through the vegetation patch in the present study is shorter than observed for large-aspect-ratio channels in other studies. Outside the lateral edge of the vegetation patch, negative velocity gradient (\(\partial \overline{u}/\partial z < 0\)) and a local velocity maximum are observed in the vertical profile of the longitudinal velocity in the near-bed region, corresponding to the negative Reynolds stress (\(- \overline{{u^{\prime}w^{\prime}}} < 0\)) at the same location. Assuming coherent vortices to be the dominant factor influencing the mean flow field, an improved Spalart–Allmaras turbulence model is developed. The model improvement is based on an enhanced turbulence length scale accounting for coherent vortices due to the effect of the porous vegetation canopy and channel bed. This particular flow characteristic is more profound in the case of high vegetation density due to the stronger momentum exchange of horizontal coherent vortices. Numerical simulations confirmed the local maximum velocity and negative gradient in the velocity profile due to the presence of vegetation and bed friction. This in turn supports the physical interpretation of the flow processes in the partly obstructed channel with vegetation patch. In addition, the vertical profile of the longitudinal velocity can also be explained by the vertical behavior of the horizontal coherent vortices based on a theoretical argument.     

乌东德库区消落带土壤水分分布特征

消落带土壤水分是制约植物生长发育的关键因素,对库区消落带生态修复意义重大.本文系统分析了乌东德库区消落带修复试验区土壤含水量分布,探索了土壤水分空间变化和影响因素.结果表明:研究区5月水分含量低于5％,不同试验分区和高程区域土壤水分分布受到微地形、土壤特征、植被类型和植被覆盖度等影响.土壤水分随深度变化的主控因素包括地...