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.

We introduce a quantum trace map for an ideally triangulated hyperbolic knot complement S^3∖K. The map assigns a quantum operator to each element of Kauffmann Skein module of the 3-manifold. The quantum operator lives in a module generated by products of quantized edge parameters of the ideal triangulation modulo some equivalence relations determined by gluing equations. Combining the quantum map with a state-integral model of SL(2,C) Chern-Simons theory, one can define perturbative invariants of knot K in the knot complement whose leading part is determined by its complex hyperbolic length. We then conjecture that the perturbative invariants determine an asymptotic expansion of the Jones polynomial for a link composed of K and K. We propose the explicit quantum trace map for figure-eight knot complement and confirm the length conjecture up to the second order in the asymptotic expansion both numerically and analytically.

Based on the double shockwave approximation procedure, a local Riemann solver for strongly nonlinear equations of state (EOS) such as the Jones-Wilkins-Lee (JWL) EOS is presented, which has suppressed successfully numerical oscillation caused by high-density ratio and high-pressure ratio across the interface between explosion products and air. The real ghost fluid method (RGFM) and the level set method have been used for converting multi-medium flows into pure flows and for implicitly tracking the interface, respectively. A fifth order finite difference weighted essentially non-oscillatory (WENO) scheme and a third order TVD Runge-Kutta method are utilized for the spatial discretization and the time advance, respectively. An enclosed-type MPI-based parallel methodology for the RGFM procedure on a uniform structured mesh is presented to realize the parallelization of three-dimensional {\color{red}(3D)} air explosion. The overall process of 3D air explosion of both TNT and aluminized explosives has been successfully simulated.

We study BPS spectra of D-branes on local Calabi-Yau threefolds O(−p)⊕O(p−2)→P^1 with p=0,1, corresponding to C^3/Z_2 and the resolved conifold. Nonabelianization for exponential networks is applied to compute directly unframed BPS indices counting states with D2 and D0 brane charges. Known results on these BPS spectra are correctly reproduced by computing new types of BPS invariants of 3d-5d BPS states, encoded by nonabelianization, through their wall-crossing. We also develop the notion of exponential BPS graphs for the simplest toric examples, and show that they encode both the quiver and the potential associated to the Calabi-Yau via geometric engineering.

Consider an age-dependent, single-species branching process defined by a progeny number distribution, and a lifetime distribution associated with each independent particle. In this paper, we focus on the associated inverse problem where one wishes to formally solve for the progeny number distribution or the lifetime distribution that defines the Bellman-Harris branching process. We derive results for the existence and uniqueness (the identifiability) of these two distributions given one of two types of information: the extinction time probability of the entire process (extinction time distribution), or the distribution of the total number of particles at one fixed time. We demonstrate that perfect knowledge of the distribution of extinction times allows us to formally determine either the progeny number distribution or the lifetime distribution. Furthermore, we show that these constructions are unique. We then consider “data” consisting of a perfectly known total number distribution given at one specific time. For a process with known progeny number distribution and exponentially distributed lifetimes, we show that the rate parameter is identifiable. For general lifetime distributions, we also show that the progeny distribution is globally unique. Our results are presented through four theorems, each describing the constructions in the four distinct cases.