Well-balanced high-order centered schemes on unstructured meshes for shallow water equations with fixed and mobile bed     

Alberto Canestrelli  Michael Dumbser  Annunziato Siviglia  Eleuterio F. Toro 《Advances in water resources》2010

In this paper, we study the numerical approximation of the two-dimensional morphodynamic model governed by the shallow water equations and bed-load transport following a coupled solution strategy. The resulting system of governing equations contains non-conservative products and it is solved simultaneously within each time step. The numerical solution is obtained using a new high-order accurate centered scheme of the finite volume type on unstructured meshes, which is an extension of the one-dimensional PRICE-C scheme recently proposed in Canestrelli et al. (2009) [5]. The resulting first-order accurate centered method is then extended to high order of accuracy in space via a high order WENO reconstruction technique and in time via a local continuous space–time Galerkin predictor method. The scheme is applied to the shallow water equations and the well-balanced properties of the method are investigated. Finally, we apply the new scheme to different test cases with both fixed and movable bed. An attractive future of the proposed method is that it is particularly suitable for engineering applications since it allows practitioners to adopt the most suitable sediment transport formula which better fits the field data.  相似文献   

Numerical modelling of coupled surface and subsurface flow systems     

Prince Chidyagwai  Beatrice Rivière 《Advances in water resources》2010

This paper presents and compares several numerical solutions of the coupled system of Navier–Stokes and Darcy equations. The schemes are based on combinations of the finite element method and the discontinuous Galerkin method. Accuracy and robustness of the methods are investigated for heterogeneous porous media. The importance of local mass conservation for filtration problems is also discussed.  相似文献   

A SANISAND model with anisotropic elasticity     

Ali Lashkari 《Soil Dynamics and Earthquake Engineering》2010

Numerous experimental studies indicate that as a result of shear stress, the elastic behavior of granular media becomes both non-linear and anisotropic. This paper presents a simple constitutive model for sands with respect to anisotropic elasticity. To this aim, using the concept of second order fabric tensor, a simplified elasticity theory is presented which is capable of considering the effect of induced anisotropy on the elastic response. SANISAND is the name used for a family of simple anisotropic sand models developed in the framework of critical state soil mechanics and bounding surface plasticity. An existing SANISAND model is modified in order to include the proposed anisotropic elasticity. The modified model simulations are compared with those obtained from the other members of this family. It is shown that considering anisotropic elasticity can take part in explanation of drastic loss of mean principal stress when sand is subjected to reverse loading in dilative branch of behavior and as a result, improve the liquefaction simulations.  相似文献   

基于内变量和张量函数表示定理的本构方程   总被引：2，自引：2，他引：0  

陈明祥  汪碧飞 《岩土力学》2010,31(2):397-402

针对各向同性材料,基于张量函数表示定理,建立了本构关系的张量不变性表示,其中,3个不可约基张量取决于应力的0～2次幂,且相互正交,3个系数由塑性应变增量和应力的不变量表示。基于塑性应变增量的不变量定义内变量,本构关系归结为确定内变量的演化。使用张量函数表示定理,给出了内变量演化方程的一般表达式,它取决于应力不变量的增量,因而与主轴旋转无关。讨论了如何根据试验资料和引入适当的假定,确定具体的演化方程。通过与塑性势理论和多重屈服面理论进行比较,表明所建模型是这些理论的最一般表示,且简捷直观、使用方便。  相似文献   

饱和地基上条形弹性基础的摇摆振动     

马晓华  蔡袁强  徐长节 《岩土力学》2010,31(7):2164-2172

基于Biot动力方程,研究了饱和均质弹性半空间上弹性条形基础的摇摆振动问题。通过Fourier积分变换求解了饱和土的动力控制方程,然后结合基础底部为混合边界的条件得到了弹性条形基础的摇摆振动对偶积分方程,利用正交多项式将对偶积分方程转化为求解一组线性代数方程组,同时利用复合Simpson法则,得到了动力柔度系数的表达式,通过算例得出了不同参数时地基动力柔度系数随无量纲频率的关系曲线。  相似文献   

大尺度圆柱墩群周围的波流场的数值模拟   总被引：5，自引：0，他引：5  

李玉成  刘德良  陈兵  李林普 《海洋通报》2005,24(2):1-12

本文对波流共同作用下大尺度圆柱墩群周围的波流场进行了数值研究。利用波浪弥散关系的迭代计算求得波向与流向的夹角以及波浪的相对频率。流场通过求解浅水环流方程得到,波浪场通过求解含流的缓坡方程得到,通过二者的迭代计算得到大尺度圆柱墩群周围的波流场的耦合解。用有限元法建立了数值模型,并将本文的计算数据与试验数据以及其他学者计算数据进行了比较,结果较为合理。  相似文献   

泥石流的二维数学模型   总被引：5，自引：2，他引：3  

王纯祥  白世伟  江崎哲郎  三谷泰浩 《岩土力学》2007,28(6):1237-1241

泥石流是在重力作用下,由砂粒石块和水等组成的固液混合物,是一种发生于山区的复杂的地质灾害现象。泥石流主要是由暴雨诱发引起的,它沿着复杂的三维地形高速流动,具有流体流动的特性。为了模拟泥石流的运动规律,预测降雨诱发的泥石流的到达距离和泛滥范围,减少和避免泥石流引起的灾害,把泥石和雨水组成的固液混合物假定为遵循均匀、连续、不可压缩的、非定常的牛顿流体运动规律。基于质量守恒方程和Naiver-stokes方程,采用深度积分方法,推导出了一个模拟泥石流运动的二维数学模型。所有方程式可用有限差分法来求解。结合GIS,该模型可用于预测泥石流的流动距离和泛滥范围,以及泛滥范围内的危险房屋和路段,也可以用于泥石流灾害的风险性分析。  相似文献   

“98.7”暴雨β中尺度低涡生成发展结构演变:双向四重嵌套网格模拟   总被引：18，自引：0，他引：18  

程麟生  冯伍虎 《气象学报》2003,61(4):385-395

用非静力中尺度模式MM5.V3全物理湿过程和双向四重嵌套网格技术及四维资料同化分析张弛方法,对1998年7月20～22日(简称"98.7")鄂东武汉沿江地区特大暴雨进行了36 h数值模拟,结果大大改进了原二重和三重嵌套网格的模拟,特别是水平分辨率为2 km的嵌套子域D04能更正确地模拟出暴雨的落区及雨强中心,进而揭示了强暴雨β中尺度系统发生和发展的结构及演变. 其主要结果包括(1)β中尺度切变线在鄂东沿江低空强烈发展及辐合中心的出现与其β中尺度低涡的形成和发展直接关联.(2)β中尺度切变线强烈发展的垂直结构强辐合层和强辐散层复式迭置并与强上升运动耦合发展;强涡度层和强位涡度层与强辐合层互伴发展;低空湿位温中心与中空饱和水汽带共存.(3)β中尺度低涡生成的垂直结构散度和上升运动均呈双支柱状发展;涡度和位涡度均呈单支柱状发展;高湿能柱呈双支耦合发展,水汽通道呈阶梯斜升状.(4)β中尺度低涡发展的垂直结构V字型散度柱和上升运动柱互耦发展;涡度和位涡度呈双支柱状;双支高湿能柱强烈发展,阶梯斜升水汽通道变宽增厚.至此低涡发展达最强,其结构具有典型性.模拟结果还指出,发展时空分辨率更高的多重嵌套网格模拟技术和应用四维资料同化方法,将有助于更细致的了解β中尺度强对流系统发生和发展的结构及其演变,并能进一步提高对暴雨落区及雨强的预报水平.  相似文献   

GENERALIZED WAVE EQUATION FINITE ELEMENT METHOD FOR SOLVING TWO-DIMENSIONAL TIDAL WAVES     

吴伣康  赵汉中 《中国海洋湖沼学报》1992,10(4):301-312

The study of tidal circulation has a long history . The numerical simulation of tidal flow has been developed greatly with the development of computer techniques in the past two decades. The generalized wave equation finite-element method is a relatively new numerical model for studying shallow water flow . This method was used to simulate tidal waves of the Gulf of St. Lawrence in Canada . The very good agreement of the numerical results with the field data indicated that the model is an effective and promising numerical method for solving two-dimensional tidal wave problems .  相似文献   

MESOSCALE NONLINEAR GRAVITATIONAL WAVES AND THEIR INTERACTION IN VERTICAL SHEAR FLOW^*          下载免费PDF全文

ZHOU Weican  ZHU Lihu  ZOU Lanjun 《Acta Meteorologica Sinica》2004,18(4):491-504

Interaction equations of two nonlinear gravitational waves in baroclinic atmosphere are presented via multi-scale perturbation method,which can be classified into coupling nonlinear Schrodinger equations.In particular,the interaction course of two nonlinear gravitational waves of basic flow in vertical linear and quadratic shear is illustrated.Numerical calculation displays that wave amplitude enlarges and wave width narrows when two solitary gravitational waves meet and chase;that basic flow with single shear is more beneficial than that with quadratic shear to the interaction of two nonlinear wave packets;and that the interaction of two wave packets makes wave shape change more greatly and energy more dispersive,which contributes to the occurrence of changeable weather.Therefore,one of the probable mechanisms for the appearance of strong convection weather is the interaction between mesoscale nonlinear gravitational waves.  相似文献