Thi-Thao-Phuong Hoang

Thi-Thao-Phuong Hoang, Caroline Japhet, Michel Kern et al.

In this paper we are interested in the "fast path" fracture and we aim to use global-in-time, nonoverlapping domain decomposition methods to model flow and transport problems in a porous medium containing such a fracture. We consider a reduced model in which the fracture is treated as an interface between the two subdomains. Two domain decomposition methods are considered: one uses the time-dependent SteklovPoincar{é} operator and the other uses optimized Schwarz waveform relaxation (OSWR) based on Ventcell transmission conditions. For each method, a mixed formulation of an interface problem on the space-time interface is derived, and different time grids are employed to adapt to different time scales in the subdomains and in the fracture. Demonstrations of the well-posedness of the Ventcell subdomain problems is given for the mixed formulation. An analysis for the convergence factor of the OSWR algorithm is given in the case with fractures to compute the optimized parameters. Numerical results for two-dimensional problems with strong heterogeneities are presented to illustrate the performance of the two methods.

Thi-Thao-Phuong Hoang, Lili Ju, Zhu Wang

The paper is concerned with overlapping domain decomposition and exponential time differencing for the diffusion equation discretized in space by cell-centered finite differences. Two localized exponential time differencing methods are proposed to solve the fully discrete problem: the first method is based on Schwarz iteration applied at each time step and involves solving stationary problems in the subdomains at each iteration, while the second method is based on the Schwarz waveform relaxation algorithm in which time-dependent subdomain problems are solved at each iteration. The convergence of the associated iterative solutions to the corresponding fully discrete multidomain solution and to the exact semi-discrete solution is rigorously proved. Numerical experiments are carried out to confirm theoretical results and to compare the performance of the two methods.

Thi-Thao-Phuong Hoang, Caroline Japhet, Michel Kern et al.

This paper is concerned with the numerical solution of porous-media flow and transport problems , i. e. heterogeneous, advection-diffusion problems. Its aim is to investigate numerical schemes for these problems in which different time steps can be used in different parts of the domain. Global-in-time, non-overlapping domain-decomposition methods are coupled with operator splitting making possible the different treatment of the advection and diffusion terms. Two domain-decomposition methods are considered: one uses the time-dependent Steklov--Poincar{é} operator and the other uses optimized Schwarz waveform relaxation (OSWR) based on Robin transmission conditions. For each method, a mixed formulation of an interface problem on the space-time interface is derived, and different time grids are employed to adapt to different time scales in the subdomains. A generalized Neumann-Neumann preconditioner is proposed for the first method. To illustrate the two methods numerical results for two-dimensional problems with strong heterogeneities are presented. These include both academic problems and more realistic prototypes for simulations for the underground storage of nuclear waste.

Thi-Thao-Phuong Hoang, Lili Ju, Wei Leng et al.

In this paper we present and analyze a general framework for constructing high order explicit local time stepping (LTS) methods for hyperbolic conservation laws. In particular, we consider the model problem discretized by Runge-Kutta discontinuous Galerkin (RKDG) methods and design LTS algorithms based on strong stability preserving Runge-Kutta (SSP-RK) schemes, that allow spatially variable time step sizes to be used for time integrations in different regions. The proposed algorithms are of predictor-corrector type, in which the interface information along the time direction is first predicted based on the SSP-RK approximations and Taylor expansions, and then the fluxes over the region of interface are corrected to conserve mass exactly at each time step. Following the proposed framework, we detail the corresponding LTS schemes with accuracy up to the fourth order, and prove their conservation property and nonlinear stability for the scalar conservation laws. Numerical experiments are also presented to demonstrate excellent performance of the proposed LTS algorithms.

Thanh Hai Ong, Duc Cam Hai Vo, Thi-Thao-Phuong Hoang

The paper is concerned with the derivation and analysis of nonoverlapping domain decomposition for heterogeneous, anisotropic diffusion problems discretized by the finite element cell-centered (FECC) scheme. Differently from the standard finite element method (FEM), the FECC method involves only cell unknowns and satisfies local conservation of fluxes by using a technique of dual mesh and multipoint flux approximations to construct the discrete gradient operator. Consequently, if the domain is decomposed into nonoverlapping subdomains, the transmission conditions (on the interfaces between subdomains) associated with the FECC scheme are different from those of the standard FEM. However, the substructuring procedure as well as the Neumann-Neumann type preconditioner can be adapted to the domain decomposition-based FECC method naturally. Convergence analysis of a preconditioned iterative algorithm, namely the Dirichlet-Neumann to Neumann-Neumann algorithm, associated with the discrete FECC interface problem is the main focus of this work. Two dimensional numerical results for two subdomains with conforming meshes demonstrate that the preconditioned iterative algorithm converges independently of the mesh size and the coefficient jump.