m

F Nous contacter


0

CEMRACS  | enregistrements trouvés : 98

O

-A +A

P Q

Déposez votre fichier ici pour le déplacer vers cet enregistrement.

Recently, an important research activity on mean field games (MFGs for short) has been initiated by the pioneering works of Lasry and Lions: it aims at studying the asymptotic behavior of stochastic differential games (Nash equilibria) as the number $n$ of agents tends to infinity. The field is now rapidly growing in several directions, including stochastic optimal control, analysis of PDEs, calculus of variations, numerical analysis and computing, and the potential applications to economics and social sciences are numerous.
In the limit when $n \to +\infty$, a given agent feels the presence of the others through the statistical distribution of the states. Assuming that the perturbations of a single agent's strategy does not influence the statistical states distribution, the latter acts as a parameter in the control problem to be solved by each agent. When the dynamics of the agents are independent stochastic processes, MFGs naturally lead to a coupled system of two partial differential equations (PDEs for short), a forward Fokker-Planck equation and a backward Hamilton-Jacobi-Bellman equation.
The latter system of PDEs has closed form solutions in very few cases only. Therefore, numerical simulation are crucial in order to address applications. The present mini-course will be devoted to numerical methods that can be used to approximate the systems of PDEs.
The numerical schemes that will be presented rely basically on monotone approximations of the Hamiltonian and on a suitable weak formulation of the Fokker-Planck equation.
These schemes have several important features:

- The discrete problem has the same structure as the continous one, so existence, energy estimates, and possibly uniqueness can be obtained with the same kind of arguments

- Monotonicity guarantees the stability of the scheme: it is robust in the deterministic limit

- convergence to classical or weak solutions can be proved

Finally, there are particular cases named variational MFGS in which the system of PDEs can be seen as the optimality conditions of some optimal control problem driven by a PDE. In such cases, augmented Lagrangian methods can be used for solving the discrete nonlinear system. The mini-course will be orgamized as follows

1. Introduction to the system of PDEs and its interpretation. Uniqueness of classical solutions.

2. Monotone finite difference schemes

3. Examples of applications

4. Variational MFG and related algorithms for solving the discrete system of nonlinear equations
Recently, an important research activity on mean field games (MFGs for short) has been initiated by the pioneering works of Lasry and Lions: it aims at studying the asymptotic behavior of stochastic differential games (Nash equilibria) as the number $n$ of agents tends to infinity. The field is now rapidly growing in several directions, including stochastic optimal control, analysis of PDEs, calculus of variations, numerical analysis and ...

49K20 ; 49N70 ; 35K40 ; 35K55 ; 35Q84 ; 65K10 ; 65M06 ; 65M12 ; 91A23 ; 91A15

Déposez votre fichier ici pour le déplacer vers cet enregistrement.

Post-edited  Interview au CIRM : Yvon Maday
Maday, Yvon (Personne interviewée) | CIRM (Editeur )

Le CIRM : écrin estival du CEMRACS depuis 20 ans !

Déposez votre fichier ici pour le déplacer vers cet enregistrement.

At the end of the 70', Littlejohn [1, 2, 3] shed new light on what is called the Gyro-Kinetic Approximation. His approach incorporated high-level mathematical concepts from Hamiltonian Mechanics, Differential Geometry and Symplectic Geometry into a physical affordable theory in order to clarify what has been done for years in the domain. This theory has been being widely used to deduce the numerical methods for Tokamak and Stellarator simulation. Yet, it was formal from the mathematical point of view and not directly accessible for mathematicians.
This talk will present a mathematically rigorous version of the theory. The way to set out this Gyro-Kinetic Approximation consists of the building of a change of coordinates that decouples the Hamiltonian dynamical system satisfied by the characteristics of charged particles submitted to a strong magnetic field into a part that concerns the fast oscillation induced by the magnetic field and a other part that describes a slower dynamics.
This building is made of two steps. The goal of the first one, so-called "Darboux Algorithm", is to give to the Poisson Matrix (associated to the Hamiltonian system) a form that would achieve the goal of decoupling if the Hamiltonian function does not depend on one given variable. Then the second change of variables (which is in fact a succession of several ones), so-called "Lie Algorithm", is to remove the given variable from the Hamiltonian function without changing the form of the Poisson Matrix.
(Notice that, beside this Geometrical Gyro-Kinetic Approximation Theory, an alternative approach, based on Asymptotic Analysis and Homogenization Methods was developed in Frenod and Sonnendrücker [5, 6, 7], Frenod, Raviart and Sonnendrücker [4], Golse and Saint-Raymond [9] and Ghendrih, Hauray and Nouri [8].)
At the end of the 70', Littlejohn [1, 2, 3] shed new light on what is called the Gyro-Kinetic Approximation. His approach incorporated high-level mathematical concepts from Hamiltonian Mechanics, Differential Geometry and Symplectic Geometry into a physical affordable theory in order to clarify what has been done for years in the domain. This theory has been being widely used to deduce the numerical methods for Tokamak and Stellarator s...

70H05 ; 82D10 ; 58Z05 ; 58J37 ; 58J45 ; 58D10

Déposez votre fichier ici pour le déplacer vers cet enregistrement.

Post-edited  Darcy problem and crowd motion modeling
Maury, Bertrand (Auteur de la Conférence) | CIRM (Editeur )

We describe here formal analogies between the Darcy equations, that describe the flow of a viscous fluid in a porous medium, and some problems arising from the handing of congestion in crowd motion models.
At the microscopic level, individuals are identified to rigid discs, and the dual handling of the non overlapping constraint leads to discrete Darcy-like equations with a unilateral constraint that involves the velocities and interaction pressures, and that are set on the contact network. At the macroscopic level, a similar problem is obtained, that is set on the congested zone.
We emphasize the differences between the two settings: at the macroscopic level, a straight use of the maximum principle shows that congestion actually favors evacuation, which is in contradiction with experimental evidence. On the contrary, in the microscopic setting, the very particular structure of the discrete differential operators makes it possible to reproduce observed "Stop and Go waves", and the so called "Faster is Slower" effect.
We describe here formal analogies between the Darcy equations, that describe the flow of a viscous fluid in a porous medium, and some problems arising from the handing of congestion in crowd motion models.
At the microscopic level, individuals are identified to rigid discs, and the dual handling of the non overlapping constraint leads to discrete Darcy-like equations with a unilateral constraint that involves the velocities and interaction ...

34A60 ; 34D20 ; 35F31 ; 35R70 ; 70E50 ; 70E55

Déposez votre fichier ici pour le déplacer vers cet enregistrement.

Many physical phenomena deal with a fluid interacting with a moving rigid or deformable structure. These kinds of problems have a lot of important applications, for instance, in aeroelasticity, biomechanics, hydroelasticity, sedimentation, etc. From the analytical point of view as well as from the numerical point of view they have been studied extensively over the past years. We will mainly focus on viscous fluid interacting with an elastic structure. The purpose of the present lecture is to present an overview of some of the mathematical and numerical difficulties that may be encountered when dealing with fluid-structure interaction problems such as the geometrical nonlinearities or the added mass effect and how one can deal with these difficulties.
Many physical phenomena deal with a fluid interacting with a moving rigid or deformable structure. These kinds of problems have a lot of important applications, for instance, in aeroelasticity, biomechanics, hydroelasticity, sedimentation, etc. From the analytical point of view as well as from the numerical point of view they have been studied extensively over the past years. We will mainly focus on viscous fluid interacting with an elastic ...

74S05 ; 76M10 ; 74F10 ; 76D05

Déposez votre fichier ici pour le déplacer vers cet enregistrement.

Reduced MHD models in Tokamak geometry are convenient simplifications of full MHD and are fundamental for the numerical simulation of MHD stability in Tokamaks. This presentation will address the mathematical well-posedness and the justification of the such models.
The first result is a systematic design of hierachies of well-posed reduced MHD models. Here well-posed means that the system is endowed with a physically sound energy identity and that existence of a weak solution can be proved. Some of these models will be detailed.
The second result is perhaps more important for applications. It provides understanding on the fact the the growth rate of linear instabilities of the initial (non reduced) model is lower bounded by the growth rate of linear instabilities of the reduced model.
This work has been done with Rémy Sart.
Reduced MHD models in Tokamak geometry are convenient simplifications of full MHD and are fundamental for the numerical simulation of MHD stability in Tokamaks. This presentation will address the mathematical well-posedness and the justification of the such models.
The first result is a systematic design of hierachies of well-posed reduced MHD models. Here well-posed means that the system is endowed with a physically sound energy ...

76W05 ; 35L65 ; 65M60 ; 35Q30

Déposez votre fichier ici pour le déplacer vers cet enregistrement.

A simple, robust and accurate HLLC-type Riemann solver for two-phase 7-equation type models is built. It involves 4 waves per phase, i.e. the three conventional right- and left-facing and contact waves, augmented by an extra "interfacial" wave. Inspired by the Discrete Equations Method (Abgrall and Saurel, 2003), this wave speed $u_I$ is assumed function only of the piecewise constant initial data. Therefore it is computed easily from these initial data. The same is done for the interfacial pressure $P_I$. Interfacial variables $u_I$ and $P_I$ are thus local constants in the Riemann problem. Thanks to this property there is no difficulty to express the non-conservative system of partial differential equations in local conservative form. With the conventional HLLC wave speed estimates and the extra interfacial speed $u_I$, the four-waves Riemann problem for each phase is solved following the same strategy as in Toro et al. (1994) for the Euler equations. As $u_I$ and $P_I$ are functions only of the Riemann problem initial data, the two-phase Riemann problem consists in two independent Riemann problems with 4 waves only. Moreover, it is shown that these solvers are entropy producing. The method is easy to code and very robust. Its accuracy is validated against exact solutions as well as experimental data.
A simple, robust and accurate HLLC-type Riemann solver for two-phase 7-equation type models is built. It involves 4 waves per phase, i.e. the three conventional right- and left-facing and contact waves, augmented by an extra "interfacial" wave. Inspired by the Discrete Equations Method (Abgrall and Saurel, 2003), this wave speed $u_I$ is assumed function only of the piecewise constant initial data. Therefore it is computed easily from these ...

76Mxx ; 76TXX

Déposez votre fichier ici pour le déplacer vers cet enregistrement.
Déposez votre fichier ici pour le déplacer vers cet enregistrement.

The momentum transport in a fusion device such as a tokamak has been in a scope of the interest during last decade. Indeed, it is tightly related to the plasma rotation and therefore its stabilization, which in its turn is essential for the confinement improvement. The intrinsic rotation, i.e. the part of the rotation occurring without any external torque is one of the possible sources of plasma stabilization.
The modern gyrokinetic theory [3] is an ubiquitous theoretical framework for lowfrequency fusion plasma description. In this work we are using the field theory formulation of the modern gyrokinetics [1]. The main attention is focussed on derivation of the momentum conservation law via the Noether method, which allows to connect symmetries of the system with conserved quantities by means of the infinitesimal space-time translations and rotations.
Such an approach allows to consistently keep the gyrokinetic dynamical reduction effects into account and therefore leads towards a complete momentum transport equation.
Elucidating the role of the gyrokinetic polarization is one of the main results of this work. We show that the terms resulting from each step of the dynamical reduction (guiding-center and gyrocenter) should be consistently taken into account in order to establish physical meaning of the transported quantity. The present work [2] generalizes previous result obtained in [4] by taking into the account purely geometrical contributions into the radial polarization.
The momentum transport in a fusion device such as a tokamak has been in a scope of the interest during last decade. Indeed, it is tightly related to the plasma rotation and therefore its stabilization, which in its turn is essential for the confinement improvement. The intrinsic rotation, i.e. the part of the rotation occurring without any external torque is one of the possible sources of plasma stabilization.
The modern gyrokinetic theory [3] ...

82D10 ; 82C40 ; 35L65 ; 35Q83 ; 70S10

Déposez votre fichier ici pour le déplacer vers cet enregistrement.

This lecture will present a short overview on kinetic MHD. The advantages and drawbacks of kinetic versus fluid modelling will be summarized. Various techniques to implement kinetic effects in the fluid description will be introduced with increasing complexity: bi-fluid effects, gyroaverage fields, Landau closures. Hybrid formulations, which combine fluid and kinetic approaches will be presented. It will be shown that these formulations raise several difficulties, including inconsistent ordering and choice of representation. The non linear dynamics of an internal kink mode in a tokamak will be used as a test bed for the various formulations. It will be shown that bi-fluid effects can explain to some extent fast plasma relaxations (reconnection), but cannot address kinetic instabilities due to energetic particles. Some results of hybrid codes will be shown. Recent developments and perspectives will be given in conclusion.
This lecture will present a short overview on kinetic MHD. The advantages and drawbacks of kinetic versus fluid modelling will be summarized. Various techniques to implement kinetic effects in the fluid description will be introduced with increasing complexity: bi-fluid effects, gyroaverage fields, Landau closures. Hybrid formulations, which combine fluid and kinetic approaches will be presented. It will be shown that these formulations raise ...

82D10 ; 76W05

Déposez votre fichier ici pour le déplacer vers cet enregistrement.
Déposez votre fichier ici pour le déplacer vers cet enregistrement.
Déposez votre fichier ici pour le déplacer vers cet enregistrement.
Déposez votre fichier ici pour le déplacer vers cet enregistrement.
Déposez votre fichier ici pour le déplacer vers cet enregistrement.

We analyse how reverting Random Number Generator can be efficiently used to save memory in solving dynamic programming equation. For SDEs, it takes the form of forward and backward Euler scheme. Surprisingly the error induced by time reversion is of order 1.

60H10 ; 60H15 ; 60H30 ; 65C10

Déposez votre fichier ici pour le déplacer vers cet enregistrement.

The mathematical framework of variational inequalities is a powerful tool to model problems arising in mechanics such as elasto-plasticity where the physical laws change when some state variables reach a certain threshold [1]. Somehow, it is not surprising that the models used in the literature for the hysteresis effect of non-linear elasto-plastic oscillators submitted to random vibrations [2] are equivalent to (finite dimensional) stochastic variational inequalities (SVIs) [3]. This presentation concerns (a) cycle properties of a SVI modeling an elasto-perfectly-plastic oscillator excited by a white noise together with an application to the risk of failure [4,5]. (b) a set of Backward Kolmogorov equations for computing means, moments and correlation [6]. (c) free boundary value problems and HJB equations for the control of SVIs. For engineering applications, it is related to the problem of critical excitation [7]. This point concerns what we are doing during the CEMRACS research project. (d) (if time permits) on-going research on the modeling of a moving plate on turbulent convection [8]. This is a mixture of joint works and / or discussions with, amongst others, A. Bensoussan, L. Borsoi, C. Feau, M. Huang, M. Laurière, G. Stadler, J. Wylie, J. Zhang and J.Q. Zhong.
The mathematical framework of variational inequalities is a powerful tool to model problems arising in mechanics such as elasto-plasticity where the physical laws change when some state variables reach a certain threshold [1]. Somehow, it is not surprising that the models used in the literature for the hysteresis effect of non-linear elasto-plastic oscillators submitted to random vibrations [2] are equivalent to (finite dimensional) stochastic ...

74H50 ; 35R60 ; 60H10 ; 60H30 ; 74C05

Déposez votre fichier ici pour le déplacer vers cet enregistrement.

The theory of mean field type control (or control of MacKean-Vlasov) aims at describing the behaviour of a large number of agents using a common feedback control and interacting through some mean field term. The solution to this type of control problem can be seen as a collaborative optimum. We will present the system of partial differential equations (PDE) arising in this setting: a forward Fokker-Planck equation and a backward Hamilton-Jacobi-Bellman equation. They describe respectively the evolution of the distribution of the agents' states and the evolution of the value function. Since it comes from a control problem, this PDE system differs in general from the one arising in mean field games.
Recently, this kind of model has been applied to crowd dynamics. More precisely, in this talk we will be interested in modeling congestion effects: the agents move but try to avoid very crowded regions. One way to take into account such effects is to let the cost of displacement increase in the regions where the density of agents is large. The cost may depend on the density in a non-local or in a local way. We will present one class of models for each case and study the associated PDE systems. The first one has classical solutions whereas the second one has weak solutions. Numerical results based on the Newton algorithm and the Augmented Lagrangian method will be presented.
This is joint work with Yves Achdou.
The theory of mean field type control (or control of MacKean-Vlasov) aims at describing the behaviour of a large number of agents using a common feedback control and interacting through some mean field term. The solution to this type of control problem can be seen as a collaborative optimum. We will present the system of partial differential equations (PDE) arising in this setting: a forward Fokker-Planck equation and a backward Hamilto...

35K40 ; 35K55 ; 35K65 ; 35D30 ; 49N70 ; 49K20 ; 65K10 ; 65M06

Déposez votre fichier ici pour le déplacer vers cet enregistrement.

In this work, we consider the discretization of some nonlinear Fokker-Planck-Kolmogorov equations. The scheme we propose preserves the non-negativity of the solution, conserves the mass and, as the discretization parameters tend to zero, has limit measure-valued trajectories which are shown to solve the equation. This convergence result is proved by assuming only that the coefficients are continuous and satisfy a suitable linear growth property with respect to the space variable. In particular, under these assumptions, we obtain a new proof of existence of solutions for such equations.
We apply our results to several examples, including Mean Field Games systems and variations of the Hughes model for pedestrian dynamics.
In this work, we consider the discretization of some nonlinear Fokker-Planck-Kolmogorov equations. The scheme we propose preserves the non-negativity of the solution, conserves the mass and, as the discretization parameters tend to zero, has limit measure-valued trajectories which are shown to solve the equation. This convergence result is proved by assuming only that the coefficients are continuous and satisfy a suitable linear growth property ...

35K55 ; 35Q84 ; 60H15 ; 60H30

Déposez votre fichier ici pour le déplacer vers cet enregistrement.

Multi angle  Project evaluation under uncertainty
Zubelli, Jorge P. (Auteur de la Conférence) | CIRM (Editeur )

Industrial strategic decisions have evolved tremendously in the last decades towards a higher degree of quantitative analysis. Such decisions require taking into account a large number of uncertain variables and volatile scenarios, much like financial market investments. Furthermore, they can be evaluated by comparing to portfolios of investments in financial assets such as in stocks, derivatives and commodity futures. This revolution led to the development of a new field of managerial science known as Real Options.
The use of Real Option techniques incorporates also the value of flexibility and gives a broader view of many business decisions that brings in techniques from quantitative finance and risk management. Such techniques are now part of the decision making process of many corporations and require a substantial amount of mathematical background. Yet, there has been substantial debate concerning the use of risk neutral pricing and hedging arguments to the context of project evaluation. We discuss some alternatives to risk neutral pricing that could be suitable to evaluation of projects in a realistic context with special attention to projects dependent on commodities and non-hedgeable uncertainties. More precisely, we make use of a variant of the hedged Monte-Carlo method of Potters, Bouchaud and Sestovic to tackle strategic decisions. Furthermore, we extend this to different investor risk profiles. This is joint work with Edgardo Brigatti, Felipe Macias, and Max O. de Souza.
If time allows we shall also discuss the situation when the historical data for the project evaluation is very limited and we can make use of certain symmetries of the problem to perform (with good estimates) a nonintrusive stratified resampling of the data. This is joint work with E. Gobet and G. Liu.
Industrial strategic decisions have evolved tremendously in the last decades towards a higher degree of quantitative analysis. Such decisions require taking into account a large number of uncertain variables and volatile scenarios, much like financial market investments. Furthermore, they can be evaluated by comparing to portfolios of investments in financial assets such as in stocks, derivatives and commodity futures. This revolution led to the ...

91B26 ; 91B06 ; 91B30 ; 91B24

Déposez votre fichier ici pour le déplacer vers cet enregistrement.

The valuation of American options (a widespread type of financial contract) requires the numerical solution of an optimal stopping problem. Numerical methods for such problems have been widely investigated. Monte-Carlo methods are based on the implementation of dynamic programming principles coupled with regression techniques. In lower dimension, one can choose to tackle the related free boundary PDE with deterministic schemes.
Pricing of American options will therefore be inevitably heavier than the one of European options, which only requires the computation of a (linear) expectation. The calibration (fitting) of a stochastic model to market quotes for American options is therefore an a priori demanding task. Yet, often this cannot be avoided: on exchange markets one is typically provided only with market quotes for American options on single stocks (as opposed to large stock indexes - e.g. S&P500 - for which large amounts of liquid European options are typically available).
In this talk, we show how one can derive (approximate, but accurate enough) explicit formulas - therefore replacing other numerical methods, at least in a low-dimensional case - based on asymptotic calculus for diffusions.
More precisely: based on a suitable representation of the PDE free boundary, we derive an approximation of this boundary close to final time that refines the expansions known so far in the literature. Via the early premium formula, this allows to derive semi-closed expressions for the price of the American put/call. The final product is a calibration recipe of a Dupire's local volatility to American option data.
Based on joint work with Pierre Henry-Labordère.
The valuation of American options (a widespread type of financial contract) requires the numerical solution of an optimal stopping problem. Numerical methods for such problems have been widely investigated. Monte-Carlo methods are based on the implementation of dynamic programming principles coupled with regression techniques. In lower dimension, one can choose to tackle the related free boundary PDE with deterministic schemes.
Pricing of ...

93E20 ; 91G60

Z