Concentration phenomena in nonlinear partial differential equations.

Lead Research Organisation: University of Bath
Department Name: Mathematical Sciences

Abstract

Broadly speaking, the area of my project is partial differential equations (PDEs). This is the branch of Mathematics which uses the tools of calculus to model phenomena in nature. Indeed, many laws in Physics, Biology, Economics, Social Studies, can be formulated as PDEs.

The study of PDEs is a very broad field within Mathematics and can encompass both theoretical and more applied perspectives. For instance, Euler equations are a set of PDEs that describe how the velocity, pressure and density of a moving fluid are related. These equations neglect the effects of the viscosity which are included in the Navier-Stokes equations. A solution of the Euler equations is therefore only an approximation to a real fluids model. For some problems, like the lift of a thin airfoil at low angle of attack, a solution of the Euler equations provides a good model of reality. For other problems, like the growth of the boundary layer on a flat plate, the Euler equations do not properly model the problem.

My main interest is for the purely mathematical aspects of the study of PDEs. The typical questions that arise in the study of PDEs include: Do solutions of a given equation (theoretically) exist? (If not, our model is not capturing something essential.) Are they stable under perturbations of the initial data? (If not, they may be difficult or impossible to observe in nature.) Do they have some inherent symmetry that reflects the underlying physical or biological phenomena being modeled? (Nature is intrinsically economical, and often the 'simplest' solutions have the most symmetry.) Do the solutions vary smoothly over time and space, or are abrupt changes possible (what mathematicians refer to as formation of singularities in PDEs)?

In this proposal I will address all these questions for specific non-linear PDEs, with main emphasis on the last one: the mathematical analysis of formation of singularities. In many models, static or dynamic in nature, governed by non-linear PDEs, one observes the formation of singularities or some form of concentration of their solutions, as the time-variable or some parameter approaches a limit value. This happens when solutions become concentrated on lower-dimensional sets, or some expressions dependent on the solution become arbitrarily large. From a PDEs' point of view, this phenomenon reflects lack of compactness in the variational formulation of the problem or loss of regularity in the solution set, which is usually related with relevant episodes of the modeled event. Think of the explosion of some substance triggered by a chemical reaction or the appearance of fractures in planes or bridges.

We propose the construction of solutions with singularities for some significant non-linear PDEs, such as for Euler equations for incompressible inviscid fluids, for Ginzburg-Landau model in superconductivity, for sine-Gordon equations, for Keller-Segel model in chemotaxis and for the prescribed mean curvature problem. My aim is to elaborate new refined gluing techniques to carry out these constructions and to derive precise descriptions on why, where and how formation of singularities takes place. My results will be of interest not only in Mathematical Analysis, but also in Geometric Flows, Geometric Partial Differential Equations and Boundary Value Problems for Nonlinear PDE's.

Planned Impact

Partial Differential Equations (PDEs) are key equations for anyone wanting to use mathematics to solve real life problems. Mathematical descriptions of continuous systems are typically phrased in terms of rates of change, or derivatives. Many laws in Physics, Biology, Ecomonics, Social Studies, can be formulated as PDEs. To mention some examples: Newton's laws of motion, non-linear diffusion equations that describe density fluctuations in a material undergoing diffusion, parabolic equation of dissipation type to describe movement of prices in economics, the Fokker-Planck equation for the evolution of individuals in a population, among others. Another interesting feature of PDEs is their universal applicability and their flexibility in allowing us to model and understand changes in different systems.
My project focuses on the theoretical analysis of important classes of PDEs, and on the formation of possible singularities in their solutions. A singularity is a location in space and time where some relevant quantity related to the equation becomes infinity, in a way that does not depend on the chosen coordinate system. Detecting singularities has an important impact in the applications the PDEs are modeling.
As an example, in aerodynamics, a mayor issue is the optimal design of a wing or an airfoil. The designer seeks to optimize the geometric shape of a configuration taking into account the trade-offs between aerodynamic performance, structure weight, and the requirement for internal volume to contain fuel and payload. The subtlety and complexity of fluid flow is such that it is unlikely that repeated trials in an interactive analysis and design procedure can lead to a truly optimum design. Progress toward automatic design has been restricted by the extreme computing costs that might be incurred from brute force numerical optimization. However, useful design methods have been devised for various simplified models, such as two-dimensional wings in inviscid flows, governed by Euler equations, which are part of my research proposal. Theoretical studies of formation of singularity in the vorticity of the fluid's velocity, and of its dynamics, are fundamental to optimize the shape of a wing or an airfoil. Indeed, the aerodynamic properties of the airfoil are linked to the vortices shed at the trailing edge, and understanding these vortices is the key to quantifying lift.
Another example, in mathematical biology, is the movement of cells by chemotaxis in response to a chemical stimulus, which is widely observed in various biological situations, as morphogenesis, bacterial self-organization and inflammatory processes. The first and most studied mathematical model of chemotaxis is the Keller-Segel system of PDEs, one of the topics of my project, whose analysis has helped to understand certain characteristics of bacterial chemotaxis. Even after 40 years of research in this problem, very little is known on formation of singularities, here represented by the cells density becoming arbitrarily large.
My research will also have a strong human resource impact component. I will employ and train a highly skilled post-doctoral research associate (PDRA) who will gain new training and practice in cutting edge energy critical heat equations and in the analysis of blow-up solutions for such problems. Furthermore, the research project will benefit two doctoral students and train them similarly. These doctoral students are funded respectively by a grant from the Royal Society and by EPSRC. The PDRA and students will be part of a much larger community of mathematicians and benefit, through the attendance to seminars and other activities (such as the Integrated Think Tanks organized by the EPSRC Centre for Doctoral Training in Statistical Applied Mathematics SAMBa based at the University of Bath), of a broad exposure to applied, computational and industrial mathematics.
 
Description The main objective of this research project is to analyze some specific non-linear Partial Differential Equations (PDEs), which are used to model phenomena in fluids, chemotaxis, geometry and related problems. The common feature of these PDEs is the formation of singularities in their solutions: when the time-variable or a parameter in the model approaches a limit value, solutions become highly concentrated on lower-dimensional sets (like points, curves, etc), losing smoothness and approaching a singular limit.
One of the fundamental equations in fluids are Euler equations. They describe the evolution of the velocity of an incompressible fluid. The vorticity of the fluid measures its local spinning motion, it is mathematically described as the curl of the fluid's velocity and its value is very high around swirls. Existence of solutions to Euler equations with vorticity highly concentrated around points in the plane is well known. The corresponding problem in dimension 3 is to find a smooth solution to Euler equations whose vorticity is high in tubular neighborhood of curves, or filaments. Only few rigorous results are known, like the vortex rings. One key finding of this research (Work Project 1) is the construction and the precise description of solutions with vorticity concentrated around families of evolving helices: these helices evolve in time following the binormal flow (they move in the direction of the binormal vector, with speed proportional to their curvature) and do not change their form during the evolution (see preprint https://arxiv.org/abs/2007.00606, accepted for pubblication in the journal Calculus of Variations and PDEs in March 2022). Another key finding of this line of research is the construction of traveling and rotation solutions of the generalized inviscid surface-quasi-geostrophic equations (see preprint https://arxiv.org/pdf/2008.12911.pdf, accepted for publication in the Transactions of the American Mathematical Society). These equations are models used in atmospheric and oceanographic flows that describe the balance between the Coriolis force and horizontal pressure forces in winds.
With Shrish Parmeshwar, the PDRA associated to this project, we have found the law of evolution of two vortex pairs for the Euler equation in dimension 2 that travel in opposite direction. We are able to describe in detail the long time behaviour of this configuration and to construct rigorously such solution. We are completing the details of this work in the form of a paper.
A third key finding is related to the construction of shrinking sphere and bubble towers for the critical heat equation (WorkProject 2). It was conjectured that in dimensions greater of equal to 7, the critical heat equation had solutions whose energy became infinity for large values of the time variable (the so-called infinite time blow-up). We could provide this construction and we discovered the existence of solutions with the form of a superposition with alternate signs of singularly scaled Aubin-Talenti solitons (Preprint https://arxiv.org/abs/1905.13500, accepted for publication to Analysis and PDE). There exists also a type II finite time blow-up solution with singularity taking place along an (n-4)-dimensional shrinking sphere (Preprint https://arxiv.org/pdf/2006.00716.pdf). Analogous phenomena have been previously detected for a modified version of 3-dimension Euler equations by Terence Tao. Another key finding related to WorkProject 2 is the construction of infinite-time multiple blow-up for a smooth perturbation of the Yamabe flow on a smooth compact Riemannian manifold, under the validity of the positive mass theorem. The result can be found in the form of a preprint in arXiv:2106.09220.
Key findings in singularity formation associated to geometry related problems are a new doubling construction for the critical Yamabe equation in the whole space. This procedure allows to find sign-changing solutions to Yamabe equation, giving rise to solutions with maximal rank rigidity in odd dimensions. (Preprint https://arxiv.org/abs/2001.03548).
Exploitation Route The techniques developed in the findings described above are new: we developed an original approach to deal with the outer regime far away from the region of high vorticity. This approach has the potential to be applied in related problems, like the long-time dynamics of a pair of vortex-pairs traveling in opposite directions (this is a joint project with the PDRA), or the leapfrogging of two parallel vortex rings. This phenomenon has only been observed in experiments but never detected in a mathematically rigorous way.
Sectors Education

URL https://arxiv.org/abs/2007.0060
 
Description I include the Academic Impact as part of my Pathway to Impact since it forms part of my critical pathway to Economic and Societal Impact. The reason is that it includes reaching academic disciplines beyond my own. A central part of this project is the training of one early career researcher, the PDRA. Training a PDRA, who is an already formed mathematician, is a quite challenging experience, since it is not just bringing on a highly skilled researcher, but committing to the promotion and success of the person's career. I was very happy to hire Shrish Parmeshwar. The interaction with Shrish Parmeshwar is producing striking and innovative ideas in the analysis of long-time dynamics for vortex pairs. I am training two PhD students, who are working on topics closely related to this project. One is studying finite-time blow-up phenomena in the classical Keller-Segel model for chemotaxis-driven motion of cells or bacteria. This topic is strongly related to my Workpackage 2. The other student is studying the dynamics of vortices in the Gross-Pitaevskii equation. The Gross-Pitaevskii equation (GPE, named after Eugene P. Gross and Lev Petrovich Pitaevskii) describes the ground state of a quantum system of identical bosons using the Hartree-Fock approximation and the pseudopotential interaction model. is the Dirac delta-function. This equation is a significant model often used in statistical physics, and has a number of applications from dynamics of a Bose-Einstein condensate to the excitations of the gas cloud.
First Year Of Impact 2020
Sector Education
Impact Types Cultural,Societal,Economic

 
Description One World PDE Seminar 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Study participants or study members
Results and Impact With the mass cancellation of in-person conferences and seminars due to the coronavirus pandemic, research communities need novel ways to stay connected. Inspired by the One World Probability Seminar, the One World PDE Seminar aims to provide such a venue for the PDE community, accessible to as many researchers as possible.
Since March 2020, Miles Wheeler, Angela Pistoia and I have been organizing two 45-minutes talks every week, with a summer break in August, delivered by leading experts in different aspects of the theory of Partial Differential Equations. We have had an average of 150 participants all over the world, with a maximum of 460 participants.
1163 people have signed up for our mailing list, and we have 277 subscribers to our YouTube channel. On YouTube we have 10.6k views totaling 886.2 hours, for an average view duration about 5 minutes. These numbers are unusually high for this type of events in Mathematics.
https://people.bath.ac.uk/mw2319/owpde/
https://www.youtube.com/channel/UCuAjF7px9Tu5huvUgn48SSw/videos
Year(s) Of Engagement Activity 2020
URL http://people.bath.ac.uk/mw2319/owpde/