arxiv: v3 [math.oc] 20 Dec 2018

Transcription

1 Repreenttion nd Stbility Anlyi of PDE-ODE Coupled Sytem Amritm D Schin Shivkumr Siep Weilnd Mtthew Peet riv: v3 mthoc 2 Dec 218 Deprtment of Electricl Engineering Eindhoven Univerity of Technology 56 MB The Netherlnd e-mil: {mdweilnd}@tuenl School for Engineering of Mtter Trnport nd Energy Arion Stte Univerity Tempe AZ USA e-mil: {hivk8 mpeet}@uedu Abtrct: In thi work we preent clble Liner Mtri Inequlity LMI bed frmework to verify the tbility of et of liner Prtil Differentil Eqution PDE in one ptil dimenion coupled with et of Ordinry Differentil Eqution ODE vi input-output bed interconnection Our pproch etend the newly developed tte pce repreenttion nd tbility nlyi of coupled PDE tht llow prmetriing the Lypunov function on L 2 with multiplier nd integrl opertor uing polynomil kernel of emi-eprble cl In prticulr under rbitrry well-poed boundry condition we define the liner opertor inequlitie on R n L 2 nd ct the tbility condition feibility problem contrined by LMI In thi frmework no dicretition or pproimtion i required to verify the tbility condition of PDE-ODE coupled ytem The developed lgorithm h been implemented in MATLAB the tbility of emple PDE-ODE coupled ytem re verified Keyword: Prtil Differentil Eqution Ordinry Differentil eqution Stbility Lypunov Function Liner Mtri Inequlity Sum of Squre 1 INTRODUCTION Coupled Ordinry Differentil EqutionODE-Prtil Differentil EqutionPDE re ueful to model ytem tht involve tte vrying over time well tte vrying over both pce nd time Interction mong ytem tht re governed by ODE nd PDE pper in multi-cle modeling of fluid flow for intnce in Qurteroni nd Veneini 23 nd Stinner et l 214 the flow tructure re the function of pce nd time; however they re driven by globl flow pttern tht i only function of time They lo frequently pper in thermo-mechnicl ytemeg Tng nd ie21 nd chemicl rection eg Chritofide 212 A lrge cl of thee ppliction involve liner ODE nd PDE tht re mutully coupled vi input nd output illutrted in Fig 1 ut wt PDE yt ODE vt Fig 1 Interconnection between PDE-ODE model In thi pper the liner ODE re of the following form: ẋt AtBvt 1 wt Ct For b the liner one dimenionl PDE re decribed in the following form: żt A ta 1 t 2 A 2 2 t yt C 1 t B c col tbt t 2 B 1 ut C tc b t t d b Here t R no nd t R np re the tte of the ODE nd the PDE repectively The model 1-2 re connected by the ignl colvtwt R mo R po nd colutyt R mp R pp In 1 ABC re contnt mtrice of uitble dimenion In 2 A A 1 A 2 B 1 C C b re mtri vlued function of pproprite dimenion nd C 1 i contnt mtri B c R 2np 4np i required to hve full row rnk The interconnection between 1-2 i decribed by mtri L R momp popp tht h full row rnk nd relte the interconnection ignl with the help of the following lgebric eqution: vt ut L wt 3 yt

2 In thi pper we develop clble pproch to verify the tbility of the coupled ytem in 1-2 tht doe not depend on ny dicretition or pproimtion Regrding verifying tbility mny method ued in the nlyi of Ditributed Prmeter SytemDPS rely on either ppliction-pecific nlyticl pprocheeg Nicie et l 29 or projection of the entire ytem onto finite-dimenionl vector pce hown in Rvindrn 2 D et l 218 nd the reference therein Once projected to finite-dimenionl ubpce the reduced order model i treted et of ODE nd there re mny well-tudied pproche to tet for tbility nd deign controller However dicretition technique for pproimting n infinite-dimenionl ytem with finitedimenionl model re prone to numericl intbilitie nd other trunction error Thi reult in lrge ytem mtrice which often mke them computtionlly demnding Moreover thee finite-dimenionl pproimtion of the DPS do not necerily repreent the behvior of the originl ytem to required level of ccurcy thu the propertie uch tbility robutne re not ect if the finite dimenionl pproimtion replce the originl model Another pproch frequently ued in tbilition nd control of coupled PDE-ODE ytem i the bcktepping method eg Krtic nd Smyhlyev 28 Krtic 29 nd Suto nd Krtic 21 The bcktepping technique doe not require ny dicretition or pproimtionofthemodel;howeveritdoenotllowutoutilie Lypunov theory to verify tbility Alo thee method re neither generic nor clble for generl cl of PDE- ODE coupled ytem On the other hnd uing the etenion of Lypunov theory in infinite-dimenionl pce by Dtko 197 mny ttempt were mde to ue Sum-of-Squre SOS optimition method for contructing Lypunov function for PDE Some of the notble work include but not limited to Ppchritodoulou nd Peet 26 Fridmn nd Orlov 29 nd Vlmorbid et l 216 However thee work cnnot del with generl cl of coupled PDE tht involve rection trnport nd diffuion term From our lb numerou problem relted to the nlyi of tbility nd robutne nd controller-oberver ynthei re ddreed in Meyer nd Peet 215 Ghlwt nd Peet 217 etc Recently Peet 218 Shivkumr nd Peet 218 hve developed novel SOS frmework for repreenttion tbility nd robutne nlyi for coupled PDE with rbitrry boundry condition Thee reult provide unifying frmework nd computtionlly clble method for the nlyi of lrge cl of PDE The contribution of thi pper involve etending the repreenttion nd the tbility nlyi of coupled PDE to the ce of liner PDE-ODE coupled ytem To thi end we contruct qudrtic Lypunov function of the form V P {PQQ SR 1R 2} P {PQQ SR 1R 2} i cl of elf-djoint liner opertor on R no L 2 The criterion for the tbility of the PDE-ODE coupled ytem i to erch for poitive P {PQ1Q 2SR 1R 2}tht retrict the time derivtive of the Lypunov function to be negtive We prmetrie the multiplier nd integrl opertor of P {PQ1Q 2SR 1R 2}with polynomil kernel of emieprble cl in 4 We how tht the poitivity contrint on thee opertor cn be re-formulted clble Liner Mtri InequlitieLMI tht etblih the tbility condition for PDE-ODE coupled ytem P {PQQ SR 1 R 2 } : P Qd 4 Q S R 1d R 2d The orgnition of the pper i follow Section 2 provide preliminry dicuion on the nottion model of the PDE-ODE coupled ytem nd the Lypunov tbility theory In Section 3 we define the cl of opertor P {PQ1Q 2SR 1R 2}nd derive it poitivity condition well the compoition nd djoint opertion of thee opertor In Section 4 the PDE-ODE coupled model h been repreented uing the me cl of opertor In Section 5 the equivlent LMI condition hve been derived to verify the tbility of PDE-ODE coupled ytem In Section 6 we numericlly illutrte the methodology to verify the tbility of emple PDE-ODE coupled ytem in MATLAB At lt Section 7 provide ome concluion nd direction for future reerch 21 Nottion 2 PRELIMINARIES For convenience we denote : nd 2 : 2 We ue Z to denote field of integer We ue S m R m m to denote the ymmetric mtrice We define the pce of qure integrble R m -vlued function on L m 2 L m 2 i equipped with the inner product y L 2 yd For denoting the inner product in R n we ue y R n y The Sobolov pce i defined by W qn : { L n 2 k L n 2 for ll k q} For n inner product pce k opertor P : i clled poitive if for ll we hve P We ue P to indicte tht P i poitive opertor We y tht P : i coercive if there eit ome ǫ > uch tht P ǫ 2 for ll 22 PDE-ODE Coupled Sytem In thi ection we elborte on the decription of the coupled liner PDE-ODE ytemthe ytem i repreented n interconnection between two ytem; one of them i governed by et of prtil differentil eqution PDE while the other i governed by et of ordinry differentil eqution ODE In the following ubection we briefly decribe the correponding PDE nd ODE model A ODE Model: The ODE model h the following tte pce repreenttion: ẋt AtBvt wt Ct 5 t R no re the tte vrible vt R mo wt R po re the interconnection ignl ABC re mtrice of uitble dimenion

3 B PDE Model: An opertor bed nottion cn be utilied to repreent the PDE model in the following btrct tte pce form: żt AtBut 6 yt Ct For function : b R R np u : R R mp y : R R pp the ytem opertora : L np 2 b D A L np 2 b B : L np Rmp 2 b nd C : Lnp 2 b re Rpp pecified by the following definition: A : A A 1 A 2 Bu : B 1 u C : C 1 b C C b d 7 with it domin D A :{ W 22 : B c b } b :col b b In 7 A A 1 A 2 R np np re mtri vlued function B 1 R np mp C C b R mp np re lo mtri vlued function C 1 R mp 4np i contnt mtri B c R 2np 4np h full row rnk C PDE-ODE Interconnection: A regulr interconnection between the ODE model in 5 nd the PDE model in 6 i pecified by the following lgebric contrint on the vrible pir vw nd uy: vt ut L1 L 2 L 3 L 4 }{{} L wt yt 8 9 The mtri L R momp popp h full row rnk nd i prtitioned to relte interconnection ignl of ODE nd PDE individully Subtituting 9 in 5 nd 6 we obtin the following utonomou tte-pce model ẋt ABL1 C BL 2 C żt BL 3 C ABL 4 C }{{} Ā t t 1 with Ā : L np Rno 2 b D Rno A R no L np 2 b We define : R no L np 2 b with inner product 1 1 : L np 2 b 11 nd norm Lypunov Stbility Theorem 1 2 R no 2 2 L np 12 2 b The necery nd ufficient condition for eponentil tbility of liner infinite dimenionl ytem h been provided in Curtin nd Zwrt 1995 The me cn be etended Theorem 1 tht provide condition to verify the eponentil tbility of PDE-ODE coupled ytem Theorem 1 Suppoe Ā : D Rno A generte trongly continuou emigroup nd there eit αβγ > nd P : uch tht α 2 P β 2 nd PĀ Ā P γ 2 13 for ll R no D A Then 1 i eponentilly tble 3 PARAMETRIZATION OF THE CLASS OF OPERATOR P {PQ1Q 2SR 1R 2} Definition 2 For mtri P R no no nd bounded polynomil function Q 1 : b R no np Q 2 : b R np no S : b R np np nd R 1 R 2 : b b R np np with R no nd L np 2 b the cl of opertor P {PQ1Q 2SR 1R 2}re defined P {PQ1Q 2SR 1R 2} : Q 2 S P Q 1d R 1d R 2d 14 In the reminder of thi ection we invetigte three propertieoftheclofopertorp {PQ1Q 2SR 1R 2}Firt we provide ufficient condition on the poitivity of the opertor Subequently n efficient method to compute compoition nd djoint of the opertor of the form P {PQ1Q 2SR 1R 2}hbeen providedthe compoition nd the djoint opertion on the opertor will llow u to obtin clble repreention of the PDE-ODE coupled ytem With the help of the poitivity condition on the opertor the Lypunov tbility condition cn be reformulted LMI 31 Poitivity of P {PQQ SR 1R 2} Opertor Theorem 3 For ny function Z R d1 n Z R d2 n b uppoe there eit mtri T uch tht P T 11 Q T 12 Z S Z T 22 Z T 13 Zd R 1 Z T 23 ZZ T 42 Z Z θt 33 Zθdθ Z θt 44 Zθdθ T 14 Zd Z θt 43 Zθdθ R 2 Z T 32 ZZ T 24 Z Z θt 33 Zθdθ Z θt 34 Zθdθ Z θt 44 Zθdθ 15

4 with T 11 T 12 T 13 T 14 T T 21 T 22 T 23 T 24 T 31 T 32 T 33 T 34 T 41 T 42 T 43 T 44 Then the opertor P {PQQ SR 1R 2} defined in 14 i poitive ie P {PQQ SR 1R 2} Proof Since T we cn define qure root of T U U 1 U 2 U 3 U 4 T 11 T 12 T 13 T 14 T T 21 T 22 T 23 T 24 T 31 T 32 T 33 T U U 34 T 41 T 42 T 43 T 44 U1 U 1 U1 U 2 U1 U 3 U 1 U 4 U2 U 2 U2 U 3 U2 U 4 U3 U 3 U3 U 4 U4 U 4 For 1 R m nd 2 L n 2blet udefinev : U 1 1 Ψ 2 Ψ 2 U 2 Z U 3 Z 2 d U 4 Z 2 d Then vv R m L n 2 b U 1 1 U 1 1 R m U 1 1 Ψ 2 R m Ψ 2 U 1 1 L n 2 b Ψ 2 Ψ 2 L n 2 b 1 P {PQQ SR 1R 2 2} Hence P {PQQ SR 1R 2} 32 Compoition of The Opertor 1 2 R m L n 2 b In thi ubection we how tht the compoition of two opertoroftheclp {PQ1Q 2SR 1R 2}lotketheme tructure nd provide n efficient method to clculte it Lemm 4 For ny mtrice AP R m m nd bounded function B 1 Q 1 : b R m n B 2 Q 2 : b R n m DS : b R n n C i R i : b b R n n with i {12} the following identity hold P {AB1B 2DC 1C 2}P {PQ1Q 2SR 1R 2} P {ˆP ˆQ1ˆQ 2ŜˆR 1ˆR 2} 16 ˆP AP B 1 Q 2 d ˆQ 1 AQ 1 B 1 S B 1 R 2 d ˆQ 2 B 2 P DQ 2 C 2 Q 2 d B 1 R 1 d C 1 Q 2 d Ŝ DS ˆR 1 B 2 Q 1 DR 1 C 1 S C 1 θr 2 θdθ C 2 θr 1 θdθ C 1 θr 1 θdθ ˆR 2 B 2 Q 1 DR 2 C 2 S C 1 θr 2 θdθ C 2 θr 2 θdθ C 2 θr 1 θdθ 17 Proof The proof cn be found in the Appendi C Nottion 1 For convenience we y ˆP ˆQ1 ˆQ 2 Ŝ ˆR 1 ˆR 2 AB 1 B 2 DC 1 C 2 18 PQ 1 Q 2 SR 1 R 2 19 if ˆP ˆQ1 ˆQ 2 Ŝ ˆR 1 ˆR 2 re relted to AB 1 B 2 DC 1 C 2 nd PQ 1 Q 2 SR 1 R 2 ccording to 17 in Lemm 4 33 Adjoint of The Opertor Now we lo dicu the djoint of n opertor of the cl P {PQ1Q 2SR 1R 2}nd it i typiclly denoted by P {PQ 1Q 2SR 1R 2} Lemm 5 For ny mtrice P R m m nd bounded function Q 1 : b R m n Q 2 : b R n m S : b R n n R 1 R 2 : b b R n n the following identity hold for ny R m y L n 2b P {PQ1Q 2SR 1R 2}y R m L n 2 b P{PQ1Q2SR1R2} y 2 R m L n 2 b P{PQ P 1Q 2SR 1R 2} {ˆPˆQ 1 ˆQ with 2ŜˆR 1ˆR 2} ˆP P Ŝ S ˆQ 1 Q 2 ˆR1 R2 ˆQ 2 Q 1 ˆR2 R1 21 Proof The proof cn be found in the Appendi D Nottion 2 For convenience we y ˆP ˆQ1 ˆQ 2 Ŝ ˆR 1 ˆR 2 PQ 1 Q 2 SR 1 R 2 22 if ˆP ˆQ1 ˆQ 2 Ŝ ˆR 1 ˆR 2 re relted to PQ 1 Q 2 SR 1 R 2 ccording to 21 in Lemm 5 4 FUNDAMENTAL IDENTITIES In thi ection we utilie the cl of opertor defined by P {PQ1Q 2SR 1R 2}to uniquely repreent the ytem 1 in term of the vrible nd replcing Lemm 6 Suppoe R no nd W 22 b nd B c colb b with B c R 2np 4np h full row rnk Then P {IG1G 2} 23

5 with P {IG3G 4} 24 G 1 B I 25 G 2 B 26 G 3 B b I 27 G 4 B b 28 B I IB f 29 B b IB f 3 B f B c B d 1 b I B c 31 I I I b I B d I 32 I Proof For detiled proof we refer to Peet 218 p 2-3 Lemm 7 Suppoe R no nd W 22 b nd B c colb b with B c R 2np 4np h full row rnk Then ẋ ABL1 C BL 2 C 33 ż BL 3 C ABL 4 C implie ẋ P ż {UV1V 2YW 1W 2} 34 U Λ 1 Λ 2 Y Π 1 Π 2 A 2 B 1 L 4 C 3 V 1 Ω 1 Ω 2 BL 2 C 3 W 1 Ψ 1 Ψ 2 V 2 Σ 1 Σ 2 W 2 Υ 1 Υ 2 with Λ1 Ω 1 Σ 1 Π 1 Ψ 1 Υ 1 ABL 1 CBL 2 C B 1 L 3 CA B 1 L 4 C IG 1 G 2 Λ2 Ω 2 Σ 2 Π 2 Ψ 2 Υ 2 BL 2 C b A 1 B 1 L 4 C b IG 3 G 4 b I C 3 C 1 B d B c C 1 35 I Moreover the definition of G 1 G 2 G 3 G 4 re given ccording to Lemm 6 Proof For detiled proof ee Appendi A The bove identitie imply tht the originl tte cn be fullyrecontructedbythe fundmentltte eepeet 218 Moreover the new repreenttion of the ytem dynmic in term of i independent of the boundry contrint In the ubequent ection thee identitie re utilied for deriving tbility condition 5 LYAPUNOV STABILITY CONDITION In order to verify the tbility of the PDE-ODE coupled ytem we hve propoed Lypunov function of the form V P {PQQ SR 1R 2} uch tht P {PQQ SR 1R 2} i elf-djoint ie P {PQQ SR 1R 2} P {PQQ SR 1R 2} 51 Derivtive of The Lypunov Function 36 Theorem 8 Suppoe R no nd W 22 b nd B c colb b with B c R 2np 4np h full row rnk Moreover uppoe there eit mtri P R no no nd bounded polynomil function Q : b R no np S : b R np np R 1 R 2 : b b R np np nd ABL1 C BL Ā 2 C 37 BL 3 C ABL 4 C Then Ā P {PQQ SR 1R 2}Ā P {PQQ SR 1R 2} P { ˆKˆLˆL ˆM ˆN 1 ˆN 2} ˆK K K ˆM M M ˆL L1 L 1 ˆL L 2 L 2 ˆN1 N 1 N 1 ˆN2 N 2 N 2 K L1 L 2 M N 1 N 2 KL1 L 2 MN 1 N 2 KL1 L 2 MN 1 N 2 IG1 G 2 PQQ SR 1 R 2 38 U 1 V 1 V 2 YW 1 W 2 39 Moreover G 1 G 2 re defined in Lemm 6 nd U V 1 V 2 Y W 1 W 2 re defined in Lemm 7 Proof The proof of Theorem 8 cn be eily hown by ubtituting Lemm 6-7 in Condition of Eponentil Stbility In order to verify the eponentil tbility of the PDE- ODE coupled ytem we erch for the opertor P {PQQ SR 1R 2} ubject to P { ˆKˆLˆL ˆM ˆN 1 ˆN 2} Definition 9 The cone of poitive opertor P tht re given in the Definition 2 with polynomil multiplier nd kernel ocited with degree d 1 nd d 2 re pecified by Φ d1d 2 : {PQQ SR 1 R 2 : the poitivity condition i tified ccording to Theorem 3} 4 Theorem 1 Suppoe there eit ǫ > d 1 d 2 Z with P R no nd bounded polynomil function Q : b R no np S : b R np np R 1 R 2 : b b R np np uch tht P ǫiqq S ǫir 1 R 2 Φd1d 2

6 ˆK ˆL ˆL MN 1 N 2 ǫ IT1 T 2 Φd1d 2 41 with ˆK K K ˆM M M ˆL L1 L 1 ˆL L 2 L 2 ˆN1 N 1 N 1 ˆN2 N 2 N 2 K L1 L 2 M N 1 N 2 KL1 L 2 MN 1 N 2 KL1 L 2 MN 1 N 2 IG1 G 2 PQQ SR 1 R 2 U 1 V 1 V 2 YW 1 W 2 IT1 T 2 IG1 G 2 IG 1 G 2 42 G 1 G 2 re defined in Lemm 6 nd U V 1 V 2 Y W 1 W 2 re defined in Lemm 7 Then the ytem in 1 i eponentilly tble Proof We hve defined the Lypunov function V P {PQQ SR 1R 2} P {P ǫiqq S ǫir 1R 2} ǫ 2 ǫ 2 43 Thi how the trict poitivity of P {PQQ SR 1R 2} Furthermore bed on Theorem 8 the eponentil tbility in Theorem 1 cn be reformulted ǫ P { ˆKˆLˆL ˆN 1 ˆN 2} P { ˆKˆLˆL ˆN 1 ˆN 2} ǫ 2 ǫ 2 P { ˆKˆLˆL ˆN 1 ˆN 2} P {IG1G 2} ǫ 2 P { ˆKˆL1ˆL 2 ˆM ˆN 1 ˆN 2} P {IG1G 2} ǫ P {IT1T 2} ǫ 2 ǫ Thi prove the eponentil tbility ccording to Theorem 1 2 In Theorem 3 we cn chooe Z to be vector of monomil of degree d 1 nd Z to be vector of monomil of degree d 2 Hence the contrint PQQ SR 1 R 2 Φ d1d 2 cn be viewed n LMI contrint A reult the verifiction of tbility for PDE-ODE coupled ytem mount to feibility tet of tifying the LMI contrint relted to 41 in Theorem 1 6 NUMERICAL IMPLEMENTATIONS In thi ection we illutrte the developed methodology for verifying the tbility of PDE-ODE coupled ytem The method h been implemented in MATLAB R uing n dpted verion of SOSTOOLS ee Ppchritodoulou et l ODE Coupled with Diffuion Eqution Firt we tudy the boundry controlled thermo-mechnicl proce lumped mechnicl ytem i driven by converting therml energy to mechnicl work In Tng nd ie 21 uch ytem i modeled finite dimenionl ODE model with the ctutor dynmic tht re governed by the diffuion eqution A cloed loop repreenttion of uch controlled PDE-ODE coupled ytem i given in Tng nd ie 21 Ẋt 3twt 44 w t t w t 45 w t 46 w1t 47 Bed on the definition in Section 3 we obtin A 3B 1C for the ODE ytem For the PDE model A A 1 A 2 1 The input nd output opertor re B 1 C 1 1 with C C b The interconnection mtri L 1 1 The developed methodology prove the eponentil tbility of thi ytem which h been verified nlyticlly in Tng nd ie ODE Coupled with Diffuion-Rection PDE The net emple i tken from chemicl procee due to cler ditinction between microcopic procee nd mcrocopic procee model often involve PDE governing mcrocopic dynmic coupled with ODE governing microcopic dynmic For eothermic or endothermic rection procee we conider generl cl of forced diffuion-rection in the ptil domin 1 which i repreented by the following coupled PDE żt λt tb 1 ut 48 In Vlmorbid et l 214 it h been hown tht if the boundry condition re choen to be t 1t then the unforced ytem given in 48ie B 1 i tble if nd only if λ < π Similrly in Vlmorbid et l 216 it h been hown tht for the boundry condition t 1t the ytem i tble if nd only if λ 2467 Such diffuion-rection procee in 48 re coupled with finite dimenionl ytem tht typiclly repreent

7 moleculr trnport in the microcopic cle Specificlly we conider tht the effect of thee moleculr trnport pplie uniform in-domin input to the PDE over the entire ptil domin ie B 1 1 b On the other hnd the output flu of the chemicl rection t ie t drive the moleculr trnport nd i conidered to be the input to the ODE model I The interconnection mtri L The ODE model I i given in 49 The tbility of 49 cn be verified by checking the eigenvlue ẋt vt wt t t 49 The developed methodology i ble to verify the tbility of the coupled ytem for λ < π with boundry condition t 1t nd b λ 2467 with boundry condition t 1t 7 CONCLUSIONS In thi pper we hve preented novel methodology to verify tbility for lrge cl of mutully coupled liner PDE-ODE ytem in one ptil dimenion In prticulr by uing poitive-definite mtrice to prmetrie cone of Lypunov function we hve hown tht the tbility of liner PDE-ODE coupled ytem mount to tifying pecific Liner Mtri Inequlity LMI tht cn be efficiently olved uing conve optimition We hve provided clble vrition of Sum of Squre SOS lgorithm tht cn be utilied to verify the tbility lrge cl of PDE-ODE coupled ytem efficiently In future the reult of tbility nlyi for PDE-ODE coupled ytem cn be etended to tudy the robutne propertie nd H nlyi of PDE-ODE coupled ytem Etending the current frmework for PDE in more thn one ptil dimenion nd nonliner PDE re potentil direction of future reerch Furthermore viewing the coupling of ODE nd PDE feedbck interconnection the problem of yntheiing finite dimenionl controller nd oberver re of gret prcticl importnce REFERENCES Chritofide PD 212 Nonliner nd robut control of PDE ytem: Method nd ppliction to trnportrection procee Springer Science & Buine Medi Curtin RF nd Zwrt HJ 1995 An Introduction to Infinite-dimenionl Liner Sytem Theory Springer- Verlg New York D A Weilnd S nd Ipichino L 218 Model pproimtion of thermo-fluidic diffuion procee in ptilly interconnected tructure In 218 Europen Control Conference ECC IEEE Dtko R 197 Etending theorem of AM Lipunov to Hilbert pce Journl of Mthemticl nlyi nd ppliction Fridmn E nd Orlov Y 29 Eponentil tbility of liner ditributed prmeter ytem with time-vrying dely Automtic Ghlwt A nd Peet MM 217 A conve um-ofqure pproch to nlyi tte feedbck nd output feedbck control of prbolic pde IEEE Trnction on Automtic Control doi: 1119/TAC Krtic M 29 Compenting tring pde in the ctution or ening pth of n untble ode In Americn Control Conference 29 ACC 9 IEEE Krtic M nd Smyhlyev A 28 Bcktepping boundry control for firt-order hyperbolic pde nd ppliction to ytem with ctutor nd enor dely Sytem & Control Letter Meyer E nd Peet MM 215 Stbility nlyi of prbolic liner pde with two ptil dimenion uing lypunov method nd o In Deciion nd Control CDC 215 IEEE 54th Annul Conference on IEEE Nicie S Vlein J nd Fridmn E 29 Stbility of the het nd of the wve eqution with boundry timevrying dely Dicrete nd Continuou Dynmicl Sytem- Serie S 23 Ppchritodoulou A Anderon J Vlmorbid G Prjn S Seiler P nd Prrilo PA 213 SOS- TOOLS: Sum of qure optimition toolbo for MAT- LAB Ppchritodoulou A nd Peet MM 26 On the nlyi of ytem decribed by cle of prtil differentil eqution In Proceeding of the 45th IEEE Conference on Deciion nd Control doi: 1119/CDC Peet MM 218 A new tte-pce repreenttion of lypunov tbility for coupled pde nd clble tbility nlyi in the o frmework riv preprint riv: Qurteroni A nd Veneini A 23 Anlyi of geometricl multicle model bed on the coupling of ode nd pde for blood flow imultion Multicle Modeling & Simultion Rvindrn SS 2 A reduced-order pproch for optiml control of fluid uing proper orthogonl decompoition Interntionl journl for numericl method in fluid Shivkumr S nd Peet MM 218 Computing Input- Output Propertie of Coupled PDE ytem riv e- print riv: Stinner C Surulecu C nd Winkler M 214 Globl wek olution in pde-ode ytem modeling multicle cncer cell invion SIAM Journl on Mthemticl Anlyi Suto GA nd Krtic M 21 Control of pde ode ccde with neumnn interconnection Journl of the Frnklin Intitute Tng S nd ie C 21 Stbilition of coupled pde-ode ytem by boundry control In 49th IEEE Conference on Deciion nd Control CDC doi:1119/cdc Vlmorbid G Ahmdi M nd Ppchritodoulou A 214 Semi-definite progrmming nd functionl in-

8 equlitie for ditributed prmeter ytem In 53rd IEEE Conference on Deciion nd Control doi:1119/cdc Vlmorbid G Ahmdi M nd Ppchritodoulou A 216 Stbility nlyi for cl of prtil differentil eqution vi emidefinite progrmming IEEE Trnction on Automtic Control Appendi A PROOF OF LEMMA 7 Lemm 7 Suppoe R no nd W 22 b nd B c colb b with B c R 2np 4np h full row rnk Then ẋ ABL1 C BL 2 C ż BL 3 C ABL 4 C implie ẋ P ż {UV1V 2YW 1W 2} A1 A2 U Λ 1 Λ 2 Y Π 1 Π 2 A 2 B 1 L 4 C 3 V 1 Ω 1 Ω 2 BL 2 C 3 W 1 Ψ 1 Ψ 2 V 2 Σ 1 Σ 2 W 2 Υ 1 Υ 2 with Λ1 Ω 1 Σ 1 Π 1 Ψ 1 Υ 1 ABL 1 CBL 2 C B 1 L 3 CA B 1 L 4 C IG 1 G 2 Λ2 Ω 2 Σ 2 Π 2 Ψ 2 Υ 2 BL 2 C b A 1 B 1 L 4 C b IG 3 G 4 b I C 3 C 1 B d B c C 1 A3 I Moreover the definition of G 1 G 2 G 3 G 4 re given ccording to Lemm 6 Proof Firt of ll uing 8 we cn how tht C 1 b C 3 d A4 b I C 3 C 1 B d B c C 1 A5 I Now we ue the definition ABC in 7 nd re-write 33 ẋ ABL1 C BL 2 C ż BL 3 C A BL 4 C BL2 C b BL2 C 3 A 1 BL 4 C b A 2 BL 4 C 3 Bed on 7 A : A A 1 : A 1 A 2 : A 2 C : b C d A6 C b : b nd C 3 : C b d C 3 d Rewriting in the form of P {PQ1Q 2SR 1R 2} ẋ ż P {ABL1CBL 2C B 1L 3CA B 1L 4C } P {BL2C b A 1B 1L 4C b } P {BL2C 3A 2B 1L 4C 3} A7 Now ubtituting Lemm 6 in A8 we obtin ẋ P ż {Λ1Ω 1Σ 1Π 1Ψ 1Υ 1} P {Λ2Ω 2Σ 2Π 2Ψ 2Υ 2} P {BL2C 3A 2B 1L 4C 3} A8 P {Λ1Ω 1Σ 1Π 1Ψ 1Υ 1} i defined P {ABL1CBL 2C B 1L 3CA B 1L 4C }P {IG1G 2} Similrly P {Λ2Ω 2Σ 2Π 2Ψ 2Υ 2} i defined P {BL2C b A 1B 1L 4C b }P {IG3G 4} A reult P {Λ1Ω 1Σ 1Π 1Ψ 1Υ 1} P {Λ2Ω 2Σ 2Π 2Ψ 2Υ 2} P {BL2C 3A 2B 1L 4C 3} P {UV1V 2YW 1W 2} Appendi B PROOF OF THEOREM 8 Theorem 8 Suppoe R no nd W 22 b nd B c colb b with B c R 2np 4np h full row rnk Moreover uppoe there eit mtri P R no no nd bounded polynomil function Q : b R no np S : b R np np R 1 R 2 : b b R np np nd ABL1 C BL Ā 2 C B1 BL 3 C ABL 4 C Then Ā P {PQQ SR 1R 2}Ā P {PQQ SR 1R 2} P { ˆKˆLˆL ˆM ˆN 1 ˆN 2} ˆK K K ˆM M M ˆL L1 L 1 ˆL L 2 L 2 ˆN1 N 1 N 1 ˆN2 N 2 N 2 K L1 L 2 M N 1 N 2 KL1 L 2 MN 1 N 2 KL1 L 2 MN 1 N 2 IG1 G 2 B2

9 PQQ SR 1 R 2 U 1 V 1 V 2 YW 1 W 2 B3 Moreover G 1 G 2 re defined in Lemm 6 nd U V 1 V 2 Y W 1 W 2 re defined in Lemm 7 Proof The time derivtive of the Lypunov function i P {PQQ SR 1R 2}Ā Ā P {PQQ SR 1R 2} P {U1V 1V 2YW 1W 2} P {PQQ SR 1R 2} }{{} :Γ 1 P {PQQ SR 1R 2}P {U1V 1V 2YW 1W 2} } {{ } :Γ 2 B4 Now tking two individul term eprtely nd ubtituting the identity 23 we obtin Γ 2 P {IG 1G 2} P {PQQ SR 1R 2}P {U1V 1V 2YW 1W 2} B5 P {KL1L 2MN 1N 2} Similrly Γ 1 P {U1V 1V 2YW 1W 2} P {PQQ SR 1R 2}P {IG1G 2} P {U 1V 1V 2YW 1W 2} B6 P {PQQ SR 1R 2}P {IG1G 2} P {KL 1L 2MN 1N 2} Now uing the djoint opertion of the opertor we obtin P {KL1L 2MN 1N 2} P {KL 1L 2MN 1N 2} B7 P { ˆKˆLˆL ˆM ˆN 1 ˆN 2} with P { ˆKˆLˆL ˆM ˆN 1 ˆN 2} P {KL 1L 2MN 1N 2} P {KL 1L 2MN 1N 2} Thi complete the proof Appendi C COMPOSITION OF THE OPERATORS Lemm 4 For ny mtrice AP R m m nd bounded function B 1 Q 1 : b R m n B 2 Q 2 : b R n m DS : b R n n C i R i : b b R n n with i {12} the following identity hold P {AB1B 2DC 1C 2}P {PQ1Q 2SR 1R 2} P {ˆPˆQ1 ˆQ 2ŜˆR 1ˆR 2} ˆP AP B 1 Q 2 d ˆQ 1 AQ 1 B 1 S B 1 R 2 d ˆQ 2 B 2 P DQ 2 C 2 Q 2 d B 1 R 1 d C 1 Q 2 d Ŝ DS ˆR 1 B 2 Q 1 DR 1 C 1 S C 1 θr 2 θdθ C 2 θr 1 θdθ C 1 θr 1 θdθ ˆR 2 B 2 Q 1 DR 2 C 2 S Proof Suppoe C 1 θr 2 θdθ C 2 θr 1 θdθ P {AB1B 2DC 1C 2} P {ˆPˆQ1 ˆQ 2ŜˆR 1ˆR 2} P {PQ1Q 2SR 1R 2} 1 C 2 θr 2 θdθ Let 1 y1 P {PQ1Q 2SR 1R 2} 2 y 2 y 1 P 1 Q 1 2 d nd y 2 Q 2 1 S 2 R 1 2 d R 2 2 d Then P {AB1B 2DC 1C 2}P {PQ1Q 2SR 1R 2} P {AB1B 2DC 1C 2} y1 1 2 y Ay 1 B 1y 2 d nd 2 B 2 y 1 Dy 2 C 1y 2 d C 2y 2 d

10 Finding the compoition i tright-forwrd lgebric opertion We do thi by epnding ech term eprtely Firtly B 1 y 2 d B 1 Q 2 1 S 2 R 1 2 d R 2 2 d d Seperting the term involving 1 nd 2 we obtin 1 AP B 1 Q 2 } {{ } ˆP ˆQ 1 AQ 1 B 1 S 1 B 1 R 2 d ˆQ 1 2 d B 1 R 1 d Net we epnd the term of 2 nd pecificlly collect the term hving 1 We obtin B 2 P DQ 2 C 1 Q 2 d 1 ˆQ 2 1 Net grouping the term hving 2 we obtin DS 2 Ŝ 2 Similrly we cn group the term involving 2 θ B 2 Q 1 DR 1 C 1 S C 1 θr 1 θdθ C 2 θr 1 θdθ ˆR 1 2 d 2 d C 1 θr 2 θdθ nd b B 2 Q 1 DR 2 C 2 S θ C 2 θr 2 θdθ C 2 θr 1 θdθ ˆR 2 2 d Thi complete the proof 2 d C 1 θr 2 θdθ P{PQ1Q2SR1R2} y R m L n 2 b P {PQ 1Q 2SR 1R 2} P {ˆPˆQ 1 ˆQ 2ŜˆR 1ˆR 2} with D1 ˆP P Ŝ S ˆQ 1 Q 2 ˆR1 R2 ˆQ 2 Q 1 ˆR2 R1 D2 Proof We ue the fct tht for ny clr we hve 1 y1 Let nd y 2 y 2 Then P {PQ1 Q 2 SR 1 R 2 }y 1 Py 1 R m L n 2 b 1 Q 1y 2d 2 Sy 2d y 1 P 1 2 R 2y 2dd 1 Q 2 y 2d y 2 S 2d y 1 ˆP 1 y 2 R 1 2dd 1 ˆQ 1y 2d y 2 Ŝ2d y2 ˆR 2 2dd P { ˆP ˆQ1 ˆQ2ŜˆR1ˆR2} y P {PQ 1 Q 2 SR 1 R 2 }y ˆP P ˆQ 1 Q 2 ˆQ 2 Q 1 Thi complete the proof y 2 Q 2 1d 2 R 1y 2dd y 2Q 1 1d y 2 R 2 2dd R m L n 2 b R m L n 2 b y 2 ˆQ 2 1d y 2 ˆR 1 2dd Ŝ S ˆR1 R 2 ˆR2 R 1 Appendi D ADJOINT OF THE OPERATORS Lemm 5 For ny mtrice P R m m nd bounded function Q 1 : b R m n Q 2 : b R n m S : b R n n R 1 R 2 : b b R n n the following identity hold for ny R m y L n 2 b P {PQ1Q 2SR 1R 2}y R m L n 2 b