Archivio 2021 170 seminari

Optimization problems subject to PDE constraints form a mathematical tool that can be applied to a wide range of scientific processes, including fluid flow control, medical imaging, biological and chemical processes, and many others. These problems involve minimizing some function arising from a physical objective, while obeying a system of PDEs which describe the process. It is necessary to obtain accurate solutions to such problems within a reasonable CPU time, in particular for time-dependent problems, for which the “all-at-once” solution can lead to extremely large linear systems. In this talk we consider Krylov subspace methods to solve such systems, accelerated by fast and robust preconditioning strategies. A key consideration is which time-stepping scheme to apply — much work to date has focused on the backward Euler scheme, as this method is stable and the resulting systems are amenable to existing preconditioners, however this leads to linear systems of even larger dimension than those obtained when using other (higher-order) methods. We will summarise some recent advances in addressing this challenge, including a new preconditioner for the more difficult linear systems obtained from a Crank-Nicolson discretization, and a Newton-Krylov method for nonlinear PDE-constrained optimization. At the end of the talk we plan to discuss some recent developments in the preconditioning of multiple saddle-point systems, specifically positive definite preconditioners which may be applied within MINRES, which may find considerable utility for solving optimization problems as well as other applications. This talk is based on work with Stefan Güttel (University of Manchester), Santolo Leveque (University of Edinburgh), and Andreas Potschka (TU Clausthal).
Le classi di Jordan sono state introdotte da Borho e Kraft nel loro studio delle sheet per algebre di Lie semisemplici. Sono le classi di equivalenza di elementi in un'algebra di Lie che hanno stessa decomposizione di Jordan, o, equivalentemente di elementi che hanno stabilizzatori (per l'azione aggiunta) coniugati tra loro. Sono localmente chiuse, irriducibili, lisce, e le loro chiusure danno luogo ad una stratificazione finita. La stessa costruzione può essere adattata per definire le classi di Jordan in gruppi algebrici riduttivi: la stratificazione che ne risulta compare nello studio di Lusztig dei fasci carattere. In collaborazione con Ambrosio ed Esposito abbiamo osservato che localmente le chiusure di classi di Jordan nel gruppo si comportano come chiusure di classi di Jordan in un'opportuna algebra di Lie. Un analogo di classe di Jordan per algebre di Lie Z_2-graduate è stato introdotto da Tauvel e Yu e le chiusure sono state studiate da Bulois ed Hivert: si perdono alcune delle caratteristiche dei casi precedenti ma il quadro complessivo è ancora chiaro. Motivato dallo studio della modalità per azioni di gruppi, Popov ha recentemente introdotto le classi di Jordan anche per algebre di Lie ciclicamente graduate. In collaborazione con Esposito e Santi abbiamo fornito una descrizione geometrica locale delle loro chiusure, mostrando in particolare che anche in questo caso la chiusura delle classi di Jordan è un'unione di classi. Con una serie di esempi mostreremo affinità e divergenze tra i vari contesti e le situazioni nelle quali la partizione in classi di Jordan ha un ruolo importante.
Ottobre
dal giorno
11/10/2021
al giorno
13/10/2021
Daniela di Serafino
Sparse Approximations with Interior Point Methods
analisi numerica
Large-scale optimization problems that seek sparse solutions have become ubiquitous. They are routinely solved with various specialized first-order methods. Although such methods are often fast, they usually struggle with not-so-well conditioned problems. In this talk, specialized variants of an interior point-proximal method of multipliers are proposed and analyzed for problems of this class. Computational experience on a variety of problems, namely, multi-period portfolio optimization, classification of data coming from functional Magnetic Resonance Imaging, restoration of images corrupted by Poisson noise, and classification via regularized logistic regression, provides substantial evidence that interior point methods, equipped with suitable linear algebra, can offer a noticeable advantage over first-order approaches. This is joint work with V. De Simone and M. Viola (University of Campania "L. Vanvitelli", Italy) and with J. Gondzio and S. Pougkakiotis (University of Edinburgh, UK).
Ottobre
dal giorno
11/10/2021
al giorno
13/10/2021
Silvia Gazzola
Hybrid projection methods for large-scale linear inverse problems
analisi matematica
Inverse problems are ubiquitous in many areas of Science and Engineering and, once discretized, they lead to ill-conditioned linear systems, often of huge dimensions: regularization consists in replacing the original system by a nearby problem with better numerical properties, in order to find a meaningful approximation of its solution. After briefly surveying some standard regularization methods, both iterative (such as many Krylov methods) and direct (such as Tikhonov method), this talk will introduce a recent class hybrid projection methods, which merge an iterative and a direct approach to regularization. In particular, strategies for choosing the regularization parameter and the regularization matrix will be emphasized, eventually leading to the computation of approximate solutions of Tikhonov problems involving a regularization term expressed in a p-norm.
The three-dimensional reconstruction of an object is an interesting topic with many applications in different fields and has attracted several researchers. The applications range goes from the biomedical 3D reconstruction of human tissues to the approximation of the surface of astronomical objects, from archeology for the digitization of artistic works to the recent development of 3D printing. The first being interested in this problem were some opticians in the Fifties-Sixties. Afterwards, B.K.P. Horn first formulated the Shape-from-Shading (SfS) problem for a single gray-level image of the object. The goal was to get the 3D surface represented in the input image solving a partial differential equation or a variational problem. This problem gave rise to an expansion in the field of mathematics and some researchers tried to prove the well-posedness in the framework of weak solutions. The first works of Lions, Rouy and Tourin in the early 90s inserted the SfS problem in the context of the viscosity solutions frameworks, hence in a much more theoretical area. In this seminar I will start dealing with the orthographic SfS problem with Lambertian reflectance model, the classical and simplest setup for this ill-posed problem that can be modeled by first order Hamilton-Jacobi equations. During the seminar I will briefly introduce some notions of Hamilton-Jacobi equations, viscosity solutions and other ingredients necessary to understand the problem in a general setting. I will continue exploring some non-Lambertian reflectance models and we will see how it is possible to derive a well-posed problem adding information in a natural way. Finally, I will talk about the more recent Shape-from-Polarization problem and the advantages of it with respect to the SfS.
Settembre
dal giorno
16/09/2021
al giorno
18/09/2021
Michiel Bertsch
Nonlinear diffusion and the dynamics of glasslike fluids
analisi matematica
Solutions of degenerate nonlinear diffusion equations such as the porous media equation have the property of finite speed of propagation of disturbances. In this talk we review a class of more degenerate equations characterized by vanishing speed of propagation. We explain how an example of such an equation arises as the singular limit of a parabolic equation describing the dynamics of glasslike fluids. This naturally leads to the problem of finding asymptotic formulas for the velocities of diffuse free boundaries. This is work in progress with Roberto Benzi (Dept. Physics, Univ. of Roma Tor Vergata) and Francesco Deangelis (GSSI, L'Aquila).
Settembre
dal giorno
16/09/2021
al giorno
18/09/2021
Thierry Coulhon
Regularizing effects of non-linear semigroups
analisi matematica
We shall present a summary of a forthcoming book with Daniel Hauer (University of Sydney). The aim of this monograph is to introduce natural and simple functional analytic methods to deduce $L^1-L^\infty$ regularization estimates on nonlinear semigroups from natural Gagliardo-Nirenberg inequalities satisfied by their infinitesimal generator. This enables one to treat in an optimal and unified way a wealth of examples, including the p-Laplace operator, the porous medium operator, as well as variations and combinations of them.
Settembre
dal giorno
16/09/2021
al giorno
18/09/2021
Italo Capuzzo Dolcetta
Invariant cones for linear elliptic systems with gradient couplings
analisi matematica
I will discuss firs the validity of the weak Maximum Principle (wMP) for vector functions u = (u1, .., um) satisfying systems of the form Au + Cu ≥ 0 in a bounded open set Ω of Rn where A is a diagonal matrix of linear degenerate second order elliptic operators and C is a cooperative matrix. Next some counterexamples to the validity of (wMP) are discussed when non diagonal couplings in first order partial derivatives of the ui appear in the system. In this more general setting I will show, through a suitable reduction to a nonlinear scalar equation of Bellman type, that some algebraic condition on the structure of gradient couplings and a cooperativity condition on the matrix of zero order couplings guarantee the existence of invariant cones in the sense of Weinberger.
Settembre
dal giorno
16/09/2021
al giorno
18/09/2021
Luigi Ambrosio
On the Gamma-convergence of action functionals induced by gradients of convex functions
analisi matematica
We investigate the asymptotic behaviour, with respect to Gamma-convergence, of action functionals made up of a kinetic term and a term induced by the the gradient of a convex function. The initial motivation has been the derivation of continuous models from particle systems, but the result has an independent interest. Joint works with A.Baradat, Y.Brenier and C.Brena.
Settembre
dal giorno
16/09/2021
al giorno
18/09/2021
Zoltan Balogh
SHARP ISOPERIMETRIC AND SOBOLEV INEQUALITIES IN SPACES WITH NONNEGATIVE RICCI CURVATURE
analisi matematica
By using optimal mass transport theory we prove a sharp isoperimetric inequal- ity in CD(0, N ) metric measure spaces assuming an asymptotic volume growth at infinity. Our result extends recently proven isoperimetric inequalities for normed spaces and Riemannian manifolds to a nonsmooth framework. As applications of the isoperimetric inequality, we establish Sobolev and Rayleigh-Faber-Krahn inequalities with explicit sharp constants in Riemannian manifolds with nonnegative Ricci curvature; here we use appropriate symmetrization techniques and optimal volume non-collapsing properties. The equality cases in the latter inequalities are also characterized by stating that sufficiently smooth, nonzero extremal functions exist if and only if the Riemannian manifold is isometric to the Euclidean space.
Settembre
dal giorno
16/09/2021
al giorno
18/09/2021
Xavier Cabre
Stable solutions to semilinear elliptic equations are smooth up to dimension 9
analisi matematica
The regularity of stable solutions to semilinear elliptic PDEs has been studied since the 1970's. In dimensions 10 and higher, there exist singular stable energy solutions. In this talk I will describe a recent work in collaboration with Figalli, Ros-Oton, and Serra, where we prove that stable solutions are smooth up to the optimal dimension 9. This answers to an open problem posed by Brezis in the mid-nineties concerning the regularity of extremal solutions to Gelfand-type problems.
Settembre
dal giorno
16/09/2021
al giorno
18/09/2021
Benoit Perthame
Structured equations in biology
analisi matematica
Models arising in biology are often written in terms of Ordinary Differential Equations. The celebrated paper of Kermack-McKendrick (1927), founding mathematical epidemiology, showed the necessity to include parameters in order to describe the state of the individuals as time elapsed after infection. During the 70s, many mathematical studies where developed when equations are structured by age, size or a physiological trait. The talk will present some structured equations, show that a universal relative entropy structure is available in the linear case. In the nonlinear cases it might be that periodic solutions occur, which can be interpreted, e.g., as network activity in the neuroscience. When the equations are conservation laws, a variant of the Monge-Kantorovich distance also gives a general control of solutions.
Settembre
dal giorno
16/09/2021
al giorno
18/09/2021
Pierre Pansu
L1 analysis and related geometry
analisi matematica
Although the Laplacian cannot be inverted on L1, the exterior differential can sometimes be inverted on L1. This discovery, due to Bourgain, Brezis et al. in the early 2000's, can be explained in geometric terms. Such results have generalizations in a Heisenberg group setting. This is joint work with A. Baldi, B. Franchi and F. Tripaldi.
Settembre
dal giorno
16/09/2021
al giorno
18/09/2021
Francesco Serra Cassano
Variational convergences for integral functionals and PDE depending on vector fields
analisi matematica
Settembre
dal giorno
16/09/2021
al giorno
18/09/2021
Juan Manfredi
NATURAL p-MEANS FOR THE p-LAPLACIAN IN EUCLIDEAN SPACE AND THE HEISENBERG GROUP
analisi matematica
We prove uniform convergence in Lipschitz domains in $\mathbb{R}^n$ (and in $C^{1,1}$ domains in the Heisenberg group) of approximations to $p$-harmonic functions obtained using the natural $p$-means introduced by Ishiwata, Magnanini, and Wadade in 2017. This talk is based in joint work with Andr\'as Domokos and Diego Ricciotti (Sacramento) and Bianca Stroffolini (Naples)
Settembre
dal giorno
16/09/2021
al giorno
18/09/2021
Cristian Gutierrez
L∞-estimates in optimal transport for non quadratic costs
analisi matematica
This research originates from recent results by M. Goldman and F. Otto concerning regularity of optimal transport maps for the quadratic cost. We consider cost functions having the form c(x, y) = h(x − y), where h is positively homogeneous of degree p ≥ 2 and h ∈ C2(Rn). A mapping T : Rn → Rn is c-monotone if c(Tx,x) + c(Ty,y) ≤ c(Tx,y) + c(Ty,x). Using Green’s representation formulas, if T is c-monotone, we prove local L∞- estimates of Tx−x in terms of Lp-averages of Tx−x. From this we deduce estimates for the interpolating maps between T and Id, and when T is optimal, L∞-estimates of T −1x − x. As a consequence of the technique, we also obtain a.e. differentiability of monotone maps. This is joint work with Annamaria Montanari.
Settembre
dal giorno
16/09/2021
al giorno
18/09/2021
Irina Markina
On exceptional families of measures
analisi matematica
It is an ongoing project to study collections of measures that are negligible in a sense of ``modules". The idea is originated in complex analysis as ``a conformal module of a family of curves" in looking for an invariant object under conformal transformations on the complex plane. The notion is closely related to the potential theory, certain capacity, and Hausdorff measure. Later the definition of the module was successfully applied to the nonlinear potential theory and quasiconformal analysis in a wider sense in Euclidean spaces. B. Fuglede, by studying the completion of functional spaces, generalized the notion of the module of a family of curves to the module of a family of measures. The arc length of a curve was thought of as a measure. A collection of measures is exceptional if the corresponding module vanishes. In the talk, I will remind examples of exceptional measures in Euclidean space. We aim to find exceptional families of measures on Carnot groups, related to geometric objects such as "intrinsic graphs". It leads to the notion of a Grassmannian on specific Carnot groups.
Settembre
dal giorno
16/09/2021
al giorno
18/09/2021
Ermanno Lanconelli
Polarity measures and their rigidity properties.
analisi matematica
Let $\mu$ be a measure concentrated on a domain $ D \subset \mathbb {R}^N$ , and let $ x_0 \in D$. Denote by $ \Gamma$ the fundamental solution of the Laplacian, and by $\Gamma_{\mu}$ the Newtonian potential of $\mu$. We say that $\mu$ is a polarity measure for $D$ at $x_0$ if $\Gamma_{\mu} = \Gamma (x_0 - x)$ for every $x$ in the complementary of $D$. If we also have $\Gamma_{\mu} < \Gamma (x_0 - x)$ for every $x \in D$ then we say that $\mu$ is a strong polarity measure for $D$ at $x_0$. In the present talk we first recall the following results: A. Every sufficiently smooth domain supports a polarity measure at an arbitrarily given point. B. Every strong polarity measure characterizes its supporting domain. Then we show how to extend A and B to the general context of the hypoellipitic semi-elliptic linear second order PDEs. All the results we present have been obtained in collaboration with Giovanni Cupini.
Settembre
dal giorno
16/09/2021
al giorno
18/09/2021
Fulvio Ricci
Recent developments on maximal restriction operators
analisi matematica
Suppose that, for a surface $S\subset\bR^n$ with (weighted) surface measure $\sigma$ and for some $p,q$ with $p\in(1,2)$, the Fourier restriction operator $\cR:f\longmapsto \widehat f_S$ satisfies the inequality $$ \|\mathcal R f\|_{L^q(S,\sigma)}\le C\|f\|_{L^p(\mathbb R^n)}\ ,\qquad \forall f\in\cS(\mathbb R^n)\ . $$ Then extendability of $\mathcal R$ to all of $L^p(\mathbb R^n)$ indicates, heuristically, that, for general $f\in L^p(\mathbb R^n)$, $\widehat f$ can be assigned values on $S$, despite the fact that it is only defined a.e. The notion of ``maximal restriction operator'' has been introduced in a paper of 2019 by D.~M\"uller, J.~Wright and myself, for the purpose of giving measure-theoretic ground to this statement. In this talk I give a precise presentation of the problem, the improvements of our original result by various authors and some of the open problems.
Settembre
dal giorno
16/09/2021
al giorno
18/09/2021
Raul Serapioni
Regular and irregular solutions of degenerate elliptic equations: a glance at vintage mathematics.
analisi matematica
The problem of Holder regularity of a variational solutions u = u(x) of a degenerate uniformly elliptic second order equations as (1) \sum_{i=1}^n Di(w(x)Diu(x)) = 0; x\in\Omega\subset R^n has been addressed since the beginning of the seventies. Now it is well known that if w belongs to the Muckenhaupt class A_2 then variational solutions of (1) are Holder continuous. On the other side the necessity of the assumption w\in A_2, or of similar structural assumptions on the weights, is far from being well understood. The simpler question of the necessity/sufficiency of quantitative assumptions on w and 1/w, even if better understood, is not yet completely settled.
Let G be a permutation group acting on a finite set Omega. A subset B of Omega is called a base for G if the pointwise stabilizer of B in G is trivial. In the 19th century, bounding the order of a finite primitive permutation group G was a problem that attracted a lot of attention. Early investigations of bases then arose because such a problem reduces to that of bounding the minimal size of a base of G. Some other far- reaching applications across Pure Mathematics led the study of the base size to be a crucial area of current research in permutation groups. In the first part of the talk, we will investigate some of these applications and review some results about base size. We will present a recent improvement of a famous estimation due to Liebeck that estimates the base size of a primitive permutation group in terms of its degree. In the second part of the talk, we will define the concept of irredundant bases of G and the concept of IBIS groups. Whereas bases of minimal size have been well studied, irredundant bases and IBIS groups have not yet received a similar degree of attention. Indeed, Cameron and Fon-Der-Flaas, already in 1995, defined such groups and proposed to classify some meaningful families. But only this year, a systematic investigation of primitive permutation IBIS groups has been started. We will discuss how we reduced the classification of primitive IBIS groups to the almost simple groups and affine groups. Eventually, we will conclude by mentioning recent advances towards a complete classification.
The dynamics of an inviscid and incompressible fluid flow on a Riemannian manifold is governed by the Euler equations. Recently, Tao launched a programme to address the global existence problem for the Euler and Navier Stokes equations based on the concept of universality. Inspired by this proposal, we show that the stationary Euler equations exhibit several universality features, In the sense that, any non-autonomous flow on a compact manifold can be extended to a smooth stationary solution of the Euler equations on some Riemannian manifold of possibly higher dimension. These results can be viewed as lending support to the intuition that solutions to the Euler equations can be extremely complicated in nature. A key point in the proof is looking at the h-principle in contact geometry through a contact mirror, unveiled by Sullivan, Etnyre and Ghrist more than two decades ago. We end up this talk addressing an apparently different question: What kind of physics might be non-computational? Using the former universality result, we can establish the Turing completeness of the steady Euler flows, i.e., there exist solutions that encode a universal Turing machine and, in particular, these solutions have undecidable trajectories. But, in view of the increase of dimension yielded by our proof, the question is: can this be done in dimension 3? We will prove the existence of Turing complete fluid flows on a 3-dimensional geometric domain. Our novel strategy uses the computational power of symbolic dynamics and the contact mirror again. This talk is based on joint work with Robert Cardona and Fran Presas (arXiv:1911.01963 and arXiv:2012.12828).
Numerical first-order methods are the most suitable choice for solving large-scale nonlinear optimization problems which model many real life applications. Among these approaches, gradient methods have widely proved their effectiveness in solving challenging unconstrained and constrained problems arising in machine learning, compressive sensing, image processing and other areas. These methods became extremely popular since the work by Barzilai and Borwein  (BB) (1988), which showed how a suitable choice of the steplength can significantly accelerate the classical Steepest Descent method. It is well-known that the performance of gradient methods based on the BB steplength does not depend on the decrease of the objective function at each iteration but relies on the relationship between the steplengths used and the eigenvalues of the average Hessian matrix; hence BB based methods are also denoted as Spectral Gradient methods. The first part of this seminar will be devoted to a review of spectral gradient methods for unconstrained optimization while the second part will focus on recent advances on the extension of these methods to the solution of large nonlinear systems of equations, the so-called Spectral Residual methods. These methods are derivative-free, low-cost per iteration and are particularly suitable when the Jacobian matrix of the residual function is not available analytically or its computation is not relatively easy. In this framework, numerical experience will be presented on sequences of nonlinear systems arising from rolling contact models which play a central role in many important applications, such as rolling bearings and wheel-rail interaction.