Past PDE CDT Lunchtime Seminar

We investigate the long-time behaviour of solutions to a nonlocal partial differential equation on smooth Riemannian manifolds of bounded sectional curvature. The equation models self-collective behaviour with intrinsic interactions that are modeled by an interaction potential. Without the diffusion term, we consider attractive interaction potentials and establish sufficient conditions for a consensus state to form asymptotically. In addition, we quantify the approach to consensus, by deriving a convergence rate for the diameter of the solution’s support. With the diffusion term, the attractive interaction and the diffusion compete. We provide the conditions of the attractive interaction for each part to win.

The Camassa--Holm (CH) equation is a nonlocal equation that manifests supercritical behaviour in ``wave-breaking" and non-uniqueness. In this talk, I will discuss the existence of global (dissipative weak martingale) solutions to the CH equation with multiplicative, gradient type noise, derived as an inviscid limit. The goal of the talk is twofold. The stochastic CH equation will be used to illustrate aspects of a stochastic compactness and renormalisation method which is popularly used to derive well-posedness and continuous dependence results in SPDEs. I shall also discuss how a lack of temporal compactness introduces fundamental difficulties in the case of the stochastic CH equation.

This talk is based on joint works with L. Galimbert and H. Holden, both at NTNU, and with K.H. Karlsen at the University of Oslo. 

We will show asymptotic analysis for hydrodynamic system, as Navier-Stokes-Fourier system, as a useful tool in in the situation when certain parameters in the system – called characteristic numbers – vanish or become infinite. The choice of proper scaling, namely proper system of reference units, the parameters determining the behaviour of the system under consideration allow to eliminate unwanted or unimportant for particular phenomena modes of motion. The main goal of many studies devoted to asymptotic analysis of various physical systems is to derive a simplified set of equations - simpler for mathematical or numerical analysis. Such systems may be derived in a very formal way, however we will concentrate on rigorous mathematical analysis. I will concentrate on low Mach number limits with so called ill-prepared data and I will present some results which concerns passage from compressible to incompressible models of fluid flow emphasising difficulties characteristic for particular problems. In particular we will discuss Navier-Stokes-Fourier system on varying domains, a multi-scale problem for viscous heat-conducting fluids in fast rotation and the incompressible limit of compressible finitely extensible nonlinear bead-spring chain models for dilute polymeric fluids.

This talk is motivated by work on the existence theory for viscoelasticity of Kelvin-Voigt type with non-convex stored energies (joint with K. Koumatos (U. of Sussex), C. Lattanzio and S. Spirito (U. of LAquila)), which shows propagation of H1-regularity for the deformation gradient of weak solutions for semiconvex stored energies. It turns out that weak solutions with deformation gradient in H1 are in fact unique in two-space dimensions, providing a striking analogy to corresponding results in the theory of 2D Euler equations with bounded vorticity.

While weak solutions still exist for initial data in L2, oscillations on the deformation gradi- ent can now persist and propagate in time. This can be seen via a counterexample indicating that for non-monotone stress-strain relations in 1-d oscillations of the strain lead to solutions with sustained oscillations. The existence of sustained oscillations in hyperbolic-parabolic system is then studied in several examples motivated by viscoelasticity and thermoviscoelas- ticity. Sufficient conditions for persistent oscillations are developed for linear problems, and examples in some nonlinear systems of interest. In several space dimensions oscillatory exam- ples are associated with lack of rank-one convexity of the stored energy. Nonlinear examples in models with thermal effects are also developed.

The formulation of Mean Field Games (MFG) via partial differential equations typically requires continuous differentiability of the Hamiltonian in order to determine the advective term in the Kolmogorov--Fokker--Planck equation for the density of players. However, in many cases of practical interest, the underlying optimal control problem may exhibit bang-bang controls, which typically lead to nondifferentiable Hamiltonians. In this talk we will present results on the analysis and numerical approximation of stationary second-order MFG systems for the general case of convex, Lipschitz, but possibly nondifferentiable Hamiltonians. In particular, we will propose a generalization of the MFG system as a Partial Differential Inclusion (PDI) based on interpreting the derivative of the Hamiltonian in terms of subdifferentials of convex functions. We present results that guarantee the existence of unique weak solutions to the stationary MFG PDI under a monotonicity condition similar to one that has been considered previously by Lasry and Lions. Moreover, we will propose a monotone finite element discretization of the weak formulation of the MFG PDI, and present results that confirm the strong H^1-norm convergence of the approximations to the value function and strong L^q-norm convergence of the approximations to the density function. The performance of the numerical method will be illustrated in experiments featuring nonsmooth solutions. This talk is based on joint work with my supervisor Iain Smears.

Irreversible processes are accompanied by an increase in the internal entropy of a continuum, and as such the entropy production function is fundamental in determining the overall state of the system. In this talk, it will be shown that the entropy production function can be utilized for a variational analysis of certain dissipative continua in two different ways. Firstly, a novel unified Lagrangian-Hamiltonian formalism is constructed giving phase space extra structure, and applied to the study of fluid flow and brittle fracture.  Secondly, a maximum entropy production principle is presented for simple bodies and its implications to the study of fluid flow discussed. 

Second-order regularity results are established for solutions to elliptic equations and systems with the principal part having a Uhlenbeck structure and square-integrable right-hand sides. Both local and global estimates are obtained. The latter apply to solutions to homogeneous Dirichlet problems under minimal regularity assumptions on the boundary of the domain. In particular, if the domain is convex, no regularity of its boundary is needed. A critical step in the approach is a sharp pointwise inequality for the involved elliptic operator. This talk is based on joint investigations with A.Kh.Balci, L.Diening, and V.Maz'ya.

Localisation is a powerful tool in proving and analysing various geometric inequalities, including isoperimertic inequality in the context of metric measure spaces. Its multi-dimensional generalisation is linked to optimal transport of vector measures and vector-valued Lipschitz maps. I shall present recent developments in this area: a partial affirmative answer to a conjecture of Klartag concerning partitions associated to Lipschitz maps on Euclidean space, and a negative answer to another conjecture of his concerning mass-balance condition for absolutely continuous vector measures. During the course of the talk I shall also discuss an intriguing notion of ghost subspaces related to the above mentioned partitions. 

In General Relativity, an impulsive gravitational wave is a localized and singular solution of the 

Einstein equations modeling the spacetime distortions created by a strongly gravitating source.
I will present a comprehensive theory allowing for ternary interactions of such impulsive gravitational waves in translation-symmetry, offering the first examples of such an interaction.  

The proof combines new techniques from harmonic analysis, Lorentzian geometry, and hyperbolic PDEs that are helpful to treat highly anisotropic low-regularity questions beyond the considered problem.  

This is joint work with Jonathan Luk.

On a smooth Riemannian manifold, the uniqueness of a geodesic given initial conditions follows from standard ODE theory. This is known to fail in the setting of RCD(K,N) spaces (metric measure spaces satisfying a synthetic notion of Ricci curvature bounded below) through an example of Cheeger-Colding. Strengthening the assumption a little, one may ask if two geodesics which agree for a definite amount of time must continue on the same trajectory. In this talk, I will show that this is true for RCD(K,N) spaces. In doing so, I will generalize a well-known result of Colding-Naber concerning the Hölder continuity of small balls along geodesics to this setting.

We derive quantitatively the weak and strong Harnack inequality for kinetic Fokker--Planck type equations with a non-local diffusion operator for the full range of the non-locality exponents in (0,1).  This implies Hölder continuity.  We give novel proofs on the boundedness of the bilinear form associated to the non-local operator and on the construction of a geometric covering accounting for the non-locality to obtain the Harnack inequalities.  Our results apply to the inhomogeneous Boltzmann equation in the non-cutoff case.

I will discuss the regularity of solutions to quasilinear systems satisfying a Legendre-Hadamard ellipticity condition. For such systems it is known that weak solutions may which fail to be C^1 in any neighbourhood, so we cannot expect a general regularity theory. However if we assume an a-priori regularity condition of the solutions we can rule out such counterexamples. Focusing on solutions to Euler-Lagrange systems, I will present an improved regularity results for solutions whose gradient satisfies a suitable BMO / VMO condition. Ideas behind the proof will be presented in the interior case, and global consequences will also be discussed.

We will consider the viscous Burgers driven by a localised one-dimensional control. The problem is considered in a bounded domain and is supplemented with the Dirichlet boundary condition. We will prove that any solution of the equation in question can be exponentially stabilised. Combining this result with an earlier result on local exact controllability we will show global exact controllability by a localised control. This is a joint work with A. Shirikyan.

We construct smooth flows with surgery that approximate weak mean curvature flows with only spherical and neck-pinch singularities. This is achieved by combining the recent work of Choi-Haslhofer-Hershkovits, and Choi-Haslhofer-Hershkovits-White, establishing canonical neighbourhoods of such singularities, with suitable barriers to flows with surgery. A limiting argument is then used to control these approximating flows. We demonstrate an application of this surgery flow by improving the entropy bound on the low-entropy Schoenflies conjecture.

In this talk I will review some recent results obtained in collaboration with E. Runa and A. Kerschbaum on the one-dimensionality of the minimizers
of a family of continuous local/nonlocal interaction functionals in general dimension. Such functionals have a local term, typically the perimeter or its Modica-Mortola approximation, which penalizes interfaces, and a nonlocal term favouring oscillations which are high in frequency and in amplitude. The competition between the two terms is expected by experiments and simulations to give rise to periodic patterns at equilibrium. Functionals of this type are used  to model pattern formation, either in material science or in biology. The difficulty in proving the emergence of such structures is due to the fact that the functionals are symmetric with respect to permutation of coordinates, while in more than one space dimensions minimizers are one-dimesnional, thus losing the symmetry property of the functionals. We will present new techniques and results showing that for two classes of functionals (used to model generalized anti-ferromagnetic systems, respectively  colloidal suspensions), both in sharp interface and in diffuse interface models, minimizers are one-dimensional and periodic, in general dimension and also while imposing a nontrivial volume constraint.

The integral of mean curvature squared is a conformal invariant that measures the distance from a given immersion to the standard embedding of a round sphere. Following work of Robert Bryant who showed that all Willmore spheres in the 3-sphere are conformally minimal, Robert Kusner proposed in the early 1980s to use the Willmore energy to obtain an “optimal” sphere eversion, called the min-max sphere eversion.

We will present a method due to Tristan Rivière that permits to tackle a wide variety of min-max problems, including ones about the Willmore energy. An important step to solve Kusner’s conjecture is to determine the Morse index of branched Willmore spheres, and we show that the Morse index of conformally minimal branched Willmore spheres is equal to the index of a canonically associated matrix whose dimension is equal to the number of ends of the dual minimal surface.

Elastic strings are among the simplest one-dimensional continuum bodies and have a rich mechanical and mathematical theory dating back to the derivation of their equations of motion by Euler and Lagrange. In classical treatments, the string is either completely extensible (tensile force produces elongation) or completely inextensible (every segment has a fixed length, regardless of the motion). However, common experience is that a string can be stretched (is extensible), and after a certain amount of tensile force is applied the stretch of the string is maximized (becomes inextensible). In this talk, we discuss a model for these stretch-limited elastic strings, in what way they model elastic behavior, the well-posedness and asymptotic stability of certain simple motions, and (many) open questions.

In this talk, we discuss the feasibility of algorithms based on deep artificial neural networks (DNN) for the solution of high-dimensional PDEs, such as those arising from stochastic control and games. In the first part, we show that in certain cases, DNNs can break the curse of dimensionality in representing high-dimensional value functions of stochastic control problems. We then exploit policy iteration to reduce the associated nonlinear PDEs into a sequence of linear PDEs, which are then further approximated via a multilayer feedforward neural network ansatz. We establish that in suitable settings the numerical solutions and their derivatives converge globally, and further demonstrate that this convergence is superlinear, by interpreting the algorithm as an inexact Newton iteration. Numerical experiments on Zermelo's navigation problem and on consensus control of interacting particle systems are presented to demonstrate the effectiveness of the method. This is joint work with Kazufumi Ito, Christoph Reisinger and Wolfgang Stockinger.

I will discuss regularity properties for solutions of linear second order non-uniformly elliptic equations in divergence form. Assuming certain integrability conditions on the coefficient field, we obtain local boundedness and validity of Harnack inequality. The assumed integrability assumptions are sharp and improve upon classical results due to Trudinger from the 1970s.

As an application of the local boundedness result, we deduce a quenched invariance principle for random walks among random degenerate conductances. If time permits I will discuss further regularity results for nonlinear non-uniformly elliptic variational problems.

We consider the homogenization of a Stokes system in a domain having many small random holes. This model mainly arises from problems of solid-fluid interaction (e.g. the flow of a viscous and incompressible fluid through a porous medium). We aim at the rigorous derivation of the homogenization limit both in the Brinkmann regime and in the one of Darcy’s law. In particular, we focus on holes that are distributed according to probability measures that allow for overlapping and clustering phenomena.

We explore one facet of an old problem: the approximation of hyperbolic conservation laws by viscous counterparts. While qualitative convergence results are well-known, quantitative rates for the inviscid limit are less common. In this talk, we consider the simplest case: a one-dimensional scalar strictly-convex conservation law started from "generic" smooth initial data. Using a matched asymptotic expansion, we quantitatively control the inviscid limit up to the time of first shock. We conclude that the inviscid limit has a universal character near the first shock. This is joint work with Sanchit Chaturvedi.

Phase transitions are present in a wide array of systems ranging from traffic to machine learning algorithms. In this talk, we will relate the concept of phase transitions to the convexity properties of the associated thermodynamic energy. Motivated by noisy stochastic gradient descent in supervised learning, we will consider the problem of understanding the thermodynamic limit of exchangeable weakly interacting diffusions (AKA propagation of chaos) from an energetic perspective. The strategy will be to exploit the 2-Wasserstein gradient flow structure associated with the thermodynamic energy in the infinite particle setting. Using this perspective, we will show how the convexity properties of the thermodynamic energy affects the homogenization limit or the stability of the log-Sobolev inequality.

A topology optimization problem for Stokes flow finds the optimal material distribution of a fluid in Stokes flow that minimizes the fluid’s power dissipation under a volume constraint. In 2003, T. Borrvall and J. Petersson [1] formulated a nonconvex optimization problem for this objective. They proved the existence of minimizers in the infinite-dimensional setting and showed that a suitably chosen finite element method will converge in a weak(-*) sense to an unspecified solution. In this talk, we will extend and refine their numerical analysis. In particular, we will show that there exist finite element functions, satisfying the necessary first-order conditions of optimality, that converge strongly to each isolated local minimizer of the problem.

/

Fully nonlinear PDEs involving the eigenvalues of matrix-valued differential operators (such as the Hessian) have been the subject of intensive study over the last few decades, since the seminal work of Caffarelli, Kohn, Nirenberg and Spruck. In this talk I will discuss some recent joint work with Luc Nguyen on the regularity theory for a large class of these equations, with a particular emphasis on a special case known as the sigma_k-Yamabe equation, which arises in conformal geometry. 

[1] T. Borrvall, J. Petersson, Topology optimization of fluids in Stokes flow, International Journal for Numerical Methods in Fluids 41 (1) (2003) 77–107. doi:10.1002/fld.426.

We analyze the steady motion of a viscous incompressible fluid in a two- and three-dimensional channel containing an obstacle through the Navier-Stokes equations under different types of boundary conditions. In the 2D case we take constant non-homogeneous Dirichlet boundary data in a (virtual) square containing the obstacle, and emphasize the connection between the appearance of lift and the unique solvability of Navier-Stokes equations. In the 3D case we consider mixed boundary conditions: the inflow is given by a fairly general datum and the flow is assumed to satisfy a constant traction boundary condition on the outlet. In the absence of external forcing, explicit bounds on the inflow velocity guaranteeing existence and uniqueness of such steady motion are provided after estimating some Sobolev embedding constants and constructing a suitable solenoidal extension of the inlet velocity. In the 3D case, this solenoidal extension is built through the Bogovskii operator and explicit bounds on its Dirichlet norm (in terms of the geometric parameters of the obstacle) are found by solving a variational problem involving the infinity-Laplacian.

The talk accounts for results obtained in collaboration with Filippo Gazzola and Ilaria Fragalà (both at Politecnico di Milano).

Solving kinetic or related models with high-dimensional random parameters has been a challenging problem. In this talk, we will discuss how to employ the bi-fidelity stochastic collocation and choose efficient low-fidelity models in order to solve a class of multi-scale kinetic equations with uncertainties, including the Boltzmann equation, linear transport and the Vlasov-Poisson equation. In addition, some error analysis for the bi-fidelity method based on these PDEs will be presented. Finally, several numerical examples are shown to validate the efficiency and accuracy of the proposed method.