Hydrodynamics and contaminant transport on a degraded bed at a 90-degree channel confluence

Channel confluences at which two channels merge have an important effect on momentum exchange and contaminant diffusion in both natural rivers and artificial canals. In this study, a three-dimensional numerical model, which is based on the Reynolds Averaged Navier–Stokes equations and Reynolds Stress Turbulence model, is applied to simulate and compare flow patterns and contaminant transport processes for different bed morphologies. The results clearly show that the distribution of contaminant concentrations is mainly controlled by the shear layer and two counter-rotating helical cells, which in turn are affected by the discharge ratio and the bed morphology. As the discharge ratio increases, the shear flow moves to the outer bank and the counter-clockwise tributary helical cell caused by flow deflection is enlarged, leading the mixing happens near the outer bank and the mixing layer distorted. The bed morphology can induce shrinkage of the separation zone and increase of the clockwise main channel helical cell, which is initiated by the interaction between the tributary helical cell and the main channel flow and strengthened by the deep scour hole. The bed morphology can also affect the distortion direction of the mixing layer. Both a large discharge ratio and the bed morphology could lead to an increase in mixing intensity.     

A model of vertical distribution of suspended matter in an open channel flow

A quasi-stationary model of vertical distribution of concentration of suspended particular matter in the bottom layer of 1D open channel with a sloped bottom and varying free surface slope is discussed. The model proceeds from the balance between the turbulent diffusion and settling with the buoyancy flux effects on the medium turbulence neglected. The model outcome is formulated in the form of an analytic formula for the vertical distribution of concentration. It is shown that the derived formula embraces two basic types of vertical distribution of concentration, one with a monotonic decrease of concentration gradient and the other with a gradient maximum (lutocline) located at some distance from the bottom. The first distribution type realizes for a relatively large settling velocity or low intensity of turbulence and the second type for a small settling velocity or high intensity of turbulence. The skill of the model to mimic realistic situations is demonstrated on data measured in the Jiaojiang Estuary (China).     

Curvature-induced secondary flow in submarine channels

The curvature-driven secondary flow in sinuous submarine channels has been a subject of considerable interest and controversy. Here, results from numerical model studies involving saline flow in laboratory-scale channels are presented. A 3D finite volume model of density and turbidity currents is used and simulations are run with different inflow discharges and channel-axis slopes. The simulation results show strong influence of bend wave length, channel gradient, confinement and cross sectional shape on the structure of secondary flow in submarine channels. Major findings are: (i) reversal of secondary flow in submarine channels is strongly associated with a tight bend characterized by a smaller wave length to width ratio or larger wave number, (ii) for the same inflow condition and planform characteristics, a trapezoidal channel cross section is more favorable to secondary flow reversal than a rectangular cross section, (iii) lateral convection resulting from the interaction between in-channel and overbank flows leads to the reversal of secondary flow in an unconfined channel at a lower channel slope than in a confined channel with the same dimensions, (iv) flow discharge has only nominal effect on the secondary flow in submarine channels.     

Influence of the secondary motions on pollutant mixing in a meandering open channel flow

This paper presents large eddy simulation of turbulent flow in a meandering open channel with smooth wall and rectangular cross-section. The Reynolds number based on the channel height is 40,000 and the aspect ratio of the cross-section is 4.48. The depth-averaged mean stream-wise velocity agree well to experimental measurements. In this specific case, two interacting cells are formed that swap from one bend to the other. Transport and mixing of a pollutant is analysed using three different positions of release, e.g. on the inner bank, on the outer bank and on the centre of the cross section. The obtained depth-average mean concentration profiles are reasonably consistent with available experimental data. The role of the secondary motions in the mixing processes is the main focus of the discussion. It is found that the mixing when the scalar is released on the centre of the cross-section is stronger and faster than the mixing of the scalar released on the sides. When the position of release is close to a bank side, the mixing is weaker and a clear concentration of scalar close to the corresponding side-wall can be observed in both cases.     

应用GIS模拟城市大气污染物浓度分布

从点源大气污染扩散模型出发,借助GIS的栅格与数据库运算功能,结合自已设计的推算流程,对兰州市大气污染物SO2的空间分布情况进行了模拟。经检验,该方法能够较准确的反映出兰州大气污染物的空间分布特征,较应用实测值进行插值的方法简便,具有较好的应用前景。     

Modeling the buoyant flow of heated water discharged from surface and submerged side outfalls in shallow and deep water with a cross flow

In this study, a three-dimensional model was used to numerically study the buoyant flow, along with its mixing characteristics, of heated water discharged from the surface and submerged side outfalls in shallow and deep water with a cross flow. Hydraulic experimental data were used to evaluate the applicability of the model. The simulation results agree well with the experimental results, particularly, the jet trajectories, the dimensions of the recirculating zone, and the distribution of the dimensionless excess temperature. The level of accuracy of the simulation results of the present study is nearly identical to that of the results conducted by McGuirk and Rodi (1978). If the heated water is discharged into shallow water where the momentum flux ratio and the discharge densimetric Froude number are high, the submerged discharge method is better than the surface discharge method in terms of the scale of the recirculating zone and the minimum dilution. In deep water, where the momentum flux ratio and discharge densimetric Froude number are low, however, the submerged discharge method had few advantages. In shallow water, the discharge jet is deflected by the ambient cross flow, while forcing the ambient flow to bend towards the far bank for the full depth. For a submerged discharge in shallow water, the recirculating zone is the largest in the lowest layer but becomes smaller in the upper layer. As the water depth increases, the ambient flow goes over the jet and diminishes the blocking effect, thereby decreasing the bending of the jet.     

An ecological risk assessment model for a pulsed contaminant emission into a wetland channel flow

As a continuation of the modelling on ecological degradation and hydraulic dispersion of pollutant emission into an idealized two-dimensional free-surface wetland flow (Zeng, L., Chen, G.Q., 2009b. Ecological degradation and hydraulic dispersion of contaminant in wetland. Ecol. Model., doi:10.1016/j.ecolmodel.2009.10.024), an ecological risk assessment model for the typical case of a pulsed contaminant emission into a realistic three-dimensional wetland channel flow is presented in this paper for the fate of cross-sectional mean concentration under environmental dispersion. An environmental dispersion model for the mean concentration is devised as an extension of Taylor’s classical analysis on dispersion in fluid flows. The velocity distribution and the environmental dispersivity in the fully developed steady flow through the wetland is found and illustrated with limiting cases covering various known solutions for the porous media flow between parallel plates, flow in a shallow wetland, sweeping flow in a densely vegetated wetland, and single phase flow in a channel. Obtained by Aris’s method of moments, the environmental dispersivity is shown characterized with multi-scale asymptotic time variations with stem dominated stage, transitional stage, and width-depth-stem dominated stage. Based on the solution for the evolution of contaminant cloud in the wetland channel flow, critical length and duration of the contaminant cloud with concentration beyond given environmental standard level are concretely illustrated for typical pollutant constituents in wastewater emission. Under the same emission intensity and environmental standard, the duration of contaminant cloud in the wetland channel is revealed shorter than that in a free surface wetland, due to the lateral effect.     

Mass-based depth and velocity scales for gravity currents and related flows

Gravity driven flows on inclines can be caused by cold, saline or turbid inflows into water bodies. Another example are cold downslope winds, which are caused by cooling of the atmosphere at the lower boundary. In a well-known contribution, Ellison and Turner (ET) investigated such flows by making use of earlier work on free shear flows by Morton, Taylor and Turner (MTT). Their entrainment relation is compared here with a spread relation based on a diffusion model for jets by Prandtl. This diffusion approach is suitable for forced plumes on an incline, but only when the channel topography is uniform, and the flow remains supercritical. A second aspect considered here is that the structure of ET’s entrainment relation, and their shallow water equations, agrees with the one for open channel flows, but their depth and velocity scales are those for free shear flows, and derived from the velocity field. Conversely, the depth of an open channel flow is the vertical extent of the excess mass of the liquid phase, and the average velocity is the (known) discharge divided by the depth. As an alternative to ET’s parameterization, two sets of flow scales similar to those of open channel flows are outlined for gravity currents in unstratified environments. The common feature of the two sets is that the velocity scale is derived by dividing the buoyancy flux by the excess pressure at the bottom. The difference between them is the way the volume flux is accounted for, which—unlike in open channel flows—generally increases in the streamwise direction. The relations between the three sets of scales are established here for gravity currents by allowing for a constant co-flow in the upper layer. The actual ratios of the three width, velocity, and buoyancy scales are evaluated from available experimental data on gravity currents, and from field data on katabatic winds. A corresponding study for free shear flows is referred to. Finally, a comparison of mass-based scales with a number of other flow scales is carried out for available data on a two-layer flow over an obstacle. Mass-based flow scales can also be used for other types of flows, such as self-aerated flows on spillways, water jets in air, or bubble plumes.     

QUAL2E模型在大沽河干流青岛段水质模拟中的应用

采用QUAL2E模型对大沽河干流青岛段的水质进行了模拟和预测。针对大沽河的具体情况,选用BOD5、COD和氮作为模拟预测指标,用实验模拟方法、模型率定法并参考相关文献确定了BOD耗氧系数k1、BOD复氧系数k2、BOD沉降系数k3、COD耗氧系数和弥散系数等水质参数,并对模拟结果进行了验证,表明预测值和实测值的相关性较好;对BOD5、k1、k2和Q(流量)进行了灵敏度分析,结果表明对大沽河DO浓度影响敏感的参数依次是:Q、k2、BOD、k1,即流量Q是模型最敏感的参数,说明河流的水力学参数对DO影响较大。     

外源性污染对太湖梅梁湾水质影响的定量化

通过太湖的水量水质数学模型,模拟了梅梁湾内主要的入湖污染源———直湖港、武进港和梁溪河排入的污染物质的迁移、转化规律;分析了梅梁湾中梅园、小湾里、闾江口、拖山4个监测点水质浓度受直湖港、武进港和梁溪河排污影响程度的大小;并分别建立了监测点污染物浓度与排污口排污量的响应关系曲线和响应关系表达式。通过这些结论可以方便的定量的计算出在排污口排污量发生变化的情况下监测点水质浓度的变化量,为合理控制外源污染物质的入湖量提供了技术支持。     

A σ-coordinate transport model coupled with rotational boussinesq-type equations

This paper describes a σ-coordinate scalar transport model coupled with a Boussinesq-type hydrodynamic model. The Boussinesq model has the ability to calculate both three-dimensional velocity distributions and the water surface motion. To capture ‘dispersion’ processes in open channel flow, horizontal vorticity effects induced by a bottom shear stress are included in the Boussinesq model. Thus, a reasonable representation of vertical flow structure can be captured in shallow and wavy flow fields. To solve the coupled Boussinesq and scalar transport system, a finite-volume method, based on a Godunov-type scheme with the HLL Riemann solver, is employed. Basic advection and advection–diffusion numerical tests in a non-rectangular domain were carried out and the computed results show good agreement with analytic solutions. With quantitative comparisons of dispersion experiments in an open channel, it is verified that the proposed coupled model is appropriate for both near and far field scalar transport predictions. From numerical simulations in the surf zone, physically reasonable results showing expected vertical variation are obtained.     

Experimental study of recirculating flows generated by lateral shock waves in very large channels

Due to the lack of data on hydraulic-jump dynamics in very large channels, the present paper describes the main characteristics of the velocity field and turbulence in a large rectangular channel with a width of 4 m. Although a hydraulic jump is always treated as a wave that is transversal to the channel wall, in the case of this study it has a trapezoidal front shape, first starting from a point at the sidewalls and then developing downstream in an oblique manner, finally giving rise to a trapezoidal shape. The oblique wave front may be regarded as a lateral shockwave that arises from a perturbation at a certain point of the lateral wall and travels obliquely toward the centreline of the channel. The experimental work was carried out at the Coastal Engineering Laboratory of the Water Engineering and Chemistry Department of the Technical University of Bari (Italy). In addition to the hydraulic jump formation, a large recirculating flow zone starts to develop from the separating point of the lateral shock wave and a separate boundary layer occurs. Intensive measurements of the streamwise and spanwise flow velocity components along one-half width of the channel were taken using a bidimensional Acoustic Doppler Velocimeter (ADV). The water surface elevation was obtained by means of an ultrasonic profiler. Velocity vectors, transversal velocity profiles, turbulence intensities and Reynolds shear stresses were all investigated. The experimental results of the separated boundary layer were compared with numerical predictions and related work presented in literature and showed good agreement. The transversal velocity profiles indicated the presence of adverse pressure gradient zones and the law of the wall appears to govern the region around the separated boundary layer.     

A stochastic Lagrangian micromixing model for the dispersion of reactive scalars in turbulent flows: role of concentration fluctuations and improvements to the conserved scalar theory under non-homogeneous conditions

Several reaction schemes, based on the conserved scalar theory, are implemented within a stochastic Lagrangian micromixing model to simulate the dispersion of reactive scalars in turbulent flows. In particular, the formulation of the reaction-dominated limit (RDL) reaction scheme is here extended to improve the model performance under non-homogeneous conditions (NHRDL scheme). The validation of the stochastic model is obtained by comparison with the available measurements of reactive pollutant concentrations in a grid-generated turbulent flow. This test case describes the dispersion of two atmospheric reactant species (NO and O₃) and their reaction product (NO₂) in an unbounded turbulent flow. Model inter-comparisons are also assessed, by considering the results of state-of-the-art models for pollutant dispersion. The present validation shows that RDL reaction scheme provides a systematic overestimation (relative error of ca. 85% around the centreline) in computing the local reactant consumption/production rate, whereas the NHRDL scheme drastically reduces this gap (relative error lower than 5% around the centreline). In terms of NO₂ production (or reactant consumption), neglecting concentration fluctuations determines overestimations of the product mean of around 100% and a NO₂ local production of one order of magnitude higher than the reference simulation. In terms of standard deviations, the concentration fluctuations of both the passive and reactive scalars are generally of the same order of magnitude or up to 1 or 2 orders of magnitudes higher than the corresponding ensemble mean values, except for the background reactant close to the plume edges. The study highlights the importance of modelling pollutant reactions depending on the instantaneous instead of the mean concentrations of the reactants, thus quantifying the role of the turbulent fluctuations of concentration, in terms of scalar statistics (mean, standard deviation, intensity of fluctuations, skewness and kurtosis of concentration, segregation coefficient, simulated reaction rate). This stochastic particle method represents an efficient numerical technique to solve the convection–diffusion equation for reactive scalars and involves several application fields: micro-scale air quality (urban and street-canyon scales), accidental releases, impact of odours, water quality and fluid flow industrial processes (e.g. combustion).     

漓江上游纳污最适总量及水质变化规律研究

为确保 2 0 0 0 - 2 0 10年桂林漓江上游进入城区水质能保持国家II类水质标准进行纳污总量控制及其水质变化规律研究。采用系数法对 2 0 0 0 - 2 0 10年陆源排污总量进行科学预测 ,按河流功能区划要求 ,首次利用动态数学模型计算河流在不同流量、不同纳污总量条件下保持II类水质的最适纳污值 ,探讨了漓江水质变化规律。为环境管理打下基础     

Two-phase SPH modelling of advective diffusion processes

This paper deals with a two-dimensional numerical model based on the smoothed particle hydrodynamics (SPH) technique for the evaluation of the concentration field of pollutants in water. A SPH model is formulated to solve the fickian diffusion equation applied to pollutants with the same density as the water. A lagrangian SPH formalism of the advective diffusion equation is also developed for pollutant-water, taking into account the effects of molecular diffusion and natural advection induced by differences between the fluid densities. These equations are coupled with the fluid mechanics equations. Attention is paid to the numerical aspects involved in the solution procedure and to the optimization of the model parameters. Environmental engineering problems concerning diffusion and natural advection phenomena occur in the presence of a pollutant in still water. Numerical tests referring to a strip and a bubble of contaminant in a water tank with different initial concentration laws have been carried out. The results obtained by the proposed SPH models are compared with other available SPH formulations, showing an overall better agreement with standard analytical solutions in terms of spatial evolution of the concentration values. Capabilities and limits of the proposed SPH models to simulate advective diffusion phenomena for a wide range of density ratios are discussed.     

Deterministic particle tracking simulation of pollutant discharges in rivers and estuaries

A two-dimensional deterministic particle tracking model, in which the anisotropic-dispersive process is described by a particle strength exchange scheme, is established for the simulation of pollutant transport in vertically well-mixed rivers and estuaries. By simulating two benchmark problems with analytic solutions, the PSE scheme is shown to be accurate even if the anisotropic ratio of dispersion coefficients is very high. Further simulations of two specific problems concerning the optimal effluent discharge location and procedure are presented. The major conclusion is that in a tidal estuary with a relatively large fresh-water flow, setting the discharge position at the mixing center and making the discharge rate proportional to flow speed may minimize the peaks of concentration.     

机动车污染模拟和地理信息系统支持

机动车污染是城市大气污染物的重要来源，对街区机动车污染排放量的预测和排放模拟对了解城市街区环境污染状况有重要价值。文章中将ＣＡＲ模式的模拟计算和ＧＩＳ集成在一起，通过ＶｉｓｕａｌＢａｓｉｃ进行了一个样例开发，其中ＧＩＳ功能由ＭａｐＸ控件实现。     

Improved formulae of velocity distributions along the vertical and transverse directions in natural rivers with the sidewall effect

The distribution of flow velocity is a basis for the research into the transport of flow and sediment in natural rivers. Characteristics of velocity distribution in narrow-deep natural rivers are different from those in wide-shallow open channels, and the effect of sidewalls on the distribution of flow velocity is considerable, which leads to a large transverse gradient of the depth-averaged velocity, with the maximum velocity occurring below the water surface. Based on the Reynolds equation of turbulence flow and the analysis of the features of velocity distribution in natural rivers, improved formulae with two empirical parameters α and β have been proposed for the velocity distributions along the vertical and transverse directions, with the effect of sidewalls being considered, through solving the definite solution problem by the method of variable separation. The proposed formulae were validated fully through comparisons between the calculated and measured velocity profiles and depth-averaged velocities at several sections in the Yangtze and Baitarani Rivers, with close agreement between them being obtained. The formula of velocity distribution along the transverse direction in natural rivers with the sidewall effect was also compared with previous studies, and the calculation accuracy of this formula at a section with a narrow-deep geometry was higher than the accuracy of the previous equations. It is confirmed that the proposed formulae can reproduce well the distribution characteristics of flow velocity along the vertical and transverse directions in narrow-deep natural rivers, with a more wide application in practice.     

Effects of secondary current and stratification on suspension concentration in an open channel flow

In this work, a mathematical model on concentration distribution is developed for a steady, uniform open channel turbulent flow laden with sediments by incorporating the effect of secondary current through velocity distribution together with the stratification effect due to presence of sediments. The effect of particle-particle interaction at reference level and the effect of incipient motion probability, non-ceasing probability and pick-up probability of the sediment particles at reference concentration are taken into account. The proposed model is compared with the Rouse equation as well as verified with existing experimental data. Good agreement between computed value and experimental data indicates that secondary current influences the suspension of particles significantly. The direction and magnitude (strength) of secondary current lead to different patterns of concentration distribution and theoretical analysis shows that type II profile (where maximum concentration appears at significant height above channel bed surface) always corresponds to upward direction and greater magnitude of secondary current.