• 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.
  • 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.
  • 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}.
  • 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$.
  • 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)$.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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$.
  • We prove the SBV regularity of the characteristic speed of the scalar hyperbolic conservation law and SBV-like regularity of the eigenvalue functions of the Jacobian matrix of flux function for general systems of conservation laws. More precisely, for the equation u_t + f(u)_x = 0, \quad u : \R^+ \times \R \to \Omega \subset \R^N, we only assume the flux $f$ is $C^2$ function in the scalar case (N=1) and Jacobian matrix $Df$ has distinct real eigenvalues in the system case $(N\geq 2)$. Using the modification of the main decay estimate in Lau and localization method applied in \cite{R}, we show that for the scalar equation $f'(u)$ belongs to SBV, and for system of conservation laws the scalar measure \[\big(D_u \lambda_i(u) \cdot r_i(u) \big) \big(l_i(u) \cdot u_x \big)] has no Cantor part, where $\lambda_i$, $r_i$, $l_i$ are the $i$-th eigenvalue, $i$-th right eigenvector and $i$-th left eigenvector of the matrix $Df$.
  • We prove that if $t \mapsto u(t) \in \mathrm {BV}(\R)$ is the entropy solution to a $N \times N$ strictly hyperbolic system of conservation laws with genuinely nonlinear characteristic fields \[ u_t + f(u)_x = 0, \] then up to a countable set of times $\{t_n\}_{n \in \mathbb N}$ the function $u(t)$ is in $\mathrm {SBV}$, i.e. its distributional derivative $u_x$ is a measure with no Cantorian part. The proof is based on the decomposition of $u_x(t)$ into waves belonging to the characteristic families \[ u(t) = \sum_{i=1}^N v_i(t) \tilde r_i(t), \quad v_i(t) \in \mathcal M(\R), \ \tilde r_i(t) \in \mathrm R^N, \] and the balance of the continuous/jump part of the measures $v_i$ in regions bounded by characteristics. To this aim, a new interaction measure $\mu_{i,\jump}$ is introduced, controlling the creation of atoms in the measure $v_i(t)$. The main argument of the proof is that for all $t$ where the Cantorian part of $v_i$ is not 0, either the Glimm functional has a downward jump, or there is a cancellation of waves or the measure $\mu_{i,\mathrm{jump}}$ is positive.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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$.
  • 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.
  • 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.
  • 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.