首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
This paper presents a computational fluid dynamics model for simulation of twodimensional (2-D) water flow, sediment transport, bank failure processes, and the subsequent channel pattern changes. Effects of secondary currents at channel bends are included in the modified momentum conservation equation of water flow. An improved bank failure model is applied to calculate bank failure due to riverbed erosion, and to simulate lateral migration and planform changes of alluvial channels. The water flow model has been validated using laboratory measurements of flow in consecutive bends designed by the authors, in addition to flume test data from the literature.  相似文献   

2.
Numerical simulation of channel pattern changes Part I: Mathematical model   总被引:1,自引:0,他引:1  
This paper presents a computational fluid dynamics model for simulation of two- dimensional (2-D) water flow, sediment transport, bank failure processes, and the subsequent channel pattern changes. Effects of secondary currents at channel bends are included in the modified momentum conservation equation of water flow. An improved bank failure model is applied to calculate bank failure due to riverbed erosion, and to simulate lateral migration and planform changes of alluvial channels. The water flow model has been validated using laboratory measurements of flow in consecutive bends designed by the authors, in addition to flume test data from the literature.  相似文献   

3.
Weiming Wu 《Ocean Dynamics》2014,64(7):1061-1071
A 3-D shallow-water flow model has been developed to simulate the flow in coastal vegetated waters with short waves. The model adopts the 3-D phase-averaged shallow-water flow equations with radiation stresses induced by short waves. It solves the governing equations using an implicit finite volume method based on quadtree rectangular mesh in the horizontal plane and stretching mesh in the vertical direction. The flow model is coupled with a spectral wave deformation model called CMS-Wave. The wave model solves the spectral wave-action balance equation and provides wave characteristics to the flow model. The model considers the effects of vegetation on currents and waves by including the drag and inertia forces of vegetation in the momentum equations and the wave energy loss due to vegetation resistance in the wave-action balance equation. The model has been tested using several sets of laboratory experiments, including steady flows in a straight channel with submerged vegetation and in a compound channel with vegetated floodplain and random waves through a vegetated channel and on a vegetated beach slope. The calculated water levels, current velocities, and wave heights are in general good agreement with the measured data.  相似文献   

4.
1 INTRODUCTIONFor many hydraulic engineering problems, the analysis of flow and bed level variations in openchannels is a fundamental prerequisite. forcal methOds fOr alluvial rivers are well develoPednowadays as far as onediInensional descriPtions are concemed. A cOmPrhensive analysis of Ihe wellknown models is Presented by Habersack(l998). HOwever, for a number of Problems such as channelwidening, flow pattem close to sPuds and etc. a more deailed knowledge of the bed level behavio…  相似文献   

5.
We propose an improvement of the overland‐flow parameterization in a distributed hydrological model, which uses a constant horizontal grid resolution and employs the kinematic wave approximation for both hillslope and river channel flow. The standard parameterization lacks any channel flow characteristics for rivers, which results in reduced river flow velocities for streams narrower than the horizontal grid resolution. Moreover, the surface areas, through which these wider model rivers may exchange water with the subsurface, are larger than the real river channels potentially leading to unrealistic vertical flows. We propose an approximation of the subscale channel flow by scaling Manning's roughness in the kinematic wave formulation via a relationship between river width and grid cell size, following a simplified version of the Barré de Saint‐Venant equations (Manning–Strickler equations). The too large exchange areas between model rivers and the subsurface are compensated by a grid resolution‐dependent scaling of the infiltration/exfiltration rate across river beds. We test both scaling approaches in the integrated hydrological model ParFlow. An empirical relation is used for estimating the true river width from the mean annual discharge. Our simulations show that the scaling of the roughness coefficient and the hydraulic conductivity effectively corrects overland flow velocities calculated on the coarse grid leading to a better representation of flood waves in the river channels.  相似文献   

6.
The helical flow significantly affects the flow, sediment transport and morphological evolution in curved channels. A semi-empirical formula is proposed to determine the cross-stream distribution of the helical flow intensity in the developed regions of a channel bend. It is then used to evaluate the dispersion terms in the depth-averaged 2-D momentum equations and suspended-load transport equation as well as the bed-load transport angle, thus enhancing the depth-averaged 2-D model to account for the effect of helical flow. The tests in several experimental and field cases show that the enhanced depth-averaged 2-D model can much more reasonably predict the shifting of main flow from inner bank to outer bank, the erosion along outer bank and deposition along inner bank than the depth-averaged 2-D model without considering this effect.  相似文献   

7.
1 wrsoooCTIoxThe Yeuow mver crtes a huge amoun of sedimcht and the noods often cause raPid and severeerosinn and dePOsihon. The channl bed of the YelOw mver often exPeriences degIadation in the mainchannel during fioods. In some cases vigorous erosion uP to l0 meters takes place in a shOrt Period oftiIn. Such phenomenon usually occurs in the ndddie reaChs of the Yelow mver and its tributaries suchas the Wdse mVer the Beiluohe mver etc. For examPle, th6 hyPenconcentraed nood in July l…  相似文献   

8.
The behaviour of river waves is described using a simplified dimensionless form of the momentum equation in conjunction with the continuity equation. Three dimensionless parameters were derived based on a quantitative linear analysis. These parameters, which depend on the Froude number of the steady uniform flow and the geometric characteristics of the river, permit quantification of the influence of inertia and pressure in the momentum equation. It was found that dynamic and diffusion waves occur mainly on gentle channel slopes and the transition between them is characterized by the Froude number. On the other hand, the kinematic wave has a wide range of applications. If the channel slope is greater than 1%, the kinematic wave is particularly suitable for describing the hydraulics of flow. Since slopes in natural channel networks are often greater than 1%, an analytical solution of the linearized kinematic wave equation with lateral inflow uniformly distributed along the channel is desirable and was therefore derived. The analytical solution was then implemented in a channel routing module of an existing simple rainfall–runoff model. The results obtained using the analytical solution compared well with those obtained from a non‐linear kinematic wave model. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

9.
I.INTRODUCTIONTwo-dimensionalnumericalmodelisaPOwerfoltoolforengineersandriVermanagerstopredictfloodhydxaulics,identifyareasofinundation,anddesignoptionsforfloodcontrollingstructures.SomespecialproblemswithheavilysedimentladenflowriVershouldbecarefullyconsideredforthenumericalmodeldesigning;1.theplaneformofariVerisusuallybraidedanditsmainchannelshiftsoften.Themainchannelandbarreplaceeachotherseveraltimesinonefloodevent.Atagivenlocationthewaterdepthmaychangefromover10meterstoseveralcenhm…  相似文献   

10.
Various parameters such as bed and bank materials, shape and irregularity of the section, vegetation, river meanders, plan of the river path etc. affect the flow hydraulic resistance. In open channel hydraulics the effects of all these parameters are generally considered as the roughness coefficient. The Manning’s equation is one of the most practical equations to flow resistance analysis, in which the surface roughness is defined by Manning coefficient. Since many parameters are effective on the value of this coefficient, in this research study it was tried to define the roughness coefficient somehow that it be able to dynamically change with different river and hydraulic conditions. The collected data in Karun River (Iran) for two periods were used as the case study. It is shown that the accuracy of model predictions for water surface elevations were improved more than 13% in error estimation in comparison with the corresponding results obtained for a constant roughness coefficient. The roughness coefficient (n) for Karun River was also estimated using the empirical method proposed by Cowan for two different dry and wet periods. These values were then successfully compared with the average corresponding roughness coefficients calculated by the numerical model for those periods.  相似文献   

11.
V. P. Singh 《水文研究》1995,9(7):783-796
Error equations for the kinematic wave and diffusion wave approximations with lateral inflow neglected in the momentum equation are derived under simplified conditions for space-independent flows. These equations specify error as a function of time in the flow hydrograph. The kinematic wave, diffusion wave and dynamic wave solutions are parameterized through a dimensionless parameter γ which is dependent on the initial conditions. This parameter reflects the effect of initial flow depth, channel-bed slope, lateral inflow, infiltration and channel roughness when the initial condition is non-vanishing; it reflects the effect of bed slope, channel roughness and acceleration due to gravity when the initial condition is vanishing. The error equations are found to be the Riccati equation. The structure of the error equations in the case when the momentum equation neglects lateral inflow is different from that when the lateral inflow is included.  相似文献   

12.
V. P. Singh 《水文研究》1994,8(4):311-326
Error equations for the kinematic wave and diffusion wave approximations with lateral inflow neglected in the momentum equation are derived under simplified conditions for space-independent flows. These equations specify error as a function of time in the flow hydrograph. The kinematic wave, diffusion wave and dynamic wave solutions are parameterized through a dimensionless parameter γ which is dependent on the initial conditions. This parameter reflects the effect of initial flow depth, channel-bed slope, lateral inflow and channel roughness when the initial condition is non-vanishing; and it reflects the effect of bed slope, channel roughness and acceleration due to gravity when the initial condition is vanishing. The error equations are found to be the Riccati equation. The structure of the error equations in the case when the momentum equation neglects lateral inflow is different from that when the lateral inflow is included.  相似文献   

13.
Mountainous areas are characterized by steep slopes and rocky landforms, with hydrological conditions varying rapidly from upstream to downstream, creating variable interactions between groundwater and surface water. In this study, mechanisms of groundwater–surface water interactions within a headwater catchment of the North China Plain were assessed along the stream length and during different seasons, using hydrochemical and stable isotope data, and groundwater residence times estimated using chlorofluorocarbons. These tracers indicate that the river is gaining, due to groundwater discharge in the headwater catchment both in the dry and rainy seasons. Residence time estimation of groundwater using chlorofluorocarbons data reveals that groundwater flow in the shallow sedimentary aquifer is dominated by the binary mixing of water approximating a piston flow model along 2 flow paths: old water, carried by a regional flow system along the direction of river flow, along with young water, which enters the river through local flow systems from hilly areas adjacent to the river valley (particularly during the rainy season). The larger mixing ratio of young water from lateral groundwater recharge and return flow of irrigation during the rainy season result in higher ion concentrations in groundwater than in the dry season. The binary mixing model showed that the ratio of young water versus total groundwater ranged from 0.88 to 0.22 and 1.0 to 0.74 in the upper and lower reaches, respectively. In the middle reach, meandering stream morphology allows some loss of river water back into the aquifer, leading to increasing estimates of the ratio of young water (from 0.22 to 1). This is also explained by declining groundwater levels near the river, due to groundwater extraction for agricultural irrigation. The switch from a greater predominance of regional flow in the dry season, to more localized groundwater flow paths in the wet season is an important groundwater–surface water interactions mechanism, with important catchment management implications.  相似文献   

14.
15.
针对现有的河道水流洪水演算模型只能模拟单一变量(流量或水位)的问题,以水流连续方程和河段蓄水量的两种不同表达形式(蓄水量等于平均过水断面面积与河段长乘积,蓄水量等于河段平均流量与传播时间的乘积)为基础,对马斯京根模型进行了通用性改进,提出了双变量耦合通用演算模型.选取了四大水系(包括内陆河流和入海河流)的16个河段汛期洪水资料进行模型检验,模型验证考虑了地理范围、不同的河段特征和水力特征、洪水量级等因素,全面地检验了模型结构的合理性和模拟实际洪水的有效性.将双变量耦合通用演算模型与传统的马斯京根法进行了效果比较,结果表明双变量耦合通用演算模型的模拟精度高于马斯京根法,模拟效果比马斯京根法稳定一些,而且具有较好的通用性.  相似文献   

16.
Mathematical modelling of overland flow is a critical task in simulating transport of water, sediment and other pollutants from land surfaces to receiving waters. In this paper, an overland flow routing method is developed based on the Saint‐Venant equations using a discretized hillslope system for areas with high roughness and steep slope. Under these conditions, the momentum equation reduces to a unique relationship between the flow depth and discharge. A hillslope is treated as a system divided into several subplanes. A set of first‐order non‐linear differential equations for subsequent subplanes are solved analytically using Chezy's formula in lieu of the momentum equation. Comparison of the analytical solution of the first‐order non‐linear ordinary differential equations and a numerical solution using the Runge‐Kutta method shows a relative error of 0·3%. Using runoff data reported in the literature, comparison between the new approach and a numerical solution of the full Saint‐Venant equations showed a close agreement. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

17.
An approach to represent drying and wetting processes in a three-dimensional finite element sigma coordinate model is described. This approach makes use of capillaries in dry areas, which can connect to the nearby wet areas. The time marching of the mass conservation equation is modified by introducing a “size factor” coefficient and a water level diffusion term. Therefore, the fictitious water level of the dry nodes can fluctuate with the adjacent wet nodes. This eliminates the artificial pressure gradient appearing in some drying and wetting approaches in the partially wet (transition) elements. This approach results in a null momentum computation at the dry areas, which can guarantee numerical stability and satisfy the mass and momentum conservation. The approach has been applied in a hypothetical case and a real case in Xiamen Estuary, China, with satisfactory results.  相似文献   

18.
This study numerically investigates effects of cutting riparian vegetation on flow characteristics by using a two-dimensional numerical model. The numerical model is based on depth-averaging the time- and volume-averaged Navier–Stokes equation with turbulent effects determined by the standard kε turbulence model. Drag forces exerted by the flow on vegetation are considered by adding source terms into momentum equations. In a rectangular channel and compound channel with vegetation along one side, numerical predictions show are in good agreement with those of previous studies. Five cutting scenarios, including the original, cutting along the main channel side, cutting along the bank side, alternative cutting, and reducing vegetative density, are analyzed in this study. The influences of the cutting scenarios on hydrodynamic behaviors are evaluated via numerical simulations. Simulation results suggest that cutting along the main channel side is the most effective scenario for reducing water depth and flow velocities.  相似文献   

19.
The basic aim of this paper is to formulate rigorous conservation equations for mass, momentum, energy and entropy for a watershed organized around the channel network. The approach adopted is based on the subdivision of the whole watershed into smaller discrete units, called representative elementary watersheds (REW), and the formulation of conservation equations for these REWs. The REW as a spatial domain is divided into five different subregions: (1) unsaturated zone; (2) saturated zone; (3) concentrated overland flow; (4) saturated overland flow; and (5) channel reach. These subregions all occupy separate volumina. Within the REW, the subregions interact with each other, with the atmosphere on top and with the groundwater or impermeable strata at the bottom, and are characterized by typical flow time scales.The balance equations are derived for water, solid and air phases in the unsaturated zone, water and solid phases in the saturated zone and only the water phase in the two overland flow zones and the channel. In this way REW-scale balance equations, and respective exchange terms for mass, momentum, energy and entropy between neighbouring subregions and phases, are obtained. Averaging of the balance equations over time allows to keep the theory general such that the hydrologic system can be studied over a range of time scales. Finally, the entropy inequality for the entire watershed as an ensemble of subregions is derived as constraint-type relationship for the development of constitutive relationships, which are necessary for the closure of the problem. The exploitation of the second law and the derivation of constitutive equations for specific types of watersheds will be the subject of a subsequent paper.  相似文献   

20.
A depth-averaged 2-D numerical model for unsteady flow, salinity and cohesive sediment transport in estuaries is established using the finite volume method on the non-staggered, curvilinear grid. The convection terms are discretized by upwind schemes, the diffusion terms are by the central difference scheme, and the time derivative terms are by the three-time-level implicit scheme. The coupling of flow velocity and water level in the 2-D shallow water equations is achieved by the SIMPLEC algorithm with the Rhie and Chow's momentum interpolation method. The sediment model calculates the settling, deposition, erosion and transport of cohesive sediment, taking into account the influence of sediment size, sediment concentration, salinity and turbulence intensity on the flocculation of cohesive sediment. The flow model is first tested against the measurement data in the Tokyo Bay and San Francisco Bay, showing good agreements. And then, the entire model of flow, salinity and sediment transport is verified in the Gironde Estuary. The water elevation, flow velocity, salinity and sediment concentration are well predicted.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号