MAST-2D diffusive model for flood prediction on domains with triangular Delaunay unstructured meshes     

C. Aricò  M. Sinagra  L. Begnudelli  T. Tucciarelli 《Advances in water resources》2011,34(11):1427-1449

A new methodology for the solution of the 2D diffusive shallow water equations over Delaunay unstructured triangular meshes is presented. Before developing the new algorithm, the following question is addressed: it is worth developing and using a simplified shallow water model, when well established algorithms for the solution of the complete one do exist?The governing Partial Differential Equations are discretized using a procedure similar to the linear conforming Finite Element Galerkin scheme, with a different flux formulation and a special flux treatment that requires Delaunay triangulation but entire solution monotonicity. A simple mesh adjustment is suggested, that attains the Delaunay condition for all the triangle sides without changing the original nodes location and also maintains the internal boundaries. The original governing system is solved applying a fractional time step procedure, that solves consecutively a convective prediction system and a diffusive correction system. The non linear components of the problem are concentrated in the prediction step, while the correction step leads to the solution of a linear system of the order of the number of computational cells. A semi-analytical procedure is applied for the solution of the prediction step. The discretized formulation of the governing equations allows to handle also wetting and drying processes without any additional specific treatment. Local energy dissipations, mainly the effect of vertical walls and hydraulic jumps, can be easily included in the model.Several numerical experiments have been carried out in order to test (1) the stability of the proposed model with regard to the size of the Courant number and to the mesh irregularity, (2) its computational performance, (3) the convergence order by means of mesh refinement. The model results are also compared with the results obtained by a fully dynamic model. Finally, the application to a real field case with a Venturi channel is presented.  相似文献   

基于Midas/GTS的FLAC3D边坡建模技术及工程应用     

钟志辉  刘祚秋  杨光华  张玉成 《地震工程学报》2011,33(Z1)

探讨了如何基于Midas/GTS来实现用FLAC3D建立复杂边坡模型。首先在Midas/GTS中建立几何模型和划分网格,然后根据两种软件的网格数据形式进行转换,把Midas/GTS的网格模型导入FLAC3D中。对于二维边坡问题,认为采用四边形单元比三角形单元的计算精度更高,在模型较为复杂的情况下建议在Midas/GTS中采用四边形+三角形的方式生成混合网格,然后导入FLAC3D中形成brick+wedge形式的混合网格进行分析,以保证计算精度。最后将上述建模方法应用于某边坡工程的稳定性分析,验证了该建模方法的可行性和实用性。  相似文献   

关于机载LiDAR点云数据形态学滤波的几点思考     

董保根  秦志远  朱传新  徐验兵 《测绘科学》2013,38(4)

形态学滤波是一种典型的机载LiDAR点云数据滤波方法,然而滤波过程中还存在某些容易被忽视的细节.本文在回顾现有各种形态学滤波算法的基础上,分析了特殊地形下点云数据内插存在的问题,提出了改进的内插算法;深入探讨了两种不同数据组织形式的滤波算法,指出基于不规则离散点和基于规则格网两种滤波方法各自的优势及不足;最后通过实验数据对两种方法在滤波性能和效率上给出了定性分析和定量评价.  相似文献   

Matlab环境下瑞利波有限差分正演与曲线绘制   总被引：1，自引：0，他引：1  

杨天春  朱自强  周勇 《煤田地质与勘探》2010,38(1):62

以瑞利波频散方程为出发点的Abo-Zena传递矩阵等方法,只能研究层状介质中瑞利波的传播特性,对于非层状介质,传递矩阵无能为力。因此,不得不考虑有限元、有限差分等方法。根据弹性动力学方程,采用交错网格有限差分方法对均匀弹性半空间介质进行全波场模拟,并在Matlab环境下实现编程计算,再现了瑞利波在近地表的传播状态,实现了地震剖面曲线的正振幅充填问题,从而获得了类似野外地震勘探的剖面记录。模拟结果表明,采用一阶差分格式所得到的地震记录存在较严重的数值频散,在采用有限差分法研究瑞利波"之"字形成因时,须尽量压制数值频散现象,才能获得更接近实际情况的频散曲线。   相似文献   

A two-dimensional morphological model based on a next generation circulation solver I: Formulation and validation     

Cui-Ping Kuang  Si-Yu Chen  Yu Zhang  Jie GuYi Pan  Jing Huang 《Coastal Engineering》2012,59(1):1-13

  相似文献   

电性参数分块连续变化二维MT有限元数值模拟   总被引：4，自引：2，他引：2       下载免费PDF全文

刘云  王绪本 《地球物理学报》2012,55(6):2079-2086

为了易于模拟野外复杂地形和地下任意形状地电体模型,将有限元单元网格设计为三角单元;并考虑到野外实际勘探中,地球介质的电性参数均是连续变化的情况,单元内的场值和电性参数被设计为双线性变化;推导出二维起伏地形条件下大地电磁法有限元数值模拟算法;根据单元节点主场值和线性插值形函数间的关系,计算出单元节点的辅助场值;在二维起伏地形情况下,定义TE、TM模式视电阻率和阻抗相位.4个模型的计算的结果与解析法的均方根误差小于1%,地形模拟与前人的计算结果相符,模拟倾斜界面异常体,能有效的反映出其异常形态.  相似文献   

Mesh type tradeoffs in 2D hydrodynamic modeling of flooding with a Godunov-based flow solver     

《Advances in water resources》2014

The effect of mesh type on the accuracy and computational demands of a two-dimensional Godunov-type flood inundation model is critically examined. Cartesian grids, constrained and unconstrained triangular grids, constrained quadrilateral grids, and mixed meshes are considered, with and without local time stepping (LTS), to determine the approach that maximizes computational efficiency defined as accuracy relative to computational effort. A mixed-mesh numerical scheme is introduced so all grids are processed by the same solver. Analysis focuses on a wide range of dam-break type test cases, where Godunov-type flood models have proven very successful. Results show that different mesh types excel under different circumstances. Cartesian grids are 2–3 times more efficient with relatively simple terrain features such as rectilinear channels that call for a uniform grid resolution, while unstructured grids are about twice as efficient in complex domains with irregular terrain features that call for localized refinements. The superior efficiency of locally refined, unstructured grids in complex terrain is attributable to LTS; the locally refined unstructured grid becomes less efficient using global time stepping. These results point to mesh-type tradeoffs that should be considered in flood modeling applications. A mixed mesh model formulation with LTS is recommended as a general purpose solver because the mesh type can be adapted to maximize computational efficiency.  相似文献   

Petrology of phlogopite-sapphirine-bearing Al-Mg granulites from Ihouhaouene, In Ouzzal, Hoggar, Algeria: an example of phlogopite stability at high temperature   总被引：1，自引：0，他引：1  

H. MOURI  M. GUIRAUD  B. J. HENSEN 《Journal of Metamorphic Geology》1996,14(6):725-738

Al-Mg granulites, with cordierite, garnet, sapphirine, orthopyroxene, sillimanite, spinel, phlogopite, K-feldspar, plagioclase and variable quartz from Ihouhaouene (In Ouzzal, Algeria), display a range of decompression textures involving the breakdown of orthopyroxene and sillimanite, and of garnet. The succession of parageneses suggests that the P–T–t evolution corresponds to decompression with cooling from peak conditions of about 950°C and 10 kbar. This decompression path is obtained from the paragenetic analysis in the FMAS system. However, according to current KFMASH grids, this P–T–t evolution should take place outside the stability field of phlogopite+quartz; yet this assemblage is probably stable during most of the P-T evolution, notably during peak metamorphism. This discrepancy is interpreted as the effect of the high content of F in phlogopite which should shift its stability limit towards higher temperature. The consequences of this shift on the phase relationships in the KFeMASH system are investigated and it is concluded that a topological inversion could exist in the F-bearing system.  相似文献   

Calculation of Well Index for Nonconventional Wells on Arbitrary Grids     

Christian Wolfsteiner  Louis J. Durlofsky  Khalid Aziz 《Computational Geosciences》2003,7(1):61-82

Wells are seldom modeled explicitly in large scale finite difference reservoir simulations. Instead, the well is coupled to the reservoir through the use of a well index, which relates wellbore flow rate and pressure to grid block quantities. The use of an accurate well index is essential for the detailed modeling of nonconventional wells; i.e., wells with an arbitrary trajectory or multiple branches. The determination of a well index for such problems is complicated, particularly when the simulation grid is irregular or unstructured. In this work, a general framework for the calculation of accurate well indices for general nonconventional wells on arbitrary grids is presented and applied. The method entails the use of an accurate semianalytical well model based on Green's functions as a reference single phase flow solution. This result is coupled with a finite difference calculation to provide an accurate well index for each grid block containing a well segment. The method is demonstrated on a number of homogeneous example cases involving deviated, horizontal and multilateral wells oriented skew to the grid. Both Cartesian and globally unstructured multiblock grids are considered. In all these cases, the method is shown to provide results that are considerably more accurate compared to results using standard procedures. The method is also applied to heterogeneous problems involving horizontal wells, where it is shown to be capable of approximating the effects of subgrid heterogeneity in coarse finite difference models.  相似文献   

Mass-conservative finite volume methods on 2-D unstructured grids for the Richards’ equation     

Gianmarco Manzini  Stefano Ferraris   《Advances in water resources》2004,27(12):215

The solution to the 2-D time-dependent unsaturated flow equation is numerically approximated by a second-order accurate cell-centered finite-volume discretization on unstructured grids. The approximation method is based on a vertex-centered Least Squares linear reconstruction of the solution gradients at mesh edges.A Taylor series development in time of the water content dependent variable in a finite-difference framework guarantees that the proposed finite volume method is mass conservative. A Picard iterative scheme solves at each time step the resulting non-linear algebraic problem. The performance of the method is assessed on five different test cases and implementing four distinct soil constitutive relationships. The first test case deals with a column infiltration problem. It shows the capability of providing a mass-conservative behavior. The second test case verifies the numerical approximation by comparison with an analytical mixed saturated–unsaturated solution. In this case, the water drains from a fully saturated portion of a 1-D column. The third and fourth test cases illustrate the performance of the approximation scheme on sharp soil heterogeneities on 1-D and 2-D multi-layered infiltration problems. The 2-D case shows the passage of an abrupt infiltration front across a curved interface between two layers. Finally, the fifth test case compares the numerical results with an analytical solution that is developed for a 2-D heterogeneous soil with a source term representing plant roots. This last test case illustrates the formal second-order accuracy of the method in the numerical approximation of the pressure head.  相似文献