Second order sucient optimality condition for a. nonlinear elliptic control problem. Eduardo Casas and Fredi Troltzsch and Andreas Unger

Size: px
Start display at page:

Download "Second order sucient optimality condition for a. nonlinear elliptic control problem. Eduardo Casas and Fredi Troltzsch and Andreas Unger"

Transcription

1 Second order sucient optimality condition for a nonlinear elliptic control problem by Eduardo Casas and Fredi Troltzsch and Andreas Unger Technical University of Chemnitz wickau Faculty of Mathematics 97 Chemnitz Germany Introduction Optimal control problems governed by nonlinear elliptic partial dierential equations have already been considered by many authors. We refer only to the papers of Casas [], [3], to the book of Tiba [4], and to the references therein. Meanwhile, the existence of optimal controls and rst order necessary optimality conditions of the maximum{ principle type are well investigated. It is known that in the case of nonlinear equations the rst order conditions are not in general sucient for optimality. In this paper we are going to derive a second{order sucient optimality condition for a class of semilinear elliptic boundary control problems. For parabolic boundary control problems second{ order sucient conditions were established in the papers by Goldberg and Troltzsch [6]. It is more or less obvious that these conditions can be transferred by the same technique to elliptic problems. However, a comparison of the results in [6] with second order conditions for optimization problems in spaces of nite dimension reveals that the gap between the sucient optimality conditions in [6] and corresponding second order necessary optimality conditions is quite large. Taking into account the active set of optimal controls, this gap can be partially closed. The known diculty in the theory of sucient optimality conditions for extremal problems in function spaces is the so{called two{norm discrepancy. This problem was resolved, for instance, by Ioe [8] and Maurer []. In recent years further considerations have shown that some weaker sucient optimality conditions can be established for abstract optimization problems in Banach spaces and for optimal control problems governed by ordinary dierential equations. We refer, for instance, to Dontchev et al. [4] and Malanowski []. The main idea is to introduce a third norm taking into account the measure of positivity of some terms occuring in the variational inequality. We will use this idea to derive of a second{order sucient condition in the elliptic case.

2 Formulation of the optimal control problem We consider the optimal control problem to minimize F (w; u) '(w(x))dx (u(x); w(x))ds x () subject to and to the constraints on the control u 4w(x) w(x) in (x) b(w(x); u(x)) on u U ad fv L ()j < a v(x) b < a. e. on g: (4) U ad is a bounded, convex, and closed subset of L (). The partial dierential equation is considered in the following weak sense. Denition A function w W () is said to be a weak solution of the partial dierential equation (), (3), if for all v W () the equation is valid. @x i wv! dx To make this well dened we impose the following assumptions: b(w; u)vds x (5) (A) The function ' : IR! IR is of class C ;. The functions : IR IR! IR and b : IRIR! IR are twice continuously dierentiable with respect to both arguments. Furthermore, the function b is monotone decreasing and all partial derivatives of and b up to the second order are globally Lipschitz{continuous. This property of the function b guarantees the existence of a unique weak solution of the partial dierential equation (), (3) for each xed u U ad (cf. Kinderlehrer and Stampacchia [9]). (A) The set IR n is assumed to be a bounded domain with a sucently smooth boundary such that the weak solution of the partial dierential equation (), (3) can be shown to be an element of W p () with some xed p > n. Then known embedding theorems ensure that the solution is continuous on the closure of. (A3) Moreover, we suppose that the following estimate is valid for all admissible controls u ; u U ad and the corresponding solutions of the partial dierential equation w and w : kw w k W p () cku u k Lp() Cku u k L(): (6)

3 For instance, suppose that is a bounded domain with a C {boundary. Using results of Groger [7], the following regularity result can be proved Lemma The weak solution w of the boundary value problem 4w w g on (8) is an element of W p (), provided the function g is an element of L p (). Let w i denote the solution of the equation (7), (8) corresponding to g i L p () for i ;. Then the estimate kw w k W p () ckg g k Lp() is valid, where c does not depend on the choice of g i. Proof: We will outline only the main idea of the proof. For the details we refer to the paper of Groger [7]. The domain is bounded and has a C {boundary. Therefore, it can be shown that the dierential operator represented by the equation (7), (8) is a continuous mapping from Wp () onto the dual space of W q (), pq. In the known way, any g L p() can be regarded as an element of (Wq ()). The statement of the lemma is now a conclusion of the properties of the dierential operator. The uniqueness of the solution in W () ensures also the invertibility of the dierential operator between Wp () and (Wq ()). We can now apply the Banach theorem to derive the continuity of the inverse operator. This lemma allows us to ensure the assumed regularity of the weak solution. We proceed in the following way. The unique solution w W () of the state equation (), (3) exists due to the assumptions on b and. The trace theorem shows that the trace wj is in W (). This space is embedded in L p () with some p > because of the embedding theorem. Owing to the Lipschitz property of b and to u L (), we can identify the function b(w(); u()) with a function g L p (). The regularity result ensures that the solution w belongs to Wp (). We repeat the considerations starting with p instead of. In this way the solution w of the state equation is seen to be an element of Wp () with some p > p. Repeating this bootstrapping argument suciently often, we construct a sequence fp i g, where p i > p i. It can be shown (see Eppler and Unger [5]) that there exists a constant c > such that p i > p i c for all i. So we are sure to arrive at p > n after a nite number of steps. The assumed estimate follows from the estimate in the regularity result and the monotonicity of the function b with respect to the rst variable. As a further corollary we have that the norm of all feasible states is uniformly bounded in Wp () since of U ad is bounded. If we assume in addition that the cost functional is convex and the function b is linear both with respect to the control u, the existence of a solution of the optimal control problem follows by standard methods (cf. for instance Eppler and Unger [5]). However, we shall not rely on this assumption. We just suppose that a xed reference control u is given. Let us introduce the Langrange{function for our optimal control problem in the following way: L(w; u; y) '(w)dx (w; u)ds x 3 wy i b(w; u)yds x : (9)

4 This function is well dened and twice continuously dierentiable with respect to the space C() \ W () L (). Standard considerations apply to derive the rst order necessary optimality conditions at u. They can be written in the equivalent form (adjoint equation) and L w (w ; u ; y ) L u (w ; u ; y )(u u ) 8u U ad (variational inequality). Subscripts denote associated partial derivatives. In detail, the adjoint equation reads 4y y ' w (w ) b w(w ; u )y w (w ; u ) on : () The variational inequality admitts the form ( u (w ; u ) b u (w ; u )y)(u u )ds x 8u U ad : () To simplify the notation we denote in the sequel the pair (w; u) by v. Denition The pair v (w; u) is called an admissible pair, if u belongs to U ad and w is the weak solution of (), (3) corresponding to u. Two norms of v corresponding to the spaces V W () L () and V Wp () L () are given by kvk (kwk W () kuk L ()) and kvk maxfkwk W p (); kuk L()g; respectively. The simplest second order sucient optimality condition which can be transferred to our problem from the parabolic case can be formulated as follows: Assume that for a pair v (w ; u ) the rst order necessary conditions are satised. Let y be the associated adjoint state. Suppose that the second order derivative of the Lagrange{function with respect to w and u fulls the estimate L vv (v ; y )[(h; ); (h; )] k(h; )k (3) for all (h; ), where uu with u U ad and h is the weak solution of the corresponding linearized equation 4h h b w(v )h b u (v ) on : (5) In this condition, the second order derivative of the Lagrange{function has to full the estimate with respect to all admissible controls u. From the theory of optimization in spaces of nite dimensions weaker second order sucient conditions are known. More precisely, the estimate of the type (3) has to be fullled only with respect to all u from some subset of the admissible set. Using ideas of Dontchev et al. [4] and Malanowski [] we will derive such a sucient condition in the next sections. 4

5 3 Motivation of a second order sucient optimality condition Let us consider the mathematical program to minimize subject to x IR k f(x) h(x) g(x) ; wheref, h, and g are smooth functions. For this problem the following second order sucient condition is known (cf., for instance, Spellucci [3]): Theorem Let x be a admissible point satisfying the rst order necessary optimality condition with associated Lagrange{multipliers and. The point x is a strong local minimizer, if an > exists such that z r f(x ) mx i for all z satisfying the conditions i r g i (x ) kx j j r h j (x ) A z z T z (6) and z T rh i (x ) ; (7) z T rg i (x ) ; if g i (x ) ; (8) z T rg i (x ) ; if i > : (9) Compairing this condition with the sucient condition (3) for the elliptic control problem we observe: The estimate (6) corresponds to the estimate (3). The equation (7) can be viewed as a representation of the partial dierential equation in the sucient condition for the control problem. The inequality (8) means that ~x x z is admissible with respect to the linearized inequality constraints for suciently small. In the context of the control problem this corresponds to restrictions on in the sucient condition for the elliptic control problem. However, there is no corresponding term for the condition (9). It is this additional condition, which we shall add to the sucient condition. Usually, necessary optimality conditions for optimal control problems have to be satised for all u from some control set U ad. To simplify the presentation, let this set be described by U ad fu L () j u(x) a.e. on g: Introducing formally a Lagrange{multiplier function with respect to the inequality constraints, the variational inequality reads ( u (w ; u ) b u (w ; u )y )(u u )ds x 5

6 for all u from the whole space L (). Thus, the left factor under the integral sign has to vanish identically on. In this way, we can identify with u (w ; u ) b u (w ; u )y in some sense. The formal multiplier function plays in this case the same rule as the vector in theorem. A suitable interpretation of the condition (9) on z is now given in the case of a optimal control problem by u(x) u (x) if u (w (x); u (x)) b u (w (x); u (x))y (x) > : In the next section we will modify this intuitive additional condition and prove the suciency of the corresponding second{order condition. 4 Second{order sucient optimality condition In what follows let v (w ; u ) be admissible for the optimal control problem. We do not assume that v achieves the global minimum of the optimal control problem. Suppose that the rst order necessary optimality conditions are fullled at v with the associated adjoint state y. This is the optimality system: (variational inequality), (adjoint equation), and ( u (v ) b u (v )y )(u u )ds x 8u U ad () 4y y ' w (v b w(v )y w (v ) on 4w b(w ; u ) on (state equation). It is well known that () holds i f u (v (x)) b u (v (x))y (x)g(u(x) u (x)) almost everywhere on for all u() U ad. We dene a set of positivity: Denition 3 For " > the set of positivity " denotes the measurable subset of such that for all u U ad holds for almost all x ". f u (v (x)) b u (v (x))y (x)g(u(x) u (x)) "ju(x) u (x)j () Remark: Owing to the assumptions on, the adjoint state y is continuous on. Thus, the function u (v (x))b u (v (x))y (x) is measurable on, and the denition make sense. It is possible that the set of positivity is of measure zero. 6

7 Corollary The estimate " f u (v (x)) b u (v (x))y (x)g(u(x) u (x))ds x "ku u k L (") holds true for all u U ad. Before we start to discuss second order sucient optimality conditions we derive some useful technical results. At rst we remind of the continuity of the adjoint state y. Therefore, the assumptions on the functions standing in the objective and in the partial dierential equations allows us to prove the following statements. Lemma The second order derivative of the Lagrange{function at (v ; y ) with respect to v fulls the estimates and jl vv (v ; y )[(h ; ); (h ; )]j ck(h ; )k k(h ; )k () jl vv (v ; y )[(h; ); (h; )] L vv (v; y )[(h; ); (h; )]j ckv v k k(h; )k (3) for all (h; ); (h i ; i ) W () L () and all v W p () L ().. Proof: The second order derivative of the Lagrange{function is given by: L vv (v ; y )[(h ; ); (h ; )] ' ww (w )h h dx ( ww (v )h h wu (v )h uw (v ) h uu (v ) )ds x (b ww (v )h h b wu (v )h b uw (v ) h b uu (v ) ) y ds x : The control u is bounded in the sense of L (), and the state w is bounded in the sense of C(). Therefore, all second order derivatives under the integral sign are uniformly bounded in the sense of L (). The adjoint state y is continuous on. For that reason we are able to estimate jl vv (v ; y )[(h ; ); (h ; )]j To prove the second estimate we consider h h dx (h h h h )ds x (h h h h )ds x A ck(h ; )k k(h ; )k : L vv (v ; y )[(h; ); (h; )] L vv (v; y )[(h; ); (h; )] ' ww (w )h dx ( ww (v )h wu (v )h uu (v ) )ds x (b ww (v )h b wu (v )h b uu (v ) )y ds x 7

8 ' ww (w)h dx ( ww (v)h wu (v)h uu (v) )ds x (b ww (v)h b wu (v)h b uu (v) )y ds x (' ww (w ) ' ww (w))h dx (( ww (v ) ww (v))h ( wu (v ) wu (v))h ( uu (v ) uu (v)) )ds x ((b ww (v ) b ww (v))h (b wu (v ) b wu (v))h (b uu (v ) b uu (v)) )y ds x The second order derivatives of ',, and b are assumed to be Lipschitz{continuous. Therefore, jl vv (v ; y )[(h; ); (h; )] L vv (v; y )[(h; ); (h; )]j jw wjh dx (jv vjh jv vj jhj jj jv vj )ds x (jv vjh jv vj jhj jj jv vj )y ds x A Now the second estimate follows immediately. The next results are concerned with properties of remainder terms. We will use the following notations: r b (v ; v) b(v) b(v ) b w (v )(w w ) b u (v )(u u ) (rst order remainder term) r L (v ; y ; v) L(v; y ) L(v ; y ) L v (v v ) L vv(v ; y )[v v ; v v ] Lemma 3 The estimate is valid for all v W p () L (). (second order remainder term) kr b (v ; v)k L () kv vk ckv v k (4) Proof: We have kr(v b ; v)k L () (b(v) b(v ) b v (v )(v v )) b v (v s(v v ))(v v )ds b v (v )(v v ) (b v (v s(v v )) b v (v )))(v v )ds 8 A A ds x ds x

9 @ cjs(v v )jjv v jds A by the Lipschitz property of b v c jv v j 4 ds x ckv v k C()L() ckv v k jv v j ds x : ds x jv v j ds x So we nish with kr b (v ; ~v)k ck~v v k k~v v k : Lemma 4 The second order remainder term r L (v ; y ; v) satises the estimate Proof: We have r L (v ; y ; v) jr L (v ; y ; v)j kv v k ckv v k : (5) L(v; y ) L(v ; y ) L v (v ; y )(v v ) L vv(v ; y )[v v ; v v ] L v (v s(v v ); y )(v v )ds L v (v ; y )(v v ) L vv(v ; y )[v v ; v v s L vv (v (v v ); y )[v v ; v v ]d L(v ; y )(v v ) A ds L v (v ; y )(v v ) L vv(v ; y )[v v ; v v ] s L vv (v (v v ); y )[v v ; v v ]d ds L vv(v ; y )[v v ; v v ] s L vv (v (v v ); y )[v v ; v v ] L vv (v ; y )[v v ; v v ]d ds Invoking the previous lemma we conclude jr L (v ; y ; v)j c s kv v k kv v k d ds ckv v k kv v k : 9

10 Now we are able to establish a strenghened second order sucient optimality condition. Theorem Let v (w ; u ) be admissible for the optimal control problem and full the rst order necessary optimality condition with associated adjoint state y. Choose " > and the corresponding set of positivity ". Suppose the existence of a > such that L vv (v ; y )[(h; z); (h; z)] k(h; z)k (6) holds for all z u u, where u U ad and u(x) u (x) for almost all x ", and h is the solution of the linearized equation Then there exist > and > such that for all v U ad with kv v k : 4h h b w(v )h b u (v )z on : (8) F (v) F (v ) kv v k (9) Remark: The solution h of the linearized equation (7), (8) exists, as b w (v) is nonpositive. Proof: Let v (u; w) be an admissible pair. Suppose that v fullls the assumptions of the theorem. We have F (v) L(v; y ) since v is admissible. By means of a Taylor{expansion of the Lagrange{function at v with respect to v we get L(v; y ) L(v ; y ) L u (v ; y )(u u ) L w (v ; y )(w w )) The admissibility of v yields L vv(v ; y )[v v ; v v ] r L (v ; y ; v): L(v ; y ) F (v ): Moreover, the necessary optimality conditions are fullled at v with adjoint state y. Thus, the third term is equal to zero. On the other hand, the second term is nonnegative due to the variational inequality (). However, we get even more. The denition of " yields: L u (v ; y )(u u ) ( u (v ) b u (v )y )(u u )ds x ( u (v ) b u (v )y )(u u )ds x " "ku u k L ("):

11 Continuing, we estimate F (v) F (v ) L u (v ; y )(u u ) L w (v ; y )(w w )) L vv(v ; y )[v v ; v v ] r L (v ; y ; v) F (v ) "ku u k L (") r L (v ; y ; v) L vv(v ; y )[v v ; v v ]: Next, we discuss the second derivative of the Lagrange{function. exploit (6). To this aim, a new control ~u is introduced by Here, we intend to ~u(x) ( u (x) on " u(x) on n " Let ~w be the state corresponding to the control ~u, i. e. b(w; ~u): Let ~v denote the pair ( ~w; ~u). Inserting these notations, we obtain by means of (): L vv (v ; y )[v v ; v v ] L vv (v ; y )[v ~v ~v v ; v ~v ~v v ] L vv (v ; y )[v ~v; v ~v] L vv (v ; y )[~v v ; ~v v ] L vv (v ; y )[~v v ; v ~v] L vv (v ; y )[~v v ; ~v v ] ck~v vk ~ck~v vk k~v v k : Moreover, we introduce w t as the solution of the linearized state equation at v : 4w t w b(v ) b w (v )(w t w ) b u (v )(~u u ): Setting v t (w t ; ~u) and invoking (), we continue by L vv (v ; y )[~v v ; ~v v ] L vv (v ; y )[~v v t v t v ; ~v v t v t v ] L vv (v ; y )[~v v t ; ~v v t ] L vv (v ; y )[v t v ; v t v ] L vv (v ; y )[v t v ; ~v v t ] L vv (v ; y )[v t v ; v t v ] ck~v v t k ~ck~v v t k kv t v k : (3) For the rst term of (3) the estimate (6) applies. To handle the second and the third term we consider rst the dierence ~v v t. By denition we have k~v v t k (k ~w w t k W k~u () ~uk L () k ~w w t k W (), and ~w w t fulls the equation 4( ~w w t ) ( ~w w t ~w w t b( ~w; ~u) b(w ; u ) b w (w ; u )(w t w ) b u (w ; u )(~u u ):

12 For the right hand side of the boundary condition we get b( ~w; ~u) b(w ; u ) b w (w ; u )(w t w ) b u (w ; u )(~u u ) b( ~w; ~u) b(w ; u ) b w (w ; u )(w t ~w ~w w ) b u (w ; u )(~u u ) b( ~w; ~u) b(w ; u ) b w (w ; u )( ~w w ) b u (w ; u )(~u u ) b w (w ; u )( ~w w t ) r b (v ; ~v) b w (w ; u )( ~w w t ): Thus the boundary condition is equivalent to Lemma ~w w t b w (w ; u )( ~w w t ) r b (v ; ~v): k~v v t k k ~w w t k W () ckr b (v ; ~v)k : (3) Inserting (3) we continue the estimation of (3) by L vv (v ; y )[~v v ; ~v v ] L vv (v ; y )[v t v ; v t v ] ck~v v t k ~ck~v v t k kv t v k kv t v k ckr b (v ; ~v)k ~ckr b (v ; ~v)k kv t v k kv t ~v ~v v k ckr b (v ; ~v)k ~ckr b (v ; ~v)k kv t ~v ~v v k k~v v k kv t ~vk k~v v k ckr b (v ; ~v)k ~ckr b (v ; ~v)k kv t ~vk ~ckr b (v ; ~v)k k~v v k k~v v k ckr(v b ; ~v)k ~ckr(v b ; ~v)k k~v v k! k~v v k ckrb(v ; ~v)k ~ckrb(v ; ~v)k : k~v v k k~v v k By (4), L vv (v ; y )[~v v ; ~v v ] k~v v k ck~v v k ~ck~v v k If v is suciently close to v with respect to the L {norm, then k~v v k is suciently close to zero, too. Therefore, the term in the brackets is greater or equal to ~ > for all admissible v such that kv v k. In this way we conclude that L vv (v ; y )[~v v ; ~v v ] ~k~v v k (3) holds for all admissible v with kv v k and all corresponding pairs ~v. Further, L vv (v ; y )[v v ; v v ] L vv (v ; y )[~v v ; ~v v ] ck~v vk ~ck~v vk k~v v k ~k~v v k ck~v vk ~ck~v vk k~v v k ~k~v v v v k ck~v vk ~ck~v vk k~v v v v k ~kv v k ~kv v k k~v vk ck~v vk ~ck~v vk kv v k ~kv v k ck~v vk ~ck~v vk kv v k

13 Now we are able to complete the estimation of the cost functional. We have F (v) F (v ) "ku u k L (") r L (v ; y ; v) L vv(v ; y )[v v ; v v ] F (v ) "ku u k L (") r L (v ; y ; v) ~kv v k ck~v vk ~ck~v vk kv v k F (v ) "ku u k L (") r L (v ; y ; v) ~kv v k ck~v vk ~ck~v vk kv v k : The next term under consideration is k~v vk kv v k. Youngs inequality implies k~v vk kv v k k~v vk kv v k for all >. Moreover, we have by construction k~v vk ck~u uk ; where the term on the right hand side is uniformly bounded. estimation with Thus, we continue the F (v) F (v ) "ku u k L (") r L (v ; y ; v) ~kv v k ck~v vk ~ck~v vk kv v k F (v ) "ku u k L (") r L (v ; y ; v) ~kv v k ck~u uk ~c k~u uk kv v k F (v ) "ku u k L (") c k~u uk! ~ kv v k ~c jrl (v ; y ; v)j kv v k The general assumptions imply by (5) that jr L (v ; y ; v)j ckv v k kv v k : Therefore, for kv v k 3 and suciently small, the term in the brackets is greater or equal to > for all v with kv v k 3. Hence F (v) F (v ) "ku u k L (") c k~u uk! ~ kv v k ~c jrl(v ; y ; v)j kv v k F (v ) "ku u k L (") c k~u uk kv v k ~u diers from u only on ", where ~u is equal to u. Therefore the norm k~u uk is equal to ku uk L (") 3ku u k L ("). Now we take 3 suciently small such that 3

14 " c( ) 3 >. Then we are allowed to omit the terms with ku u k L (") and arrive at F (v) F (v ) kv v k (33) for all admissible v with kv v k 3. Re{formulating this result we arrive at Theorem 3 Suppose that the assumptions of Theorem hold true. Then u is a locally optimal control in the sense of L (). 5 Remarks The results of the proceding sections ensure local optimality of the control u in a suciently small neighbourhood in L (). If jumps of u cannot be excluded, all functions in this neighbourhood must have jumps at the same position. This is too strong for many applications. Aiming to weaken this, we consider now the optimal control problem under stronger assumptions. We shall show that in this case the assumptions of Theorem are sucient for local optimality in the sense of L p (), where p < is suciently large (cf. section ). To do so we suppose the following additional properties of and b : The functions and b full the assumptions (A), (A), and (A3) from section and have the form: (A4) b(w; u) b (w) b (w)u (34) (w; u) (w) (w)u u ; (35) where L () with (x) almost everywhere on. We assume that the functions b i and i and all of their derivatives up to the second order are Lipschitz continuous. The crucial point in the proof of Theorem were the estimates of the rst and second order remainder terms of b and the Lagrange{function L. In general, these estimates are valid only with respect to the norm in W p () L (). For that reason Theorem 3 states only local optimality in the sense of L (). The additional assumptions (34) and (35) allow us to improve the estimates of the remainder terms. We set for convenience kvk p p k(w; u)k p p kwk p W p () kukp L p() : Lemma 5 The rst order remainder term of b fullls the estimate for all v (w; u) W p () L (). kr b (v ; v)k L () ckv v k kv v k p (36) Proof: We have kr b (v ; v)k L () (b(v) b(v ) b v (v )(v v )) ds fb v (v s(v v ) b v (v )g(v v )ds A ds x : (37) 4

15 Making use of the special structure of b we can rewrite the integrand at almost all x in the following way jb v (v (x) s(v(x) v (x)) b v (v (x))g(v(x) v (x))j jfb w (w (x) s(w(x) w (x))) b w (w (x))g(w(x) w (x)) fb w (w (x) s(w(x) w (x)))(u (x) s(u(x) u (x))) b w (w (x))u (x)g (w(x) w (x)) fb (w (x) s(w(x) w (x))) b (w (x))g(u(x) u (x)))j: The functions b and b and the corresponding derivatives are assumed to be Lipschitz{ continuous. Using u L (), we can estimate jb v (v (x) s(v(x) v (x)) b v (v (x))g(v(x) v (x))j c (js(w(x) w (x))(w(x) w (x))j sj(u(x) u (x))(w(x) w (x))j ju (x)jjs(w(x) w (x))(w(x) w (x))j) c (js(w(x) w (x))(w(x) w (x))j sj(u(x) u (x))(w(x) w (x))j) : Inserting this estimate in (37) we get kr b (v ; v)k L () c c fjs(w(x) w (x))(w(x) w (x))j sj(u(x) u (x))(w(x) w (x))jgds fj(w(x) w (x))j j(u(x) u (x))(w(x) w (x))jg ds x (jw w j 4 jw w j 3 ju u j jw w j ju u j )ds x ckw w k C() ckw w k W p ()kv v k ckv v k p kv v k (jw w j jw w jju u j ju u j )ds x The statement of the Lemma is now a simple conclusion. To prove estimates of the second order remainder term of the Lagrange{ function we have used in the general case the Lipschitz{argument of Lemma. This general result cannot be transfered to our special case replacing the {norm by the p{norm. However, we are able to show Lemma 6 The estimate jl vv (v (vv ); y )[vv ; vv ]L vv (v ; y )[vv ; vv ]j ckvv k p kvv k (38) holds for all admissible pairs v W p () L () and (; ). Proof: By denition we have L vv (v ; y )[v v ; v v ] A ds x 5

16 ' ww (w )(w w ) dx ( ww (v )(w w ) wu (v )(w w )(u u ) uu (v )(u u ) ds x (b ww (v )(w w ) b wu (v )(w w )(u u ) b uu (v )(u u ) ) y ds x : (34) and (35) imply: L vv (v ; y )[v v ; v v ] ' ww (w )(w w ) dx (( ww (w ) ww (w )u )(w w ) w (w )(w w )(u u ) (u u ) )ds x y (b ww (w )(w w ) b ww (w )u (w w ) b wu (w )(w w )(u u ))ds x Analogously we have with w w (w w ) and u u (u u ) L vv (v ; y )[v v ; v v ] ' ww (w )(w w ) dx (( ww (w ) ww (w )u )(w w ) w (w )(w w )(u u ) (u u ) )ds x y (b ww (w )(w w ) b ww (w )u (w w ) b wu (w )(w w )(u u ))ds x Using the Lipschitz{continuity of the partial derivatives of ',, and b, we can estimate jl vv (v ; y )[v v ; v v ] L vv (v ; y )[v v ; v v ]j c c (jw w j(w w ) ju u j(w w ) jw w jj(u u )(w w )j)ds x jw w j(w w ) dx (jw w j(w w ) ju u j(w w ) jw w jj(u u )(w w )j)ds x jw w j(w w ) dx ckw w k jw w jjw w jdx ckv v k p kv v k : ((w w ) ju u jjw w j j(u u )(w w )j)ds x A 6

17 In particular, the embedding result kwk C() ckwk W p () was used. Therefore, the statement of the Lemma holds true. The next statement is an immediate consequence. Lemma 7 The second order remainder term r L (v ; y ; v) satises the estimate for all admissible v W p () L (). jr L (v ; y ; v)j ckv v k p kv v k (39) Proof: The proof is along the lines of the proof of Lemma 4. We only have to use Lemma 6 instead of Lemma for the Lipschitz{estimate. Summarizing up we arrive at Theorem 4 Let the general assumptions (A), (A), and (A3) and the assumption (A4) of this section be fullled. Suppose further that the assumptions of Theorem are satised. Then u is locally optimal in the sense of L p (). Proof: The proof of local optimality of u in the sense ofl () was mainly based on Wp (){regularity of the state w, on dierentiability of the nonlinearities in C()L (), and on the estimates of the remainder terms of b and L. Wp (){regularity is preserved for controls u L p (). Due to the strengthened assumptions on the appearance of u, the dierentiability is guaranteed in C() L p (), too. Therefore, to prove the statement of Theorem 4, we can proceed along the lines of the proof of Theorem. The only dierence is the use of estimates according to Lemma 5 and Lemma 7 instead of Lemma 3 and Lemma 4. References [] Robert A. Adams. Sobolev Spaces, volume 65 of Pure and Applied Mathematics. Academic Press Inc., 978. [] E. Casas. Pontryagin's Principle for Optimal Control Problems Governed by Semilinear Elliptic Equations, volume 8 of International Series of Numerical Mathematics, pages 97{4. Birkhauser Verlag Basel, 994. [3] E. Casas and J. Yong. Maximum principle for state-constrained optimal control problems governed by quasilinear elliptic equations. Dierential and Integral equations, 8():{8, January 995. [4] A. L. Dontchev, W. W. Hager, A. B. Poore, and B. Yang. Optimality, stability, and convergence in nonlinear control. Appl. Math. Optim., 3(3):97{36, May/June 995. [5] K. Eppler and A. Unger. Boundary control of semilinear elliptic equations { existence of optimal solutions. to appear, 995. [6] H. Goldberg and F. Troltzsch. Second{order sucient optimality conditions for a class of nonlinear parabolic boundary control problems. SIAM J. Control and Optimization, 3(4):7{5, July

18 [7] K. Groger. A W ;p {estimate for solutions to mixed boundary value problems for second order elliptic dierential equations. Math. Ann., 83:679{687, 989. [8] A. D. Ioe. Necessary and sucient conditions for a local minimum 3: Second order conditions and augmented duality. SIAM J. Control and Optimization, 7:66{88, 979. [9] D. Kinderlehrer and G. Stampacchia. An Introduction to variational inequalities and their applications. Academic Press, New York, 98. [] O. A. Ladyzhenskaya and N. N. Ural'tseva. Linear and Quasilinear Elliptic Equations. Academic Press, New York/London, 968. [] K. Malanowski. Sucient optimality conditions in optimal control. Technical report, Systems Research Institute, Polish Academy of Sciences, 994. [] H. Maurer. First and second order sucient optimality conditions in mathematical programming and optimal control. Mathematical Programming Study, 4:63{77, 98. [3] P. Spellucci. Numerische Verfahren der nichtlinearen Optimierung. ISNM Lehrbuch. Birkhauser, Basel, Boston, Berlin, 993. [4] D. Tiba. Lectures on the Optimal Control of Elliptic Equations, volume 3 of Lecture Notes. University of Jyvaskyla, Department of Mathematics,

LIPSCHITZ STABILITY OF SOLUTIONS OF LINEAR-QUADRATIC PARABOLIC CONTROL PROBLEMS WITH RESPECT TO PERTURBATIONS. FREDI TR OLTZSCH y

LIPSCHITZ STABILITY OF SOLUTIONS OF LINEAR-QUADRATIC PARABOLIC CONTROL PROBLEMS WITH RESPECT TO PERTURBATIONS. FREDI TR OLTZSCH y LIPSCHIT STABILITY OF SOLUTIONS OF LINEAR-QUADRATIC PARABOLIC CONTROL PROBLEMS WITH RESPECT TO PERTURBATIONS FREDI TR OLTSCH y Abstract. We consider a class of control problems governed by a linear parabolic

More information

SECOND ORDER SUFFICIENT OPTIMALITY CONDITIONS FOR SEMILINEAR ELLIPTIC EQUATIONS

SECOND ORDER SUFFICIENT OPTIMALITY CONDITIONS FOR SEMILINEAR ELLIPTIC EQUATIONS SECOND ORDER SUFFICIENT OPTIMALITY CONDITIONS FOR SOME STATE{CONSTRAINED CONTROL PROBLEMS OF SEMILINEAR ELLIPTIC EQUATIONS EDUARDO CASAS y, FREDI TR OLTSCH z, AND ANDREAS UNGER z Abstract. This paper deals

More information

REGULAR LAGRANGE MULTIPLIERS FOR CONTROL PROBLEMS WITH MIXED POINTWISE CONTROL-STATE CONSTRAINTS

REGULAR LAGRANGE MULTIPLIERS FOR CONTROL PROBLEMS WITH MIXED POINTWISE CONTROL-STATE CONSTRAINTS REGULAR LAGRANGE MULTIPLIERS FOR CONTROL PROBLEMS WITH MIXED POINTWISE CONTROL-STATE CONSTRAINTS fredi tröltzsch 1 Abstract. A class of quadratic optimization problems in Hilbert spaces is considered,

More information

Adaptive methods for control problems with finite-dimensional control space

Adaptive methods for control problems with finite-dimensional control space Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy

More information

quantitative information on the error caused by using the solution of the linear problem to describe the response of the elastic material on a corner

quantitative information on the error caused by using the solution of the linear problem to describe the response of the elastic material on a corner Quantitative Justication of Linearization in Nonlinear Hencky Material Problems 1 Weimin Han and Hong-ci Huang 3 Abstract. The classical linear elasticity theory is based on the assumption that the size

More information

Error estimates for the finite-element approximation of a semilinear elliptic control problem

Error estimates for the finite-element approximation of a semilinear elliptic control problem Error estimates for the finite-element approximation of a semilinear elliptic control problem by Eduardo Casas 1 and Fredi röltzsch 2 1 Departamento de Matemática Aplicada y Ciencias de la Computación

More information

Pointwise convergence rate for nonlinear conservation. Eitan Tadmor and Tao Tang

Pointwise convergence rate for nonlinear conservation. Eitan Tadmor and Tao Tang Pointwise convergence rate for nonlinear conservation laws Eitan Tadmor and Tao Tang Abstract. We introduce a new method to obtain pointwise error estimates for vanishing viscosity and nite dierence approximations

More information

Hamburger Beiträge zur Angewandten Mathematik

Hamburger Beiträge zur Angewandten Mathematik Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael

More information

SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS

SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation

More information

ON THE LAGRANGE-NEWTON-SQP METHOD FOR THE OPTIMAL CONTROL OF SEMILINEAR PARABOLIC EQUATIONS. FREDI TR OLTZSCH y

ON THE LAGRANGE-NEWTON-SQP METHOD FOR THE OPTIMAL CONTROL OF SEMILINEAR PARABOLIC EQUATIONS. FREDI TR OLTZSCH y ON THE LAGRANGE-NEWTON-SP METHOD FOR THE OPTIMAL CONTROL OF SEMILINEAR PARABOLIC EUATIONS FREDI TR OLTSCH y Abstract. A class of Lagrange-Newton-SP methods is investigated for optimal control problems

More information

A Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract

A Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract A Finite Element Method for an Ill-Posed Problem W. Lucht Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D-699 Halle, Germany Abstract For an ill-posed problem which has its origin

More information

QUASI-UNIFORMLY POSITIVE OPERATORS IN KREIN SPACE. Denitizable operators in Krein spaces have spectral properties similar to those

QUASI-UNIFORMLY POSITIVE OPERATORS IN KREIN SPACE. Denitizable operators in Krein spaces have spectral properties similar to those QUASI-UNIFORMLY POSITIVE OPERATORS IN KREIN SPACE BRANKO CURGUS and BRANKO NAJMAN Denitizable operators in Krein spaces have spectral properties similar to those of selfadjoint operators in Hilbert spaces.

More information

system of PenroseFife type Olaf Klein Weierstrass Institute for Applied Analysis and Stochastics Mohrenstrasse 39 D10117 Berlin Germany

system of PenroseFife type Olaf Klein Weierstrass Institute for Applied Analysis and Stochastics Mohrenstrasse 39 D10117 Berlin Germany A class of time discrete schemes for a phaseeld system of PenroseFife type Olaf Klein Weierstrass Institute for Applied Analysis and Stochastics Mohrenstrasse 39 D7 Berlin Germany email: klein@wias-berlin.de

More information

Math 497 R1 Winter 2018 Navier-Stokes Regularity

Math 497 R1 Winter 2018 Navier-Stokes Regularity Math 497 R Winter 208 Navier-Stokes Regularity Lecture : Sobolev Spaces and Newtonian Potentials Xinwei Yu Jan. 0, 208 Based on..2 of []. Some ne properties of Sobolev spaces, and basics of Newtonian potentials.

More information

ON TRIVIAL GRADIENT YOUNG MEASURES BAISHENG YAN Abstract. We give a condition on a closed set K of real nm matrices which ensures that any W 1 p -grad

ON TRIVIAL GRADIENT YOUNG MEASURES BAISHENG YAN Abstract. We give a condition on a closed set K of real nm matrices which ensures that any W 1 p -grad ON TRIVIAL GRAIENT YOUNG MEASURES BAISHENG YAN Abstract. We give a condition on a closed set K of real nm matrices which ensures that any W 1 p -gradient Young measure supported on K must be trivial the

More information

Lecture 2: Review of Prerequisites. Table of contents

Lecture 2: Review of Prerequisites. Table of contents Math 348 Fall 217 Lecture 2: Review of Prerequisites Disclaimer. As we have a textbook, this lecture note is for guidance and supplement only. It should not be relied on when preparing for exams. In this

More information

This method is introduced by the author in [4] in the case of the single obstacle problem (zero-obstacle). In that case it is enough to consider the v

This method is introduced by the author in [4] in the case of the single obstacle problem (zero-obstacle). In that case it is enough to consider the v Remarks on W 2;p -Solutions of Bilateral Obstacle Problems Srdjan Stojanovic Department of Mathematical Sciences University of Cincinnati Cincinnati, OH 4522-0025 November 30, 995 Abstract The well known

More information

Technische Universität Graz

Technische Universität Graz Technische Universität Graz Stability of the Laplace single layer boundary integral operator in Sobolev spaces O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2016/2 Technische

More information

Bifurcation from the rst eigenvalue of some nonlinear elliptic operators in Banach spaces

Bifurcation from the rst eigenvalue of some nonlinear elliptic operators in Banach spaces Nonlinear Analysis 42 (2000) 561 572 www.elsevier.nl/locate/na Bifurcation from the rst eigenvalue of some nonlinear elliptic operators in Banach spaces Pavel Drabek a;, Nikos M. Stavrakakis b a Department

More information

SECOND-ORDER SUFFICIENT OPTIMALITY CONDITIONS FOR THE OPTIMAL CONTROL OF NAVIER-STOKES EQUATIONS

SECOND-ORDER SUFFICIENT OPTIMALITY CONDITIONS FOR THE OPTIMAL CONTROL OF NAVIER-STOKES EQUATIONS 1 SECOND-ORDER SUFFICIENT OPTIMALITY CONDITIONS FOR THE OPTIMAL CONTROL OF NAVIER-STOKES EQUATIONS fredi tröltzsch and daniel wachsmuth 1 Abstract. In this paper sufficient optimality conditions are established

More information

We consider the application of an SQP method to an optimal control. programming problems have to be treated. The solution of the constrained

We consider the application of an SQP method to an optimal control. programming problems have to be treated. The solution of the constrained Optimal Control: Theory, Algorithms, and Applications, pp. 1-2 W. W. Hager and P. M. Pardalos, Editors c1997 Kluwer Academic Publishers B.V. On a SQP-Multigrid Technique for Nonlinear Parabolic Boundary

More information

system of equations. In particular, we give a complete characterization of the Q-superlinear

system of equations. In particular, we give a complete characterization of the Q-superlinear INEXACT NEWTON METHODS FOR SEMISMOOTH EQUATIONS WITH APPLICATIONS TO VARIATIONAL INEQUALITY PROBLEMS Francisco Facchinei 1, Andreas Fischer 2 and Christian Kanzow 3 1 Dipartimento di Informatica e Sistemistica

More information

Abstract. A front tracking method is used to construct weak solutions to

Abstract. A front tracking method is used to construct weak solutions to A Front Tracking Method for Conservation Laws with Boundary Conditions K. Hvistendahl Karlsen, K.{A. Lie, and N. H. Risebro Abstract. A front tracking method is used to construct weak solutions to scalar

More information

On second order sufficient optimality conditions for quasilinear elliptic boundary control problems

On second order sufficient optimality conditions for quasilinear elliptic boundary control problems On second order sufficient optimality conditions for quasilinear elliptic boundary control problems Vili Dhamo Technische Universität Berlin Joint work with Eduardo Casas Workshop on PDE Constrained Optimization

More information

2 Ben Schweizer As Hocking describes in [5], a discussion focused on the following question: Can one assume that the dynamic contact angle is a consta

2 Ben Schweizer As Hocking describes in [5], a discussion focused on the following question: Can one assume that the dynamic contact angle is a consta A well-posed model for dynamic contact angles Ben Schweizer February 27, 1998 Abstract We consider uid systems with a free boundary and with a point of contact of the free boundary with a solid wall. We

More information

WEAK LOWER SEMI-CONTINUITY OF THE OPTIMAL VALUE FUNCTION AND APPLICATIONS TO WORST-CASE ROBUST OPTIMAL CONTROL PROBLEMS

WEAK LOWER SEMI-CONTINUITY OF THE OPTIMAL VALUE FUNCTION AND APPLICATIONS TO WORST-CASE ROBUST OPTIMAL CONTROL PROBLEMS WEAK LOWER SEMI-CONTINUITY OF THE OPTIMAL VALUE FUNCTION AND APPLICATIONS TO WORST-CASE ROBUST OPTIMAL CONTROL PROBLEMS ROLAND HERZOG AND FRANK SCHMIDT Abstract. Sufficient conditions ensuring weak lower

More information

THE STOKES SYSTEM R.E. SHOWALTER

THE STOKES SYSTEM R.E. SHOWALTER THE STOKES SYSTEM R.E. SHOWALTER Contents 1. Stokes System 1 Stokes System 2 2. The Weak Solution of the Stokes System 3 3. The Strong Solution 4 4. The Normal Trace 6 5. The Mixed Problem 7 6. The Navier-Stokes

More information

t x 0.25

t x 0.25 Journal of ELECTRICAL ENGINEERING, VOL. 52, NO. /s, 2, 48{52 COMPARISON OF BROYDEN AND NEWTON METHODS FOR SOLVING NONLINEAR PARABOLIC EQUATIONS Ivan Cimrak If the time discretization of a nonlinear parabolic

More information

1. Introduction The nonlinear complementarity problem (NCP) is to nd a point x 2 IR n such that hx; F (x)i = ; x 2 IR n + ; F (x) 2 IRn + ; where F is

1. Introduction The nonlinear complementarity problem (NCP) is to nd a point x 2 IR n such that hx; F (x)i = ; x 2 IR n + ; F (x) 2 IRn + ; where F is New NCP-Functions and Their Properties 3 by Christian Kanzow y, Nobuo Yamashita z and Masao Fukushima z y University of Hamburg, Institute of Applied Mathematics, Bundesstrasse 55, D-2146 Hamburg, Germany,

More information

Rearrangements and polar factorisation of countably degenerate functions G.R. Burton, School of Mathematical Sciences, University of Bath, Claverton D

Rearrangements and polar factorisation of countably degenerate functions G.R. Burton, School of Mathematical Sciences, University of Bath, Claverton D Rearrangements and polar factorisation of countably degenerate functions G.R. Burton, School of Mathematical Sciences, University of Bath, Claverton Down, Bath BA2 7AY, U.K. R.J. Douglas, Isaac Newton

More information

PDE Constrained Optimization selected Proofs

PDE Constrained Optimization selected Proofs PDE Constrained Optimization selected Proofs Jeff Snider jeff@snider.com George Mason University Department of Mathematical Sciences April, 214 Outline 1 Prelim 2 Thms 3.9 3.11 3 Thm 3.12 4 Thm 3.13 5

More information

Chapter 12. Partial di erential equations Di erential operators in R n. The gradient and Jacobian. Divergence and rotation

Chapter 12. Partial di erential equations Di erential operators in R n. The gradient and Jacobian. Divergence and rotation Chapter 12 Partial di erential equations 12.1 Di erential operators in R n The gradient and Jacobian We recall the definition of the gradient of a scalar function f : R n! R, as @f grad f = rf =,..., @f

More information

Weak solutions for some quasilinear elliptic equations by the sub-supersolution method

Weak solutions for some quasilinear elliptic equations by the sub-supersolution method Nonlinear Analysis 4 (000) 995 100 www.elsevier.nl/locate/na Weak solutions for some quasilinear elliptic equations by the sub-supersolution method Manuel Delgado, Antonio Suarez Departamento de Ecuaciones

More information

Variational Formulations

Variational Formulations Chapter 2 Variational Formulations In this chapter we will derive a variational (or weak) formulation of the elliptic boundary value problem (1.4). We will discuss all fundamental theoretical results that

More information

ARCHIVUM MATHEMATICUM (BRNO) Tomus 32 (1996), 13 { 27. ON THE OSCILLATION OF AN mth ORDER PERTURBED NONLINEAR DIFFERENCE EQUATION

ARCHIVUM MATHEMATICUM (BRNO) Tomus 32 (1996), 13 { 27. ON THE OSCILLATION OF AN mth ORDER PERTURBED NONLINEAR DIFFERENCE EQUATION ARCHIVUM MATHEMATICUM (BRNO) Tomus 32 (996), 3 { 27 ON THE OSCILLATION OF AN mth ORDER PERTURBED NONLINEAR DIFFERENCE EQUATION P. J. Y. Wong and R. P. Agarwal Abstract. We oer sucient conditions for the

More information

2 BAISHENG YAN When L =,it is easily seen that the set K = coincides with the set of conformal matrices, that is, K = = fr j 0 R 2 SO(n)g: Weakly L-qu

2 BAISHENG YAN When L =,it is easily seen that the set K = coincides with the set of conformal matrices, that is, K = = fr j 0 R 2 SO(n)g: Weakly L-qu RELAXATION AND ATTAINMENT RESULTS FOR AN INTEGRAL FUNCTIONAL WITH UNBOUNDED ENERGY-WELL BAISHENG YAN Abstract. Consider functional I(u) = R jjdujn ; L det Duj dx whose energy-well consists of matrices

More information

A. RÖSCH AND F. TRÖLTZSCH

A. RÖSCH AND F. TRÖLTZSCH ON REGULARITY OF SOLUTIONS AND LAGRANGE MULTIPLIERS OF OPTIMAL CONTROL PROBLEMS FOR SEMILINEAR EQUATIONS WITH MIXED POINTWISE CONTROL-STATE CONSTRAINTS A. RÖSCH AND F. TRÖLTZSCH Abstract. A class of nonlinear

More information

PARAMETER IDENTIFICATION IN THE FREQUENCY DOMAIN. H.T. Banks and Yun Wang. Center for Research in Scientic Computation

PARAMETER IDENTIFICATION IN THE FREQUENCY DOMAIN. H.T. Banks and Yun Wang. Center for Research in Scientic Computation PARAMETER IDENTIFICATION IN THE FREQUENCY DOMAIN H.T. Banks and Yun Wang Center for Research in Scientic Computation North Carolina State University Raleigh, NC 7695-805 Revised: March 1993 Abstract In

More information

LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,

LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov, LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical

More information

1 The Stokes System. ρ + (ρv) = ρ g(x), and the conservation of momentum has the form. ρ v (λ 1 + µ 1 ) ( v) µ 1 v + p = ρ f(x) in Ω.

1 The Stokes System. ρ + (ρv) = ρ g(x), and the conservation of momentum has the form. ρ v (λ 1 + µ 1 ) ( v) µ 1 v + p = ρ f(x) in Ω. 1 The Stokes System The motion of a (possibly compressible) homogeneous fluid is described by its density ρ(x, t), pressure p(x, t) and velocity v(x, t). Assume that the fluid is barotropic, i.e., the

More information

Department of Mathematics. University of Notre Dame. Abstract. In this paper we study the following reaction-diusion equation u t =u+

Department of Mathematics. University of Notre Dame. Abstract. In this paper we study the following reaction-diusion equation u t =u+ Semilinear Parabolic Equations with Prescribed Energy by Bei Hu and Hong-Ming Yin Department of Mathematics University of Notre Dame Notre Dame, IN 46556, USA. Abstract In this paper we study the following

More information

Nonlinear Analysis. Global solution curves for boundary value problems, with linear part at resonance

Nonlinear Analysis. Global solution curves for boundary value problems, with linear part at resonance Nonlinear Analysis 71 (29) 2456 2467 Contents lists available at ScienceDirect Nonlinear Analysis journal homepage: www.elsevier.com/locate/na Global solution curves for boundary value problems, with linear

More information

Priority Program 1253

Priority Program 1253 Deutsche Forschungsgemeinschaft Priority Program 1253 Optimization with Partial Differential Equations Klaus Deckelnick and Michael Hinze A note on the approximation of elliptic control problems with bang-bang

More information

ON SCALAR CONSERVATION LAWS WITH POINT SOURCE AND STEFAN DIEHL

ON SCALAR CONSERVATION LAWS WITH POINT SOURCE AND STEFAN DIEHL ON SCALAR CONSERVATION LAWS WITH POINT SOURCE AND DISCONTINUOUS FLUX FUNCTION STEFAN DIEHL Abstract. The conservation law studied is @u(x;t) + @ F u(x; t); x @t @x = s(t)(x), where u is a f(u); x > 0 concentration,

More information

A Priori Error Analysis for Space-Time Finite Element Discretization of Parabolic Optimal Control Problems

A Priori Error Analysis for Space-Time Finite Element Discretization of Parabolic Optimal Control Problems A Priori Error Analysis for Space-Time Finite Element Discretization of Parabolic Optimal Control Problems D. Meidner and B. Vexler Abstract In this article we discuss a priori error estimates for Galerkin

More information

Necessary conditions for convergence rates of regularizations of optimal control problems

Necessary conditions for convergence rates of regularizations of optimal control problems Necessary conditions for convergence rates of regularizations of optimal control problems Daniel Wachsmuth and Gerd Wachsmuth Johann Radon Institute for Computational and Applied Mathematics RICAM), Austrian

More information

FIXED POINT METHODS IN NONLINEAR ANALYSIS

FIXED POINT METHODS IN NONLINEAR ANALYSIS FIXED POINT METHODS IN NONLINEAR ANALYSIS ZACHARY SMITH Abstract. In this paper we present a selection of fixed point theorems with applications in nonlinear analysis. We begin with the Banach fixed point

More information

2 Garrett: `A Good Spectral Theorem' 1. von Neumann algebras, density theorem The commutant of a subring S of a ring R is S 0 = fr 2 R : rs = sr; 8s 2

2 Garrett: `A Good Spectral Theorem' 1. von Neumann algebras, density theorem The commutant of a subring S of a ring R is S 0 = fr 2 R : rs = sr; 8s 2 1 A Good Spectral Theorem c1996, Paul Garrett, garrett@math.umn.edu version February 12, 1996 1 Measurable Hilbert bundles Measurable Banach bundles Direct integrals of Hilbert spaces Trivializing Hilbert

More information

Congurations of periodic orbits for equations with delayed positive feedback

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

WALTER LITTMAN BO LIU kind of problem is related to the deformation of a membrane with strings attached to several sides of the polygon. In other word

WALTER LITTMAN BO LIU kind of problem is related to the deformation of a membrane with strings attached to several sides of the polygon. In other word THE REGULARITY AND SINGULARITY OF SOLUTIONS OF CERTAIN ELLIPTIC PROBLEMS ON POLYGONAL DOMAINS Walter Littman Bo Liu School of Math., University of Minnesota, Minneapolis, MN55455 Abstract. The regularity

More information

EXISTENCE OF REGULAR LAGRANGE MULTIPLIERS FOR A NONLINEAR ELLIPTIC OPTIMAL CONTROL PROBLEM WITH POINTWISE CONTROL-STATE CONSTRAINTS

EXISTENCE OF REGULAR LAGRANGE MULTIPLIERS FOR A NONLINEAR ELLIPTIC OPTIMAL CONTROL PROBLEM WITH POINTWISE CONTROL-STATE CONSTRAINTS EXISTENCE OF REGULAR LAGRANGE MULTIPLIERS FOR A NONLINEAR ELLIPTIC OPTIMAL CONTROL PROBLEM WITH POINTWISE CONTROL-STATE CONSTRAINTS A. RÖSCH AND F. TRÖLTZSCH Abstract. A class of optimal control problems

More information

Differentiability of Implicit Functions Beyond the Implicit Function Theorem

Differentiability of Implicit Functions Beyond the Implicit Function Theorem Preprint Differentiability of Implicit Functions Beyond the Implicit Function Theorem Gerd November 19, 2012 Research Group Numerical Mathematics (Partial Differential Equations) Abstract Typically, the

More information

LECTURE 15 + C+F. = A 11 x 1x1 +2A 12 x 1x2 + A 22 x 2x2 + B 1 x 1 + B 2 x 2. xi y 2 = ~y 2 (x 1 ;x 2 ) x 2 = ~x 2 (y 1 ;y 2 1

LECTURE 15 + C+F. = A 11 x 1x1 +2A 12 x 1x2 + A 22 x 2x2 + B 1 x 1 + B 2 x 2. xi y 2 = ~y 2 (x 1 ;x 2 ) x 2 = ~x 2 (y 1 ;y 2  1 LECTURE 5 Characteristics and the Classication of Second Order Linear PDEs Let us now consider the case of a general second order linear PDE in two variables; (5.) where (5.) 0 P i;j A ij xix j + P i,

More information

MARY ANN HORN be decoupled into three wave equations. Thus, we would hope that results, analogous to those available for the wave equation, would hold

MARY ANN HORN be decoupled into three wave equations. Thus, we would hope that results, analogous to those available for the wave equation, would hold Journal of Mathematical Systems, Estimation, and Control Vol. 8, No. 2, 1998, pp. 1{11 c 1998 Birkhauser-Boston Sharp Trace Regularity for the Solutions of the Equations of Dynamic Elasticity Mary Ann

More information

Contents: 1. Minimization. 2. The theorem of Lions-Stampacchia for variational inequalities. 3. Γ -Convergence. 4. Duality mapping.

Contents: 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 information

i=1 α i. Given an m-times continuously

i=1 α i. Given an m-times continuously 1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable

More information

1 Constrained Optimization

1 Constrained Optimization 1 Constrained Optimization Let B be an M N matrix, a linear operator from the space IR N to IR M with adjoint B T. For (column vectors) x IR N and y IR M we have x B T y = Bx y. This vanishes for all y

More information

Weak Formulation of Elliptic BVP s

Weak Formulation of Elliptic BVP s Weak Formulation of Elliptic BVP s There are a large number of problems of physical interest that can be formulated in the abstract setting in which the Lax-Milgram lemma is applied to an equation expressed

More information

84 G. Schmidt (1.3) T : '! u := uj ; where u is the solution of (1.1{2). It is well known (cf. [7], [1]) that the linear operator T : H ( )! e H ( ) i

84 G. Schmidt (1.3) T : '! u := uj ; where u is the solution of (1.1{2). It is well known (cf. [7], [1]) that the linear operator T : H ( )! e H ( ) i Numer. Math. 69: 83{11 (1994) Numerische Mathematik c Springer-Verlag 1994 Electronic Edition Boundary element discretization of Poincare{Steklov operators Gunther Schmidt Institut fur Angewandte Analysis

More information

ON THE ASYMPTOTIC STABILITY IN TERMS OF TWO MEASURES FOR FUNCTIONAL DIFFERENTIAL EQUATIONS. G. Makay

ON THE ASYMPTOTIC STABILITY IN TERMS OF TWO MEASURES FOR FUNCTIONAL DIFFERENTIAL EQUATIONS. G. Makay ON THE ASYMPTOTIC STABILITY IN TERMS OF TWO MEASURES FOR FUNCTIONAL DIFFERENTIAL EQUATIONS G. Makay Student in Mathematics, University of Szeged, Szeged, H-6726, Hungary Key words and phrases: Lyapunov

More information

Obstacle problems and isotonicity

Obstacle problems and isotonicity Obstacle problems and isotonicity Thomas I. Seidman Revised version for NA-TMA: NA-D-06-00007R1+ [June 6, 2006] Abstract For variational inequalities of an abstract obstacle type, a comparison principle

More information

COMBINED EFFECTS FOR A STATIONARY PROBLEM WITH INDEFINITE NONLINEARITIES AND LACK OF COMPACTNESS

COMBINED EFFECTS FOR A STATIONARY PROBLEM WITH INDEFINITE NONLINEARITIES AND LACK OF COMPACTNESS Dynamic Systems and Applications 22 (203) 37-384 COMBINED EFFECTS FOR A STATIONARY PROBLEM WITH INDEFINITE NONLINEARITIES AND LACK OF COMPACTNESS VICENŢIU D. RĂDULESCU Simion Stoilow Mathematics Institute

More information

FLUX IDENTIFICATION IN CONSERVATION LAWS 2 It is quite natural to formulate this problem more or less like an optimal control problem: for any functio

FLUX IDENTIFICATION IN CONSERVATION LAWS 2 It is quite natural to formulate this problem more or less like an optimal control problem: for any functio CONVERGENCE RESULTS FOR THE FLUX IDENTIFICATION IN A SCALAR CONSERVATION LAW FRANCOIS JAMES y AND MAURICIO SEP ULVEDA z Abstract. Here we study an inverse problem for a quasilinear hyperbolic equation.

More information

Second Order Elliptic PDE

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

Entrance Exam, Real Analysis September 1, 2017 Solve exactly 6 out of the 8 problems

Entrance Exam, Real Analysis September 1, 2017 Solve exactly 6 out of the 8 problems September, 27 Solve exactly 6 out of the 8 problems. Prove by denition (in ɛ δ language) that f(x) = + x 2 is uniformly continuous in (, ). Is f(x) uniformly continuous in (, )? Prove your conclusion.

More information

The Erwin Schrodinger International Pasteurgasse 6/7. Institute for Mathematical Physics A-1090 Wien, Austria

The Erwin Schrodinger International Pasteurgasse 6/7. Institute for Mathematical Physics A-1090 Wien, Austria ESI The Erwin Schrodinger International Pasteurgasse 6/7 Institute for Mathematical Physics A-1090 Wien, Austria On the Point Spectrum of Dirence Schrodinger Operators Vladimir Buslaev Alexander Fedotov

More information

A DUALITY ALGORITHM FOR THE OBSTACLE PROBLEM

A DUALITY ALGORITHM FOR THE OBSTACLE PROBLEM Ann. Acad. Rom. Sci. Ser. Math. Appl. ISSN 2066-6594 Vol. 5, No. -2 / 203 A DUALITY ALGORITHM FOR THE OBSTACLE PROBLEM Diana Merluşcă Abstract We consider the obstacle problem in Sobolev spaces, of order

More information

Spurious Chaotic Solutions of Dierential. Equations. Sigitas Keras. September Department of Applied Mathematics and Theoretical Physics

Spurious Chaotic Solutions of Dierential. Equations. Sigitas Keras. September Department of Applied Mathematics and Theoretical Physics UNIVERSITY OF CAMBRIDGE Numerical Analysis Reports Spurious Chaotic Solutions of Dierential Equations Sigitas Keras DAMTP 994/NA6 September 994 Department of Applied Mathematics and Theoretical Physics

More information

Existence and Uniqueness

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

Memoirs on Dierential Equations and Mathematical Physics

Memoirs on Dierential Equations and Mathematical Physics Memoirs on Dierential Equations and Mathematical Physics Volume 20, 2000, 113{126 T. Chantladze, N. Kandelaki, A. Lomtatidze, and D. Ugulava THE P-LAPLACIAN AND CONNECTED PAIRS OF FUNCTIONS Abstract. The

More information

Methods for a Class of Convex. Functions. Stephen M. Robinson WP April 1996

Methods for a Class of Convex. Functions. Stephen M. Robinson WP April 1996 Working Paper Linear Convergence of Epsilon-Subgradient Descent Methods for a Class of Convex Functions Stephen M. Robinson WP-96-041 April 1996 IIASA International Institute for Applied Systems Analysis

More information

XIAO-CHUAN CAI AND MAKSYMILIAN DRYJA. strongly elliptic equations discretized by the nite element methods.

XIAO-CHUAN CAI AND MAKSYMILIAN DRYJA. strongly elliptic equations discretized by the nite element methods. Contemporary Mathematics Volume 00, 0000 Domain Decomposition Methods for Monotone Nonlinear Elliptic Problems XIAO-CHUAN CAI AND MAKSYMILIAN DRYJA Abstract. In this paper, we study several overlapping

More information

Techinical Proofs for Nonlinear Learning using Local Coordinate Coding

Techinical Proofs for Nonlinear Learning using Local Coordinate Coding Techinical Proofs for Nonlinear Learning using Local Coordinate Coding 1 Notations and Main Results Denition 1.1 (Lipschitz Smoothness) A function f(x) on R d is (α, β, p)-lipschitz smooth with respect

More information

Neumann problem: p u = juj N() : p 2 u; in ; u n = ; on : Robin problem: R() : Steklov problem: S() : p u = juj p 2 u; in ; p 2 u jruj n + jujp 2 u =

Neumann problem: p u = juj N() : p 2 u; in ; u n = ; on : Robin problem: R() : Steklov problem: S() : p u = juj p 2 u; in ; p 2 u jruj n + jujp 2 u = EIGENVALUE PROBLEMS FOR THE p-laplacian AN L^E DEPARTMENT OF MATHEMATICS, UNIVERSITY OF UTAH 55 SOUTH 4 EAST SALT LAKE CITY, UT 842, USA Abstract. We study nonlinear eigenvalue problems for the p-laplace

More information

on! 0, 1 and 2 In the Zienkiewicz-Zhu SPR p 1 and p 2 are obtained by solving the locally discrete least-squares p

on! 0, 1 and 2 In the Zienkiewicz-Zhu SPR p 1 and p 2 are obtained by solving the locally discrete least-squares p Analysis of a Class of Superconvergence Patch Recovery Techniques for Linear and Bilinear Finite Elements Bo Li Zhimin Zhang y Abstract Mathematical proofs are presented for the derivative superconvergence

More information

EXISTENCE OF NONTRIVIAL SOLUTIONS FOR A QUASILINEAR SCHRÖDINGER EQUATIONS WITH SIGN-CHANGING POTENTIAL

EXISTENCE OF NONTRIVIAL SOLUTIONS FOR A QUASILINEAR SCHRÖDINGER EQUATIONS WITH SIGN-CHANGING POTENTIAL Electronic Journal of Differential Equations, Vol. 2014 (2014), No. 05, pp. 1 8. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu EXISTENCE OF NONTRIVIAL

More information

K. Krumbiegel I. Neitzel A. Rösch

K. Krumbiegel I. Neitzel A. Rösch SUFFICIENT OPTIMALITY CONDITIONS FOR THE MOREAU-YOSIDA-TYPE REGULARIZATION CONCEPT APPLIED TO SEMILINEAR ELLIPTIC OPTIMAL CONTROL PROBLEMS WITH POINTWISE STATE CONSTRAINTS K. Krumbiegel I. Neitzel A. Rösch

More information

Space-time Finite Element Methods for Parabolic Evolution Problems

Space-time Finite Element Methods for Parabolic Evolution Problems Space-time Finite Element Methods for Parabolic Evolution Problems with Variable Coefficients Ulrich Langer, Martin Neumüller, Andreas Schafelner Johannes Kepler University, Linz Doctoral Program Computational

More information

OPTIMALITY CONDITIONS AND ERROR ANALYSIS OF SEMILINEAR ELLIPTIC CONTROL PROBLEMS WITH L 1 COST FUNCTIONAL

OPTIMALITY CONDITIONS AND ERROR ANALYSIS OF SEMILINEAR ELLIPTIC CONTROL PROBLEMS WITH L 1 COST FUNCTIONAL OPTIMALITY CONDITIONS AND ERROR ANALYSIS OF SEMILINEAR ELLIPTIC CONTROL PROBLEMS WITH L 1 COST FUNCTIONAL EDUARDO CASAS, ROLAND HERZOG, AND GERD WACHSMUTH Abstract. Semilinear elliptic optimal control

More information

Existence of minimizers for the pure displacement problem in nonlinear elasticity

Existence of minimizers for the pure displacement problem in nonlinear elasticity Existence of minimizers for the pure displacement problem in nonlinear elasticity Cristinel Mardare Université Pierre et Marie Curie - Paris 6, Laboratoire Jacques-Louis Lions, Paris, F-75005 France Abstract

More information

Linear Algebra. 1.1 Introduction to vectors 1.2 Lengths and dot products. January 28th, 2013 Math 301. Monday, January 28, 13

Linear Algebra. 1.1 Introduction to vectors 1.2 Lengths and dot products. January 28th, 2013 Math 301. Monday, January 28, 13 Linear Algebra 1.1 Introduction to vectors 1.2 Lengths and dot products January 28th, 2013 Math 301 Notation for linear systems 12w +4x + 23y +9z =0 2u + v +5w 2x +2y +8z =1 5u + v 6w +2x +4y z =6 8u 4v

More information

On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities

On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities Heiko Berninger, Ralf Kornhuber, and Oliver Sander FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)

More information

Pathwise Construction of Stochastic Integrals

Pathwise Construction of Stochastic Integrals Pathwise Construction of Stochastic Integrals Marcel Nutz First version: August 14, 211. This version: June 12, 212. Abstract We propose a method to construct the stochastic integral simultaneously under

More information

Realization of set functions as cut functions of graphs and hypergraphs

Realization of set functions as cut functions of graphs and hypergraphs Discrete Mathematics 226 (2001) 199 210 www.elsevier.com/locate/disc Realization of set functions as cut functions of graphs and hypergraphs Satoru Fujishige a;, Sachin B. Patkar b a Division of Systems

More information

The continuity method

The continuity method The continuity method The method of continuity is used in conjunction with a priori estimates to prove the existence of suitably regular solutions to elliptic partial differential equations. One crucial

More information

Maximum principle for the fractional diusion equations and its applications

Maximum principle for the fractional diusion equations and its applications Maximum principle for the fractional diusion equations and its applications Yuri Luchko Department of Mathematics, Physics, and Chemistry Beuth Technical University of Applied Sciences Berlin Berlin, Germany

More information

Pseudo-monotonicity and degenerate elliptic operators of second order

Pseudo-monotonicity and degenerate elliptic operators of second order 2002-Fez conference on Partial Differential Equations, Electronic Journal of Differential Equations, Conference 09, 2002, pp 9 24. http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu

More information

NONTRIVIAL SOLUTIONS TO INTEGRAL AND DIFFERENTIAL EQUATIONS

NONTRIVIAL SOLUTIONS TO INTEGRAL AND DIFFERENTIAL EQUATIONS Fixed Point Theory, Volume 9, No. 1, 28, 3-16 http://www.math.ubbcluj.ro/ nodeacj/sfptcj.html NONTRIVIAL SOLUTIONS TO INTEGRAL AND DIFFERENTIAL EQUATIONS GIOVANNI ANELLO Department of Mathematics University

More information

LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)

LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,

More information

Boundary Layer Solutions to Singularly Perturbed Problems via the Implicit Function Theorem

Boundary Layer Solutions to Singularly Perturbed Problems via the Implicit Function Theorem Boundary Layer Solutions to Singularly Perturbed Problems via the Implicit Function Theorem Oleh Omel chenko, and Lutz Recke Department of Mathematics, Humboldt University of Berlin, Unter den Linden 6,

More information

Economics 204 Fall 2011 Problem Set 2 Suggested Solutions

Economics 204 Fall 2011 Problem Set 2 Suggested Solutions Economics 24 Fall 211 Problem Set 2 Suggested Solutions 1. Determine whether the following sets are open, closed, both or neither under the topology induced by the usual metric. (Hint: think about limit

More information

EXISTENCE OF SOLUTIONS FOR CROSS CRITICAL EXPONENTIAL N-LAPLACIAN SYSTEMS

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

Lecture Notes on PDEs

Lecture Notes on PDEs Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential

More information

ARE202A, Fall Contents

ARE202A, Fall Contents ARE202A, Fall 2005 LECTURE #2: WED, NOV 6, 2005 PRINT DATE: NOVEMBER 2, 2005 (NPP2) Contents 5. Nonlinear Programming Problems and the Kuhn Tucker conditions (cont) 5.2. Necessary and sucient conditions

More information

Scientiae Mathematicae Japonicae Online, Vol. 5, (2001), Ryo Ikehata Λ and Tokio Matsuyama y

Scientiae Mathematicae Japonicae Online, Vol. 5, (2001), Ryo Ikehata Λ and Tokio Matsuyama y Scientiae Mathematicae Japonicae Online, Vol. 5, (2), 7 26 7 L 2 -BEHAVIOUR OF SOLUTIONS TO THE LINEAR HEAT AND WAVE EQUATIONS IN EXTERIOR DOMAINS Ryo Ikehata Λ and Tokio Matsuyama y Received November

More information

Optimal Control of Partial Differential Equations I+II

Optimal Control of Partial Differential Equations I+II About the lecture: Optimal Control of Partial Differential Equations I+II Prof. Dr. H. J. Pesch Winter semester 2011/12, summer semester 2012 (seminar), winter semester 2012/13 University of Bayreuth Preliminary

More information

epub WU Institutional Repository

epub WU Institutional Repository epub WU Institutional Repository Josef Leydold A Faber-Krahn-type Inequality for Regular Trees Working Paper Original Citation: Leydold, Josef (996) A Faber-Krahn-type Inequality for Regular Trees. Preprint

More information

Abstract. Previous characterizations of iss-stability are shown to generalize without change to the

Abstract. Previous characterizations of iss-stability are shown to generalize without change to the On Characterizations of Input-to-State Stability with Respect to Compact Sets Eduardo D. Sontag and Yuan Wang Department of Mathematics, Rutgers University, New Brunswick, NJ 08903, USA Department of Mathematics,

More information

GRONWALL'S INEQUALITY FOR SYSTEMS OF PARTIAL DIFFERENTIAL EQUATIONS IN TWO INDEPENDENT VARIABLES

GRONWALL'S INEQUALITY FOR SYSTEMS OF PARTIAL DIFFERENTIAL EQUATIONS IN TWO INDEPENDENT VARIABLES PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 33, Number I, May 1972 GRONWALL'S INEQUALITY FOR SYSTEMS OF PARTIAL DIFFERENTIAL EQUATIONS IN TWO INDEPENDENT VARIABLES DONALD R. SNOW Abstract.

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