BOUNDARY CONTROLLABILITY IN PROBLEMS JOHN E. LAGNESE. linear hyperbolic systems having piecewise constant coecients in

Size: px
Start display at page:

Download "BOUNDARY CONTROLLABILITY IN PROBLEMS JOHN E. LAGNESE. linear hyperbolic systems having piecewise constant coecients in"

Transcription

1 ESAIM: Control, Optimisation and Calculus of Variations URL: BOUNDARY CONTROLLABILITY IN PROBLEMS OF TRANSMISSION FOR A CLASS OF SECOND ORDER HYPERBOLIC SYSTEMS JOHN E. LAGNESE Abstract. We consider transmission problems for general second order linear hyperbolic systems having piecewise constant coecients in a bounded, open connected set with smooth boundary and controlled through the Dirichlet boundary condition. It is proved that such a system is exactly controllable in an appropriate function space provided the interfaces where the coecients have a jump discontinuity are all star-shaped with respect to one and the same point and the coecients satisfy a certain monotonicity condition.. Introduction and main results Let, be bounded, open, connected sets in R n with smooth boundaries ; and ;, respectively, such that, and set = n, whose boundary is ; =;[ ;. Let m and n be positive integers. For i k = ::: m,forj ` = ::: n and for =, let be real numbers having the symmetry = a k`ij = i k = ::: m j ` = ::: n: (.) When m = n we further assume that The = a jik` = a ij`k = i j k ` = ::: n: (.) are also assumed to satisfy the following ellipticity condition: if m 6= n then ij k` c ij ij c > 8( ij ) R mn (.3) if m = n then ij k` c ij ij c > 8( ij ) S n (.4) where S n denotes the set of all real symmetric n n matrices. In (.3), (.4) and in what follows, a repeated Roman index indicates summation with respect to that index over an appropriate range (either to m or to n, as the context suggests), unless otherwise noted. Fix T>and let W = B w. w m C A W = B w. w m C A Department of Mathematics, Georgetown University, Washington, DC 57 USA. Received by the journal April 5, 997. Accepted for publication October 3, 997. c Societe de Mathematiques Appliquees et Industrielles. Typeset by L A TEX.

2 344 JOHN E. LAGNESE be vector functions dened on the cylinders = ( T)and = ( T), respectively. Weset=; ( T) and =; ( T) and consider the following problem of transmission: ( w i ; w k `j = in w i ; (.5) w k `j = in i = ::: m where _ = =t and where a subscript following a comma indicates dierentiation with respect to the corresponding spatial variable w i = w i W = U on (.6) w k ` j = w k ` j on i = ::: m (.7) where = ( ::: n ) is the unit normal to ; pointing into (and towards the exterior of ), W j t= = _ W j t= = in, = : (.8) The function U in (.6) is viewed as a control function which ischosen to eect the evolution of the dynamics. More precisely, we consider below the reachability problem, which is to determine the set R T = f(w (T ) W (T ) _W (T ) _W (T )) : U Ug where U is an appropriate space of controls. The question of boundary controllability in problems of transmission has been considered by several authors. In particular, J.-L. Lions considered the system (.5) - (.8) in the special case of two wave equations, that is, he considered the problem ( w ; aw = in (.9) w ; a w = in where a > and is the ordinary Laplacian in R n, w = u on (.) w w = w a = a w on (.) w j t= = _w j t= = in, = : (.) For this problem the following result is proved in [5, Chapter VI]. Theorem.. Assume that ; is star-shaped with respect to some point x and set ;(x )=fx ;j (x ; x ) >g, (x )=;(x ) ( T), where is the unit outer normal to ;. Let R T = f(w (T ) w (T ) _w (T ) _w (T )) : u L () u=on n(x )g: If a >a and T>T(x ):=R(x )= p a,where R(x ) = max x jx ; x j, then R T = H V where H = L ( ) L ( ) V = f(' ' ) H ( ) H ( ): ' j ; = ' j ; = ' j ; g:

3 BOUNDARY CONTROLLABILITY IN PROBLEMS OF TRANSMISSION 345 The main purpose of this paper is prove the following generalization of Theorem. to the system (.5) - (.8). Let us dene ja j = 8 >< >: kl ij inf m 6= n 6=( ij )R mn ij ij kl ij inf m = n: 6=( ij )S n ij ij (.3) Theorem.. Assume that ; is star-shaped with respect to some point x and let ;(x ), (x ) and R(x ) be asintheorem.. Let R T =f(w (T ) W (T ) _W (T ) _W (T )):U [L ()] m U= on n(x )g: If and if 8 < : ij k` ij k` 8( ij ) R mn m6= n then R T = H V where ij k` ij k` 8( ij ) S n m= n T>T(x )= 8 >< >: R(x ) p ja j m 6= n p R(x ) p m = n: ja j H =[L ( )] m [L ( )] m V = f( ) [H ( )] m [H ( )] m : j ; = j ; = j ; g: (.4) (.5) Remark.3. Whether or not the monotonicity condition (.4) is necessary for exact controllability in dimension n is unknown, even in the case of ordinary wave equations considered in Theorem.. It is known that no such condition is necessary for exact boundary controllability of transmission problems in dimension one or, more generally, for exact controllability of - dimensional networks in R 3 [3, Chapter IV]. Remark.4. Set H ;( )=f' H ( ): 'j ; =g: The space V is a closed subspace of [H ( )] m [H ; ( )] m and V may be identied with the restriction to V of elements in the product space [(H ( )) ] m [(H ; ( )) ] m (cf. [5, Remark 5., p. 376]. Remark.5. With only minor modications in the proof, Theorem. can be extended to the situation involving and p nested open sets! :::! p with Set! i! i+ i = ::: p;! p =: =! i =! i n! i; i = ::: p: Setting ; i :=! i, one should assume that ; i is star-shaped with respect to a point x, i = ::: p;, and that the parameters, = ::: p,

4 346 JOHN E. LAGNESE associated with satisfy appropriate monotonicity conditions with respect to. Example.6. (Wave equations) Suppose that m =. In this case the system (.5) may be written ( w ; a j` w j` in w ; a j` w j` in where a j` = a`j.ifwe set A =(a j` ), the condition (.4) with m =<n becomes the requirement that the matrix A ; A be positive semidenite. The control time T (x ) is given by R(x )= p ja j,where A ja j = inf 6=R n jj : In particular, if A = a I,whereI is the identity matrix in R n,werecover Theorem.. Example.7. (Linear isotropic elastodynamics) Here n = m = 3 and the system of equations is (we do not use summation convention in this example) 8 >< >: w i ; w i ; ( + ) w i ; w i ; ( + ) 3X 3X j= j= w j ji = in w j ji = in i = 3 where is the (constant) mass density per unit of reference volume and, are the Lame parameters for the elastic material occupying the region.we assume that > and : (.6) By a change of notation we may assume that = =. Then the nonzero coecients are given by a ijij a iiii = + i = 3 = a ijji = a iijj = i 6= j = 3 and therefore for ( ij ) S 3 wehave X 3X a ijk` ij k` =( + ) ii + = 3X i= i= ii! + 3X i j= i6=j 3X i j= ii jj + ij 3X i j= 3X i j= i6=j ij: ( ij + ji ) In the same way, the monotonicity condition (.4) (with strict inequality) will be satised if and : The control time T (x )isgiven by p R(x )= p ja j =R(x )= p, which is what is expected since the two wave speeds associated with the region are p and p +.

5 BOUNDARY CONTROLLABILITY IN PROBLEMS OF TRANSMISSION 347 Let us briey recount other works on boundary controllability in problems of transmission or, more generally, in interface problems. These works are for the most part disjoint from the present paper. For one-dimensional interface problems (for example, string or beam equations on one-dimensional graphs in R 3 ), one may refer to [3] for a comprehensive analysis of the reachability problem utilizing controls supported at the nodes of the graph. These systems, however, do not generally fall within the framework of the problem (.5) - (.8). For higher dimensional systems, one may cite[], where boundary controllability in transmissions problems for the Reissner-Mindlin system of thin plate theory is considered. This system consists of two coupled second order equations with constant coecients but with lower order terms as well, so they too fall outside of the scope of the present work. There are also several works dealing with boundary controllability in problems of transmission or in interface problems for the wave equation in the singular case. If, for example, ; \ ; 6= in the problem (.9) - (.) above, the solution will generally suer a loss of regularity in a neighborhood of points belonging to this set. In this case the proof technique for Theorem. given in [5] fails that is the reason for the assumption. In the singular case, we can refer to two recent papers of S. Nicaise [6], [7], where the question of exact boundary controllability for the wave equation on two-dimensional networks in R n is treated. In these papers it is shown, among other things, that a result analogous to Theorem. may be obtained provided that, in addition to boundary controls, controls are also imposed near the singular vertices. The reader is referred to the cited papers for the precise statements of Nicaise's results. References to other works related to the topic of the present paper may be found in the bibliographies of the books and papers cited above.. A priori estimates Let H and V be as in Theorem.. For k set H k () = [H k ( )] m [H k ( )] m : Let := ( ), :=( ) be elements of V and consider the bilinear form ( ) = X = ' k ` i j dx V: This form is continuous on V V in the topology induced by H () and is symmetric by virtue of (.). If m 6= n, is also coercive onv because of the ellipticity assumption (.3) and j ; =, that is, for some k>, ( ) kkk H () 8 V: On the other hand, if m = n, because of the symmetry assumptions (.) we may write ( ) = X = " k`( )" ij ( ) dx

6 348 JOHN E. LAGNESE where " k`() = (' k ` + '` k ): If V then the functions ( ' i (x) x ' i (x)= ' i (x) x i = ::: m are in H () and therefore the following Korn's inequality holds [8, Theorem.]: X = " ij ( )" ij ( ) dx It follows from (.) and (.4) that ( ) c o X = X = " ij ( )" ij ( ) dx kkk H () ' i j ' i j dx: (.) 8 V Therefore in all cases [( )] = denes a norm on V equivalent to the one induced by H (). We retopologize V with this norm and let V denote the dual space of V. With H identied with its dual space we have the dense and continuous embeddings V H V. Let us consider the homogeneous boundary value problem ( 'i ; ' k `j = in ' i ; (.) ' k `j = in i = ::: m = on (.3) ' i = ' i ' k ` j = ' k ` j on i = ::: m (.4) j t= = _ j t= = in, = : (.5) Assume that := ( ) V, := ( ) H. According to standard variational theory (.) - (.5) has a unique solution in the following sense: there is a unique function := ( ) C([ T] V ) \ C ([ T] H) \ C ([ T] V ) such that h i V + ( ) = 8 V <t T and which satises (.5), where h i V denotes the pairing between elements V and V in the V ; V duality. Furthermore, the mapping S(t) : ( ) 7! ((t) (t)) _ is a denes a unitary group on V H, that is k(t)k V + k _ (t)k H = k k V + k k H t : (.6) Let A denote the canonical isomorphism of V onto V, so that ( ) = ha i V, 8 V, and set D(A) =f V : A Hg k k D(A) = ka k H : If ( )D(A)V, then C([ T] D(A))\C ([ T] V )\C ([ T] H) and S(t) is a unitary group on D(A) V. Because of our assumptions on the regions and, an element D(A) has the regularity

7 BOUNDARY CONTROLLABILITY IN PROBLEMS OF TRANSMISSION 349 H () \ V and A =;( ' k j` ' k j`). Moreover, satises both of the transmission conditions in (.4) in the trace sense and the norm kk D(A) is equivalent to that induced on D(A) by H (). These results are standard and may beproved by the usual methods. Lemma.. Suppose that ( ) D(A) V, and let h [W ()] n. Then the solution of (.) - (.5) satises the following identity: X = _' i (h r' i ) dx + = T (div h)[ _' i _' i ; + + (h ) + (h ) 'i ' k ` ' i p h p j dxdt 'k 'i ' k ` ' i j ] dxdt j ` d h ( ; )' i j ' k ` ; ' i 'k Proof. The proof utilizes classical multipliers. Form = and integrate by parts. One has ' i (h r' i ) dxdt = + ; ' k j ` d: (.7) (' i ; ' k `j )h r' i dxdt (.8) _' i (h r' i ) dx (div h) _' i _' i dxdt ; T (h ) _' i _' i d (.9) where = ( T) and is the unit normal vector pointing into the exterior of. In addition, ( ' k `j )h r' i dxdt = ; ' k `' i p h p j dxdt ; From the symmetry property (.)we obtain ' k `' i p h p j d ' k `' i pj h p dxdt: a ijk` ' k `' i pj = (aijk` ' k `' i j ) p

8 35 JOHN E. LAGNESE and therefore ( ' k `j )h r' i dxdt = ' k `' i p h p j d ; ; ' k `' i p h p j dxdt + We then obtain from (.8) - (.) the identity = _' i (h r' i ) dx + T + (h )( _' i _' i ; (div h)[ _' i _' i ; (h ) ' k `' i j d (div h) ' k ` ' i j dxdt: (.) ' k `' i p h p j dxdt ' k `' i j ) d+ We sum (.) over =, noting the relations ' i j = ' i j ' i j j = ' i j ' k ` ' i j ] dxdt ' k `' i p h p j d: i = ::: m j = ::: n (.) and observing that along the boundary region ; that is common to and the unit normal vector pointing out of points into.we obtain X ( = _' i (h r' i ) dx + + = T + (div h)( _' i _' i ; (h ) 'k ' k `' i p h p j dxdt ' k ` ' i j ) dxdt j ` d 'i (h )( ' k `' i j ; ' k `' i j ) d + ( ' k `' i p ; ' k `' i p )h p j d: (.) To complete the proof we need to show that the last two integrals in (.) agree with the last integral in (.7). Consider the last integral in (.). By choosing suitable extensions of and wemay assume that ' i H () for = and i = ::: n. We have ( ' k `' i p ; ' k `' i p )h p j = ' k ` j h p (' i p ; ' i p )=(h ) ' k ` j 'i ; ' i (.3)

9 on, since BOUNDARY CONTROLLABILITY IN PROBLEMS OF TRANSMISSION 35 (' i p ; ' i p )j ; = 'i It also follows from (.4) that on and ; ' i p i = ::: m p = ::: n: (.4) ' k `' i j = ' k ` ' i j + 'i ' k `' i j = 'i ' k ` ' i j ; ; ' i j ; ' i j : If the last expression is subtracted from the preceding one we obtain ' k `' i j ; ' k `' i j =( ; )' k `' i j + ' k ` j 'i ; ' i on (.5) where we have used the symmetry property (.) and the second transmission condition in (.4). Itfollows from (.3) and (.5) that (h )( ' k `' i j ; ' k `' i j ) d + ( ' k `' i p h p j ; ' k `' i p h p j ) d = We use (.4) once again to obtain on ( ; )' k `' i j (h )( ; )' k `' i j d: (.6) =( ; )' k `' i j +( ; 'i )' k ` j =( ; )' k `' i j + ' k ` j 'i ; ' k ` j 'i ; ' i ; ' i =( ; )' k `' i j + 'i (' k ` ; ' k `) j =( ; )' k `' i j + 'k ; ' k 'i Now substitute (.7) into (.6). ; ' i ; ' i ; ' i j `: (.7) Theorem.. Let ( ) V H. Then the solution of (.) - (.5) satises 'k for some constant C. 'i j ` d C(T + )(k k V + k k H ) (.8)

10 35 JOHN E. LAGNESE Proof. By density it suces to prove (.8) for ( ) D(A) V.Since ; and ; are smooth, there is a vector eld h [C ()] n such that h = on ; and h = on ; (see Komornik [, Lemma.]). We utilize this vector eld in (.7) and observe that the left side of (.7) is bounded above in absolute value by C o (k (T _ )k H +k(t )k V +k k H +k k V )+C T sup(k(t)k V +k (t)k _ H) [ T ] in view of (.6). C(T +)(k k H + k k V ) Remark.3. If m 6= n is follows from (.8) and (.3) that ' i L () i = ::: m: (.9) If m = n, then (.8) and (.4) imply that ' i j + ' j i L () i j = ::: n so, in particular, i ' i = L () for each i = ::: n. Since nx 'i j + ' nx j 'i nx ' i i = + i j i j= i= i j= ' j and the second sum is in L (), it follows that (.9) holds true also in this case. Theorem.4. Assume that ; is star-shaped with respect to some point x and suppose that the coecients satisfy (.4). Let ( ) V H. Then the solution of (.) - (.5) satises (T ; T (x ))(k k H + k k V ) (h ) 'k (x ) where T (x ) is given by (.5) and h(x) =x ; x. 'i j ` d (.) Proof. Because of (.8) it suces to prove (.) for ( ) D(A) V. We apply (.7) with h(x) =x ; x. Then h on;. On the other hand, assumption (.4) and the standing assumption (.) (or (.)) imply that ( ; )' i j ' k ` + Therefore X = _' i (h r' i ) dx T 'i + n +(; n ) ' k ` ' j dxdt + (x ) (h ) ; ' i 'k _' i _' i dxdt 'k ; ' k j ` : ' k ` ' i p h p j dxdt 'i j ` d: (.)

11 Now form BOUNDARY CONTROLLABILITY IN PROBLEMS OF TRANSMISSION 353 = X = and integrate by parts to obtain X T _' i ' i ( = ; = (' i ; (_' i _' i ; ' k `j )' i dxdt ' k `' i j ) dxdt ( ' k `' i j ; ' k `' i j ) d: (.) The right side vanishes because of the transmission conditions (.4). Multiply (.) by (n ; )= and add the result from (.) to obtain X = _' i (h r' i + n ; ' i ) dx + (_' i _' i + (x ) (h ) T ' k ` ' j ) dxdt 'k 'i ) j ` d: (.3) One has X = (_' i _' i + ' k ` ' j ) dxdt = T (k k V + k k H) (.4) and (dropping for the moment summation convention over repeated Roman indices) X mx = i= X mx _' i (h r' i + n ; ' i) dx K = i= _' i dx + 4K (h r' i +(n ; )' i ) dx for any K>. From Komornik [, Lemma 3.] we have (hr' i +(n;)' i ) dx (n;) Hence X mx = i= _' i (h r' i + n ; ' i) dx X mx K = i= ; (h)' i d+4 _' i dx + K (hr' i ) dx: (h r' i ) dx

12 354 JOHN E. LAGNESE where we have utilized (.3) and (.4). Nowchoose K = K(x )tobeaconstant such that (returning to summation convention over repeated Roman indices) X = (h r' i )(h r' i ) dx K We then obtain X = X K _' i (h r' i + n ; ' i) dx = X = ' k l ' i j dx 8 V: (_' i _' i + ' k `' i j ) dx = K (k k V + k k H): (.5) Use of (.4) and (.5) in (.3) yields the inequality (T ; K)(k k V + k k H) (h ) 'k (x ) 'i j ` d: One may choose K = 8 >< >: R(x ) p ja j m 6= n p R(x ) p m = n: ja j In fact, if m 6= n one has X = (R(x )) ja j (h r' i )(h r' i ) dx (R(x )) X = X = a ijkl ' k `' i j dx (R(x )) ja j while if m = n we have from (.) X = (h r' i )(h r' i ) dx (R(x )) (R(x )) ja j (R(x )) ja j = (R(x )) ja j X = ' i j ' i j dx X = X = X = X = a ijkl ' k `' i j dx " ij ( )" ij ( ) dx " kl ( )" ij ( ) dx " kl ( )" ij ( ) dx ' k `' i j dx:

13 Set BOUNDARY CONTROLLABILITY IN PROBLEMS OF TRANSMISSION Proof of Theorem. U = fu [L ()] m : U = onn(x )g: We rst make precise the sense in which (.5) - (.8) is solvable when U U. Proceeding formally for the moment, let D(A) and form = X = (w i ; w k `j )' i dx: If we integrate by parts and use (.) and properties of we obtain = X = (w i ' i dx ; = ; The last equation is the same as ' k `j w i dx ;(x ) ( W ) H +(W A) H = ; Let H and set = A ;. One has ;(x ) ' k ` j u i d; ' k ` j u i d; U = ;(x ) and therefore there is a B L(U H)such that ; ;(x ) Then (3.) may be written which is the same as B u. u m C A : ' k ` j u i d;: (3.) CkUk U kk D(A) = CkUk U k k H ' k ` j u k d;=(bu ) H 8 H: ( W A ; ) H +(W ) H =(BU ) H 8 H W + AW = ABU in (D(A)) (3.) where (D(A)) is the dual of D(A) with respect to H and A is an isomorphism of H onto (D(A)) dened by ha i D(A) =( A) H 8 H D(A): Since ABU L ( T (D(A)) )itfollows that given initial data W() = W V _W() = W (D(A)) (3.3) the problem (3.), (3.3) has a unique solution with W C([ T] V ), _W C([ T] (D(A)) ). Further, because of Theorem. and utilizing a \lifting theorem" of Lasiecka and Triggiani [4] we have that (W W ) H V =) W C([ T] H) \ C ([ T] V ) and, moreover, k(w _W)k C([ T ] HV ) C T [k(w W )k HV + kuk U ]:

14 356 JOHN E. LAGNESE In particular, the problem (.5) - (.8) has a unique solution in the sense of (3.) and W() = _W() = : (3.4) To complete the proof we apply the Hilbert uniqueness method [5]. Let =( ) be the solution of the problem ( 'i ; ' k `j = in ' i ; ' k `j = in i = ::: m = on ' i = ' i ' k ` j = ' k ` j on i = ::: m j t=t = V _j t=t = H: According to Theorems. and., for T>T(x ) one may dene on V H an equivalent norm kk F by " k( 'k = 'i )k F = j ` d# : (x ) Now letw =(W W ) be the solution of (.5) - (.8) in which wechoose From = 8 < U = : ; on (x ) on n(x ): X = we obtain by integration by parts (w i ; =h _ W(T ) i V ; (W(T ) ) H ; which may be rewritten h( W(T _ ) ;W(T )) ( )i V H = 'k (x ) w k `j )' i dx (x ) 'i ' k ` j d 'i j ` d =k( )k F : Therefore for T > T(x ) the map : ( ) 7! ( _ W(T ) ;W(T )) is the canonical isomorphism of F onto F = V H. Hence, given any (W ;W ) V H one may choose ( ) V H such that that is, ( )=(W ;W ) W(T )=W _W(T )=W : This proves that R T = H V if T>T(x ).

15 BOUNDARY CONTROLLABILITY IN PROBLEMS OF TRANSMISSION 357 References [] V. Komornik: Exact Controllability and Stabilization: The Multiplier Method,Masson, Paris, 994. [] J. E. Lagnese: Boundary controllability in transmission problems for thin plates, in Dierential Equations, Dynamical Systems, and Control Science, K.D.Elworthy, W. N. Everitt and E. B. Lee, Eds., Lecture Notes in Pure and Applied Mathematics, vol. 5, Marcel Dekker, New York, 994, 64{658. [3] J. Lagnese, G. Leugering, G. Schmidt: Modeling, Analysis and Control of Dynamic Elastic Multi-Link Structures, Birkhauser, Boston, 994. [4] I. Lasiecka, R. Triggiani: A lifting theorem for the time regularity of solutions of abstract equations with unbounded operators and applications to hyperbolic equations, Proc. AMS 4 (988), 745{755. [5] J.-L. Lions: Contr^olabilite Exacte, Perturbations et Stabilisation de Systemes Distribues: Tome, Contr^olabilite Exacte, Coll. RMA, vol. 8, Masson, Paris, 988. [6] S. Nicaise: Boundary exact controllability of interface problems with singularities I: addition of the coecients of singularities, SIAM J. Control & Opt. 34 (996), 5{ 53. [7] S. Nicaise: Boundary exact controllability of interface problems with singularities II: addition of internal controls, SIAM J. Control & Opt. 35 (997), 585{63. [8] O. A. Oleinik, A. S. Shamaev, G. A. Yosian: Mathematical Problems in Elasticity and Homogenization, Studies in Mathematics and Its Applications, vol. 6, North-Holland, Amsterdam, 99.

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

Stabilization and Controllability for the Transmission Wave Equation

Stabilization and Controllability for the Transmission Wave Equation 1900 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 46, NO. 12, DECEMBER 2001 Stabilization Controllability for the Transmission Wave Equation Weijiu Liu Abstract In this paper, we address the problem of

More information

Bilinear spatial control of the velocity term in a Kirchhoff plate equation

Bilinear spatial control of the velocity term in a Kirchhoff plate equation Electronic Journal of Differential Equations, Vol. 1(1), No. 7, pp. 1 15. ISSN: 17-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) Bilinear spatial control

More information

Controllability of linear PDEs (I): The wave equation

Controllability 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 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

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

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

Introduction. Statement of the Problem Let be a bounded open subset of R with suciently smooth boundary, and >. In this work, we will study the exact

Introduction. Statement of the Problem Let be a bounded open subset of R with suciently smooth boundary, and >. In this work, we will study the exact Exact Controllability of a hermoelastic System with Control in the hermal Component Only George Avalos January 8, 998 Abstract In this work we give a result of exact controllability for a thermoelastic

More information

Spectrum and Exact Controllability of a Hybrid System of Elasticity.

Spectrum and Exact Controllability of a Hybrid System of Elasticity. Spectrum and Exact Controllability of a Hybrid System of Elasticity. D. Mercier, January 16, 28 Abstract We consider the exact controllability of a hybrid system consisting of an elastic beam, clamped

More information

The Mortar Wavelet Method Silvia Bertoluzza Valerie Perrier y October 29, 1999 Abstract This paper deals with the construction of wavelet approximatio

The Mortar Wavelet Method Silvia Bertoluzza Valerie Perrier y October 29, 1999 Abstract This paper deals with the construction of wavelet approximatio The Mortar Wavelet Method Silvia Bertoluzza Valerie Perrier y October 9, 1999 Abstract This paper deals with the construction of wavelet approximation spaces, in the framework of the Mortar method. We

More information

2.3 Variational form of boundary value problems

2.3 Variational form of boundary value problems 2.3. VARIATIONAL FORM OF BOUNDARY VALUE PROBLEMS 21 2.3 Variational form of boundary value problems Let X be a separable Hilbert space with an inner product (, ) and norm. We identify X with its dual X.

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

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

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

Lecture Note III: Least-Squares Method

Lecture Note III: Least-Squares Method Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,

More information

Nonlinear stabilization via a linear observability

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

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

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

More information

Strong Stabilization of the System of Linear Elasticity by a Dirichlet Boundary Feedback

Strong Stabilization of the System of Linear Elasticity by a Dirichlet Boundary Feedback To appear in IMA J. Appl. Math. Strong Stabilization of the System of Linear Elasticity by a Dirichlet Boundary Feedback Wei-Jiu Liu and Miroslav Krstić Department of AMES University of California at San

More information

Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems

Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Elliptic boundary value problems often occur as the Euler equations of variational problems the latter

More information

arxiv: v1 [math.na] 27 Jan 2016

arxiv: v1 [math.na] 27 Jan 2016 Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università

More information

Exponential stability of abstract evolution equations with time delay feedback

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

Optimal shape and position of the support for the internal exact control of a string

Optimal shape and position of the support for the internal exact control of a string Optimal shape and position of the support for the internal exact control of a string Francisco Periago Abstract In this paper, we consider the problem of optimizing the shape and position of the support

More information

QUADRATIC RATE OF CONVERGENCE FOR CURVATURE DEPENDENT SMOOTH INTERFACES: A SIMPLE PROOF 1

QUADRATIC RATE OF CONVERGENCE FOR CURVATURE DEPENDENT SMOOTH INTERFACES: A SIMPLE PROOF 1 QUADRATIC RATE OF CONVERGENCE FOR CURVATURE DEPENDENT SMOOTH INTERFACES: A SIMPLE PROOF R. H. NOCHETTO Department of Mathematics and Institute for Physical Science and Technology University of Maryland,

More information

DYNAMIC WEIGHT FUNCTIONS FOR A MOVING CRACK II. SHEAR LOADING. University of Bath. Bath BA2 7AY, U.K. University of Cambridge

DYNAMIC WEIGHT FUNCTIONS FOR A MOVING CRACK II. SHEAR LOADING. University of Bath. Bath BA2 7AY, U.K. University of Cambridge DYNAMIC WEIGHT FUNCTIONS FOR A MOVING CRACK II. SHEAR LOADING A.B. Movchan and J.R. Willis 2 School of Mathematical Sciences University of Bath Bath BA2 7AY, U.K. 2 University of Cambridge Department of

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

LEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.

LEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates. LEGENDRE POLYNOMIALS AND APPLICATIONS We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.. Legendre equation: series solutions The Legendre equation is

More information

Asymptotic Behavior of a Hyperbolic-parabolic Coupled System Arising in Fluid-structure Interaction

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

ON SOME ELLIPTIC PROBLEMS IN UNBOUNDED DOMAINS

ON SOME ELLIPTIC PROBLEMS IN UNBOUNDED DOMAINS Chin. Ann. Math.??B(?), 200?, 1 20 DOI: 10.1007/s11401-007-0001-x ON SOME ELLIPTIC PROBLEMS IN UNBOUNDED DOMAINS Michel CHIPOT Abstract We present a method allowing to obtain existence of a solution for

More information

Strong stabilization of the system of linear elasticity by a Dirichlet boundary feedback

Strong stabilization of the system of linear elasticity by a Dirichlet boundary feedback IMA Journal of Applied Mathematics (2000) 65, 109 121 Strong stabilization of the system of linear elasticity by a Dirichlet boundary feedback WEI-JIU LIU AND MIROSLAV KRSTIĆ Department of AMES, University

More information

Finite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product

Finite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )

More information

1. Subspaces A subset M of Hilbert space H is a subspace of it is closed under the operation of forming linear combinations;i.e.,

1. Subspaces A subset M of Hilbert space H is a subspace of it is closed under the operation of forming linear combinations;i.e., Abstract Hilbert Space Results We have learned a little about the Hilbert spaces L U and and we have at least defined H 1 U and the scale of Hilbert spaces H p U. Now we are going to develop additional

More information

Math 413/513 Chapter 6 (from Friedberg, Insel, & Spence)

Math 413/513 Chapter 6 (from Friedberg, Insel, & Spence) Math 413/513 Chapter 6 (from Friedberg, Insel, & Spence) David Glickenstein December 7, 2015 1 Inner product spaces In this chapter, we will only consider the elds R and C. De nition 1 Let V be a vector

More information

Stability of implicit extrapolation methods. Abstract. Multilevel methods are generally based on a splitting of the solution

Stability of implicit extrapolation methods. Abstract. Multilevel methods are generally based on a splitting of the solution Contemporary Mathematics Volume 00, 0000 Stability of implicit extrapolation methods Abstract. Multilevel methods are generally based on a splitting of the solution space associated with a nested sequence

More information

EXACT NEUMANN BOUNDARY CONTROLLABILITY FOR SECOND ORDER HYPERBOLIC EQUATIONS

EXACT NEUMANN BOUNDARY CONTROLLABILITY FOR SECOND ORDER HYPERBOLIC EQUATIONS Published in Colloq. Math. 76 998, 7-4. EXAC NEUMANN BOUNDARY CONROLLABILIY FOR SECOND ORDER HYPERBOLIC EUAIONS BY Weijiu Liu and Graham H. Williams ABSRAC Using HUM, we study the problem of exact controllability

More information

Stanford Mathematics Department Math 205A Lecture Supplement #4 Borel Regular & Radon Measures

Stanford Mathematics Department Math 205A Lecture Supplement #4 Borel Regular & Radon Measures 2 1 Borel Regular Measures We now state and prove an important regularity property of Borel regular outer measures: Stanford Mathematics Department Math 205A Lecture Supplement #4 Borel Regular & Radon

More information

Linear Algebra. Session 12

Linear Algebra. Session 12 Linear Algebra. Session 12 Dr. Marco A Roque Sol 08/01/2017 Example 12.1 Find the constant function that is the least squares fit to the following data x 0 1 2 3 f(x) 1 0 1 2 Solution c = 1 c = 0 f (x)

More information

ON THE ARITHMETIC-GEOMETRIC MEAN INEQUALITY AND ITS RELATIONSHIP TO LINEAR PROGRAMMING, BAHMAN KALANTARI

ON THE ARITHMETIC-GEOMETRIC MEAN INEQUALITY AND ITS RELATIONSHIP TO LINEAR PROGRAMMING, BAHMAN KALANTARI ON THE ARITHMETIC-GEOMETRIC MEAN INEQUALITY AND ITS RELATIONSHIP TO LINEAR PROGRAMMING, MATRIX SCALING, AND GORDAN'S THEOREM BAHMAN KALANTARI Abstract. It is a classical inequality that the minimum of

More information

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

Chapter 8 Integral Operators

Chapter 8 Integral Operators Chapter 8 Integral Operators In our development of metrics, norms, inner products, and operator theory in Chapters 1 7 we only tangentially considered topics that involved the use of Lebesgue measure,

More information

Abstract. Jacobi curves are far going generalizations of the spaces of \Jacobi

Abstract. Jacobi curves are far going generalizations of the spaces of \Jacobi Principal Invariants of Jacobi Curves Andrei Agrachev 1 and Igor Zelenko 2 1 S.I.S.S.A., Via Beirut 2-4, 34013 Trieste, Italy and Steklov Mathematical Institute, ul. Gubkina 8, 117966 Moscow, Russia; email:

More information

arxiv: v1 [math.na] 29 Feb 2016

arxiv: v1 [math.na] 29 Feb 2016 EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)

More information

Simple Examples on Rectangular Domains

Simple Examples on Rectangular Domains 84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation

More information

arxiv: v1 [math.na] 27 Jan 2016

arxiv: v1 [math.na] 27 Jan 2016 Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università

More information

Elliptic Kirchhoff equations

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

SEMIGROUP WELL-POSEDNESS OF MULTILAYER MEAD-MARKUS PLATE WITH SHEAR DAMPING

SEMIGROUP WELL-POSEDNESS OF MULTILAYER MEAD-MARKUS PLATE WITH SHEAR DAMPING SEMIGROUP WELL-POSEDNESS OF MULTILAYER MEAD-MARKUS PLATE WITH SHEAR DAMPING Scott W. Hansen Department of Mathematics Iowa State University Ames, IA 50011 shansen@iastate.edu Abstract A multilayer plate

More information

Variational Principles for Equilibrium Physical Systems

Variational Principles for Equilibrium Physical Systems Variational Principles for Equilibrium Physical Systems 1. Variational Principles One way of deriving the governing equations for a physical system is the express the relevant conservation statements and

More information

Weak Convergence Methods for Energy Minimization

Weak Convergence Methods for Energy Minimization Weak Convergence Methods for Energy Minimization Bo Li Department of Mathematics University of California, San Diego E-mail: bli@math.ucsd.edu June 3, 2007 Introduction This compact set of notes present

More information

1 Introduction We study innite-dimensional systems that can formally be written in the familiar form x 0 (t) = Ax(t) + Bu(t); y(t) = Cx(t) + Du(t); t

1 Introduction We study innite-dimensional systems that can formally be written in the familiar form x 0 (t) = Ax(t) + Bu(t); y(t) = Cx(t) + Du(t); t Quadratic Optimal Control of Regular Well-Posed Linear Systems, with Applications to Parabolic Equations Olof J. Staans Abo Akademi University Department of Mathematics FIN-20500 Abo, Finland Olof.Staans@abo.

More information

Real Analysis Qualifying Exam May 14th 2016

Real Analysis Qualifying Exam May 14th 2016 Real Analysis Qualifying Exam May 4th 26 Solve 8 out of 2 problems. () Prove the Banach contraction principle: Let T be a mapping from a complete metric space X into itself such that d(tx,ty) apple qd(x,

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

Real Analysis Notes. Thomas Goller

Real Analysis Notes. Thomas Goller Real Analysis Notes Thomas Goller September 4, 2011 Contents 1 Abstract Measure Spaces 2 1.1 Basic Definitions........................... 2 1.2 Measurable Functions........................ 2 1.3 Integration..............................

More information

The Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria

The Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria ESI The Erwin Schrodinger International Boltzmanngasse 9 Institute for Mathematical Physics A-1090 Wien, Austria Noncommutative Contact Algebras Hideki Omori Yoshiaki Maeda Naoya Miyazaki Akira Yoshioka

More information

Bilinear Stochastic Elliptic Equations

Bilinear Stochastic Elliptic Equations Bilinear Stochastic Elliptic Equations S. V. Lototsky and B. L. Rozovskii Contents 1. Introduction (209. 2. Weighted Chaos Spaces (210. 3. Abstract Elliptic Equations (212. 4. Elliptic SPDEs of the Full

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

NONLINEAR FREDHOLM ALTERNATIVE FOR THE p-laplacian UNDER NONHOMOGENEOUS NEUMANN BOUNDARY CONDITION

NONLINEAR FREDHOLM ALTERNATIVE FOR THE p-laplacian UNDER NONHOMOGENEOUS NEUMANN BOUNDARY CONDITION Electronic Journal of Differential Equations, Vol. 2016 (2016), No. 210, pp. 1 7. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu NONLINEAR FREDHOLM ALTERNATIVE FOR THE p-laplacian

More information

Introduction to Real Analysis Alternative Chapter 1

Introduction to Real Analysis Alternative Chapter 1 Christopher Heil Introduction to Real Analysis Alternative Chapter 1 A Primer on Norms and Banach Spaces Last Updated: March 10, 2018 c 2018 by Christopher Heil Chapter 1 A Primer on Norms and Banach Spaces

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

1 Vectors. Notes for Bindel, Spring 2017 Numerical Analysis (CS 4220)

1 Vectors. Notes for Bindel, Spring 2017 Numerical Analysis (CS 4220) Notes for 2017-01-30 Most of mathematics is best learned by doing. Linear algebra is no exception. You have had a previous class in which you learned the basics of linear algebra, and you will have plenty

More information

Vector Space Basics. 1 Abstract Vector Spaces. 1. (commutativity of vector addition) u + v = v + u. 2. (associativity of vector addition)

Vector Space Basics. 1 Abstract Vector Spaces. 1. (commutativity of vector addition) u + v = v + u. 2. (associativity of vector addition) Vector Space Basics (Remark: these notes are highly formal and may be a useful reference to some students however I am also posting Ray Heitmann's notes to Canvas for students interested in a direct computational

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

ON SUM OF SQUARES DECOMPOSITION FOR A BIQUADRATIC MATRIX FUNCTION

ON SUM OF SQUARES DECOMPOSITION FOR A BIQUADRATIC MATRIX FUNCTION Annales Univ. Sci. Budapest., Sect. Comp. 33 (2010) 273-284 ON SUM OF SQUARES DECOMPOSITION FOR A BIQUADRATIC MATRIX FUNCTION L. László (Budapest, Hungary) Dedicated to Professor Ferenc Schipp on his 70th

More information

Small zeros of hermitian forms over quaternion algebras. Lenny Fukshansky Claremont McKenna College & IHES (joint work with W. K.

Small zeros of hermitian forms over quaternion algebras. Lenny Fukshansky Claremont McKenna College & IHES (joint work with W. K. Small zeros of hermitian forms over quaternion algebras Lenny Fukshansky Claremont McKenna College & IHES (joint work with W. K. Chan) Institut de Mathématiques de Jussieu October 21, 2010 1 Cassels Theorem

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

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

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

The Mortar Boundary Element Method

The Mortar Boundary Element Method The Mortar Boundary Element Method A Thesis submitted for the degree of Doctor of Philosophy by Martin Healey School of Information Systems, Computing and Mathematics Brunel University March 2010 Abstract

More information

ON CERTAIN CLASSES OF CURVE SINGULARITIES WITH REDUCED TANGENT CONE

ON CERTAIN CLASSES OF CURVE SINGULARITIES WITH REDUCED TANGENT CONE ON CERTAIN CLASSES OF CURVE SINGULARITIES WITH REDUCED TANGENT CONE Alessandro De Paris Università degli studi di Napoli Federico II Dipartimento di Matematica e Applicazioni R. Caccioppoli Complesso Monte

More information

12 CHAPTER 1. PRELIMINARIES Lemma 1.3 (Cauchy-Schwarz inequality) Let (; ) be an inner product in < n. Then for all x; y 2 < n we have j(x; y)j (x; x)

12 CHAPTER 1. PRELIMINARIES Lemma 1.3 (Cauchy-Schwarz inequality) Let (; ) be an inner product in < n. Then for all x; y 2 < n we have j(x; y)j (x; x) 1.4. INNER PRODUCTS,VECTOR NORMS, AND MATRIX NORMS 11 The estimate ^ is unbiased, but E(^ 2 ) = n?1 n 2 and is thus biased. An unbiased estimate is ^ 2 = 1 (x i? ^) 2 : n? 1 In x?? we show that the linear

More information

Boundary Stabilization of Coupled Plate Equations

Boundary Stabilization of Coupled Plate Equations Palestine Journal of Mathematics Vol. 2(2) (2013), 233 242 Palestine Polytechnic University-PPU 2013 Boundary Stabilization of Coupled Plate Equations Sabeur Mansouri Communicated by Kais Ammari MSC 2010:

More information

Conservative Control Systems Described by the Schrödinger Equation

Conservative Control Systems Described by the Schrödinger Equation Conservative Control Systems Described by the Schrödinger Equation Salah E. Rebiai Abstract An important subclass of well-posed linear systems is formed by the conservative systems. A conservative system

More information

using the Hamiltonian constellations from the packing theory, i.e., the optimal sphere packing points. However, in [11] it is shown that the upper bou

using the Hamiltonian constellations from the packing theory, i.e., the optimal sphere packing points. However, in [11] it is shown that the upper bou Some 2 2 Unitary Space-Time Codes from Sphere Packing Theory with Optimal Diversity Product of Code Size 6 Haiquan Wang Genyuan Wang Xiang-Gen Xia Abstract In this correspondence, we propose some new designs

More information

An explicit nite element method for convection-dominated compressible viscous Stokes system with inow boundary

An explicit nite element method for convection-dominated compressible viscous Stokes system with inow boundary Journal of Computational and Applied Mathematics 156 (2003) 319 343 www.elsevier.com/locate/cam An explicit nite element method for convection-dominated compressible viscous Stokes system with inow boundary

More information

Finite difference method for elliptic problems: I

Finite difference method for elliptic problems: I Finite difference method for elliptic problems: I Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen

More information

A Mixed Nonconforming Finite Element for Linear Elasticity

A Mixed Nonconforming Finite Element for Linear Elasticity A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,

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

A priori error analysis of the BEM with graded meshes for the electric eld integral equation on polyhedral surfaces

A priori error analysis of the BEM with graded meshes for the electric eld integral equation on polyhedral surfaces A priori error analysis of the BEM with graded meshes for the electric eld integral equation on polyhedral surfaces A. Bespalov S. Nicaise Abstract The Galerkin boundary element discretisations of the

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

Convergence of A Galerkin Method for 2-D Discontinuous Euler Flows Jian-Guo Liu 1 Institute for Physical Science andtechnology and Department of Mathe

Convergence of A Galerkin Method for 2-D Discontinuous Euler Flows Jian-Guo Liu 1 Institute for Physical Science andtechnology and Department of Mathe Convergence of A Galerkin Method for 2-D Discontinuous Euler Flows Jian-Guo Liu 1 Institute for Physical Science andtechnology and Department of Mathematics University of Maryland College Park, MD 2742

More information

1 Math 241A-B Homework Problem List for F2015 and W2016

1 Math 241A-B Homework Problem List for F2015 and W2016 1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let

More information

Appendix A Functional Analysis

Appendix A Functional Analysis Appendix A Functional Analysis A.1 Metric Spaces, Banach Spaces, and Hilbert Spaces Definition A.1. Metric space. Let X be a set. A map d : X X R is called metric on X if for all x,y,z X it is i) d(x,y)

More information

Lecture 1 Operator spaces and their duality. David Blecher, University of Houston

Lecture 1 Operator spaces and their duality. David Blecher, University of Houston Lecture 1 Operator spaces and their duality David Blecher, University of Houston July 28, 2006 1 I. Introduction. In noncommutative analysis we replace scalar valued functions by operators. functions operators

More information

Elementary 2-Group Character Codes. Abstract. In this correspondence we describe a class of codes over GF (q),

Elementary 2-Group Character Codes. Abstract. In this correspondence we describe a class of codes over GF (q), Elementary 2-Group Character Codes Cunsheng Ding 1, David Kohel 2, and San Ling Abstract In this correspondence we describe a class of codes over GF (q), where q is a power of an odd prime. These codes

More information

Functional Limits and Continuity

Functional Limits and Continuity Chapter 4 Functional Limits and Continuity 4.1 Discussion: Examples of Dirichlet and Thomae Although it is common practice in calculus courses to discuss continuity before differentiation, historically

More information

Functional Analysis. James Emery. Edit: 8/7/15

Functional Analysis. James Emery. Edit: 8/7/15 Functional Analysis James Emery Edit: 8/7/15 Contents 0.1 Green s functions in Ordinary Differential Equations...... 2 0.2 Integral Equations........................ 2 0.2.1 Fredholm Equations...................

More information

MAT 578 FUNCTIONAL ANALYSIS EXERCISES

MAT 578 FUNCTIONAL ANALYSIS EXERCISES MAT 578 FUNCTIONAL ANALYSIS EXERCISES JOHN QUIGG Exercise 1. Prove that if A is bounded in a topological vector space, then for every neighborhood V of 0 there exists c > 0 such that tv A for all t > c.

More information

10. Smooth Varieties. 82 Andreas Gathmann

10. Smooth Varieties. 82 Andreas Gathmann 82 Andreas Gathmann 10. Smooth Varieties Let a be a point on a variety X. In the last chapter we have introduced the tangent cone C a X as a way to study X locally around a (see Construction 9.20). It

More information

The WENO Method for Non-Equidistant Meshes

The WENO Method for Non-Equidistant Meshes The WENO Method for Non-Equidistant Meshes Philip Rupp September 11, 01, Berlin Contents 1 Introduction 1.1 Settings and Conditions...................... The WENO Schemes 4.1 The Interpolation Problem.....................

More information

CIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen

CIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall

More information

Generalized Korn s inequality and conformal Killing vectors

Generalized Korn s inequality and conformal Killing vectors Generalized Korn s inequality and conformal Killing vectors arxiv:gr-qc/0505022v1 4 May 2005 Sergio Dain Max-Planck-Institut für Gravitationsphysik Am Mühlenberg 1 14476 Golm Germany February 4, 2008 Abstract

More information

Micro-support of sheaves

Micro-support of sheaves Micro-support of sheaves Vincent Humilière 17/01/14 The microlocal theory of sheaves and in particular the denition of the micro-support is due to Kashiwara and Schapira (the main reference is their book

More information

CS 468. Differential Geometry for Computer Science. Lecture 13 Tensors and Exterior Calculus

CS 468. Differential Geometry for Computer Science. Lecture 13 Tensors and Exterior Calculus CS 468 Differential Geometry for Computer Science Lecture 13 Tensors and Exterior Calculus Outline Linear and multilinear algebra with an inner product Tensor bundles over a surface Symmetric and alternating

More information

with Applications to Elasticity and Compressible Flow Daoqi Yang March 20, 1997 Abstract

with Applications to Elasticity and Compressible Flow Daoqi Yang March 20, 1997 Abstract Stabilized Schemes for Mixed Finite Element Methods with Applications to Elasticity and Compressible Flow Problems Daoqi Yang March 20, 1997 Abstract Stabilized iterative schemes for mixed nite element

More information

Linear Algebra, 4th day, Thursday 7/1/04 REU Info:

Linear Algebra, 4th day, Thursday 7/1/04 REU Info: Linear Algebra, 4th day, Thursday 7/1/04 REU 004. Info http//people.cs.uchicago.edu/laci/reu04. Instructor Laszlo Babai Scribe Nick Gurski 1 Linear maps We shall study the notion of maps between vector

More information

ABOUT SOME TYPES OF BOUNDARY VALUE PROBLEMS WITH INTERFACES

ABOUT SOME TYPES OF BOUNDARY VALUE PROBLEMS WITH INTERFACES 13 Kragujevac J. Math. 3 27) 13 26. ABOUT SOME TYPES OF BOUNDARY VALUE PROBLEMS WITH INTERFACES Boško S. Jovanović University of Belgrade, Faculty of Mathematics, Studentski trg 16, 11 Belgrade, Serbia

More information

Problem Set 9 Due: In class Tuesday, Nov. 27 Late papers will be accepted until 12:00 on Thursday (at the beginning of class).

Problem Set 9 Due: In class Tuesday, Nov. 27 Late papers will be accepted until 12:00 on Thursday (at the beginning of class). Math 3, Fall Jerry L. Kazdan Problem Set 9 Due In class Tuesday, Nov. 7 Late papers will be accepted until on Thursday (at the beginning of class).. Suppose that is an eigenvalue of an n n matrix A and

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

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

LECTURE 5, FRIDAY

LECTURE 5, FRIDAY LECTURE 5, FRIDAY 20.02.04 FRANZ LEMMERMEYER Before we start with the arithmetic of elliptic curves, let us talk a little bit about multiplicities, tangents, and singular points. 1. Tangents How do we

More information

5.3. Compactness results for Dirichlet spaces Compactness results for manifolds 39 References 39 sec:intro 1. Introduction We mean by a spectr

5.3. Compactness results for Dirichlet spaces Compactness results for manifolds 39 References 39 sec:intro 1. Introduction We mean by a spectr CONVERGENCE OF SPECTRAL STRUCTURES: A FUNCTIONAL ANALYTIC THEORY AND ITS APPLICATIONS TO SPECTRAL GEOMETRY KAZUHIRO KUWAE AND TAKASHI SHIOYA Abstract. The purpose of this paper is to present a functional

More information

Classes of Linear Operators Vol. I

Classes of Linear Operators Vol. I Classes of Linear Operators Vol. I Israel Gohberg Seymour Goldberg Marinus A. Kaashoek Birkhäuser Verlag Basel Boston Berlin TABLE OF CONTENTS VOLUME I Preface Table of Contents of Volume I Table of Contents

More information