In this paper, we discuss a general multiscale model reduction framework based on multiscale finite element methods. We give a brief overview of related multiscale methods. Due to page limitations, the overview focuses on a few related methods and is not intended to be comprehensive. We present a general adaptive multiscale model reduction framework, the Generalized Multiscale Finite Element Method. Besides the method's basic outline, we discuss some important ingredients needed for the method's success. We also discuss several applications. The proposed method allows performing local model reduction in the presence of high contrast and no scale separation.

No abstract uploaded!

High order path-conservative schemes have been developed for solving nonconservative hyperbolic systems. Recently, it has been observed that this approach may have some computational issues and shortcomings. In this paper, a modification to the high order path-conservative scheme in is proposed to improve its computational performance and to overcome some of the shortcomings. This modification is based on the high order finite volume WENO scheme with subcell resolution and it uses an exact Riemann solver to catch the right paths at the discontinuities. An application to one-dimensional compressible two-medium flows of nonconservative or primitive Euler equations is carried out to show the effectiveness of this new approach.

In this paper, we present L2 and negative-order norm estimates for the local discontinuous Galerkin (LDG) method solving variable coefficient Schrödinger equations. For these special solutions the LDG method exhibits ‘‘hidden accuracy", and we are able to extract it through the use of a convolution kernel that is composed of a linear combination of B-splines. This technical was initially introduced by Cockburn, Luskin, Shu, and Süli for linear hyperbolic equations and extended by Ryan et al. as a smoothness-increasing accuracy-conserving (SIAC) filter. We demonstrate that it is possible to extend the SIAC filter on Schrödinger equations. When polynomials of degree k are used, we can prove theoretically the LDG method solutions are of order k + 1, whereas the post-processed solutions that convolution with the SIAC filter are of order at least 2k. Additionally, we present numerical results to confirm that the accuracy of LDG solutions can be improved from O(hk+1) to at least O(h2k+1) for Schrödinger equations by using alternating numerical fluxes.

In this paper we derive new fractional error bounds for polynomial systems with exponents explicitly determined by the dimension of the underlying space and the number/degree of the involved polynomials. Our major result extends the existing error bounds from the system involving only a single polynomial to a general polynomial system and do not require any regularity assumptions. In this way we resolve, in particular, some open questions posed in the literature. The developed techniques are largely based on variational analysis and generalized differentiation, which allow us to establish, e.g., a nonsmooth extension of the seminal Łojasiewicz’s gradient inequality to maxima of polynomials with explicitly determined exponents. Our major applications concern quantitative Hölderian stability of solution maps for parameterized polynomial optimization problems and nonlinear complementarity systems with polynomial data as well as high-order semismooth properties of the eigenvalues of symmetric tensors.