From nonlocal to local Cahn-Hilliard equation. Stefano Melchionna Helene Ranetbauer Lara Trussardi. Uni Wien (Austria) September 18, 2018

Size: px
Start display at page:

Download "From nonlocal to local Cahn-Hilliard equation. Stefano Melchionna Helene Ranetbauer Lara Trussardi. Uni Wien (Austria) September 18, 2018"

Transcription

1 From nonlocal to local Cahn-Hilliard equation Stefano Melchionna Helene Ranetbauer Lara Trussardi Uni Wien (Austria) September 18, 2018 SFB P D ME S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 1 / 26

2 Index 1 Introduction 2 Local CH 3 Nonlocal CH 4 Main result S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 2 / 26

3 Cahn-Hilliard equation It has been proposed in 1958 by [Cahn, Hilliard 1958, Cahn 1961] It describes the process of phase separation (spinodal decomposition) in binary alloys (iron-nickel) Phase field model u [0, 1] (vs sharp interface model u {0, 1}) It has a variety of applications: image processing [Capuzzo Dolcetta, Finzi Vita, March 2002] population dynamics [Cohen, Murray 1981] formation of Saturn rings [Tremaine 2003] tumour growth [Colli, Garcke, Gilardi, Lam, Rocca, Sprekels, Scala... ] S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 3 / 26

4 Settings u R: real valued function representing the local concentration of one of the two components u = 0, u = 1: pure phases Two species A, B with concentrations c A and c B = 1 c A at each point. Then u(x, t) = 1 c A (x, t) and: u 1 pure phase B; u 0 pure phase A. Source of the picture: Wikipedia S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 4 / 26

5 Settings u R: real valued function representing the local concentration of one of the two components u = 0, u = 1: pure phases Two species A, B with concentrations c A and c B = 1 c A at each point. Then u(x, t) = 1 c A (x, t) and: u 1 pure phase B; u 0 pure phase A. Source of the picture: Wikipedia S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 4 / 26

6 Settings u R: real valued function representing the local concentration of one of the two components u = 0, u = 1: pure phases Two species A, B with concentrations c A and c B = 1 c A at each point. Then u(x, t) = 1 c A (x, t) and: u 1 pure phase B; u 0 pure phase A. Source of the picture: Wikipedia S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 4 / 26

7 Settings u R: real valued function representing the local concentration of one of the two components u = 0, u = 1: pure phases Two species A, B with concentrations c A and c B = 1 c A at each point. Then u(x, t) = 1 c A (x, t) and: u 1 pure phase B; u 0 pure phase A. Source of the picture: Wikipedia S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 4 / 26

8 Settings u R: real valued function representing the local concentration of one of the two components u = 0, u = 1: pure phases Two species A, B with concentrations c A and c B = 1 c A at each point. Then u(x, t) = 1 c A (x, t) and: u 1 pure phase B; u 0 pure phase A. Source of the picture: Wikipedia S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 4 / 26

9 Settings u R: real valued function representing the local concentration of one of the two components u = 0, u = 1: pure phases Two species A, B with concentrations c A and c B = 1 c A at each point. Then u(x, t) = 1 c A (x, t) and: u 1 pure phase B; u 0 pure phase A. Source of the picture: Wikipedia S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 4 / 26

10 Settings u R: real valued function representing the local concentration of one of the two components u = 0, u = 1: pure phases Two species A, B with concentrations c A and c B = 1 c A at each point. Then u(x, t) = 1 c A (x, t) and: u 1 pure phase B; u 0 pure phase A. Source of the picture: Wikipedia S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 4 / 26

11 Settings u R: real valued function representing the local concentration of one of the two components u = 0, u = 1: pure phases Two species A, B with concentrations c A and c B = 1 c A at each point. Then u(x, t) = 1 c A (x, t) and: u 1 pure phase B; u 0 pure phase A. Source of the picture: Wikipedia S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 4 / 26

12 Settings u R: real valued function representing the local concentration of one of the two components u = 0, u = 1: pure phases Two species A, B with concentrations c A and c B = 1 c A at each point. Then u(x, t) = 1 c A (x, t) and: u 1 pure phase B; u 0 pure phase A. Source of the picture: Wikipedia S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 4 / 26

13 Settings u R: real valued function representing the local concentration of one of the two components u = 0, u = 1: pure phases Two species A, B with concentrations c A and c B = 1 c A at each point. Then u(x, t) = 1 c A (x, t) and: u 1 pure phase B; u 0 pure phase A. Source of the picture: Wikipedia S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 4 / 26

14 Settings u R: real valued function representing the local concentration of one of the two components u = 0, u = 1: pure phases Two species A, B with concentrations c A and c B = 1 c A at each point. Then u(x, t) = 1 c A (x, t) and: u 1 pure phase B; u 0 pure phase A. Source of the picture: Wikipedia S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 4 / 26

15 Settings u R: real valued function representing the local concentration of one of the two components u = 0, u = 1: pure phases Two species A, B with concentrations c A and c B = 1 c A at each point. Then u(x, t) = 1 c A (x, t) and: u 1 pure phase B; u 0 pure phase A. Source of the picture: Wikipedia S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 4 / 26

16 Settings u R: real valued function representing the local concentration of one of the two components u = 0, u = 1: pure phases Two species A, B with concentrations c A and c B = 1 c A at each point. Then u(x, t) = 1 c A (x, t) and: u 1 pure phase B; u 0 pure phase A. Source of the picture: Wikipedia S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 4 / 26

17 Settings u R: real valued function representing the local concentration of one of the two components u = 0, u = 1: pure phases Two species A, B with concentrations c A and c B = 1 c A at each point. Then u(x, t) = 1 c A (x, t) and: u 1 pure phase B; u 0 pure phase A. Source of the picture: Wikipedia S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 4 / 26

18 Settings u R: real valued function representing the local concentration of one of the two components u = 0, u = 1: pure phases Two species A, B with concentrations c A and c B = 1 c A at each point. Then u(x, t) = 1 c A (x, t) and: u 1 pure phase B; u 0 pure phase A. Source of the picture: Wikipedia S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 4 / 26

19 Local Cahn-Hilliard Free energy functional E CH (u) = where ( τ 2 ) 2 u 2 + F(u) dx τ: small positive parameter related to the transition region thickness F: double well potential with two global minima in the pure phases u 2 : reflects intermolecular interactions (penalising the creation of interfaces) S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 5 / 26

20 Local Cahn-Hilliard Corresponding evolution problem (4 th order PDE): where u t + J CH = 0, µ: mobility (constant = 1) J CH = µ(u) v CH, v CH = δe CH(u) δu v = δe δu : chemical potential = τ 2 u + F (u) S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 6 / 26

21 Nonlocal Cahn-Hilliard Proposed by [Giacomin, Lebowitz 1997] Free energy functional E NL (u) = 1 K (x, y)(u(x) u(y)) 2 dxdy + 4 F(u(x))dx, where K (x, y): positive and symmetric convolution kernel F: double well potential with two global minima in the pure phases S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 7 / 26

22 Nonlocal Cahn-Hilliard Corresponding evolution problem (2 nd order PDE): u t + J NL = 0, J NL = µ(u) v NL, v NL = δe NL(u) δu = (K 1)u K u + F (u) S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 8 / 26

23 Nonlocal Cahn-Hilliard Corresponding evolution problem (2 nd order PDE): u t + J NL = 0, J NL = µ(u) v NL, v NL = δe NL(u) δu = (K 1)u K u + F (u) Goal: prove the convergence of solutions of the nonlocal Cahn-Hilliard equation to solutions of the local version in a periodic setting. S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 8 / 26

24 Local vs nonlocal They share fundamental features: underlying gradient flow structure, lack of comparison principle, separation from the pure phases,... Both energy functionals allow the same Γ-limit for vanishing interface thickness [Gal, Grasselli, Miranville, Rocca,... ] Pointwise convergence is of little use due to non convexity and lack of coercivity of the nonlocal energy functional E NL in H 1. There exists Γ-convergence for the energy functionals [Ponce 2004] but it is not trivial to prove convergence for solutions of the corresponding dynamic problems using [Sandier, Serfaty 2011] S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 9 / 26

25 Overview of known results and properties Eq F µ Existence Separation Long time pol non-deg Garcke 00, Temam probably false Temam 88 deg Elliott Garcke 96?? CH log non-deg Elliott Luckhaus 91 Miranville Zelik 04 Cherfils Zelik Miranville 11 deg Elliott Garcke 96? Debussche Dettori 95 pol non-deg Bates Han 04 probably false Gal Grasselli 17 deg??? NLCH log non-deg Gal Giorgini Grasselli 17 Gal Giorgini Grasselli 17 Gal Giorgini Grasselli 17 deg Gajewski Zacharias 03 Londen Petzeltova 11 Londen Petzeltova 11 S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 10 / 26

26 Hypothesis H1 d-dimensional flat torus with d 3 S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 11 / 26

27 Hypothesis H1 d-dimensional flat torus with d 3 H2 Family of convolution kernels parametrised by a parameter ε: ( K ε (x, y) = ε d 2 x y J 2 ) ε with J : R R sufficiently smooth nonnegative function with compact support and 1 J( z 2 ) z 2 dz = 1 d S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 11 / 26

28 Hypothesis H1 d-dimensional flat torus with d 3 H2 Family of convolution kernels parametrised by a parameter ε: ( K ε (x, y) = ε d 2 x y J 2 ) ε with J : R R sufficiently smooth nonnegative function with compact support and 1 J( z 2 ) z 2 dz = 1 d H3 F C 2 (R) double well potential with two global minima at 0 and 1 such that F (s) 0 for s (, a] [a, + ) with a nonnegative, and C l ( u 3 + 1) F (u) C u ( u 3 + 1) for C l, C u > 0 S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 11 / 26

29 Hypothesis H1 d-dimensional flat torus with d 3 H2 Family of convolution kernels parametrised by a parameter ε: ( K ε (x, y) = ε d 2 x y J 2 ) ε with J : R R sufficiently smooth nonnegative function with compact support and 1 J( z 2 ) z 2 dz = 1 d H3 F C 2 (R) double well potential with two global minima at 0 and 1 such that F (s) 0 for s (, a] [a, + ) with a nonnegative, and C l ( u 3 + 1) F (u) C u ( u 3 + 1) for C l, C u > 0 H4 u 0,ε L 2 () converges strongly in L 2 () to the limit u 0 H 1 () and satisfies E ε (u 0,ε ), E(u 0 ) C 0 for some constant C 0 > 0 independent of ε S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 11 / 26

30 Definition of solutions Definition (Weak solution to the nonlocal Cahn-Hilliard equation) Let ε > 0 and T > 0 be fixed. We define u ε to be a weak solution to the nonlocal Cahn-Hilliard equation on [0, T ] associated with the initial datum u 0,ε L 2 () if u ε H 1 (0, T ; (H 1 ()) ) L 2 (0, T ; H 1 ()), satisfies t u ε, ϕ (H 1 ()),H 1 () + for all ϕ H 1 (), and u ε (0) = u 0,ε. [(K ε 1)u ε K ε u ε + F (u ε )] ϕ dx = 0 S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 12 / 26

31 Definition of solutions Definition (Weak solution to the local Cahn-Hilliard equation) Let T > 0 be fixed. We define u to be a weak solution to the Cahn-Hilliard equation on [0, T ] associated with the initial datum u 0 H 1 () if u H 1 (0, T ; (H 1 ()) ) L 2 (0, T ; H 2 ()), satisfies t u, ϕ (H 1 ()),H 1 () + for all ϕ H 2 (), and u(0) = u 0. u ϕ dx F (u) ϕ dx = 0 S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 13 / 26

32 Main result Remark: existence and uniqueness of weak solutions to both problems are well known with different choices for the boundary conditions both systems have been largely studied (qualitative properties, numerical aspects, long-time behaviour, asymptotics with different kinds of boundary conditions and different potentials) S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 14 / 26

33 Main result Remark: existence and uniqueness of weak solutions to both problems are well known with different choices for the boundary conditions both systems have been largely studied (qualitative properties, numerical aspects, long-time behaviour, asymptotics with different kinds of boundary conditions and different potentials) Theorem Let u ε be a solution of the nonlocal CH with periodic boundary conditions and kernel K ε (x, y). Then u ε u in L 2 (0, T ; H 1 ()) H 1 (0, T ; H 1 ()) where u is a solution of the local CH. S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 14 / 26

34 Some comments Note that in the Neumann case CH has two boundary conditions (one for u and one for the chemical potential), while NLCH has just one (for the chemical potential v). The mobility can be annoying for the estimates on the chemical potential. S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 15 / 26

35 Some comments Note that in the Neumann case CH has two boundary conditions (one for u and one for the chemical potential), while NLCH has just one (for the chemical potential v). The mobility can be annoying for the estimates on the chemical potential. Idea: use advantage of the dynamic structure: for every fixed ε, u ε H 1 S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 15 / 26

36 Proof: uniform estimates t u ε, ϕ (H 1 ()),H 1 () + [(K ε 1)u ε K ε u ε + F (u ε )] ϕ dx = 0 Test function: ϕ = u ε 0 = 1 d 2 dt u ε 2 L 2 () + [ (K ε 1)u ε K ε u ε + F (u ε ) ] u ε dx S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 16 / 26

37 Proof: uniform estimates t u ε, ϕ (H 1 ()),H 1 () + [(K ε 1)u ε K ε u ε + F (u ε )] ϕ dx = 0 Test function: ϕ = u ε 0 = 1 d 2 dt u ε 2 L 2 () + 0 = 1 d 2 dt u ε 2 L 2 () + [ (K ε 1)u ε K ε u ε + F (u ε ) ] u ε dx [(K ε 1) u ε 2 (K ε u ε ) u ε +F (u ε ) u ε 2 ] dx }{{} 1 2 Kε(x,y) uε(x) uε(y) 2 dx dy S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 16 / 26

38 Proof: uniform estimates Change of variable x y ε =: z 1 d 2 dt u ε 2 L 2 () + 1 J( z 2 ) u ε (y + εz) u ε (y) 2 ε = F (u ε ) u ε 2 dx B 1 u ε 2 L 2 (). 2 dy dz S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 17 / 26

39 Proof: uniform estimates Change of variable x y ε =: z 1 d 2 dt u ε 2 L 2 () + 1 J( z 2 ) u ε (y + εz) u ε (y) 2 ε = F (u ε ) u ε 2 dx B 1 u ε 2 L 2 (). 2 dy dz Aim: estimate the blue term with u ε 2 Key idea: [Ponce 2004], Poincaré type inequality S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 17 / 26

40 Key idea R N with N 1, bounded domain with Lipschitz boundary, 1 p < Poincaré It exist C p > 0 s.t. f f p C p Df p, f W 1,p () Let (ρ n ) L 1 (R N ) be a sequence of radial functions satisfying: ρ n 0 a.e. in R N ρ n = 1 n 1 R N ρ n (h)dh = 0 δ > 0 lim n h >δ S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 18 / 26

41 Theorem 1 Let (ρ n ) L 1 (R N ) be a sequence of radial functions as defined before. Given δ > 0, there exists n 0 1 sufficiently large, such that ( f f p Cp ) f (x) f (y) p + δ K p,n x y p ρ n ( x y ) dx dy for every f L p () and n n 0. S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 19 / 26

42 Theorem 1 Let (ρ n ) L 1 (R N ) be a sequence of radial functions as defined before. Given δ > 0, there exists n 0 1 sufficiently large, such that ( f f p Cp ) f (x) f (y) p + δ K p,n x y p ρ n ( x y ) dx dy for every f L p () and n n 0. Observe: this formulation is stronger than Poincaré. It is easy to see that f (x) f (y) p x y p ρ n ( x y ) dx dy Df p C R N Df p S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 19 / 26

43 Theorem 2 compactness If (f n ) L p (R n ) is a bounded sequence such that f n (x) f n (y) p x y p ρ n ( x y )dxdy B, n 1 then (f n ) is relatively compact in L p. Assume that f n f in L p () then f W 1,p () if 1 < p < f BV () if p = 1 S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 20 / 26

44 Proof: uniform estimates t u ε, ϕ (H 1 ()),H 1 () + [(K ε 1)u ε K ε u ε + F (u ε )] ϕ dx = 0 Test function: ϕ = ( ) 1 U ε where ( ) 1 : (H 1 ()) H 1 () is the map assigning to every v (H 1 ()) the unique solution w of the equation w = v such that the mean value is zero, i.e. w = 0. We define U ε = u ε u ε. S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 21 / 26

45 Proof: uniform estimates t u ε, ϕ (H 1 ()),H 1 () + [(K ε 1)u ε K ε u ε + F (u ε )] ϕ dx = 0 Test function: ϕ = ( ) 1 U ε where ( ) 1 : (H 1 ()) H 1 () is the map assigning to every v (H 1 ()) the unique solution w of the equation w = v such that the mean value is zero, i.e. w = 0. We define U ε = u ε u ε. u ε u ε 2 L 2 () C p and u ε u ε 2 L 2 () C p ( ) 2 J( z 2 ) z 2 Uε (y + εz) U ε (y) dy dz B 2 ε z J( z 2 ) z 2 u ε (y + εz) u ε (y) ε z 2 dy dz, S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 21 / 26

46 Uniform estimates 1 d 2 dt u ε 2 L 2 () J( z 2 ) u ε (y + εz) u ε (y) ε 2 dy dz B 1 u ε 2 L 2 () S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 22 / 26

47 Uniform estimates 1 d 2 dt u ε 2 L 2 () J( z 2 ) u ε (y + εz) u ε (y) ε 2 dy dz B 1 u ε 2 L 2 () ( ) u ε(t ) 2 L 2 () + B 1 u ε 2 2C L 2 (0,T ;L 2 ()) 1 p 2 u ε(0) 2 L 2 () S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 22 / 26

48 Uniform estimates 1 d 2 dt u ε 2 L 2 () J( z 2 ) u ε (y + εz) u ε (y) ε 2 dy dz B 1 u ε 2 L 2 () ( ) u ε(t ) 2 L 2 () + B 1 u ε 2 2C L 2 (0,T ;L 2 ()) 1 p 2 u ε(0) 2 L 2 () T 0 J( z 2 ) u ε (y + εz) u ε (y) ε 2 dy dz dt C S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 22 / 26

49 Proof: convergence Goal: prove the limit u to be a weak solution of the local Cahn-Hilliard equation u ε u weakly in L 2 (0, T ; H 1 ()) t u ε t u weakly in L 2 (0, T ; (H 1 ()) ) u ε u strongly in C([0, T ]; L 2 ()) for some limit u L 2 (0, T ; H 1 ()) H 1 (0, T ; (H 1 ()) ) S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 23 / 26

50 Proof: convergence Goal: prove the limit u to be a weak solution of the local Cahn-Hilliard equation u ε u weakly in L 2 (0, T ; H 1 ()) t u ε t u weakly in L 2 (0, T ; (H 1 ()) ) u ε u strongly in C([0, T ]; L 2 ()) for some limit u L 2 (0, T ; H 1 ()) H 1 (0, T ; (H 1 ()) ) [Ponce 2004] u L 2 (0, T ; H 2 ()) S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 23 / 26

51 Proof: convergence Test function: ϕ C () T 0 = ( t u ε )ϕ dx dt 0 }{{} I T F (u ε ) ϕ dx dt 0 } {{ } II 1 T K ε (x, y)(u ε (x) u ε (y))( ϕ(x) ϕ(y)) dy dx dt 2 0 }{{} III (I) (II): the growth conditions and continuity on F suffice to pass to the limit S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 24 / 26

52 Proof: convergence Thanks to (H2): 1 d Hd 1 (S d 1 ) 0 J(r 2 )r d+1 dr = 1 and by using the weak convergence: 1 T J( z 2 ) z (u ε (y + εz) u ε (y)) ( ϕ(y + εz) ϕ(y)) dy dzdt ε z ε z 1 T u(y) ϕ(y)dy dt 2 Thus, the limit u satisfies T T ( t u)ϕ dx dt T u ϕ dx dt F (u) ϕ dx dt = 0. 0 S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 25 / 26

53 Outlook and open questions Proved the convergence of weak solutions of the NLCH equation to the weak solutions to the local one as the convolution kernel approximates a Dirac delta (case with periodic boundary conditions) using a compactness argument. Dirichlet / Neumann boundary conditions? S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 26 / 26

54 Outlook and open questions Proved the convergence of weak solutions of the NLCH equation to the weak solutions to the local one as the convolution kernel approximates a Dirac delta (case with periodic boundary conditions) using a compactness argument. Dirichlet / Neumann boundary conditions? Thanks for your attention SFB P D ME Melchionna S., Ranetbauer H. and T. L., From nonlocal to local Cahn-Hilliard equation, submitted (2018) S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal to local Cahn-Hilliard equation 26 / 26

THE CAHN-HILLIARD EQUATION WITH A LOGARITHMIC POTENTIAL AND DYNAMIC BOUNDARY CONDITIONS

THE CAHN-HILLIARD EQUATION WITH A LOGARITHMIC POTENTIAL AND DYNAMIC BOUNDARY CONDITIONS THE CAHN-HILLIARD EQUATION WITH A LOGARITHMIC POTENTIAL AND DYNAMIC BOUNDARY CONDITIONS Alain Miranville Université de Poitiers, France Collaborators : L. Cherfils, G. Gilardi, G.R. Goldstein, G. Schimperna,

More information

THE CAUCHY PROBLEM AND STEADY STATE SOLUTIONS FOR A NONLOCAL CAHN-HILLIARD EQUATION

THE CAUCHY PROBLEM AND STEADY STATE SOLUTIONS FOR A NONLOCAL CAHN-HILLIARD EQUATION Electronic Journal of Differential Equations, Vol. 24(24), No. 113, pp. 1 9. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) THE CAUCHY

More information

EXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN EQUATION WITH DEGENERATE MOBILITY

EXISTENCE 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 information

Well-posedness and asymptotic analysis for a Penrose-Fife type phase-field system

Well-posedness and asymptotic analysis for a Penrose-Fife type phase-field system 0-0 Well-posedness and asymptotic analysis for a Penrose-Fife type phase-field system Buona positura e analisi asintotica per un sistema di phase field di tipo Penrose-Fife Salò, 3-5 Luglio 2003 Riccarda

More information

Diffuse and sharp interfaces in Biology and Mechanics

Diffuse and sharp interfaces in Biology and Mechanics Diffuse and sharp interfaces in Biology and Mechanics E. Rocca Università degli Studi di Pavia SIMAI 2016 Politecnico of Milan, September 13 16, 2016 Supported by the FP7-IDEAS-ERC-StG Grant EntroPhase

More information

Una aproximación no local de un modelo para la formación de pilas de arena

Una aproximación no local de un modelo para la formación de pilas de arena Cabo de Gata-2007 p. 1/2 Una aproximación no local de un modelo para la formación de pilas de arena F. Andreu, J.M. Mazón, J. Rossi and J. Toledo Cabo de Gata-2007 p. 2/2 OUTLINE The sandpile model of

More information

Sobolev Spaces. Chapter 10

Sobolev Spaces. Chapter 10 Chapter 1 Sobolev Spaces We now define spaces H 1,p (R n ), known as Sobolev spaces. For u to belong to H 1,p (R n ), we require that u L p (R n ) and that u have weak derivatives of first order in L p

More information

Weak solutions for the Cahn-Hilliard equation with degenerate mobility

Weak solutions for the Cahn-Hilliard equation with degenerate mobility Archive for Rational Mechanics and Analysis manuscript No. (will be inserted by the editor) Shibin Dai Qiang Du Weak solutions for the Cahn-Hilliard equation with degenerate mobility Abstract In this paper,

More information

Numerical Approximation of Phase Field Models

Numerical Approximation of Phase Field Models Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School

More information

PDEs, Homework #3 Solutions. 1. Use Hölder s inequality to show that the solution of the heat equation

PDEs, Homework #3 Solutions. 1. Use Hölder s inequality to show that the solution of the heat equation PDEs, Homework #3 Solutions. Use Hölder s inequality to show that the solution of the heat equation u t = ku xx, u(x, = φ(x (HE goes to zero as t, if φ is continuous and bounded with φ L p for some p.

More information

A sharp diffuse interface tracking method for approximating evolving interfaces

A sharp diffuse interface tracking method for approximating evolving interfaces A sharp diffuse interface tracking method for approximating evolving interfaces Vanessa Styles and Charlie Elliott University of Sussex Overview Introduction Phase field models Double well and double obstacle

More information

(d). Why does this imply that there is no bounded extension operator E : W 1,1 (U) W 1,1 (R n )? Proof. 2 k 1. a k 1 a k

(d). Why does this imply that there is no bounded extension operator E : W 1,1 (U) W 1,1 (R n )? Proof. 2 k 1. a k 1 a k Exercise For k 0,,... let k be the rectangle in the plane (2 k, 0 + ((0, (0, and for k, 2,... let [3, 2 k ] (0, ε k. Thus is a passage connecting the room k to the room k. Let ( k0 k (. We assume ε k

More information

On some weighted fractional porous media equations

On some weighted fractional porous media equations On some weighted fractional porous media equations Gabriele Grillo Politecnico di Milano September 16 th, 2015 Anacapri Joint works with M. Muratori and F. Punzo Gabriele Grillo Weighted Fractional PME

More information

Asymptotic behavior of the degenerate p Laplacian equation on bounded domains

Asymptotic behavior of the degenerate p Laplacian equation on bounded domains Asymptotic behavior of the degenerate p Laplacian equation on bounded domains Diana Stan Instituto de Ciencias Matematicas (CSIC), Madrid, Spain UAM, September 19, 2011 Diana Stan (ICMAT & UAM) Nonlinear

More information

NONLOCAL DIFFUSION EQUATIONS

NONLOCAL DIFFUSION EQUATIONS NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,

More information

Optimal Control in Diffuse Interface Models of Tumor Growth

Optimal Control in Diffuse Interface Models of Tumor Growth Optimal Control in Diffuse Interface Models of Tumor Growth E. Rocca Università degli Studi di Pavia WIAS-Berlin May 4, 2017 joint work with Harald Garcke and Kei Fong Lam (Regensburg) Fondazione Cariplo

More information

Homogenization of Lattice Systems

Homogenization of Lattice Systems ANDREA BRAIDES Homogenization of Lattice Systems Ginzburg-Landau equations, Dislocations and Homogenization May 23, 2011, Ile de Ré From discrete to continuous energies Discrete system: with discrete variables

More information

Fractional Cahn-Hilliard equation

Fractional Cahn-Hilliard equation Fractional Cahn-Hilliard equation Goro Akagi Mathematical Institute, Tohoku University Abstract In this note, we review a recent joint work with G. Schimperna and A. Segatti (University of Pavia, Italy)

More information

In this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,

In this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots, Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical

More information

Small energy regularity for a fractional Ginzburg-Landau system

Small energy regularity for a fractional Ginzburg-Landau system Small energy regularity for a fractional Ginzburg-Landau system Yannick Sire University Aix-Marseille Work in progress with Vincent Millot (Univ. Paris 7) The fractional Ginzburg-Landau system We are interest

More information

Free energy estimates for the two-dimensional Keller-Segel model

Free energy estimates for the two-dimensional Keller-Segel model Free energy estimates for the two-dimensional Keller-Segel model dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine in collaboration with A. Blanchet (CERMICS, ENPC & Ceremade) & B.

More information

HOW TO APPROXIMATE THE HEAT EQUATION WITH NEUMANN BOUNDARY CONDITIONS BY NONLOCAL DIFFUSION PROBLEMS. 1. Introduction

HOW TO APPROXIMATE THE HEAT EQUATION WITH NEUMANN BOUNDARY CONDITIONS BY NONLOCAL DIFFUSION PROBLEMS. 1. Introduction HOW TO APPROXIMATE THE HEAT EQUATION WITH NEUMANN BOUNDARY CONDITIONS BY NONLOCAL DIFFUSION PROBLEMS CARMEN CORTAZAR, MANUEL ELGUETA, ULIO D. ROSSI, AND NOEMI WOLANSKI Abstract. We present a model for

More information

Global minimization. Chapter Upper and lower bounds

Global minimization. Chapter Upper and lower bounds Chapter 1 Global minimization The issues related to the behavior of global minimization problems along a sequence of functionals F are by now well understood, and mainly rely on the concept of -limit.

More information

Calculus of Variations. Final Examination

Calculus 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 information

Nonlocal Cahn-Hilliard equations with fractional dynamic boundary conditions

Nonlocal Cahn-Hilliard equations with fractional dynamic boundary conditions Nonlocal Cahn-Hilliard equations with fractional dynamic boundary conditions Ciprian G. Gal To cite this version: Ciprian G. Gal. Nonlocal Cahn-Hilliard equations with fractional dynamic boundary conditions.

More information

Some lecture notes for Math 6050E: PDEs, Fall 2016

Some lecture notes for Math 6050E: PDEs, Fall 2016 Some lecture notes for Math 65E: PDEs, Fall 216 Tianling Jin December 1, 216 1 Variational methods We discuss an example of the use of variational methods in obtaining existence of solutions. Theorem 1.1.

More information

On Multigrid for Phase Field

On Multigrid for Phase Field On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis

More information

Applications of the periodic unfolding method to multi-scale problems

Applications of the periodic unfolding method to multi-scale problems Applications of the periodic unfolding method to multi-scale problems Doina Cioranescu Université Paris VI Santiago de Compostela December 14, 2009 Periodic unfolding method and multi-scale problems 1/56

More information

Analysis of a herding model in social economics

Analysis of a herding model in social economics Analysis of a herding model in social economics Lara Trussardi 1 Ansgar Jüngel 1 C. Kühn 1 1 Technische Universität Wien Taormina - June 13, 2014 www.itn-strike.eu L. Trussardi, A. Jüngel, C. Kühn (TUW)

More information

Modelling of interfaces and free boundaries

Modelling 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 information

Sobolev spaces. May 18

Sobolev 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 information

BOUNDARY FLUXES FOR NON-LOCAL DIFFUSION

BOUNDARY FLUXES FOR NON-LOCAL DIFFUSION BOUNDARY FLUXES FOR NON-LOCAL DIFFUSION CARMEN CORTAZAR, MANUEL ELGUETA, JULIO D. ROSSI, AND NOEMI WOLANSKI Abstract. We study a nonlocal diffusion operator in a bounded smooth domain prescribing the flux

More information

Optimal control in diffuse interface models of tumor growth

Optimal control in diffuse interface models of tumor growth Optimal control in diffuse interface models of tumor growth E. Rocca Università degli Studi di Pavia 37ème Colloque de la Société Francophone de Biologie Théorique, Poitiers joint work with Harald Garcke

More information

arxiv: v1 [math.ap] 11 Apr 2019

arxiv: v1 [math.ap] 11 Apr 2019 Global well-posedness of the 3D Cahn-Hilliard equations Zhenbang Li Caifeng Liu arxiv:194.6191v1 [math.ap] 11 Apr 219 April 15, 219 Abstract The Cauchy problem of the Cahn-Hilliard equations is studied

More information

Diffusion on the half-line. The Dirichlet problem

Diffusion on the half-line. The Dirichlet problem Diffusion on the half-line The Dirichlet problem Consider the initial boundary value problem (IBVP) on the half line (, ): v t kv xx = v(x, ) = φ(x) v(, t) =. The solution will be obtained by the reflection

More information

with deterministic and noise terms for a general non-homogeneous Cahn-Hilliard equation Modeling and Asymptotics

with deterministic and noise terms for a general non-homogeneous Cahn-Hilliard equation Modeling and Asymptotics 12-3-2009 Modeling and Asymptotics for a general non-homogeneous Cahn-Hilliard equation with deterministic and noise terms D.C. Antonopoulou (Joint with G. Karali and G. Kossioris) Department of Applied

More information

Slow motion for the nonlocal Allen Cahn equation in n-dimensions

Slow motion for the nonlocal Allen Cahn equation in n-dimensions Slow motion for the nonlocal Allen Cahn equation in n-dimensions Ryan Murray Carnegie Mellon University Pittsburgh, PA, USA Matteo Rinaldi Carnegie Mellon University Pittsburgh, PA, USA December 4, 215

More information

THE STEINER REARRANGEMENT IN ANY CODIMENSION

THE STEINER REARRANGEMENT IN ANY CODIMENSION THE STEINER REARRANGEMENT IN ANY CODIMENSION GIUSEPPE MARIA CAPRIANI Abstract. We analyze the Steiner rearrangement in any codimension of Sobolev and BV functions. In particular, we prove a Pólya-Szegő

More information

Laplace s Equation. Chapter Mean Value Formulas

Laplace 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 information

PDEs in Image Processing, Tutorials

PDEs 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 information

Math The Laplacian. 1 Green s Identities, Fundamental Solution

Math 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 information

Institut für Mathematik

Institut für Mathematik RHEINISCH-WESTFÄLISCHE TECHNISCHE HOCHSCHULE AACHEN Institut für Mathematik Travelling Wave Solutions of the Heat Equation in Three Dimensional Cylinders with Non-Linear Dissipation on the Boundary by

More information

Dynamical systems with Gaussian and Levy noise: analytical and stochastic approaches

Dynamical systems with Gaussian and Levy noise: analytical and stochastic approaches Dynamical systems with Gaussian and Levy noise: analytical and stochastic approaches Noise is often considered as some disturbing component of the system. In particular physical situations, noise becomes

More information

UNCONDITIONALLY ENERGY STABLE SCHEME FOR CAHN-MORRAL EQUATION 3

UNCONDITIONALLY ENERGY STABLE SCHEME FOR CAHN-MORRAL EQUATION 3 UNCONDITIONALLY ENERGY STABLE TIME STEPPING SCHEME FOR CAHN-MORRAL EQUATION: APPLICATION TO MULTI-COMPONENT SPINODAL DECOMPOSITION AND OPTIMAL SPACE TILING R. TAVAKOLI Abstract. An unconditionally energy

More information

OPTIMAL POTENTIALS FOR SCHRÖDINGER OPERATORS. 1. Introduction In this paper we consider optimization problems of the form. min F (V ) : V V, (1.

OPTIMAL POTENTIALS FOR SCHRÖDINGER OPERATORS. 1. Introduction In this paper we consider optimization problems of the form. min F (V ) : V V, (1. OPTIMAL POTENTIALS FOR SCHRÖDINGER OPERATORS G. BUTTAZZO, A. GEROLIN, B. RUFFINI, AND B. VELICHKOV Abstract. We consider the Schrödinger operator + V (x) on H 0 (), where is a given domain of R d. Our

More information

2 A Model, Harmonic Map, Problem

2 A Model, Harmonic Map, Problem ELLIPTIC SYSTEMS JOHN E. HUTCHINSON Department of Mathematics School of Mathematical Sciences, A.N.U. 1 Introduction Elliptic equations model the behaviour of scalar quantities u, such as temperature or

More information

GLOBAL SOLUTION TO THE ALLEN-CAHN EQUATION WITH SINGULAR POTENTIALS AND DYNAMIC BOUNDARY CONDITIONS

GLOBAL SOLUTION TO THE ALLEN-CAHN EQUATION WITH SINGULAR POTENTIALS AND DYNAMIC BOUNDARY CONDITIONS GLOBAL SOLUTION TO THE ALLEN-CAHN EQUATION WITH SINGULAR POTENTIALS AND DYNAMIC BOUNDARY CONDITIONS Luca Calatroni Cambridge Centre for Analysis, University of Cambridge Wilberforce Road, CB3 0WA, Cambridge,

More information

Entropy-dissipation methods I: Fokker-Planck equations

Entropy-dissipation methods I: Fokker-Planck equations 1 Entropy-dissipation methods I: Fokker-Planck equations Ansgar Jüngel Vienna University of Technology, Austria www.jungel.at.vu Introduction Boltzmann equation Fokker-Planck equations Degenerate parabolic

More information

LECTURE 3: DISCRETE GRADIENT FLOWS

LECTURE 3: DISCRETE GRADIENT FLOWS LECTURE 3: DISCRETE GRADIENT FLOWS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA Tutorial: Numerical Methods for FBPs Free Boundary Problems and

More information

A generalized MBO diffusion generated motion for constrained harmonic maps

A 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 information

Pseudo-Poincaré Inequalities and Applications to Sobolev Inequalities

Pseudo-Poincaré Inequalities and Applications to Sobolev Inequalities Pseudo-Poincaré Inequalities and Applications to Sobolev Inequalities Laurent Saloff-Coste Abstract Most smoothing procedures are via averaging. Pseudo-Poincaré inequalities give a basic L p -norm control

More information

CONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence

CONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence 1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.

More information

Small Volume Fraction Limit of the Diblock Copolymer Problem: I. Sharp Interface Functional

Small Volume Fraction Limit of the Diblock Copolymer Problem: I. Sharp Interface Functional Small Volume Fraction Limit of the Diblock Copolymer Problem: I. Sharp Interface Functional Rustum Choksi Mark A. Peletier February, 2 Abstract We present the first of two articles on the small volume

More information

A CONVEX-CONCAVE ELLIPTIC PROBLEM WITH A PARAMETER ON THE BOUNDARY CONDITION

A CONVEX-CONCAVE ELLIPTIC PROBLEM WITH A PARAMETER ON THE BOUNDARY CONDITION A CONVEX-CONCAVE ELLIPTIC PROBLEM WITH A PARAMETER ON THE BOUNDARY CONDITION JORGE GARCÍA-MELIÁN, JULIO D. ROSSI AND JOSÉ C. SABINA DE LIS Abstract. In this paper we study existence and multiplicity of

More information

HOMEWORK 4 1. P45. # 1.

HOMEWORK 4 1. P45. # 1. HOMEWORK 4 SHUANGLIN SHAO P45 # Proof By the maximum principle, u(x, t x kt attains the maximum at the bottom or on the two sides When t, x kt x attains the maximum at x, ie, x When x, x kt kt attains

More information

Diffuse interface methods on graphs: Data clustering and Gamma-limits

Diffuse interface methods on graphs: Data clustering and Gamma-limits Diffuse interface methods on graphs: Data clustering and Gamma-limits Yves van Gennip joint work with Andrea Bertozzi, Jeff Brantingham, Blake Hunter Department of Mathematics, UCLA Research made possible

More information

Coupled second order singular perturbations for phase transitions

Coupled second order singular perturbations for phase transitions Coupled second order singular perturbations for phase transitions CMU 06/09/11 Ana Cristina Barroso, Margarida Baía, Milena Chermisi, JM Introduction Let Ω R d with Lipschitz boundary ( container ) and

More information

Diffuse interface models of tumor growth: optimal control and other issues

Diffuse interface models of tumor growth: optimal control and other issues Diffuse interface models of tumor growth: optimal control and other issues E. Rocca Università degli Studi di Pavia 4-6 décembre 2017 Laboratoire Jacques-Louis Lions First meeting of the French-German-Italian

More information

Interface and contact problems

Interface and contact problems Interface and contact problems Heiko Gimperlein Heriot Watt University and Maxwell Institute, Edinburgh and Universität Paderborn TU Wien May 30, 2016 H. Gimperlein (Edinburgh) Interface and contact problems

More information

Partial Differential Equations

Partial 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 information

Dirichlet s principle and well posedness of steady state solutions in peridynamics

Dirichlet s principle and well posedness of steady state solutions in peridynamics Dirichlet s principle and well posedness of steady state solutions in peridynamics Petronela Radu Work supported by NSF - DMS award 0908435 January 19, 2011 The steady state peridynamic model Consider

More information

Conservation law equations : problem set

Conservation law equations : problem set Conservation law equations : problem set Luis Silvestre For Isaac Neal and Elia Portnoy in the 2018 summer bootcamp 1 Method of characteristics For the problems in this section, assume that the solutions

More information

Stability Analysis of Stationary Solutions for the Cahn Hilliard Equation

Stability Analysis of Stationary Solutions for the Cahn Hilliard Equation Stability Analysis of Stationary Solutions for the Cahn Hilliard Equation Peter Howard, Texas A&M University University of Louisville, Oct. 19, 2007 References d = 1: Commun. Math. Phys. 269 (2007) 765

More information

Minimal Surface equations non-solvability strongly convex functional further regularity Consider minimal surface equation.

Minimal 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 information

A C 1 virtual element method for the Cahn-Hilliard equation with polygonal meshes. Marco Verani

A C 1 virtual element method for the Cahn-Hilliard equation with polygonal meshes. Marco Verani A C 1 virtual element method for the Cahn-Hilliard equation with polygonal meshes Marco Verani MOX, Department of Mathematics, Politecnico di Milano Joint work with: P. F. Antonietti (MOX Politecnico di

More information

ANALYSIS OF THE TV REGULARIZATION AND H 1 FIDELITY MODEL FOR DECOMPOSING AN IMAGE INTO CARTOON PLUS TEXTURE. C.M. Elliott and S.A.

ANALYSIS OF THE TV REGULARIZATION AND H 1 FIDELITY MODEL FOR DECOMPOSING AN IMAGE INTO CARTOON PLUS TEXTURE. C.M. Elliott and S.A. COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 6, Number 4, December 27 pp. 917 936 ANALYSIS OF THE TV REGULARIZATION AND H 1 FIDELITY MODEL FOR DECOMPOSING AN IMAGE

More information

Partial Differential Equations, 2nd Edition, L.C.Evans Chapter 5 Sobolev Spaces

Partial Differential Equations, 2nd Edition, L.C.Evans Chapter 5 Sobolev Spaces Partial Differential Equations, nd Edition, L.C.Evans Chapter 5 Sobolev Spaces Shih-Hsin Chen, Yung-Hsiang Huang 7.8.3 Abstract In these exercises always denote an open set of with smooth boundary. As

More information

Math 5520 Homework 2 Solutions

Math 5520 Homework 2 Solutions Math 552 Homework 2 Solutions March, 26. Consider the function fx) = 2x ) 3 if x, 3x ) 2 if < x 2. Determine for which k there holds f H k, 2). Find D α f for α k. Solution. We show that k = 2. The formulas

More information

NAVIER-STOKES EQUATIONS IN THIN 3D DOMAINS WITH NAVIER BOUNDARY CONDITIONS

NAVIER-STOKES EQUATIONS IN THIN 3D DOMAINS WITH NAVIER BOUNDARY CONDITIONS NAVIER-STOKES EQUATIONS IN THIN 3D DOMAINS WITH NAVIER BOUNDARY CONDITIONS DRAGOŞ IFTIMIE, GENEVIÈVE RAUGEL, AND GEORGE R. SELL Abstract. We consider the Navier-Stokes equations on a thin domain of the

More information

CAHN-HILLIARD-NAVIER-STOKES SYSTEMS WITH MOVING CONTACT LINES

CAHN-HILLIARD-NAVIER-STOKES SYSTEMS WITH MOVING CONTACT LINES CAHN-HILLIARD-NAVIER-STOKES SYSTEMS WITH MOVING CONTACT LINES C.G. Gal M Grasselli A Miranville To cite this version: C.G. Gal M Grasselli A Miranville. CAHN-HILLIARD-NAVIER-STOKES SYSTEMS WITH MOVING

More information

PREPUBLICACIONES DEL DEPARTAMENTO DE MATEMÁTICA APLICADA UNIVERSIDAD COMPLUTENSE DE MADRID MA-UCM

PREPUBLICACIONES DEL DEPARTAMENTO DE MATEMÁTICA APLICADA UNIVERSIDAD COMPLUTENSE DE MADRID MA-UCM PREPUBLICACIONES DEL DEPARTAMENTO DE MATEMÁTICA APLICADA UNIVERSIDAD COMPLUTENSE DE MADRID MA-UCM 2009-13 Extremal equilibria for monotone semigroups in ordered spaces with application to evolutionary

More information

Allen Cahn Equation in Two Spatial Dimension

Allen Cahn Equation in Two Spatial Dimension Allen Cahn Equation in Two Spatial Dimension Yoichiro Mori April 25, 216 Consider the Allen Cahn equation in two spatial dimension: ɛ u = ɛ2 u + fu) 1) where ɛ > is a small parameter and fu) is of cubic

More information

u( x) = g( y) ds y ( 1 ) U solves u = 0 in U; u = 0 on U. ( 3)

u( x) = g( y) ds y ( 1 ) U solves u = 0 in U; u = 0 on U. ( 3) M ath 5 2 7 Fall 2 0 0 9 L ecture 4 ( S ep. 6, 2 0 0 9 ) Properties and Estimates of Laplace s and Poisson s Equations In our last lecture we derived the formulas for the solutions of Poisson s equation

More information

1 Cheeger differentiation

1 Cheeger differentiation 1 Cheeger differentiation after J. Cheeger [1] and S. Keith [3] A summary written by John Mackay Abstract We construct a measurable differentiable structure on any metric measure space that is doubling

More information

Fact Sheet Functional Analysis

Fact Sheet Functional Analysis Fact Sheet Functional Analysis Literature: Hackbusch, W.: Theorie und Numerik elliptischer Differentialgleichungen. Teubner, 986. Knabner, P., Angermann, L.: Numerik partieller Differentialgleichungen.

More information

We denote the space of distributions on Ω by D ( Ω) 2.

We denote the space of distributions on Ω by D ( Ω) 2. Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study

More information

ANALYSIS OF A SCALAR PERIDYNAMIC MODEL WITH A SIGN CHANGING KERNEL. Tadele Mengesha. Qiang Du. (Communicated by the associate editor name)

ANALYSIS OF A SCALAR PERIDYNAMIC MODEL WITH A SIGN CHANGING KERNEL. Tadele Mengesha. Qiang Du. (Communicated by the associate editor name) Manuscript submitted to AIMS Journals Volume X, Number 0X, XX 200X doi:10.3934/xx.xx.xx.xx pp. X XX ANALYSIS OF A SCALAR PERIDYNAMIC MODEL WITH A SIGN CHANGING KERNEL Tadele Mengesha Department of Mathematics

More information

Lectures on. Sobolev Spaces. S. Kesavan The Institute of Mathematical Sciences, Chennai.

Lectures on. Sobolev Spaces. S. Kesavan The Institute of Mathematical Sciences, Chennai. Lectures on Sobolev Spaces S. Kesavan The Institute of Mathematical Sciences, Chennai. e-mail: kesh@imsc.res.in 2 1 Distributions In this section we will, very briefly, recall concepts from the theory

More information

Stationary solutions of curvature driven phase-field flow

Stationary solutions of curvature driven phase-field flow A.J. Vromans Stationary solutions of curvature driven phase-field flow Results for the Functionalized Cahn-Hilliard energy Master thesis, defended on August 12, 2013 Thesis advisor: prof.dr. A. Doelman

More information

Lecture No 2 Degenerate Diffusion Free boundary problems

Lecture No 2 Degenerate Diffusion Free boundary problems Lecture No 2 Degenerate Diffusion Free boundary problems Columbia University IAS summer program June, 2009 Outline We will discuss non-linear parabolic equations of slow diffusion. Our model is the porous

More information

Γ -convergence of the Allen Cahn energy with an oscillating forcing term

Γ -convergence of the Allen Cahn energy with an oscillating forcing term Interfaces and Free Boundaries 8 (2006), 47 78 Γ -convergence of the Allen Cahn energy with an oscillating forcing term N. DIRR Max-Planck Institute for Mathematics in the Sciences, Inselstr. 22, D-04103

More information

Time periodic solutions of Cahn Hilliard system with dynamic boundary conditions

Time periodic solutions of Cahn Hilliard system with dynamic boundary conditions Time periodic solutions of Cahn Hilliard system with dynamic boundary conditions arxiv:7.3697v [math.ap] Dec 7 Taishi Motoda Graduate School of Education, Kyoto University of Education Fujinomori, Fukakusa,

More information

COMPLETION OF A METRIC SPACE

COMPLETION OF A METRIC SPACE COMPLETION OF A METRIC SPACE HOW ANY INCOMPLETE METRIC SPACE CAN BE COMPLETED REBECCA AND TRACE Given any incomplete metric space (X,d), ( X, d X) a completion, with (X,d) ( X, d X) where X complete, and

More information

On the Cahn-Hilliard equation with irregular potentials and dynamic boundary conditions

On the Cahn-Hilliard equation with irregular potentials and dynamic boundary conditions On the Cahn-Hilliard equation with irregular potentials and dynamic boundary conditions Gianni Gilardi 1) e-mail: gianni.gilardi@unipv.it Alain Miranville 2) e-mail: Alain.Miranville@mathlabo.univ-poitiers.fr

More information

Convergence for periodic Fourier series

Convergence for periodic Fourier series Chapter 8 Convergence for periodic Fourier series We are now in a position to address the Fourier series hypothesis that functions can realized as the infinite sum of trigonometric functions discussed

More information

Switching, sparse and averaged control

Switching, 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 information

Adaptive evolution : a population approach Benoît Perthame

Adaptive 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 information

Solution Sheet 3. Solution Consider. with the metric. We also define a subset. and thus for any x, y X 0

Solution 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 information

A Mean Field Equation as Limit of Nonlinear Diffusions with Fractional Laplacian Operators

A Mean Field Equation as Limit of Nonlinear Diffusions with Fractional Laplacian Operators A Mean Field Equation as Limit of Nonlinear Diffusions with Fractional Laplacian Operators Sylvia Serfaty and Juan Luis Vázquez June 2012 Abstract In the limit of a nonlinear diffusion model involving

More information

Sobolev Spaces. Chapter Hölder spaces

Sobolev Spaces. Chapter Hölder spaces Chapter 2 Sobolev Spaces Sobolev spaces turn out often to be the proper setting in which to apply ideas of functional analysis to get information concerning partial differential equations. Here, we collect

More information

Bulletin of the. Iranian Mathematical Society

Bulletin of the. Iranian Mathematical Society ISSN: 1017-060X (Print) ISSN: 1735-8515 (Online) Bulletin of the Iranian Mathematical Society Vol. 42 (2016), No. 1, pp. 129 141. Title: On nonlocal elliptic system of p-kirchhoff-type in Author(s): L.

More information

On Asymptotic Variational Wave Equations

On Asymptotic Variational Wave Equations On Asymptotic Variational Wave Equations Alberto Bressan 1, Ping Zhang 2, and Yuxi Zheng 1 1 Department of Mathematics, Penn State University, PA 1682. E-mail: bressan@math.psu.edu; yzheng@math.psu.edu

More information

CAPACITIES ON METRIC SPACES

CAPACITIES ON METRIC SPACES [June 8, 2001] CAPACITIES ON METRIC SPACES V. GOL DSHTEIN AND M. TROYANOV Abstract. We discuss the potential theory related to the variational capacity and the Sobolev capacity on metric measure spaces.

More information

Boundary measures, generalized Gauss Green formulas, and mean value property in metric measure spaces

Boundary measures, generalized Gauss Green formulas, and mean value property in metric measure spaces Boundary measures, generalized Gauss Green formulas, and mean value property in metric measure spaces Niko Marola, Michele Miranda Jr, and Nageswari Shanmugalingam Contents 1 Introduction 2 2 Preliminaries

More information

OPTIMAL DISTRIBUTION OF OPPOSITELY CHARGED PHASES: PERFECT SCREENING AND OTHER PROPERTIES

OPTIMAL DISTRIBUTION OF OPPOSITELY CHARGED PHASES: PERFECT SCREENING AND OTHER PROPERTIES OPTIMAL DISTRIBUTION OF OPPOSITELY CHARGED PHASES: PERFECT SCREENING AND OTHER PROPERTIES MARCO BONACINI, HANS KNÜPFER, AND MATTHIAS RÖGER Abstract. We study the minimum energy configuration of a uniform

More information

TRAVELLING WAVES FOR A NONLOCAL DOUBLE OBSTACLE PROBLEM

TRAVELLING WAVES FOR A NONLOCAL DOUBLE OBSTACLE PROBLEM TRAVELLING WAVES FOR A NONLOCAL DOUBLE OBSTACLE PROBLEM PAUL C. FIFE, Mathematics Department, University of Utah, Salt Lake City, Utah 84112, USA Abstract Existence, uniqueness, and regularity properties

More information

GIOVANNI COMI AND MONICA TORRES

GIOVANNI COMI AND MONICA TORRES ONE-SIDED APPROXIMATION OF SETS OF FINITE PERIMETER GIOVANNI COMI AND MONICA TORRES Abstract. In this note we present a new proof of a one-sided approximation of sets of finite perimeter introduced in

More information

Convergence to equilibrium for rough differential equations

Convergence to equilibrium for rough differential equations Convergence to equilibrium for rough differential equations Samy Tindel Purdue University Barcelona GSE Summer Forum 2017 Joint work with Aurélien Deya (Nancy) and Fabien Panloup (Angers) Samy T. (Purdue)

More information

MATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Fall 2011 Professor: Jared Speck

MATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Fall 2011 Professor: Jared Speck MATH 8.52 COURSE NOTES - CLASS MEETING # 6 8.52 Introduction to PDEs, Fall 20 Professor: Jared Speck Class Meeting # 6: Laplace s and Poisson s Equations We will now study the Laplace and Poisson equations

More information

MATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Spring 2018 Professor: Jared Speck

MATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Spring 2018 Professor: Jared Speck MATH 8.52 COURSE NOTES - CLASS MEETING # 6 8.52 Introduction to PDEs, Spring 208 Professor: Jared Speck Class Meeting # 6: Laplace s and Poisson s Equations We will now study the Laplace and Poisson equations

More information