On the upstream mobility scheme for two-phase flow in porous media

When neglecting capillarity, two-phase incompressible flow in porous media is modelled as a scalar nonlinear hyperbolic conservation law. A change in the rock type results in a change of the flux function. Discretizing in one-dimensional with a finite volume method, we investigate two numerical fluxes, an extension of the Godunov flux and the upstream mobility flux, the latter being widely used in hydrogeology and petroleum engineering. Then, in the case of a changing rock type, one can give examples when the upstream mobility flux does not give the right answer.Comment: A preprint to be published in Computational Geoscience

A threshold phenomenon for embeddings of H0mH^m_0 into Orlicz spaces

We consider a sequence of positive smooth critical points of the Adams-Moser-Trudinger embedding of $H^m_0$ into Orlicz spaces. We study its concentration-compactness behavior and show that if the sequence is not precompact, then the liminf of the $H^m_0$-norms of the functions is greater than or equal to a positive geometric constant.Comment: 14 Page

A theory of L1L^1-dissipative solvers for scalar conservation laws with discontinuous flux

We propose a general framework for the study of $L^1$ contractive semigroups of solutions to conservation laws with discontinuous flux. Developing the ideas of a number of preceding works we claim that the whole admissibility issue is reduced to the selection of a family of "elementary solutions", which are certain piecewise constant stationary weak solutions. We refer to such a family as a "germ". It is well known that (CL) admits many different $L^1$ contractive semigroups, some of which reflects different physical applications. We revisit a number of the existing admissibility (or entropy) conditions and identify the germs that underly these conditions. We devote specific attention to the anishing viscosity" germ, which is a way to express the "$\Gamma$-condition" of Diehl. For any given germ, we formulate "germ-based" admissibility conditions in the form of a trace condition on the flux discontinuity line $x=0$ (in the spirit of Vol'pert) and in the form of a family of global entropy inequalities (following Kruzhkov and Carrillo). We characterize those germs that lead to the $L^1$-contraction property for the associated admissible solutions. Our approach offers a streamlined and unifying perspective on many of the known entropy conditions, making it possible to recover earlier uniqueness results under weaker conditions than before, and to provide new results for other less studied problems. Several strategies for proving the existence of admissible solutions are discussed, and existence results are given for fluxes satisfying some additional conditions. These are based on convergence results either for the vanishing viscosity method (with standard viscosity or with specific viscosities "adapted" to the choice of a germ), or for specific germ-adapted finite volume schemes

A Computer-Assisted Uniqueness Proof for a Semilinear Elliptic Boundary Value Problem

A wide variety of articles, starting with the famous paper (Gidas, Ni and Nirenberg in Commun. Math. Phys. 68, 209-243 (1979)) is devoted to the uniqueness question for the semilinear elliptic boundary value problem -{\Delta}u={\lambda}u+u^p in {\Omega}, u>0 in {\Omega}, u=0 on the boundary of {\Omega}, where {\lambda} ranges between 0 and the first Dirichlet Laplacian eigenvalue. So far, this question was settled in the case of {\Omega} being a ball and, for more general domains, in the case {\lambda}=0. In (McKenna et al. in J. Differ. Equ. 247, 2140-2162 (2009)), we proposed a computer-assisted approach to this uniqueness question, which indeed provided a proof in the case {\Omega}=(0,1)x(0,1), and p=2. Due to the high numerical complexity, we were not able in (McKenna et al. in J. Differ. Equ. 247, 2140-2162 (2009)) to treat higher values of p. Here, by a significant reduction of the complexity, we will prove uniqueness for the case p=3

Sharp constants in weighted trace inequalities on Riemannian manifolds

We establish some sharp weighted trace inequalities W^{1,2}(\rho^{1-2\sigma}, M)\hookrightarrow L^{\frac{2n}{n-2\sigma}}(\pa M) on $n+1$ dimensional compact smooth manifolds with smooth boundaries, where $\rho$ is a defining function of $M$ and $\sigma\in (0,1)$. This is stimulated by some recent work on fractional (conformal) Laplacians and related problems in conformal geometry, and also motivated by a conjecture of Aubin.Comment: 34 page

Computing the first eigenpair of the p-Laplacian via inverse iteration of sublinear supersolutions

We introduce an iterative method for computing the first eigenpair $(\lambda_{p},e_{p})$ for the $p$-Laplacian operator with homogeneous Dirichlet data as the limit of $(\mu_{q,}u_{q})$ as $q\rightarrow p^{-}$, where $u_{q}$ is the positive solution of the sublinear Lane-Emden equation $-\Delta_{p}u_{q}=\mu_{q}u_{q}^{q-1}$ with same boundary data. The method is shown to work for any smooth, bounded domain. Solutions to the Lane-Emden problem are obtained through inverse iteration of a super-solution which is derived from the solution to the torsional creep problem. Convergence of $u_{q}$ to $e_{p}$ is in the $C^{1}$-norm and the rate of convergence of $\mu_{q}$ to $\lambda_{p}$ is at least $O(p-q)$. Numerical evidence is presented.Comment: Section 5 was rewritten. Jed Brown was added as autho