基于非结构化网格的2.5-D直流电阻率自适应有限元数值模拟   总被引：3，自引：2，他引：1       下载免费PDF全文

汤井田  王飞燕  任政勇 《地球物理学报》2010,53(3):708-716

2.5-D直流电阻率有限元数值模拟中,模型的剖分及加密主要通过手动实现.另外,采用的单元类型比较规则如矩形单元等,不易实现复杂模型的模拟.为解决上述问题,文中提出了一种自适应有限元算法.算法中采用稳健的后验误差估计来自动预测下一次网格的单元尺寸,直到设定的迭代条件满足为止.另外,采用非结构化三角形单元实现了任意复杂模型的灵活剖分.基于此,利用垂直接触面模型分析和对比了不同自适应策略的效率.通过对比发现,点源附近的单元得到了加密以消除源的奇异性.另外,对于任意一种策略,有限元结果均能最终收敛到精确解.最后,模拟了两个模型：2-D单个异常体模型和2-D地形模型.  相似文献   

TWO DIMENSIONAL NUMERICAL SIMULATION OF FLOW AND GEO-MORPHOLOGICAL PROCESSES NEAR HEADLANDS BY USING UNSTRUCTURED GRID     

A.MOHAMMADIAN  M.TAJRISHI  F.LOTFIAZAD 《国际泥沙研究》2004,19(4)

1 INTRODUCTION In recent years, due to the increase in population and industrial developments, mankind has faced manyproblems associated with rivers, coastal waters and reservoirs. Some of these problems are flood control,water supply, power generation, and irrigation. In addition, making new hydraulic structures changesnatural conditions. Prediction of these changes is necessary for designing such constructions. For solutionof these problems usually an assessment of flow pattern, sedim…  相似文献   

基于非结构化网格的边界积分方程方法的断层自发破裂模拟          下载免费PDF全文

钱峰  吴葆宁  冯禧  张海明 《地球物理学报》2019,62(9):3421-3431

地震自发破裂模拟是震源动力学研究的重要内容,了解复杂的断层动力学破裂过程对深入认识震源特征和解释运动学反演结果具有重要意义.基于边界积分方程方法的破裂模拟已经被广泛使用,大多采用的是平面断层模型的结构化网格划分.由于实际的断层往往具有较为复杂的几何特征,为了更为灵活地刻画断层几何复杂性,我们建立断层模型的三角形网格离散方案,通过精确的解析解形式来计算断层各个单元之间的应力格林函数,联立滑动弱化摩擦准则和非奇异边界积分方程,对断层的自发破裂过程进行了模拟.在简单的平面断层模型下,将计算结果与前人的结果进行了对比,验证了方法的正确性与有效性.对于几种常见的复杂断层模型,例如弯折、阶跃、含障碍体断层等,我们模拟了其破裂过程并对计算结果进行了比较与分析.模拟结果表明,非结构化网格划分的边界积分方程方法能够很好地模拟平面矩形断层或由其组成的规则断层,同时也能成功地模拟具有复杂几何形状的不规则断层上的动力学破裂过程.本研究的结果显示了边界积分方程方法在模拟复杂断层系统的动力学破裂问题上具有较广阔的应用前景.  相似文献   

二维TI介质非结构网格弹性波矢量逆时偏移          下载免费PDF全文

芦永明  张剑锋  杨凯  胡锦银 《地球物理学报》2017,60(12):4776-4789

波场延拓得到的多分量波场中既包含纵波信息也包含横波信息,能否在全波场中实现纵横波的分离对各向同性和各向异性逆时偏移都有非常重要的意义.传统的散度旋度分离只适应于各向同性介质而对各向异性介质却无效.在非规则、非结构网格的弹性波数值模拟方法的基础上,发展了一种适应于各向异性介质的波场分离方法.该方法通过求解Christoffel方程,得到相角和极化角的关系,再利用群角和相角的关系,直接得到群角和极化角的关系.该方法与现存的各向异性波场分离相比,获得的计算效率改进更显著,而且存储量小.用简单各向异性模型和SEG各向异性Hess模型进行测试,都得到了较好的效果,证明了本文方法的有效性.  相似文献   

基于微增模型的海洋可控源电磁法三维非结构化矢量有限元数值模拟          下载免费PDF全文

陈汉波  李桐林  熊彬  王恒  张镕哲  李少朋 《地球物理学报》2018,61(6):2560-2577

为了分析海洋中含金属矿的储油气砂层的成分、结构对于海洋可控源电磁场的影响特征,引入多相微增等效介质模型,结合非结构化网格实现了海洋可控源电磁三维矢量有限元正演.首先对多相微增模型进行了介绍,并分析不同参数影响下的等效电导率变化特征.而后,详细推导了频率域可控源电磁法非结构化矢量有限元方程.接着,利用预处理的IDR （s）迭代算法求解线性方程组.最后设计典型的海洋地电模型,通过正演计算,验证了本文算法的正确性和有效性,同时研究了海洋储油砂层中的金属矿的含量、孔隙度、含水饱和度等参数对于可控源电磁场响应特征的影响.  相似文献   

鄱阳湖二维水动力和水质耦合数值模拟   总被引：5，自引：3，他引：2  

赖锡军  姜加虎  黄群  徐力刚 《湖泊科学》2011,23(6):893-902

针对大型通江湖泊水位变化剧烈,地形起伏多变,岸线复杂,湖泊内部窄小洪道与大面洲滩和洼地连接复杂,湖泊洲滩出露、淹没频繁交替等湖泊水情和地貌特征,基于二维浅水方程和对流扩散方程组构建了鄱阳湖二维水动力和水质耦合模拟模型.模型采用非结构网格有限体积法进行离散,以HLLC算法计算单元界面的水量、动量和物质输运通量.水陆边界通...  相似文献   

Two-dimensional scour simulations based on coupled model of shallow water equations and sediment transport on unstructured meshes     

X. Liu  B.J. Landry  M.H. García 《Coastal Engineering》2008

A two-dimensional scour model based on coupled system of shallow water equations (SWEs) and sediment transport on unstructured mesh is developed. The coupled system of hydrodynamic and morphodynamic equations is solved by finite volume method using Godunov scheme. Roe's approximate Riemann solver is used to calculate the inviscid fluxes. The use of unstructured mesh makes the model applicable to complex domains. However, it is difficult to evaluate the eigenvalues and eigenvectors of the Jacobian matrix in the global coordinate. The method proposed herein to deal with this difficulty is to transform the system into the local coordinate with one of the axes in the same direction as the interface outward normal vector. In the local coordinate system, the Jacobian matrix is simplified and the eigenvalues are analyzed using asymptotic method. Regular expansion breaks down when the flow is near critical. Uniformity of the expansion is achieved by changing the scales. Rotational invariance theorem is used to relate the interfacial fluxes in the global and local coordinate systems. Special treatment of the source term on unstructured grid makes the scheme stable and physically balanced (both mass and momentum). The method proposed in this paper for the eigen-system is very efficient comparing to iterative numerical methods. Results from the test cases show good agreement with the experiments.  相似文献   

DG-FEM solution for nonlinear wave-structure interaction using Boussinesq-type equations     

Allan P. Engsig-Karup  Jan S. Hesthaven  Harry B. Bingham  T. Warburton 《Coastal Engineering》2008

  相似文献   

A three-dimensional coastal barotropic model in generalized curvilinear grids     

鲍献文  孙文心  史峰岩 《中国海洋湖沼学报》1999,17(4):289-299

INTRODUCTIONThemethodofautomaticallygeneratinggeneralizedcurvilinearmeshesisaneffectivetoolforobtainingnumericalsolutionsofpartialdifferentialequationsinaregionwitharbitraryboundaries.Thegrids,calledadaptivemesh ,canbeadaptedtotheshapeoftheboundaryortothespatialdis tributionofthesimulatedfield .Insteadofthetraditionalequidistantgrids,automaticallyproducedcurvilineargrids,whichcanbeadaptedtotheequationsolutions’featuresandtotheirregularshapeofthelateralboundary,areused .Thenumericalcalcul…  相似文献   

The effects of resolving the Canadian Arctic Archipelago in a finite element sea ice model   总被引：1，自引：0，他引：1  

Olivier Lietaer  Thierry Fichefet  Vincent Legat 《Ocean Modelling》2008,24(3-4):140-152

Though narrow straits may have a strong influence on the large-scale sea ice mass balance, they are often crudely represented in coarse resolution sea ice models. Unstructured meshes, with their natural ability to fit boundaries and locally increase the mesh resolution, propose an alternative framework to capture the complex oceanic areas formed by coasts and islands. In this paper, we develop a finite element sea ice model to investigate the sensitivity of the Arctic sea ice cover features to the resolution of the narrow straits constituting the Canadian Arctic Archipelago. The model is a two-level dynamic-thermodynamic sea ice model, including a viscous-plastic rheology. It is run over 1979–2005, forced by daily NCEP/NCAR reanalysis data. Confronting qualitatively numerical experiments with observations shows a good agreement with satellite and buoys measurements. Due to its simple representation of the oceanic interactions, the model overestimates the sea ice extent during winter in the southernmost parts of the Arctic, while the Baffin Bay and Kara Sea remain ice-covered during summer. In order to isolate the benefits from resolving the Canadian Arctic Archipelago, a numerical experiment is performed where we artificially close the archipelago. Focusing on the large-scale sea ice thickness pattern, no significant change is found in our model, except in the close surroundings of the archipelago. However, the local and short-term influences of the ice exchanges are nonnegligible. In particular, we show that the ice volume associated to the Canadian Arctic Archipelago represents 10% of the Northern Hemisphere sea ice volume and that the annual mean ice export towards Baffin Bay amounts to 125 km³ yr⁻¹, which may play an important role on the convective overturning in the Labrador Sea.  相似文献