We present a new stabilization technique for multiscale convection diffusion problems. Stabilization for these problems has been a challenging task, especially for the case with high Peclet numbers. Our method is based on a constraint energy minimization idea and the discontinuous Petrov-Galerkin formulation. In particular, the test functions are constructed by minimizing an appropriate energy subject to certain orthogonality conditions, and are related to the trial space. The resulting test functions have a localization property, and can therefore be computed locally. We will prove the stability, and present several numerical results. Our numerical results confirm that our test space gives a good stability, in the sense that the solution error is close to the best approximation error.

In this paper we investigate staggered discontinuous Galerkin method for the Helmholtz equation with large wave number on general quadrilateral and polygonal meshes. The method is highly flexible by allowing rough grids such as the trapezoidal grids and highly distorted grids, and at the same time, is numerical flux free. Furthermore, it allows hanging nodes, which can be simply treated as additional vertices. By exploiting a modified duality argument, the stability and convergence can be proved under the condition that \kappa h is sufficiently small, where \kappa h is the wave number and \kappa h is the mesh size. Error estimates for both the scalar and vector variables in \kappa h norm are established. Several numerical experiments are tested to verify our theoretical results and to present the capability of our method for capturing singular solutions.

In this paper, we present a constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM) for parabolic equations with multiscale coefficients, arising from applications in porous media. We will present the construction of CEM-GMsFEM and rigorously analyze its convergence for the parabolic equations. The convergence rate is characterized by the coarse grid size and the eigenvalue decay of local spectral problems, but is independent of the scale length and contrast of the media. The analysis shows that the method has a first order convergence rate with respect to coarse grid size in the energy norm and second order convergence rate with respect to coarse grid size in L^2 norm under some appropriate assumptions. For the temporal discretization, finite difference techniques are used and the convergence analysis of full discrete scheme is given. Moreover, a posteriori error estimator is

In this paper, we present a staggered discontinuous Galerkin immersed boundary method (SDGIBM) for the numerical approximation of fluid-structure interaction. The immersed boundary method is used to model the fluid-structure interaction, while the fluid flow is governed by incompressible Navier-Stokes equations. One advantage of using Galerkin method over the finite difference method with immersed boundary method is that we can avoid approximations of the Dirac Delta function. Another key ingredient of our method is that our solver for incompressible Navier-Stokes equations combines the advantages of discontinuous Galerkin methods and staggered meshes, and results in many good properties, namely local and global conservations and pointwise divergence-free velocity field by a local postprocessing technique. Furthermore, energy stability is improved by a skew-symmetric discretization of the convection term. We will present numerical results to show the performance of the method.

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