首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
This paper presents the analytical properties of the sensitivity of the two-dimensional, steady-state groundwater flow equation to the flow parameters and to the boundary conditions, based on the perturbation approach. These analytical properties are used to provide guidelines for model design, model calibration and monitoring network design. The sensitivity patterns are shown to depend on the nature of both the perturbed parameter and the variable investigated. Indeed, the sensitivity of the hydraulic head to the hydraulic conductivity extends mainly in the flow direction, while the sensitivity to the recharge spreads radially. Besides, the sensitivity of the flow longitudinal velocity to the hydraulic conductivity propagates in both the longitudinal and transverse directions, whereas the sensitivity of the flow transverse velocity propagates in the diagonal directions to the flow. The analytical results are confirmed by application examples on idealized and real-world simulations. These analytical findings allow some general rules to be established for model design, model calibration and monitoring network design. In particular, the optimal location of measurement points depends on the nature of the variable of interest. Measurement network design thus proves to be problem-dependent. Moreover, adequate monitoring well network design may allow to discriminate between the possible sources of error.  相似文献   

2.
An analytical and a numerical model are used to understand the response of velocity and sediment distributions over Gaussian-shaped estuarine cross-sections to changes in tidal forcing and water depth. The estuaries considered here are characterized by strong mixing and a relatively weak along-channel density gradient. It is also examined under what conditions the fast, two-dimensional analytical flow model yields results that agree with those obtained with the more complex three-dimensional numerical model. The analytical model reproduces and explains the main velocity and sediment characteristics in large parts of the parameter space considered (average tidal velocity amplitude, 0.1–1 m s − 1 and maximum water depth, 10–60 m). Its skills are lower for along-channel residual flows if nonlinearities are moderate to high (strong tides in deep estuaries) and for transverse flows and residual sediment concentrations if the Ekman number is small (weak tides in deep estuaries). An important new aspect of the analytical model is the incorporation of tidal variations in the across-channel density gradient, causing a double circulation pattern in the transverse flow during slack tides. The gradient also leads to a new tidally rectified residual flow component via net advection of along-channel tidal momentum by the density-induced transverse tidal flow. The component features landward currents in the channel and seaward currents over the slopes and is particularly effective in deeper water. It acts jointly with components induced by horizontal density differences, Coriolis-induced tidal rectification and Stokes discharge, resulting in different along-channel residual flow regimes. The residual across-channel density gradient is crucial for the residual transverse circulation and for the residual sediment concentration. The clockwise density-induced circulation traps sediment in the fresher water over the left slope (looking up-estuary in the northern hemisphere). Model results are largely consistent with available field data of well-mixed estuaries.  相似文献   

3.
Abstract

The stability of a plane parallel shear flow with the profile U(z) = tanh z is considered in a rotating system with the axis of rotation in the z-direction. The establishment of the basic flow requires a baroclinic state, but baroclinic effects are suppressed in the stability analysis by assuming a limit of high thermal conductivity. It is shown that the strongest growing disturbance changes from a purely transverse form in the limit of vanishing rotation rate to a nearly longitudinal form as the angular velocity of rotation increases. An analytical solution of the stability equation is obtained for vanishing growth rates of the transverse form of the instability. But, in general, the solution of the problem requires numerical integrations which demonstrate that the preferred direction of the wave vector of the instability is towards the left of the direction of the mean flow.  相似文献   

4.
A 2D depth‐averaged hydrodynamic, sediment transport and bed morphology model named STREMR HySeD is presented. The depth‐averaged sediment transport equations are derived from the 3D dilute, multiphase, flow equations and are incorporated into the hydrodynamic model STREMR. The hydrodynamic model includes a two‐equation turbulence model and a correction for the mean flow due to secondary flows. The suspended sediment load can be subdivided into different size classes using the continuum (two‐fluid) approach; however, only one bed sediment size is used herein. The validation of the model is presented by comparing the suspended sediment transport module against experimental measurements and analytical solutions for the case of equilibrium sediment‐laden in a transition from a rigid bed to a porous bed where re‐suspension of sediment is prevented. On the other hand, the bed‐load sediment transport and bed evolution numerical results are compared against bed equilibrium experimental results for the case of a meander bend. A sensitivity analysis based on the correction for secondary flow on the mean flow including the effect of secondary flow on bed shear stresses direction as well as the downward acceleration effect due to gravity on transverse bed slopes is performed and discussed. In general, acceptable agreement is found when comparing the numerical results obtained with STREMR HySeD against experimental measurements and analytical solutions. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

5.
The Xinanjiang model, which is a conceptual rainfall‐runoff model and has been successfully and widely applied in humid and semi‐humid regions in China, is coupled by the physically based kinematic wave method based on a digital drainage network. The kinematic wave Xinanjiang model (KWXAJ) uses topography and land use data to simulate runoff and overland flow routing. For the modelling, the catchment is subdivided into numerous hillslopes and consists of a raster grid of flow vectors that define the water flow directions. The Xinanjiang model simulates the runoff yield in each grid cell, and the kinematic wave approach is then applied to a ranked raster network. The grid‐based rainfall‐runoff model was applied to simulate basin‐scale water discharge from an 805‐km2 catchment of the Huaihe River, China. Rainfall and discharge records were available for the years 1984, 1985, 1987, 1998 and 1999. Eight flood events were used to calibrate the model's parameters and three other flood events were used to validate the grid‐based rainfall‐runoff model. A Manning's roughness via a linear flood depth relationship was suggested in this paper for improving flood forecasting. The calibration and validation results show that this model works well. A sensitivity analysis was further performed to evaluate the variation of topography (hillslopes) and land use parameters on catchment discharge. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

6.
Motivated by field studies of the Ems estuary which show longitudinal gradients in bottom sediment concentration as high as O(0.01 kg/m4), we develop an analytical model for estuarine residual circulation based on currents from salinity gradients, turbidity gradients, and freshwater discharge. Salinity is assumed to be vertically well mixed, while the vertical concentration profile is assumed to result from a balance between a constant settling velocity and turbulent diffusive flux. Width and depth of the model estuary are held constant. Model results show that turbidity gradients enhance tidally averaged circulation upstream of the estuarine turbidity maximum (ETM), but significantly reduce residual circulation downstream, where salinity and turbidity gradients oppose each other. We apply the condition of morphodynamic equilibrium (vanishing sediment transport) and develop an analytical solution for the position of the turbidity maximum and the distribution of suspended sediment concentration (SSC) along a longitudinal axis. A sensitivity study shows great variability in the longitudinal distribution of suspended sediment with the applied salinity gradient and six model parameters: settling velocity, vertical mixing, horizontal dispersion, total sediment supply, fresh water flow, and water depth. Increasing depth and settling velocity move the ETM upstream, while increasing freshwater discharge and vertical mixing move the ETM downstream. Moreover, the longitudinal distribution of SSC is inherently asymmetric around the ETM, and depends on spatial variations in the residual current structure and the vertical profile of SSC.  相似文献   

7.
Herein a simplified quasi-two dimensional horizontal hydro-morphological mathematical model is presented. The governing equations for the quasi-2D horizontal time-depending flow field are represented by the well-known approach of interconnected cells. New discharge laws between cells are incorporated. The model is capable of predicting temporal changes in water depth, velocity distribution,sediment transport, bed elevation, as well as water and suspended sediment exchanges between main stream and flood plains. An application of the model to the middle reach of the Argentinean Parana River is presented. Satisfactory results were obtained during model calibration, validation and application.  相似文献   

8.
1 INTRODUCTION Mathematical models to analyze the one-dimensional morphological evolution of alluvial rivers, induced either by natural events or anthropic actions, have been commonly applied since the original work of de Vries (1957, 1965, 1969). In the last decades much effort was made in developing suitable 2D horizontal and 3D time-depending mathematical models to study riverbed changes. In fact, full 2D-H hydro-morphological models (Olesen, 1987; Spasojevic and Holly, 1988) as w…  相似文献   

9.
Seismoelectric coupling in an electric isotropic and elastic anisotropic medium is developed using a primary–secondary formulation. The anisotropy is of vertical transverse isotropic type and concerns only the poroelastic parameters. Based on our finite difference time domain algorithm, we solve the seismoelectric response to an explosive source. The seismic wavefields are computed as the primary field. The electric field is then obtained as a secondary field by solving the Poisson equation for the electric potential. To test our numerical algorithm, we compared our seismoelectric numerical results with analytical results obtained from Pride's equation. The comparison shows that the numerical solution gives a good approximation to the analytical solution. We then simulate the seismoelectric wavefields in different models. Simulated results show that four types of seismic waves are generated in anisotropic poroelastic medium. These are the fast and slow longitudinal waves and two separable transverse waves. All of these seismic waves generate coseismic electric fields in a homogenous anisotropic poroelastic medium. The tortuosity has an effect on the propagation of the slow longitudinal wave. The snapshot of the slow longitudinal wave has an oval shape when the tortuosity is anisotropic, whereas it has a circular shape when the tortuosity is isotropic. In terms of the Thomsen parameters, the radiation anisotropy of the fast longitudinal wave is more sensitive to the value of ε, while the radiation anisotropy of the transverse wave is more sensitive to the value of δ.  相似文献   

10.
Shallow water depths on steep slopes of as much as fifty per cent can be measured easily by weighing a light flume and the water it contains. Because water accelerates along the flume, a good approximation of the steady state depth is obtained when the recording balance is fixed to its bottom end. From the unit discharge and the depth, and not from measurements of the surface velocity, the Darcy-Weisbach friction coefficient can be calculated. The present results show that this friction coefficient is larger in thin sheet flows than that calculated from the equation for rough turbulent flow. This latter could fit at a Reynolds Number of 50,000. When the regime is laminar (Re < 2,440) the Darcy-Weisbach friction coefficient always exceeds the theoretical value of 96/Re. The great relative depth of standing and travelling waves could account for this discrepancy together with turbulence and wake formation around bottom grains. Herein it is assumed that a regime can prevail where a laminar superlayer glides over a turbulent sublayer in the vicinity of bottom grains, because the ratio of the surface velocity to the mean velocity can greatly exceed 1.5, especially on steep slopes. Until photographs of the streamlines are taken, no statement about flow regimes in supercritical sheet flow can be made.  相似文献   

11.
The present paper focuses on the governing equations for the sensitivity of the variables to the parameters in flow models that can be described by one-dimensional scalar, hyperbolic conservation laws. The sensitivity is shown to obey a hyperbolic, scalar conservation law. The sensitivity is a conserved scalar except in the case of discontinuous flow solutions, where an extra, point source term must be added to the equations in order to enforce conservation. The propagation speed of the sensitivity waves being identical to that of the conserved variable in the original conservation law, the system of conservation laws formed by the original hyperbolic equation and the equation satisfied by the sensitivity is linearly degenerate. A consequence on the solution of the Riemann problem is that rarefaction waves for the variable of the original equation result in vacuum regions for the sensitivity. The numerical solution of the hyperbolic conservation law for the sensitivity by finite volume methods requires the implementation of a specific shock detection procedure. A set of necessary conditions is defined for the discretisation of the source term in the sensitivity equation. An application to the one-dimensional kinematic wave equation shows that the proposed numerical technique allows analytical solutions to be reproduced correctly. The computational examples show that first-order numerical schemes do not yield satisfactory numerical solutions in the neighbourhood of moving shocks and that higher-order schemes, such as the MUSCL scheme, should be used for sharp transients.  相似文献   

12.
Many simplifications are used in modeling surface runoff over a uniform slope. A very common simplification is to determine the infiltration rate independent of the overland flow depth and to combine it afterward with the kinematic-wave equation to determine the overland flow depth. Another simplication is to replace the spatially variable infiltration rates along the slope i(x, t) due to the water depth variations h(x,t) with an infiltration rate that is determined at a certain location along the slope. The aim of this study is to evaluate the errors induced by these simplications on predicted infiltration rates, overland flow depths, and total runoff volume. The error analysis is accomplished by comparing a simplified model with a model where the interaction between the overland flow depth and infiltration rate is counted. In this model, the infiltration rate is assumed to vary along the slope with the overland flow depth, even for homogeneous soil profiles. The kinematic-wave equation with interactive infiltration rate, calculated along the slopy by Richard's equation, are then solved by a finite difference scheme for a 100-m-long uniform slope. In the first error analysis, we study the effect of combining an ‘exact’ and ‘approximate’ one-dimensional infiltration rate with the kinematic-wave equation for three different soil surface roughness coefficients. The terms ‘exact’ and ‘approximate’ stand for the solution of Richard's equation with and without using the overland flow depth in the boundary condition, respectively. The simulations showed that higher infiltration rates and lower overland flow depths are obtained during the rising stage of the hydrograph when overland flow depth is used in the upper boundary condition of the one-dimensional Richard's equation. During the recession period, the simplified model predicts lower infiltration rates and higher overland flow depths. The absolute relative errors between the ‘exact’ and ‘approximate’ solutions are positively correlated to the overland flow depths which increase with the soil surface roughness coefficient. For this error analysis, the relative errors in surface runoff volume per unit slope width throughout the storm are much smaller than the relative errors in momentary overland flow depths and discharges due to the alternate signs of the deviations along the rising and falling stages. In the second error analysis, when the spatially variable infiltration rate along the slope i(x, t) is replaced in the kinematic-wave equation by i(t), calculated at the slope outlet, the overland flow depth is underestimated during the rising stage of the hydrograph and overestimated during the falling stage. The deviations during the rising stage are much smaller than the deviations during the falling stage, but they are of a longer duration. This occurs because the solution with i(x, t) recognizes that part of the slope becomes dry after rainfall stops, while overland flow still exists with i(t) determined at the slope outlet. As obtained for the first error analysis, the relative errors in surface runoff volume per unit slope width are also much smaller than the relative errors in momentary overland flow depths and discharges. The relation between the errors in overland flow depth and discharge to different mathematical simplifications enables to evaluate whether certain simplifications are justified or more computational efforts should be used.  相似文献   

13.
The conditions under which the Saint Venant equations system for unsteady open channel flow, as an initial–boundary value problem, becomes self‐similar are investigated by utilizing one‐parameter Lie group of point scaling transformations. One of the advantages of this methodology is that the self‐similarity conditions due to the initial and boundary conditions can also be investigated thoroughly in addition to the conditions due to the governing equation. The obtained self‐similarity conditions are compared with the scaling relationships that are derived through the Froude similitude. It is shown that the initial–boundary value problem of a one‐dimensional unsteady open channel flow process in a prototype domain can be self‐similar with that of several different scaled domains. However, the values of all the flow variables (at specified time and space) under different scaled domains can be upscaled to the same values in the prototype domain (at the corresponding time and space), as shown in this study. Distortion in scales of different space dimensions has been implemented extensively in physical hydraulic modelling, mainly because of cost, space and time limitations. Unlike the traditional approach, the distinction is made between the longitudinal–horizontal and transverse–horizontal length scales in this study. The scaled domain obtained by the proposed approach, when scaling ratios of channel width and water depth are equal, is particularly important for the similarity of flow characteristics in a cross‐section because the width‐to‐depth ratio and the inclination angles of the banks are conserved in a cross‐section. It is also shown that the scaling ratio of the roughness coefficient under distorted channel conditions depends on that of hydraulic radius and longitudinal length. The proposed scaling relations obtained by the Lie group scaling approach may provide additional spatial, temporal and economical flexibility in setting up physical hydraulic models. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

14.
The rate of biodegradation in contaminated aquifers depends to a large extent on dispersive mixing processes that are now generally accepted to result from spatial variations in the velocity field. It has been shown, however, that transient flow fields can also contribute to dispersive mixing. The influence of transient flow on biodegrading contaminants is particularly important since it can enhance mixing with electron acceptors, further promoting the reactive process. Using numerical simulations, the effect of transient flow on the behavior of a biodegradable contaminant is evaluated here both with respect to the development of apparently large horizontal transverse dispersion and also with respect to enhanced mixing between the substrate (electron donor) and electron acceptor. The numerical model BIO3D, which solves for advective-dispersive transport coupled with Monod-type biodegradation of substrates in the presence of an electron acceptor, was used for the simulations. The model was applied in a two-dimensional plan view mode considering a single substrate. Transient flow fields were found to yield larger apparent transverse dispersion because the longitudinal dispersivity also acts transverse to the mean flow direction. In the reactive case, the transient flow field increases substrate-oxygen mixing, which in turn enhances the overall rate of biodegradation. The results suggest that in the case of moderate changes of flow directions, a steady-state flow field can be justified, thereby avoiding the higher computational costs of a fully transient simulation. The use of a higher transverse horizontal dispersivity in a steady flow field can, under these conditions, adequately forecast plume development.  相似文献   

15.
Sensitivity analysis of the hydrological behaviour of basins has mainly focused on the correlation between streamflow and climate, ignoring the uncertainty of future climate and not utilizing complex hydrological models. However, groundwater storage is affected by climatic change and human activities. The streamflow of many basins is primarily sourced from the natural discharge of aquifers in upstream regions. The correlation between streamflow and groundwater storage has not been thoroughly discussed. In this study, the storage–discharge sensitivity of 22 basins in Taiwan was investigated by means of daily streamflow and rainfall data obtained over more than 30 years. The relationship between storage and discharge variance was evaluated using low‐flow recession analysis and a water balance equation that ignores the influence of rainfall and evapotranspiration. Based on the obtained storage–discharge sensitivity, this study explored whether the water storage and discharge behaviour of the studied basins is susceptible to climate change or human activities and discusses the regional differences in storage–discharge sensitivity. The results showed that the average storage–discharge sensitivities were 0.056 and 0.162 mm?1 in the northern and southern regions of Taiwan, respectively. In the central and eastern regions, the values were both 0.020 mm?1. The storage–discharge sensitivity was very high in the southern region. The regional differences in storage–discharge sensitivity with similar climate conditions are primarily due to differences in aquifer properties. Based on the recession curve, other factors responsible for these differences include land utilization, land coverage, and rainfall patterns during dry and wet seasons. These factors lead to differences in groundwater recharge and thus to regional differences in storage–discharge sensitivity.  相似文献   

16.
The occurrence of preferential flow in the subsurface has often been shown in field experiments. However, preferential flow is rarely included in models simulating the hydrological response at the catchment scale. If it is considered, preferential flow parameters are typically determined at the plot scale and then transferred to larger-scale simulations. Here, we successfully used the optimization algorithm DiffeRential Evolution Adaptive Metropolis (DREAM) to calibrate a 3D physics-based dual-permeability model directly at the catchment scale. In order to keep computational costs of the optimization routine at a reasonable level, we limited the number of parameters to be calibrated to the ones that had been shown before to be most influential for the simulation of discharge. We also calibrated parameters of the matrix domain and the macropore domain with a fixed parameter ratio between soil layers instead of calibrating every layer separately. These ratios reflected observed depth profiles of soil hydraulic properties at our study site. The dual-permeability parameter sets identified during calibration were able to simulate observed discharge time series satisfactorily but did not outperform a calibrated single-domain reference model scenario. Saturated hydraulic conductivities of the macropore domain were calibrated such that they became very similar to matrix saturated hydraulic conductivities, thereby effectively removing the effect of macropores. This suggests that the incorporation of vertical preferential flow as represented by the dual-permeability approach was not relevant for reproducing the hydrometric response reasonably well in the studied catchment. We also tested the scale-invariance of the calibrated dual-permeability parameter sets by using the parameter sets performing best at catchment scale to simulate plot-scale bromide depth profiles obtained from tracer irrigation experiments. This parameter transfer proved to be not successful, indicating that soil hydraulic parameters are scale-variant, independent of the direction of parameter transfer.  相似文献   

17.
A three-dimensional stochastic Lagrangian particle tracking sediment transport model is developed to solve the discrete advection-dispersion equation using a combination of empirical dispersion equations.The performance of three widely-used longitudinal dispersion coefficient equations was examined to select one of them as the primary dispersion equation term in the developed model. Also, a conditional empirical equation was used to consider the effect of vertical dispersion term in top layers n...  相似文献   

18.
An analytical solution is developed for flow in a manifold. The interest is primarily for trickle irrigation laterals, but the solution has broader applications including those for which pressure increases in the direction of flow and for intake manifolds. Both velocity head losses and variable discharge along the manifold are considered in the fundamental analysis. The appropriate second order, nonlinear equation is solved for two flow regimes, laminar and fully turbulent. Results indicate that for most trickle irrigation laterals the velocity head loss is negligible, but for an example from a chemical processing system the effect is important.  相似文献   

19.
Antidunes are fluvial bedforms that form in rivers with supercritical flows. The water surface over antidunes is strongly in phase with the bed surface, and the water surface is amplified to produce large surface waves. Many experimental studies have addressed antidunes; however, the shapes of three-dimensional antidunes in a wide channel with alternate bars have not yet been appropriately understood. In this study, we experimentally investigated the streamwise and transverse length scales of antidunes under conditions with a large width–depth ratio. Our experimental results provide evidence for the coevolution of antidunes and free alternate bars, and show for the first time that the development of free bars greatly alters the three-dimensional shape of water surface waves over antidunes. In the absence of free bars in a wide channel, multiple longitudinal wave trains form, and the number of wave trains counted in the transverse direction increases with increases in the width–depth ratio. However, the presence of free bars affects the local flow characteristics, resulting in a decrease of the number of wave trains in the transverse direction. Therefore, we propose a simple model for predicting the reduction in the number of wave trains by combining two previous theories for antidunes and free bars. Results obtained by the model were found to largely agree with experimental observations. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd  相似文献   

20.
为研究双轴各向异性介质多分量感应测井响应特征,本文基于三重傅里叶变换,推导任意方向偶极子源的谱域电磁场解析式;采用围线积分方法,自适应截断积分区间,结合谱域电磁场周期特性,实现三重傅里叶变换的精确快速积分;进而,针对双轴各向异性倾斜地层,模拟研究不同纵横向各向异性条件多分量感应测井响应特征.结果表明:利用谱域内电磁场在周向的周期特性简化解析式,可将计算速度提高4倍;自适应截断积分区间方法保证了计算精度,并极大地减少了积分节点数.对于倾斜双轴各向异性介质,倾角较大时,共面分量可反映地层横向各向异性,同轴分量可反映地层纵向各向异性;倾角较小时,同轴分量可反映地层横向各向异性,共面分量可反映地层纵向各向异性.  相似文献   

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

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