A natural metric in 2-manifold surfaces is to use geodesic distance. If a 2-manifold surface is represented by a triangle mesh T , the geodesic metric on T can be computed exactly using computational geometry methods. Previous work for establishing the geodesic metric on T only supports using half-edge data structures; i.e., each edge e in T is split into two halves (he1, he2) and each half-edge corresponds to one of two faces incident to e. In this paper, we prove that the exact-geodesic structures on two halfedges of e can be merged into one structure associated with e. Four merits are achieved based on the properties which are studied in this paper: (1) Existing CAD systems that use edge-based data structures can directly add the geodesic distance function without changing the kernel to a half-edge data structure; (2) To find the geodesic path from inquiry points to the source, the MMP algorithm can be run in an onthe- fly fashion such that the inquiry points are covered by correct wedges; (3) The MMP algorithm is sped up by pruning unnecessary wedges during the wedge propagation process; (4) The storage of the MMP algorithm is reduced since fewer wedges need to be stored in an edge-based data structure. Experimental results show that when compared to the classic half-edge data structure, the edge-based implementation of the MMP algorithm reduces running time by 44% and storage by 29% on average.

In this paper, we analyze the Lax-Wendroff discontinuous Galerkin (LWDG) method for solving linear conservation laws. The method was originally proposed by Guo et al., where they applied local discontinuous Galerkin (LDG) techniques to approximate high order spatial derivatives in the Lax-Wendroff time discretization. We show that, under the standard CFL condition $\tau\leq \lambda h$ (where $\tau$ and $h$ are the time step and the maximum element length respectively and $\lambda>0$ is a constant) and uniform or non-increasing time steps, the second order schemes with piecewise linear elements and the third order schemes with arbitrary piecewise polynomial elements are stable in the $L^2$ norm. The specific type of stability may differ with different choices of numerical fluxes. Our stability analysis includes multidimensional problems with divergence-free coefficients. Besides solving the equation itself, the LWDG method also gives approximations to its time derivative simultaneously. We obtain optimal error estimates for both the solution $u$ and its first order time derivative $u_t$ in one dimension, and numerical examples are given to validate our analysis.

We further investigate the high order positivity-preserving discontinuous Galerkin (DG) methods for linear hyperbolic and radiative transfer equations developed in \cite{yuan}. The DG methods in \cite{yuan} can maintain positivity and high order accuracy, but they rely both on the scaling limiter in \cite{zhang1} and a rotational limiter, the latter may alter cell averages of the unmodulated DG scheme, thereby affecting conservation. Even though a Lax-Wendroff type theorem is proved in \cite{yuan}, guaranteeing convergence to weak solutions with correct shock speed when such rotational limiter is applied, it would still be desirable if a conservative DG method without changing the cell averages can be obtained which has both high order accuracy and positivity-preserving capability. In this paper, we develop and analyze such a DG method for both linear hyperbolic equations and radiative transfer equations. In the one-dimensional case, the method uses traditional DG space $P^k$ of piecewise polynomials of degree at most $k$. A key result is proved that the unmodulated DG solver in this case can maintain positivity of the cell average if the inflow boundary value and the source term are both positive, therefore the positivity-preserving framework in \cite{zhang1} can be used to obtain a high order conservative positivity-preserving DG scheme. Unfortunately, in two-dimensions this is no longer the case. We show that the unmodulated DG solver based either on $P^k$ or $Q^k$ spaces (piecewise $k$-th degree polynomials or piecewise tensor-product $k$-th degree polynomials) could generate negative cell averages. We augment the DG space with additional functions so that the positivity of cell averages from the unmodulated DG solver can be restored, thereby leading to high order conservative positivity-preserving DG scheme based on these augmented DG spaces following the framework in \cite{zhang1}. Computational results are provided to demonstrate the good performance of our DG schemes.

The inflow problem of full compressible NavierStokes equations is considered on the half-line (0,+\infty). First, we give the existence (or nonexistence) of the boundary layer solution to the inflow problem when the right end state (0,+\infty) belongs to the subsonic, transonic, and supersonic regions, respectively. Then the asymptotic stability of not only the single contact wave but also the superposition of the subsonic boundary layer solution, the contact wave, and the rarefaction wave to the inflow problem are investigated under some smallness conditions. Note that the amplitude of the rarefaction wave is not necessarily small. The proofs are given by the elementary energy method.

Let $f(x)=\sum_{n=0}^{\infty}\frac{1}{n!}q^{n(n-1)/2}x^n$ ($0<q<1$) be the deformed exponential function. It is known that the zeros of $f(x)$ are real and form a negative decreasing sequence $(x_k)$ ($k\ge 1$). We investigate the complete asymptotic expansion for $x_{k}$ and prove that for any $n\ge1$, as $k\to \infty$, \begin{align*} x_k=-kq^{1-k}\Big(1+\sum_{i=1}^{n}C_i(q)k^{-1-i}+o(k^{-1-n})\Big), \end{align*} where $C_i(q)$ are some $q$ series which can be determined recursively. We show that each $C_{i}(q)\in \mathbb{Q}[A_0,A_1,A_2]$, where $A_{i}=\sum_{m=1}^{\infty}m^i\sigma(m)q^m$ and $\sigma(m)$ denotes the sum of positive divisors of $m$. When writing $C_{i}$ as a polynomial in $A_0, A_1$ and $A_2$, we find explicit formulas for the coefficients of the linear terms by using Bernoulli numbers. Moreover, we also prove that $C_{i}(q)\in \mathbb{Q}[E_2,E_4,E_6]$, where $E_2$, $E_4$ and $E_6$ are the classical Eisenstein series of weight 2, 4 and 6, respectively.