In our earlier work, we constructed uniformly high order accurate discontinuous Galerkin (DG) and finite volume schemes which preserve positivity of density and pressure for the Euler equations of compressible gas dynamics. In this paper, we present an extension of this framework to construct positivity-preserving high order essentially non-oscillatory (ENO) and weighted essentially non-oscillatory (WENO) finite difference schemes for compressible Euler equations. General equations of state and source terms are also discussed. Numerical tests of the fifth order finite difference WENO scheme are reported to demonstrate the good behavior of such schemes.

No abstract uploaded!

No abstract uploaded!

We quantized the full Einstein equations in a globally hyperbolic spacetime $N=N^{n+1}$, $n\ge 3$, and found solutions of the resulting hyperbolic equation in a fiber bundle $E$ which can be expressed as a product of spatial eigenfunctions (eigendistributions) and temporal eigenfunctions. The spatial eigenfunctions form a basis in an appropriate Hilbert space while the temporal eigenfunctions are solutions to a second order ordinary differential equation in $\R[]_+$. In case $n\ge 17$ and provided the cosmological constant $\Lam$ is negative the temporal eigenfunctions are eigenfunctions of a self-adjoint operator $\hat H_0$ such that the eigenvalues are countable and the eigenfunctions form an orthonormal basis of a Hilbert space.

Shape registration has a wide range of applications in geometric modeling, medical imaging, and computer vision. This paper focuses on the registration of the genus-3 vestibular systems and studies the geometric differences between the normal and Adolescent Idiopathic Scoliosis (AIS) groups. The non-trivial topology of the VS poses great technical challenges to the geometric analysis. To tackle these challenges, we present an effective and practical solution to register the vestibular systems. We first extract six geodesic landmarks for the VS, which are stable, intrinsic, and insensitive to the VS’s resolution and tessellation. Moreover, they are highly consistent regardless of the AIS and normal groups. The detected geodesic landmarks partition the VS into three patches, a topological annulus and two topological disks. For each pair of patches of the AIS subject and the control, we compute a bijective map using the holomorphic 1-form and harmonic map techniques. With a carefully designed boundary condition, the three individual maps can be glued in a seamless manner so that the resulting registration is a homeomorphism with exact landmark matching. Our method is robust, automatic and efficient. It takes only a few seconds on a low-end PC, which significantly outperforms the non-rigid ICP algorithm. We conducted a student’s t-test on the test data. Computational results show that using the mean curvature measure E_H, our method can distinguish the AIS subjects and the normal subjects.