In this paper, we apply local discontinuous Galerkin methods to the compressible wormhole propagation. Optimal error estimates for the pressure, velocity, porosity and concentration in different norms are established on non-uniform grids. Numerical experiments are presented to verify the theoretical analysis and show the good performance of the LDG scheme for compressible wormhole propagation.

Mortar methods are widely used techniques for discretizations of partial differential equations and preconditioners for the algebraic systems resulting from the discretizations. For problems with high contrast and multiple scales, the standard mortar spaces are not robust, and some enrichments are necessary in order to obtain an efficient and robust mortar space. In this paper, we consider a class of flow problems in high contrast heterogeneous media, and develop a systematic approach to obtain an enriched multiscale mortar space. Our approach is based on the constructions of local multiscale basis functions. The multiscale basis functions are constructed from local problems by following the framework of the Generalized Multiscale Finite Element Method (GMsFEM). In particular, we first create a local snapshot space. Then we select the dominated modes within the snapshot space using an appropriate Proper Orthogonal Decomposition (POD) technique. These multiscale basis functions show better accuracy than polynomial basis for multiscale problems. Using the proposed multiscale mortar space, we will construct a multiscale finite element method to solve the flow problem on a coarse grid and a preconditioning technique for the fine scale discretization of the flow problem. In particular, we develop a multiscale mortar mixed finite element method using the mortar space. In addition, we will design a two-level additive preconditioner and a two-level hybrid preconditioner based on the proposed mortar space for the iterative method applied to the fine scale discretization of the flow problem. We present several numerical examples to demonstrate

We consider an inverse problem of recovering a time-dependent factor of an unknown source on some subboundary for a diffusion equation with time fractional derivative by nonlocal measurement data. Such fractional-order equations describe anomalous diffusion of some contaminants in heterogeneous media such as soil and model the contamination process from an unknown source located on a part of the boundary of the concerned domain. For this inverse problem, we firstly establish the well-posedness in some Sobolev space. Then we propose two regularizing schemes in order to reconstruct an unknown boundary source stably in terms of the noisy measurement data. The first regularizing scheme is based on an integral equation of the second kind which an unknown boundary source solves, and we prove a convergence rate of regularized solutions with a suitable choice strategy of the regularizing parameter. The second regularizing scheme relies directly on discretization by the radial basis function for the initial-boundary value problem for fractional diffusion equation, and we carry out numerical tests, which show the validity of our proposed regularizing scheme.

No abstract uploaded!

We present a regularized finite difference method for the logarithmic Schr ̈odingerequation (LogSE) and establish its error bound. Due to the blowup of the logarithmic nonlinearity, i.e., ln ρ→ −\infty when ρ→0+v with ρ = |u|^2 being the density and u being the complex-valuedwave function or order parameter, there are significant difficulties in designing numerical methodsand establishing their error bounds for the LogSE. In order to suppress the roundoff error and toavoid blowup, a regularized LogSE (RLogSE) is proposed with a small regularization parameter 0 < ε << 1 and linear convergence is established between the solutions of RLogSE and LogSE interm of ε. Then a semi-implicit finite difference method is presented for discretizing the RLogSEand error estimates are established in terms of the mesh sizehand time stepτas well as the smallregularization parameterε. Finally numerical results are reported to illustrate our error bounds.