Simulation of electromagnetic wave propagation in metamaterials leads to more complicated time domain Maxwell's equations than the standard Maxwell's equations in free space. In this paper, we develop and analyze a non-dissipative discontinuous Galerkin (DG) method for solving the Maxwell's equations in Drude metamaterials. Previous discontinuous Galerkin methods in the literature for electromagnetic wave propagation in metamaterials were either non-dissipative but sub-optimal, or dissipative and optimal. Our method uses a different and simple choice of numerical fluxes, achieving provable non-dissipative stability and optimal error estimates simultaneously. We prove the stability and optimal error estimates for both semi- and fully discrete DG schemes, with the leap-frog time discretization for the fully discrete case. Numerical results are given to demonstrate that the DG method can solve metamaterial Maxwell's equations effectively.

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 consider an online basis enrichment mixed generalized multiscale method with oversampling, for solving flow problems in highly heterogeneous porous media. This is an exten-sion of the online mixed generalized multiscale method [6]. The multiscale online basis functions are computed by solving a Neumann problem in an over-sampled domain, instead of a standard neighborhood of a coarse face. We are motivated by the restricted domain decomposition method. Extensive numerical experiments are presented to demonstrate the performance of our methods for both steady-state flow, and two-phase flow and transport problems.

The well-posedness of the three space dimensional Prandtl equations is studied under some constraint on its flow structure. Together with the instability result given in [28], it gives an almost necessary and sufficient structural condition for the stability of the three-dimensional Prandtl equations. It reveals that the classical Burgers equation plays an important role in determining this type of flow with special structure, that avoids the appearance of the secondary flow, an unstabilizing factor in the three-dimensional Prandtl boundary layers. And the sufficiency of the monotonicity condition on the tangential velocity field for the existence of solutions to the Prandtl boundary layer equations is illustrated in the three-dimensional setting. Moreover, it is shown that this structured flow is linearly stable for any smooth three-dimensional perturbation.

Characterization of homogeneous polynomials with isolated critical point at the origin follows from a study of complex geometry. Yau previously proposed a Numerical Characterization Conjecture. A step forward in solving this Conjecture, the Granville-Lin-Yau Conjecture was formulated, with a sharp estimate that counts the number of positive integral points in ndimensional (n≥3) real right-angled simplices with vertices whose distance to the origin are at least n-1. The estimate was proven for n≤6 but has a counterexample for n = 7. In this project we come up with an idea of forming a new sharp estimate conjecture where we need the distances of the vertices to be n. We have proved this new sharp estimate conjecture for n≤7 and are in the process of proving the general n case.