No abstract uploaded!

This paper analyzes the random fluctuations obtained by a heterogeneous multi-scale first-order finite element method applied to solve elliptic equations with a random potential. Several multiscale numerical algorithms have been shown to correctly capture the homogenized limit of solutions of elliptic equations with coefficients modeled as stationary and ergodic random fields. Because theoretical results are available in the continuum setting for such equations, we consider here the case of a secondorder elliptic equations with random potential in two dimensions of space. We show that the random fluctuations of such solutions are correctly estimated by the heterogeneous multi-scale algorithm when appropriate fine-scale problems are solved on subsets that cover the whole computational domain. However, when the fine-scale problems are solved over patches that do not cover the entire domain, the random fluctuations may or may not be estimated accurately. In the case of random potentials with short-range interactions, the variance of the random fluctuations is amplified as the inverse of the fraction of the medium covered by the patches. In the case of random potentials with long-range interactions, however, such an amplification does not occur and random fluctuations are correctly captured independent of the (macroscopic) size of the patches. These results are consistent with those obtained in [9] for more general equations in the one-dimensional setting and provide indications on the loss in accuracy that results from using coarser, and hence computationally less intensive, algorithms.

Global geometrical optics method is a new semi-classical approach for the high frequency linear waves proposed by the author in [13]. In this paper, we rederive it in a more concise way. It is shown that the right candidate of solution ansatz for the high frequency wave equations is the extended WKB function, other than the WKB function used in the classical geometrical optics approximation. A new and main contribution of this paper is an interface analysis for the Helmholtz equation when the incident wave is of extended WKB-type. We derive asymptotic expressions for the reflected and/or transmitted propagating waves in the general case. These expressions are valid even when the incident rays include caustic points.

No abstract uploaded!

We propose second order accurate discontinuous Galerkin (DG) schemes which satisfy a strict maximum principle for general nonlinear convection-diffusion equations on unstructured triangular meshes. Motivated by genuinely high order maximum-principle-satisfying DG schemes for hyperbolic conservation laws, we prove that under suitable time step restriction for forward Euler time stepping, for general nonlinear convection-diffusion equations, the same scaling limiter coupled with second order DG methods preserves the physical bounds indicated by the initial condition while maintaining uniform second order accuracy. Similar to the purely convection cases, the limiters are mass conservative and easy to implement. Strong stability preserving (SSP) high order time discretizations will keep the maximum principle. Following the idea in earlier work, we extend the schemes to two-dimensional convection-diffusion equations on triangular meshes. There are no geometric constraints on the mesh such as angle acuteness. Numerical results including incompressible Navier-Stokes equations are presented to validate and demonstrate the effectiveness of the numerical methods.