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.  相似文献   

“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.  相似文献   

A Free Boundary Problem Related to the Location of Volcanic Gas Sources     

Jesús?Ildefonso DíazEmail author  Giorgio?Talenti 《Pure and Applied Geophysics》2004,161(7):1509-1517

A mathematical model of general gas emitting systems is derived, and a sample of relevant mathematical results is offered. The present paper indicates that shallow subsurface gas sources in typical volcanic areas can be located if appropriate physico-chemical measurements are made on the Earths surface and put to use.  相似文献   

On the use of STF-tensors in celestial mechanics     

Torsten Hartmann  Michael H. Soffel  Thomas Kioustelidis 《Celestial Mechanics and Dynamical Astronomy》1994,60(1):139-159

The purpose of this article is to emphasize the usefulness of STF-tensors in celestial mechanics. Using STF-mass multipole moments and Cartesian coordinates the derivations of equations of motion, the interaction- and tidal-potentials for an isolated system ofN arbitrarily shaped and composed, purely gravitationally interacting bodies are particularly simple. Using simple relations between STF-tensors and spherical harmonics it is shown how all Cartesian formulas can be converted easily into the usual spherical representations. Some computational aspects of STF-tensors and spherical harmonics are discussed. A list of useful formulas for STF-tensors is provided.  相似文献   

A unified formulation for the three-dimensional shallow water equations using orthogonal co-ordinates: theory and application     

Herman W. J. Kernkamp  Henri A. H. Petit  Herman Gerritsen  Erik D. de Goede 《Ocean Dynamics》2005,55(3-4):351-369

In this paper, the formulations of the primitive equations for shallow water flow in various horizontal co-ordinate systems and the associated finite difference grid options used in shallow water flow modelling are reviewed. It is observed that horizontal co-ordinate transformations do not affect the chosen co-ordinate system and representation in the vertical, and are the same for the three- and two-dimensional cases. A systematic derivation of the equations in tensor notation is presented, resulting in a unified formulation for the shallow water equations that covers all orthogonal horizontal grid types of practical interest. This includes spherical curvilinear orthogonal co-ordinate systems on the globe. Computational efficiency can be achieved in a single computer code. Furthermore, a single numerical algorithmic code implementation satisfies. All co-ordinate system specific metrics are determined as part of a computer-aided model grid design, which supports all four orthogonal grid types. Existing intuitive grid design and visual interpretation is conserved by appropriate conformal mappings, which conserve spherical orthogonality in planar representation. A spherical curvilinear co-ordinate solution of wind driven steady channel flow applying a strongly distorted grid is shown to give good agreement with a regular spherical co-ordinate model approach and the solution based on a β-plane approximation. Especially designed spherical curvilinear boundary fitted model grids are shown for typhoon surge propagation in the South China Sea and for ocean-driven flows through Malacca Straits. By using spherical curvilinear grids the number of grid points in these single model grid applications is reduced by a factor of 50–100 in comparison with regular spherical grids that have the same horizontal resolution in the area of interest. The spherical curvilinear approach combines the advantages of the various grid approaches, while the overall computational effort remains acceptable for very large model domains.  相似文献   

流体-结构相互作用下用微分求积法求弹性板上的动水压力场          下载免费PDF全文

张立翔  郭亚昆 《水科学进展》2006,17(6):767-773

采用微分求积法(DQ)在流体-结构相互作用(FSI)框架下求解原变量形式的Navier-Stokes(N-S)方程和板的弯曲振动方程,获取作用在弹性板上的动水压力场.分析对象是4块矩形板组成的箱型流道,其中两块为弹性薄板,另外两块为刚性板.为比较FSI对动水压力场的影响,数值结果给出了1.5m×0.5m×0.2m流道、来流雷诺数分别为000和10000、相应DQ网格为25×25×25和29×29×29考虑FSI与否的压力分布.数值计算结果表明,所建立的考虑FSI的DQ法数值求解原变量形式N-S方程和弹性板振动方程,获取弹性板上的动水压力场是可行的,FSI效应对动水压力场的影响是明显的.  相似文献   

The contribution of Very Long Baseline Interferometry to ITRF2005     

Markus Vennebusch  Sarah Böckmann  Axel Nothnagel 《Journal of Geodesy》2007,81(6-8):553-564

The contribution of the International VLBI Service for Geodesy and Astrometry (IVS) to the ITRF2005 (International Terrestrial Reference Frame 2005) has been computed by the IVS Analysis Coordinator’s office at the Geodetic Institute of the University of Bonn, Germany. For this purpose the IVS Analysis Centres (ACs) provided datum-free normal equation matrices in Solution INdependent EXchange (SINEX) format for each 24 h observing session to be combined on a session-by-session basis by a stacking procedure. In this process, common sets of parameters, transformed to identical reference epochs and a prioris, and especially representative relative weights have been taken into account for each session. In order to assess the quality of the combined IVS files, Earth orientation parameters (EOPs) and scaling factors have been derived from the combined normal equation matrices. The agreement of the EOPs of the combined normal equation matrices with those of the individual ACs in terms of weighted root mean square (WRMS) is in the range of 50–60 μas for the two polar motion components and about 3 μs for UT1−UTC. External comparisons with International GNSS Serive (IGS) polar motion components is at the level of 130–170 μas and 21 μs/day for length of day (LOD). The scale of the terrestrial reference frame realized through the IVS SINEX files agrees with ITRF2000 at the level of 0.2 ppb.  相似文献   

直线方程在工程测量中的运用     

翟警卫 《测绘与空间地理信息》2008,31(6)

通过实例,介绍一种在工程测量中从某一已知点上放样直线或直线上任意点确定其平面位置和高程位置的方法.此方法既可提高工作效率和放样精度,又能满足工程测量的需要.  相似文献