Adaptive evolution : a population approach Benoît Perthame
|
|
- Roland May
- 5 years ago
- Views:
Transcription
1 Adaptive evolution : a population approach Benoît Perthame
2 Adaptive dynamic : selection principle d dt n(x, t) = n(x, t)r( x, ϱ(t) ), ϱ(t) = R d n(x, t)dx. given x, n(x) = ϱ δ(x x), R( x, ϱ) = 0, ϱ( x). They are stable by perturbation of the weight ϱ (strong topology) d ϱ(t) = ϱ(t)r( x, ϱ(t)). dt But they are unstable by approximation in measures (weak topology), and by mutation (structural)...
3 Adaptive dynamic : mutations Off-springs undergo small mutations that change slightly the trait t n(x, t) n = n(x, t)r( x, ϱ(t) ), ϱ(t) = R d n(x, t)dx. t n(x, t) = n(x, t)r( x, ϱ(t) ) + ϱ(t) = R d n(x, t)dx. M(x, y)b(y)n(y, t)dy,
4 Adaptive dynamic : mutations Off-springs undergo small mutations that change slightly the trait t n(x, t) n = n(x, t)r( x, ϱ(t) ), ϱ(t) = R d n(x, t)dx. We assume that mutations are RARE and introduce a scale ε for small mutations ε t n(x, t) ε2 n = n(x, t)r ( x, ϱ(t) ), ϱ(t) = R d n(x, t)dx.
5 Adaptive dynamic : mutations Off-springs undergo small mutations that change slightly the trait t n(x, t) n = n(x, t)r( x, ϱ(t) ), ϱ(t) = R d n(x, t)dx. We assume that mutations are RARE and introduce a scale ε for small mutations ε t n ε(x, t) ε 2 n ε = n ε (x, t)r ( x, ϱ ε (t) ), ϱ ε (t) = R d n ε (x, t)dx.
6 Population model of adaptive dynamics : mutations This is not far from Fisher/KPP equation for invasion fronts/chemical reaction ε t n ε(x, t) ε 2 n ε = n ε (x, t) ( 1 n ε (x, t) ),. {n ε = 1} {n ε = 0} WKB, large deviations, level sets, geometric motion G. Barles, L. C. Evans, W. Fleming, P. E. Souganidis, S. Osher, J. Sethian...
7 Population model of adaptive dynamics : mutations This is not far from Fisher/KPP equation for invasion fronts/chemical reaction ε t n ε(x, t) ε 2 n ε = n ε (x, t) ( 1 n ε (x, t) ),. {n ε = 1} {n ε = 0} in the limit n(x, t)(1 n(x, t) = 0.
8 Population model of adaptive dynamics : mutations The situation is very different for the nonlocal equation ε t n ε(x, t) ε 2 n ε = n ε (x, t)r ( x, ϱ ε (t) ), ϱ ε (t) = R d n ε(x, t)dx.
9 Population model of adaptive dynamics : mutations The situation is very different for the nonlocal equation ε t n ε(x, t) ε 2 n ε = n ε (x, t)r ( x, ϱ ε (t) ), ϱ ε (t) = n ε(x, t)dx. R d In the limit one can expect 0 = n(x, t)r ( x, ϱ(t) ), n(x, t) = ϱδ Γ(t), Γ(t) { R(, ϱ(t)) = 0 }.
10 Asymptotic method Question. What tools to describe Dirac concentrations in PDEs? n ε (x) = ϱ (2πε) d/2e x x 2 /(2ε) ε 0 ϱδ(x x) n ε (x) = e ( x x 2 +ε ln O(ε))/(2ε) ε 0 ϱδ(x x) More generally (Hopf-Cole/WKB) with the conditions n ε (x) = e ϕ ε(x)/ε ε 0 ϱδ(x x) ϕ ε ε 0 ϕ, max x ϕ(x) = 0 = ϕ( x)
11 Asymptotic method Theorem Suppose x R, R x > 0, R ϱ < 0. Then, for subsequences n ε (x, t) ϱ(t)δ(x = x(t)), ϱ ε ϱ(t) = n(x, t)dx, εk 0 εk 0 Can one give a law for the dynamics of x(t)?
12 Asymptotic method Theorem Suppose x R, R x > 0, R ϱ < 0. Then, for subsequences n ε (x, t) ϱ(t)δ(x = x(t)), ϱ ε ϱ(t) = n(x, t)dx, εk 0 εk 0 and the fittest trait x(t) is characterised by the Eikonal equation with constraints t ϕ(x, t) = R( x, ϱ(t) ) + ϕ(x, t) 2 max x ϕ(x, t) = 0 = ϕ(t, x(t)) Definition This situation is called monomorphism Difficulty Solutions to H.-J. eq. are not smooth
13 Asymptotic method Theorem Suppose x R, R x > 0, R ϱ < 0. Then, for subsequences n ε (x, t) ϱ(t)δ(x = x(t)), ϱ ε ϱ(t) = n(x, t)dx, εk 0 εk 0 and the fittest trait x(t) is characterised by the Eikonal equation with constraints However t ϕ(x, t) = R( x, ϱ(t) ) + ϕ(x, t) 2 max x ϕ(x, t) = 0 = ϕ(t, x(t)) t ϕ( x(t), t ) = 0 = x ϕ( x(t), t ) R ( x(t), ϱ(t) ) = 0 (Pessimism Principle)
14 Asymptotic method This problem should be understood as follows max x ϱ(t) ϕ(x, t) = 0, t is a constraint, is a Lagrange multiplier.. This is not an obstacle problem! 0!!"! -0.05!!"# -0.1!!"$ -0.15!!"% -0.2!!"& -0.25!!"' -0.3!!"( !!")!"!!"#!"$!"%!"&!"'!"(!")!"*!"+ #"!
15 Asymptotic method Theorem In R d, set n ε (x, t) = e ϕ ε(x,t)/ε. After extraction, ϕ ε εk 0 ϕ (locally uniformly), ϱ ε(t) εk 0 ϱ(t) And n ε (x, t) εk 0 t ϕ(x, t) = R( x, ϱ(t) ) + ϕ(x, t) ) 2 max x ϕ(x, t) = 0 ( = ϕ(t, x(t)) n(x, t) weakly in measures, supp ( n(t) ) {ϕ(t) = 0} ).
16 Asymptotic method Proof 1. ϱ ε (t) is BV (and converges after extraction) and its limit ϱ(t) is non-decreasing 2. Because n ε (x, t) = e ϕ ε(x,t)/ε we have t ϕ ε(x, t) = R ( x, ϱ ε (t) ) + ϕ ε (x, t) 2 ε ϕ ε and ϕ ε is Lipschitz continous in x (difficulty in t) (gives the H.-J. equation in viscosity sense) 3. ϱ m n ε (x, t)dx ϱ M (gives the constraint)
17 Asymptotic method Conclusion : n ε L ϕ ε C
18 Asymptotic method Theorem (G. Barles, BP) Uniqueness With reasonable assumptions there exist a unique lipschitz continuous solution ( ϱ, ϕ) to the constraint H.-J. equation t ϕ(x, t) = b(x) ϱ(t)d(x) + ϕ 2, max x ϕ(x, t) = 0 ( = ϕ(t, x(t)) ) Open question Extend uniqueness to ϕ(x, t) = b(x) t 1 + ϱ(t) ϱ(t)d(x) + ϕ 2.
19 Asymptotic method Proof of uniqueness : the difficulty The L contraction property is lost! Define M(t) := max x [ϕ 1(x, t) ϕ 2 (x, t)] d dt M(t) R(x M(t), ϱ 1 (t)) R(x M (t), ϱ 2 (t)) C ϱ 1 (t) ϱ 2 (t) But the constraint cannot be used here to control ϱ 1 (t) ϱ 2 (t) by M(t).
20 Asymptotic method Proof of uniqueness : idea. Work on Define ψ(t) := ϕ(x, t) + d(x) R(x, ϱ) = b(x) d(x)ϱ t 0 ϱ(s)ds. t ψ(x, t) = b(x) + ψ d(x) t 0 ϱ(s)ds 2, max x ϕ(x, t) = 0 ( = ϕ(t, x(t)) ) M(t) := max x [ψ 1(x, t) ψ 2 (x, t)] d dt M(t) p M d(x M ) t 0 ϱ 1(s)ds 2 p M d(x M ) t 0 ϱ 2(s)ds 2
21 Asymptotic method t d dt M(t) p M d(x M ) ϱ 1(s)ds 2 p M d(x M ) ϱ 2(s)ds Use that solutions are Lipschitz and the specific form of R t d dt M(t) C ϱ 1(s)ds 0 But we may choose ϕ(t, x 1 ) = 0 and get t 0 ϱ 2(s)ds M(t) ψ 1 (t, x 1 ) ψ 2 (t, x 1 ) d(x 1 ) [ t ϱ 1(s)ds 0 The opposite inequality holds true similarly and thus d dt M(t) CM(t). t t 0 ϱ 2(s)ds ]
22 Asymptotic method Numerical tests : b(x) =.5 + x Direct simulation (1500 points) H.-J. solution (200 points)
23 Asymptotic method Numerical tests : b(x) =.5 + x(2 x) Direct simulation (1500 points) H.-J. solution (200 points)
24 Asymptotic method Numerical tests : min(.45 + x. 2, x) Direct simulation (1500 points) H.-J. solution (200 points)
25 Survival threshold Next ingredient is the notion of survival threshold. n(t, x) t ε n(t, x) = n(t, x) R ( x, [n(t)] ) ε nn(t, x) ε Motivated by Population really vanishes ; some traits are not represented The notion of individual is somehow included in the parameter n A because n(t, x) really vanishes at a level related to n compatibility with Monte-Carlo simulations
26 Survival threshold Next ingredient is the notion of survival threshold. n(t, x) t ε n(t, x) = n(t, x) R ( x, [n(t)] ) ε nn(t, x) ε Motivated by Population really vanishes ; some traits are not represented The notion of individual is somehow included in the parameter n A similar notion represents demographic stochasticity compatibility with Monte-Carlo simulations
27 Survival threshold n(t, x) n(t, x) ε n(t, x) = R ( x, [n(t)] ) t ε and choose a threshold of the form n n(t, x) ε n = e ϕ st/ε, ϕ st < 0 (constants). The formal constrained H.-J. equation is a free boundary problem ϕ t = ϕ 2 + R in Ω(t) := {(x, t), s.t. ϕ > ϕ st }, ϕ = in Ω c, ϕ ϕ st in Ω. Open questions : prove it rigorously ; other scales
28 n(t, x) ε n(t, x) = t n(t, x) R ( x, [n(t)] ) ε n n(t, x) ε n = e ϕ st/ε, ϕ st < 0 (constants), Formal derivation n ε = e ϕ ε/ε, then ϕ ε t ε ϕ ε ϕ ε 2 = R ( x, [n(t)] ) n n(t, x) ϕ ε t ε ϕ ε ϕ ε 2 = R ( x, [n(t)] ) e ϕst ϕε 2ε when ϕ st < ϕ ε then e ϕ st ϕ ε 2ε when ϕ st > ϕ ε then e ϕ st ϕ ε 2ε 0 (disappears) i.e. ϕ ε
29 Survival threshold Theorem Fix R(x) 0 then the free boundary problem ϕ ε ε 0 ϕ(t, x) ϕ t = ϕ 2 + R(x) in Ω(t) := {(x, t), s.t. ϕ > ϕ st }, ϕ = in Ω c, ϕ ϕ st in Ω. characterized by one of the equivalent statements it is the minimal solution the Dirichlet boundary condition should be satisfied ϕ = ϕ st on Ω(t). ϕ = is a truncation to of the global solution in R d.
30 Survival threshold When R(x) changes sign. The previous truncation formula is wrong The additional Dirichlet boundary condition is not enough, one should maybe impose also a state constraint boundary condition in Ω(t) The semi-relaxed limits can be compared to two relatively close functions of optimal control type (given by an iterative cleaning ). This implies at least that In opposition to the case R 0, the solution is changed drastically (see numerics) The limit does not depend on the specific power n
31 Survival threshold n(t, x) t ε n(t, x) = n(t, x) R(x) ε nn(t, x) ε Related to another asymptotic (other scales) : Bernouilli problem (see Lorz, Markowich, BP) n(x) + n(x) = R(x) 0, x Ω n(x) = 0 x Ω, n ν = n x Ω.
32 Numerical results Effect of the survival threshold
33 Numerical results trait x
34 Model with two nutrients : no survival threshold ϱ 1 (t) = ψ 1 (x)n(x, t)dx, ϱ 2 (t) = See Champagnat and Jabin ψ 2 (x)n(x, t)dx Direct simulation (1500 points) H.-J. solution (200 points) density n phase ϕ
35 Numerical results And the dynamics looks like
36 Open questions Uniqueness for a general R(x, ϱ) Case of multiple nutrients (See Champagnat and Jabin) R := R(x, ϱ 1, ϱ 2,..., ϱ I ), ϱ i = Survival threshold (R(x, ϱ), other scales) Explain branching ψ i (x)n(x, t)dx
Adaptive evolution : a population approach Benoît Perthame
Adaptive evolution : a population approach Benoît Perthame 30 20 t 50 t 10 x 0 1 0.5 0 0.5 1 x 0 1 0.5 0 0.5 1 Motivation : population adaptive evolution Darwin (1809-1882) On the origin of species (1859)
More informationAdaptive evolution : A population point of view Benoît Perthame
Adaptive evolution : A population point of view Benoît Perthame 30 20 t 50 t 10 x 0 1 0.5 0 0.5 1 x 0 1 0.5 0 0.5 1 Population formalism n(x, t) t }{{} variation of individuals = birth with mutations {}}{
More informationAdaptive evolution : a population approach Benoît Perthame
Adaptive evolution : a population approach Benoît Perthame 30 20 t 50 t 10 x 0 1 0.5 0 0.5 1 x 0 1 0.5 0 0.5 1 OUTLINE OF THE LECTURE DIRECT COMPETITION AND POLYMORPHIC CONCENTRATIONS I. Direct competition
More informationAdaptive evolution : a population approach Benoît Perthame
Adaptive evolution : a population approach Benoît Perthame OUTLINE OF THE LECTURE Interaction between a physiological trait and space I. Space and physiological trait II. Selection of dispersal (bounded
More informationDeterministic and stochastic approaches for evolutionary branching in adaptive dynamics
Deterministic and stochastic approaches for evolutionary branching in adaptive dynamics Nicolas Champagnat (Inria Nancy, Institut Elie Cartan de Lorraine (IECL)) 57th Annual AustMS Meeting, Sydney, 2nd
More informationSingular limits for reaction-diffusion equations with fractional Laplacian and local or nonlocal nonlinearity
Singular limits for reaction-diffusion equations with fractional Laplacian and local or nonlocal nonlinearity Sylvie Méléard, Sepideh Mirrahimi September 2, 214 Abstract We perform an asymptotic analysis
More informationAdaptive dynamics in an individual-based, multi-resource chemostat model
Adaptive dynamics in an individual-based, multi-resource chemostat model Nicolas Champagnat (INRIA Nancy) Pierre-Emmanuel Jabin (Univ. Maryland) Sylvie Méléard (Ecole Polytechnique) 6th European Congress
More informationA PDE approach to studying evolutionary and spatial dynamics in cancer cell populations
A PDE approach to studying evolutionary and spatial dynamics in cancer cell populations Tommaso Lorenzi LJLL, 30th June 2017 T. Lorenzi (University of St Andrews) June 30th 2 / 43 Synopsis Object of study
More informationExistence of viscosity solutions for a nonlocal equation modelling polymer
Existence of viscosity solutions for a nonlocal modelling Institut de Recherche Mathématique de Rennes CNRS UMR 6625 INSA de Rennes, France Joint work with P. Cardaliaguet (Univ. Paris-Dauphine) and A.
More informationNew Physical Principle for Monte-Carlo simulations
EJTP 6, No. 21 (2009) 9 20 Electronic Journal of Theoretical Physics New Physical Principle for Monte-Carlo simulations Michail Zak Jet Propulsion Laboratory California Institute of Technology, Advance
More informationTransport Equations in Biology; Adaptive Dynamics
Transport Equations in Biology; Adaptive Dynamics Matt Becker University of Maryland, College Park Applied Mathematics, Statistics, and Scientific Computation March 26, 2014 Becker Cancer RIT University
More informationConcentration in Lotka-Volterra parabolic or integral equations: a general convergence result
Concentration in Lotka-Volterra parabolic or integral equations: a general convergence result Guy Barles, Sepideh Mirrahimi, Benoît Perthame To cite this version: Guy Barles, Sepideh Mirrahimi, Benoît
More informationEvolution of species trait through resource competition
Evolution of species trait through resource competition Sepideh Mirrahimi, Benoît Perthame, Joe Yuichiro Wakano To cite this version: Sepideh Mirrahimi, Benoît Perthame, Joe Yuichiro Wakano. Evolution
More informationCompactness in Ginzburg-Landau energy by kinetic averaging
Compactness in Ginzburg-Landau energy by kinetic averaging Pierre-Emmanuel Jabin École Normale Supérieure de Paris AND Benoît Perthame École Normale Supérieure de Paris Abstract We consider a Ginzburg-Landau
More informationREPEATED GAMES FOR NON-LINEAR PARABOLIC INTEGRO-DIFFERENTIAL EQUATIONS AND INTEGRAL CURVATURE FLOWS
REPEATED GAMES FOR NON-LINEAR PARABOLIC INTEGRO-DIFFERENTIAL EQUATIONS AND INTEGRAL CURVATURE FLOWS CYRIL IMBERT AND SYLVIA SERFATY Abstract. The main purpose of this paper is to approximate several non-local
More informationFrequency functions, monotonicity formulas, and the thin obstacle problem
Frequency functions, monotonicity formulas, and the thin obstacle problem IMA - University of Minnesota March 4, 2013 Thank you for the invitation! In this talk we will present an overview of the parabolic
More informationApplied Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books.
Applied Analysis APPM 44: Final exam 1:3pm 4:pm, Dec. 14, 29. Closed books. Problem 1: 2p Set I = [, 1]. Prove that there is a continuous function u on I such that 1 ux 1 x sin ut 2 dt = cosx, x I. Define
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationFiltered scheme and error estimate for first order Hamilton-Jacobi equations
and error estimate for first order Hamilton-Jacobi equations Olivier Bokanowski 1 Maurizio Falcone 2 2 1 Laboratoire Jacques-Louis Lions, Université Paris-Diderot (Paris 7) 2 SAPIENZA - Università di Roma
More informationEvolution in prey-predator systems
Darwin: Le hasard et l évolution S. Méléard, Ecole Polytechnique, Centre de Mathématiques CMAP, École Polytechnique Appliquées in collaboration with N. Loeuille, C. Hauzy and S. Méléard 02/04/2013, Aussois
More informationPropagation d interfaces avec termes non locaux
Propagation d interfaces avec termes non locaux P. Cardaliaguet Univ. Brest Janvier 2008 Joint works with G. Barles (Tours), O. Alvarez (Rouen), O. Ley (Tours), R. Monneau (CERMICS), A. Monteillet (Brest).
More informationKramers formula for chemical reactions in the context of Wasserstein gradient flows. Michael Herrmann. Mathematical Institute, University of Oxford
eport no. OxPDE-/8 Kramers formula for chemical reactions in the context of Wasserstein gradient flows by Michael Herrmann Mathematical Institute, University of Oxford & Barbara Niethammer Mathematical
More informationModelling of interfaces and free boundaries
University of Regensburg Regensburg, March 2009 Outline 1 Introduction 2 Obstacle problems 3 Stefan problem 4 Shape optimization Introduction What is a free boundary problem? Solve a partial differential
More informationHomogenization of a Hele-Shaw-type problem in periodic time-dependent med
Homogenization of a Hele-Shaw-type problem in periodic time-dependent media University of Tokyo npozar@ms.u-tokyo.ac.jp KIAS, Seoul, November 30, 2012 Hele-Shaw problem Model of the pressure-driven }{{}
More informationEntropy and Relative Entropy
Entropy and Relative Entropy Joshua Ballew University of Maryland October 24, 2012 Outline Hyperbolic PDEs Entropy/Entropy Flux Pairs Relative Entropy Weak-Strong Uniqueness Weak-Strong Uniqueness for
More informationAPPLICATIONS OF DIFFERENTIABILITY IN R n.
APPLICATIONS OF DIFFERENTIABILITY IN R n. MATANIA BEN-ARTZI April 2015 Functions here are defined on a subset T R n and take values in R m, where m can be smaller, equal or greater than n. The (open) ball
More informationMotivation Power curvature flow Large exponent limit Analogues & applications. Qing Liu. Fukuoka University. Joint work with Prof.
On Large Exponent Behavior of Power Curvature Flow Arising in Image Processing Qing Liu Fukuoka University Joint work with Prof. Naoki Yamada Mathematics and Phenomena in Miyazaki 2017 University of Miyazaki
More informationSobolev spaces. May 18
Sobolev spaces May 18 2015 1 Weak derivatives The purpose of these notes is to give a very basic introduction to Sobolev spaces. More extensive treatments can e.g. be found in the classical references
More informationInput to state Stability
Input to state Stability Mini course, Universität Stuttgart, November 2004 Lars Grüne, Mathematisches Institut, Universität Bayreuth Part III: Lyapunov functions and quantitative aspects ISS Consider with
More informationParameter Dependent Quasi-Linear Parabolic Equations
CADERNOS DE MATEMÁTICA 4, 39 33 October (23) ARTIGO NÚMERO SMA#79 Parameter Dependent Quasi-Linear Parabolic Equations Cláudia Buttarello Gentile Departamento de Matemática, Universidade Federal de São
More informationA Hamilton-Jacobi approach for a model of population structured by space and trait
A Hamilton-Jacobi approach for a model of population structured by space and trait Emeric Bouin Sepideh Mirrahimi March 26, 2014 Abstract We study a non-local parabolic Lotka-Volterra type equation describing
More informationFlux-limited solutions for quasi-convex Hamilton-Jacobi equations on networks
Flux-limited solutions for quasi-convex Hamilton-Jacobi equations on networks C. Imbert and R. Monneau June 24, 2014 Abstract We study Hamilton-Jacobi equations on networks in the case where Hamiltonians
More informationLucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche
Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second
More informationREMARKS ON THE VANISHING OBSTACLE LIMIT FOR A 3D VISCOUS INCOMPRESSIBLE FLUID
REMARKS ON THE VANISHING OBSTACLE LIMIT FOR A 3D VISCOUS INCOMPRESSIBLE FLUID DRAGOŞ IFTIMIE AND JAMES P. KELLIHER Abstract. In [Math. Ann. 336 (2006), 449-489] the authors consider the two dimensional
More informationSébastien Chaumont a a Institut Élie Cartan, Université Henri Poincaré Nancy I, B. P. 239, Vandoeuvre-lès-Nancy Cedex, France. 1.
A strong comparison result for viscosity solutions to Hamilton-Jacobi-Bellman equations with Dirichlet condition on a non-smooth boundary and application to parabolic problems Sébastien Chaumont a a Institut
More informationAnalysis of PDE models of molecular motors. Benoît Perthame based on works with P. E. Souganidis
Analysis of PDE models of molecular motors Benoît Perthame based on works with P. E. Souganidis Asymmetry and motion Except human constructions movement always uses asymmetry Asymmetry and motion See M.
More informationA Model of Optimal Portfolio Selection under. Liquidity Risk and Price Impact
A Model of Optimal Portfolio Selection under Liquidity Risk and Price Impact Huyên PHAM Workshop on PDE and Mathematical Finance KTH, Stockholm, August 15, 2005 Laboratoire de Probabilités et Modèles Aléatoires
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationWELL-POSEDNESS FOR HYPERBOLIC PROBLEMS (0.2)
WELL-POSEDNESS FOR HYPERBOLIC PROBLEMS We will use the familiar Hilbert spaces H = L 2 (Ω) and V = H 1 (Ω). We consider the Cauchy problem (.1) c u = ( 2 t c )u = f L 2 ((, T ) Ω) on [, T ] Ω u() = u H
More informationarxiv: v1 [math.ap] 11 Jan 2013
The Hele-Shaw asymptotics for mechanical models of tumor growth Benoît Perthame Fernando Quirós Juan Luis Vázquez arxiv:1301.2612v1 [math.ap] 11 Jan 2013 January 15, 2013 Abstract Models of tumor growth,
More informationContinuum Modeling of Transportation Networks with Differential Equations
with Differential Equations King Abdullah University of Science and Technology Thuwal, KSA Examples of transportation networks The Silk Road Examples of transportation networks Painting by Latifa Echakhch
More informationOn the infinity Laplace operator
On the infinity Laplace operator Petri Juutinen Köln, July 2008 The infinity Laplace equation Gunnar Aronsson (1960 s): variational problems of the form S(u, Ω) = ess sup H (x, u(x), Du(x)). (1) x Ω The
More informationMATH 819 FALL We considered solutions of this equation on the domain Ū, where
MATH 89 FALL. The D linear wave equation weak solutions We have considered the initial value problem for the wave equation in one space dimension: (a) (b) (c) u tt u xx = f(x, t) u(x, ) = g(x), u t (x,
More informationControllability of linear PDEs (I): The wave equation
Controllability of linear PDEs (I): The wave equation M. González-Burgos IMUS, Universidad de Sevilla Doc Course, Course 2, Sevilla, 2018 Contents 1 Introduction. Statement of the problem 2 Distributed
More informationSINGULARITY DETECTION AND PROCESSING WITH WAVELETS
SINGULARITY DETECTION AND PROCESSING WITH WAVELETS Stephane Mallat and Wen Liang Hwang Courant Institute of Mathematical Sciences New York University, New York, NY 10012 Technical Report March 1991 Abstract
More informationSolutions in the sense of distributions. Solutions in the sense of distributions Definition, non uniqueness
Solutions in the sense of distributions Definition, non uniqueness 1. Notion of distributions In order to build weak solutions to the Hopf equation, we need to define derivatives of non smooth functions,
More informationOn the local well-posedness of compressible viscous flows with bounded density
On the local well-posedness of compressible viscous flows with bounded density Marius Paicu University of Bordeaux joint work with Raphaël Danchin and Francesco Fanelli Mathflows 2018, Porquerolles September
More informationStochastic optimal control with rough paths
Stochastic optimal control with rough paths Paul Gassiat TU Berlin Stochastic processes and their statistics in Finance, Okinawa, October 28, 2013 Joint work with Joscha Diehl and Peter Friz Introduction
More informationExample 1. Hamilton-Jacobi equation. In particular, the eikonal equation. for some n( x) > 0 in Ω. Here 1 / 2
Oct. 1 0 Viscosity S olutions In this lecture we take a glimpse of the viscosity solution theory for linear and nonlinear PDEs. From our experience we know that even for linear equations, the existence
More informationWeak solutions to Fokker-Planck equations and Mean Field Games
Weak solutions to Fokker-Planck equations and Mean Field Games Alessio Porretta Universita di Roma Tor Vergata LJLL, Paris VI March 6, 2015 Outlines of the talk Brief description of the Mean Field Games
More information1. Statement of the problem.
218 Нелинейные граничные задачи 18, 218-229 (2008 c 2008. M. A. Borodin THE STEFAN PROBLEM The Stefan problem in its classical statement is a mathematical model of the process of propagation of heat in
More informationSwitching, sparse and averaged control
Switching, sparse and averaged control Enrique Zuazua Ikerbasque & BCAM Basque Center for Applied Mathematics Bilbao - Basque Country- Spain zuazua@bcamath.org http://www.bcamath.org/zuazua/ WG-BCAM, February
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationDerivation of the time dependent Gross-Pitaevskii equation for the dynamics of the Bose-Einstein condensate
Derivation of the time dependent Gross-Pitaevskii equation for the dynamics of the Bose-Einstein condensate László Erdős University of Munich Evian-les-Bains, 2007 June Joint work with B. Schlein and H.T.
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationA Central Limit Theorem for Fleming-Viot Particle Systems Application to the Adaptive Multilevel Splitting Algorithm
A Central Limit Theorem for Fleming-Viot Particle Systems Application to the Algorithm F. Cérou 1,2 B. Delyon 2 A. Guyader 3 M. Rousset 1,2 1 Inria Rennes Bretagne Atlantique 2 IRMAR, Université de Rennes
More informationViscosity Iterative Approximating the Common Fixed Points of Non-expansive Semigroups in Banach Spaces
Viscosity Iterative Approximating the Common Fixed Points of Non-expansive Semigroups in Banach Spaces YUAN-HENG WANG Zhejiang Normal University Department of Mathematics Yingbing Road 688, 321004 Jinhua
More informationarxiv: v1 [math.ap] 10 Apr 2013
QUASI-STATIC EVOLUTION AND CONGESTED CROWD TRANSPORT DAMON ALEXANDER, INWON KIM, AND YAO YAO arxiv:1304.3072v1 [math.ap] 10 Apr 2013 Abstract. We consider the relationship between Hele-Shaw evolution with
More informationNumerical Methods of Applied Mathematics -- II Spring 2009
MIT OpenCourseWare http://ocw.mit.edu 18.336 Numerical Methods of Applied Mathematics -- II Spring 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationA class of Hamilton-Jacobi equations with constraint: uniqueness and constructive approach
A class of Hamilton-Jacobi equations with constraint: uniqueness and constructive approach Sepideh Mirrahimi, Jean-Michel Roquejoffre To cite this version: Sepideh Mirrahimi, Jean-Michel Roquejoffre. A
More informationThere are five problems. Solve four of the five problems. Each problem is worth 25 points. A sheet of convenient formulae is provided.
Preliminary Examination (Solutions): Partial Differential Equations, 1 AM - 1 PM, Jan. 18, 16, oom Discovery Learning Center (DLC) Bechtel Collaboratory. Student ID: There are five problems. Solve four
More informationGeometric projection of stochastic differential equations
Geometric projection of stochastic differential equations John Armstrong (King s College London) Damiano Brigo (Imperial) August 9, 2018 Idea: Projection Idea: Projection Projection gives a method of systematically
More informationb i (µ, x, s) ei ϕ(x) µ s (dx) ds (2) i=1
NONLINEAR EVOLTION EQATIONS FOR MEASRES ON INFINITE DIMENSIONAL SPACES V.I. Bogachev 1, G. Da Prato 2, M. Röckner 3, S.V. Shaposhnikov 1 The goal of this work is to prove the existence of a solution to
More informationPoisson random measure: motivation
: motivation The Lévy measure provides the expected number of jumps by time unit, i.e. in a time interval of the form: [t, t + 1], and of a certain size Example: ν([1, )) is the expected number of jumps
More informationRisk Averse Shape Optimization
Risk Averse Shape Optimization Sergio Conti 2 Martin Pach 1 Martin Rumpf 2 Rüdiger Schultz 1 1 Department of Mathematics University Duisburg-Essen 2 Rheinische Friedrich-Wilhelms-Universität Bonn Workshop
More informationA Necessary and Sufficient Condition for the Continuity of Local Minima of Parabolic Variational Integrals with Linear Growth
A Necessary and Sufficient Condition for the Continuity of Local Minima of Parabolic Variational Integrals with Linear Growth E. DiBenedetto 1 U. Gianazza 2 C. Klaus 1 1 Vanderbilt University, USA 2 Università
More informationMTH 847: PDE I (Fall 2017) Exam 2,
MTH 847: PDE I (Fall 2017) Exam 2, 2017.12.8 Name: Standard exam rules apply: You are not allowed to give or receive help from other students. All electronic devices must be turned off for the duration
More informationA Variational Analysis of a Gauged Nonlinear Schrödinger Equation
A Variational Analysis of a Gauged Nonlinear Schrödinger Equation Alessio Pomponio, joint work with David Ruiz Dipartimento di Meccanica, Matematica e Management, Politecnico di Bari Variational and Topological
More informationTopological complexity of motion planning algorithms
Topological complexity of motion planning algorithms Mark Grant mark.grant@ed.ac.uk School of Mathematics - University of Edinburgh 31st March 2011 Mark Grant (University of Edinburgh) TC of MPAs 31st
More informationChapter 2: First Order DE 2.6 Exact DE and Integrating Fa
Chapter 2: First Order DE 2.6 Exact DE and Integrating Factor First Order DE Recall the general form of the First Order DEs (FODE): dy dx = f(x, y) (1) (In this section x is the independent variable; not
More informationON A MODEL OF A POPULATION WITH VARIABLE MOTILITY
ON A MODEL OF A POPULATION WITH VARIABLE MOTILITY OLGA TURANOVA Abstract. We study a reaction-diffusion equation with a nonlocal reaction term that models a population with variable motility. We establish
More informationViscosity Solutions for Dummies (including Economists)
Viscosity Solutions for Dummies (including Economists) Online Appendix to Income and Wealth Distribution in Macroeconomics: A Continuous-Time Approach written by Benjamin Moll August 13, 2017 1 Viscosity
More informationThe Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition
The Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition Sukjung Hwang CMAC, Yonsei University Collaboration with M. Dindos and M. Mitrea The 1st Meeting of
More informationHomogenization of a 1D pursuit law with delay
Homogenization of a 1D pursuit law with delay Jeremy Firozaly (Joint work with R. Monneau) Paris-Est University June 2, 2016 J. Firozaly (Paris-Est University) Homogenization of a delay equation June 2,
More informationA junction condition by specified homogenization
A junction condition by specified homogenization G. Galise, C. Imbert, R. Monneau June 19, 2014 Abstract Given a coercive Hamiltonian which is quasi-convex with respect to the gradient variable and periodic
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationn 2 xi = x i. x i 2. r r ; i r 2 + V ( r) V ( r) = 0 r > 0. ( 1 1 ) a r n 1 ( 1 2) V( r) = b ln r + c n = 2 b r n 2 + c n 3 ( 1 3)
Sep. 7 The L aplace/ P oisson Equations: Explicit Formulas In this lecture we study the properties of the Laplace equation and the Poisson equation with Dirichlet boundary conditions through explicit representations
More informationCoarsening process in the 2d voter model
Alessandro Tartaglia (LPTHE) Coarsening in the 2d voter model May 8, 2015 1 / 34 Coarsening process in the 2d voter model Alessandro Tartaglia LPTHE, Université Pierre et Marie Curie alessandro.tartaglia91@gmail.com
More informationWELL-POSEDNESS BY NOISE FOR SCALAR CONSERVATION LAWS. 1. Introduction
WELL-POSEDNESS BY NOISE FOR SCALAR CONSERVATION LAWS BENJAMIN GESS AND MARIO MAURELLI Abstract. We consider stochastic scalar conservation laws with spatially inhomogeneous flux. The regularity of the
More informationHomogeniza*ons in Perforated Domain. Ki Ahm Lee Seoul Na*onal University
Homogeniza*ons in Perforated Domain Ki Ahm Lee Seoul Na*onal University Outline 1. Perforated Domain 2. Neumann Problems (joint work with Minha Yoo; interes*ng discussion with Li Ming Yeh) 3. Dirichlet
More informationGENETICS. On the Fixation Process of a Beneficial Mutation in a Variable Environment
GENETICS Supporting Information http://www.genetics.org/content/suppl/2/6/6/genetics..24297.dc On the Fixation Process of a Beneficial Mutation in a Variable Environment Hildegard Uecker and Joachim Hermisson
More informationA path integral approach to the Langevin equation
A path integral approach to the Langevin equation - Ashok Das Reference: A path integral approach to the Langevin equation, A. Das, S. Panda and J. R. L. Santos, arxiv:1411.0256 (to be published in Int.
More informationIntroduction to Random Diffusions
Introduction to Random Diffusions The main reason to study random diffusions is that this class of processes combines two key features of modern probability theory. On the one hand they are semi-martingales
More informationMath 220A - Fall 2002 Homework 5 Solutions
Math 0A - Fall 00 Homework 5 Solutions. Consider the initial-value problem for the hyperbolic equation u tt + u xt 0u xx 0 < x 0 u t (x, 0) ψ(x). Use energy methods to show that the domain of dependence
More informationSTRUCTURAL OPTIMIZATION BY THE LEVEL SET METHOD
Shape optimization 1 G. Allaire STRUCTURAL OPTIMIZATION BY THE LEVEL SET METHOD G. ALLAIRE CMAP, Ecole Polytechnique with F. JOUVE and A.-M. TOADER 1. Introduction 2. Setting of the problem 3. Shape differentiation
More informationLiquid crystal flows in two dimensions
Liquid crystal flows in two dimensions Fanghua Lin Junyu Lin Changyou Wang Abstract The paper is concerned with a simplified hydrodynamic equation, proposed by Ericksen and Leslie, modeling the flow of
More informationEffective dynamics of many-body quantum systems
Effective dynamics of many-body quantum systems László Erdős University of Munich Grenoble, May 30, 2006 A l occassion de soixantiéme anniversaire de Yves Colin de Verdiére Joint with B. Schlein and H.-T.
More informationMalliavin Calculus in Finance
Malliavin Calculus in Finance Peter K. Friz 1 Greeks and the logarithmic derivative trick Model an underlying assent by a Markov process with values in R m with dynamics described by the SDE dx t = b(x
More informationNONTRIVIAL SOLUTIONS TO INTEGRAL AND DIFFERENTIAL EQUATIONS
Fixed Point Theory, Volume 9, No. 1, 28, 3-16 http://www.math.ubbcluj.ro/ nodeacj/sfptcj.html NONTRIVIAL SOLUTIONS TO INTEGRAL AND DIFFERENTIAL EQUATIONS GIOVANNI ANELLO Department of Mathematics University
More informationRegularity of local minimizers of the interaction energy via obstacle problems
Regularity of local minimizers of the interaction energy via obstacle problems J. A. Carrillo, M. G. Delgadino, A. Mellet September 22, 2014 Abstract The repulsion strength at the origin for repulsive/attractive
More informationMathematical modelling of collective behavior
Mathematical modelling of collective behavior Young-Pil Choi Fakultät für Mathematik Technische Universität München This talk is based on joint works with José A. Carrillo, Maxime Hauray, and Samir Salem
More informationMetric Spaces and Topology
Chapter 2 Metric Spaces and Topology From an engineering perspective, the most important way to construct a topology on a set is to define the topology in terms of a metric on the set. This approach underlies
More informationNumerical Solutions of Geometric Partial Differential Equations. Adam Oberman McGill University
Numerical Solutions of Geometric Partial Differential Equations Adam Oberman McGill University Sample Equations and Schemes Fully Nonlinear Pucci Equation! "#+ "#* "#) "#( "#$ "#' "#& "#% "#! "!!!"#$ "
More informationSpaces with Ricci curvature bounded from below
Spaces with Ricci curvature bounded from below Nicola Gigli February 23, 2015 Topics 1) On the definition of spaces with Ricci curvature bounded from below 2) Analytic properties of RCD(K, N) spaces 3)
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More informationDiscretization of Stochastic Differential Systems With Singular Coefficients Part II
Discretization of Stochastic Differential Systems With Singular Coefficients Part II Denis Talay, INRIA Sophia Antipolis joint works with Mireille Bossy, Nicolas Champagnat, Sylvain Maire, Miguel Martinez,
More informationAn Introduction to Malliavin calculus and its applications
An Introduction to Malliavin calculus and its applications Lecture 3: Clark-Ocone formula David Nualart Department of Mathematics Kansas University University of Wyoming Summer School 214 David Nualart
More informationCalculus of Variations. Final Examination
Université Paris-Saclay M AMS and Optimization January 18th, 018 Calculus of Variations Final Examination Duration : 3h ; all kind of paper documents (notes, books...) are authorized. The total score of
More informationThe semi-geostrophic equations - a model for large-scale atmospheric flows
The semi-geostrophic equations - a model for large-scale atmospheric flows Beatrice Pelloni, University of Reading with M. Cullen (Met Office), D. Gilbert, T. Kuna INI - MFE Dec 2013 Introduction - Motivation
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More information