A spatially and temporally adaptive solution of Richards’ equation     

Cass T. Miller  Chandra Abhishek  Matthew W. Farthing 《Advances in water resources》2006

Efficient, robust simulation of groundwater flow in the unsaturated zone remains computationally expensive, especially for problems characterized by sharp fronts in both space and time. Standard approaches that employ uniform spatial and temporal discretizations for the numerical solution of these problems lead to inefficient and expensive simulations. In this work, we solve Richards’ equation using adaptive methods in both space and time. Spatial adaption is based upon a coarse grid solve and a gradient error indicator using a fixed-order approximation. Temporal adaption is accomplished using variable order, variable step size approximations based upon the backward difference formulas up to fifth order. Since the advantages of similar adaptive methods in time are now established, we evaluate our method by comparison with a uniform spatial discretization that is adaptive in time for four different one-dimensional test problems. The numerical results demonstrate that the proposed method provides a robust and efficient alternative to standard approaches for simulating variably saturated flow in one spatial dimension.  相似文献   

Integrated surface–groundwater flow modeling: A free-surface overland flow boundary condition in a parallel groundwater flow model     

Stefan J. Kollet  Reed M. Maxwell 《Advances in water resources》2006

Interactions between surface and groundwater are a key component of the hydrologic budget on the watershed scale. Models that honor these interactions are commonly based on the conductance concept that presumes a distinct interface at the land surface, separating the surface from the subsurface domain. These types of models link the subsurface and surface domains via an exchange flux that depends upon the magnitude and direction of the hydraulic gradient across the interface and a proportionality constant (a measure of the hydraulic connectivity). Because experimental evidence of such a distinct interface is often lacking in field systems, there is a need for a more general coupled modeling approach.  相似文献   

Influence of uncertain boundary conditions and model structure on flood inundation predictions   总被引：2，自引：0，他引：2  

Florian Pappenberger  Patrick Matgen  Keith J. Beven  Jean-Baptiste Henry  Laurent Pfister  Paul Fraipont de 《Advances in water resources》2006

In this study, the GLUE methodology is applied to establish the sensitivity of flood inundation predictions to uncertainty of the upstream boundary condition and bridges within the modelled region. An understanding of such uncertainties is essential to improve flood forecasting and floodplain mapping. The model has been evaluated on a large data set. This paper shows uncertainty of the upstream boundary can have significant impact on the model results, exceeding the importance of model parameter uncertainty in some areas. However, this depends on the hydraulic conditions in the reach e.g. internal boundary conditions and, for example, the amount of backwater within the modelled region. The type of bridge implementation can have local effects, which is strongly influenced by the bridge geometry (in this case the area of the culvert). However, the type of bridge will not merely influence the model performance within the region of the structure, but also other evaluation criteria such as the travel time. This also highlights the difficulties in establishing which parameters have to be more closely examined in order to achieve better fits. In this study no parameter set or model implementation that fulfils all evaluation criteria could be established. We propose four different approaches to this problem: closer investigation of anomalies; introduction of local parameters; increasing the size of acceptable error bounds; and resorting to local model evaluation. Moreover, we show that it can be advantageous to decouple the classification into behavioural and non-behavioural model data/parameter sets from the calculation of uncertainty bounds.  相似文献   

Pore scale modeling of immiscible and miscible fluid flows using smoothed particle hydrodynamics     

Alexandre M. Tartakovsky  Paul Meakin 《Advances in water resources》2006

A numerical model based on smoothed particle hydrodynamics (SPH) was developed and used to simulate immiscible and miscible fluid flows in porous media and to study effects of pore scale heterogeneity and anisotropy on such flows.  相似文献   

A bed‐load transport model for rough turbulent open‐channel flows on plane beds     

Athol D. Abrahams  Peng Gao 《地球表面变化过程与地形》2006,31(7):910-928

Data from flume studies are used to develop a model for predicting bed‐load transport rates in rough turbulent two‐dimensional open‐channel flows moving well sorted non‐cohesive sediments over plane mobile beds. The object is not to predict transport rates in natural channel flows but rather to provide a standard against which measured bed‐load transport rates influenced by factors such as bed forms, bed armouring, or limited sediment availability may be compared in order to assess the impact of these factors on bed‐load transport rates. The model is based on a revised version of Bagnold's basic energy equation i_bs_b = e_bω, where i_b is the immersed bed‐load transport rate, ω is flow power per unit area, e_b is the efficiency coefficient, and s_b is the stress coefficient defined as the ratio of the tangential bed shear stress caused by grain collisions and fluid drag to the immersed weight of the bed load. Expressions are developed for s_b and e_b in terms of G, a normalized measure of sediment transport stage, and these expressions are substituted into the revised energy equation to obtain the bed‐load transport equation i_b = ω G ^3·4. This equation applies regardless of the mode of bed‐load transport (i.e. saltation or sheet flow) and reduces to i_b = ω where G approaches 1 in the sheet‐flow regime. That i_b = ω does not mean that all the available power is dissipated in transporting the bed load. Rather, it reflects the fact that i_b is a transport rate that must be multiplied by s_b to become a work rate before it can be compared with ω. It follows that the proportion of ω that is dissipated in the transport of bed load is i_bs_b/ω, which is approximately 0·6 when i_b = ω. It is suggested that this remarkably high transport efficiency is achieved in sheet flow (1) because the ratio of grain‐to‐grain to grain‐to‐bed collisions increases with bed shear stress, and (2) because on average much more momentum is lost in a grain‐to‐bed collision than in a grain‐to‐grain one. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

雅鲁藏布江缝合带内构造岩片划分及形成——以仲巴-桑桑段为例   总被引：1，自引：2，他引：1  

吴新国  吕继东 《大地构造与成矿学》2006,30(3):320-325

构造岩片作为造山带中非史密斯地层地区的填图单位有其重要的地质意义,本文将雅鲁藏布江缝合带分为6大构造岩片,其中蛇绿混杂岩岩片又分4个次级岩片,进一步阐述了缝合带内不同单元具有不同岩性和不同构造型式,是一条复杂的构造带,并发现俯冲板块一侧的不同构造单元具有自南而北,由脆性变形到低温高压变质、脆-韧性变形,超基性岩的深层次韧性变形的分布规律,大多数岩片的分布是通过多重逆冲推覆型或走滑型的断层将其从地壳深部挤出定位的。雅鲁藏布江缝合带内的构造岩片特征,反映了冈瓦那大陆北缘自三叠纪末期到第四纪经历了裂陷、扩张、多岛洋、闭合、俯冲、碰撞、褶皱、推覆造山的全过程。  相似文献   

裂隙岩质边坡地下水渗流数值模拟     

刘建军  卢正  江峰 《地下水》2006,28(1):26-28

本文从反倾向岩质高边坡的渗流特征和渗流普遍规律出发,根据渗流力学知识,建立数学模型,通过渗流数值模拟计算分析了渗透压力、水头变化规律以及对岩质边坡稳定性的影响.  相似文献   

基于GIS的滑坡、泥石流灾害危险性区划关键问题研究   总被引：20，自引：3，他引：20  

丛威青  潘懋  李铁锋  吴自兴  吕广宪 《地学前缘》2006,13(1):185-190

随着GIS技术的引入,滑坡、泥石流灾害危险性区划的效率和准确性得以大大提高。依据工程地质类比原则,在灾害学理论指导下,结合专家打分、层次分析、人工神经网络、信息量、Logistic回归、统计量等模型方法,以MAPGIS软件为平台,利用C++语言开发了滑坡、泥石流灾害危险性区划评价分析系统;并重点探讨了GIS支持下的滑坡、泥石流灾害危险性区划过程中的因子分析、模型选取、模型复合、单元划分、系统集成、结果评价等关键问题,建立了一整套基于GIS的滑坡、泥石流灾害评价方法体系。应用该系统对长江三峡库区和辽宁省鞍山市分别开展了滑坡、泥石流灾害危险性区划研究,取得了较好的效果。  相似文献   

姑婆山复式岩体的基本特征及其与成矿作用的关系   总被引：5，自引：0，他引：5  

王建辉  蒋柏昌  孔祥伟  唐奎保 《矿产与地质》2006,20(1):47-53

姑婆山复式岩体由晚侏罗世里松岩体、姑婆山岩体及早白垩世新路岩体组成。复式岩体为壳幔混合源成因。据岩石地球化学判断和分析,姑婆山岩体和新路岩体为含锡岩体;整个姑婆山复式岩体为含钨岩体,里松岩体和姑婆山岩体为含稀有、稀土岩体。上述矿产的形成与岩石中挥发份F、B及放射性U、T h等元素的高含量有关。  相似文献   

Mantle heat flow and thermal structure of the northern block of Southern Granulite Terrain, India   总被引：1，自引：0，他引：1  

Ajay Manglik   《Journal of Geodynamics》2006,41(5):510-519

Continental shield regions are normally characterized by low-to-moderate mantle heat flow. Archaean Dharwar craton of the Indian continental shield also follows the similar global pattern. However, some recent studies have inferred significantly higher mantle heat flow for the Proterozoic northern block of Southern Granulite Terrain (SGT) in the immediate vicinity of the Dharwar craton by assuming that the radiogenic elements depleted exposed granulites constitute the 45-km-thick crust. In this study, we use four-layered model of the crustal structure revealed by integrated geophysical studies along a geo-transect in this region to estimate the mantle heat flow. The results indicate that: (i) the mantle heat flow of the northern block of SGT is 17 ± 2 mW/m², supporting the global pattern, and (ii) the lateral variability of 10–12 mW/m² in the surface heat flow within the block is of crustal origin. In terms of temperature, the Moho beneath the eastern Salem–Namakkal region appears to be at 80–100 °C higher temperature than that beneath the western Avinashi region.  相似文献