Efficient DEM‐based overland flow routing using integrated recursive algorithms          下载免费PDF全文

Pin‐Chun Huang  Kwan Tun Lee 《水文研究》2017,31(5):1007-1017

Two primary concerns in performing watershed overland flow routing are the numerical instability and computational efficiency. The stability of executing an explicit scheme has to be maintained by observing the Courant–Friedrich–Lewy criterion, which is adopted to confirm that the numerical marching speed is larger than the wave celerity. Moreover, there is another criterion of time step devised in previous studies to avoid back‐and‐forth refluxing between adjacent grids. The situation of refluxing usually occurs on flat regions. In light of this, the selection of a small time increment to honor both restrictions simultaneously is believed to decrease the computational efficiency in performing overland flow routing. This study aims at creating a robust algorithm to relax both restrictions. The proposed algorithm was first implemented on a one‐dimensional overland plane to evaluate the accuracy of the numerical result by comparing it with an analytical solution. Then, the algorithm was further applied to a watershed for 2D runoff simulations. The results show that the proposed integrated algorithm can provide an accurate runoff simulation and achieve satisfactory performance in terms of computational speed.  相似文献   

Application of the numerical manifold method for stress wave propagation across rock masses     

G. F. Zhao  X. B. Zhao  J. B. Zhu 《国际地质力学数值与分析法杂志》2014,38(1):92-110

In this paper, the numerical manifold method (NMM) is extended to study wave propagation across rock masses. First, improvements to the system equations, contact treatment, and boundary conditions of the NMM are performed, where new system equations are derived based on the Newmark assumption of the space–time relationship, the edge‐to‐edge contact treatment is further developed for the NMM to handle stress wave propagation across discontinuities, and the viscous non‐reflection boundary condition is derived based on the energy minimisation principle. After the modification, numerical comparisons between the original and improved NMM are presented. The results show that the original system equations result in artificial numerical damping, which can be overcome by the Newmark system equations. Meanwhile, the original contact scheme suffers some calculation problems when modelling stress wave propagation across a discontinuity, which can be solved by the proposed edge‐to‐edge contact scheme. Subsequently, the influence of the mesh size and time step on the improved NMM for stress wave propagation is studied. Finally, 2D wave propagation is modelled, and the model's results are in good agreement with the analytical solution. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

全球百年海表面温度年际和年代际变化特征分析   总被引：1，自引：0，他引：1  

张锦婷  刘秦玉  武术 《海洋学报》2010,32(4):24-31

利用100 a(1903—2002年)HadiSST的逐月资料,将全球海表面温度异常(SSTA)作为整体进行经验正交函数分解(EOF分解),提取了控制各大洋SSTA的主导模态和各大洋之间的联合模态,分析了相应的空间分布和时间序列。研究结果表明:SSTA变化最剧烈的海区是赤道中东太平洋、西北大西洋湾流海区和北太平洋黑潮延伸体海区。热带太平洋厄尔尼诺和南方涛动是主导模态并具有2~7 a周期的年际变化;SSTA变化第二主模态和第三主模态都是以大约70 a为周期的年代际变化为主的跨大洋联合模态。第二主模态的空间分布主要表现为中纬度北太平洋和北大西洋反位相、热带太平洋与大西洋反位相的太平洋-大西洋双偶极子型分布。SSTA变化的第三模态主要呈现南北半球海洋反位相的特征,为北太平洋-北大西洋-南大洋联合模态。第四模态基本上是反映各个不同海域特有的局地海洋-大气相互作用模态,该模态的时间序列具有1~4 a周期的年际变化和约9 a周期的年代际变化。  相似文献   

Analytical method of radiation by a wave energy device with dual rectangular floating bodies     

WANG Wensheng  YOU Yage  ZHANG Yunqiu  WU Bijun  YE Yin 《海洋学报(英文版)》2010,29(3):106-115

The system with one floating rectangular body on the free surface and one submerged rectangular body has been applied to a wave energy conversion device in water of finite depth. The radiation problem by this device on a plane incident wave is solved by the use of an eigenfunction expansion method, and a new analytical expression for the radiation velocity potential is obtained. The wave excitation force is calculated via the known incident wave potential and the radiation potential with a theorem of Haskind employed. To verify the correctness of this method, an example is computed respectively through the bound element method and analytical method. Results show that two numerical methods. are in good agreement, which shows that the present method is applicable. In addition, the trends of hydrodynamic coefficients and wave force are analyzed under different conditions by use of the present analytical method.  相似文献   

Predicting the dependency of a degree of saturation on void ratio and suction using effective stress principle for unsaturated soils     

David Mašín 《国际地质力学数值与分析法杂志》2010,34(1):73-90

The paper presents an approach to predicting variation of a degree of saturation in unsaturated soils with void ratio and suction. The approach is based on the effective stress principle for unsaturated soils and several underlying assumptions. It focuses on the main drying and wetting processes and does not incorporate the effects of hydraulic hysteresis. It leads to the dependency of water retention curve (WRC) on void ratio, which does not require any material parameters apart from the parameters specifying WRC for the reference void ratio. Its validity is demonstrated by comparing predictions with the experimental data on four different soils taken over from the literature. Good correlation between the measured and predicted behaviour indirectly supports applicability of the effective stress principle for unsaturated soils. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

Cracking risk of partially saturated porous media—Part I: Microporoelasticity model     

Bernhard Pichler  Luc Dormieux 《国际地质力学数值与分析法杂志》2010,34(2):135-157

Drying of deformable porous media results in their shrinkage, and it may cause cracking provided that shrinkage deformations are hindered by kinematic constraints. This is the motivation to develop a thermodynamics‐based microporoelasticity model for the assessment of cracking risk in partially saturated porous geomaterials. The study refers to 3D representative volume elements of porous media, including a two‐scale double‐porosity material with a pore network comprising (at the mesoscale) 3D mesocracks in the form of oblate spheroids, and (at the microscale) spherical micropores of different sizes. Surface tensions prevailing in all interfaces between solid, liquid, and gaseous matters are taken into account. To establish a thermodynamics‐based crack propagation criterion for a two‐scale double‐porosity material, the potential energy of the solid is derived, accounting—in particular—for mesocrack geometry changes (main original contribution) and for effective micropore pressures, which depend (due to surface tensions) on the pore radius. Differentiating the potential energy with respect to crack density parameter yields the thermodynamical driving force for crack propagation, which is shown to be governed by an effective macrostrain. It is found that drying‐related stresses in partially saturated mesocracks reduce the cracking risk. The drying‐related effective underpressures in spherical micropores, in turn, result in a tensile eigenstress of the matrix in which the mesocracks are embedded. This way, micropores increase the mesocracking risk. Model application to the assessment of cracking risk during drying of argillite is the topic of the companion paper (Part II). Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

Effect of bed surface roughness on longitudinal dispersion in artificial open channels     

Marcus Schulz  Jan Priegnitz  Jörg Klasmeier  Stefan Heller  Stefan Meinecke  Michael Feibicke 《水文研究》2012,26(2):272-280

In several empirical and modelling studies on river hydraulics, dispersion was negatively correlated to surface roughness. In this study, it was aimed to investigate the influence of surface roughness on longitudinal dispersion under controlled conditions. In artificial flow channels with a length of 104 m, tracer experiments with variations in channel bed material were performed. By use of measured tracer breakthrough curves, average flow velocity, mean longitudinal dispersion, and mean longitudinal dispersivity were calculated. Longitudinal dispersion coefficients ranged from 0·018 m² s^?1 in channels with smooth bed surface up to 0·209 m² s^?1 in channels with coarse gravel as bed material. Longitudinal dispersion was linearly related to mean flow velocity. Accordingly, longitudinal dispersivities ranged between 0·152 ± 0·017 m in channels with smooth bed surface and 0·584 ± 0·015 m in identical channels with a coarse gravel substrate. Grain size and surface roughness of the channel bed were found to correlate positively to longitudinal dispersion. This finding contradicts several existing relations between surface roughness and dispersion. Future studies should include further variation in surface roughness to derive a better‐founded empirical equation forecasting longitudinal dispersion from surface roughness. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

Difference-frequency ultrasound generation from microbubbles under dual-frequency excitation     

马青玉  邱媛媛  黄蓓  章东  龚秀芬 《海洋学报(英文版)》2010,29(9)

The difference-frequency (DF) ultrasound generated by using parametric effect promises to improve detection depth owing to its low attenuation, which is beneficial for deep tissue imaging. With ultrasound contrast agents infusion, the harmonic components scattered from the microbubbles, including DF, can be generated due to the nonlinear vibration. A theoretical study on the DF generation from microbubbles under the dual-frequency excitation is proposed in formula based on the solution of the RPNNP equation. The optimisation of the DF generation is discussed associated with the applied acoustic pressure, frequency, and the microbubble size. Experiments are performed to validate the theoretical predictions by using a dual-frequency signal to excite microbubbles. Both the numerical and experimental results demonstrate that the optimised DF ultrasound can be achieved as the difference frequency is close to the resonance frequency of the microbubble and improve the contrast-to-tissue ratio in imaging.  相似文献   

The behavior of magneto-acoustic-gravity waves near the cusp resonance in a lossless,compressible, isothermal,stratified, electrically conducting and uniformly magnetized atmosphere. II: The valve effect     

Leon P. J. Kamp 《地球物理与天体物理流体动力学》2013,107(4):257-271

Abstract

In the present paper the behavior of (internal) magneto-acoustic-gravity waves near the cusp resonance in a lossless, compressible, isothermal, stratified, electrically conducting atmosphere that is permeated by a uniform, nearly horizontal magnetic field is re-addressed (Kamp, 1989). The previously analyzed linear conversion of long acoustic-gravity waves into short magneto-acoustic waves that carry off the energy from the resonance region along the magnetic field, is re-analyzed with boundary layer techniques that are based on the smallness of the vertical component of the magnetic field. More specifically the existence of the so-called valve effect for the generated magneto-acoustic mode near the critical level is explicitly demonstrated and shown to be governed by two rivalling effects.  相似文献   

Monitoring over‐winter soil water dynamics in a freezing and snow‐covered environment using a thermally insulated tensiometer     

Yukiyoshi Iwata  Tomoyoshi Hirota 《水文研究》2005,19(15):3013-3019

Water potential below a frozen soil layer was continuously monitored over an entire winter period (using thermally insulated tensiometers sheltered in a heated chamber) along with other soil, snow and atmospheric variables. In early winter, the freezing front advanced under a thin snow cover, inducing upward soil water flow in the underlying unfrozen soil. The freezing front started to retreat when the snow cover became thick enough to insulate the soil, resulting in the reversal of the flow direction in the unfrozen zone. These data provide a clear illustration of soil water dynamics, which have rarely been monitored with a tensiometer. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献