共查询到20条相似文献，搜索用时 15 毫秒
1.
《Advances in water resources》2005,28(8):779-792
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. 相似文献
2.
Olivier Gourgue Richard ComblenJonathan Lambrechts Tuomas KärnäVincent Legat Eric Deleersnijder 《Advances in water resources》2009
We present a flux-limiting wetting–drying approach for finite-element discretizations of the shallow-water equations using discontinuous linear elements for the elevation. The key ingredient of the method is the use of limiters for generalized nodal fluxes. This method is implemented into the Second-generation Louvain-la-Neuve Ice-ocean Model (SLIM), and is verified against standard test cases. The method is further applied to the wetting and drying of sand banks in the Scheldt Estuary, which is located in northern Belgium and the southern Netherlands. The results obtained for both the benchmarks and the realistic problem illustrate the accuracy of the method in describing the hydrodynamics in the vicinity of dry areas. In particular, the method strictly conserves mass, and there is no transport through dry areas. 相似文献
3.
An important part in the numerical simulation of tsunami and storm surge events is the accurate modeling of flooding and the appearance of dry areas when the water recedes. This paper proposes a new algorithm to model inundation events with piecewise linear Runge–Kutta discontinuous Galerkin approximations applied to the shallow water equations. This study is restricted to the one-dimensional case and shows a detailed analysis and the corresponding numerical treatment of the inundation problem.The main feature is a velocity based “limiting” of the momentum distribution in each cell, which prevents instabilities in case of wetting or drying situations. Additional limiting of the fluid depth ensures its positivity while preserving local mass conservation. A special flux modification in cells located at the wet/dry interface leads to a well-balanced method, which maintains the steady state at rest. The discontinuous Galerkin scheme is formulated in a nodal form using a Lagrange basis. The proposed wetting and drying treatment is verified with several numerical simulations. These test cases demonstrate the well-balancing property of the method and its stability in case of rapid transition of the wet/dry interface. We also verify the conservation of mass and investigate the convergence characteristics of the scheme. 相似文献
4.
Details are given of the development and application of a numerical model for predicting free-surface flows in estuarine and coastal basins using the finite volume method. Both second- and third-order accurate and oscillation free explicit numerical schemes have been used to solve the shallow water equations. The model deploys an unstructured triangular mesh and incorporates two types of mesh layouts, namely the ‘cell centred’ and ‘mesh vertex’ layouts, and provides a powerful mesh generator in which a user can adjust the mesh-size distribution interactively to create a desirable mesh. The quality of mesh has been shown to have a major impact on the overall performance of the numerical model.The model has been applied to simulate two-dimensional dam break flows for which transient water level distributions measured within a laboratory flume were available. In total 12 model runs were undertaken to test the model for various flow conditions. These conditions include: (1) different bed slopes (ranging from zero to 0.8%), (2) different upstream and downstream water level conditions, and (3) initially wet and dry bed conditions, downstream of the dam. Detailed comparisons have been made between model predicted and measured water levels and good agreement achieved between both sets of results. The model was then used to predict water level and velocity distributions in a real estuary, i.e. the Ribble Estuary, where the bed level varies rapidly at certain locations. In order to model the whole estuary, a 1-D numerical model has also been used to model the upper part of the estuary and this model was linked dynamically to the 2-D model. Findings from this application are given in detail. 相似文献
5.
During seismic wave propagation on a free surface, a strong material contrast boundary develops in response to interference by P- and S- waves to create a surfacewave phenomenon. To accurately determine the effects of this interface on surface-wave propagation, the boundary conditions must be accurately modeled. In this paper, we present a numerical approach based on the dynamic poroelasticity for a space–time-domain staggeredgrid finite-difference simulation in porous media that contain a free-surface boundary. We propose a generalized stess mirror formulation of the free-surface boundary for solids and fluids in porous media for the grid mesh on which lays the free-surface plane. Its analog is that used for elastic media, which is suitable for precise and stable Rayleigh-type surface-wave modeling. The results of our analysis of first kind of Rayleigh (R1) waves obtained by this model demonstrate that the discretization of the mesh in a similar way to that for elastic media can realize stable numerical solutions with acceptable precision. We present numerical examples demonstrating the efficiency and accuracy of our proposed method. 相似文献
6.
Yiping Wu Shuguang Liu Wende Yan Jiangzhou Xia Wenhua Xiang Kelin Wang Qiao Luo Wei Fu Wenping Yuan 《Stochastic Environmental Research and Risk Assessment (SERRA)》2016,30(1):225-235
Anthropogenic activities have altered the climate and led to changes in the water cycle. Understanding the climate change and hydrological responses is critical to derive adaptive strategies for sustainable water resources management. In this study, we diagnosed the trends of primary climate elements and hydrological components during the past half century (1960–2009) for the humid Xiangjiang River Basin in central-south China at multiple temporal and spatial scales. The air temperature trend demonstrated an overall warming climate but with a quicker pace in recent years; however, the wind speed reduced significantly in the early period, and this downtrend had largely disappeared after the mid-1990s. Under such a shifting climate, the hydrological responses were not monotonic during the past 50 years: the evapotranspiration behaved in a decreasing trend in the early 35 years (1960–1994), followed by an uptrend in the later period (1995–2009). The stepwise analysis of soil water content and baseflow demonstrated a wetting trend followed by a drying one but with a steeper slope, indicating an accelerated drying trend which may cause a concern in stream water availability especially in the dry season. Spatial trend analysis showed that some areas experienced a downtrend (drying) in the dry season, but most areas had an uptrend (wetting) in the wet season for the whole study period. Overall, the analyses of temporal and spatial changes are useful for decision makers to deal with the continuing changes in climate and hydrology. This study also highlighted the necessity of climate change studies at multiple temporal and spatial scales. 相似文献
7.
8.
A two‐dimensional (2D) finite‐difference shallow water model based on a second‐order hybrid type of total variation diminishing (TVD) approximate solver with a MUSCL limiter function was developed to model flooding and inundation problems where the evolution of the drying and wetting interface is numerically challenging. Both a minimum positive depth (MPD) scheme and a non‐MPD scheme were employed to handle the advancement of drying and wetting fronts. We used several model problems to verify the model, including a dam break in a slope channel, a dam break flooding over a triangular obstacle, an idealized circular dam‐break, and a tide flow over a mound. Computed results agreed well with the experiment data and other numerical results available. The model was then applied to simulate the dam breaking and flooding of Hsindien Creek, Taiwan, with the detailed river basin topography. Computed flooding scenarios show reasonable flow characteristics. Though the average speed of flooding is 6–7 m s?1, which corresponds to the subcritical flow condition (Fr < 1), the local maximum speed of flooding is 14·12 m s?1, which corresponds to the supercritical flow condition (Fr ≈ 1·31). It is necessary to conduct some kind of comparison of the numerical results with measurements/experiments in further studies. Nevertheless, the model exhibits its capability to capture the essential features of dam‐break flows with drying and wetting fronts. It also exhibits the potential to provide the basis for computationally efficient flood routing and warning information. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
9.
Soil moisture states,lateral flow,and streamflow generation in a semi‐arid,snowmelt‐driven catchment
Hydraulic connectivity on hillslopes and the existence of preferred soil moisture states in a catchment have important controls on runoff generation. In this study we investigate the relationships between soil moisture patterns, lateral hillslope flow, and streamflow generation in a semi‐arid, snowmelt‐driven catchment. We identify five soil moisture conditions that occur during a year and present a conceptual model based on field studies and computer simulations of how streamflow is generated with respect to the soil moisture conditions. The five soil moisture conditions are (1) a summer dry period, (2) a transitional fall wetting period, (3) a winter wet, low‐flux period, (4) a spring wet, high‐flux period, and (5) a transitional late‐spring drying period. Transitions between the periods are driven by changes in the water balance between rain, snow, snowmelt and evapotranspiration. Low rates of water input to the soil during the winter allow dry soil regions to persist at the soil–bedrock interface, which act as barriers to lateral flow. Once the dry‐soil flow barriers are wetted, whole‐slope hydraulic connectivity is established, lateral flow can occur, and upland soils are in direct connection with the near‐stream soil moisture. This whole‐slope connectivity can alter near‐stream hydraulics and modify the delivery of water, pressure, and solutes to the stream. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
10.
The basic idea of subgrid models is the use of available high-resolution bathymetric data at subgrid level in computations that are performed on relatively coarse grids allowing large time steps. For that purpose, an algorithm that correctly represents the precise mass balance in regions where wetting and drying occur was derived by Casulli (Int J Numer Method Fluids 60:391–408, 2009) and Casulli and Stelling (Int J Numer Method Fluids 67:441–449, 2010). Computational grid cells are permitted to be wet, partially wet, or dry, and no drying threshold is needed. Based on the subgrid technique, practical applications involving various scenarios were implemented including an operational forecast model for water level, salinity, and temperature of the Elbe Estuary in Germany. The grid generation procedure allows a detailed boundary fitting at subgrid level. The computational grid is made of flow-aligned quadrilaterals including few triangles where necessary. User-defined grid subdivision at subgrid level allows a correct representation of the volume up to measurement accuracy. Bottom friction requires a particular treatment. Based on the conveyance approach, an appropriate empirical correction was worked out. The aforementioned features make the subgrid technique very efficient, robust, and accurate. Comparison of predicted water levels with the comparatively highly resolved classical unstructured grid model shows very good agreement. The speedup in computational performance due to the use of the subgrid technique is about a factor of 20. A typical daily forecast can be carried out in less than 10 min on a standard PC-like hardware. The subgrid technique is therefore a promising framework to perform accurate temporal and spatial large-scale simulations of coastal and estuarine flow and transport processes at low computational cost. 相似文献
11.
Quadtree-adaptive tsunami modelling 总被引：3，自引：0，他引：3
Stéphane Popinet 《Ocean Dynamics》2011,61(9):1261-1285
12.
The partially decomposed organic layer (duff: F and H layers) of the forest floor is an important boundary between the soil and atmospheric processes. Here we use both empirical data and a three‐dimensional coupled heat and water budget model to explain the duff hydrological hillslope shift between very brief wet periods when lateral flow in the duff and infiltration into the mineral soil occur and dry periods when evaporative flow dominates and both lateral and mineral soil flow are not important. The duff moisture transitions from wet to dry periods were the result of low lateral flow which moves liquid and water vapour only centimetres to metres, very rapidly and mostly in the H layer immediately after precipitation. During wet periods, the net lateral fluxes were negative on divergent areas and positive on convergent areas of the hillslope, leading to a net moisture loss in divergent areas and a net gain in convergent areas. The response to lateral flow in the H layer was more rapid than in the F layer. The transition from the lateral downwards flow to mineral soil to evaporative control was within approximately 48 h of precipitation. Canopy species and aspect were important with lodgepole pine, southwest aspect and 4‐cm deep duff controlled by evaporative processes while Engelmann spruce, northeast aspect and 30‐cm duff were more controlled by hillslope redistribution processes. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
13.
为了研究浅层膨胀土的细观结构在湿干循环过程中的演化规律,研制了可以进行CT扫描的CT-固结仪。在不同压力作用下对南阳膨胀土的原状样及其重塑土样进行多组多次湿干循环试验,并对每次烘干和增湿稳定后的试样进行CT扫描,跟踪观察湿干循环过程中试样的细观结构变化。结果表明:在湿干循环中,土样的膨胀量和收缩量都比较大,经历三次湿干循环之后,相同饱和度下土样体积基本不再变化;增湿和干燥都能使原状膨胀土及其重塑土的结构损伤,第一次烘干之后,土样即产生明显裂隙;随着湿干循环次数的增加,土样裂隙继续开展,土样结构损伤发生累积;上部荷载对裂隙开展有一定抑制作用,并影响裂隙的发育和裂隙网络的形状;原状膨胀土的裂隙围绕第一次形成的主裂隙发育,而重塑膨胀土的裂隙呈龟背状或辐射状;原状膨胀土及其重塑土在湿态(饱和度Sr=85%)时的CT数ME值均较干态(饱和度Sr=25%)时的ME值大;湿态和干态的ME值均随湿干循环次数呈近似线性变化,前者的坡度比较平缓,后者则较陡;基于CT数据定义了膨胀土的细观结构参数,提出定量描述湿干循环过程中膨胀土细观结构参数与试样饱和度、湿干循环次数及所受荷载的数学表达式,预测结果与试验资料比较接近。 相似文献
14.
The effects of floodplain vegetation on river planform have been investigated for a medium‐sized river using a 2D morphodynamic model with submodels for flow resistance and plant colonization. The flow resistance was divided into a resistance exerted by the soil and a resistance exerted by the plants. In this way it was possible to reproduce both the decrease in bed shear stress, reducing the sediment transport capacity of the flow within the plants, and the increase in hydraulic resistance, reducing the flow velocities. Colonization by plants was obtained by instantaneously assigning vegetation to the areas that became dry at low water stages. This colonization presents a step forward in the modelling of bank accretion. Bank erosion was related to bed degradation at adjacent wet cells. Bank advance and retreat were reproduced as drying and wetting of the computational cells at the channel margins. The model was applied to a hypothetical case with the same characteristics as the Allier River (France). The river was allowed to develop its own geometry starting from a straight, uniform, channel. Different vegetation densities produced different planforms. With bare floodplains, the river always developed a braided planform, even if the discharge was constant and below bankfull. With the highest vegetation density (grass) the flow concentrated in a single channel and formed incipient meanders. Lower vegetation density (pioneer vegetation) led to a transitional planform, with a low degree of braiding and distinguishable incipient meanders. The results comply with flume experiments and field observations reported in the literature. 相似文献
15.
A periodic vertical movement of the groundwater table results in a subsequent cyclic response of the water content and pressure profiles in the vadose zone. The sequence of periodic wetting and drying processes can be affected by hysteresis effects in this zone. A one-dimensional saturated/unsaturated flow model based on Richards’ equation and the Mualem (Soil Sci. 137 (1984) 283) hysteresis model is formulated which can take into account multi-cycle hysteresis effects in the relation between capillary pressure and water content. The numerical integration of the unsaturated flow equation is based on a Galerkin-type finite element method. The flow domain is discretised by finite elements with linear shape functions. Simulations start with static water content and pressure profiles, which correspond to either a boundary drying or wetting retention curve. To facilitate the numerical solution of the hysteretic case an implicit non-iterative procedure was chosen for the solution of the nonlinear differential equation. Laboratory experiments were performed with a vertical sand column by imposing a high frequency periodic pressure head at the lower end of the column. The total water volume in the column, and the periodic water content profile averaged over time were measured. The boundary drying and wetting curves of the relation between water content and capillary pressure were determined by independent experiments. The simulations of the experimental conditions show a clear effect of the hysteresis phenomenon on the water content profile. The simulations with hysteresis agree well with the measurements. Computed dimensionless water content profiles are presented for different oscillation frequencies with and without consideration of hysteresis. 相似文献
16.
K.A. Culligan D. Wildenschild B.S.B. Christensen W.G. Gray M.L. Rivers 《Advances in water resources》2006
Studies of NAPL dissolution in porous media have demonstrated that measurement of saturation alone is insufficient to describe the rate of dissolution. Quantification of the NAPL–water interfacial area provides a measure of the expected area available for mass transfer and will likely be a primary determinant of NAPL removal efficiency. To measure the interfacial area, we have used a synchrotron-based CMT technique to obtain high-resolution 3D images of flow in a Soltrol–water–glass bead system. The interfacial area is found to increase as the wetting phase saturation decreases, reach a maximum, and then decrease as the wetting phase saturation goes to zero. These results are compared to previous findings for an air–water–glass bead study; The Soltrol–water interfacial areas were found to peak at similar saturations as those measured for the air–water system (20–35% saturation range), however, the peak values were in some cases almost twice as high for the oil-water system. We believe that the observed differences between the air–water and oil–water systems to a large degree can be explained by the differences in interfacial tensions for the two systems. 相似文献
17.
M. ALIPARAST 《国际泥沙研究》2009,24(1):99-107
A numerical model based upon a second-order upwind ceil-center f＇mite volume method on unstructured triangular grids is developed for solving shallow water equations. The assumption of a small depth downstream instead of a dry bed situation changes the wave structure and the propagation speed of the front which leads to incorrect results. The use of Harten-Lax-vau Leer （HLL） allows handling of wet/dry treatment. By usage of the HLL approximate Riemann solver, also it make possible to handle discontinuous solutions. As the assumption of a very small depth downstream oftbe dam can change the nature of the dam break flow problem which leads to incorrect results, the HLL approximate Riemann solver is used for the computation of inviscid flux functions, which makes it possible to handle discontinuous solutions. A multidimensional slope-limiting technique is applied to achieve second-order spatial accuracy and to prevent spurious oscillations. To alleviate the problems associated with numerical instabilities due to small water depths near a wet/dry boundary, the friction source terms are treated in a fully implicit way. A third-order Runge-Kutta method is used for the time integration of semi-discrete equations. The developed numerical model has been applied to several test cases as well as to real flows. The tests are tested in two cases： oblique hydraulic jump and experimental dam break in converging-diverging flume. Numerical tests proved the robustness and accuracy of the model. The model has been applied for simulation of dam break analysis of Torogh in Iran. And finally the results have been used in preparing EAP （Emergency Action Plan）. 相似文献
18.
The no-slip boundary condition has usually been assumed to hold for the Reynolds equations (local cubic law) for fluid flow through rough-walled fractures. However, its validity for non-wetting fluid flow, such as prevails in fractured oil reservoirs, has been questioned. A series of experiments with a rough-walled fracture with mean aperture of 760 μm finds a higher flow rate for non-wetting fluid than wetting fluid. A modified Reynolds equation with a slip boundary condition is derived for non-wetting fluid flow through rough-walled fractures. Comparison of the modified Reynolds equation predictions with experimental results confirms that slip was a plausible explanation for a higher flow rate. The amount by which the flow rate for non-wetting fluids exceeds that for wetting fluids is found to depend highly on, and increase with, the degree to which the flowing fluid was non-wetted to thin immobile films on the surfaces. Numerical studies using the modified Reynolds equation indicate that the flow rate of non-wetting fluid became higher than that of wetting fluid as the roughness of the fracture increases. As the aperture becomes smaller, the flow rate ratio of non-wetting fluid to wetting fluid becomes large, leading to the endpoint relative permeability for the non-wetting fluid to exceed 1. The experimental and numerical studies clearly show that as the aperture of the fracture became less than a few hundred microns, the modified Reynolds equation with slip boundary conditions provides a better model for flow of a non-wetting fluid through rough-walled fractures. 相似文献
19.
A robust well-balanced finite volume model for shallow water flows with wetting and drying over irregular terrain 总被引：2，自引：0，他引：2
An unstructured Godunov-type finite volume model is developed for the numerical simulation of geometrically challenging two-dimensional shallow water flows with wetting and drying over convoluted topography. In the framework of sloping bottom model, a modified formulation of shallow water equations is used to preserve mass conservation during flooding and recession. The key ingredient of the model is the use of this combination of the sloping bottom model and the modified shallow water equations to provide a robust technique for wet/dry fronts tracking and, together with centered discretization of the bed slope source term, to exactly preserve the static flow on irregular topographies. The variable reconstruction technique ensures nonnegative reconstructed water depth and reasonable reconstructed velocity, and the friction terms are solved by semi-implicit scheme that does not invert the direction of velocity components. The robustness and accuracy of the proposed model are assessed by comparing numerical and reference results of extensive test cases. Moreover, the results of a dam-break flooding over real topography are presented to show the capability of the model on field-scale application. 相似文献
20.
A new solution of transient confined–unconfined flow driven by a pumping well is developed and compared to previous approximate solutions of Moench and Prickett [Moench AF, Prickett TA. Radial flow in an infinite aquifer undergoing conversion from artesian to water table conditions. Water Resour Res 1972;8:494–9] and Hu and Chen [Hu L, Chen C. Analytical methods for transient flow to a well in a confined–unconfined aquifer. Ground Water 2008;46(4):642–6]. The problem is rewritten in dimensionless form with the Boltzmann transform. The nonlinear equation for flow in the unconfined zone is solved with the Runge–Kutta method. Position of the conversion interface is determined with an iteration scheme. This study shows that the confined–unconfined flow depends on three dimensionless parameters that represent the confined–unconfined storativity ratio (aD), the ratio of the initial hydraulic head over the aquifer thickness (fi), and the dimensionless pumping rate (qD). The rate of expansion of the unconfined zone increases with qD, but decreases with aD and fi. Differences between the two previous approximate solutions and the new solution of this study are observable in the estimated position of the conversion interface and the drawdown–time curves. The new solution can be applied to estimate the time for confined–unconfined conversion to occur (critical conversion time), and the time when the pumping well becomes dry (critical drying time). The critical conversion time is found to be very sensitive to the initial hydraulic head. The critical drying time is often much larger than the critical conversion time and may never be observed during a finite pumping period. 相似文献