No abstract uploaded!

We review some important ideas in the design and analysis of robust overlapping domain decomposition algorithms for high-contrast multiscale problems.

We use a suite of cosmological hydrodynamic simulations, run by two fixed grid codes, in the context of LCDM model to investigate the properties of solenoidal and dilatational, i.e, curl and compressive, motions of the intergalactic medium (IGM), and the impact of numerical viscosity on turbulence. The codes differ only in the spatial difference discretization, one is a second-order TVD scheme and the other is a fifth-order positivity-preserving WENO scheme. We find that (1) The vortical motion grows rapidly since $z=2$, and reaches $\sim 10 km/s -75 km/s$, depending on the gas density, at low redshifts. The small-scale compressive ratio $r_{CS}$ drops from $0.8$ to $0.5$ during the same period. (2) Power spectra of the solenoidal velocity possess two regimes, $\propto k^{-0.89}$ and $\propto k^{-2.02}$, while the total velocity and dilatational velocity follow the scaling $k^{-1.88}$ and $k^{-2.20}$ respectively all through the turbulent scale $\sim 2-0.2 h^{-1}Mpc$. The IGM turbulence may contain two different phases, the supersonic phase and the post-supersonic phase, simultaneously along the hierarchical structure formation history. (3) The non-thermal pressure support, measured by the vortical kinetic energy, is comparable with the thermal pressure for $\rho_b \sim 10$, or $T \sim 10^5 K$ at $z=0.25$. The deviation of the baryon fraction from the cosmic mean shows a positive correlation with the turbulence pressure support. (4) A relatively higher numerical viscosity would artificially dissipate both the compressive and vortical motions in the IGM more seriously, resulting in less developed density fluctuation and vorticity, and leading to remarkably shortened turbulence scale. Shocks in regions out of the clusters are outstandingly suppressed by the numerical viscosity since $z=2$, which could directly cause the different levels of turbulence in the IGM between the codes. We also discuss the impact of ingredients not included in our simulation.

In [H. Guo, Q. Zhang, J. Wang, Applied Mathematics and Computation, 259 (2015), 88-105], a nonconservative local discontinuous Galerkin (LDG) method for both flow and transport equations was introduced for the one-dimensional coupled system of compressible miscible displacement problem. In this paper, we will continue our effort and develop a conservative LDG method for the problem in two space dimensions. Optimal error estimates in $L^{\infty}(0, T; L^{2})$ norm for not only the solution itself but also the auxiliary variables will be derived. The main difficulty is how to treat the inter-element discontinuities of two independent solution variables (one from the flow equation and the other from the transport equation) at cell interfaces. Numerical experiments will be given to confirm the accuracy and efficiency of the scheme.

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