SUPPRESSION OF THE SECONDARY FLOWS IN A BEND PIPE USING NAVIER-STOKES SOLVER AND EVOLUTIONARY ALGORITHMS

Hydrodynamic optimization design of the bend pipe from pump using the Navier-Stokes solver and evolutionary algorithms was conducted. The minimization of the total pressure loss of the bend pipe was chosen as the design object in order to obtain the uniform exit flows through suppressing the secondary flows. The 3-D Navier-Stokes solver was applied to evaluate the hydrodynamic performance of the bendpipe flows. A 7th-order Bezier curve was used to parameterize the meridional section and elliptic representation was adopted to represent the cross-section profiles of the bend pipe. Evolutionary algorithms were applied in optimization. The obtained results show that the designed bend pipe shape has much more uniform exit flows compared with the initial one and much weaker secondary flows, and that the evolutionary algorithms and CFD technique are the powerful optimization tools for the fluid machinery design.     

A SPLIT-CHARACTERISTIC FINITE ELEMENT MODEL FOR 1-D UNSTEADY FLOWS

An efficient and accurate solution algorithm was proposed for 1-D unsteady flow problems widely existing in hydraulic engineering. Based on the split-characteristic finite element method, the numerical model with the Saint-Venant equations of 1-D unsteady flows was established. The assembled finite element equations were solved with the tri-diagonal matrix algorithm. In the semi-implicit and explicit scheme, the critical time step of the method was dependent on the space step and flow velocity, not on the wave celerity. The method was used to eliminate the restriction due to the wave celerity for the computational analysis of unsteady open-channel flows. The model was verified by the experimental data and theoretical solution and also applied to the simulation of the flow in practical river networks. It shows that the numerical method has high efficiency and accuracy and can be used to simulate 1-D steady flows, and unsteady flows with shock waves or flood waves. Compared with other numerical methods, the algorithm of this method is simpler with higher accuracy, less dissipation, higher computation efficiency and less computer storage.     

NUMERICAL SIMULATION OF TWO-DIMENSIONAL OVERTOPPING AGAINST SEAWALLS ARMORED WITH ARTIFICIAL UNITS IN REGULAR WAVES

A new mathematical model for the overtopping against seawalls armored with artificial units in regular waves was established. The 2-D numerical wave flume, based on the Reynolds Averaged Navier-Stokes (RANS) equations and the standard k-ε turbulence model, was developed to simulate the turbulent flows with the free surface, in which the Volume Of Fluid (VOF) method was used to handle the large deformation of the free surface and the relaxation approach of combined wave generation and absorbing was implemented. In order to consider the effects of energy dissipation due to the armors on a slope seawall, a porous media model was proposed and implemented in the numerical wave flume. A series of physical model experiments were carried out in the same condition of the numerical simulation to determine the drag coefficient in the porous media model in terms of the overtopping discharge. Compared the computational value of overtopping over the seawall with the experimental data, the values of the effective drag coefficient was calibrated for the layers of blocks at different locations along the seawalls.     

AN OPERATOR-SPLITTING ALGORITHM FOR ADVECTION-DIFFUSION-REACTION EQUATION

An operator-splitting algorithm for three-dimensional advection-diffusion-reactionequation is presented.The method of characteristics is adopted for the pure advection operator,the explicit difference scheme is used for diffusion,and a prediction-correction scheme is em-ployed for reaction.The condition for stability of the algorithm is analysed.Severall inear andnonlinear examples are illustrated to test the convergence and accuracy of the numerical proce-dure,and satisfactory agreements between computed and analytical solutions are achieved.Dueto its simplicity,stability,and validity for both one-and two-dimensional problems,the success-ful algorithm can be used to numerical simulations of viscous fluid flows,the transport of pollu-tants and sedimentations in reservoirs,lakes,rivers,estuaries and other environments,cooling-problems in heat or nuclear power plants,etc.     

The mechanical energy equation for total flow in open channels

The mechanical energy equation is a fundamental equation of a 1-D mathematical model in Hydraulics and Engineering Fluid Mechanics. This equation for the total flow used to be deduced by extending the Bernoulli＇s equation for the ideal fluid in the streamline to a stream tube, and then revised by considering the viscous effect and integrated on the cross section. This derivation is not rigorous and the effect of turbulence is not considered. In this paper, the energy equation for the total flow is derived by using the Navier-Stokes equations in Fluid Mechanics, the results are as follows：（1） A new energy equation for steady channel flows of incompressible homogeneous liquid is obtained, which includes the variation of the turbulent kinetic energy along the channel, the formula for the mechanical energy loss of the total flow can be determined directly in the deduction process.（2） The theoretical solution of the velocity field for laminar flows in a rectangular open channel is obtained and the mechanical energy loss in the energy equation is calculated. The variations of the coefficient of the mechanical energy loss against the Reynolds number and the width-depth ratio are obtained.（3） The turbulent flow in a rectangular open channel is simulated using 3-D Reynolds averaged equations closed by the Reynolds stress model（RSM）, and the variations of the coefficient of the mechanical energy loss against the Reynolds number and the width-depth ratio are discussed.     

NUMERICAL SIMULATION OF 3-D TURBULENT FLOWS OF PLUNGE POOL AND ENERGY DISSIPATION ANALYSIS

A 3-D numerical model for simulating the complicated turbulent flows was developed and the code was made. A numerical calculation of the plunge pool of Laxiwa project in China was carried on. Those 3-D distributions of velocity,pressure on bottom wall,turbulence kinetic energy and turbulence energy dissipation rate are revealed in detail and the detailed flow patterns of plunge pool were shown. By studying on the characteristics of turbulent diffusion and energy dissipation,the calculated results show that the major turbulence and energy dissipation are taken place near the axis of water jet. The calculated results also indicated that the calculated maximum impact pressures on the bottom wall of the plunge pool have a good agreement with those obtained by physical hydraulic model test.     

NUMERICAL EXPERIMENTS ABOUT GROWTH OF THE ENERGY OF DISTURBANCE IN A BAROCLINIC FLOW

This article is concerned with the growth of energy of disturbances in a baroclinic flow within a finite time period. The implicit difference scheme was applied to the linearized vorticity equation, and the disturbance energy was computed for three kinds of vertical shears. It turns out that all the disturbance energy rapidly increases initially, and during the succeeding period there are several stages of growth and decay of energy of disturbances, and from a certain time on, all the disturbance energy begins to decrease.     

SIMULATION OF SEDIMENT-LADEN FLOW BY DEPTH-AVERAGED MODEL BASED ON UNSTRUCTURED COLLOCATED GRID

The 2-D depth-averaged mathematical model for sediment-laden flows has been widely used in river control and other related engineering problems, and now it is usually solved on structured grids. Since the natural river is usually very complicated in plane boundary, and unstructured grids are more attractive in solving the problems with complicated domains, the following questions about solving 2-D depth-averaged model were discussed in this article: (1) a modified Bowyer algorithm was suggested to generate unstructured grids for natural rivers, (2) the Finite Volume Method (FVM) is employed to discretize the governing equations of the 2-D depth-averaged model and an implicit scheme was suggested with unstructured collocated grids, (3) the observed hydrological data of the Chenglingji Reach in the Yangtze River are used for verification of the presented method. It seems that the suggested numerical scheme works very well, and the simulation results of both hydraulic characteristics and river bed deformation are in good agreement with the observed ones.     

NUMERICAL SIMULATION OF TWO-DIMENSIONAL DAM-BREAK FLOWS IN CURVED CHANNELS

Two-dimensional transient dam-break flows in a river with bends were theoretically studied. The river was modeled as a curved channel with a constant width and a flat bottom. The water was assumed to be an incompressible and homogeneous fluid. A channel-fitted orthogonal curvilinear coordinate system was established and the corresponding two-dimensional shallow-water equations were derived for this system. The governing equations with well-posed initial and boundary conditions were numerically solved in a rectangular domain by use of the Godunov-type finite-difference scheme, which can capture the hydraulic jump of dam-break flows. The comparison between the obtained numerical results and the experimental data of Miller and Chaudry in a semicircle channel shows the validity of the present numerical scheme. The mathematical model and the numerical method were applied to the dam-break flows in channels with various curvatures. Based on the numerical results, the influence of river curvatures on the dam-break flows was analyzed in details.     

Numerical Model for Flow Motion with Vegetation

A set of governing equations for turbulent flows in vegetated area were derived with the assumption that vegetation is of straight and rigid cylinder. The effect of vegetation on flow motion was represented by additional inertial and drag forces. The new model was validated by available experimental data for open channel flows passing through vegetated areas with different vegetation size, density and distribution. Numerical results are in good agreement with the experimental data. Finally, the flow around a supposed isolated vegetated pile was simulated and the effects of vegetation density on the wake flow were discussed. It is found that the presence of vegetation, even at a very low density, has the pronounced influence on the dissipation of flow energy, both inside the vegetation domain and outside it in the wake flow region.     

Call for Papers

正As the eldest and most distributed professional journal of water industry,Water Wastewater Engineering(W WE)gained very wide acceptance from the colleagues working in this field in this country since start publication in 1964.Now it is administrated by the Ministry of Construction of the People’s Republic of China and sponsored by the China Civil Engineering Society etc.     

库水波动下库岸边坡地下水浸润线解析解及误差分析

解析法是库岸边坡地下水浸润线计算中便于实际应用的方法,但该方法须基于若干假定并对潜水运动基本方程线性化后才能求解。针对各假定和线性化过程建立不同的地下水渗流数学模型,用解析法和有限元法解答上述数学模型,分析各误差大小及其规律。结果表明：库岸垂直处理带来的误差要小于方程线性化处理和不考虑非饱和渗流带来的误差,而且它们都随渗透系数变大而减小;在三峡库区库水调度情况下,假定库水位等速变化带来的误差对大多岸坡而言可以忽略;计算库岸边坡地下水浸润线时,解析法只适用于水位变化幅度相比含水层厚度较小,且几何边界规则、岩土结构简单、岩土体渗透性较好时的情况。     

Changes in Summer Irrigated Crop Area and Water Use in Southeastern Turkey from 1993 to 2002: Implications for Current and Future Water Resources

Changes in summer irrigated cropland acreage and related water use are estimated from satellite remote sensing and ancillary data in semi-arid Southeastern Turkey where traditionally dry agricultural lands are being rapidly transformed into irrigated fields with the help of water from the Euphrates-Tigris Rivers. An image classification methodology based on thresholding of Landsat NDVI images from the peak summer period reveals that the total area of summer irrigated crops has increased three-fold (from 35,000 ha to over 100,000) in the Harran Plain between 1993 and 2002. Coupled analysis of annual irrigated crop area from remote sensing and potential evapotranspiration based estimates of irrigation water requirements for cotton indicate a corresponding increase in agricultural water use from about 370 million cubic meters to over one billion cubic meters, a volume in accordance with the state estimates. These estimates have important implications for understanding the rapid changes in current agricultural withdrawals in Southeastern Turkey and form a quantitative basis for exploring the changes in future water demands in the region. For example, expansion of irrigated lands have led to a steady decrease in potential evaporation due to increased roughness and decreased humidity deficit in the Harran Plain. Assuming that the changes in future evaporation conditions will be of similar nature, water use for irrigation is expected to decrease over 40 percent in future irrigation sites. Incorporating this decrease in overall planning of the irrigation projects currently under construction should lead to improved management, and by extension, sustainability of water resources in the region.     

浅谈优化Vray渲染在制作三维场景的应用研究

用Vray三维渲染引擎来渲染三维场景,其质量和速度是由一些影响因素和参数设置决定的,不同的参数和因素,可以得出不同的三维渲染效果.因此,分析、优化三维渲染质量和速度的影响因素,正确、合理的选择参数对提高Vray最后渲染阶段的质量和速度非常重要.     

固结、渗流和沉积控制方程研究

本文从两相流体动力学的守恒方程出发，考虑固相、液相的相对运动，推导了可以考虑场的不均匀性和流体压缩性的统一形式的控制方程，并对方程所基于的物理概念和相互联系进行了必要的分析。通过对方程适当形式的表达，可导出用于大、小变形固结、渗流和沉积的简化方程，显示了在不同领域的三个理论存在一定的统一性。方程推导过程并未对本构关系加以限制，适用范围较广。依据实际问题的要求，作了不同形式的简化。     

四川省及重庆市中型水电站效益分析

该文分析了四川省及重庆市几座中型水电站的调查材料表明：四川省及重庆市中型水电站不仅在川渝主网和地方电网中起着很大的作用，并具有巨大的社会效益，同时还有综合利用效益，说明修建中型水电站的优越性．     

Site modelling methods for detecting hydrologic alteration of flood frequency and flood duration in the floodplain below the carlyle dam,lower Kaskaskia River,Illinois, USA

Reports concerning the influence of dams on river hydrology vary among researchers, interest groups and government agencies. These often contradicting statements may occur because changes in hydrology caused by dams are distinct for each dam and river watershed. The objective of this research was to use site specific techniques to determine if the 1967 installation of the Carlyle Dam, lower Kaskaskia River, Illinois, altered flood frequency and duration within the forested floodplain located below the dam. Results indicated a decrease in flood duration and frequency, and a decrease in annual flood frequency variation at a site 6.4 km below the dam. Pre‐dam versus post‐dam differences in flood frequency and duration at the site 32.2 km below the dam were related to climate rather than dam effects. Although dam impacts are a concern, this research shows that distance downstream from the dam and downstream tributary and watershed characteristics should be considered before assuming that the dam has changed hydrologic parameters for portions of rivers. This research also indicates that areas of the lower Kaskaskia River may still maintain hydrologic ecological integrity, and could be targeted for restoration and adaptive management purposes. Hydrologic modelling combined with river gage and on‐site well measurement techniques presented in this study could provide detailed flood frequency and duration information for land use, sociological and geomorphological questions in focus areas within river floodplains. Copyright © 2008 John Wiley & Sons, Ltd.     

Statistical tools for thermal regime characterization at segment river scale: Case study of the Ste‐Marguerite River

Suitable thermal fish habitats are constrained by both maximum and minimum temperature tolerances. A multivariate and geostatistical approach was developed to estimate stream thermal characteristics at the river segment scale. Data from 22 temperature‐monitoring stations during summer 2007 were used to estimate monthly maximum temperature as well as thermal characteristics such as the number of events, the cumulative degree–days and the associated duration over specific temperature thresholds of 19 and 21°C. The probability of exceeding these temperature thresholds has also been interpolated. The methodology relies on the construction of a multivariate space using physiographic and hydrological characteristics of gauging stations as inputs in a canonical correlation analysis (CCA). A geostatistical interpolation technique, ordinary kriging, was subsequently used to perform interpolation in the physiographical space constructed using CCA. Results from this study were obtained for thermal characteristics estimated into two different interpolation spaces: (1) a 7 metrics space, and (2) an 8 metrics space. Cross‐validation technique has been performed and satisfactory results were obtained. Kriging thermal characteristics (magnitude and duration) into the 7 metric space for a 19°C threshold exceedance leads to best results with Relative Root Mean Square Error (RRMSE) ranging between 9.66 and 15.08%. The study shows that kriging in a multivariate space is a promising tool for water resources managers, especially in cases where risk mapping for lethal or sub‐lethal temperature thresholds may be required for a specific fish species. Copyright © 2010 John Wiley & Sons, Ltd.     

掺MgO混凝土的宏细观试验研究

利用MgO膨胀剂的延迟微膨胀效应能显著提高大体积混凝土的力学性能,产生的预压应力还能有效提高大体积混凝土抵抗温度开裂的能力,但目前对其微观机理研究不多。采用宏观与微观相结合的方式,研究了不同活性指数、不同掺量的MgO混凝土的力学性能,利用SEM/EDS微观分析手段揭示了MgO对混凝土性能影响的作用机理。试验结果表明：MgO取代部分水泥掺入混凝土中,混凝土的力学性能降低;掺量为4%~5%时,活性指数为100s的MgO混凝土的膨胀率大于活性指数为50 s的MgO混凝土;MgO混凝土的膨胀特性与Ca,Si,Mg,Al等元素的分布情况有关,掺入的MgO将改变界面区Ca,Si,Mg元素的富集特性,Ca元素和Mg元素更易于在界面区富集,导致混凝土宏观力学性能的降低。     

基于MIDAS软件探讨施工工序对深基坑稳定性的影响

基于武汉地铁名都站深基坑工程与水文地质勘察资料,建立了三维工程地质仿真计算模型。依据名都站深基坑开挖支护方案,利用有限元软件MIDAS的摩尔-库仑本构模型,对武汉地铁名都站深基坑施工过程中,不同施工工序下每步开挖之后,基坑的变形情况做仿真模拟计算,并对比分析其对基坑稳定性的影响。结果表明：基坑围护结构的变形情况与施工工序有很大的关系,施工工序越合理,上部土体变形越小,反之越大。因此,在基坑工程的施工过程中,要注重施工步骤的划分与合理安排,从而保证基坑工程的安全。研究成果为以后武汉地区基坑工程的施工设计提供了值得借鉴的经验。