An arbitrary high-order discontinuous Galerkin method for elastic waves on unstructured meshes – I. The two-dimensional isotropic case with external source terms     

Martin Käser  Michael Dumbser 《Geophysical Journal International》2006,166(2):855-877

  相似文献   

电性源短偏移距瞬变电磁地形影响特征分析     

齐彦福  李貅  孙乃泉  戚志鹏  周建美  张明晶  邢涛 《吉林大学学报(地球科学版)》2022,52(1):247-260

电性源短偏移距瞬变电磁法是进行山区矿产资源勘探的一种有效手段,然而地形起伏对观测数据影响严重,给瞬变电磁数据精细化解释带来很大困难。为此,本文基于三维非结构时间域有限元正演模拟方法开展电性源短偏移距瞬变电磁地形影响特征分析。首先利用非结构网格对起伏地形进行精细剖分,结合矢量有限元和后退欧拉离散技术实现复杂地质模型快速正演。然后通过模拟大量典型起伏地表模型,系统分析了三维地形响应的空间分布特征,结果表明地形对电性源短偏移距瞬变电磁响应影响严重,而且影响特征十分复杂。在此基础上,我们进一步分析了复杂地形对地下目标体探测的影响,结果发现复杂地形与地下目标体相互耦合导致无法对地下目标体进行有效分辨,这给地下目标体的准确探测带来了巨大困难。  相似文献   

适用于移动终端的全球离散格网快速生成方法     

陈浩  许长辉  宋现锋  徐江  潘苗苗 《测绘科学》2015,(12):63-67

针对全球离散格网依据剖分方法和投影方式的不同,在实现方式、剖分速度、格网均匀性等方面呈现出较大差异的问题,该文顾及移动终端用户对高精度位置应用日益增长的需求,以及顾及移动终端计算性能较低的特点,分析了球面格网生成方法、格网向椭球面投影方式和几何变形特点。通过分析椭球面格网与球面格网间变形差异,该文认为球面格网代替椭球面格网在一定程度上能够满足不同精度的格网剖分需求。通过建立用户精度指标与剖分层次对应关系,提高了格网剖分的计算效率,可应用于移动终端等计算性能较低的设备。  相似文献   

An Evaluation of Publicly Available Global Bathymetry Grids     

K. M. Marks  W. H. F. Smith 《Marine Geophysical Researches》2006,27(1):19-34

We evaluate the strengths and weaknesses of six publicly available global bathymetry grids: DBDB2 (Digital Bathymetric Data Base; an ongoing project of the Naval Research Laboratory), ETOPO2 (Earth Topography; National Geophysical Data Center, 2001, ETOPO2 Global 2’ Elevations [CD-ROM]. Boulder, Colorado, USA: U.S. Department of Commerce, National Oceanic and Atmospheric Administration), GEBCO (General Bathymetric Charts of the Oceans; British Oceanographic Data Centre, 2003, Centenary Edition of the GEBCO Digital Atlas [CD-ROM] Published on behalf of the Intergovernmental Oceanographic Commission and the International Hydrographic Organization Liverpool, UK), GINA (Geographic Information Network of Alaska; Lindquist et al., 2004), Smith and Sandwell (1997), and S2004 (Smith, unpublished). The Smith and Sandwell grid, derived from satellite altimetry and ship data combined, provides high resolution mapping of the seafloor, even in remote regions. DBDB2, ETOPO2, GINA, and S2004 merge additional datasets with the Smith and Sandwell grid; but moving from a pixel to grid registration attenuates short wavelengths (<20 km) in the ETOPO2 and DBDB2 solutions. Short wavelengths in the GINA grid are also attenuated, but the cause is not known. ETOPO2 anomalies are offset to the northeast, due to a misregistration in both latitude and longitude. The GEBCO grid is interpolated from 500 m contours that were digitized from paper charts at 1:10 million scale, so it is artificially smooth; yet new efforts have captured additional information from shallow water contours on navigational charts. The S2004 grid merges the Smith and Sandwell grid with GEBCO over shallow depths and polar regions, and so is intended to capture the best of both products. Our evaluation makes the choice of which bathymetry grid to use a more informed one. The U.S. government right to retain a non-exclusive royalty-free license in and to any copyright is acknowledged.  相似文献   

Gravitational mode calculation of basins discretized by orthogonal curvilinear grids     

Gian Mario BeltramiRiccardo Briganti 《Ocean Engineering》2003,30(7):833-853

This paper presents a simple and straightforward method for carrying out the direct numerical solution of the eigenvalue problem associated to the homogeneous linear shallow-water equations expressed using orthogonal curvilinear coordinates, when ‘adiabatic’ boundary conditions apply. These equations, together with the boundary conditions, define a self-adjoint problem in the continuum. The method presented here, which is thought for calculating the 2-D theoretical gravity modes of both natural and artificial basins, relies on a change of basis of the dependent variable vector. This preliminary transformation makes it, in fact, possible to formulate two different numerical approaches which guarantee the self-adjoint property of the discrete form of the system consisting of the governing equations and the boundary conditions. The method is tested using a square and a fully circular domain, both of which allow comparisons with well-known analytical and numerical solutions. Discretizing the physical domain of a fully circular basin by a cylindrical coordinate grid makes it possible to show the actual efficiency of the method in calculating the theoretical gravity modes of basins discretized by a boundary-following coordinate grid which allows laterally variable resolution.  相似文献   

Diagnoses of vertical transport in a three-dimensional finite element model of the tidal circulation around an island     

Laurent White  Eric Deleersnijder   《Estuarine, Coastal and Shelf Science》2007,74(4):655

A three-dimensional finite element model is used to investigate the formation of shallow-water eddies in the wake of Rattray Island (Great Barrier Reef, Australia). Field measurements and visual observations show that stable eddies develop in the lee of the island at rising and falling tides. The water turbidity downstream of the island suggests the existence of strong upwelling that would be responsible for carrying bed sediments up to the sea surface. We first propose to look at the upwelling velocity and then use the theory of the age to diagnose vertical transport. The water age is defined as the time elapsed since particles of water left the sea bottom, where the age is prescribed to be zero. Two versions of this diagnosis are considered. Although the model predicts upwelling within the eddies, it is not sufficiently intense to account for vertical transport throughout the water column during the life span of the eddies. As mesh resolution increases, this upwelling does not intensify. However, strong upwelling is then resolved off the island's tips, which is confirmed by the results obtained with the age. This study also shows that the finite element method, together with unstructured meshes, performs well for representing three-dimensional flow past an island.  相似文献   

Numerical modeling of wind waves generated by tropical cyclones using moving grids   总被引：1，自引：0，他引：1  

Hendrik L. Tolman  Jose-Henrique G.M. Alves 《Ocean Modelling》2005,9(4):305-323

A version of the WAVEWATCH III wave model featuring a continuously moving spatial grid is presented. The new model option/version is intended for research into wind waves generated by tropical cyclones in deep water away from the coast. The main advantage of such an approach is that the cyclones can be modeled with spatial grids that cover much smaller areas than conventional fixed grids, making model runs with high spatial resolution more economically feasible. The model modifications necessary are fairly trivial. Most complications occur due to the Garden Sprinkler effect (GSE) and methods used to mitigate it. The basic testing of the model is performed using idealized wind fields consisting of a Rankine vortex. The model is also applied to hurricane Lili in the Gulf of Mexico in October 2002. The latter application shows that the moving grid approach provides a natural way to deal with hurricane wind fields that have a high-resolution in space, but a low resolution in time. Although the new model version is originally intended for tropical cyclones, it is suitable for high-resolution modeling of waves due to any moving weather pattern.  相似文献   

A Finite Volume Method with Unstructured Triangular Grids for Numerical Modeling of Tidal Current   总被引：1，自引：1，他引：1  

史宏达  刘臻 《中国海洋工程》2005,19(4):693-700

The finite volume method （FVM） has many advantages in 2-D shallow water numerical simulation. In this study, the finite volume method is used with unstructured triangular grids to simulate the tidal currents. The Roe scheme is applied in the calculation of the intercell numerical flux, and the MUSCL method is introduced to improve its accuracy. The time integral is a two-step scheme of forecast and revision. For the verification of the present method, the Stoker＇s problem is calculated and the result is compared with the mathematically analytic solutions. The comparison indicates that the method is feasible. A sea area of a port is used as an example to test the method established here. The result shows that the present computational method is satisfactory, and it could be applied to the engineering fields.  相似文献   

Horizontal Pressure Gradient Errors of the Monterey Bay Sigma Coordinate Ocean Model with Various Grids     

Le Ngoc Ly  Lin Jiang 《Journal of Oceanography》1999,55(1):87-97

A coastal ocean -coordinate model of Monterey Bay (MOB) with realistic bottom topography and coastlines is developed using the Princeton Ocean Model (POM) and grid generation technique (GGT) to study the horizontal pressure gradient errors associated with the MOB steep topography. The submarine canyon in MOB features some of the steepest topography encountered anywhere in the world oceans. The MOB grids are designed using the EAGEAL View and GENIE⁺⁺ grid generation systems. A grid package developed by Ly and Luong (1993) is used in this study to couple grids to the model. The MOB model is tested with both orthogonal and curvilinear nearly-orthogonal (CNO) grids. The CNO grid has horizontal resolution which varies from 300 m to 2 km, while the resolution of the orthogonal grid is uniform with x = 1.25 km and y = 1.38 km. These grids cover a domain of 180 × 160 km with the same number of grid points of 131 × 131. Vertical resolutions of 25, 35 and 45 vertical sigma levels are tested. The error in the MOB are evaluated in terms of mean kinetic energy and velocity against various grids, vertical, horizontal resolution and distributions, and bottom topography smoothing. Simulations with various grids show that GGT can be used as another tool in reducing -coordinate errors in coastal ocean modeling besides increasing resolution and smoothing bottom topography. Topographical smoothing not only reduces topographic slope, but changes realistic topography. A CNO grid with a high grid density packed along steep slopes and Monterey Submarine Canyon reduces the errors by 40% compared to a rectangular grid with the same number of grid points. The CNO grid is more efficient than the rectangular grid, since it has most of its grids over water. The simulations show that the presented MOB -coordinate model can be used with a confidence regarding horizontal pressure gradient error.  相似文献   

一种面向虚拟地球的海面动态可视化优化方法     

黄吴蒙  陈静 《测绘学报》2016,45(Z1):135-143

现有的虚拟地球中的海面可视化方法主要采用投影网格的方式组织海面格网,这种特殊的网格组织方式,导致其在反映不同海域差异性特征时存在缺陷。基于全球离散格网的海面可视化方法由于能够与虚拟地球离散的空间剖分相匹配,有效弥补了投影网格法的缺陷,因此更符合虚拟地球海面仿真应用的需要。然而现有的离散格网法存在绘制效率差、场景加载慢、需要修补格网缝隙等问题,限制了其应用发展。对此,本文在现有离散格网法基础上进行优化:首先,在数据结构上对传统等经纬度离散格网进行扩展,设计了一种面向GPU绘制的多尺度海面网格模型来组织管理海面格网;其次,为了实现风场驱动下的海浪动态绘制,在多尺度海面网格模型的基础上提出了一种支持实时风场更新的海浪动态绘制方法;同时,考虑到格网缝隙修补对系统效率的影响,本文针对海面格网特点并结合GPU技术提出了一种高效的海面格网缝隙修补方法。最后通过对比试验验证了本文方法的可行性和有效性。试验结果表明,本文方法绘制效率稳定,加载速度快,且能弥补现有方法在功能上的缺失,因此应用范围更加广泛。  相似文献