• ### Vanishing Viscosity Solutions for Conservation Laws with Regulated Flux(1805.01766)

May 4, 2018 math.AP
In this paper we introduce a concept of "regulated function" $v(t,x)$ of two variables, which reduces to the classical definition when $v$ is independent of $t$. We then consider a scalar conservation law of the form $u_t+F(v(t,x),u)_x=0$, where $F$ is smooth and $v$ is a regulated function, possibly discontinuous w.r.t.both $t$ and $x$. By adding a small viscosity, one obtains a well posed parabolic equation. As the viscous term goes to zero, the uniqueness of the vanishing viscosity limit is proved, relying on comparison estimates for solutions to the corresponding Hamilton--Jacobi equation. As an application, we obtain the existence and uniqueness of solutions for a class of $2\times2$ triangular systems of conservation laws with hyperbolic degeneracy.
• ### On the Optimal Shape of Tree Roots and Branches(1803.01042)

March 2, 2018 math.OC
This paper introduces two classes of variational problems, determining optimal shapes for tree roots and branches. Given a measure $\mu$, describing the distribution of leaves, we introduce a sunlight functional $\S(\mu)$ computing the total amount of light captured by the leaves. On the other hand, given a measure $\mu$ describing the distribution of root hair cells, we consider a harvest functional $\H(\mu)$ computing the total amount of water and nutrients gathered by the roots. In both cases, we seek to maximize these functionals subject to a ramified transportation cost, for transporting nutrients from the roots to the trunk and from the trunk to the leaves. The main results establish various properties of these functionals, and the existence of optimal distributions. In particular, we prove the upper semicontinuity of $\S$ and $\H$, together with a priori estimates on the support of optimal distributions.
• ### On finite time BV blow-up for the p-system(1710.03382)

Oct. 10, 2017 math.AP
The paper studies the possible blowup of the total variation for entropy weak solutions of the p-system, modeling isentropic gas dynamics. It is assumed that the density remains uniformly positive, while the initial data can have arbitrarily large total variation (measured in terms of Riemann invariants). Two main results are proved. (I) If the total variation blows up in finite time, then the solution must contain an infinite number of large shocks in a neighborhood of some point in the $t$-$x$ plane. (II) Piecewise smooth approximate solutions can be constructed whose total variation blows up in finite time. For these solutions the strength of waves emerging from each interaction is exact, while rarefaction waves satisfy the natural decay estimates stemming from the assumption of genuine nonlinearity.
• ### A Model of Controlled Growth(1608.08645)

July 29, 2017 math.AP
We consider a free boundary problem for a system of PDEs, modeling the growth of a biological tissue. A morphogen, controlling volume growth, is produced by specific cells and then diffused and absorbed throughout the domain. The geometric shape of the growing tissue is determined by the instantaneous minimization of an elastic deformation energy, subject to a constraint on the volumetric growth. For an initial domain with $C^{2,\alpha}$ boundary, our main result establishes the local existence and uniqueness of a classical solution, up to a rigid motion.
• ### On the Controllability of Lagrangian Systems by Active Constraints(1702.04824)

Feb. 16, 2017 math-ph, math.MP, math.CA, math.OC
We consider a mechanical system which is controlled by means of moving constraints. Namely, we assume that some of the coordinates can be directly assigned as functions of time by means of frictionless constraints. This leads to a system of ODE's whose right hand side depends quadratically on the time derivative of the control. In this paper we introduce a simplified dynamics, described by a differential inclusion. We prove that every trajectory of the differential inclusion can be uniformly approximated by a trajectory of the original system, on a sufficiently large time interval, starting at rest. Under a somewhat stronger assumption, we show this second trajectory reaches exactly the same terminal point.
• ### Growth Models for Tree Stems and Vines(1701.06680)

Jan. 23, 2017 math.OC
The paper introduces a PDE model for the growth of a tree stem or a vine. The equations describe the elongation due to cell growth, and the response to gravity and to external obstacles. An additional term accounts for the tendency of a vine to curl around branches of other plants. When obstacles are present, the model takes the form of a differential inclusion with unilateral constraints. At each time t, a cone of admissible reactions is determined by the minimization of an elastic deformation energy. The main theorem shows that local solutions exist and can be prolonged globally in time, except when a specific "breakdown configuration" is reached. Approximate solutions are constructed by an operator-splitting technique. Some numerical simulations are provided at the end of the paper.
• ### A Stochastic Model of Optimal Debt Management and Bankruptcy(1609.05983)

Sept. 23, 2016 math.OC
A problem of optimal debt management is modeled as a noncooperative game between a borrower and a pool of lenders, in infinite time horizon with exponential discount. The yearly income of the borrower is governed by a stochastic process. When the debt-to-income ratio $x(t)$ reaches a given size $x^*$, bankruptcy instantly occurs. The interest rate charged by the risk-neutral lenders is precisely determined in order to compensate for this possible loss of their investment. For a given bankruptcy threshold $x^*$, existence and properties of optimal feedback strategies for the borrower are studied, in a stochastic framework as well as in a limit deterministic setting. The paper also analyzes how the expected total cost to the borrower changes, depending on different values of $x^*$, changes, depending on different values of $x^*$?.
• ### Existence of optima and equilibria for traffic flow on networks(1211.1355)

March 24, 2016 math.AP
This paper is concerned with a conservation law model of traffic flow on a network of roads, where each driver chooses his own departure time in order to minimize the sum of a departure cost and an arrival cost. The model includes various groups of drivers, with different origins and destinations and having different cost functions. Under a natural set of assumptions, two main results are proved: (i) the existence of a globally optimal solution, minimizing the sum of the costs to all drivers, and (ii) the existence of a Nash equilibrium solution, where no driver can lower his own cost by changing his departure time or the route taken to reach destination. In the case of Nash solutions, all departure rates are uniformly bounded and have compact support.
• ### The Riemann Solver for Traffic Flow at an Intersection with Buffer of Vanishing Size(1512.03269)

Dec. 10, 2015 math.AP
The paper examines the model of traffic flow at an intersection introduced in [2], containing a buffer with limited size. As the size of the buffer approach zero, it is proved that the solution of the Riemann problem with buffer converges to a self-similar solution described by a specific Limit Riemann Solver (LRS). Remarkably, this new Riemann Solver depends Lipschitz continuously on all parameters.
• ### Lipschitz Metrics for a Class of Nonlinear Wave Equations(1506.06310)

June 21, 2015 math.AP
The nonlinear wave equation $u_{tt}-c(u)(c(u)u_x)_x=0$ determines a flow of conservative solutions taking values in the space $H^1(\mathbb{R})$. However, this flow is not continuous w.r.t. the natural $H^1$ distance. Aim of this paper is to construct a new metric which renders the flow uniformly Lipschitz continuous on bounded subsets of $H^1(\mathbb{R})$. For this purpose, $H^1$ is given the structure of a Finsler manifold, where the norm of tangent vectors is defined in terms of an optimal transportation problem. For paths of piecewise smooth solutions, one can carefully estimate how the weighted length grows in time. By the generic regularity result proved in [7], these piecewise regular paths are dense and can be used to construct a geodesic distance with the desired Lipschitz property.
• ### Structurally Stable Singularities for a Nonlinear Wave Equation(1503.08807)

March 30, 2015 math.AP
For the nonlinear wave equation $u_{tt} - c(u)\big(c(u) u_x\big)_x~=~0$, it is well known that solutions can develop singularities in finite time. For an open dense set of initial data, the present paper provides a detailed asymptotic description of the solution in a neighborhood of each singular point, where $|u_x|\to\infty$. The different structure of conservative and dissipative solutions is analyzed.
• ### Generic Regularity of Conservative Solutions to a Nonlinear Wave Equation(1502.02611)

Feb. 9, 2015 math.AP
The paper is concerned with conservative solutions to the nonlinear wave equation $u_{tt} - c(u)\big(c(u) u_x\big)_x = 0$. For an open dense set of $C^3$ initial data, we prove that the solution is piecewise smooth in the $t$-$x$ plane, while the gradient $u_x$ can blow up along finitely many characteristic curves. The analysis is based on a variable transformation introduced in [7], which reduces the equation to a semilinear system with smooth coefficients, followed by an application of Thom's transversality theorem.
• ### No BV bounds for approximate solutions to p-system with general pressure law(1502.02616)

Feb. 9, 2015 math.AP
For the p-system with large BV initial data, an assumption introduced in [3] by Bakhvalov guarantees the global existence of entropy weak solutions with uniformly bounded total variation. The present paper provides a partial converse to this result. Whenever Bakhvalov's condition does not hold, we show that there exist front tracking approximate solutions, with uniformly positive density, whose total variation becomes arbitrarily large. The construction extends the arguments in [4] to a general class of pressure laws.
• ### Unique Conservative Solutions to a Variational Wave Equation(1411.2012)

Nov. 7, 2014 math.AP
Relying on the analysis of characteristics, we prove the uniqueness of conservative solutions to the variational wave equation $u_{tt}-c(u) (c(u)u_x)_x=0$. Given a solution $u(t,x)$, even if the wave speed $c(u)$ is only H\"older continuous in the $t$-$x$ plane, one can still define forward and backward characteristics in a unique way. Using a new set of independent variables $X,Y$, constant along characteristics, we prove that $t,x,u$, together with other variables, satisfy a semilinear system with smooth coefficients. From the uniqueness of the solution to this semilinear system, one obtains the uniqueness of conservative solutions to the Cauchy problem for the wave equation with general initial data $u(0,\cdot)\in H^1(\mathbb{R})$, $u_t(0,\cdot)\in L^2(\mathbb{R})$.
• ### Representation of Dissipative Solutions to a Nonlinear Variational Wave Equation(1407.1220)

July 4, 2014 math.AP
The paper introduces a new way to construct dissipative solutions to a second order variational wave equation. By a variable transformation, from the nonlinear PDE one obtains a semilinear hyperbolic system with sources. In contrast with the conservative case, here the source terms are discontinuous and the discontinuities are not always crossed transversally. Solutions to the semilinear system are obtained by an approximation argument, relying on Kolmogorov's compactness theorem. Reverting to the original variables, one recovers a solution to the nonlinear wave equation where the total energy is a monotone decreasing function of time.
• ### Uniqueness of Conservative Solutions to the Camassa-Holm Equation via Characteristics(1401.0312)

Jan. 1, 2014 math.AP
The paper provides a direct proof the uniqueness of solutions to the Camassa-Holm equation, based on characteristics. Given a conservative solution $u=u(t,x)$, an equation is introduced which singles out a unique characteristic curve through each initial point. By studying the evolution of the quantities $u$ and $v= 2\arctan u_x$ along each characteristic, it is proved that the Cauchy problem with general initial data $u_0\in H^1(\mathbb{R})$ has a unique solution, globally in time.
• ### Non-existence and Non-uniqueness for Multidimensional Sticky Particle Systems(1312.1636)

Dec. 5, 2013 math.AP, math-ph, math.MP
The paper is concerned with sticky weak solutions to the equations of pressureless gases in two or more space dimensions. Various initial data are constructed, showing that the Cauchy problem can have (i) two distinct sticky solutions, or (ii) no sticky solution, not even locally in time. In both cases the initial density is smooth with compact support, while the initial velocity field is continuous.
• ### Lack of BV Bounds for Approximate Solutions to the p-system with Large Data(1310.7058)

Oct. 26, 2013 math.AP
We consider front tracking approximate solutions to the p-system of isentropic gas dynamics. At interaction times, the outgoing wave fronts have the same strength as in the exact solution of the Riemann problem, but some error is allowed in their speed. For large BV initial data, we construct examples showing that the total variation of these approximate solutions can become arbitrarily large, or even blow up in finite time. This happens even if the density of the gas remains uniformly positive.
• ### On the Control of Non Holonomic Systems by Active Constraints(1208.4335)

Aug. 21, 2012 math.DS
The paper is concerned with mechanical systems which are controlled by implementing a number of time-dependent, frictionless holonomic constraints. The main novelty is due to the presence of additional non-holonomic constraints. We develop a general framework to analyze these problems, deriving the equations of motion and studying the continuity properties of the "control-to-trajectory" maps. Various geometric characterizations are provided, in order that the equations be affine w.r.t. the time derivative of the control. In this case the system is fit for jumps, and the evolution is well defined also in connection with discontinuous control functions. The classical Roller Racer provides an example where the non-affine dependence of the equations on the derivative of the control is due only to the non-holonomic constraint. This is a case where the presence of quadratic terms in the equations can be used for controllability purposes.
• ### Lipschitz metric for the Hunter-Saxton equation(0904.3615)

April 23, 2009 math.AP
We study stability of solutions of the Cauchy problem for the Hunter-Saxton equation $u_t+uu_x=\frac14(\int_{-\infty}^xu_x^2 dx-\int_{x}^\infty u_x^2 dx)$ with initial data $u_0$. In particular, we derive a new Lipschitz metric $d_\D$ with the property that for two solutions $u$ and $v$ of the equation we have $d_\D(u(t),v(t))\le e^{Ct} d_\D(u_0,v_0)$.
• ### Nearly Optimal Patchy Feedbacks for Minimization Problems with Free Terminal Time(0708.3897)

Oct. 31, 2007 math.CA, math.OC
The paper is concerned with a general optimization problem for a nonlinear control system, in the presence of a running cost and a terminal cost, with free terminal time. We prove the existence of a patchy feedback whose trajectories are all nearly optimal solutions, with pre-assigned accuracy.
• ### Nearly Time Optimal Stabilizing Patchy Feedbacks(math/0512531)

Dec. 22, 2005 math.CA
We consider the time optimal stabilization problem for a nonlinear control system $\dot x=f(x,u)$. Let $\tau(y)$ be the minimum time needed to steer the system from the state $y\in\R^n$ to the origin, and call $\A(T)$ the set of initial states that can be steered to the origin in time $\tau(y)\leq T$. Given any $\ve>0$, in this paper we construct a patchy feedback $u=U(x)$ such that every solution of $\dot x=f(x, U(x))$, $x(0)=y\in \A(T)$ reaches an $\ve$-neighborhood of the origin within time $\tau(y)+\ve$.
• ### Infinite Horizon Noncooperative Differential Games(math/0505212)

May 11, 2005 math.AP
For a non-cooperative differential game, the value functions of the various players satisfy a system of Hamilton-Jacobi equations. In the present paper, we consider a class of infinite-horizon games with nonlinear costs exponentially discounted in time. By the analysis of the value functions, we establish the existence of Nash equilibrium solutions in feedback form and provide results and counterexamples on their uniqueness and stability.
• ### An Optimal Transportation Metric for Solutions of the Camassa-Holm Equation(math/0504450)

May 4, 2005 math.AP
In this paper we construct a global, continuous flow of solutions to the Camassa-Holm equation on the entire space $H^1$. Our solutions are conservative, in the sense that the total energy $\int (u^2+u_x^2) dx$ remains a.e. constant in time. Our new approach is based on a distance functional $J(u,v)$, defined in terms of an optimal transportation problem, which satisfies ${d\over dt} J(u(t), v(t))\leq \kappa\cdot J(u(t),v(t))$ for every couple of solutions. Using this new distance functional, we can construct arbitrary solutions as the uniform limit of multi-peakon solutions, and prove a general uniqueness result.
• ### On Asymptotic Variational Wave Equations(math/0502124)

Feb. 7, 2005 math.AP
We investigate the equation $(u_t + (f(u))_x)_x = f''(u) (u_x)^2/2$ where $f(u)$ is a given smooth function. Typically $f(u)= u^2/2$ or $u^3/3$. This equation models unidirectional and weakly nonlinear waves for the variational wave equation $u_{tt} - c(u) (c(u)u_x)_x =0$ which models some liquid crystals with a natural sinusoidal $c$. The equation itself is also the Euler-Lagrange equation of a variational problem. Two natural classes of solutions can be associated with this equation. A conservative solution will preserve its energy in time, while a dissipative weak solution loses energy at the time when singularities appear. Conservative solutions are globally defined, forward and backward in time, and preserve interesting geometric features, such as the Hamiltonian structure. On the other hand, dissipative solutions appear to be more natural from the physical point of view. We establish the well-posedness of the Cauchy problem within the class of conservative solutions, for initial data having finite energy and assuming that the flux function $f$ has Lipschitz continuous second-order derivative. In the case where $f$ is convex, the Cauchy problem is well-posed also within the class of dissipative solutions. However, when $f$ is not convex, we show that the dissipative solutions do not depend continuously on the initial data.