Convergence of a large time-step scheme for Mean Curvature Motion
|
|
- Shanna Hicks
- 5 years ago
- Views:
Transcription
1 Outline Convergence of a large time-step scheme for Mean Curvature Motion M. Falcone Dipartimento di Matematica SAPIENZA, Università di Roma in collaboration with E. Carlini (SAPIENZA) and R. Ferretti (Roma Tre) College de France, June 2, 29
2 Outline Outline Introduction
3 Outline Introduction Introduction
4 Level Set Method The Level Set method has had a great success for the analysis of front propagation problems for its capability to handle many different physical phenomena within the same theoretical framework. One can use it for isotropic and anisotropic front propagation, for merging different fronts, for Mean Curvature Motion (MCM) and other situations when the velocity depends on some geometrical properties of the front. It allows to develop the analysis also after the on set of singularities.
5 Level Set Method Our unknown is the representation function u : R n [,T] R and looking at the -level set of u we can back to the front, i.e. Γ t {x : u(x,t) = } The model equation corresponding to the LS method is { ut + c(x) u(x) = x R n [,T] u(x) = u (x) x R n where u must be a representation function for the front (i.e. u >, x R n \ Ω u(x) = x R n () u(x) < x Ω The front Γ = Ω.
6 Front propagation and minimum time problem Let the velocity of the front in the normal direction is a given c : R n R Minimum time problem Target set =Ω Dynamics { ẏ(t) = c(y)a, a B(,) Minimum time function y() = x (2) T(x) inf{t R + : y x (t;a(t)) Ω }
7 Monotone evolution T( ) is the unique viscosity solution of max a B(,) {c(y)a T(x)} = Rn \ Ω with the Dirichlet condition T(x) = on Ω If c does not change sign the evolution is monotone (increasing or decreasing) and we have the following link u(x,t) = T(x) t So we can solve the stationary problem and get any front Γ t, for t >.
8 More General Models In the standard model the normal velocity c : R n R is given, but the same approach applies to other scalar velocities c(x,t) isotropic growth with time varying velocity c(x,η) anisotropic growth, cristal growth c(x,k(x)) Mean Curvature Motion c(x) obtained by convolution (dislocation dynamics)
9 Outline Introduction Introduction
10 SL schemes for MCM u t (x,t) = div u(x,) = u (x) Representation formula (Soner Touzi): ( ) Du(x,t) Du(x,t) Du(x, t) u(x,t) = E{u (y(x,t,t))}, Du { dy(x,t,s) = 2P(y,t,s)dW (s) y(x,t,) = x P(y,t,s) = ( u 2 ) x 2 u x u x2 Du 2 u x u x2 ux 2
11 Soner Touzi formula between t and t + t : u(x,t + t) = E{u(y(x,t + t, t),t)} Brownian dimension reduction 2PdW = 2 Du ( ux2 u x )( ux2 dw Du u ) x dw 2 = σdŵ Du where σ(x,t) = ( 2 ux2 Du u x )
12 we can replace the stochastic Cauchy problem by { dy(x,t,s) = σ(y,t,s)dŵ (s) y(x,t,) = x. i.e. with a -dimensional brownian motion in the direction tangent to the curve.
13 Euler scheme for SDE with Time-discretization { yk+ = y k + 2σ(y,t k,) Ŵk y = x. P( Ŵk = ± t) = 2. u t (x,t n+ ) = 2 u t(x + 2σ(x,t n,) t,t n ) + Fully discrete scheme Du u n+ j = u t(x 2σ(x,t n,) t,t n ). ( I[u n ](x j + σ n j t) + I[u n ](x j σ n j t) )
14 Some references MCM via viscosity solution techniques Evans-Spruck, J. Diff. Geom., 99 Evans - Soner-Souganidis, Comm. Pure Appl. Math, 992,... Kohn-Serfaty, Comm. Pure Appl. Math., 26 Soner-Touzi, J. Eur. Math. Soc, 22 and Ann. Prob. 23 Numerical methods Osher-Sethian, JCP, 988 Merriman-Bence-Osher, AMS LN, 993 Crandall-Lions, Numer. Math., 996 Barles-Georgelin, SIAM J. Num. Anal., 995 Catté-Dibos-Koepfler, SIAM J. Num. Anal., 995 F. Ferretti, ENUMATH 2 Carlini- F. -Ferretti, JCP 25 and Interfaces and Free Boundaries, 29 (submitted)
15 Outline Introduction Introduction
16 Modified scheme with threshold C x s u n+ j = 2 u n+ j = 4 ui n i D(j) [ ] I[u n ](x j + σj n t) + I[u n ](x j σj n t) if D n j > C xs if D n j C xs Consistency error (case Dj n > C xs ) ( ) ( x r x q s τ x, t = O + O t t 2 ) + O( t 2) + O( t)
17 Consistency for small gradients Let us consider the case D n j C xs. Case a: (x j,t n ) is such that Du(x j,t n ) =. This is a stardad computation based on the consistency with the heat equation. Case b : (x j,t n ) is such that Du(x j,t n ) and D j [w] C x s. By the lower semicontinuity of F we have that for any ε > there exists a δ (ε ) such that F(Du,D 2 u)(y,s) F(Du,D 2 u)(x,t) ε for any (y,s) B δ (x,t). (3)
18 Consistency for small gradients The upper semicontinuity of F implies that for any ε 2 > there exists a δ 2 (ε 2 ) such that F(Du,D 2 u)(y,s) F(Du,D 2 u)(x,t)+ε 2 for any (y,s) B δ2 (x,t). (4) The interesting case is when for ε max(ε,ε 2 ) there exists (y,s) B δ (x,t) B δ2 (x,t) such that Du(y,s) =. For (x,t) = (x j,t n ), we can apply both (3) and (4) getting F(Du,D 2 u)(x j,t n ) F(Du,D 2 u)(x j,t n ) + 2ε (5)
19 Consistency for small gradients In fact, we have = liminf t F(Du,D 2 u)(x j,t n ) 2 D 2 u(y,s) + ε ε = u(x j,t n) H(w;j) t + ε = limsup t u(x j,t n) H(w;j) t 2 D 2 u(y,s) + ε F(Du,D 2 u)(x j,t n ) + 2ε and, since ε and ε 2 are arbitrary, we obtain the result. + ε
20 Monotonicity The standard scheme presented at the beginning is not monotone even local monotone interpolation operators (linear, bilinear). We modify it in order to satisfy a relaxed monotonicity property. u n+ j u n j t = ρ 2 ( 2 I[un ](x j + σ n j ρ)+ 2 I[un ](x j σ n j ρ) un j ) u n+ j = 4 i D(j) In compact form we write u n i if D n j C xs. u n+ j = H ρ (u n ;j).
21 Monotonicity Since for this scheme H ρ we have H ρ (u n ;j) u n i 4 tl I[u n ],G Cρ x s+ we need to compensate the negative bound and we add a small viscosity P H ρ (u n ;j) = H ρ (u n ;j) + t W x i D(j) un i 4un j ρ x, if D n s x 2 j > C xs u n+ j = 4 i D(j) where W is a positive constant. u n i if D n j C xs.
22 Monotonicity Differentiating H ρ we get: H ρ (u n ;j) uj n t 4W t ρ2 ρ x +s H ρ (u n ;j) = ψ i (x j + σj n ρ) for i S(j) \ (D(j) {j}) u n i H ρ (u n ;j) u n i 4 tl I[u n ],G W t + Cρ xs+ ρ x +s Then H ρ (u n ;j) is monotone if t 4W t ρ2 ρ x +s 4L I[u n ],G < W. C for i D(j). (6) (7)
23 Weak monotonicity property Theorem Under the above assumption, the scheme H ρ satisfies H ρ H ρ (η;j) + o( t).
24 Outline Introduction Introduction
25 Let us consider now a general scheme that is supposed to approximate our HJ equation. Its abstract form on a lattice will be given by u n+ j = S Dt (u n ;j) for j Z 2 n =,...,N u j = u (x j ) for j Z 2 (HJ Dt ) where S Dt : B(G x ) R and B(D) is the space of the bounded functions defined on D.
26 The scheme has to satisfy the following conditions: A - Invariance with respect to the addition of constants For any k R and j Z 2, S Dt (v + k;j) = S Dt (v;j) + k (A) A2 - Weak consistency Let us define F(Dφ,D 2 φ)(x,t) = liminf (y,s) (x,t) F(Dφ,D2 φ)(y,s), F(Dφ,D 2 φ)(x,t) = limsup F(Dφ,D 2 φ)(y,s). (y,s) (x,t)
27 Consistency The consistency assumption (A2) requires φ t (x,t) + F(Dφ,D 2 φ)(x,t) φ(x j,t n ) S Dt (φ n ;j) liminf (x j,t n) (x,t) t t φ(x j,t n ) S Dt (φ n ;j) lim sup (x j,t n) (x,t) t t φ t (x,t) + F(Dφ,D 2 φ)(x,t) where φ C (R 2 (,T]) and φ n = (φ(x j,t n )) xj G x. Note that F(Dφ,D 2 φ), F(Dφ,D 2 φ)(x,t) are respectively lower and upper semicontinuous extension of F.
28 Note that if F is continuous, then the liminf and the limsup must coincide, and the definition reduces to the usual definition of consistency. A3 - Generalized monotonicity v j φ n j for j Z 2 implies S Dt (v;j) S Dt (φ n ;j)+o(dt) where v B(G x ) and S Dt is a (possibly different) scheme weakly consistent.
29 Then, consider u n = (uj n) j Z 2 with un j solution of (HJ Dt ) and its piecewise constant (in time) interpolation u Dt defined as: { u Dt I[u n ](x) if t [t n,t n+ ), (x,t) = u (x) if t [, t), where I[ ] : B(G x ) R is a general interpolation operator I[u n ](x) ψ l (x)ul n (8) l I(x) where ψ l (x) are basis functions in R 2 and I(x) is the set of indices corresponding to the vertices of the cell containing x.
30 Let us note that we can always couple the choice of Dt and Dx according to Dt = CDx γ, with C a positive constant, in order to deal with a unique discretization parameter. We assume γ (this parameter has to be tuned to guarantee monotonicity). The interpolation operator I[ ] has to verify a relaxed monotonicity property: if v j η j for any j I(x) then I[v](x) I[η](x)+o(Dt) (A4) with v B(G x ) and η = (f (x j )) xj G x,where f (x) is a smooth function.
31 Moreover I[ ] satisfies I[η](x) f (x) = o(dt) for any x R 2. (A5) Theorem Assume (A) (A5) and let u(x,t) be the unique viscosity solution of (HJ). Then u Dt (x,t) u(x,t) locally uniformly on R 2 [,T] as Dt.
32 Outline Introduction Introduction
33 TEST: a shrinking circle t = O( x 3 4), s = 4. x t L order L order Errors for the SL scheme
34 TEST: Development of non empty interior Fattening: evolution of the level curves u =.95,,.5
35 The torus evolving into a sphere, R < r
36 The torus collapsing in a circle, R > r
37 Dumb-bell: topology change in R 3
38 Crandall-Lions scheme In the scheme proposed by such problems are avoided by replacing the matrix Θ in the equation of MCM by the following one: with ǫ >. Θ ǫ (Du) = I Du Du Du 2 + ǫ,
39 Crandall-Lions scheme cntd Denoting by e j (j =,2) the canonical base of R 2, we can write the CL scheme as: = H CL (u n ;j) where u n+ j H CL (u n ;j) = u n j + t ρ 2 (I[un ](x j + ρθ ǫ (D j [u n ])e )+ +I[u n ](x j + ρθ ǫ (D j [u n ])e 2 ) + +I[u n ](x j ρθ ǫ (D j [u n ])e ) + + I[u n ](x j ρθ ǫ (D j [u n ])e 2 ) uj n )
40 The scheme for which convergence is proved is u n+ j = H CL (u n ;j) + tk ui n 4uj n. ρ x i D(j)
41 Kohn-Serfaty scheme It has been proved that we can write ( ) Du(x,t) { } div Du(x, t) = min a T D 2 u(x,t)a, (9) Du(x, t) a S,a Du= By penalization we can write min {a max T D 2 u(x,t)a ε a Du,aT D 2 u(x,t)a + ε } a Du a S
42 Kohn-Serfaty scheme We write the scheme as u n+ j = H KS (u n ;j), where { ( H KS (w;j) = min max I[w] x j + ) ( 2 t a,i[w] x j )} 2 t a a S The average of the (SL) scheme is replaced by a min-max operator. This scheme is more expensive and difficult to extend to higher dimension.
43 Test : Evolution of a circle (shrinking) Errors for the SL scheme x t order order CPU s s s
44 Test : Evolution of a circle (shrinking) Errors for the min max scheme x t order order CPU s s m29s ms
45 Outline Introduction Introduction
46 Cuves evolving in the space Let us consider the evolution of a curve C in R 3. There are two ways to handle the problem: The curve is described as the intersection of two surfaces This leads to a system of HJ equation (Osher at alia) The curve is replaced by an ε-tube centered at the curve C. We study the evolution of the surface and get back to the curve in the limit for ε. Following the second characterization, the (SL) scheme has been extended to codimension-2 problems.
47 Cuves evolving in the space { u t = F(D 2 u,du) R 3 [, ) u(x,) = 2d(x, C)2 where F(A,p) = inf ν N(p) {trace[apν ]} and N(p) = {ν S 2 : P ν p = }, and P ν = νν T Soner Touzi formula u(x,t) = inf ν U {E{u(y ν(x,t,t),)}} with U = {ν : [,T] S 2 R 3 : ν Du(x,t) = }
48 The generalized characteristics curves solve { dy ν (x,t,s) = 2ν(s)dŴ(s) y ν (x,t,) = x L.Ambrosio, M.Soner, Level Set Approach to Mean Curvature Flow in Arbitrary Codimension, J.Differential Geometry, 43 (996),
49 MCM cod 2: Numerical approximation Time-discrete scheme u t (x,t n+ ) = 2 inf ν n b U Fully-discrete scheme { u t (x + 2 tν n,t n ) + u t (x } 2 tν n,t n ) u n j = min ν n R 3{ 2 I[un ](x + 2 tν n ) + 2 I[un ](x 2 tν n ) + + ( Dn j νn ) 2 ǫ + ( νn ) 2 ǫ 2 }.
50 Evolution of ǫ helical surface,ǫ =.8.
51 Evolution of two linked circles in R 3 (ǫ-tube).
52 Optimal Trajectory Algorithm s := x jmin ν n j point of the curve = argmin{u n (s j + 2 sνj n )} j =,...,ĵ νj n Û s j+ = s j + 2 sνj n s sĵ ǫ j =,...,ĵ
53 Evolution of a helix in R 3.
54 Evolution of two linked circles in R 3.
Fast Marching Methods for Front Propagation Lecture 1/3
Outline for Front Propagation Lecture 1/3 M. Falcone Dipartimento di Matematica SAPIENZA, Università di Roma Autumn School Introduction to Numerical Methods for Moving Boundaries ENSTA, Paris, November,
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 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 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 informationMean Field Games on networks
Mean Field Games on networks Claudio Marchi Università di Padova joint works with: S. Cacace (Rome) and F. Camilli (Rome) C. Marchi (Univ. of Padova) Mean Field Games on networks Roma, June 14 th, 2017
More informationGAME-THEORETIC SCHEMES FOR GENERALIZED CURVATURE FLOWS IN THE PLANE
Sixth Mississippi State Conference on Differential Equations and Computational Simulations, Electronic Journal of Differential Equations, Conference 5 (7), pp. 5 64. ISSN: 7-669. URL: http://ejde.math.txstate.edu
More informationAn adaptive RBF-based Semi-Lagrangian scheme for HJB equations
An adaptive RBF-based Semi-Lagrangian scheme for HJB equations Roberto Ferretti Università degli Studi Roma Tre Dipartimento di Matematica e Fisica Numerical methods for Hamilton Jacobi equations in optimal
More informationStochastic Homogenization for Reaction-Diffusion Equations
Stochastic Homogenization for Reaction-Diffusion Equations Jessica Lin McGill University Joint Work with Andrej Zlatoš June 18, 2018 Motivation: Forest Fires ç ç ç ç ç ç ç ç ç ç Motivation: Forest Fires
More informationConvergence of a first order scheme for a non local eikonal equation
Convergence of a first order scheme for a non local eikonal equation O. Alvarez, E. Carlini, R. Monneau, E. Rouy March 22, 2005 Abstract We prove the convergence of a first order finite difference scheme
More informationGLOBAL EXISTENCE RESULTS AND UNIQUENESS FOR DISLOCATION EQUATIONS
GLOBAL EXISTENCE RESULTS AND UNIQUENESS FOR DISLOCATION EQUATIONS GUY BARLES, PIERRE CARDALIAGUET, OLIVIER LEY & RÉGIS MONNEAU Abstract. We are interested in nonlocal Eikonal Equations arising in the study
More informationConvergence of a Generalized Fast Marching Method for an Eikonal equation with a Velocity Changing Sign
Convergence of a Generalized Fast Marching Method for an Eikonal equation with a Velocity Changing Sign E. Carlini, M. Falcone, N. Forcadel, R. Monneau December 28, 2007 Abstract We present a new Fast
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 informationOn level-set approach to motion of manifolds of arbitrary codimension
Interfaces and Free Boundaries 5 2003, 417 458 On level-set approach to motion of manifolds of arbitrary codimension DEJAN SLEPČEV Department of Mathematics, University of Toronto, Toronto, ON M5S 3G3,
More informationA generalized MBO diffusion generated motion for constrained harmonic maps
A generalized MBO diffusion generated motion for constrained harmonic maps Dong Wang Department of Mathematics, University of Utah Joint work with Braxton Osting (U. Utah) Workshop on Modeling and Simulation
More informationRough Burgers-like equations with multiplicative noise
Rough Burgers-like equations with multiplicative noise Martin Hairer Hendrik Weber Mathematics Institute University of Warwick Bielefeld, 3.11.21 Burgers-like equation Aim: Existence/Uniqueness for du
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationECOLE POLYTECHNIQUE. Approximation of the anisotropic mean curvature ow CENTRE DE MATHÉMATIQUES APPLIQUÉES UMR CNRS A. Chambolle, M.
ECOLE POLYTECHNIQUE CENTRE DE MATHÉMATIQUES APPLIQUÉES UMR CNRS 7641 91128 PALAISEAU CEDEX (FRANCE). Tél: 01 69 33 41 50. Fax: 01 69 33 30 11 http://www.cmap.polytechnique.fr/ Approximation of the anisotropic
More informationContinuous dependence estimates for the ergodic problem with an application to homogenization
Continuous dependence estimates for the ergodic problem with an application to homogenization Claudio Marchi Bayreuth, September 12 th, 2013 C. Marchi (Università di Padova) Continuous dependence Bayreuth,
More informationMean field games and related models
Mean field games and related models Fabio Camilli SBAI-Dipartimento di Scienze di Base e Applicate per l Ingegneria Facoltà di Ingegneria Civile ed Industriale Email: Camilli@dmmm.uniroma1.it Web page:
More informationThe Hopf - Lax solution for state dependent Hamilton - Jacobi equations
The Hopf - Lax solution for state dependent Hamilton - Jacobi equations Italo Capuzzo Dolcetta Dipartimento di Matematica, Università di Roma - La Sapienza 1 Introduction Consider the Hamilton - Jacobi
More informationA Narrow-Stencil Finite Difference Method for Hamilton-Jacobi-Bellman Equations
A Narrow-Stencil Finite Difference Method for Hamilton-Jacobi-Bellman Equations Xiaobing Feng Department of Mathematics The University of Tennessee, Knoxville, U.S.A. Linz, November 23, 2016 Collaborators
More informationPresenter: Noriyoshi Fukaya
Y. Martel, F. Merle, and T.-P. Tsai, Stability and Asymptotic Stability in the Energy Space of the Sum of N Solitons for Subcritical gkdv Equations, Comm. Math. Phys. 31 (00), 347-373. Presenter: Noriyoshi
More informationPDEs in Image Processing, Tutorials
PDEs in Image Processing, Tutorials Markus Grasmair Vienna, Winter Term 2010 2011 Direct Methods Let X be a topological space and R: X R {+ } some functional. following definitions: The mapping R is lower
More informationDynamic Blocking Problems for Models of Fire Propagation
Dynamic Blocking Problems for Models of Fire Propagation Alberto Bressan Department of Mathematics, Penn State University bressan@math.psu.edu Alberto Bressan (Penn State) Dynamic Blocking Problems 1 /
More informationZdzislaw Brzeźniak. Department of Mathematics University of York. joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York)
Navier-Stokes equations with constrained L 2 energy of the solution Zdzislaw Brzeźniak Department of Mathematics University of York joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York) LMS
More informationEXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN EQUATION WITH DEGENERATE MOBILITY
Electronic Journal of Differential Equations, Vol. 216 216), No. 329, pp. 1 22. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN
More informationMath The Laplacian. 1 Green s Identities, Fundamental Solution
Math. 209 The Laplacian Green s Identities, Fundamental Solution Let be a bounded open set in R n, n 2, with smooth boundary. The fact that the boundary is smooth means that at each point x the external
More informationGlobal existence results for eikonal equation with BV initial data
Global existence results for eikonal equation with BV initial data. Boudjerada 1, A. El Hajj 2 July 8, 2014 Abstract In this paper, we study a local and a non-local eikonal equations in one dimensional
More informationMulti-dimensional Stochastic Singular Control Via Dynkin Game and Dirichlet Form
Multi-dimensional Stochastic Singular Control Via Dynkin Game and Dirichlet Form Yipeng Yang * Under the supervision of Dr. Michael Taksar Department of Mathematics University of Missouri-Columbia Oct
More informationHomogenization and error estimates of free boundary velocities in periodic media
Homogenization and error estimates of free boundary velocities in periodic media Inwon C. Kim October 7, 2011 Abstract In this note I describe a recent result ([14]-[15]) on homogenization and error estimates
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 225 Estimates of the second-order derivatives for solutions to the two-phase parabolic problem
More informationThe main motivation of this paper comes from the following, rather surprising, result of Ecker and Huisken [13]: for any initial data u 0 W 1,
Quasilinear parabolic equations, unbounded solutions and geometrical equations II. Uniqueness without growth conditions and applications to the mean curvature flow in IR 2 Guy Barles, Samuel Biton and
More informationADJOINT METHODS FOR OBSTACLE PROBLEMS AND WEAKLY COUPLED SYSTEMS OF PDE. 1. Introduction
ADJOINT METHODS FOR OBSTACLE PROBLEMS AND WEAKLY COPLED SYSTEMS OF PDE F. CAGNETTI, D. GOMES, AND H.V. TRAN Abstract. The adjoint method, recently introduced by Evans, is used to study obstacle problems,
More informationMinimal Surface equations non-solvability strongly convex functional further regularity Consider minimal surface equation.
Lecture 7 Minimal Surface equations non-solvability strongly convex functional further regularity Consider minimal surface equation div + u = ϕ on ) = 0 in The solution is a critical point or the minimizer
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
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 informationON SELECTION OF SOLUTIONS TO VECTORIAL HAMILTON-JACOBI SYSTEM
ON SELECTION OF SOLUTIONS TO VECTORIAL HAMILTON-JACOBI SYSTEM BAISHENG YAN Abstract. In this paper, we discuss some principles on the selection of the solutions of a special vectorial Hamilton-Jacobi system
More informationMEAN CURVATURE FLOW OF ENTIRE GRAPHS EVOLVING AWAY FROM THE HEAT FLOW
MEAN CURVATURE FLOW OF ENTIRE GRAPHS EVOLVING AWAY FROM THE HEAT FLOW GREGORY DRUGAN AND XUAN HIEN NGUYEN Abstract. We present two initial graphs over the entire R n, n 2 for which the mean curvature flow
More informationA Semi-Lagrangian Scheme for the Open Table Problem in Granular Matter Theory
A Semi-Lagrangian Scheme for the Open Table Problem in Granular Matter Theory M. Falcone and S. Finzi Vita Abstract We introduce and analyze a new scheme for the approximation of the two-layer model proposed
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 informationSingular Diffusion Equations With Nonuniform Driving Force. Y. Giga University of Tokyo (Joint work with M.-H. Giga) July 2009
Singular Diffusion Equations With Nonuniform Driving Force Y. Giga University of Tokyo (Joint work with M.-H. Giga) July 2009 1 Contents 0. Introduction 1. Typical Problems 2. Variational Characterization
More informationConvergence of Finite Volumes schemes for an elliptic-hyperbolic system with boundary conditions
Convergence of Finite Volumes schemes for an elliptic-hyperbolic system with boundary conditions Marie Hélène Vignal UMPA, E.N.S. Lyon 46 Allée d Italie 69364 Lyon, Cedex 07, France abstract. We are interested
More informationZdzislaw Brzeźniak. Department of Mathematics University of York. joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York)
Navier-Stokes equations with constrained L 2 energy of the solution Zdzislaw Brzeźniak Department of Mathematics University of York joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York) Stochastic
More informationSolution Sheet 3. Solution Consider. with the metric. We also define a subset. and thus for any x, y X 0
Solution Sheet Throughout this sheet denotes a domain of R n with sufficiently smooth boundary. 1. Let 1 p
More informationDetermination of thin elastic inclusions from boundary measurements.
Determination of thin elastic inclusions from boundary measurements. Elena Beretta in collaboration with E. Francini, S. Vessella, E. Kim and J. Lee September 7, 2010 E. Beretta (Università di Roma La
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 informationSelf-intersection local time for Gaussian processes
Self-intersection local Olga Izyumtseva, olaizyumtseva@yahoo.com (in collaboration with Andrey Dorogovtsev, adoro@imath.kiev.ua) Department of theory of random processes Institute of mathematics Ukrainian
More informationOptimal stopping time formulation of adaptive image filtering
Optimal stopping time formulation of adaptive image filtering I. Capuzzo Dolcetta, R. Ferretti 19.04.2000 Abstract This paper presents an approach to image filtering based on an optimal stopping time problem
More informationCentre for Mathematics and Its Applications The Australian National University Canberra, ACT 0200 Australia. 1. Introduction
ON LOCALLY CONVEX HYPERSURFACES WITH BOUNDARY Neil S. Trudinger Xu-Jia Wang Centre for Mathematics and Its Applications The Australian National University Canberra, ACT 0200 Australia Abstract. In this
More informationExistence of weak solutions for general nonlocal and nonlinear second-order parabolic equations
Existence of weak solutions for general nonlocal and nonlinear second-order parabolic equations Guy Barles, Pierre Cardaliaguet, Olivier Ley, Aurélien Monteillet To cite this version: Guy Barles, Pierre
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationApplications of the compensated compactness method on hyperbolic conservation systems
Applications of the compensated compactness method on hyperbolic conservation systems Yunguang Lu Department of Mathematics National University of Colombia e-mail:ylu@unal.edu.co ALAMMI 2009 In this talk,
More informationLYAPUNOV STABILITY OF CLOSED SETS IN IMPULSIVE SEMIDYNAMICAL SYSTEMS
Electronic Journal of Differential Equations, Vol. 2010(2010, No. 78, pp. 1 18. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu LYAPUNOV STABILITY
More informationΓ-convergence of functionals on divergence-free fields
Γ-convergence of functionals on divergence-free fields N. Ansini Section de Mathématiques EPFL 05 Lausanne Switzerland A. Garroni Dip. di Matematica Univ. di Roma La Sapienza P.le A. Moro 2 0085 Rome,
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 informationComputing High Frequency Waves By the Level Set Method
Computing High Frequency Waves By the Level Set Method Hailiang Liu Department of Mathematics Iowa State University Collaborators: Li-Tien Cheng (UCSD), Stanley Osher (UCLA) Shi Jin (UW-Madison), Richard
More informationUNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE
Surveys in Mathematics and its Applications ISSN 1842-6298 (electronic), 1843-7265 (print) Volume 5 (2010), 275 284 UNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE Iuliana Carmen Bărbăcioru Abstract.
More informationThe HJB-POD approach for infinite dimensional control problems
The HJB-POD approach for infinite dimensional control problems M. Falcone works in collaboration with A. Alla, D. Kalise and S. Volkwein Università di Roma La Sapienza OCERTO Workshop Cortona, June 22,
More informationBrownian Motion. 1 Definition Brownian Motion Wiener measure... 3
Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................
More informationA finite element level set method for anisotropic mean curvature flow with space dependent weight
A finite element level set method for anisotropic mean curvature flow with space dependent weight Klaus Deckelnick and Gerhard Dziuk Centre for Mathematical Analysis and Its Applications, School of Mathematical
More informationA Ginzburg-Landau Type Problem for Nematics with Highly Anisotropic Elastic Term
A Ginzburg-Landau Type Problem for Nematics with Highly Anisotropic Elastic Term Peter Sternberg In collaboration with Dmitry Golovaty (Akron) and Raghav Venkatraman (Indiana) Department of Mathematics
More informationL -uniqueness of Schrödinger operators on a Riemannian manifold
L -uniqueness of Schrödinger operators on a Riemannian manifold Ludovic Dan Lemle Abstract. The main purpose of this paper is to study L -uniqueness of Schrödinger operators and generalized Schrödinger
More informationSome asymptotic properties of solutions for Burgers equation in L p (R)
ARMA manuscript No. (will be inserted by the editor) Some asymptotic properties of solutions for Burgers equation in L p (R) PAULO R. ZINGANO Abstract We discuss time asymptotic properties of solutions
More informationOn continuous time contract theory
Ecole Polytechnique, France Journée de rentrée du CMAP, 3 octobre, 218 Outline 1 2 Semimartingale measures on the canonical space Random horizon 2nd order backward SDEs (Static) Principal-Agent Problem
More informationGENERAL EXISTENCE OF SOLUTIONS TO DYNAMIC PROGRAMMING PRINCIPLE. 1. Introduction
GENERAL EXISTENCE OF SOLUTIONS TO DYNAMIC PROGRAMMING PRINCIPLE QING LIU AND ARMIN SCHIKORRA Abstract. We provide an alternative approach to the existence of solutions to dynamic programming equations
More informationWeak convergence and Brownian Motion. (telegram style notes) P.J.C. Spreij
Weak convergence and Brownian Motion (telegram style notes) P.J.C. Spreij this version: December 8, 2006 1 The space C[0, ) In this section we summarize some facts concerning the space C[0, ) of real
More informationOn Stopping Times and Impulse Control with Constraint
On Stopping Times and Impulse Control with Constraint Jose Luis Menaldi Based on joint papers with M. Robin (216, 217) Department of Mathematics Wayne State University Detroit, Michigan 4822, USA (e-mail:
More informationAdaptive evolution : a population approach Benoît Perthame
Adaptive evolution : a population approach Benoît Perthame 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).
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationCOMPARISON PRINCIPLES FOR CONSTRAINED SUBHARMONICS PH.D. COURSE - SPRING 2019 UNIVERSITÀ DI MILANO
COMPARISON PRINCIPLES FOR CONSTRAINED SUBHARMONICS PH.D. COURSE - SPRING 2019 UNIVERSITÀ DI MILANO KEVIN R. PAYNE 1. Introduction Constant coefficient differential inequalities and inclusions, constraint
More informationThuong Nguyen. SADCO Internal Review Metting
Asymptotic behavior of singularly perturbed control system: non-periodic setting Thuong Nguyen (Joint work with A. Siconolfi) SADCO Internal Review Metting Rome, Nov 10-12, 2014 Thuong Nguyen (Roma Sapienza)
More informationWeek 2 Notes, Math 865, Tanveer
Week 2 Notes, Math 865, Tanveer 1. Incompressible constant density equations in different forms Recall we derived the Navier-Stokes equation for incompressible constant density, i.e. homogeneous flows:
More informationExtremal Solutions of Differential Inclusions via Baire Category: a Dual Approach
Extremal Solutions of Differential Inclusions via Baire Category: a Dual Approach Alberto Bressan Department of Mathematics, Penn State University University Park, Pa 1682, USA e-mail: bressan@mathpsuedu
More informationSYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS. Massimo Grosi Filomena Pacella S. L. Yadava. 1. Introduction
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 21, 2003, 211 226 SYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS Massimo Grosi Filomena Pacella S.
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 informationMATH 6605: SUMMARY LECTURE NOTES
MATH 6605: SUMMARY LECTURE NOTES These notes summarize the lectures on weak convergence of stochastic processes. If you see any typos, please let me know. 1. Construction of Stochastic rocesses A stochastic
More informationMATH 220: MIDTERM OCTOBER 29, 2015
MATH 22: MIDTERM OCTOBER 29, 25 This is a closed book, closed notes, no electronic devices exam. There are 5 problems. Solve Problems -3 and one of Problems 4 and 5. Write your solutions to problems and
More informationLocal semiconvexity of Kantorovich potentials on non-compact manifolds
Local semiconvexity of Kantorovich potentials on non-compact manifolds Alessio Figalli, Nicola Gigli Abstract We prove that any Kantorovich potential for the cost function c = d / on a Riemannian manifold
More informationPATH FUNCTIONALS OVER WASSERSTEIN SPACES. Giuseppe Buttazzo. Dipartimento di Matematica Università di Pisa.
PATH FUNCTIONALS OVER WASSERSTEIN SPACES Giuseppe Buttazzo Dipartimento di Matematica Università di Pisa buttazzo@dm.unipi.it http://cvgmt.sns.it ENS Ker-Lann October 21-23, 2004 Several natural structures
More informationOn the approximation of the principal eigenvalue for a class of nonlinear elliptic operators
On the approximation of the principal eigenvalue for a class of nonlinear elliptic operators Fabio Camilli ("Sapienza" Università di Roma) joint work with I.Birindelli ("Sapienza") I.Capuzzo Dolcetta ("Sapienza")
More informationLOCAL BEHAVIOR AND GLOBAL EXISTENCE OF POSITIVE SOLUTIONS OF au λ u u λ. COMPORTEMENT LOCAL ET EXISTENCE GLOBALE DES SOLUTIONS POSITIVES DE au λ u u λ
Ann. I. H. Poincaré AN 19, 6 (2002) 889 901 2002 Éditions scientifiques et médicales Elsevier SAS. All rights reserved S0294-1449(02)00105-1/FLA LOCAL BEHAVIOR AND GLOBAL EXISTENCE OF POSITIVE SOLUTIONS
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 informationarxiv: v3 [math.ds] 9 Nov 2012
Positive expansive flows Alfonso Artigue November 13, 2018 arxiv:1210.3202v3 [math.ds] 9 Nov 2012 Abstract We show that every positive expansive flow on a compact metric space consists of a finite number
More informationShape from Shading with discontinuous image brightness
Shape from Shading with discontinuous image brightness Fabio Camilli Dip. di Matematica Pura e Applicata, Università dell Aquila (Italy), camilli@ing.univaq.it Emmanuel Prados UCLA Vision Lab. (USA), eprados@cs.ucla.edu
More informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
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 informationLipschitz continuity for solutions of Hamilton-Jacobi equation with Ornstein-Uhlenbeck operator
Lipschitz continuity for solutions of Hamilton-Jacobi equation with Ornstein-Uhlenbeck operator Thi Tuyen Nguyen Ph.D student of University of Rennes 1 Joint work with: Prof. E. Chasseigne(University of
More informationOptimal Transportation. Nonlinear Partial Differential Equations
Optimal Transportation and Nonlinear Partial Differential Equations Neil S. Trudinger Centre of Mathematics and its Applications Australian National University 26th Brazilian Mathematical Colloquium 2007
More informationA Nonlinear PDE in Mathematical Finance
A Nonlinear PDE in Mathematical Finance Sergio Polidoro Dipartimento di Matematica, Università di Bologna, Piazza di Porta S. Donato 5, 40127 Bologna (Italy) polidoro@dm.unibo.it Summary. We study a non
More informationA stochastic representation for the level set equations
A stochastic representation for the level set equations H. Mete Soner Koç University Department of Mathematics Rumelifener Yolu, Sariyer Istanbul 80910, Turkey msoner@ku.edu.tr Nizar Touzi CREST 15 Bd
More informationA CONNECTION BETWEEN A GENERAL CLASS OF SUPERPARABOLIC FUNCTIONS AND SUPERSOLUTIONS
A CONNECTION BETWEEN A GENERAL CLASS OF SUPERPARABOLIC FUNCTIONS AND SUPERSOLUTIONS RIIKKA KORTE, TUOMO KUUSI, AND MIKKO PARVIAINEN Abstract. We show to a general class of parabolic equations that every
More informationSome issues on Electrical Impedance Tomography with complex coefficient
Some issues on Electrical Impedance Tomography with complex coefficient Elisa Francini (Università di Firenze) in collaboration with Elena Beretta and Sergio Vessella E. Francini (Università di Firenze)
More informationReaction-Diffusion Equations In Narrow Tubes and Wave Front P
Outlines Reaction-Diffusion Equations In Narrow Tubes and Wave Front Propagation University of Maryland, College Park USA Outline of Part I Outlines Real Life Examples Description of the Problem and Main
More informationFormulation of the problem
TOPICAL PROBLEMS OF FLUID MECHANICS DOI: https://doi.org/.43/tpfm.27. NOTE ON THE PROBLEM OF DISSIPATIVE MEASURE-VALUED SOLUTIONS TO THE COMPRESSIBLE NON-NEWTONIAN SYSTEM H. Al Baba, 2, M. Caggio, B. Ducomet
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 informationWeek 6 Notes, Math 865, Tanveer
Week 6 Notes, Math 865, Tanveer. Energy Methods for Euler and Navier-Stokes Equation We will consider this week basic energy estimates. These are estimates on the L 2 spatial norms of the solution u(x,
More informationRichard F. Bass Krzysztof Burdzy University of Washington
ON DOMAIN MONOTONICITY OF THE NEUMANN HEAT KERNEL Richard F. Bass Krzysztof Burdzy University of Washington Abstract. Some examples are given of convex domains for which domain monotonicity of the Neumann
More informationClosed-Loop Impulse Control of Oscillating Systems
Closed-Loop Impulse Control of Oscillating Systems A. N. Daryin and A. B. Kurzhanski Moscow State (Lomonosov) University Faculty of Computational Mathematics and Cybernetics Periodic Control Systems, 2007
More informationIntegro-differential equations: Regularity theory and Pohozaev identities
Integro-differential equations: Regularity theory and Pohozaev identities Xavier Ros Oton Departament Matemàtica Aplicada I, Universitat Politècnica de Catalunya PhD Thesis Advisor: Xavier Cabré Xavier
More informationSYMMETRY OF POSITIVE SOLUTIONS OF SOME NONLINEAR EQUATIONS. M. Grossi S. Kesavan F. Pacella M. Ramaswamy. 1. Introduction
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 12, 1998, 47 59 SYMMETRY OF POSITIVE SOLUTIONS OF SOME NONLINEAR EQUATIONS M. Grossi S. Kesavan F. Pacella M. Ramaswamy
More information