Stability of nonlinear locally damped partial differential equations: the continuous and discretized problems. Part I
|
|
- Lora Hopkins
- 5 years ago
- Views:
Transcription
1 . Stability of nonlinear locally damped partial differential equations: the continuous and discretized problems. Part I Fatiha Alabau-Boussouira 1 Emmanuel Trélat 2 1 Univ. de Lorraine, LMAM 2 Univ. Paris 6 (LJLL) and IUF 33th International Summer School of Automatic Control, Grenoble, september
2 Outline of the first lecture 1 General motivations for dissipative vibrating systems 2 Dissipative systems in finite dimensions 3 Infinite dimensional dissipative systems 4 The scalar case 5 The optimal-weight convexity method 6 Lower energy estimates: an energy comparison principle 7 Examples of decay rates and optimality 8 Vectorial case: parameter-dependent upper estimates
3 General motivations for dissipative vibrating systems Sommaire 1 General motivations for dissipative vibrating systems 2 Dissipative systems in finite dimensions 3 Infinite dimensional dissipative systems 4 The scalar case 5 The optimal-weight convexity method 6 Lower energy estimates: an energy comparison principle 7 Examples of decay rates and optimality 8 Vectorial case: parameter-dependent upper estimates
4 General motivations for dissipative vibrating systems Many physical phenomenon s, such as the propagation of waves in a string, a membrane or a plate are modelized by finite or infinite dimensional vibrating systems. In applications, the engineers want to reduce to zero the vibrations of the solutions by auto-regulation. This is performed through the implementation of appropriate feedbacks on selected regions of the physical device. The auto-regulated system is a damped system and the vibrations of its solutions are in general measured through their energies.
5 General motivations for dissipative vibrating systems We consider reversible dissipative systems in both finite and infinite dimensions in autonomous or non autonomous form, meaning that the damping effect can be nonlocal. The common feature for dissipative systems is : one can associate a (physical) energy to them moreover this energy is decaying as time increases. Damping (feedback) = the energy is decaying
6 General motivations for dissipative vibrating systems We consider reversible dissipative systems in both finite and infinite dimensions in autonomous or non autonomous form, meaning that the damping effect can be nonlocal. The common feature for dissipative systems is : one can associate a (physical) energy to them moreover this energy is decaying as time increases. Damping (feedback) = the energy is decaying
7 General motivations for dissipative vibrating systems We consider reversible dissipative systems in both finite and infinite dimensions in autonomous or non autonomous form, meaning that the damping effect can be nonlocal. The common feature for dissipative systems is : one can associate a (physical) energy to them moreover this energy is decaying as time increases. Damping (feedback) = the energy is decaying
8 General motivations for dissipative vibrating systems We consider reversible dissipative systems in both finite and infinite dimensions in autonomous or non autonomous form, meaning that the damping effect can be nonlocal. The common feature for dissipative systems is : one can associate a (physical) energy to them moreover this energy is decaying as time increases. Damping (feedback) = the energy is decaying
9 General motivations for dissipative vibrating systems We consider reversible dissipative systems in both finite and infinite dimensions in autonomous or non autonomous form, meaning that the damping effect can be nonlocal. The common feature for dissipative systems is : one can associate a (physical) energy to them moreover this energy is decaying as time increases. Damping (feedback) = the energy is decaying
10 General motivations for dissipative vibrating systems The energy measures the vibrations of the device, thus it is important to determine its asymptotic behavior at large time. Questions: Is the damping sufficient to induce strong stabilization (the energy goes to 0 at )? This is known to hold in general. References: Dafermos, Ball, Slemrod... in the 70 s, beggining of 80 s What about decay rates and how they depend on the feedback properties? Unified approach between feedbacks types? Optimality of the decay rates? Links/differences between the finite and infinite dimensions? Discretization of systems of infinite dimensional vibrating systems and links between the continuous and discretized systems?
11 General motivations for dissipative vibrating systems The energy measures the vibrations of the device, thus it is important to determine its asymptotic behavior at large time. Questions: Is the damping sufficient to induce strong stabilization (the energy goes to 0 at )? This is known to hold in general. References: Dafermos, Ball, Slemrod... in the 70 s, beggining of 80 s What about decay rates and how they depend on the feedback properties? Unified approach between feedbacks types? Optimality of the decay rates? Links/differences between the finite and infinite dimensions? Discretization of systems of infinite dimensional vibrating systems and links between the continuous and discretized systems?
12 Dissipative systems in finite dimensions Sommaire 1 General motivations for dissipative vibrating systems 2 Dissipative systems in finite dimensions 3 Infinite dimensional dissipative systems 4 The scalar case 5 The optimal-weight convexity method 6 Lower energy estimates: an energy comparison principle 7 Examples of decay rates and optimality 8 Vectorial case: parameter-dependent upper estimates
13 Dissipative systems in finite dimensions The scalar case : where Dissipative systems in finite dimensions u + ν u + f (u) + ρ(u ) = 0. ν > 0 u is a scalar unknown f is locally lipschitz continuous and sf (s) 0 for all s R ρ is monotone increasing, ρ(0) = 0. Multiply the equation by u. This gives the so-called dissipation relation: where ( 1 2( u (t) 2 + ν u(t) 2) + F(u(t))) = u (t)ρ(u (t))., F (u) = u 0 f (s) ds.
14 Dissipative systems in finite dimensions The scalar case : where Dissipative systems in finite dimensions u + ν u + f (u) + ρ(u ) = 0. ν > 0 u is a scalar unknown f is locally lipschitz continuous and sf (s) 0 for all s R ρ is monotone increasing, ρ(0) = 0. Multiply the equation by u. This gives the so-called dissipation relation: where ( 1 2( u (t) 2 + ν u(t) 2) + F(u(t))) = u (t)ρ(u (t))., F (u) = u 0 f (s) ds.
15 Dissipative systems in finite dimensions Define the energy E(t) = 1 2 ( u (t) 2 + ν u(t) 2) + F(u(t)), Kinetic energy Then the dissipation relation becomes potential energy E (t) = u (t)ρ(u (t)). Since sρ(s) 0 for all s R, we have E (t) = u (t)ρ(u (t)) 0 t 0. = the energy is a Lyapunov function, however it is not a strict Lyapunov function, since the dissipation relation does not involve the potential part of the energy.
16 Dissipative systems in finite dimensions The vectorial case: where the unknown u R n and u + Au + f (u) + Bρ(u ) = 0. A is a symmetric, positive definite matrix B = diag(b i ) 1 i n, with b i 0 i {1,..., n} (f (u)) i = f (u i ), (ρ(u )) i = ρ(u i ) Take the euclidian scalar product of the above equality with u, and denote by the euclidian norm. This gives n n u, u + Au, u + (F(u i )) = b i ρ(u i )u i. where F = f as for the scalar case. i=1 i=1
17 Dissipative systems in finite dimensions The vectorial case: where the unknown u R n and u + Au + f (u) + Bρ(u ) = 0. A is a symmetric, positive definite matrix B = diag(b i ) 1 i n, with b i 0 i {1,..., n} (f (u)) i = f (u i ), (ρ(u )) i = ρ(u i ) Take the euclidian scalar product of the above equality with u, and denote by the euclidian norm. This gives n n u, u + Au, u + (F(u i )) = b i ρ(u i )u i. where F = f as for the scalar case. i=1 i=1
18 Dissipative systems in finite dimensions Thus, for this example the energy is given by E(t) = 1 2 ( u 2 + Au, u + n ) (F(u i )), and we have the dissipation relation (under the same hypothesis on ρ than in the scalar case) i=1 n E (t) = b i ρ(u i )u i 0 t 0. i=1 Thus, the energy is again a Lyapunov function, but not a strict one.
19 Dissipative systems in finite dimensions The motivation for vectorial system comes from the semi-discretization of wave-type equations. Let us consider a frictional dissipative wave equation in the one-dimensional space domain (0, 1). t 2u 2 x u + f (u) + b(x)ρ(x, t u) = 0, 0 < t, 0 < x < 1, u(t, x) = 0, for x = 0, x = 1, 0 < t, u(0, x) = u 0 (x), t u(0, x) = u 1 (x), 0 < x < 1. (1) We assume that this system is dissipative, i.e. b 0 is the damping coefficient (locally supported in general), ρ is monotone nondecreasing with respect to the second variable and ρ(., 0) = 0. Additional hypotheses on f to guarantee existence for all t 0.
20 Dissipative systems in finite dimensions A semi-discretization of the above equation in space, with for instance a uniform mesh x i = i h for i = 0,..., n + 1 with a parameter of discretization h = 1/(n + 1) gives the finite dimensional system u i u i+1 + u i 1 2u i h 2 + f (u i ) + b i ρ i (u i ) = 0, 0 < t, i = 1,..., n, u 0 (t) = u n+1 (t) = 0, 0 < t, u i (0) = u i,0, t u i (0) = u i,1, i = 1,..., n, (2) where u i is a function of t which stands for an approximation of the solution u at point x i b i = b(x i ), ρ i (s) = ρ(x i, s) for all s R.
21 Dissipative systems in finite dimensions A semi-discretization of the above equation in space, with for instance a uniform mesh x i = i h for i = 0,..., n + 1 with a parameter of discretization h = 1/(n + 1) gives the finite dimensional system u i u i+1 + u i 1 2u i h 2 + f (u i ) + b i ρ i (u i ) = 0, 0 < t, i = 1,..., n, u 0 (t) = u n+1 (t) = 0, 0 < t, u i (0) = u i,0, t u i (0) = u i,1, i = 1,..., n, (2) where u i is a function of t which stands for an approximation of the solution u at point x i b i = b(x i ), ρ i (s) = ρ(x i, s) for all s R.
22 Dissipative systems in finite dimensions A semi-discretization of the above equation in space, with for instance a uniform mesh x i = i h for i = 0,..., n + 1 with a parameter of discretization h = 1/(n + 1) gives the finite dimensional system u i u i+1 + u i 1 2u i h 2 + f (u i ) + b i ρ i (u i ) = 0, 0 < t, i = 1,..., n, u 0 (t) = u n+1 (t) = 0, 0 < t, u i (0) = u i,0, t u i (0) = u i,1, i = 1,..., n, (2) where u i is a function of t which stands for an approximation of the solution u at point x i b i = b(x i ), ρ i (s) = ρ(x i, s) for all s R.
23 Dissipative systems in finite dimensions Set A = h which is symmetric and positive definite. We can rewrite the semi-discretized equation as the vectorial system where the unknown u R n u + Au + f (u) + ρ(u ) = 0. In a similar way, we can consider a semi-discretization of the plate (or Petrowsky) equation.
24 Dissipative systems in finite dimensions Thus, semi-discretized infinite dimensional vibrating systems enter in the framework of finite dimensional vibrating systems. We want optimal energy decay rates for the scalar and vectorial systems. As a consequence, this would also give optimal decay rates for semi-discretized systems. But in addition, it is important for applications to trace the dependence of the estimates on the discretization parameters and to obtain numerical schemes which lead to uniform decay rates with respect to h.
25 Dissipative systems in finite dimensions Thus, semi-discretized infinite dimensional vibrating systems enter in the framework of finite dimensional vibrating systems. We want optimal energy decay rates for the scalar and vectorial systems. As a consequence, this would also give optimal decay rates for semi-discretized systems. But in addition, it is important for applications to trace the dependence of the estimates on the discretization parameters and to obtain numerical schemes which lead to uniform decay rates with respect to h.
26 Dissipative systems in finite dimensions The discrete energy is defined by E h (t) = h 2 n { (u i ) 2 (t) + ( u i+1 u i ) } 2 h i=0 and satisfies, in the linear damped, i.e. when ρ(., s) = a(.)s, the discrete dissipation relation E h(t) = h n a i (u i ) 2 (t) i=1
27 Dissipative systems in finite dimensions It is well-known that in the linear damped case, the above numerical scheme leads to nonuniform exponential decay rates with respect to h. These results are strongly related to observability estimates for the undamped equation, which are not uniform with respect to h. These results are due to Glowinski-Li-Lions, Glowinski, Infante-Zuazua, Tebou-Zuazua... This phenomenon is due to high frequency numerical spurious oscillations. Thus in the nonlinear damped case, the first step is to derive optimal energy decay rates which are h-dependent then check how to modify the numerical scheme to derive optimal energy decay rates that are uniform with respect to h.
28 Infinite dimensional dissipative systems Sommaire 1 General motivations for dissipative vibrating systems 2 Dissipative systems in finite dimensions 3 Infinite dimensional dissipative systems 4 The scalar case 5 The optimal-weight convexity method 6 Lower energy estimates: an energy comparison principle 7 Examples of decay rates and optimality 8 Vectorial case: parameter-dependent upper estimates
29 Infinite dimensional dissipative systems Infinite dimensional systems We consider the following abstract equation u + Au + feedback operator[u] = 0 (u, u )(0) = (u 0, u 1 ) Here A stands for an unbounded linear operator in an Hilbert space H, which is closed, coercive self-adjoint with dense domain in H.
30 Infinite dimensional dissipative systems As for the finite dimensional case, one can associate an energy to the solutions and this energy decays through time. The same questions than for the finite dimensional case arise, except now that there is a space dependence and the problem is set-up in an infinite dimensional framework (Sobolev spaces) Moreover the feedback can be locally or boundary supported and in general the support of the feedback coefficient has to satisfy geometric conditions. We shall come back later on infinite dimensional systems.
31 Infinite dimensional dissipative systems As for the finite dimensional case, one can associate an energy to the solutions and this energy decays through time. The same questions than for the finite dimensional case arise, except now that there is a space dependence and the problem is set-up in an infinite dimensional framework (Sobolev spaces) Moreover the feedback can be locally or boundary supported and in general the support of the feedback coefficient has to satisfy geometric conditions. We shall come back later on infinite dimensional systems.
32 The scalar case Sommaire 1 General motivations for dissipative vibrating systems 2 Dissipative systems in finite dimensions 3 Infinite dimensional dissipative systems 4 The scalar case 5 The optimal-weight convexity method 6 Lower energy estimates: an energy comparison principle 7 Examples of decay rates and optimality 8 Vectorial case: parameter-dependent upper estimates
33 The scalar case Let us start by the scalar case example: u + ν u + f (u) + ρ(u ) = 0. where ν > 0 u is a scalar unknown f is locally lipschitz continuous, sf (s) 0 for all s R µ > 0 such that 0 F (s) µsf (s), s R, with F(u) = u ρ is monotone increasing, ρ(0) = 0 0 f (s) ds
34 The scalar case We assume that the feedback satisfies the assumption ρ C(R), is monotone increasing ρ(0) = 0 a strictly increasing odd function g such that (A1) c s ρ(s) C s, s 1 c g( s ) ρ(s) C g 1 ( s ), s 1 r 0 > 0 such that g C 1 ([0, r 0 ]), g(0) = g (0) = 0 where g 1 denotes the inverse function of g and where c, C are positive constants. Remark The function g has an arbitrary sublinear growth the interesting case close to 0. It captures the behavior of the feedback close to 0.
35 The scalar case Remark If g (0) 0 then it is as for the linear feedback case. So it will not be considered here. The assumption of linear growth at infinity is made for the sake of simplicity. It can be removed using a well-known strong stabilization result (based on Lasalle invariance principle) as follows E(t) 0 as t. Thus for sufficiently large time u (t) r 0.
36 The scalar case The energy is defined as E(t) = 1 2 ( u (t) 2 + ν u(t) 2) + F(u(t)), where F (u) = u 0 f (s) ds. Recall the dissipation relation: E (t) = u (t)ρ(u (t)). where E is the energy of the solution.
37 The scalar case We suppose f 0, ν = 1. Then E is given by E(t) = 1 2 ( u (t) 2 + ν u(t) 2). It is well-known that when ρ(s) = s p 1 s with p > 1, E(t) decays as t 2/(p 1) and this decay rate is optimal: the result is due to Haraux Idea of the proof: E is a Lyapunov function, but not a strict one, one misses estimates for potential energy Introduce the function, which is indeed a modified energy V ε (t) = E(t) + ε u(t) p 1 u(t)u (t)
38 The scalar case We shall prove that it is a Lyapunov function, which is equivalent to the original one for sufficiently small ε, uniformly with respect to t. Let us first check that V ε is equivalent to E as above stated. Noting that since E is non increasing, we have u(t) 2 2ν 1 E(t) 2ν 1 E(0) for all t 0, so that V ε (t) E(t) ε(2ν 1 p 1 u(t) 2 E(0)) + ε u (t) ε max(2ν 1 E(0)) p 1, 1)E(t) εc E(0) E(t) t 0 where C E(0) > 1. Therefore we have (1 C E(0) ε)e(t) V ε (t) (1 C E(0) ε)e(t) t 0 Hence V ε is equivalent to E for sufficiently small ε, uniformly with respect to t.
39 The scalar case The additional term ε u(t) p 1 u(t)u (t) is build to "create" a potential term which is negative when one differentiates it with respect to t, so that the modified energy is a strict Lyapunov function Let us precise what it means (not with all the details but just the idea):
40 The scalar case When one differentiates V ε, two good terms appear: u (t) p+1 which is due to the contribution E and thus to the dissipation relation εν u(t) p+1 which is the important term which appears in the expression ε u p 1 uu The other terms can be absorbed in this two dominant terms, so that V ε(t) c 1 u (t) p+1 εc 2 u(t) p+1 where c 1 > 0, c 2 > 0. Hence we have where c > 0. V ε(t) cv ε (p+1)/2 (t)
41 The scalar case Thus since p > 1, we deduce that (V (1 p)/2) ε ) (t) c(p 1)/2 > 0 so that for ε sufficiently small, and t sufficiently large ( V ε (t) V (1 p)/2 ε ) 2/(p 1) (0) + c 3 t CE(0) t 2/(p 1) so that since V ε E uniformly in t for ε sufficiently small, we have for t sufficiently large E(t) C E(0) t 2/(p 1)
42 The scalar case an alternative proof based on integral inequalities Haraux 1978, Komornik Indeed one can show that there exists C 0 > 0 (not depending on E 0 such that T S E (p 1)/2 (t)e(t) dt C 0 E(S) 0 S T. where p > 1. One can show that if E is a nonnegative, non increasing, absolutely continuous function satisfying the above inequality then there exists C E( 0) such that Remark E(t) C E(0) t 2/(p 1) for sufficiently large t. If p = 1 in the above integral inequality, then one can show that E is exponentially decaying at infinity.
43 The scalar case Let us give a proof for p = 1 (linear damping case), that E is exponentially decaying Set (this is well-defined thanks to our integral inequality) k(t) = From our assumption, we have Differentiating k, we get t E(s) ds for t 0 k(t) C 0 E(t) t 0 Thus we have k (t) = E(t) k(t)/c 0 t 0 k(t) k(0)e t/c 0 C 0 E(0)e t/c 0 t 0
44 The scalar case Using non negativity of E and its decaying property we have C 0 E(t) t t C 0 E(s) ds k(t C 0 ) C 0 E(0)e (t C 0)/C 0 t 0, We finally get E(t) E(0)e 1 t/c 0 t 0
45 The scalar case For p > 1, we also write a differential inequality for a suitable function k p, more precisely for k p (t) = the differential inequality being t E (p+1)/2 (s) ds k p(t) = E (p+1)/2 (t) C (p+1)/2 0 k (p+1)/2 p (t) t 0 One can get the announced power-like decay as for its proof for the Lyapunov function V ε.
46 The scalar case What about proving a lower estimate for the energy? there exists c > 0 such that E = u p+1 ce (p+1)/2 = after integration (E (1 p)/2 ) c = E(t) c 1 t 2/(p 1) for sufficiently large t. So this proves the optimality of the upper estimate.
47 The scalar case Can this type of results be extended for general feedbacks, i.e. with "arbitrary" growth close to 0? How to build an appropriate a suitable weight for integral inequalities for general feedbacks? How to deduce a simple computable optimal decay rate from the weighted integral inequality if we get one? And further can it be extended to infinite dimensional systems, nonlocal dissipation...?
48 The optimal-weight convexity method Sommaire 1 General motivations for dissipative vibrating systems 2 Dissipative systems in finite dimensions 3 Infinite dimensional dissipative systems 4 The scalar case 5 The optimal-weight convexity method 6 Lower energy estimates: an energy comparison principle 7 Examples of decay rates and optimality 8 Vectorial case: parameter-dependent upper estimates
49 The optimal-weight convexity method Assume that the function g that measures the behavior of the feedback ρ close to 0 is an odd, strictly increasing smooth function on R. We assume that g (0) = 0, otherwise it is easy to show that E decays exponentially to 0 at (linear damping case). We set H(s) = sg( s) s 0
50 The optimal-weight convexity method Weighted integral inequalities Theorem (A.-B. 2005, 2010) We make the above hypotheses and further assume that H is strictly convex on a right neighborhood of 0, denoted by [0, r 2 0 ], r 0 > 0. Then there is a nonnegative smooth strictly increasing weight function w such that E satisfies the following integral inequality T S w(e)e dt C 0 E(S) 0 S T, Proof Set θ = min ( 1, 1 ). 2 µ
51 The optimal-weight convexity method Let for the moment w be a nonnegative C 1 and strictly increasing function defined from [0, r0 2 ) onto [0, + ). w is going to be an optimal-weight function to be determined later on We multiply the left hand side of the equation u + νu + f (u) + ρ(u ) = 0 by w(e(t))u(t) and integrate the resulting equation on [S, T ]. Since E is nonincreasing, w is nondecreasing and thanks to our assumption on f, this gives
52 The optimal-weight convexity method continued. T T θ Ew(E) dt S S 1 T 2 νθ 4 T S T S S w(e) u 2 dt 1 2 w (E) E u u dt 1 2 w(e) u 2 dt + 1 4νθ T S T S w(e)ρ(u )u+ [ w(e)u u ] T S w(e) ρ(u ) 2 dt+ w(e) u 2 dt + 1 ν E(S)w(E(S)), 0 S T. where θ is defined as above.
53 The optimal-weight convexity method continued. Thus, we have T S Ew(E) dt 2 θ 1 2νθ 2 T S T S w(e) u 2 dt+ w(e) ρ(u ) 2 dt + 2 νθ E(S)w(E(S)), 0 S T.
54 The optimal-weight convexity method The next steps will be devoted to control the weighted integrals of the linear kinetic and nonlinear kinetic energies in the right hand side of the above inequality.
55 The optimal-weight convexity method We first remark that thanks to our hypotheses on ρ, we have (up to the positive constants c and C, which may change) { c s ρ(s) C s, s r 0, c g( s ) ρ(s) C g 1 ( s ), s r 0, We set Ĥ(x) = { H(x), if x [0, r 2 0 ], +, if x R [0, r 2 0 ], We define Ĥ as the convex conjugate function of Ĥ (also called the Legendre transform), i.e. Ĥ (y) = sup(xy Ĥ(x)) We recall that the following inequality, called Young s inequality AB Ĥ(A) + Ĥ (B) A 0, B 0
56 The optimal-weight convexity method First step: estimate of the linear kinetic energy For u r 0, we have H( u 2 ) = u g( u ) 1 c u ρ(u ). This, together with Young s inequality imply for u r 0 w(e) u 2 Ĥ ( w(e) ) + H( u 2 ) Ĥ ( w(e) ) + 1 c u ρ(u ). On the other hand, we have w(e) u 2 1 c w(e)u ρ(u ), for u r 0.
57 The optimal-weight convexity method continued. Combining the above two inequalities and the dissipation relation, we obtain T S w(e) u 2 dt T S Ĥ ( w(e) ) dt + 1 c [ 1 + w(e(s)) ] E(S), 0 S T.
58 The optimal-weight convexity method continued. Second step 2: estimate of the nonlinear kinetic energy For u r 0, we have, thanks to Young s inequality w(e) ρ(u ) 2 C 2 Ĥ ( w(e) ) + ( ρ(u ) H 2 ) C Ĥ ( w(e) ) + 1 C u ρ(u ), for u r 0. On the other hand, we have w(e) ρ(u ) 2 Cw(E)u ρ(u ), for u r 0.
59 The optimal-weight convexity method continued. Combining the above two inequalities and the dissipation relations, as above, we have T S w(e) ρ(u ) 2 dt T C 2 Ĥ ( w(e) ) dt + C [ 1 + w(e(s)) ] E(S), 0 S T. S Using the above estimates, we obtain the estimate
60 The optimal-weight convexity method Proof. (continued) T S T Ew(E) dt β Ĥ ( w(e) ) dt+ S 2 ( + 2 νθ θc + C 2νθ )[ 1 + w(e(s)) ] E(S), 0 S T, where β can be easily computed and does not depend on w. One has to be cautious, β is not only chosen as the constant appearing in the above right hand side, it also should be chosen in relation with E(0) such that the weight function s w(s) is defined in the range [0, E(0)) (indeed we choose it even with a stronger criterium below for technical reasons).
61 The optimal-weight convexity method continued. We can now choose the weight by requesting We define a function L by Then the weight w should satisfy βĥ ( w(e) ) = 1 2 Ew(E) L(y) = Ĥ (y) y that is, if L is invertible L(w(E)) = E/(2β) w(.) = L 1 (./2β).
62 The optimal-weight convexity method continued. Assume for the moment that w is well-defined and satisfies the desired properties (non negativity, strictly increasing...). We saw before that the energy E satisfies T S T Ew(E) dt β Ĥ ( w(e) ) dt+ S 2 ( + 2 νθ θc + C 2νθ )[ 1 + w(e(s)) ] E(S), 0 S T, With this choice of weight function, we deduce that T S Ew(E) dt C 0 E(S), 0 S T where C 0 > 0, so that we prove that E satisfies a generalized weighted inequality.
63 The optimal-weight convexity method We now check that the optimal-weight w is indeed well-defined by this way Proposition (A.-B. 2005) Let g be a given odd, strictly increasing C 1 function from R to R such that g (0) = 0. We assume that there exists r 0 > 0 such that the function H is strictly convex on [0, r0 2 ]. Then the function L defined by Ĥ (y), if y (0, + ), L(y) = y 0, if y = 0, is the strictly increasing continuous onto function from [0, + ) on [0, r0 2 ) given by:
64 The optimal-weight convexity method Proposition (continued) (H ) 1 (y) H((H ) 1 (y)), if y [0, H (r0 2 )], L(y) = y r0 2 H(r 0 2), if y [H (r0 2 ), + ). y The weight function w is thus uniquely determined as w(s) = L 1 ( s 2β ) s [0, 2βr 2 0 ). where β > 0 is a suitable constant (independent on w) which satisfies β > E(0)/(2r 2 0 ). Remark Note that for general dampings the inverse of L is not defined on all R, and the further requested computations are not explicit, contrarily to the linear or polynomial case.
65 The optimal-weight convexity method Upper energy estimates from the weighted integral inequality Hence E is a nonnegative, nonincreasing continuous function satisfying a weighted integral inequality. In the polynomial case, it leads to a polynomial decay rate. The situation is more tricky in the general case. We still assume that H is strictly convex on [0, r0 2 ]. Let us introduce the following important function (it will be also useful to classify the feedbacks and for optimality) Λ H on (0, r0 2 ] defined by Λ H (x) = H(x) xh (x).
66 The optimal-weight convexity method Theorem (A.-B. JDE 2010) Let H be a given strictly convex C 1 function from [0, r0 2 ] to R such that H(0) = H (0) = 0, where r 0 > 0 is sufficiently small. We define Ĥ, L and Λ H as before. Assume E be a given nonincreasing, absolutely continuous, nonnegative real function defined on [0, + ), C 0 > 0 be a fixed real number and β > 0 a given real number such that E satisfies the nonlinear Gronwall inequality T S under the condition E(t)L 1 ( E(t) 2β ) dt C 0E(S), 0 S T. 0 < E(0) 2L(H (r 2 0 )) β,
67 The optimal-weight convexity method Theorem (continued) Then, if lim sup x 0 + Λ H (x) < 1, E decays at infinity as follows: E(t) 2β(H ) 1( DC 0 ), for sufficiently large t t where D is a positive constant which does not depend on E(0). Otherwise, we have a general decay rate (see A.-B. 2005). Remark If g(s) = s p 1 s, then H(s) = s (p+1)/2 so that H (s) = 1 (p + 1)s(p 1)/2 2 Thus we recover the upper estimate E(t) C(E(0))t 2/(p 1).
68 The optimal-weight convexity method Moreover one can show that the weight function is of the form w(s) = Cs (p 1)/2 that is we recover as a peculiar example the weight introduced by Haraux-Komornik. If g is such that lim sup Λ H = 1, this means that the feedback is "close" to a linear feedback as for instance g(x) = x(ln( 1 x )) p where p > 0, x close to 0 + We will come back later on this criterium. Under this criterium, we have a simple upper estimate. Next questions: Determine a sharp lower estimate Optimality of the upper estimate?
69 The optimal-weight convexity method Moreover one can show that the weight function is of the form w(s) = Cs (p 1)/2 that is we recover as a peculiar example the weight introduced by Haraux-Komornik. If g is such that lim sup Λ H = 1, this means that the feedback is "close" to a linear feedback as for instance g(x) = x(ln( 1 x )) p where p > 0, x close to 0 + We will come back later on this criterium. Under this criterium, we have a simple upper estimate. Next questions: Determine a sharp lower estimate Optimality of the upper estimate?
70 The optimal-weight convexity method Moreover one can show that the weight function is of the form w(s) = Cs (p 1)/2 that is we recover as a peculiar example the weight introduced by Haraux-Komornik. If g is such that lim sup Λ H = 1, this means that the feedback is "close" to a linear feedback as for instance g(x) = x(ln( 1 x )) p where p > 0, x close to 0 + We will come back later on this criterium. Under this criterium, we have a simple upper estimate. Next questions: Determine a sharp lower estimate Optimality of the upper estimate?
71 Lower energy estimates: an energy comparison principle Sommaire 1 General motivations for dissipative vibrating systems 2 Dissipative systems in finite dimensions 3 Infinite dimensional dissipative systems 4 The scalar case 5 The optimal-weight convexity method 6 Lower energy estimates: an energy comparison principle 7 Examples of decay rates and optimality 8 Vectorial case: parameter-dependent upper estimates
72 Lower energy estimates: an energy comparison principle Lower energy estimates: energy comparison principle We prove Lemma (A.-B. JDE 2010) Assume that f is a continuous and locally Lipschitz function on R which satisfies the above assumptions, and that ρ = g satisfies (A1). Moreover assume that H is increasing and H(0) = 0. Let u be a solution of the scalar ode and E be its energy. Then the following lower estimate holds 1 2 v 2 (t) E(t), t 0, where v is the solution the ODE v + g(v) = 0, v(0) = 2E(0).
73 Lower energy estimates: an energy comparison principle Proof. Thanks to the dissipation relation, and to our assumptions on g, we have E (t) u (t)g(u (t)) = H ( (u ) 2), t 0. Hence, thanks to the ode satisfied by v, we have ( v 2 2 E ) (t) = H ( (u ) 2) H ( (v(t)) 2) H(2E(t)) H ( (v(t)) 2), t 0, and v 2 (0) = 2E(0).
74 Lower energy estimates: an energy comparison principle Proof. Since H is strictly increasing on R, we deduce easily by comparison principles for ODE s that the stated lower energy estimate holds. This an energy comparison principle: we compare the energy of the nonlinear harmonic oscillator u + ν u + f (u) + ρ(u ) = 0. which is a second order ODE, to the energy of the first order ODE v + g(v) = 0
75 Lower energy estimates: an energy comparison principle Theorem (Optimality Theorem, A.-B. JDE 2010) Assume that f is a continuous and locally Lipschitz function on R which satisfies the above hypotheses, and that ρ = g satisfies (A1). We assume that H is strictly convex on [0, r 2 0 ]. Let (u 0, u 1 ) R 2, satisfying 0 < u 1 + u 0 be given, u be the corresponding solution and E be its energy. Moreover assume that either 0 < lim inf x 0 or that there exists µ > 0 such that ( H(µ x) z1 0 < lim inf x 0 µ x x Λ H(x) lim sup Λ H (x) < 1, x 0 1 ) H(y) dy for a certain z 1 (0, z 0 ] (arbitrary)., and lim sup Λ H (x) < 1, x 0
76 Lower energy estimates: an energy comparison principle Theorem (continued) Then the energy of solution satisfies the estimate E(t) = O(v 2 (t)) = O ((H ) 1 ( 1 ) t ), uniformly for large time The proof relies on a key comparison lemma that relies on convexity properties and allows up to compare the lower estimate v 2 (t) which is an "energy-type" estimate, since v 2 is the energy associated to the ode v + g(v) to the upper time-pointwise estimate (H ) 1 (Cste/t)
77 Lower energy estimates: an energy comparison principle Remarks The function Λ H introduces a "classification" of the nonlinearity of the feedbacks lim x 0 Λ H (x) = 0 for g(x) = e 1/x, x > 0 and more generally for very degenerate feedbacks (converging to 0 exponentially for instance) lim x 0 Λ H = 2/(p + 1) for g(x) = x p 1 x and more generally this lim inf (0, 1) for polynomial-logartithmic behavior close to 0 lim x 0 Λ H (x) = 1 for g(x) = x ln(x) 1, x > 0 and more generally for feedbacks which "are close" to a linear behavior at the origin. for feedbacks for which this lim sup < 1, all the solutions have the same asymptotic behavior at infinity
78 Lower energy estimates: an energy comparison principle The condition limsup x 0 +Λ H (x) < 1 for optimality excludes the feedbacks which are "close" to a linear behavior at the origin (see the 3rd example above). Is this sharp? Consider a linear case, that is u + ν u + υu = 0. where ν > 0, υ > 0 are given constants and u is a scalar unknown. Then if υ 2 4ν > 0 there are two linearly independent solutions, namely u 1 (t) = e υ υ 2 4ν 2, t 0, and u 2 (t) = e υ+ υ 2 4ν 2, t 0,
79 Lower energy estimates: an energy comparison principle The energies E i (t) = 1 2 ( u i (t) 2 + ν u i (t) 2), i = 1, 2, t 0, of these two solutions decay exponentially as t goes to, but at different rates and their ratio satisfies E 2 (t) lim t E 1 t) =. This implies different lower and upper bounds and different decay rates depending on the initial data Same behavior for υ 2 4ν = 0, If υ 2 4ν < 0, all the energies of all the solutions decay at the same speed e υt as t goes to. A similar situation (with more cases) arises in the linear vectorial case.
80 Lower energy estimates: an energy comparison principle We conjecture that a similar situation, that is two branches of solutions with different asymptotic behavior may arise for nonlinear feedbacks close to linear growth at the origin that is for feedbacks for which. limsup x 0 +Λ H (x) = 1
81 Examples of decay rates and optimality Sommaire 1 General motivations for dissipative vibrating systems 2 Dissipative systems in finite dimensions 3 Infinite dimensional dissipative systems 4 The scalar case 5 The optimal-weight convexity method 6 Lower energy estimates: an energy comparison principle 7 Examples of decay rates and optimality 8 Vectorial case: parameter-dependent upper estimates
82 Examples of decay rates and optimality Examples of decay rates Example 1 (polynomial case): let g be given by g(x) = x p where p > 1 on (0, r 0 ]. Then E(t) Cβ E(0) t 2 p 1, (3) for t sufficiently large and for all (u 0, u 1 ) in R 2. Moreover this estimate is optimal. Example 2 (exponential case): let g be given by g(x) = e 1 x 2 on (0, r 0 ]. Then E(t) Cβ E(0) (ln(t)) 1, for large t. Moreover, this estimate is optimal.
83 Examples of decay rates and optimality Example 3 (polynomial-logarithmic, close to linear): let g be given by g(x) = x(ln( 1 x )) p where p > 0. Then E(t) C β E(0) e 2( p t DT 0 ) 1/(p+1) t 1/(p+1) (4) for t sufficiently large. Optimality cannot be asserted. Example 4 (faster than any polynomial less than exponential) : vskip 2mm let g be given by g(x) = e (ln( 1 x ))p, 1 < p < 2, x [0, r 0 ]. Then E(t) C β E(0) e 2(ln(t))1/p This estimate is optimal.
84 Vectorial case: parameter-dependent upper estimates Sommaire 1 General motivations for dissipative vibrating systems 2 Dissipative systems in finite dimensions 3 Infinite dimensional dissipative systems 4 The scalar case 5 The optimal-weight convexity method 6 Lower energy estimates: an energy comparison principle 7 Examples of decay rates and optimality 8 Vectorial case: parameter-dependent upper estimates
85 Vectorial case: parameter-dependent upper estimates Similar result with sharp energy decay rates, and optimality results for the vectorial case. These optimal estimates depends on the dimension of the system. Moreover these results applied to semi-discretized PDE s such as the wave or plate equations and give optimal rates of decay. Therefore when applied to semi-discretized PDE s, they are not uniform with respect to the discretization parameter.
86 Vectorial case: parameter-dependent upper estimates We consider the vectorial case: u + Au + f (u) + Bρ(u ) = 0. where the unknown u R n and A is a symmetric, positive definite matrix B = diag(b i ) 1 i n, with b i 0 i {1,..., n} (f (u)) i = f (u i ), (ρ(u )) i = ρ(u i ) ρ : v = (v 1,..., v n ) R n ( g 1 (v 1 ),..., g n (v n ) ). f is a vectorial function of the form: f (u) = (f 1 (u 1 ),..., f n (u n )), u = (u 1,..., u n ) R n.
87 Vectorial case: parameter-dependent upper estimates Recall that the energy is given by E(t) = 1 2 ( u 2 + Au, u + n ) (F(u i )), i=1 and the dissipation relation by n E (t) = b i ρ(u i )u i 0 t 0. i=1
88 Vectorial case: parameter-dependent upper estimates Theorem (A.-B. 2010) Assume that f is continuous and locally Lipschitz on R n and satisfies some growth conditions (as for the scalar case) the functions g i (components of the feedback ρ) satisfies similar assumptions than in the optimality theorem for the scalar case with cg( s ) g i (s) k 1 g( s ) and g(s) k 2 g 1 ( s ) for s 1 and i {1,..., n} the function H(s) = sg( s) is strictly convex on [0, r0 2]
89 Vectorial case: parameter-dependent upper estimates Theorem (continued) Then If 0 < lim inf x 0 then the following estimates hold Λ H(x) lim sup Λ H (x) < 1, x v 2 (t) E(t) C 2 v 2 (t), for t sufficiently large, where C 2 is a positive constant and v is the solution of the ordinary differential equation: v (t) + n k 1 g(v(t)) = 0, v(0) = 2E(0), t 0. Or equivalently, we have E(t) = O(v 2 (t)) = O ((H ) 1 ( D ) t ) uniformly with respect to t, (5) where D is a positive constant.
90 Vectorial case: parameter-dependent upper estimates Theorem (continued) If there exists µ > 0 and z 1 (0, E(0)] such that ( H(µ x) z1 0 < lim inf x 0 µ x x 1 ) H(y) dy, and lim sup Λ H (x) < 1, x 0 then the above estimates hold.
91 Vectorial case: parameter-dependent upper estimates From the above theorem, we deduce easily the following corollary. Theorem (Semi-discretized PDE s) Assume that the vectorial functions defined by f i (.) = f (x i,.) and ρ(.) = ρ(x i,.)for i = 1,... n satisfy the above hypotheses, where the points x i, i = 1,..., n denotes the discretization points. Then the above estimates hold. Remark The decay rate given in the above theorem depends on n and thus on the discretization parameter h. They are not uniform in h. In particular, for prescribed initial data, the constant β behaves as C/h where C depends on the initial data but not on h. Upper estimates independent on h and sharper conditions on the nonlinearity: work in progress with Emmanuel Trélat and Yannick Privat (see course by Emmanuel).
92 Vectorial case: parameter-dependent upper estimates Hence in the finite dimensional case, we have a simple and optimal energy decay rate for finite dimensional systems, without prescribed growth assumptions on the feedbacks it applies in particular to finite dimensional systems coming out from semi-discretization of infinite dimensional nonlinear vibrating systems. this, for fixed discretization parameters, since our estimates in the vectorial case depend on n, i.e. equivalently on h. see Emmanuel s talk for extension to non diagonal matrix B (related to discretized nonlinear damping) Next lecture: infinite dimensional systems
Stability of nonlinear locally damped partial differential equations: the continuous and discretized problems. Part II
. Stability of nonlinear locally damped partial differential equations: the continuous and discretized problems. Part II Fatiha Alabau-Boussouira 1 Emmanuel Trélat 2 1 Univ. de Lorraine, LMAM 2 Univ. Paris
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationOBSERVABILITY INEQUALITY AND DECAY RATE FOR WAVE EQUATIONS WITH NONLINEAR BOUNDARY CONDITIONS
Electronic Journal of Differential Equations, Vol. 27 (27, No. 6, pp. 2. ISSN: 72-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu OBSERVABILITY INEQUALITY AND DECAY RATE FOR WAVE EQUATIONS
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationControllability of linear PDEs (I): The wave equation
Controllability of linear PDEs (I): The wave equation M. González-Burgos IMUS, Universidad de Sevilla Doc Course, Course 2, Sevilla, 2018 Contents 1 Introduction. Statement of the problem 2 Distributed
More informationThe Dirichlet s P rinciple. In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation:
Oct. 1 The Dirichlet s P rinciple In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation: 1. Dirichlet s Principle. u = in, u = g on. ( 1 ) If we multiply
More informationControl, Stabilization and Numerics for Partial Differential Equations
Paris-Sud, Orsay, December 06 Control, Stabilization and Numerics for Partial Differential Equations Enrique Zuazua Universidad Autónoma 28049 Madrid, Spain enrique.zuazua@uam.es http://www.uam.es/enrique.zuazua
More informationLyapunov Stability Theory
Lyapunov Stability Theory Peter Al Hokayem and Eduardo Gallestey March 16, 2015 1 Introduction In this lecture we consider the stability of equilibrium points of autonomous nonlinear systems, both in continuous
More informationMildly degenerate Kirchhoff equations with weak dissipation: global existence and time decay
arxiv:93.273v [math.ap] 6 Mar 29 Mildly degenerate Kirchhoff equations with weak dissipation: global existence and time decay Marina Ghisi Università degli Studi di Pisa Dipartimento di Matematica Leonida
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More informationDecay Rates for Dissipative Wave equations
Published in Ricerche di Matematica 48 (1999), 61 75. Decay Rates for Dissipative Wave equations Wei-Jiu Liu Department of Applied Mechanics and Engineering Sciences University of California at San Diego
More informationConvergence Rate of Nonlinear Switched Systems
Convergence Rate of Nonlinear Switched Systems Philippe JOUAN and Saïd NACIRI arxiv:1511.01737v1 [math.oc] 5 Nov 2015 January 23, 2018 Abstract This paper is concerned with the convergence rate of the
More informationEquilibria with a nontrivial nodal set and the dynamics of parabolic equations on symmetric domains
Equilibria with a nontrivial nodal set and the dynamics of parabolic equations on symmetric domains J. Földes Department of Mathematics, Univerité Libre de Bruxelles 1050 Brussels, Belgium P. Poláčik School
More informationON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM
ON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM OLEG ZUBELEVICH DEPARTMENT OF MATHEMATICS THE BUDGET AND TREASURY ACADEMY OF THE MINISTRY OF FINANCE OF THE RUSSIAN FEDERATION 7, ZLATOUSTINSKY MALIY PER.,
More informationStabilization of the wave equation with localized Kelvin-Voigt damping
Stabilization of the wave equation with localized Kelvin-Voigt damping Louis Tebou Florida International University Miami SEARCDE University of Memphis October 11-12, 2014 Louis Tebou (FIU, Miami) Stabilization...
More informationGLOBAL EXISTENCE AND ENERGY DECAY OF SOLUTIONS TO A PETROVSKY EQUATION WITH GENERAL NONLINEAR DISSIPATION AND SOURCE TERM
Georgian Mathematical Journal Volume 3 (26), Number 3, 397 4 GLOBAL EXITENCE AND ENERGY DECAY OF OLUTION TO A PETROVKY EQUATION WITH GENERAL NONLINEAR DIIPATION AND OURCE TERM NOUR-EDDINE AMROUN AND ABBE
More informationRégularité des équations de Hamilton-Jacobi du premier ordre et applications aux jeux à champ moyen
Régularité des équations de Hamilton-Jacobi du premier ordre et applications aux jeux à champ moyen Daniela Tonon en collaboration avec P. Cardaliaguet et A. Porretta CEREMADE, Université Paris-Dauphine,
More informationEstimates for the resolvent and spectral gaps for non self-adjoint operators. University Bordeaux
for the resolvent and spectral gaps for non self-adjoint operators 1 / 29 Estimates for the resolvent and spectral gaps for non self-adjoint operators Vesselin Petkov University Bordeaux Mathematics Days
More informationOn Chern-Simons-Schrödinger equations including a vortex point
On Chern-Simons-Schrödinger equations including a vortex point Alessio Pomponio Dipartimento di Meccanica, Matematica e Management, Politecnico di Bari Workshop in Nonlinear PDEs Brussels, September 7
More informationIntroduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems
p. 1/5 Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems p. 2/5 Time-varying Systems ẋ = f(t, x) f(t, x) is piecewise continuous in t and locally Lipschitz in x for all t
More informationMeasure and Integration: Solutions of CW2
Measure and Integration: s of CW2 Fall 206 [G. Holzegel] December 9, 206 Problem of Sheet 5 a) Left (f n ) and (g n ) be sequences of integrable functions with f n (x) f (x) and g n (x) g (x) for almost
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationh(x) lim H(x) = lim Since h is nondecreasing then h(x) 0 for all x, and if h is discontinuous at a point x then H(x) > 0. Denote
Real Variables, Fall 4 Problem set 4 Solution suggestions Exercise. Let f be of bounded variation on [a, b]. Show that for each c (a, b), lim x c f(x) and lim x c f(x) exist. Prove that a monotone function
More informationStability of an abstract wave equation with delay and a Kelvin Voigt damping
Stability of an abstract wave equation with delay and a Kelvin Voigt damping University of Monastir/UPSAY/LMV-UVSQ Joint work with Serge Nicaise and Cristina Pignotti Outline 1 Problem The idea Stability
More informationStability of Stochastic Differential Equations
Lyapunov stability theory for ODEs s Stability of Stochastic Differential Equations Part 1: Introduction Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH December 2010
More informationGlobal Solutions for a Nonlinear Wave Equation with the p-laplacian Operator
Global Solutions for a Nonlinear Wave Equation with the p-laplacian Operator Hongjun Gao Institute of Applied Physics and Computational Mathematics 188 Beijing, China To Fu Ma Departamento de Matemática
More informationOn 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 informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
More informationLecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.
Lecture 4 Chapter 4: Lyapunov Stability Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 4 p. 1/86 Autonomous Systems Consider the autonomous system ẋ
More informationNonlinear Analysis 71 (2009) Contents lists available at ScienceDirect. Nonlinear Analysis. journal homepage:
Nonlinear Analysis 71 2009 2744 2752 Contents lists available at ScienceDirect Nonlinear Analysis journal homepage: www.elsevier.com/locate/na A nonlinear inequality and applications N.S. Hoang A.G. Ramm
More informationCUTOFF RESOLVENT ESTIMATES AND THE SEMILINEAR SCHRÖDINGER EQUATION
CUTOFF RESOLVENT ESTIMATES AND THE SEMILINEAR SCHRÖDINGER EQUATION HANS CHRISTIANSON Abstract. This paper shows how abstract resolvent estimates imply local smoothing for solutions to the Schrödinger equation.
More informationNonlinear Systems Theory
Nonlinear Systems Theory Matthew M. Peet Arizona State University Lecture 2: Nonlinear Systems Theory Overview Our next goal is to extend LMI s and optimization to nonlinear systems analysis. Today we
More informationContinuous Functions on Metric Spaces
Continuous Functions on Metric Spaces Math 201A, Fall 2016 1 Continuous functions Definition 1. Let (X, d X ) and (Y, d Y ) be metric spaces. A function f : X Y is continuous at a X if for every ɛ > 0
More informationAnalysis Finite and Infinite Sets The Real Numbers The Cantor Set
Analysis Finite and Infinite Sets Definition. An initial segment is {n N n n 0 }. Definition. A finite set can be put into one-to-one correspondence with an initial segment. The empty set is also considered
More informationExistence and uniqueness of solutions for nonlinear ODEs
Chapter 4 Existence and uniqueness of solutions for nonlinear ODEs In this chapter we consider the existence and uniqueness of solutions for the initial value problem for general nonlinear ODEs. Recall
More informationAn asymptotic ratio characterization of input-to-state stability
1 An asymptotic ratio characterization of input-to-state stability Daniel Liberzon and Hyungbo Shim Abstract For continuous-time nonlinear systems with inputs, we introduce the notion of an asymptotic
More informationStationary Kirchhoff equations with powers by Emmanuel Hebey (Université de Cergy-Pontoise)
Stationary Kirchhoff equations with powers by Emmanuel Hebey (Université de Cergy-Pontoise) Lectures at the Riemann center at Varese, at the SNS Pise, at Paris 13 and at the university of Nice. June 2017
More informationObstacle Problems Involving The Fractional Laplacian
Obstacle Problems Involving The Fractional Laplacian Donatella Danielli and Sandro Salsa January 27, 2017 1 Introduction Obstacle problems involving a fractional power of the Laplace operator appear in
More information2. Dual space is essential for the concept of gradient which, in turn, leads to the variational analysis of Lagrange multipliers.
Chapter 3 Duality in Banach Space Modern optimization theory largely centers around the interplay of a normed vector space and its corresponding dual. The notion of duality is important for the following
More informationTHE WAVE EQUATION. d = 1: D Alembert s formula We begin with the initial value problem in 1 space dimension { u = utt u xx = 0, in (0, ) R, (2)
THE WAVE EQUATION () The free wave equation takes the form u := ( t x )u = 0, u : R t R d x R In the literature, the operator := t x is called the D Alembertian on R +d. Later we shall also consider the
More informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
More informationElliptic Kirchhoff equations
Elliptic Kirchhoff equations David ARCOYA Universidad de Granada Sevilla, 8-IX-2015 Workshop on Recent Advances in PDEs: Analysis, Numerics and Control In honor of Enrique Fernández-Cara for his 60th birthday
More informationON THE PATHWISE UNIQUENESS OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS
PORTUGALIAE MATHEMATICA Vol. 55 Fasc. 4 1998 ON THE PATHWISE UNIQUENESS OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS C. Sonoc Abstract: A sufficient condition for uniqueness of solutions of ordinary
More information1 The Observability Canonical Form
NONLINEAR OBSERVERS AND SEPARATION PRINCIPLE 1 The Observability Canonical Form In this Chapter we discuss the design of observers for nonlinear systems modelled by equations of the form ẋ = f(x, u) (1)
More informationStabilization for the Wave Equation with Variable Coefficients and Balakrishnan-Taylor Damping. Tae Gab Ha
TAIWANEE JOURNAL OF MATHEMATIC Vol. xx, No. x, pp., xx 0xx DOI: 0.650/tjm/788 This paper is available online at http://journal.tms.org.tw tabilization for the Wave Equation with Variable Coefficients and
More informationIntroduction LECTURE 1
LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in
More informationAsymptotic Behavior for Semi-Linear Wave Equation with Weak Damping
Int. Journal of Math. Analysis, Vol. 7, 2013, no. 15, 713-718 HIKARI Ltd, www.m-hikari.com Asymptotic Behavior for Semi-Linear Wave Equation with Weak Damping Ducival Carvalho Pereira State University
More informationNonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems. p. 1/1
Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems p. 1/1 p. 2/1 Converse Lyapunov Theorem Exponential Stability Let x = 0 be an exponentially stable equilibrium
More informationInfluence of damping on hyperbolic equations with parabolic degeneracy
University of New Orleans ScholarWorks@UNO Mathematics Faculty Publications Department of Mathematics 1-1-212 Influence of damping on hyperbolic equations with parabolic degeneracy Katarzyna Saxton Loyola
More informationMODULUS OF CONTINUITY OF THE DIRICHLET SOLUTIONS
MODULUS OF CONTINUITY OF THE DIRICHLET SOLUTIONS HIROAKI AIKAWA Abstract. Let D be a bounded domain in R n with n 2. For a function f on D we denote by H D f the Dirichlet solution, for the Laplacian,
More informationExponential stability of abstract evolution equations with time delay feedback
Exponential stability of abstract evolution equations with time delay feedback Cristina Pignotti University of L Aquila Cortona, June 22, 2016 Cristina Pignotti (L Aquila) Abstract evolutions equations
More informationSharp estimates of bounded solutions to some semilinear second order dissipative equations
Sharp estimates of ounded solutions to some semilinear second order dissipative equations Cyrine Fitouri & Alain Haraux Astract. Let H, V e two real Hilert spaces such that V H with continuous and dense
More informationAUTOMATIC CONTROL. Andrea M. Zanchettin, PhD Spring Semester, Introduction to Automatic Control & Linear systems (time domain)
1 AUTOMATIC CONTROL Andrea M. Zanchettin, PhD Spring Semester, 2018 Introduction to Automatic Control & Linear systems (time domain) 2 What is automatic control? From Wikipedia Control theory is an interdisciplinary
More informationAn introduction to some aspects of functional analysis
An introduction to some aspects of functional analysis Stephen Semmes Rice University Abstract These informal notes deal with some very basic objects in functional analysis, including norms and seminorms
More informationNonlinear Control Lecture 5: Stability Analysis II
Nonlinear Control Lecture 5: Stability Analysis II Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 Farzaneh Abdollahi Nonlinear Control Lecture 5 1/41
More informationCongurations of periodic orbits for equations with delayed positive feedback
Congurations of periodic orbits for equations with delayed positive feedback Dedicated to Professor Tibor Krisztin on the occasion of his 60th birthday Gabriella Vas 1 MTA-SZTE Analysis and Stochastics
More informationBLOW-UP OF SOLUTIONS FOR VISCOELASTIC EQUATIONS OF KIRCHHOFF TYPE WITH ARBITRARY POSITIVE INITIAL ENERGY
Electronic Journal of Differential Equations, Vol. 6 6, No. 33, pp. 8. ISSN: 7-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu BLOW-UP OF SOLUTIONS FOR VISCOELASTIC EQUATIONS OF KIRCHHOFF
More information3 Integration and Expectation
3 Integration and Expectation 3.1 Construction of the Lebesgue Integral Let (, F, µ) be a measure space (not necessarily a probability space). Our objective will be to define the Lebesgue integral R fdµ
More informationExistence and Uniqueness
Chapter 3 Existence and Uniqueness An intellect which at a certain moment would know all forces that set nature in motion, and all positions of all items of which nature is composed, if this intellect
More informationin Bounded Domains Ariane Trescases CMLA, ENS Cachan
CMLA, ENS Cachan Joint work with Yan GUO, Chanwoo KIM and Daniela TONON International Conference on Nonlinear Analysis: Boundary Phenomena for Evolutionnary PDE Academia Sinica December 21, 214 Outline
More informationInput to state Stability
Input to state Stability Mini course, Universität Stuttgart, November 2004 Lars Grüne, Mathematisches Institut, Universität Bayreuth Part III: Lyapunov functions and quantitative aspects ISS Consider with
More informationMicrolocal Methods in X-ray Tomography
Microlocal Methods in X-ray Tomography Plamen Stefanov Purdue University Lecture I: Euclidean X-ray tomography Mini Course, Fields Institute, 2012 Plamen Stefanov (Purdue University ) Microlocal Methods
More informationContents: 1. Minimization. 2. The theorem of Lions-Stampacchia for variational inequalities. 3. Γ -Convergence. 4. Duality mapping.
Minimization Contents: 1. Minimization. 2. The theorem of Lions-Stampacchia for variational inequalities. 3. Γ -Convergence. 4. Duality mapping. 1 Minimization A Topological Result. Let S be a topological
More informationBoundedly complete weak-cauchy basic sequences in Banach spaces with the PCP
Journal of Functional Analysis 253 (2007) 772 781 www.elsevier.com/locate/jfa Note Boundedly complete weak-cauchy basic sequences in Banach spaces with the PCP Haskell Rosenthal Department of Mathematics,
More informationSPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS
SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties
More informationMinimal periods of semilinear evolution equations with Lipschitz nonlinearity
Minimal periods of semilinear evolution equations with Lipschitz nonlinearity James C. Robinson a Alejandro Vidal-López b a Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. b Departamento
More informationORDINARY DIFFERENTIAL EQUATIONS
ORDINARY DIFFERENTIAL EQUATIONS GABRIEL NAGY Mathematics Department, Michigan State University, East Lansing, MI, 48824. JANUARY 3, 25 Summary. This is an introduction to ordinary differential equations.
More informationOutline. Input to state Stability. Nonlinear Realization. Recall: _ Space. _ Space: Space of all piecewise continuous functions
Outline Input to state Stability Motivation for Input to State Stability (ISS) ISS Lyapunov function. Stability theorems. M. Sami Fadali Professor EBME University of Nevada, Reno 1 2 Recall: _ Space _
More informationThe harmonic map flow
Chapter 2 The harmonic map flow 2.1 Definition of the flow The harmonic map flow was introduced by Eells-Sampson in 1964; their work could be considered the start of the field of geometric flows. The flow
More informationA generic property of families of Lagrangian systems
Annals of Mathematics, 167 (2008), 1099 1108 A generic property of families of Lagrangian systems By Patrick Bernard and Gonzalo Contreras * Abstract We prove that a generic Lagrangian has finitely many
More informationConvergence and sharp thresholds for propagation in nonlinear diffusion problems
J. Eur. Math. Soc. 12, 279 312 c European Mathematical Society 2010 DOI 10.4171/JEMS/198 Yihong Du Hiroshi Matano Convergence and sharp thresholds for propagation in nonlinear diffusion problems Received
More informationSecond order forward-backward dynamical systems for monotone inclusion problems
Second order forward-backward dynamical systems for monotone inclusion problems Radu Ioan Boţ Ernö Robert Csetnek March 6, 25 Abstract. We begin by considering second order dynamical systems of the from
More informationLECTURE 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 informationExistence and uniqueness of solutions for a diffusion model of host parasite dynamics
J. Math. Anal. Appl. 279 (23) 463 474 www.elsevier.com/locate/jmaa Existence and uniqueness of solutions for a diffusion model of host parasite dynamics Michel Langlais a and Fabio Augusto Milner b,,1
More informationAn Operator Theoretical Approach to Nonlocal Differential Equations
An Operator Theoretical Approach to Nonlocal Differential Equations Joshua Lee Padgett Department of Mathematics and Statistics Texas Tech University Analysis Seminar November 27, 2017 Joshua Lee Padgett
More informationScalar Conservation Laws and First Order Equations Introduction. Consider equations of the form. (1) u t + q(u) x =0, x R, t > 0.
Scalar Conservation Laws and First Order Equations Introduction. Consider equations of the form (1) u t + q(u) x =, x R, t >. In general, u = u(x, t) represents the density or the concentration of a physical
More informationStabilization of persistently excited linear systems
Stabilization of persistently excited linear systems Yacine Chitour Laboratoire des signaux et systèmes & Université Paris-Sud, Orsay Exposé LJLL Paris, 28/9/2012 Stabilization & intermittent control Consider
More informationẋ = f(x, y), ẏ = g(x, y), (x, y) D, can only have periodic solutions if (f,g) changes sign in D or if (f,g)=0in D.
4 Periodic Solutions We have shown that in the case of an autonomous equation the periodic solutions correspond with closed orbits in phase-space. Autonomous two-dimensional systems with phase-space R
More informationEntropy-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 informationFunctional Analysis. Franck Sueur Metric spaces Definitions Completeness Compactness Separability...
Functional Analysis Franck Sueur 2018-2019 Contents 1 Metric spaces 1 1.1 Definitions........................................ 1 1.2 Completeness...................................... 3 1.3 Compactness......................................
More informationOrdinary Differential Equations II
Ordinary Differential Equations II February 9 217 Linearization of an autonomous system We consider the system (1) x = f(x) near a fixed point x. As usual f C 1. Without loss of generality we assume x
More informationSeptember Math Course: First Order Derivative
September Math Course: First Order Derivative Arina Nikandrova Functions Function y = f (x), where x is either be a scalar or a vector of several variables (x,..., x n ), can be thought of as a rule which
More informationA TWO PARAMETERS AMBROSETTI PRODI PROBLEM*
PORTUGALIAE MATHEMATICA Vol. 53 Fasc. 3 1996 A TWO PARAMETERS AMBROSETTI PRODI PROBLEM* C. De Coster** and P. Habets 1 Introduction The study of the Ambrosetti Prodi problem has started with the paper
More informationAsymptotic Behavior of a Hyperbolic-parabolic Coupled System Arising in Fluid-structure Interaction
International Series of Numerical Mathematics, Vol. 154, 445 455 c 2006 Birkhäuser Verlag Basel/Switzerland Asymptotic Behavior of a Hyperbolic-parabolic Coupled System Arising in Fluid-structure Interaction
More informationQualitative behavior of global solutions to some nonlinear fourth order differential equations
Qualitative behavior of global solutions to some nonlinear fourth order differential equations Elvise BERCHIO - Alberto FERRERO - Filippo GAZZOLA - Paschalis KARAGEORGIS Abstract We study global solutions
More informationPropagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R
Propagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R P. Poláčik School of Mathematics, University of Minnesota Minneapolis, MN 55455 Abstract We consider semilinear
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 information(2m)-TH MEAN BEHAVIOR OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS UNDER PARAMETRIC PERTURBATIONS
(2m)-TH MEAN BEHAVIOR OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS UNDER PARAMETRIC PERTURBATIONS Svetlana Janković and Miljana Jovanović Faculty of Science, Department of Mathematics, University
More informationForchheimer Equations in Porous Media - Part III
Forchheimer Equations in Porous Media - Part III Luan Hoang, Akif Ibragimov Department of Mathematics and Statistics, Texas Tech niversity http://www.math.umn.edu/ lhoang/ luan.hoang@ttu.edu Applied Mathematics
More informationAsymptotic Convergence of the Steepest Descent Method for the Exponential Penalty in Linear Programming
Journal of Convex Analysis Volume 2 (1995), No.1/2, 145 152 Asymptotic Convergence of the Steepest Descent Method for the Exponential Penalty in Linear Programming R. Cominetti 1 Universidad de Chile,
More informationEXISTENCE OF SOLUTIONS FOR CROSS CRITICAL EXPONENTIAL N-LAPLACIAN SYSTEMS
Electronic Journal of Differential Equations, Vol. 2014 (2014), o. 28, pp. 1 10. ISS: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu EXISTECE OF SOLUTIOS
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 informationLecture 1: Entropy, convexity, and matrix scaling CSE 599S: Entropy optimality, Winter 2016 Instructor: James R. Lee Last updated: January 24, 2016
Lecture 1: Entropy, convexity, and matrix scaling CSE 599S: Entropy optimality, Winter 2016 Instructor: James R. Lee Last updated: January 24, 2016 1 Entropy Since this course is about entropy maximization,
More informationTakens embedding theorem for infinite-dimensional dynamical systems
Takens embedding theorem for infinite-dimensional dynamical systems James C. Robinson Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. E-mail: jcr@maths.warwick.ac.uk Abstract. Takens
More informationfor all subintervals I J. If the same is true for the dyadic subintervals I D J only, we will write ϕ BMO d (J). In fact, the following is true
3 ohn Nirenberg inequality, Part I A function ϕ L () belongs to the space BMO() if sup ϕ(s) ϕ I I I < for all subintervals I If the same is true for the dyadic subintervals I D only, we will write ϕ BMO
More informationNotes on uniform convergence
Notes on uniform convergence Erik Wahlén erik.wahlen@math.lu.se January 17, 2012 1 Numerical sequences We begin by recalling some properties of numerical sequences. By a numerical sequence we simply mean
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 informationPeriodic solutions of weakly coupled superlinear systems
Periodic solutions of weakly coupled superlinear systems Alessandro Fonda and Andrea Sfecci Abstract By the use of a higher dimensional version of the Poincaré Birkhoff theorem, we are able to generalize
More informationObservability for deterministic systems and high-gain observers
Observability for deterministic systems and high-gain observers design. Part 1. March 29, 2011 Introduction and problem description Definition of observability Consequences of instantaneous observability
More informationON PARABOLIC HARNACK INEQUALITY
ON PARABOLIC HARNACK INEQUALITY JIAXIN HU Abstract. We show that the parabolic Harnack inequality is equivalent to the near-diagonal lower bound of the Dirichlet heat kernel on any ball in a metric measure-energy
More information