Robustness of numerical methods has attracted an increasing interest in the community of computational fluid dynamics. One mathematical aspect of robustness for numerical methods is the positivity-preserving property. At high Mach numbers or for flows near vacuum, solving the conservative Euler equations may generate negative density or internal energy numerically, which may lead to nonlinear instability and crash of the code. This difficulty is particularly profound for high order methods, for multi-material flows and for problems with moving meshes, such as the Lagrangian methods. In this paper, we construct both first order and uniformly high order accurate conservative Lagrangian schemes which preserve positivity of physically positive variables such as density and internal energy in the simulation of compressible multi-material flows with general equations of state (EOS). We first develop a positivity-preserving approximate Riemann solver for the Lagrangian scheme solving compressible Euler equations with both ideal and non ideal EOS. Then we design a class of high order positivity-preserving and conservative Lagrangian schemes by using the essentially non-oscillatory (ENO) reconstruction, the strong stability preserving (SSP) high order time discretizations and the positivity-preserving scaling limiter which can be proven to maintain conservation and uniformly high order accuracy and is easy to implement. One-dimensional and two-dimensional numerical tests for the positivity-preserving Lagrangian schemes are provided to demonstrate the effectiveness of these methods.

No abstract uploaded!

No abstract uploaded!

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.

Kronheimer-Mrowka recently proved that the Dehn twist along a 3-sphere in the neck of K3#K3 is not smoothly isotopic to the identity. This provides a new example of self-diffeomorphisms on 4-manifolds that are isotopic to the identity in the topological category but not smoothly so. (The first such examples were given by Ruberman.) In this paper, we use the Pin(2)-equivariant Bauer-Furuta invariant to show that this Dehn twist is not smoothly isotopic to the identity even after a single stabilization (connected summing with the identity map on S2×S2). This gives the first example of exotic phenomena on simply connected smooth 4-manifolds that do not disappear after a single stabilization.