-
Two Nitsche-based mixed finite element discretizations for the seepage problem in Richards' equation
Authors:
Federico Gatti,
Andrea Bressan,
Alessio Fumagalli,
Domenico Gallipoli,
Leonardo Maria Lalicata,
Simone Pittaluga,
Lorenzo Tamellini
Abstract:
This paper proposes two algorithms to impose seepage boundary conditions in the context of Richards' equation for groundwater flows in unsaturated media. Seepage conditions are non-linear boundary conditions, that can be formulated as a set of unilateral constraints on both the pressure head and the water flux at the ground surface, together with a complementarity condition: these conditions in pr…
▽ More
This paper proposes two algorithms to impose seepage boundary conditions in the context of Richards' equation for groundwater flows in unsaturated media. Seepage conditions are non-linear boundary conditions, that can be formulated as a set of unilateral constraints on both the pressure head and the water flux at the ground surface, together with a complementarity condition: these conditions in practice require switching between Neumann and Dirichlet boundary conditions on unknown portions on the boundary. Upon realizing the similarities of these conditions with unilateral contact problems in mechanics, we take inspiration from that literature to propose two approaches: the first method relies on a strongly consistent penalization term, whereas the second one is obtained by an hybridization approach, in which the value of the pressure on the surface is treated as a separate set of unknowns. The flow problem is discretized in mixed form with div-conforming elements so that the water mass is preserved. Numerical experiments show the validity of the proposed strategy in handling the seepage boundary conditions on geometries with increasing complexity.
△ Less
Submitted 10 July, 2024;
originally announced July 2024.
-
Intermediate Domains for Scalar Conservation Laws
Authors:
Fabio Ancona,
Alberto Bressan,
Elio Marconi,
Luca Talamini
Abstract:
For a scalar conservation law with strictly convex flux, by Oleinik's estimates the total variation of a solution with initial data $\overline{u}\in \bf{L}^\infty(\mathbb R)$ decays like $t^{-1}$. This paper introduces a class of intermediate domains $\mathcal P_α$, $0<α<1$, such that for $\overline u\in \mathcal P_α$ a faster decay rate is achieved:…
▽ More
For a scalar conservation law with strictly convex flux, by Oleinik's estimates the total variation of a solution with initial data $\overline{u}\in \bf{L}^\infty(\mathbb R)$ decays like $t^{-1}$. This paper introduces a class of intermediate domains $\mathcal P_α$, $0<α<1$, such that for $\overline u\in \mathcal P_α$ a faster decay rate is achieved: $\mathrm{Tot.Var.}\bigl\{ u(t,\cdot)\bigr\}\sim t^{α-1}$. A key ingredient of the analysis is a ``Fourier-type" decomposition of $\overline u$ into components which oscillate more and more rapidly. The results aim at extending the theory of fractional domains for analytic semigroups to an entirely nonlinear setting.
△ Less
Submitted 16 April, 2024;
originally announced April 2024.
-
Generic Properties of Conjugate Points in Optimal Control Problems
Authors:
Alberto Bressan,
Marco Mazzola,
Khai T. Nguyen
Abstract:
The first part of the paper studies a class of optimal control problems in Bolza form, where the dynamics is linear w.r.t.~the control function. A necessary condition is derived, for the optimality of a trajectory which starts at a conjugate point. The second part is concerned with a classical problem in the Calculus of Variations, with free terminal point. For a generic terminal cost…
▽ More
The first part of the paper studies a class of optimal control problems in Bolza form, where the dynamics is linear w.r.t.~the control function. A necessary condition is derived, for the optimality of a trajectory which starts at a conjugate point. The second part is concerned with a classical problem in the Calculus of Variations, with free terminal point. For a generic terminal cost $ψ\in \C^4(\mathbb{R}^n)$, applying the previous necessary condition we show that the set of conjugate points is contained in the image of an $(n-2)$-dimensional manifold, and has locally bounded $(n-2)$-dimensional Hausdorff measure.
△ Less
Submitted 2 April, 2024;
originally announced April 2024.
-
Parallelization in time by diagonalization
Authors:
Andrea Bressan,
Alen Kushova,
Gabriele Loli,
Monica Montardini,
Giancarlo Sangalli,
Mattia Tani
Abstract:
This is a review of preconditioning techniques based on fast-diagonalization methods for space-time isogeometric discretization of the heat equation. Three formulation are considered: the Galerkin approach, a discrete least-square and a continuous least square. For each formulation the heat differential operator is written as a sum of terms that are kronecker products of uni-variate operators. The…
▽ More
This is a review of preconditioning techniques based on fast-diagonalization methods for space-time isogeometric discretization of the heat equation. Three formulation are considered: the Galerkin approach, a discrete least-square and a continuous least square. For each formulation the heat differential operator is written as a sum of terms that are kronecker products of uni-variate operators. These are used to speed-up the application of the operator in iterative solvers and to construct a suitable preconditioner. Contrary to the fast-diagonalization technique for the Laplace equation where all uni-variate operators acting on the same direction can be simultaneously diagonalized in the case of the heat equation this is not possible. Luckily this can be done up to an additional term that has low rank allowing for the utilization of arrow-head like factorization or inversion by Sherman-Morrison formula. The proposed preconditioners work extremely well on the parametric domain and, when the domain is parametrized or when the equation coefficients are not constant, they can be adapted and retain good performance characteristics.
△ Less
Submitted 30 November, 2023;
originally announced March 2024.
-
Energy conservative isogeometric techniques for the wave equation
Authors:
Andrea Bressan,
Annalisa Buffa,
Alen Kushova,
Rafael Vázquez
Abstract:
We analyze the wave equation in mixed form, with periodic and/or Dirichlet homogeneous boundary conditions, and nonconstant coefficients that depend on the spatial variable. For the discretization, the weak form of the second equation is replaced by a strong form, written in terms of a projection operator. The system of equations is discretized with B-splines forming a De Rham complex along with s…
▽ More
We analyze the wave equation in mixed form, with periodic and/or Dirichlet homogeneous boundary conditions, and nonconstant coefficients that depend on the spatial variable. For the discretization, the weak form of the second equation is replaced by a strong form, written in terms of a projection operator. The system of equations is discretized with B-splines forming a De Rham complex along with suitable commutative projectors for the approximation of the second equation. The discrete scheme is energy conservative when discretized in time with a conservative method such as Crank-Nicolson. We propose a convergence analysis of the method to study the dependence with respect to the mesh size $h$, with focus on the consistency error. Numerical results show optimal convergence of the error in energy norm, and a relative error in energy conservation for long-time simulations of the order of machine precision.
△ Less
Submitted 30 November, 2023;
originally announced November 2023.
-
Space-time least squares approximation for Schrödinger equation and efficient solver
Authors:
Andrea Bressan,
Alen Kushova,
Giancarlo Sangalli,
Mattia Tani
Abstract:
In this work we present a space-time least squares isogeometric discretization of the Schrödinger equation and propose a preconditioner for the arising linear system in the parametric domain. Exploiting the tensor product structure of the basis functions, the preconditioner is written as the sum of Kronecker products of matrices. Thanks to an extension to the classical Fast Diagonalization method,…
▽ More
In this work we present a space-time least squares isogeometric discretization of the Schrödinger equation and propose a preconditioner for the arising linear system in the parametric domain. Exploiting the tensor product structure of the basis functions, the preconditioner is written as the sum of Kronecker products of matrices. Thanks to an extension to the classical Fast Diagonalization method, the application of the preconditioner is efficient and robust w.r.t. the polynomial degree of the spline space. The time required for the application is almost proportional to the number of degrees-of-freedom, for a serial execution.
△ Less
Submitted 30 November, 2023;
originally announced November 2023.
-
One Dimensional Hyperbolic Conservation Laws: Past and Future
Authors:
Alberto Bressan
Abstract:
Aim of these notes is provide a brief review of the current well-posedness theory for hyperbolic systems of conservation laws in one space dimension, also pointing out open problems and possible research directions. They supplement the slides of the short course given by the author in Erice, May 2023, available at: sites.google.com/view/erice23/speakers-and-slides.
Aim of these notes is provide a brief review of the current well-posedness theory for hyperbolic systems of conservation laws in one space dimension, also pointing out open problems and possible research directions. They supplement the slides of the short course given by the author in Erice, May 2023, available at: sites.google.com/view/erice23/speakers-and-slides.
△ Less
Submitted 25 October, 2023;
originally announced October 2023.
-
Generic Singularities for 2D Pressureless Flow
Authors:
Alberto Bressan,
Geng Chen,
Shoujun Huang
Abstract:
In this paper, we consider the Cauchy problem for pressureless gases in two space dimensions with generic smooth initial data (density and velocity). These equations give rise to singular curves, where the mass has positive density w.r.t.~1-dimensional Hausdorff measure. We observe that the system of equations describing these singular curves is not hyperbolic. For analytic data, local solutions a…
▽ More
In this paper, we consider the Cauchy problem for pressureless gases in two space dimensions with generic smooth initial data (density and velocity). These equations give rise to singular curves, where the mass has positive density w.r.t.~1-dimensional Hausdorff measure. We observe that the system of equations describing these singular curves is not hyperbolic. For analytic data, local solutions are constructed using a version of the Cauchy-Kovalevskaya theorem. We then study the interaction of two singular curves, in generic position. Finally, for a generic initial velocity field, we investigate the asymptotic structure of the smooth solution up to the first time when a singularity is formed.
△ Less
Submitted 21 July, 2023;
originally announced July 2023.
-
Optimal Solutions for a Class of Set-Valued Evolution Problems
Authors:
Stefano Bianchini,
Alberto Bressan,
Maria Teresa Chiri
Abstract:
The paper is concerned with a class of optimization problems for moving sets $t\mapstoΩ(t)\subset\mathbb{R}^2$, motivated by the control of invasive biological populations. Assuming that the initial contaminated set $Ω_0$ is convex, we prove that a strategy is optimal if an only if at each given time $t\in [0,T]$ the control is active along the portion of the boundary $\partial Ω(t)$ where the cur…
▽ More
The paper is concerned with a class of optimization problems for moving sets $t\mapstoΩ(t)\subset\mathbb{R}^2$, motivated by the control of invasive biological populations. Assuming that the initial contaminated set $Ω_0$ is convex, we prove that a strategy is optimal if an only if at each given time $t\in [0,T]$ the control is active along the portion of the boundary $\partial Ω(t)$ where the curvature is maximal. In particular, this implies that $Ω(t)$ is convex for all $t\geq 0$. The proof relies on the analysis of a one-step constrained optimization problem, obtained by a time discretization.
△ Less
Submitted 7 July, 2023;
originally announced July 2023.
-
A remark on the uniqueness of solutions to hyperbolic conservation laws
Authors:
Alberto Bressan,
Camillo De Lellis
Abstract:
Given a strictly hyperbolic $n\times n$ system of conservation laws, it is well known that there exists a unique Lipschitz semigroup of weak solutions, defined on a domain of functions with small total variation, which are limits of vanishing viscosity approximations. Aim of this note is to prove that every weak solution taking values in the domain of the semigroup, and whose shocks satisfy the Li…
▽ More
Given a strictly hyperbolic $n\times n$ system of conservation laws, it is well known that there exists a unique Lipschitz semigroup of weak solutions, defined on a domain of functions with small total variation, which are limits of vanishing viscosity approximations. Aim of this note is to prove that every weak solution taking values in the domain of the semigroup, and whose shocks satisfy the Liu admissibility conditions, actually coincides with a semigroup trajectory.
△ Less
Submitted 26 May, 2023;
originally announced May 2023.
-
Unique Solutions to Hyperbolic Conservation Laws with a Strictly Convex Entropy
Authors:
Alberto Bressan,
Graziano Guerra
Abstract:
Consider a strictly hyperbolic $n\times n$ system of conservation laws, where each characteristic field is either genuinely nonlinear or linearly degenerate. In this standard setting, it is well known that there exists a Lipschitz semigroup of weak solutions, defined on a domain of functions with small total variation. If the system admits a strictly convex entropy, we give a short proof that ever…
▽ More
Consider a strictly hyperbolic $n\times n$ system of conservation laws, where each characteristic field is either genuinely nonlinear or linearly degenerate. In this standard setting, it is well known that there exists a Lipschitz semigroup of weak solutions, defined on a domain of functions with small total variation. If the system admits a strictly convex entropy, we give a short proof that every entropy weak solution taking values within the domain of the semigroup coincides with a semigroup trajectory. The result shows that the assumptions of ``Tame Variation" or ``Tame Oscillation", previously used to achieve uniqueness, can be removed in the presence of a strictly convex entropy.
△ Less
Submitted 18 May, 2023;
originally announced May 2023.
-
Controlled Traveling Profiles for Models of Invasive Biological Species
Authors:
Alberto Bressan,
Minyan Zhang
Abstract:
We consider a family of controlled reaction-diffusion equations, describing the spatial spreading of an invasive biological species. For a given propagation speed $c\in{I\!\!R}$, we seek a control with minimum cost, which achieves a traveling profile with speed $c$. For various nonlinear models, the existence of a (possibly measure valued) optimal control is proved, together with necessary conditi…
▽ More
We consider a family of controlled reaction-diffusion equations, describing the spatial spreading of an invasive biological species. For a given propagation speed $c\in{I\!\!R}$, we seek a control with minimum cost, which achieves a traveling profile with speed $c$. For various nonlinear models, the existence of a (possibly measure valued) optimal control is proved, together with necessary conditions for optimality. In the last section we study a case where the wave speed cannot be modified by any control with finite cost.
The present analysis is motivated by the recent results in arXiv:2201.01723 and arXiv:2108.09321, showing how a control problem for a reaction-diffusion equation can be approximated by a simpler problem of optimal control of a moving set.
△ Less
Submitted 19 February, 2023;
originally announced February 2023.
-
Generic Properties of First Order Mean Field Games
Authors:
Alberto Bressan,
Khai T. Nguyen
Abstract:
We consider a class of deterministic mean field games, where the state associated with each player evolves according to an ODE which is linear w.r.t. the control. Existence, uniqueness, and stability of solutions are studied from the point of view of generic theory. Within a suitable topological space of dynamics and cost functionals, we prove that, for nearly all mean field games(in the Baire cat…
▽ More
We consider a class of deterministic mean field games, where the state associated with each player evolves according to an ODE which is linear w.r.t. the control. Existence, uniqueness, and stability of solutions are studied from the point of view of generic theory. Within a suitable topological space of dynamics and cost functionals, we prove that, for nearly all mean field games(in the Baire category sense) the best reply map is single valued for a.e. player. As a consequence, the mean field game admits a strong (not randomized) solution. Examples are given of open sets of games admitting a single solution, and other open sets admitting multiple solutions. Further examples show the existence of an open set of MFG having a unique solution which is asymptotically stable w.r.t. the best reply map, and another open set of MFG having a unique solution which is unstable. We conclude with an example of a MFG with terminal constraints which does not have any solution, not even in the mild sense with randomized strategies.
△ Less
Submitted 26 October, 2022;
originally announced October 2022.
-
Shock interactions for the Burgers-Hilbert Equation
Authors:
Alberto Bressan,
Sondre T. Galtung,
Katrin Grunert,
Khai T. Nguyen
Abstract:
This paper provides an asymptotic description of a solution to the Burgers-Hilbert equation in a neighborhood of a point where two shocks interact. The solution is obtained as the sum of a function with $H^2$ regularity away from the shocks plus a corrector term having an asymptotic behavior like |x|ln|x| close to each shock. A key step in the analysis is the construction of piecewise smooth solut…
▽ More
This paper provides an asymptotic description of a solution to the Burgers-Hilbert equation in a neighborhood of a point where two shocks interact. The solution is obtained as the sum of a function with $H^2$ regularity away from the shocks plus a corrector term having an asymptotic behavior like |x|ln|x| close to each shock. A key step in the analysis is the construction of piecewise smooth solutions with a single shock for a general class of initial data.
△ Less
Submitted 5 April, 2022;
originally announced April 2022.
-
Optimal Control of Moving Sets
Authors:
Alberto Bressan,
Maria Teresa Chiri,
Najmeh Salehi
Abstract:
Motivated by the control of invasive biological populations, we consider a class of optimization problems for moving sets $t\mapsto Ω(t)\subset\mathbb{R}^2$. Given an initial set $Ω_0$, the goal is to minimize the area of the contaminated set $Ω(t)$ over time, plus a cost related to the control effort. Here the control function is the inward normal speed along the boundary $\partial Ω(t)$. We prov…
▽ More
Motivated by the control of invasive biological populations, we consider a class of optimization problems for moving sets $t\mapsto Ω(t)\subset\mathbb{R}^2$. Given an initial set $Ω_0$, the goal is to minimize the area of the contaminated set $Ω(t)$ over time, plus a cost related to the control effort. Here the control function is the inward normal speed along the boundary $\partial Ω(t)$. We prove the existence of optimal solutions, within a class of sets with finite perimeter. Necessary conditions for optimality are then derived, in the form of a Pontryagin maximum principle. Additional optimality conditions show that the sets $Ω(t)$ cannot have certain types of outward or inward corners. Finally, some explicit solutions are presented.
△ Less
Submitted 5 January, 2022;
originally announced January 2022.
-
On the Optimal Control of Propagation Fronts
Authors:
Alberto Bressan,
Maria Teresa Chiri,
Najmeh Salehi
Abstract:
We consider a controlled reaction-diffusion equation, motivated by a pest eradication problem. Our goal is to derive a simpler model, describing the controlled evolution of a contaminated set. In this direction, the first part of the paper studies the optimal control of 1-dimensional traveling wave profiles. Using Stokes' formula, explicit solutions are obtained, which in some cases require measur…
▽ More
We consider a controlled reaction-diffusion equation, motivated by a pest eradication problem. Our goal is to derive a simpler model, describing the controlled evolution of a contaminated set. In this direction, the first part of the paper studies the optimal control of 1-dimensional traveling wave profiles. Using Stokes' formula, explicit solutions are obtained, which in some cases require measure-valued optimal controls. In the last section we introduce a family of optimization problems for a moving set. We show how these can be derived from the original parabolic problems, by taking a sharp interface limit.
△ Less
Submitted 20 August, 2021;
originally announced August 2021.
-
Optimal Shapes for Tree Roots
Authors:
Alberto Bressan,
Sondre T. Galtung,
Qing Sun
Abstract:
The paper studies a class of variational problems, modeling optimal shapes for tree roots. Given a measure $μ$ describing the distribution of root hair cells, we seek to maximize a harvest functional $\mathcal{H}$, computing the total amount of water and nutrients gathered by the roots, subject to a cost for transporting these nutrients from the roots to the trunk. Earlier papers had established t…
▽ More
The paper studies a class of variational problems, modeling optimal shapes for tree roots. Given a measure $μ$ describing the distribution of root hair cells, we seek to maximize a harvest functional $\mathcal{H}$, computing the total amount of water and nutrients gathered by the roots, subject to a cost for transporting these nutrients from the roots to the trunk. Earlier papers had established the existence of an optimal measure, and a priori bounds. Here we derive necessary conditions for optimality. Moreover, in space dimension $d=2$, we prove that the support of an optimal measure is nowhere dense.
△ Less
Submitted 3 May, 2022; v1 submitted 11 August, 2021;
originally announced August 2021.
-
Diffusion Approximations of Markovian Solutions to Discontinuous ODEs
Authors:
Alberto Bressan,
Marco Mazzola,
Khai T. Nguyen
Abstract:
In a companion paper, the authors have characterized all deterministic semigroups, and all Markov semigroups, whose trajectories are Carathe'odory solutions to a given ODE x'=f(x), with f possibly discontinuous. The present paper establishes two approximation results. Namely, every deterministic semigroup can be obtained as the pointwise limit of the flows generated by a sequence of ODEs $x'=f_n(x…
▽ More
In a companion paper, the authors have characterized all deterministic semigroups, and all Markov semigroups, whose trajectories are Carathe'odory solutions to a given ODE x'=f(x), with f possibly discontinuous. The present paper establishes two approximation results. Namely, every deterministic semigroup can be obtained as the pointwise limit of the flows generated by a sequence of ODEs $x'=f_n(x) with smooth right hand sides. Moreover, every Markov semigroup can be obtained as limit of a sequence of diffusion processes with smooth drifts and with diffusion coefficients approaching zero.
△ Less
Submitted 5 July, 2021;
originally announced July 2021.
-
On the Regularity of Optimal Dynamic Blocking Strategies
Authors:
Alberto Bressan,
Maria Teresa Chiri
Abstract:
The paper studies a dynamic blocking problem, motivated by a model of optimal fire confinement. While the fire can expand with unit speed in all directions, barriers are constructed in real time. An optimal strategy is sought, minimizing the total value of the burned region, plus a construction cost. It is well known that optimal barriers exists. In general, they are a countable union of compact,…
▽ More
The paper studies a dynamic blocking problem, motivated by a model of optimal fire confinement. While the fire can expand with unit speed in all directions, barriers are constructed in real time. An optimal strategy is sought, minimizing the total value of the burned region, plus a construction cost. It is well known that optimal barriers exists. In general, they are a countable union of compact, connected, rectifiable sets. The main result of the present paper shows that optimal barriers are nowhere dense. The proof relies on new estimates on the reachable sets and on optimal trajectories for the fire, solving a minimum time problem in the presence of obstacles.
△ Less
Submitted 1 December, 2020;
originally announced December 2020.
-
Entropy Admissibility of the Limit Solution for a Nonlocal Model of Traffic Flow
Authors:
Alberto Bressan,
Wen Shen
Abstract:
We consider a conservation law model of traffic flow, where the velocity of each car depends on a weighted average of the traffic density $ρ$ ahead. The averaging kernel is of exponential type: $w_\varepsilon(s)=\varepsilon^{-1} e^{-s/\varepsilon}$. For any decreasing velocity function $v$, we prove that, as $\varepsilon\to 0$, the limit of solutions to the nonlocal equation coincides with the uni…
▽ More
We consider a conservation law model of traffic flow, where the velocity of each car depends on a weighted average of the traffic density $ρ$ ahead. The averaging kernel is of exponential type: $w_\varepsilon(s)=\varepsilon^{-1} e^{-s/\varepsilon}$. For any decreasing velocity function $v$, we prove that, as $\varepsilon\to 0$, the limit of solutions to the nonlocal equation coincides with the unique entropy-admissible solution to the scalar conservation law $ρ_t + (ρv(ρ))_x=0$.
△ Less
Submitted 10 November, 2020;
originally announced November 2020.
-
A posteriori Error Estimates for Numerical Solutions to Hyperbolic Conservation Laws
Authors:
Alberto Bressan,
Maria Teresa Chiri,
Wen Shen
Abstract:
The paper is concerned with a posteriori error bounds for a wide class of numerical schemes, for $n\times n$ hyperbolic conservation laws in one space dimension. These estimates are achieved by a "post-processing algorithm", checking that the numerical solution retains small total variation, and computing its oscillation on suitable subdomains. The results apply, in particular, to solutions obtain…
▽ More
The paper is concerned with a posteriori error bounds for a wide class of numerical schemes, for $n\times n$ hyperbolic conservation laws in one space dimension. These estimates are achieved by a "post-processing algorithm", checking that the numerical solution retains small total variation, and computing its oscillation on suitable subdomains. The results apply, in particular, to solutions obtained by the Godunov or the Lax-Friedrichs scheme, backward Euler approximations, and the method of periodic smoothing. Some numerical implementations are presented.
△ Less
Submitted 1 October, 2020;
originally announced October 2020.
-
Numerical Study of Non-uniqueness for 2D Compressible Isentropic Euler Equations
Authors:
Alberto Bressan,
Yi Jiang,
Hailiang Liu
Abstract:
In this paper, we numerically study a class of solutions with spiraling singularities in vorticity for two-dimensional, inviscid, compressible Euler systems, where the initial data have an algebraic singularity in vorticity at the origin. These are different from the multi-dimensional Riemann problems widely studied in the literature. Our computations provide numerical evidence of the existence of…
▽ More
In this paper, we numerically study a class of solutions with spiraling singularities in vorticity for two-dimensional, inviscid, compressible Euler systems, where the initial data have an algebraic singularity in vorticity at the origin. These are different from the multi-dimensional Riemann problems widely studied in the literature. Our computations provide numerical evidence of the existence of initial value problems with multiple solutions, thus revealing a fundamental obstruction toward the well-posedness of the governing equations. The compressible Euler equations are solved using the positivity-preserving discontinuous Galerkin method.
△ Less
Submitted 27 August, 2021; v1 submitted 20 September, 2020;
originally announced September 2020.
-
Markovian Solutions to Discontinuous ODEs
Authors:
Alberto Bressan,
Marco Mazzola,
Khai T. Nguyen
Abstract:
Given a possibly discontinuous, bounded function $f:\mathbb{R}\mapsto\mathbb{R}$, we consider the set of generalized flows, obtained by assigning a probability measure on the set of Carathéodory solutions to the ODE ~$\dot x = f(x)$. The paper provides a complete characterization of all such flows which have a Markov property in time. This is achieved in terms of (i) a positive, atomless measure s…
▽ More
Given a possibly discontinuous, bounded function $f:\mathbb{R}\mapsto\mathbb{R}$, we consider the set of generalized flows, obtained by assigning a probability measure on the set of Carathéodory solutions to the ODE ~$\dot x = f(x)$. The paper provides a complete characterization of all such flows which have a Markov property in time. This is achieved in terms of (i) a positive, atomless measure supported on the set $f^{-1}(0)$ where $f$ vanishes, (ii) a countable number of Poisson random variables, determining the waiting times at points in $f^{-1}(0)$, and (iii) a countable set of numbers $θ_k\in [0,1]$, describing the probability of moving up or down, at isolated points where two distinct trajectories can originate.
△ Less
Submitted 11 September, 2020;
originally announced September 2020.
-
On a Shape Optimization Problem for Tree Branches
Authors:
Alberto Bressan,
Sondre T. Galtung
Abstract:
This paper is concerned with a shape optimization problem, where the functional to be maximized describes the total sunlight collected by a distribution of tree leaves, minus the cost for transporting water and nutrient from the base of the trunk to all the leaves. In the case of 2 space dimensions, the solution is proved to be unique, and explicitly determined.
This paper is concerned with a shape optimization problem, where the functional to be maximized describes the total sunlight collected by a distribution of tree leaves, minus the cost for transporting water and nutrient from the base of the trunk to all the leaves. In the case of 2 space dimensions, the solution is proved to be unique, and explicitly determined.
△ Less
Submitted 2 June, 2020;
originally announced June 2020.
-
Soil searching by an artificial root
Authors:
Fabio Ancona,
Alberto Bressan,
Maria Teresa Chiri
Abstract:
We model an artificial root which grows in the soil for underground prospecting. Its evolution is described by a controlled system of two integro-partial differential equations: one for the growth of the body and the other for the elongation of the tip. At any given time, the angular velocity of the root is obtained by solving a minimization problem with state constraints. We prove the existence o…
▽ More
We model an artificial root which grows in the soil for underground prospecting. Its evolution is described by a controlled system of two integro-partial differential equations: one for the growth of the body and the other for the elongation of the tip. At any given time, the angular velocity of the root is obtained by solving a minimization problem with state constraints. We prove the existence of solutions to the evolution problem, up to the first time where a "breakdown configuration" is reached. Some numerical simulations are performed, to test the effectiveness of our feedback control algorithm.
△ Less
Submitted 12 May, 2020;
originally announced May 2020.
-
A Posteriori Error Estimates for Self-Similar Solutions to the Euler Equations
Authors:
Alberto Bressan,
Wen Shen
Abstract:
The main goal of this paper is to analyze a family of "simplest possible" initial data for which, as shown by numerical simulations, the incompressible Euler equations have multiple solutions. We take here a first step toward a rigorous validation of these numerical results. Namely, we consider the system of equations corresponding to a self-similar solution, restricted to a bounded domain with sm…
▽ More
The main goal of this paper is to analyze a family of "simplest possible" initial data for which, as shown by numerical simulations, the incompressible Euler equations have multiple solutions. We take here a first step toward a rigorous validation of these numerical results. Namely, we consider the system of equations corresponding to a self-similar solution, restricted to a bounded domain with smooth boundary. Given an approximate solution obtained via a finite dimensional Galerkin method, we establish a posteriori error bounds on the distance between the numerical approximation and the exact solution having the same boundary data.
△ Less
Submitted 5 February, 2020;
originally announced February 2020.
-
Variational Problems for Tree Roots and Branches
Authors:
Alberto Bressan,
Michele Palladino,
Qing Sun
Abstract:
This paper studies two classes of variational problems introduced in [7], related to the optimal shapes of tree roots and branches. Given a measure $μ$ describing the distribution of leaves, a sunlight functional $§(μ)$ computes the total amount of light captured by the leaves. For a measure $μ$ describing the distribution of root hair cells, a harvest functional $\H(μ)$ computes the total amount…
▽ More
This paper studies two classes of variational problems introduced in [7], related to the optimal shapes of tree roots and branches. Given a measure $μ$ describing the distribution of leaves, a sunlight functional $§(μ)$ computes the total amount of light captured by the leaves. For a measure $μ$ describing the distribution of root hair cells, a harvest functional $\H(μ)$ computes the total amount of water and nutrients gathered by the roots. In both cases, we seek a measure $μ$ that maximizes these functionals subject to a rami?ed transportation cost, for transporting nutrients from the roots to the trunk or from the trunk to the leaves. Compared with [7], here we do not impose any a priori bound on the total mass of the optimal measure $μ$, and more careful a priori estimates are thus required. In the unconstrained optimization problem for branches, we prove that an optimal measure exists, with bounded support and bounded total mass. In the unconstrained problem for tree roots, we prove that an optimal measure exists, with bounded support but possibly unbounded total mass. The last section of the paper analyzes how the size of the optimal tree depends on the parameters defining the various functionals.
△ Less
Submitted 13 January, 2020;
originally announced January 2020.
-
On Traffic Flow with Nonlocal Flux: a Relaxation Representation
Authors:
Alberto Bressan,
Wen Shen
Abstract:
We consider a conservation law model of traffic flow, where the velocity of each car depends on a weighted average of the traffic density $ρ$ ahead. The averaging kernel is of exponential type: $w_\varepsilon(s)=\varepsilon ^{-1} e^{-s/\varepsilon}$. By a transformation of coordinates, the problem can be reformulated as a $2\times 2$ hyperbolic system with relaxation. Uniform BV bounds on the solu…
▽ More
We consider a conservation law model of traffic flow, where the velocity of each car depends on a weighted average of the traffic density $ρ$ ahead. The averaging kernel is of exponential type: $w_\varepsilon(s)=\varepsilon ^{-1} e^{-s/\varepsilon}$. By a transformation of coordinates, the problem can be reformulated as a $2\times 2$ hyperbolic system with relaxation. Uniform BV bounds on the solution are thus obtained, independent of the scaling parameter $\varepsilon $. Letting $\varepsilon\to 0$, the limit yields a weak solution to the corresponding conservation law $ρ_t + ( ρv(ρ))_x=0$. In the case where the velocity $v(ρ)= a-bρ$ is affine, using the Hardy-Littlewood rearrangement inequality we prove that the limit is the unique entropy-admissible solution to the scalar conservation law.
△ Less
Submitted 13 December, 2019; v1 submitted 9 November, 2019;
originally announced November 2019.
-
On best constants in $L^2$ approximation
Authors:
Andrea Bressan,
Michael S. Floater,
Espen Sande
Abstract:
In this paper we provide explicit upper and lower bounds on certain $L^2$ $n$-widths, i.e., best constants in $L^2$ approximation. We further describe a numerical method to compute these $n$-widths approximately, and prove that this method is superconvergent. Based on our numerical results we formulate a conjecture on the asymptotic behaviour of the $n$-widths. Finally we describe how the numerica…
▽ More
In this paper we provide explicit upper and lower bounds on certain $L^2$ $n$-widths, i.e., best constants in $L^2$ approximation. We further describe a numerical method to compute these $n$-widths approximately, and prove that this method is superconvergent. Based on our numerical results we formulate a conjecture on the asymptotic behaviour of the $n$-widths. Finally we describe how the numerical method can be used to compute the breakpoints of the optimal spline spaces of Melkman and Micchelli, which have recently received renewed attention in the field of Isogeometric Analysis.
△ Less
Submitted 25 September, 2020; v1 submitted 30 September, 2019;
originally announced September 2019.
-
Competition Models for Plant Stems
Authors:
Alberto Bressan,
Sondre T. Galtung,
Audun Reigstad,
Johanna Ridder
Abstract:
The models introduced in this paper describe a uniform distribution of plant stems competing for sunlight. The shape of each stem, and the density of leaves, are designed in order to maximize the captured sunlight, subject to a cost for transporting water and nutrients from the root to all the leaves. Given the intensity of light, depending on the height above ground, we first solve the optimizati…
▽ More
The models introduced in this paper describe a uniform distribution of plant stems competing for sunlight. The shape of each stem, and the density of leaves, are designed in order to maximize the captured sunlight, subject to a cost for transporting water and nutrients from the root to all the leaves. Given the intensity of light, depending on the height above ground, we first solve the optimization problem determining the best possible shape for a single stem. We then study a competitive equilibrium among a large number of similar plants, where the shape of each stem is optimal given the shade produced by all others. Uniqueness of equilibria is proved by analyzing the two-point boundary value problem for a system of ODEs derived from the necessary conditions for optimality.
△ Less
Submitted 5 September, 2019;
originally announced September 2019.
-
Weighted Irrigation Plans
Authors:
Alberto Bressan,
Qing Sun
Abstract:
We model an irrigation network where lower branches must be thicker in order to support the weight of the higher ones. This leads to a countable family of ODEs, one for each branch, that must be solved by backward induction. Having introduced conditions that guarantee the existence and uniqueness of solutions, our main result establishes the lower semicontinuity of the corresponding cost functiona…
▽ More
We model an irrigation network where lower branches must be thicker in order to support the weight of the higher ones. This leads to a countable family of ODEs, one for each branch, that must be solved by backward induction. Having introduced conditions that guarantee the existence and uniqueness of solutions, our main result establishes the lower semicontinuity of the corresponding cost functional, w.r.t. pointwise convergence of the irrigation plans. In turn, this yields the existence of an optimal irrigation plan, in the presence of these additional weights.
△ Less
Submitted 5 June, 2019;
originally announced June 2019.
-
Globally Optimal Departure Rates for Several Groups of Drivers
Authors:
Alberto Bressan,
Yucong Huang
Abstract:
The first part of this paper contains a brief introduction to conservation law models of traffic flow on a network of roads. Globally optimal solutions and Nash equilibrium solutions are reviewed, with several groups of drivers sharing different cost functions.
In the second part we consider a globally optimal set of departure rates, for different groups of drivers but on a single road. Necessar…
▽ More
The first part of this paper contains a brief introduction to conservation law models of traffic flow on a network of roads. Globally optimal solutions and Nash equilibrium solutions are reviewed, with several groups of drivers sharing different cost functions.
In the second part we consider a globally optimal set of departure rates, for different groups of drivers but on a single road. Necessary conditions are proved, which lead to a practical algorithm for computing the optimal solution.
△ Less
Submitted 16 January, 2019;
originally announced January 2019.
-
Sum-factorization techniques in Isogeometric Analysis
Authors:
A. Bressan,
S. Takacs
Abstract:
The fast assembling of stiffness and mass matrices is a key issue in isogeometric analysis, particularly if the spline degree is increased. We present two algorithms based on the idea of sum factorization, one for matrix assembling and one for matrix-free methods, and study the behavior of their computational complexity in terms of the spline order $p$. Opposed to the standard approach, these algo…
▽ More
The fast assembling of stiffness and mass matrices is a key issue in isogeometric analysis, particularly if the spline degree is increased. We present two algorithms based on the idea of sum factorization, one for matrix assembling and one for matrix-free methods, and study the behavior of their computational complexity in terms of the spline order $p$. Opposed to the standard approach, these algorithms do not apply the idea element-wise, but globally or on macro-elements. If this approach is applied to Gauss quadrature, the computational complexity grows as $p^{d+2}$ instead of $p^{2d+1}$ as previously achieved.
△ Less
Submitted 2 May, 2019; v1 submitted 14 September, 2018;
originally announced September 2018.
-
Approximation in FEM, DG and IGA: A Theoretical Comparison
Authors:
Andrea Bressan,
Espen Sande
Abstract:
In this paper we compare approximation properties of degree $p$ spline spaces with different numbers of continuous derivatives. We prove that, for a given space dimension, $\smooth {p-1}$ splines provide better a priori error bounds for the approximation of functions in $H^{p+1}(0,1)$. Our result holds for all practically interesting cases when comparing $\smooth {p-1}$ splines with…
▽ More
In this paper we compare approximation properties of degree $p$ spline spaces with different numbers of continuous derivatives. We prove that, for a given space dimension, $\smooth {p-1}$ splines provide better a priori error bounds for the approximation of functions in $H^{p+1}(0,1)$. Our result holds for all practically interesting cases when comparing $\smooth {p-1}$ splines with $\smooth {-1}$ (discontinuous) splines. When comparing $\smooth {p-1}$ splines with $\smooth 0$ splines our proof covers almost all cases for $p\ge 3$, but we can not conclude anything for $p=2$. The results are generalized to the approximation of functions in $H^{q+1}(0,1)$ for $q<p$, to broken Sobolev spaces and to tensor product spaces.
△ Less
Submitted 12 February, 2019; v1 submitted 13 August, 2018;
originally announced August 2018.
-
Approximation of Sweeping Processes and Controllability for a Set Valued Evolution
Authors:
Alberto Bressan,
Marco Mazzola,
Khai T. Nguyen
Abstract:
We consider a controlled evolution problem for a set $Ω(t)\in\mathbb{R}^d$, originally motivated by a model where a dog controls a flock of sheep. Necessary conditions and sufficient conditions are given, in order that the evolution be completely controllable. Similar techniques are then applied to the approximation of a sweeping process. Under suitable assumptions, we prove that there exists a co…
▽ More
We consider a controlled evolution problem for a set $Ω(t)\in\mathbb{R}^d$, originally motivated by a model where a dog controls a flock of sheep. Necessary conditions and sufficient conditions are given, in order that the evolution be completely controllable. Similar techniques are then applied to the approximation of a sweeping process. Under suitable assumptions, we prove that there exists a control function such that the corresponding evolution of the set $Ω(t)$ is arbitrarily close to the one determined by the sweeping process.
△ Less
Submitted 23 May, 2018;
originally announced May 2018.
-
Vanishing Viscosity Solutions for Conservation Laws with Regulated Flux
Authors:
Alberto Bressan,
Graziano Guerra,
Wen Shen
Abstract:
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 parabo…
▽ More
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.
△ Less
Submitted 4 May, 2018;
originally announced May 2018.
-
Local approximation operators on box meshes
Authors:
Andrea Bressan,
Tom Lyche
Abstract:
This paper analyzes the approximation properties of spaces of piece-wise tensor product polynomials over box meshes with a focus on application to IsoGeometric Analysis (IGA). The errors are measured in Lebesgue norms. Estimates of different types are considered: local and global, with full or reduced Sobolev seminorms. Attention is also paid to the dependence on the degree and exponential converg…
▽ More
This paper analyzes the approximation properties of spaces of piece-wise tensor product polynomials over box meshes with a focus on application to IsoGeometric Analysis (IGA). The errors are measured in Lebesgue norms. Estimates of different types are considered: local and global, with full or reduced Sobolev seminorms. Attention is also paid to the dependence on the degree and exponential convergence is proved for the approximation of analytic functions.
△ Less
Submitted 20 December, 2018; v1 submitted 22 March, 2018;
originally announced March 2018.
-
On the Optimal Shape of Tree Roots and Branches
Authors:
Alberto Bressan,
Qing Sun
Abstract:
This paper introduces two classes of variational problems, determining optimal shapes for tree roots and branches. Given a measure $μ$, describing the distribution of leaves, we introduce a sunlight functional $§(μ)$ computing the total amount of light captured by the leaves. On the other hand, given a measure $μ$ describing the distribution of root hair cells, we consider a harvest functional…
▽ More
This paper introduces two classes of variational problems, determining optimal shapes for tree roots and branches. Given a measure $μ$, describing the distribution of leaves, we introduce a sunlight functional $§(μ)$ computing the total amount of light captured by the leaves. On the other hand, given a measure $μ$ describing the distribution of root hair cells, we consider a harvest functional $\H(μ)$ 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 $§$ and $\H$, together with a priori estimates on the support of optimal distributions.
△ Less
Submitted 2 March, 2018;
originally announced March 2018.
-
On finite time BV blow-up for the p-system
Authors:
Alberto Bressan,
Geng Chen,
Qingtian Zhang
Abstract:
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 sol…
▽ More
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.
△ Less
Submitted 9 October, 2017;
originally announced October 2017.
-
On the Controllability of Lagrangian Systems by Active Constraints
Authors:
Alberto Bressan,
Zipeng Wang
Abstract:
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 diffe…
▽ More
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.
△ Less
Submitted 15 February, 2017;
originally announced February 2017.
-
Growth Models for Tree Stems and Vines
Authors:
Alberto Bressan,
Michele Palladino,
Wen Shen
Abstract:
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 ea…
▽ More
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.
△ Less
Submitted 23 January, 2017;
originally announced January 2017.
-
A Stochastic Model of Optimal Debt Management and Bankruptcy
Authors:
Alberto Bressan,
Antonio Marigonda,
Khai T. Nguyen,
Michele Palladino
Abstract:
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 d…
▽ More
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^*$?.
△ Less
Submitted 23 September, 2016; v1 submitted 19 September, 2016;
originally announced September 2016.
-
A Model of Controlled Growth
Authors:
Alberto Bressan,
Marta Lewicka
Abstract:
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.…
▽ More
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,α}$ boundary, our main result establishes the local existence and uniqueness of a classical solution, up to a rigid motion.
△ Less
Submitted 29 July, 2017; v1 submitted 30 August, 2016;
originally announced August 2016.
-
The Riemann Solver for Traffic Flow at an Intersection with Buffer of Vanishing Size
Authors:
Alberto Bressan,
Anders Nordli
Abstract:
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 paramete…
▽ More
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.
△ Less
Submitted 10 December, 2015;
originally announced December 2015.
-
Lipschitz Metrics for a Class of Nonlinear Wave Equations
Authors:
Alberto Bressan,
Geng Chen
Abstract:
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 stru…
▽ More
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.
△ Less
Submitted 20 June, 2015;
originally announced June 2015.
-
Structurally Stable Singularities for a Nonlinear Wave Equation
Authors:
Alberto Bressan,
Tao Huang,
Fang Yu
Abstract:
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 an…
▽ More
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.
△ Less
Submitted 30 March, 2015;
originally announced March 2015.
-
No BV bounds for approximate solutions to p-system with general pressure law
Authors:
Alberto Bressan,
Geng Chen,
Qingtian Zhang,
Shengguo Zhu
Abstract:
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 to…
▽ More
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.
△ Less
Submitted 9 February, 2015;
originally announced February 2015.
-
Generic Regularity of Conservative Solutions to a Nonlinear Wave Equation
Authors:
Alberto Bressan,
Geng Chen
Abstract:
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…
▽ More
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.
△ Less
Submitted 9 February, 2015;
originally announced February 2015.
-
Unique Conservative Solutions to a Variational Wave Equation
Authors:
Alberto Bressan,
Geng Chen,
Qingtian Zhang
Abstract:
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ölder 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 c…
▽ More
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ölder 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})$.
△ Less
Submitted 7 November, 2014;
originally announced November 2014.
-
Representation of Dissipative Solutions to a Nonlinear Variational Wave Equation
Authors:
Alberto Bressan,
Tao Huang
Abstract:
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 ar…
▽ More
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.
△ Less
Submitted 4 July, 2014;
originally announced July 2014.