Local multiscale methods often construct multiscale basis functions in the offline stage without taking into account input parameters, such as source terms, boundary conditions, and so on. These basis functions are then used in the online stage with a specific input parameter to solve the global problem at a reduced computational cost. Recently, online approaches have been introduced, where multiscale basis functions are adaptively constructed in some regions to reduce the error significantly. In multiscale methods, it is desired to have only 12 iterations to reduce the error to a desired threshold. Using Generalized Multiscale Finite Element Framework [10], it was shown that by choosing sufficient number of offline basis functions, the error reduction can be made independent of physical parameters, such as scales and contrast. In this paper, our goal is to improve this. Using our recently proposed approach [4] and

The authors establish error estimates for recently developed finite-element methods for incompressible viscous flow in domains with no-slip boundary conditions. The methods arise by discretization of a well-posed extended Navier-Stokes dynamics for which pressure is determined from current velocity and force fields. The methods use C1 elements for velocity and C0 elements for pressure. A stability estimate is proved for a related finite-element projection method close to classical time-splitting methods of Orszag, Israeli, DeVille and Karniadakis.

In this paper, we propose a dynamically low-dimensional approximation method to solve a class of time-dependent multiscale stochastic diffusion equations. In Cheng etal. (2013) a dynamically bi-orthogonal (DyBO) method was developed to explore low-dimensional structures of stochastic partial differential equations (SPDEs) and solve them efficiently. However, when the SPDEs have multiscale features in physical space, the original DyBO method becomes expensive. To address this issue, we construct multiscale basis functions within the framework of generalized multiscale finite element method (GMsFEM) for dimension reduction in the physical space. To further improve the accuracy, we also perform online procedure to construct online adaptive basis functions. In the stochastic space, we use the generalized polynomial chaos (gPC) basis functions to represent the stochastic part of the solutions. Numerical

The main purpose of this paper is to analyze the stability and error estimates of the local discontinuous Galerkin (LDG) methods coupled with implicit-explicit (IMEX) time discretization schemes, for solving one-dimensional convection-diffusion equations with a nonlinear convection. Both Runge-Kutta and multi-step IMEX methods are considered. By the aid of the energy method, we show that the IMEX LDG schemes are unconditionally stable for the nonlinear problems, in the sense that the time-step $\dt$ is only required to be upper-bounded by a positive constant which depends on the flow velocity and the diffusion coefficient, but is independent of the mesh size $h$. We also give optimal error estimates for the IMEX LDG schemes, under the same temporal condition, if a monotone numerical flux is adopted for the convection. Numerical experiments are given to verify our main results.

In this article, we propose a numerical method based on sparse Gaussian processes (SGPs) to solve nonlinear partial differential equations (PDEs). The SGP algorithm is based on a Gaussian process (GP) method, which approximates the solution of a PDE with the maximum a posteriori probability estimator of a GP conditioned on the PDE evaluated at a finite number of sample points. The main bottleneck of the GP method lies in the inversion of a covariance matrix, whose cost grows cubically with respect to the size of samples. To improve the scalability of the GP method while retaining desirable accuracy, we draw inspiration from SGP approximations, where inducing points are introduced to summarize the information of samples. More precisely, our SGP method uses a Gaussian prior associated with a low-rank kernel generated by inducing points randomly selected from samples. In the SGP method, the size of the matrix to be inverted is proportional to the number of inducing points, which is much less than the size of the samples. The numerical experiments show that the SGP method using less than half of the uniform samples as inducing points achieves comparable accuracy to the GP method using the same number of uniform samples, which significantly reduces the computational cost. We give the existence proof for the approximation to the solution of a PDE and provide rigorous error analysis.