Operator Splitting and Local Time-Stepping Methods for Transport Problems in Fractured Porous Media  

作　　者：Phuoc-Toan Huynh Yanzhao Cao Thi-Thao-Phuong Hoang 

机构地区：[1]Department of Mathematics and Statistics,Auburn University,Auburn 36849,Alabama,USA

出　　处：《Communications in Computational Physics》2023年第10期1215-1246,共32页计算物理通讯（英文）

基　　金：partially supported by the US National Science Foundation under grant number DMS-1912626.

摘　　要：This paper is concerned with efficient numerical methods for the advectiondiffusion equation in a heterogeneous porous medium containing fractures.A dimensionally reduced fracture model is considered,in which the fracture is represented as an interface between subdomains and is assumed to have larger permeability than the surrounding area.We develop three global-in-time domain decomposition methods coupled with operator splitting for the reduced fracture model,where the advection and the diffusion are treated separately by different numerical schemes and with different time steps.Importantly,smaller time steps can be used in the fracture-interface than in the subdomains.The first two methods are based on the physical transmission conditions,while the third one is based on the optimized Schwarz waveform relaxation approach with Ventcel-Robin transmission conditions.A discrete space-time interface system is formulated for each method and is solved iteratively and globally in time.Numerical results for two-dimensional problems with various P′eclet numbers and different types of fracture are presented to illustrate and compare the convergence and accuracy in time of the proposed methods with local time stepping.

关 键 词：Domain decomposition reduced fracture model operator splitting nonconforming time grids time-dependent Steklov-Poincare operator optimized Schwarz waveform relaxation 

分 类 号：O17[理学—数学]