首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 371 毫秒
1.
D. Yu  S. N. Lane 《水文研究》2006,20(7):1567-1583
This paper develops and tests a sub‐grid‐scale wetting and drying correction for use with two‐dimensional diffusion‐wave models of urban flood inundation. The method recognizes explicitly that representations of sub‐grid‐scale topography using roughness parameters will provide an inadequate representation of the effects of structural elements on the floodplain (e.g. buildings, walls), as such elements not only act as momentum sinks, but also have mass blockage effects. The latter may dominate, especially in structurally complex urban areas. The approach developed uses high‐resolution topographic data to develop explicit parameterization of sub‐grid‐scale topographic variability to represent both the volume of a grid cell that can be occupied by the flow and the effect of that variability upon the timing and direction of the lateral fluxes. This approach is found to give significantly better prediction of fluvial flood inundation in urban areas than traditional calibration of sub‐grid‐scale effects using Manning's n. In particular, it simultaneously reduces the need to use exceptionally high values of n to represent the effects of using a coarser mesh process representation and increases the sensitivity of model predictions to variation in n. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

2.
Time nonlocal transport models such as the time fractional advection‐dispersion equation (t‐fADE) were proposed to capture well‐documented non‐Fickian dynamics for conservative solutes transport in heterogeneous media, with the underlying assumption that the time nonlocality (which means that the current concentration change is affected by previous concentration load) embedded in the physical models can release the effective dispersion coefficient from scale dependency. This assumption, however, has never been systematically examined using real data. This study fills this historical knowledge gap by capturing non‐Fickian transport (likely due to solute retention) documented in the literature (Huang et al. 1995) and observed in our laboratory from small to intermediate spatial scale using the promising, tempered t‐fADE model. Fitting exercises show that the effective dispersion coefficient in the t‐fADE, although differing subtly from the dispersion coefficient in the standard advection‐dispersion equation, increases nonlinearly with the travel distance (varying from 0.5 to 12 m) for both heterogeneous and macroscopically homogeneous sand columns. Further analysis reveals that, while solute retention in relatively immobile zones can be efficiently captured by the time nonlocal parameters in the t‐fADE, the motion‐independent solute movement in the mobile zone is affected by the spatial evolution of local velocities in the host medium, resulting in a scale‐dependent dispersion coefficient. The same result may be found for the other standard time nonlocal transport models that separate solute retention and jumps (i.e., displacement). Therefore, the t‐fADE with a constant dispersion coefficient cannot capture scale‐dependent dispersion in saturated porous media, challenging the application for stochastic hydrogeology methods in quantifying real‐world, preasymptotic transport. Hence improvements on time nonlocal models using, for example, the novel subordination approach are necessary to incorporate the spatial evolution of local velocities without adding cumbersome parameters.  相似文献   

3.
Many studies indicate that small‐scale heterogeneity and/or mobile–immobile mass exchange produce transient non‐Fickian plume behavior that is not well captured by the use of the standard, deterministic advection‐dispersion equation (ADE). An extended ADE modeling framework is presented here that is based on continuous time random walk theory. It can be used to characterize non‐Fickian transport coupled with simultaneous sequential first‐order reactions (e.g., biodegradation or radioactive decay) for multiple degrading contaminants such as chlorinated solvents, royal demolition explosive, pesticides, and radionuclides. To demonstrate this modeling framework, new transient analytical solutions are derived and are inverted in Laplace space. Closed‐form, steady‐state, multi‐species analytical solutions are also derived for non‐Fickian transport in highly heterogeneous aquifers with linear sorption–desorption and matrix diffusion for use in spreadsheets. The solutions are general enough to allow different degradation rates for the mobile and immobile zones. The transient solutions for multi‐species transport are applied to examine the effects of source remediation on the natural attenuation of downgradient plumes of both parent and degradation products in highly heterogeneous aquifers. Results for representative settings show that the use of the standard, deterministic ADE can over‐estimate cleanup rates and under‐predict the cleanup timeframe in comparison to the extended ADE analytical model. The modeling framework and calculations introduced here are also applied for a 30 year groundwater cleanup program at a site in Palm Bay, Florida. The simulated plume concentrations using the extended ADE exhibited agreement with observed long concentration tails of trichloroethene, cis 1,2 DCE, and VC that remained above cleanup goals.  相似文献   

4.
The processes of hillslope runoff and erosion are typically represented at coarse spatial resolution in catchment‐scale models due to computational limitations. Such representation typically fails to incorporate the important effects of topographic heterogeneity on runoff generation, overland flow, and soil erosion. These limitations currently undermine the application of distributed catchment models to understand the importance of thresholds and connectivity on hillslope and catchment‐scale runoff and erosion, particularly in semi‐arid environments. This paper presents a method for incorporating high‐resolution topographic data to improve sub‐grid scale parameterization of hillslope overland flow and erosion models. Results derived from simulations conducted using a kinematic wave overland flow model at 0.5 m spatial resolution are used to parameterize the depth–discharge relationship in the overland flow model when applied at 16 m resolution. The high‐resolution simulations are also used to derive a more realistic parameterization of excess flow shear stress for use in the 16 m resolution erosion model. Incorporating the sub‐grid scale parameterization in the coarse‐resolution model (16 m) leads to improved predictions of overland flow and erosion when evaluated using results derived from high‐resolution (0.5 m) model simulations. The improvement in performance is observed for a range of event magnitudes and is most notable for erosion estimates due to the non‐linear dependency between the rates of erosion and overland flow. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

5.
I. Haltas 《水文研究》2012,26(22):3448-3458
Recognizing the spatial heterogeneity of hydraulic parameters, many researchers have studied the solute transport by both groundwater and channel flow in a stochastic framework. One of the methodologies used to up‐scale the stochastic solute transport equation, from a point‐location scale to a grid scale, is the cumulant expansion method combined with the calculus for the time‐ordered exponential and the calculus for the Lie operator. When the point‐location scale transport equation is scaled up to the grid scale, using the cumulant expansion method, a new dispersion coefficient emerges in the dispersive term of the solute transport equation in addition to the molecular dispersion coefficient. This velocity driven dispersion is called ‘macrodispersion’. The macrodispersion coefficient is the integral function of the time‐ordered covariance of the random velocity field. The integral is calculated over a Lagrangian trajectory of the flow. The Lagrangian trajectory depends on the following: (i) the spatial origin of the particle; (ii) the time when the macrodispersion is calculated; and (iii) the mean velocity field along the trajectory itself. The Lagrangian trajectory is a recursive function of time because the location of the particle along the trajectory at a particular time depends on the location of the particle at the previous time. This recursive functional form of the Lagrangian trajectory makes the calculation of the macrodispersion coefficient difficult. Especially for the unsteady, spatially non‐stationary, non‐uniform flow field, the macrodispersion coefficient is a highly complex expression and, so far, calculated using numerical methods in the discrete domains. Here, an analytical method was introduced to calculate the macrodispersion coefficient in the discrete domain for the unsteady and steady, spatially non‐stationary flow cases accurately and efficiently. This study can fill the gap between the theory of the ensemble averaged solute transport model and its numerical implementations. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

6.
A new tracer experiment (referred to as MADE‐5) was conducted at the well‐known Macrodispersion Experiment (MADE) site to investigate the influence of small‐scale mass‐transfer and dispersion processes on well‐to‐well transport. The test was performed under dipole forced‐gradient flow conditions and concentrations were monitored in an extraction well and in two multilevel sampler (MLS) wells located at 6, 1.5, and 3.75 m from the source, respectively. The shape of the breakthrough curve (BTC) measured at the extraction well is strongly asymmetric showing a rapidly arriving peak and an extensive late‐time tail. The BTCs measured at seven different depths in the two MLSs are radically different from one another in terms of shape, arrival times, and magnitude of the concentration peaks. All of these characteristics indicate the presence of a complex network of preferential flow pathways controlling solute transport at the test site. Field‐experimental data were also used to evaluate two transport models: a stochastic advection‐dispersion model (ADM) based on conditional multivariate Gaussian realizations of the hydraulic conductivity field and a dual‐domain single‐rate (DDSR) mass‐transfer model based on a deterministic reconstruction of the aquifer heterogeneity. Unlike the stochastic ADM realizations, the DDSR accurately predicted the magnitude of the concentration peak and its arrival time (within a 1.5% error). For the multilevel BTCs between the injection and extraction wells, neither model reproduced the observed values, indicating that a high‐resolution characterization of the aquifer heterogeneity at the subdecimeter scale would be needed to fully capture 3D transport details.  相似文献   

7.
Numerical simulation of the acoustic wave equation is widely used to theoretically synthesize seismograms and constitutes the basis of reverse‐time migration. With finite‐difference methods, the discretization of temporal and spatial derivatives in wave equations introduces numerical grid dispersion. To reduce the grid dispersion effect, we propose to satisfy the dispersion relation for a number of uniformly distributed wavenumber points within a wavenumber range with the upper limit determined by the maximum source frequency, the grid spacing and the wave velocity. This new dispersion‐relationship‐preserving method relatively uniformly reduces the numerical dispersion over a large‐frequency range. Dispersion analysis and seismic numerical simulations demonstrate the effectiveness of the proposed method.  相似文献   

8.
1 INTRODUCTION In recent years, due to the increase in population and industrial developments, mankind has faced manyproblems associated with rivers, coastal waters and reservoirs. Some of these problems are flood control,water supply, power generation, and irrigation. In addition, making new hydraulic structures changesnatural conditions. Prediction of these changes is necessary for designing such constructions. For solutionof these problems usually an assessment of flow pattern, sedim…  相似文献   

9.
Small‐scale variations in surface moisture content were measured on a fine‐grained beach using a Delta‐T Theta probe. The resulting data set was used to examine the implications of small‐scale variability for estimating aeolian transport potential. Surface moisture measurements were collected on a 40 cm × 40 cm grid at 10 cm intervals, providing a total of 25 measurements for each grid data set. A total of 44 grid data sets were obtained from a representative set of beach sub‐environments. Measured moisture contents ranged from about 0% (dry) to 25% (saturated), by weight. The moisture content range within a grid data set was found to vary from less than 1% to almost 15%. The magnitude of within‐grid variability varied consistently with the mean moisture content of the grid sets, following an approximately normal distribution. Both very wet and very dry grid data sets exhibited little internal variability in moisture content, while intermediate moisture contents were associated with higher levels of variability. Thus, at intermediate moisture contents it was apparent that some portions of the beach surface could be dry enough to allow aeolian transport (i.e. moisture content is below the critical threshold), while adjacent portions are too wet for transport to occur. To examine the implications of this finding, cumulative distribution functions were calculated to model the relative proportions of beach surface area expected to be above or below specified threshold moisture levels (4%, 7%, and 14%). It was found that the implicit inclusion of small‐scale variability in surface moisture levels typically resulted in changes of less than 1% in the beach area available for transport, suggesting that this parameter can be ignored at larger spatial scales. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

10.
Distributed, continuous hydrologic models promote better understanding of hydrology and enable integrated hydrologic analyses by providing a more detailed picture of water transport processes across the varying landscape. However, such models are not widely used in routine modelling practices, due in part to the extensive data input requirements, computational demands, and complexity of routing algorithms. We developed a two‐dimensional continuous hydrologic model, HYSTAR, using a time‐area method within a grid‐based spatial data model with the goal of providing an alternative way to simulate spatiotemporally varied watershed‐scale hydrologic processes. The model calculates the direct runoff hydrograph by coupling a time‐area routing scheme with a dynamic rainfall excess sub‐model implemented here using a modified curve number method with an hourly time step, explicitly considering downstream ‘reinfiltration’ of routed surface runoff. Soil moisture content is determined at each time interval based on a water balance equation, and overland and channel runoff is routed on time‐area maps, representing spatial variation in hydraulic characteristics for each time interval in a storm event. Simulating runoff hydrographs does not depend on unit hydrograph theory or on solution of the Saint Venant equation, yet retains the simplicity of a unit hydrograph approach and the capability of explicitly simulating two‐dimensional flow routing. The model provided acceptable performance in predicting daily and monthly runoff for a 6‐year period for a watershed in Virginia (USA) using readily available geographic information about the watershed landscape. Spatial and temporal variability in simulated effective runoff depth and time area maps dynamically show the areas of the watershed contributing to the direct runoff hydrograph at the outlet over time, consistent with the variable source area overland flow generation mechanism. The model offers a way to simulate watershed processes and runoff hydrographs using the time‐area method, providing a simple, efficient, and sound framework that explicitly represents mechanisms of spatially and temporally varied hydrologic processes. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

11.
Considering heterogeneity in porous media pore size and connectivity is essential to predicting reactive solute transport across interfaces. However, exchange with less‐mobile porosity is rarely considered in surface water/groundwater recharge studies. Previous research indicates that a combination of pore‐fluid sampling and geoelectrical measurements can be used to quantify less‐mobile porosity exchange dynamics using the time‐varying relation between fluid and bulk electrical conductivity. For this study, we use macro‐scale (10 s of cm) advection–dispersion solute transport models linked with electrical conduction in COMSOL Multiphysics to explore less‐mobile porosity dynamics in two different types of observed sediment water interface porous media. Modeled sediment textures contrast from strongly layered streambed deposits to poorly sorted lakebed sands and cobbles. During simulated ionic tracer perturbations, a lag between fluid and bulk electrical conductivity, and the resultant hysteresis, is observed for all simulations indicating differential loading of pore spaces with tracer. Less‐mobile exchange parameters are determined graphically from these tracer time series data without the need for inverse numerical model simulation. In both sediment types, effective less‐mobile porosity exchange parameters are variable in response to changes in flow direction and fluid flux. These observed flow‐dependent effects directly impact local less‐mobile residence times and associated contact time for biogeochemical reaction. The simulations indicate that for the sediment textures explored here, less‐mobile porosity exchange is dominated by variable rates of advection through the domain, rather than diffusion of solute, for typical low‐to‐moderate rate (approximately 3–40 cm/day) hyporheic fluid fluxes. Overall, our model‐based results show that less‐mobile porosity may be expected in a range of natural hyporheic sediments and that changes in flowpath orientation and magnitude will impact less‐mobile exchange parameters. These temporal dynamics can be assessed with the geoelectrical experimental tracer method applied at laboratory and field scales.  相似文献   

12.
The numerical simulation of long‐term large‐scale (field to regional) variably saturated subsurface flow and transport remains a computational challenge, even with today's computing power. Therefore, it is appropriate to develop and use simplified models that focus on the main processes operating at the pertinent time and space scales, as long as the error introduced by the simpler model is small relative to the uncertainties associated with the spatial and temporal variation of boundary conditions and parameter values. This study investigates the effects of various model simplifications on the prediction of long‐term soil salinity and salt transport in irrigated soils. Average root‐zone salinity and cumulative annual drainage salt load were predicted for a 10‐year period using a one‐dimensional numerical flow and transport model (i.e. UNSATCHEM) that accounts for solute advection, dispersion and diffusion, and complex salt chemistry. The model uses daily values for rainfall, irrigation, and potential evapotranspiration rates. Model simulations consist of benchmark scenarios for different hypothetical cases that include shallow and deep water tables, different leaching fractions and soil gypsum content, and shallow groundwater salinity, with and without soil chemical reactions. These hypothetical benchmark simulations are compared with the results of various model simplifications that considered (i) annual average boundary conditions, (ii) coarser spatial discretization, and (iii) reducing the complexity of the salt‐soil reaction system. Based on the 10‐year simulation results, we conclude that salt transport modelling does not require daily boundary conditions, a fine spatial resolution, or complex salt chemistry. Instead, if the focus is on long‐term salinity, then a simplified modelling approach can be used, using annually averaged boundary conditions, a coarse spatial discretization, and inclusion of soil chemistry that only accounts for cation exchange and gypsum dissolution–precipitation. We also demonstrate that prediction errors due to these model simplifications may be small, when compared with effects of parameter uncertainty on model predictions. The proposed model simplifications lead to larger time steps and reduced computer simulation times by a factor of 1000. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

13.
More theoretical analysis is needed to investigate why a dual‐domain model often works better than the classical advection‐dispersion (AD) model in reproducing observed breakthrough curves for relatively homogeneous porous media, which do not contain distinct dual domains. Pore‐scale numerical experiments presented here reveal that hydrodynamics create preferential flow paths that occupy a small part of the domain but where most of the flow takes place. This creates a flow‐dependent configuration, where the total domain consists of a mobile and an immobile domain. Mass transfer limitations may result in nonequilibrium, or significant differences in concentration, between the apparent mobile and immobile zones. When the advection timescale is smaller than the diffusion timescale, the dual‐domain mass transfer (DDMT) model better captures the tailing in the breakthrough curve. Moreover, the model parameters (mobile porosity, mean solute velocity, dispersivity, and mass transfer coefficient) demonstrate nonlinear dependency on mean fluid velocity. The studied case also shows that when the Peclet number, Pe, is large enough, the mobile porosity approaches a constant, and the mass transfer coefficient can be approximated as proportional to mean fluid velocity. Based on detailed analysis at the pore scale, this paper provides a physical explanation why these model parameters vary in certain ways with Pe. In addition, to improve prediction in practical applications, we recommend conducting experiments for parameterization of the DDMT model at a velocity close to that of the relevant field sites, or over a range of velocities that may allow a better parameterization.  相似文献   

14.
Fluid flow exerts a critical impact on the convection of thermal energy in geological media, whereas heat transport in turn affects fluid properties, including fluid dynamic viscosity and density. The interplay of flow and heat transport also affects solute transport. To unravel these complex coupled flow, heat, and solute transport processes, here, we present a theory for the idealized scale‐dependent Poiseuille flow model considering a constant temperature gradient (?T) along a single fracture, where fluid dynamic viscosity connects with temperature via an exponential function. The idealized scale‐dependent model is validated based on the solutions from direct numerical simulations. We find that the hydraulic conductivity (K) of the Poiseuille flow either increases or decreases with scales depending on ?T > 0°C/m or ?T < 0°C/m, respectively. Indeed, the degree of changes in K depends on the magnitude of ?T and fracture length. The scale‐dependent model provides an alternative explanation for the well‐known scale‐dependent transport problem, for example, the dispersion coefficient increases with travel distance when ?T > 0°C/m according to the Taylor dispersion theory, because K (or equivalently flux through fractures) scales with fracture length. The proposed theory unravels intertwined interactions between flow and transport processes, which might shed light on understanding many practical geophysical problems, for example, geothermal energy exploration.  相似文献   

15.
It has been known for many years that dispersivities increase with solute displacement distance in a subsurface. The increase of dispersivities with solute travel distance results from significant variation in hydraulic properties of porous media and was identified in the literature as scale‐dependent dispersion. In this study, Laplace‐transformed analytical solutions to advection‐dispersion equations in cylindrical coordinates are derived for interpreting a divergent flow tracer test with a constant dispersivity and with a linear scale‐dependent dispersivity. Breakthrough curves obtained using the scale‐dependent dispersivity model are compared to breakthrough curves obtained from the constant dispersivity model to illustrate the salient features of scale‐dependent dispersion in a divergent flow tracer test. The analytical results reveal that the breakthrough curves at the specific location for the constant dispersivity model can produce the same shape as those from the scale‐dependent dispersivity model. This correspondence in curve shape between these two models occurs when the local dispersivity at an observation well in the scale‐dependent dispersivity model is 1·3 times greater than the constant dispersivity in the constant dispersivity model. To confirm this finding, a set of previously reported data is interpreted using both the scale‐dependent dispersivity model and the constant dispersivity model to distinguish the differences in scale dependence of estimated dispersivity from these two models. The analytical result reveals that previously reported dispersivity/distance ratios from the constant dispersivity model should be revised by multiplying these values by a factor of 1·3 for the scale‐dependent dispersion model if the dispersion process is more accurately characterized by scale‐dependent dispersion. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

16.
A comprehensive understanding of the dynamics of erosion and sedimentation in reservoirs under different management conditions is required to anticipate sedimentation issues and implement effective sediment management strategies. This paper describes a unique approach combining fluvial geomorphology tools and morphodynamic modeling for analyzing the sediment dynamics of an elongated hydropower reservoir subjected to management operations: the Génissiat Reservoir on the Rhône River. Functional sub‐reaches representative of the reservoir morphodynamics were delineated by adapting natural river segmentation methods to elongated reservoirs. The segmentation revealed the link between the spatial and temporal reservoir changes and the variability of longitudinal flow conditions during reservoir management operations. An innovative modeling strategy, incorporating the reservoir segmentation into two sediment transport codes, was implemented to simulate the dynamics of erosion and sedimentation at the reach scale during historic events. One code used a bedload approach, based on the Exner equation with a transport capacity formula, and the other used a suspended load approach based on the advection–dispersion equation. This strategy provided a fair quantification of the dynamics of erosion and sedimentation at the reach scale during different management operations. This study showed that the reservoir morphodynamics is controlled by bedload transport in upper reaches, graded suspended load transport of sand in middle reaches and suspended load transport of fine sediments in lower reaches. Eventually, it allowed a better understanding of the impact of dam management on sediment dynamics. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

17.
As a fundamental unit of the landscape, hillslopes are studied for their retention and release of water and nutrients across a wide range of ecosystems. The understanding of these near‐surface processes is relevant to issues of runoff generation, groundwater–surface water interactions, catchment export of nutrients, dissolved organic carbon, contaminants (e.g. mercury) and ultimately surface water health. We develop a 3‐D physics‐based representation of the Panola Mountain Research Watershed experimental hillslope using the TOUGH2 sub‐surface flow and transport simulator. A recent investigation of sub‐surface flow within this experimental hillslope has generated important knowledge of threshold rainfall‐runoff response and its relation to patterns of transient water table development. This work has identified components of the 3‐D sub‐surface, such as bedrock topography, that contribute to changing connectivity in saturated zones and the generation of sub‐surface stormflow. Here, we test the ability of a 3‐D hillslope model (both calibrated and uncalibrated) to simulate forested hillslope rainfall‐runoff response and internal transient sub‐surface stormflow dynamics. We also provide a transparent illustration of physics‐based model development, issues of parameterization, examples of model rejection and usefulness of data types (e.g. runoff, mean soil moisture and transient water table depth) to the model enterprise. Our simulations show the inability of an uncalibrated model based on laboratory and field characterization of soil properties and topography to successfully simulate the integrated hydrological response or the distributed water table within the soil profile. Although not an uncommon result, the failure of the field‐based characterized model to represent system behaviour is an important challenge that continues to vex scientists at many scales. We focus our attention particularly on examining the influence of bedrock permeability, soil anisotropy and drainable porosity on the development of patterns of transient groundwater and sub‐surface flow. Internal dynamics of transient water table development prove to be essential in determining appropriate model parameterization. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

18.
The MT3DMS groundwater solute transport model was modified to simulate solute transport in the unsaturated zone by incorporating the unsaturated‐zone flow (UZF1) package developed for MODFLOW. The modified MT3DMS code uses a volume‐averaged approach in which Lagrangian‐based UZF1 fluid fluxes and storage changes are mapped onto a fixed grid. Referred to as UZF‐MT3DMS, the linked model was tested against published benchmarks solved analytically as well as against other published codes, most frequently the U.S. Geological Survey's Variably‐Saturated Two‐Dimensional Flow and Transport Model. Results from a suite of test cases demonstrate that the modified code accurately simulates solute advection, dispersion, and reaction in the unsaturated zone. Two‐ and three‐dimensional simulations also were investigated to ensure unsaturated‐saturated zone interaction was simulated correctly. Because the UZF1 solution is analytical, large‐scale flow and transport investigations can be performed free from the computational and data burdens required by numerical solutions to Richards' equation. Results demonstrate that significant simulation runtime savings can be achieved with UZF‐MT3DMS, an important development when hundreds or thousands of model runs are required during parameter estimation and uncertainty analysis. Three‐dimensional variably saturated flow and transport simulations revealed UZF‐MT3DMS to have runtimes that are less than one tenth of the time required by models that rely on Richards' equation. Given its accuracy and efficiency, and the wide‐spread use of both MODFLOW and MT3DMS, the added capability of unsaturated‐zone transport in this familiar modeling framework stands to benefit a broad user‐ship.  相似文献   

19.
This study assesses hydrodynamic and morphodynamic model sensitivity and functionality in a curved channel. The sensitivity of a depth‐averaged model to user‐defined parameters (grain size, roughness, transverse bed slope effect, transport relations and secondary flow) is tested. According to the sensitivity analysis, grain size, transverse bed slope effect and sediment transport relations are critical to simulated meander bend morphodynamics. The parametrization of grain size has the most remarkable effect: field‐based grain size parametrization is necessary in a successful morphodynamic reconstruction of a meander bend. The roughness parametrization method affects the distribution of flow velocities and therefore also morphodynamics. The combined effect of various parameters needs further research. Two‐dimensional (2D) and three‐dimensional (3D) reconstructions of a natural meander bend during a flood event are assessed against field measurements of acoustic Doppler current profiler and multi‐temporal mobile laser scanning data. The depth‐averaged velocities are simulated satisfactorily (differences from acoustic Doppler current profiler velocities 5–14%) in both 2D and 3D simulations, but the advantage of the 3D hydrodynamic model is unquestionable because of its ability to model vertical and near‐bed flows. The measured and modelled near‐bed flow, however, differed notably from each other's, the reason of which was left open for future research. It was challenging to model flow direction beyond the apex. The 3D flow features, which also affected the distribution of the bed shear stress, seem not to have much effect on the predicted morphodynamics: the 2D and 3D morphodynamic reconstructions over the point bar resembled each other closely. Although common features between the modelled and measured morphological changes were also found, some specific changes that occurred were not evident in the simulation results. Our results show that short‐term, sub‐bend scale morphodynamic processes of a natural meander bend are challenging to model, which implies that they are affected by factors that have been neglected in the simulations. The modelling of short‐term morphodynamics in natural curved channel is a challenge that requires further study. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

20.
In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field.  相似文献   

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

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