No abstract uploaded!

Wormhole propagation, arising in petroleum engineering, is used to describe the distribution of acid and the increase of porosity in carbonate reservoir under dissolution of injected acid. The important physical features of porosity and acid concentration include their boundedness between $0$ and $1$, as well as the monotone increasing for porosity. How to keep these properties in the simulation is crucial to the robustness of the numerical algorithm. In this paper, we propose high-order bound-preserving discontinuous Galerkin methods to keep these important physical properties. The main technique is to introduce a new variable $r$ to replace the original acid concentration and use a consistent flux pair to deduce a ghost equation such that the positive-preserving technique can be applied on both original and deduced equations. A high-order slope limiter is used to keep a polynomial upper bound which changes over time for $r$. Moreover, the high-order accuracy is attained by the flux limiter. Numerical examples are given to demonstrate the high-order accuracy and bound-preserving property of the numerical technique.

Surface segmentation is a fundamental problem in computer graphics. It has various applications such as metamorphosis, surface matching, surface compression, 3D shape retrieval, texture mapping, etc. All orientable surfaces are Riemann surfaces, and admit conformal structures. This paper introduces a novel surface segmentation algorithm based on its conformal structure. Each segment can be conformally mapped to a planar rectangle, and the transition maps are planar translations. The segmentation is intrinsic to the surface, independent of the embedding, and consistent for surfaces with similar geometries. By using segmentation based on conformal structure, the mapping between surfaces with arbitrary topologies can be constructed explicitly. The method is rigorous, efficient and automatic. The segmentation can be applied to surface morphing, construct conformal geometry image, convert mesh to Spline

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.

No abstract uploaded!