Shape optimization is a central field in applied mathematics and engineering, with applications ranging from aerodynamics to material design. The objective is to find the optimal shape of a domain that minimizes a given criterion, such as energy or resistance, while satisfying geometric constraints like fixed volume, surface area, or minimal height. Classically, solving these optimization problems requires numerical methods that rely on shape derivatives and adjoint calculations, which are computationally expensive and difficult to adapt to complex geometries. To overcome these challenges and the inherent limitations in parallelization found in classical approaches, our objective is to find a good methodology that relies on neural networks and benefits on their specific properties (easy parallelization, automatic differentiation, and a meshless approach). This requires suitable representations for both the solution of the state equation and the shape of the domain in which the equation is defined. For this purpose, we introduce the DeepRitz method to approximate solutions to the state equation and symplectic neural networks to model the domain shape effectively. We illustrate this approach through an example: minimizing the Dirichlet energy of a domain under a volume constraint. We will then propose an extention of this work for odd dimensions.
We present the study of the non-linear stability of a class of travelling-wave solutions to the compressible pressureless Navier-Stokes system with a singular viscosity. These solutions encode the effect of congestion by connecting a congested left state to an uncongested right state. By using carefully weighted energy estimates we are able to prove the non-linear stability of viscous shock waves to this system under a small zero integral perturbation, which in particular extends previous results that do not handle the case where the viscosity is singular. This is a joint work with Muhammed Ali Mehmood from Imperial College, London.
Cette présentation vise à introduire un schéma monolithique GD/VF (Galerkin Discontinu/Volumes Finis) préservant les propriétés convexes (principes du maximum, positivité, entropie, ...) pour la résolution des systèmes de lois de conservation sur maillages non-structurés. Il est bien connu que les méthodes Galerkin discontinue (GD) nécessitent une limitation non linéaire pour éviter les oscillations parasites ou les instabilités non linéaires, susceptibles de provoquer l´arrêt anticipé du code de simulation. L'idée principale de ce travail est d'améliorer la robustesse des schémas GD tout en préservant autant que possible leur grande précision et leur résolution de sous-maille. Pour ce faire, une combinaison convexe entre un schéma GD d'ordre élevé et un schéma volumes finis (VF) d'ordre un sera localement effectué, à l'échelle des sous-mailles, où cela sera nécessaire. À cette fin, nous prouverons tout d'abord qu'il est possible de réécrire un schéma GD comme un schéma VF défini sur un sous-maillage, en introduisant des flux numériques spécifiques, qu'on appellera flux reconstruits GD. Le schéma monolithique GD/VF sera alors défini de la manière suivante : à chaque face de chaque sous-cellule seront assignés deux flux, un flux VF d'ordre un et un flux reconstruit d'ordre élevé, qui seront finalement combinés de manière convexe. L'objectif est alors de déterminer, par analyse, les coefficients de combinaison optimaux pour atteindre les propriétés souhaitées (par exemple, positivité, absence d´oscillations, inégalités d´entropie) tout en préservant la grande précision du schéma. Des résultats numériques sur divers types de problèmes hyperboliques seront présentés pour évaluer les performances de la méthode proposée.
Grâce à ce formalisme monolithique, nous tenterons de répondre à certaines questions : est-il possible d'assurer une stabilité entropique ? de quelle stabilité entropique parlons-nous (discrète, semi-discrète, de maille, de sous-maille, pour quelle entropie, ...) ? quels sont les coûts de telles stabilités (en terme de précision ou de perte d'autres propriétés) ? À quel point est-ce essentiel pour les problèmes qui nous intéressent ? Nous présenterons différents résultats numériques pour tenter de partiellement répondre à ces questions
In this seminar we present a mathematical model to describe the evolution of a city, which is determined by the interaction of two large populations of agents, workers and firms. The map of the city is described by a network with the edges representing at the same time residential areas and communication routes. The two populations compete for space while interacting through the labour market. The resulting model is described by a two population Mean-Field Game system coupled with an Optimal Transport problem. We prove existence and uniqueness of the solution and we provide some numerical tools to develop several numerical simulations. This is a joint work with Fabio Camilli (Sapienza Roma) and Luciano Marzufero (Libera Università di Bolzano).
The goal of this talk is to give an overview of two new results in the development of hybrid methods for elliptic and hyperbolic partial differential equations (PDEs). A hybrid method combines classical numerical analysis techniques (finite element method (FEM), discontinuous Galerkin (DG), ...) with tools from machine learning (ML). The first part of this talk is dedicated to a broad presentation of such ML tools, including a common framework to represent PDE approximators, be they classical or ML-based. Then, in a second part, we explain how to use a physics-informed prior to lower the error constant of the FEM while keeping the same order of accuracy. Thanks to the FEM framework applied to elliptic PDEs, we rigorously prove that our correction improves the FEM error constant by a factor depending on the prior quality. If time permits, in a third part, we discuss how to enhance the DG basis with physics-informed priors, to increase the resolution of near-equilibrium solutions to hyperbolic systems of balance laws. Once again, we rigorously prove that the error constant is improved. Numerical illustrations will be present throughout the presentation, to validate our results.
Le modèle Euler compressible bifluide ne présente pas de difficultés théoriques supplémentaires comparé au cas monofluide. Mais sa résolution numérique est notoirement plus difficile à cause du phénomène d'oscillations de pression à l'interface entre fluides. Nous présentons une approche basée sur un échantillonnage aléatoire "à la Glimm" à l'interface, qui permet de s'affranchir de ce défaut. Le schéma obtenu est applicable à des maillages non structurés, il a d'excellentes propriétés de robustesse et de convergence. Nous l'appliquons à des cas de déferlement.
Pendant l'élongation de l’axe de l'embryon de vertébré, on observe, grâce à l’imagerie live, un phénomène de turbulence cellulaire dans les différents tissus embryonnaires. Nous proposons un modèle mécanique en 2D pour modéliser la croissance des tissus pendant l'élongation de l'embryon, qui permet de retrouver ces flux turbulents à travers un rotationnel non trivial pour les vitesses des tissus. Une autre spécificité de ce modèle est que la ségrégation entre les tissus est assurée par une pression de ségrégation. Après avoir déterminé (formellement) la limite incompressible, nous étudions le comportement qualitatif à la limite et discutons d'un effet fantôme.
Plusieurs modèles physiques permettent de comprendre la dynamique des mélanges de fluides, parmi lesquels figurent les modèles dits de Baer-Nunziato. Les équations aux dérivées partielles associées à ces modèles ressemblent à celles de Navier-Stokes, avec, en plus, de nouveaux termes de relaxation.
Une stratégie pour obtenir ces modèles est l'homogénéisation : à partir d'un mélange mésoscopique, où deux fluides purs satisfaisant les équations de Navier-Stokes compressibles se répartissent l'espace, on effectue un changement d'échelle pour obtenir un mélange macroscopique, où, en chaque point de l'espace, les deux fluides peuvent coexister.
Ce problème relève de l'étude des équations de Navier-Stokes avec des données initiales fortement oscillantes. On commencera donc par expliquer certains résultats dans ce cadre de travail, en dimension un d'espace et sur le tore, d'abord pour des fluides barotropes, puis pour des gaz parfaits non barotropes. On détaillera ensuite les différentes étapes de la démonstration de l'homogénéisation.
Dynamic processes in continuum physics are modeled using time-dependent partial differential equations (PDE), which are based on the conservation of some physical quantities, such as mass, momentum and energy. Depending on the physical phenomenon under consideration, the governing equations can exhibit some mathematical structures like differential constraints, algebraic relations, physical admissible states as well as asymptotic limits and thermodynamics compatibility. An interesting class of mathematical models is provided by symmetric hyperbolic systems that intrinsically imply all the structures listed above. When passing at the discrete level, the exact satisfaction of these structural properties is not automatically guaranteed, thus Structure Preserving numerical schemes have recently emerged with the aim of exactly discretizing at least a subset of these constraints. We will investigate and present some of our research activity carried out in the framework of the development of Structure Preserving schemes, focusing on recent contributions delivered in the last three years. In particular, we will address asymptotic preserving schemes for low Mach flows, div-curl and curl-grad preserving operators for discontinuous Galerkin methods, and a novel geometric and thermodynamically compatible finite volume method for continuum mechanics.
The energy of saline gradients is a very promising source of non-intermittent renewable energy, the exploitation of which is hampered by the lack of economically viable technology. The most investigated harvesting methods rely on selective transport of ions or water molecules through semi-permeable or ion-selective membranes, which demonstrate limited power densities of the order of a few W/m2. While in the last decade single nanofluidic objects such as nanopores of nanotubes have opened up very promising prospects with power density capabilities in the kW or even MW/m2, scale-up efforts face serious issues, as concentration polarization phenomena result in a massive loss of performance.
At the LiPhy we work on a concept of nanofluidic exchanger for power generation from saline gradients, focused on designing a nanoscale flow able to harvest the power at the output of the nanopores. We will present the study of a simple exchanger made of a selective nanoslit fed by a nanofluidic assembly. One specific feature of such an exchanger relies on the non-linear ion fluxes through the nanoslit, according to the so-called 1D Poisson-Nernt Planck equations. Such an elemental brick could be massively parallelized in stackable electricity-generating layers using standard technologies of the semi-conductors industry. We demonstrate here a scheme for rationalizing the choice of the exchanger parameters, taking into account the transport properties at all scales. The simplified numerical resolution of the three-dimensional device shows that net power densities of 300 W/m2 and more can be achieved.
Orateur(e)s :
Walter Boscheri (LAMA); Camille Carvalho, (ICJ); Frédérique Charles, (LJK); Nicolae Cindea, (LMBP); Sue Claret, (LMBP); Baptiste Devyver, (IF); Martin Donati, (IF); Louis Dupaigne, (ICJ); Hugo Eulry, (UMPA); Christophe Lacave, (LAMA); Mickael Nahon, (LJK); Niami Nasr, (ICJ); Pierre-Damien Thizy, (ICJ);
Orateur(e)s :
Walter Boscheri (LAMA); Camille Carvalho, (ICJ); Frédérique Charles, (LJK); Nicolae Cindea, (LMBP); Sue Claret, (LMBP); Baptiste Devyver, (IF); Martin Donati, (IF); Louis Dupaigne, (ICJ); Hugo Eulry, (UMPA); Christophe Lacave, (LAMA); Mickael Nahon, (LJK); Niami Nasr, (ICJ); Pierre-Damien Thizy, (ICJ);
Je présenterai quelques avancées récentes sur la caractérisation des phénomènes de propagation qui apparaissent dans les équations semi-linéaires avec diffusion non locale de type Levy. Récemment différentes dichotomies entre propagation accélérée et propagation à vitesse constante en fonction des paramètres de décroissance du noyau et de l'ordre d'annulation en zéro de la non-linéarité considérée ont été obtenues. Je me concentrerai sur le cas monostable et sur une manière de contourner les difficultés liées au traitement des opérateurs de Levy généraux.