Variational and Topological methods : Theory, Applications, Numerical Simulations, and Open Problems 6-9 June 2012, Northern Arizona University

Transcription

1 Variational and Topological methods : Theory, Applications, Numerical Simulations, and Open Problems 6-9 June 22, Northern Arizona University Some methods using monotonicity for solving quasilinear parabolic equations (an introductory lecture for graduate students and postdocs) Jacques Giacomoni LMAP (UMR 542), Université de Pau et des pays de l Adour Avenue de l université, 643 Pau Cedex. (jgiacomo@univ-pau.fr) /3

2 We are first concerned by the following problem : (P t ) u t p u = f (x, u) in Q T def = (, T ) u = on Σ T = [, T ], u(, ) = u ( ) in i) is a bounded domain in IR d with smooth boundary, ii) < p <, p u def = ( u p 2 u), T >, u L () ( W,p ()). iii) f : (x, s) IR f (x, s) Caratheodory function, locally Lipschitz with respect to s uniformly in x and satisfying the growth condition : f (x, s) a s q + b a.a. x () with q >, a > and b. 2/3

3 Issues :. local existence of weak solutions, 2. regularity of weak solutions. Such problems arise in different models : ) non Newtonian flows (pseudo-plastic fluids : < p < 2, dilatant fluids : p > 2), 2) Models of flow through porous media in turbulent regimes (ex. flow through rock filled dams [Díaz-De Thelin]), 3) models in glaceology, reaction diffusion problems, petroleum extraction, climate models, etc. Bibliography :. I.J. Díaz, Nonlinear partial differential equations and free boundaries, Elliptic Equations, Research Notes in Mathematics, Vol. I, 6, Pitman London, R. Aris, The Mathematical Theory of Diffusion and Reaction in Permeable Catalysts, Volume I and II. Clarenton Press, Oxford, /3

4 Auxilary problem : def u t p u = h(x, t) in Q T = (, T ) (S t ) u = on Σ T = [, T ], u(, ) = u ( ) in with T >, q min(2, p ), u L q (), h L (, T ; L q ()). Properties of p-laplace operator ( < p < ) :. W,p () u p u W,p () isomorphim ;. monotonicity : u, v W,p () < p u + p v, u v > ; aplications : weak comparison principle, Minty-Browder theory.. variational : p u = J (u) with J(v) def = p v p dx, v W,p () ; J convex, s.c.i. p maximal monotone in L 2 (). 4/3

5 Maximal monotone operators, m-accretive operators Let H a Hilbert space and A : D(A) H P(H) (multi-valued) Definition A is maximal monotone in H if. A is monotone : u, v D(A), ψ Au, η Av, (ψ η u v) (( ) denotes the scalar product). 2. A is maximal : R(I + A) = H. Remarks :. A is maximal monotone J λ def = (I + λa) is a contraction in H for any λ >. 2. A single valued, linear A maximal monotone D(A) dense in H. Example : Let V a reflexive Banach space such that V H V (densely and continuously) and A : V V single-valued, hemicontinuous, coercive and monotone. Then A H is maximal monotone (by Minty-Browder Theorem). 5/3

6 Let X a banach space. A X X. Definition A is m-accretive in X if. (accretive) [x, y ], [x 2, y 2 ] A and λ >, x x 2 X x x 2 + λ(y y 2 ) X. 2. (maximal) R(I + A) = X. Example : p is m-accretive in L q Np () for q min(2, Np N+p ) (accretive for q ). Indeed, introduce for λ > and f L q () the energy functional E defined by E(u) def = u 2 dx + λ u p dx fu dx. 2 p E is convex and s.c.i. in L q () W,p () for q min(2, Np Np N+p ). 6/3

7 Remark : By the weak comparison principle, p is m-accretive in L (). Example 2 : Let β a maximal monotone operator in IR, q. Let β the realization of β in L q (). Then, A defined by Au def = u + β(u), is m-accretive in L q (). Other examples : porous media equation involves (β(u)) [with β() and β maximal monotone in IR] which is m-accretive in L (). Theory of maximal monotone or m-accretive operators existence of weak solutions to du dt + Au f, u() = u. (2) References : H. Brezis Opérateurs maximaux monotones et semi-groupes de contractions dans des espaces de Hilbert, V. Barbu, Non-Linear Differential Equations of Monotone Type On Banach Space. 7/3

8 Definition Let ɛ >, T >, N IN\{} and f L (, T ; X ).. An ɛ-discretization on [, T ] (denoted by D ɛ (t,, t N ; f,, f N )) of equation in (2) consists of a partition = t t t 2 t N = T and a finite sequence {f i } N i= such that t i t i < ɛ for i =,, N and N i= ti t i f (s) f i ds < ɛ. 2. An ɛ-approximate solution to (2) associated to D ɛ (t,, t N ; f,, f N ) is a piecewise constant function z : [, T ] X whose (nodal) values z i satisfies on (t i, t i ] satisfy z i z i t i t i + Az i f i, i =,, N and z() u ɛ. 8/3

9 Definition Let A be m-accretive in X.. A mild solution to the Cauchy problem (2) is a function y C([, T ]; X ) such that for any ɛ >, there is ɛ-approximate solution z on [, T ] such that sup t [,T ] y(t) z(t) ɛ and y() = u. 2. A strong solution to the Cauchy problem (2) is a function y W, ((, T ]; X ) C([, T ]; X ) such that f (t) dy dt (t) Ay(t), a.e. t (, T ), y() = u. Remark. A strong solution to (2) is unique and continuous in respect to f and u. 9/3

10 Theorem. (existence of mild solutions) Let A be m-accretive in X and u D(A) X. Then, there exists a unique mild solution to (2). Proof. (main steps) Let r, k L (, T ; X ), ɛ >, η > and u ɛ, v η be ɛ-approximate solution and η-approximate solution respectively associated to the discretisations D ɛ (iɛ, r i def D η (iη, k i def We define = η = ɛ iɛ (i )ɛ r(s) ds) and iη (i )η k(s) ds) with initial data u and v in D(A) X. ϕ(t, s) def = r(t) k(s) X (t, s) [, T ] [, T ] and for a fixed z D(A) b(t, r, k) def = u z X + v z X + t Az X + t + r(τ) X dτ + t k(τ) X dτ, t [ T, T ], and /3

11 Ψ(t, s) def = b(t s, r, k)+ s t the solution to the transport equation ϕ(t s + τ, τ)dτ if s t T, ϕ(τ, s t + τ)dτ if t s T, Ψ Ψ (t, s) + (t, s) = ϕ(t, s) (t, s) [, T ] [, T ], t s Ψ(t, ) = b(t, r, k) t [, T ], Ψ(, s) = b( s, r, k) s [, T ]. (3) Finally, Let {uɛ n } and {vη n } be the nodal functions associated to u ɛ and v η respectively. For n, m IN, u n ɛ v m η + ɛη ɛ + η (Aun ɛ Av m η ) = + ɛ ɛ + η (un ɛ v m η ) + η ɛ + η (un ɛ vη m ) ɛη ɛ + η (r n k m ), /3

12 Accretivity of A in X implies that Φ ɛ,η n,m = u n ɛ v m η X satisfies η ɛ + η Φɛ,η n,m + ɛ ɛ + η Φɛ,η n,m + Φ ɛ,η n, b(t n, r ɛ, k η ) and Φ ɛ,η,m b( s m, r ɛ, k η ), Φ ɛ,η n,m ɛη ɛ + η r n k m X, For the two last inequalities, note that from accretivity of A, we have : uɛ n z X uɛ n z + ɛ(r n + ɛ ( un ɛ + uɛ n ) Az) X u n ɛ z X + ɛ r n X + ɛ Az X, vη m z X vη m z + η(r m + η ( v η m + vη m ) Az) X v m η z X + η k m X + η Az X, 2/3

13 Clearly, Φ ɛ,η n,m Ψ ɛ,η n,m where Ψ ɛ,η n,m satisfies Ψ ɛ,η η n,m = ɛ + η Ψɛ,η n,m + ɛ ɛ + η Ψɛ,η n,m + ɛη ɛ + η r n k m X, Ψ ɛ,η n, = b(t n, r ɛ, k η ) and Ψ ɛ,η,m = b( s m, r ɛ, k η ). We have b(, r ɛ, k η ) b(, r, k) in L ( T, T ; X ) and Φ ɛ,η ϕ in L ((, T ) (, T ); X ). Then, ρ ɛ,η = Ψ ɛ,η Ψ L ([,T ] [,T ]) as (ɛ, η). We get for t [, T ] and s [, T ] u ɛ (t) v η (s) X = Φ ɛ,η (t, s) Ψ ɛ,η (t, s) Ψ(t, s) + ρ ɛ,η, (4) 3/3

14 Doing in (4) t = s, r = k, v = u : We have also u ɛ (t) u η (t) X 2 u z X + ρ ɛ,η. u ɛ (t) ũ ɛ (t) X 2 u ɛ (t) u ɛ (t ɛ) X 4 u z X + 2ɛ Az X + 2 t + 2 h(ɛ + τ) h(τ) X dτ. ɛ h(τ) X dτ From the fact that u D(A) X, {u ɛ } ɛ>, {ũ ɛ } ɛ> are Cauchy sequences in L (, T ; X ) which converge to u C([, T ]; X ) which satisfies u() = u and (by (4) with v = u ) u(t) u(s) X 2 u z X + t s Az + + max(t,s) t s h(τ) X dτ h( t s + τ) h(τ) X dτ. (5) 4/3

15 Passing to the limit in (4) with t = s we obtain t u(t) v(t) X u z X + v z X + r(τ) k(τ) X dτ, Since v D(A) X, we get t u(t) v(t) X u v X + r(s) k(s) X ds, t T. (6) which implies that u is continuous in respect to u and f. If u D(A) and h W, (, T ; X ), then u is absolutely continuous. Indeed (for s t), u(t) u(s) X u u(t s) X + h(τ) h(τ + t s) X dτ ( T ) (t s) ( Au h() X + dh(τ) dt dτ. X s 5/3

16 Remark. A generates a continuous semi-group of contractions on D(A) X, that is a family of mappings that maps D(A) X into itself with the properties : (i) S(t + s)x = S(t)S(s)x, x D(A) X, t, s. (ii) S()x = x, x D(A) X. (iii) For every x D(A) X, the function t S(t)x is continuous on [, ) and solves (2) with f =. (iv) S(t)x S(t)y X x y X, t, x, y D(A) X. Aplications : p is the infinitesimal generator of a continuous semi-group of contractions in L q () with q [2, ) and in C () (q = ) with { } D(A) = u L q () W,p () p u L q (), q [2, ]. 6/3

17 Concerning (P t ), (local) existence of weak solutions can be proved as follows :. Existence (by Schauder fixed-point theorem) and uniqueness of global solutions u k L (, T ; W,p () L (Q) such that L 2 (Q) to u k t (Pk t ) where u t p u = f k (x, u) in Q def = (, T ) u = on Σ T = [, T ], u > in Q u(, ) = u ( ) in f k (x, t) def = f (x, T k (t)) with T k (t) def = { sign(t)k if t k t if k t k 2. Existence of barrier functions or a priori estimates u k weak solution to (P t ) for T > small enough. 3. compact semi-group, maximal time interval [, T ) : blow up in finite time if T < (see [Vrabie]). 7/3

18 We can adapt the method of semi-discretization in time on the following quasilinear Barenblatt equation : { f (, t u) (P) ~p u = g in Q =], T [, u = on Σ = (, T ) Γ, u(x, ) = u (x) in where ~p u def = ) d i= x i ( xi u p i (x) 2 xi u is the anisotropic ~p(x)-laplace operator such that ~p = (p i ) i=,..,d : [p, p + ] 2d with max(, d+2 ) < p p + <, continuous in with logarithmic module of continuity. f satisfies : (f) f : (x, t) IR IR is a Carathéodory function, non-decreasing with respect to t and (f2) There exist a L 2 () and a 2 such that x a.e., t IR, f (x, t) a (x) + a 2 t, (7) (f3) there exist b L () and b 2 such that, x a.e., t IR, tf (x, t) b 2 t 2 b (x). (8) 8/3

19 Concerning g and the initial data u, we require (g) g L 2 (Q) and u W,~p(x) (). Basic properties of the spaces L p(x) () and W,~p(x) () Let h a measurable function. Set ρ p (h) def = h(x) p(x) dx and the luxembourg norm f p def = inf { ( ) } f λ >, ρ p. λ We define L p(x) () = W,~p(x) () = { } u : IR, measurable/ x u(x) p(x) L (), { } u W, () / xi u L p i (x) (), i =,.., d. 9/3

20 . (L p(x) (), f p ) is a separable uniform convex space. 2. Dual space : L p (x) () where p (x) = p(x) p(x) and Hölder inequality : f (x)g(x) dx c f p g p. (9) 3. Let (u n ) n L p(x) (). Then, lim u n = in L p(x) () if and n only if lim n ρ p (u n ) =. If u n u in L p(x) (), then u(x) p(x) dx lim inf u n (x) p(x) dx. n 4. L p(x) () L q(x) () if and only if q(x) p(x). 5. Endowed with the graph norm u,~p(x) W () = d i= x i u pi, W,~p(x) () is a reflexive and separable Banach space. 2/3

21 Morawetz Inequality Lemma Let < p <. Set s def = min(, p 2 ) and C(p) def = ( ) p 2 2 p 2 if p 2. if p < 2 Then, a, b IR, a b p ( 2 C(p) ( a p + b p ) s a p + b p 2 a + b s 2 p). () Compactness If p i (x) q i (x) for any i, then W,~q(x) () W,~p(x) (). Then, if p > 2d d+2, W,~p(x) () L 2 (). 2/3

22 Definition A solution of (P) is any function u L (, T ; W,~p(x) ()) with t u L 2 (Q) such that u(,.) = u and, for any v W,~p(x) (), a.e. in ], T [, f (., t u)vdx + a ~p (u, v) = gv dx. We show : Theorem (G, Vallet) Assume that u W,~p(x) () and the conditions (f)-(f3) and (g). Then, there exists a solution u to (P). Furthermore, u C([, T [, W,~p(x) ()). 22/3

23 Semi-discretization in time : N IN\{}, t def = T N, { f (., un u n t ) ~p u n = g n in u n = with g n def = t n t (n ) t g(s)d s on [(n ) t, n t ). u = u and u n W,~p(x) () is the unique global minimizer to E defined by E(u) = t F (, u un t ) dx + where F (., λ) = λ f (., σ)dσ. d i= p i (x) x i u p i (x) dx g n u dx 23/3

24 We define for t ((n ) t, n t ] Note that Thus, g t (t) = g n, u t (t) = u n, ũ t (t) = t (n ) t t (u n u n ) + u n. g t 2 L 2 (Q) = t N g n 2 L 2 () g 2 L 2 (Q), n= f (, t ũ t ) ~p u t = g t Energy estimates t ũ t L 2 (Q) C, in Q T =], T [. ũ t, u t are bounded in L (, T ; W,~p(x) ()), ũ t u t L 2 (Q) t t ũ t L 2 (Q) C t. 24/3

25 up to a subsequence, ũ t u in L (, T ; W,~p(x) ()), u t v in L (, T ; W,~p(x) ()), ũ t, u t u = v in L (, T ; L 2 ()) and t ũ t t u in L 2 (Q). Note that u L (, T ; W,~p(x) ()) H (, T, L 2 ()) and then u C w ([, T ], W,~p(x) ()) and u(t) W,~p(x) (), for any t [, T ]. Strong convergence by Theorem (Aubin-Simon) Consider p ], + [, q [, + ] and V, E and F three Banach spaces such that V E F. Then, if A is a bounded subset of W,p (, T, F ) and of L q (, T, V ), A is relatively compact in C([, T ], F ) and in L q (, T, E). 25/3

26 u t u in L p (, T ; W,~p(x) ()) for any p < +. Consequently, setting ~q(x) = (q i (x)) with q i (x) def = p i (x) p i (x) ~p u t ~p u in L p p (, T ; W,~q(x) ()) Next, there exists χ L 2 (Q) such that as t + (up to a subsequence) one has for any < p < +. f ( t ũ t ) χ in L 2 (, T ; L 2 ()) and χ ~p u = g in D (Q). From g t g in L 2 (Q) as t + and convexity arguments, we get s lim sup f ( t ũ t ) t ũ t dxdt + t + d s p i (x) x i u p i (x) dx + i= d i= g t u dxdt. p i (x) x i u(s) p i (x) dx 26/3

27 we have for < t << that s t = s t χ(t) u(t + t) u(t) + t g(t) u(t + t) u(t) t s t ( a ~p u(t), u(t + ) t) u(t) t where and as t + s a ~p : (u, v) lim inf t d i= χ(t) t u + lim inf t d t i= t xi u p i (x) 2 xi u xi vdx. d t i= s s t p i (x) x i u(t) p i (x) + p i (x) x i u(t) p i (x) s g(t) t u. 27/3

28 Since u C w ([, T ], W,~p(x) ()), one has s d i= lim sup t + χ(t) t u + lim inf t s d p i (x) x i u p i (x) dx + s t i= s t s f ( t ũ t ) t ũ t + d i= g(t) t u p i (x) x i u(t) p i (x) p i (x) x i u(s) p i (x). Therefore, for any Lebesgue point s of u in L p+ (, T, W,~p(x) ()), we get that s lim sup t + χ(t) t u + s d i= f ( t ũ t ) t ũ t + p i (x) x i u(s) p i (x) d i= p i (x) x i u(s) p i (x). 28/3

29 Thus, lim sup t + s f ( t ũ t ) t ũ t dxdt s χ t u dxdt. From the monotonicity of f, we know that for any v L 2 (Q s ) def where Q s = ], s[, one has [f (v) f ( t ũ t )](v t ũ t ) dxdt. () Q s Taking the lim sup as t + in (), we obtain that [f (v) χ](v t u) dxdt. Q s A classical monotonicity argument (Minty s trick!) yields for s a.e. in ], T [, χ = f ( t u) in L 2 (Q s ), that is, χ = f ( t u) in L 2 (Q) and then f ( t u) ~p u = g. 29/3

30 Lemma Assume that hypothesis (f)-(f3) are satisfied and let T >, Q =], T [, u W,~p(x) (), g L 2 (Q) and u H (Q) C w ([, T ], W,~p(x) ()) the solution of the ~p-laplace evolution equation : f (, t u) ~p u = g in Q, with u(,.) = u. Then, u C([, T [, W,~p(x) ()) and for any t [, T [, = ],t[ d i= f (, t u(σ))u(σ) dxdσ + p i (x) x i u p i (x) dx + d i= ],t[ p i (x) x i u(t, x) p i (x) g(σ) t u dxdσ. (2) 3/3