• ### Optimality of integrability estimates for advection-diffusion equations(1702.00321)

Feb. 1, 2017 math.AP
We discuss $L^p$ integrability estimates for the solution $u$ of the advection-diffusion equation $\partial_t u + \mathrm{div} (bu) = \Delta u$, where the velocity field $b \in L^r_t L^q_x$. We first summarize some classical results proving such estimates for certain ranges of the exponents $r$ and $q$. Afterwards we prove the optimality of such ranges by means of new original examples.
• ### On the structure of $L^\infty$-entropy solutions to scalar conservation laws in one-space dimension(1608.02811)

Aug. 9, 2016 math.AP
We prove that if $u$ is the entropy solution to a scalar conservation law in one space dimension, then the entropy dissipation is a measure concentrated on countably many Lipschitz curves. This result is a consequence of a detailed analysis of the structure of the characteristics. \\ In particular the characteristic curves are segments outside a countably 1-rectifiable set and the left and right traces of the solution exist in a $C^0$-sense up to the degeneracy due to the segments where $f"=0$. We prove also that the initial data is taken in a suitably strong sense and we give some counterexamples which show that these results are sharp.
• ### Renormalization for autonomous nearly incompressible BV vector fields in 2D(1412.6387)

July 11, 2019 math.AP
Given a bounded autonomous vector field $b \colon \mathbb R^d \to \mathbb R^d$, we study the uniqueness of bounded solutions to the initial value problem for the related transport equation \begin{equation*} \partial_t u + b \cdot \nabla u= 0. \end{equation*} We are interested in the case where $b$ is of class BV and it is nearly incompressible. Assuming that the ambient space has dimension $d=2$, we prove uniqueness of weak solutions to the transport equation. The starting point of the present work is the result which has been obtained in \cite{BG} (where the \emph{steady} case is treated). Our proof is based on splitting the equation onto a suitable partition of the plane: this technique was introduced in \cite{ABC1}, using the results on the structure of level sets of Lipschitz maps obtained in \cite{ABC2}. Furthermore, in order to construct the partition, we use Ambrosio's superposition principle \cite{ambrosiobv}.
• ### Quadratic interaction functional for general systems of conservation laws(1412.6408)

Dec. 19, 2014 math.AP
For the Glimm scheme approximation $u_\epsilon$ to the solution of the system of conservation laws in one space dimension \begin{equation*} u_t + f(u)_x = 0, \qquad u(0,x) = u_0(x) \in \mathbb R^n, \end{equation*} with initial data $u_0$ with small total variation, we prove a quadratic (w.r.t. $\mathrm{TV}(u_0)$) interaction estimate, which has been used in the literature for stability and convergence results. No assumptions on the structure of the flux $f$ are made (apart smoothness), and this estimate is the natural extension of the Glimm type interaction estimate for genuinely nonlinear systems. More precisely we obtain the following results: \newline - a new analysis of the interaction estimates of simple waves; \newline - a Lagrangian representation of the derivative of the solution, i.e. a map $\mathtt x(t,w)$ which follows the trajectory of each wave $w$ from its creation to its cancellation; \newline - the introduction of the characteristic interval and partition for couples of waves, representing the common history of the two waves; \item a new functional $\mathfrak Q$ controlling the variation in speed of the waves w.r.t. time. \newline This last functional is the natural extension of the Glimm functional for genuinely nonlinear systems. The main result is that the distribution $D_{tt} \mathtt x(t,w)$ is a measure with total mass $\leq \mathrm{const} \mathrm{TV}(u_0)^2$.
• ### SBV regularity for Hamilton-Jacobi equations in $\mathbb R^n$(1002.4087)

Dec. 16, 2014 math.AP
In this paper we study the regularity of viscosity solutions to the following Hamilton-Jacobi equations $$\partial_t u + H(D_{x} u)=0 \qquad \textrm{in} \Omega\subset \mathbb R\times \mathbb R^{n} .$$ In particular, under the assumption that the Hamiltonian $H\in C^2(\mathbb R^n)$ is uniformly convex, we prove that $D_{x}u$ and $\partial_t u$ belong to the class $SBV_{loc}(\Omega)$.
• ### The decomposition of optimal transportation problems with convex cost(1409.0515)

Sept. 1, 2014 math.CA
Given a positive l.s.c. convex function $\mathtt c : \mathbb R^d \to \mathbb R^d$ and an optimal transference plane $\underline{\pi}$ for the transportation problem \begin{equation*} \int \mathtt c(x'-x) \pi(dxdx'), \end{equation*} we show how the results of \cite{biadan} on the existence of a \emph{Sudakov decomposition} for norm cost $\mathtt c= |\cdot|$ can be extended to this case. More precisely, we prove that there exists a partition of $\mathbb R^d$ into a family of disjoint sets $\{S^h_\mathfrak a\}_{h,\mathfrak a}$ together with the projection $\{O^h_\mathfrak a\}_{h,\mathfrak a}$ on $\mathbb R^d$ of proper extremal faces of $\mathrm{epi}\, \mathtt c$, $h = 0,\dots,d$ and $\mathfrak a \in \mathfrak A^h \subset \mathbb R^{d-h}$, such that - $S^h_\mathfrak a$ is relatively open in its affine span, and has affine dimension $h$; \item $O^h_\mathfrak a$ has affine dimension $h$ and is parallel to $S^h_\mathfrak a$; - $\mathcal L^d(\mathbb R^d \setminus \cup_{h,\mathfrak a} S^h_\mathfrak a) = 0$, and the disintegration of $\mathcal L^d$, $\mathcal L^d = \sum_h \int \xi^h_\mathfrak a \eta^h(d\mathfrak a)$, w.r.t. $S^h_\mathfrak a$ has conditional probabilities $\xi^h_\mathfrak a \ll \mathcal H^h \llcorner_{S^h_\mathfrak a}$; - the sets $S^h_\mathfrak a$ are essentially cyclically connected and cannot be further decomposed. \end{list} The last point is used to prove the existence of an optimal transport map. The main idea is to recast the problem in $(t,x) \in [0,\infty] \times \mathbb R^d$ with an $1$-homogeneous norm $\bar{\mathtt c}(t,x) := t \mathtt c(- \frac{x}{t})$ and to extend the regularity estimates of \cite{biadan} to this case.
• ### Structure of entropy solutions to general scalar conservation laws in one space dimension(1409.0521)

Sept. 1, 2014 math.AP
In this paper, we show that the entropy solution of a scalar conservation law is - continuous outside a $1$-rectifiable set $\Xi$, - up to a $\mathcal H^1$ negligible set, for each point $(\bar t,\bar x) \in \Xi$ there exists two regions where $u$ is left/right continuous in $(\bar t,\bar x)$. We provide examples showing that these estimates are nearly optimal. In order to achieve these regularity results, we extend the wave representation of the wavefront approximate solutions to entropy solution. This representation can the interpreted as some sort of Lagrangian representation of the solution to the nonlinear scalar PDE, and implies a fine structure on the level sets of the entropy solution.
• ### Steady nearly incompressible vector fields in 2D: chain rule and renormalization(1408.2932)

Aug. 13, 2014 math.AP
Given bounded vector field $b : \mathbb R^d \to \mathbb R^d$, scalar field $u : \mathbb R^d \to \mathbb R$ and a smooth function $\beta : \mathbb R \to \mathbb R$ we study the characterization of the distribution $\mathrm{div}(\beta(u)b)$ in terms of $\mathrm{div}\, b$ and $\mathrm{div}(u b)$. In the case of $BV$ vector fields $b$ (and under some further assumptions) such characterization was obtained by L. Ambrosio, C. De Lellis and J. Mal\'y, up to an error term which is a measure concentrated on so-called \emph{tangential set} of $b$. We answer some questions posed in their paper concerning the properties of this term. In particular we construct a nearly incompressible $BV$ vector field $b$ and a bounded function $u$ for which this term is nonzero. For steady nearly incompressible vector fields $b$ (and under some further assumptions) in case when $d=2$ we provide complete characterization of $\mathrm{div}(\beta(u) b)$ in terms of $\mathrm{div}\, b$ and $\mathrm{div}(u b)$. Our approach relies on the structure of level sets of Lipschitz functions on $\mathrm R^2$ obtained by G. Alberti, S. Bianchini and G. Crippa. Extending our technique we obtain new sufficient conditions when any bounded weak solution $u$ of $\partial_t u + b \cdot \nabla u=0$ is \emph{renormalized}, i.e. also solves $\partial_t \beta(u) + b \cdot \nabla \beta(u)=0$ for any smooth function $\beta : \mathbb R \to \mathbb R$. As a consequence we obtain new uniqueness result for this equation.
• ### On Sudakov's type decomposition of transference plans with norm costs(1311.1918)

Jan. 7, 2014 math.CA
We consider the original strategy proposed by Sudakov for solving the Monge transportation problem with norm cost $|\cdot|_{D^*}$ $\min \bigg\{\int |\mathtt T(x) - x|_{D^*} d\mu(x), \ \mathtt T : \mathbb R^d \to \mathbb R^d, \ \nu = \mathtt T_\# \mu \bigg\},$ with $\mu$, $\nu$ probability measures in $\mathbb R^d$ and $\mu$ absolutely continuous w.r.t. $\mathcal L^d$. The key idea in this approach is to decompose (via disintegration of measures) the Kantorovich optimal transportation problem into a family of transportation problems in $Z_\mathfrak a\times\mathbb R^d$, where $\{Z_\mathfrak a\}_{\mathfrak a\in\mathfrak A} \subset \mathbb R^d$ are disjoint regions such that the construction of an optimal map $\mathtt T_\mathfrak a : Z_\mathfrak a \to \mathbb R^d$ is simpler than in the original problem, and then to obtain $\mathtt T$ by piecing together the maps $\mathtt T_\mathfrak a$. In this paper we show how the original idea of Sudakov can be successfully implemented. The results yield a complete characterization of the Kantorovich optimal transportation problem, whose straightforward corollary is the solution of the Monge problem in each set $Z_\mathfrak a$ and then in $\mathbb R^d$. The strategy is sufficiently powerful to be applied to other optimal transportation problems.
• ### Quadratic interaction functional for systems of conservation laws: a case study(1311.2722)

Nov. 12, 2013 math.AP
We prove a quadratic interaction estimate for wavefront approximate solutions to the triangular system of conservation laws $\begin{cases} u_t + \tilde f(u,v)_x = 0, \cr v_t - v_x = 0. \end{cases}$ This quadratic estimate has been used in the literature to prove the convergence rate of the Glimm scheme. Our aim is to extend the analysis, done for scalar conservation laws, in the presence of transversal interactions among wavefronts of different families. The proof is based on the introduction of a quadratic functional $\mathfrak Q(t)$, decreasing at every interaction, and such that its total variation in time is bounded. %cancellations it variation is controlled by the total variation growths at most of the total variation of the solution multiplied by the amount of cancellation. The study of this particular system is a key step in the proof of the quadratic interaction estimate for general systems: it requires a deep analysis of the wave structure of the solution $(u(t,x),v(t,x))$ and the reconstruction of the past history of each wavefront involved in an interaction.
• ### On a quadratic functional for scalar conservation laws(1311.2929)

Nov. 12, 2013 math.AP
We prove a quadratic interaction estimate for approximate solutions to scalar conservation laws obtained by the wavefront tracking approximation or the Glimm scheme. This quadratic estimate has been used in the literature to prove the convergence rate of the Glimm scheme. The proof is based on the introduction of a quadratic functional $\mathfrak Q(t)$, decreasing at every interaction, and such that its total variation in time is bounded. Differently from other interaction potentials present in the literature, the form of this functional is the natural extension of the original Glimm functional, and coincides with it in the genuinely nonlinear case.
• ### Global structure of admissible BV solutions to piecewise genuinely nonlinear, strictly hyperbolic conservation laws in one space dimension(1211.3526)

Dec. 2, 2012 math.AP
The paper gives an accurate description of the qualitative structure of an admissible BV solution to a strictly hyperbolic, piecewise genuinely nonlinear system of conservation laws. We prove that there are a countable set $\Theta$ which contains all interaction points and a family of countably many Lipschitz curves $\T$ such that outside $\T\cup \Theta$ $u$ is continuous, and along the curves in $\T$, u has left and right limit except for points in $\Theta$. This extends the corresponding structural result in \cite{BL,Liu1} for admissible solutions. The proof is based on approximate wave-front tracking solutions and a proper selection of discontinuity curves in the approximate solutions, which converge to curves covering the discontinuities in the exact solution $u$.
• ### SBV-like regularity for general hyperbolic systems of conservation laws(1202.2680)

Feb. 13, 2012 math.AP
• ### The Monge Problem for distance cost in geodesic spaces(1103.2796)

March 14, 2011 math.AP, math.PR
We address the Monge problem in metric spaces with a geodesic distance: (X, d) is a Polish space and dL is a geodesic Borel distance which makes (X,dL) a non branching geodesic space. We show that under the assumption that geodesics are d-continuous and locally compact, we can reduce the transport problem to 1-dimensional transport problems along geodesics. We introduce two assumptions on the transport problem {\pi} which imply that the conditional probabilities of the first marginal on each geodesic are continuous or absolutely continuous w.r.t. the 1- dimensional Hausdorff distance induced by dL. It is known that this regularity is sufficient for the construction of a transport map. We study also the dynamics of transport along the geodesic, the stability of our conditions and show that in this setting dL-cyclical monotonicity is not sufficient for optimality.
• ### A connection between viscous profiles and singular ODEs(0801.4532)

April 2, 2009 math.AP, math.CA
We deal with the viscous profiles for a class of mixed hyperbolic-parabolic systems. We focus, in particular, on the case of the compressible Navier Stokes equation in one space variable written in Eulerian coordinates. We describe the link between these profiles and a singular ordinary differential equation in the form $$dV / dt = F(V) / z (V) .$$ Here $V \in R^d$ and the function F takes values into $R^d$ and is smooth. The real valued function z is as well regular: the equation is singular in the sense that z (V) can attain the value 0.
• ### Invariant manifolds for a singular ordinary differential equation(0801.4425)

March 25, 2009 math.AP, math.CA
We study the singular ordinary differential equation $$\frac{d U}{d t} = f (U) / z (U) + g (U),$$ where $U \in R^N$, the functions $f \in R^N$ and $g \in R^N$ are of class $C^2$ and $z$ is a real valued $C^2$ function. The equation is singular in the sense that $z (U)$ can attain the value 0. We focus on the solutions of the singular ODE that belong to a small neighborhood of a point $\bar U$ such that $f (\bar U) = g (\bar U) = \vec 0$, $z (\bar U) =0$. We investigate the existence of manifolds that are locally invariant for the singular ODE and that contain orbits with a suitable prescribed asymptotic behaviour. Under suitable hypotheses on the set $\{U: z (U) = 0 \}$, we extend to the case of the singular ODE the definitions of center manifold, center stable manifold and of uniformly stable manifold. An application of our analysis concerns the study of the viscous profiles with small total variation for a class of mixed hyperbolic-parabolic systems in one space variable. Such a class includes the compressible Navier Stokes equation.
• ### Invariant Manifolds for Viscous Profiles of a Class of Mixed Hyperbolic-Parabolic Systems(0812.1184)

Dec. 5, 2008 math.AP
We are concerned with viscous profiles (travelling waves and steady solutions) for mixed hyperbolic-parabolic systems in one space variable. For a class of systems including the compressible Navier Stokes equation, these profiles satisfy a singular ordinary differential equation in the form \label{e:ab} dU / dt = F(U)/ z (U) . Here U takes values in $R^d$ and $F: R^d \to R^d$ is a regular function. The real valued function $z (U)$ is as well regular, but the equation is singular because $z (U)$ can attain the value 0. We focus on a small enough neighbourhood of a point $\bar U$ satisfying $F(\bar U) = \vec 0$, $z (\bar U) =0$. From the point of view of the applications to the study of hyperbolic-parabolic systems this means restricting to systems with small total variation. We discuss how to extend the notions of center manifold and of uniformly stable manifold. Also, we give conditions ensuring that if $z (U) \neq 0$ at $t=0$ then $z (U) \neq 0$ at every $t$. We provide an example showing that if $z(U)$ becomes zero in finite time then in general the solution $U$ of equation \eqref{e:ab} is not continuously differentiable.
• ### Asymptotic behavior of smooth solutions for partially dissipative hyperbolic systems with a convex entropy(0805.3614)

May 23, 2008 math.AP
We study the asymptotic time behavior of global smooth solutions to general entropy dissipative hyperbolic systems of balance law in m space dimensions, under the Shizuta-Kawashima condition. We show that these solutions approach constant equilibrium state in the Lp-norm at a rate O(t^(-m/2(1-1/p))), as t tends to $\infty$, for p in [min (m,2),+ \infty]. Moreover, we can show that we can approximate, with a faster order of convergence, theconservative part of the solution in terms of the linearized hyperbolic operator for m >= 2, and by a parabolic equation in the spirit of Chapman-Enskog expansion. The main tool is given by a detailed analysis of the Green function for the linearized problem.
• ### Vanishing Viscosity Solutions of Nonlinear Hyperbolic Systems(math/0111321)

Nov. 30, 2001 math.AP
We consider the Cauchy problem for a strictly hyperbolic, $n\times n$ system in one space dimension: $u_t+A(u)u_x=0$, assuming that the initial data has small total variation. We show that the solutions of the viscous approximations $u_t+A(u)u_x=\ve u_{xx}$ are defined globally in time and satisfy uniform BV estimates, independent of $\ve$. Moreover, they depend continuously on the initial data in the $\L^1$ distance, with a Lipschitz constant independent of $t,\ve$. Letting $\ve\to 0$, these viscous solutions converge to a unique limit, depending Lipschitz continuously on the initial data. In the conservative case where $A=Df$ is the Jacobian of some flux function $f:\R^n\mapsto\R^n$, the vanishing viscosity limits are precisely the unique entropy weak solutions to the system of conservation laws $u_t+f(u)_x=0$.
• ### A note on sigular limits to hyperbolic systems(math/0009012)

Sept. 1, 2000 math.AP
In this note we consider two different singular limits to hyperbolic system of conservation laws, namely the standard backward schemes for non linear semigroups and the semidiscrete scheme. Under the assumption that the rarefaction curve of the corresponding hyperbolic system are straight lines, we prove the stability of the solution and the convergence to the perturbed system to the unique solution of the limit system for initial data with small total variation.
• ### On the Stability of the Standard Riemann Semigroup(math/0007055)

July 10, 2000 math.AP
We consider the dependence of the entropic solution of a hyperbolic system of conservation laws $\{\{array}{c} u_t + f(u)_x = 0 u(0,\cdot) = u_0 \{array}$ on the flux function f. We prove that the solution in Lipschitz continuous w.r.t.~the $C^0$ norm of the derivative of the perturbation of f. We apply this result to prove the convergence of the solution of the relativistic Euler equation to the classical limit.
• ### Stability of $L^\infty$ solutions for hyperbolic systems with coinciding shocks and rarefactions(math/0006094)

June 13, 2000 math.AP
We consider a hyperbolic system of conservation laws u_t + f(u)_x = 0 and u(0,\cdot) = u_0, where each characteristic field is either linearly degenerate or genuinely nonlinear. Under the assumption of coinciding shock and rarefaction curves and the existence of a set of Riemann coordinates $w$, we prove that there exists a semigroup of solutions $u(t) = \mathcal{S}_t u_0$, defined on initial data $u_0 \in L^\infty$. The semigroup $\mathcal{S}$ is continuous w.r.t. time and the initial data $u_0$ in the $L^1_{\text{loc}}$ topology. Moreover $\mathcal{S}$ is unique and its trajectories are obtained as limits of wave front tracking approximations.