We propose an extension of the computational fluid mechanics approach to the Monge-Kantorovich mass transfer problem, which was developed by Benamou-Brenier. Our extension allows optimal transfer of unnormalized and unequal masses. We obtain a one-parameter family of simple modifications of the formulation in [4]. This leads us to a new Monge-Ampere type equation and a new Kantorovich duality formula. These can be solved efficiently by, for example, the Chambolle-Pock primal-dual algorithm. This solution to the extended mass transfer problem gives us a simple metric for computing the distance between two unnormalized densities. The L1 version of this metric was shown in [23] (which is a precursor of our work here) to have desirable properties.

No abstract uploaded!

We propose an efficient eigensolver for computing densely distributed spectra of the two-dimensional transmission eigenvalue problem (TEP), which is derived from Maxwell's equations with Tellegen media and the transverse magnetic mode. The governing equations, when discretized by the standard piecewise linear finite element method, give rise to a large-scale quadratic eigenvalue problem (QEP). Our numerical simulation shows that half of the positive eigenvalues of the QEP are densely distributed in some interval near the origin. The quadratic JacobiDavidson method with a so-called non-equivalence deflation technique is proposed to compute the dense spectrum of the QEP. Extensive numerical simulations show that our proposed method processes the convergence efficiently, even when it needs to compute more than 5000 desired eigenpairs. Numerical results also illustrate that the computed

In this paper, we extend the high order finite-difference method with subcell resolution (SR) for two-species stiff one-reaction models to multispecies and multireaction chemical reactive flows, which are significantly more difficult because of the multiple scales generated by different reactions. For reaction problems, when the reaction time scale is very small, the reaction zone scale is also small and the governing equations become very stiff. Wrong propagation speed of discontinuity may occur due to the underresolved numerical solution in both space and time. The present SR method for reactive Euler system is a fractional step method. In the convection step, any high order shock-capturing method can be used. In the reaction step, an ODE solver is applied but with certain computed flow variables in the shock region modified by the Harten subcell resolution idea. Several numerical examples of multispecies and multireaction reactive flows are performed in both one and two dimensions. Studies demonstrate that the SR method can capture the correct propagation speed of discontinuities in very coarse meshes.

Wormhole propagation, arising in petroleum engineering, is used to describe the distribution of acid and the increase of porosity in carbonate reservoir under dissolution of injected acid. The important physical features of porosity and acid concentration include their boundedness between $0$ and $1$, as well as the monotone increasing for porosity. How to keep these properties in the simulation is crucial to the robustness of the numerical algorithm. In this paper, we propose high-order bound-preserving discontinuous Galerkin methods to keep these important physical properties. The main technique is to introduce a new variable $r$ to replace the original acid concentration and use a consistent flux pair to deduce a ghost equation such that the positive-preserving technique can be applied on both original and deduced equations. A high-order slope limiter is used to keep a polynomial upper bound which changes over time for $r$. Moreover, the high-order accuracy is attained by the flux limiter. Numerical examples are given to demonstrate the high-order accuracy and bound-preserving property of the numerical technique.