No abstract uploaded!

We develop a method searching for quasi-conformal maps from point-cloud surfaces to Euclidean plane. The maps can be got by directly solving Beltrami equation with proper boundary condition or solving optimization problem of the variational formu- lation. The point integral methd (PIM) is used for discretization on point clouds. Numerical experiments suggest that the method converges as the density of points in- creases.

Positivity-preserving discontinuous Galerkin (DG) methods for solving hyperbolic conservation laws have been extensively studied in the last several years. But nearly all the developed schemes are coupled with explicit time discretizations. Explicit discretizations suffer from the constraint for the Courant-Friedrichs-Levis (CFL) number. This makes explicit methods impractical for problems involving unstructured and extremely varying meshes or long-time simulations. Instead, implicit DG schemes are often popular in practice, especially in the computational fluid dynamics (CFD) community. In this paper we develop a high-order positivity-preserving DG method with the backward Euler time discretization for conservation laws. We focus on one spatial dimension, however the result easily generalizes to multidimensional tensor product meshes and polynomial spaces. This work is based on a generalization of the positivity-preserving limiters in (X. Zhang and C.-W. Shu, Journal of Computational Physics, 229 (2010), pp.~3091--3120) and (X. Zhang and C.-W. Shu, Journal of Computational Physics, 229 (2010), pp.~8918--8934) to implicit time discretizations. Both the analysis and numerical experiments indicate that a lower bound for the CFL number is required to obtain the positivity-preserving property. The proposed scheme not only preserves the positivity of the numerical approximation without compromising the designed high-order accuracy, but also helps accelerate the convergence towards the steady-state solution and add robustness to the nonlinear solver. Numerical experiments are provided to support these conclusions.

Accurate simulation of seismic waves is of critical importance in a variety of geophysical applications. Based on recent works on staggered discontinuous Galerkin methods, we have developed a new method for the simulations of seismic waves, which has energy conservation and extremely low grid dispersion, so that it naturally provided accurate numerical simulations of wave propagation useful for geophysical applications and was a generalization of classical staggered-grid finite-difference methods. Moreover, it could handle with ease irregular surface topography and discontinuities in the subsurface models. Our new method discretized the velocity and the stress tensor on this staggered grid, with continuity imposed on different parts of the mesh. The symmetry of the stress tensor was enforced by the Lagrange multiplier technique. The resulting method was an explicit scheme, requiring the solutions of a

No abstract uploaded!