Aircraft safety modeling for time-limited dispatch

Size: px
Start display at page:

Download "Aircraft safety modeling for time-limited dispatch"

Transcription

1 Loughborough Universiy Insiuionl Reposiory Aircrf sfey modeling for ime-limied dispch This iem ws submied o Loughborough Universiy's Insiuionl Reposiory by he/n uhor. Ciion: PRESCOTT, D.R. nd ANDREWS, J.D., 005. Aircrf sfey modeling for ime-limied dispch. IN: Proceedings of he Annul Relibiliy nd Mininbiliy Symposium, RAMS'05, Virgini, USA, 4h-7h Jnury, pp [DOI:0.09/RAMS ] Addiionl Informion: This is conference pper [ c IEEE]. I is lso vilble from: hp://ieeexplore.ieee.org/ Personl use of his meril is permied. However, permission o reprin/republish his meril for dverising or promoionl purposes or for creing new collecive works for resle or redisribuion o servers or liss, or o reuse ny copyrighed componen of his work in oher works mus be obined from he IEEE. Med Record: hps://dspce.lboro.c.uk/34/3834 Publisher: c Insiue of Elecricl nd Elecronics Engineers (IEEE) Plese cie he published version.

2 This iem ws submied o Loughborough s Insiuionl Reposiory (hps://dspce.lboro.c.uk/) by he uhor nd is mde vilble under he following Creive Commons Licence condiions. For he full ex of his licence, plese go o: hp://creivecommons.org/licenses/by-nc-nd/.5/

3 Aircrf Sfey Modeling For Time-Limied Dispch Drren R. Presco, Loughborough Universiy John D. Andrews, Loughborough Universiy Key Words: Time-Limied Dispch, FADEC sysem relibiliy, Mone Crlo simulion SUMMARY & CONCLUSIONS This pper offers n lernive mehod of modeling he Time-Limied Dispch (TLD) of ircrf. Exising mehods involve he use of ful ree nlysis nd Mrkov nlysis wih vrious simplifying ssumpions. Mone Crlo simulion (MCS) is he suggesed lernive, which overcomes he problems ssocied wih he oher echniques, such s dependencies beween bsic evens (ful ree nlysis) or huge number of sysem ses (Mrkov nlysis). The resuls obined from he nlysis of simple exmple re compred for he exising modeling pproches nd MCS. MCS is seen o hve poenil dvnges, especilly when modeling TLD for lrge, full scle sysems.. BACKGROUND Inroduced o commercil rnspor ircrf bou 0 yers go, Full Auhoriy Digil Elecronic Conrol (FADEC) sysems govern engine hrus from he ime fuel meering begins o he poin of fuel shuoff. Unil he inroducion of hese elecronic engine conrol sysems hydromechnicl conrol (HMC) sysems were used. In hese pplicions of FADEC i ws o be he firs ime h pilos would hve no HMC sysems vilble s bckup in he even of n elecronic sysem filure []. FADEC sysems conin cerin level of redundncy, incorporing dul chnnel conrol sysem. This involves hving wo essenilly idenicl chnnels per engine. Ech criicl loop or funcion in he FADEC conins eiher redundn elemens or dul sysems. Despie his redundncy, he dispch crieri imposed fer he inroducion of FADEC were overly resricive, incresing he numbers of delys nd cncellions of flighs []. This ws due o he occurrence of independen fuls in more hn one chnnel. Becuse levels of relibiliy re higher for FADEC sysems hn for HMC sysems n opporuniy exised o use vilble redundncy o llow dispch wih fuls presen. This would sill llow irworhiness sndrds o be me nd lso reduce he numbers of delys nd cncellions, long wih he dded benefi of llowing beer plnning of minennce operions. This new pproch, llowing degrded redundncy dispch ws nmed ime-limied dispch (TLD).. INTRODUCTION TO TLD TLD llows he dispch of ircrf in he presence of one or more known fuls whils ssuring cerin level of sysem relibiliy. Depending on he significnce of he fuls ircrf my be dispched for differing lenghs of ime. These dispch inervls give he mximum lengh of ime h he ircrf my be dispched wih fuls presen before hose fuls mus be ddressed. There re four cegories of dispch inervl. These re: o Do No Dispch (DND), o Shor Time Dispch (), o Long Time Dispch (), o Mnufcurer/Operor Defined Dispch (MDD). The implemenion of he DND dispch cegory mens h he ircrf mus no be dispched becuse of he fuls presen nd minennce mus be underken immediely. The dispch cegory llows dispch in he shor-erm nd he cegory for relively longer ime before repirs re crried ou. The finl cegory of dispch, MDD, is reserved for fuls flling ino none of he oher hree cegories nd no ffecing he loss of hrus conrol (LOTC) re []. An upper limi for he LOTC re of 00 evens per 0 6 fligh hours (fl. hrs.) is given by he FAA for dispchble sysem configurions. The mximum verge LOTC re of he sysem mus no exceed 0 evens per 0 6 fl. hrs. This level mches h which ws chieved by he HMC sysems superceded by FADEC.. Minennce Sregies Two sregies exis h my be used o minin he FADEC sysems of ircrf when TLD is pplied. The firs of hese is minimum equipmen lis (MEL) minennce [], which would normlly be used on cegory fuls. The exc ime of occurrence of ny ful minined using his pproch mus be known. As he ful occurs coundown of he dispch ime is iniied. The ful mus be remedied, he les, by he ime he coundown reches zero. Figure illusres his process. The ful occurs ime, fer which he dispch inervl is iniied. A he end of he coundown,, he ful mus be repired, if he repir hs no lredy been implemened. The second minennce sregy used is periodic inspecion/repir (PIR), which involves checking he sysem for fuls regulr inervls. This is mos ofen used o minin fuls in he cegory. Unlike MEL minennce PIR does no require knowledge of he exc ime of occurrence of he ful. Fuls, discovered inspecions, re ssumed o occur he midpoin of consecuive inspecions []. This is considered resonble since he ful will, on verge, occur his ime, ssuming he filure res for fuls re consn wih ime nd he periodic inspecion inervl is less hn he men ime beween filures (MTBF) of RAMS /05/$ IEEE Auhorized licensed use limied o: IEEE Xplore. Downloded on Ocober 3, :30 from IEEE Xplore. Resricions pply.

4 dispch inervl Figure. MEL Minennce. A B Figure 3. Muliple Fuls dispch inervl T I f I A B 3 Figure. PIR Minennce. he sum of filure res in h cegory. The dispch inervl is hen deemed o hve begun he midpoin of he inspecions nd he llowble period of dispch fer he inspecion is clculed. This minennce sregy is illusred in Figure. I nd I represen wo consecuive inspecions. A ful, occurring ime f, is discovered I nd ssumed o occur. If he dispch inervl is ssumed o begin his ime he ircrf my be dispched unil ime, giving n llowble period of dispch of ime T I. In prcice he inspecion inervl for ful cegory mus no exceed wice he dispch inervl for fuls of h cegory. This ensures h he verge exposure o fuls does no exceed he dispch inervl. If PIR is used o minin fuls of more hn one dispch cegory siuion cn occur where ful is discovered n inspecion for fuls of differen cegory. In his cse he ful could be reed s if found for he firs ime he nex inspecion for fuls of is own dispch cegory. In his wy ful discovered inspecion could be reed s if discovered he nex inspecion for fuls []. Due o he differen pproches involved in he MEL nd PIR minennce sregies he mximum possible exposure ime of he sysem o fuls will differ. In MEL he mximum possible exposure ime is equl o he dispch inervl, bu in PIR he mximum possible exposure ime is equl o wice he dispch inervl. There is possibiliy h more hn one ful my exis wihin FADEC sysem ny one ime. In such cses hese fuls my be repired in number of differen wys. Below re some exmples of he siuions h my rise, long wih some of he minennce possibiliies. Figure 3 depics he occurrence of wo fuls, A nd B. If hese fuls were o be minined using MEL minennce dispch inervls for hese fuls would end nd respecively. Upon reching he end of he firs dispch inervl,, number of sregies re possible. A his poin ful A mus be ddressed o llow furher dispch. Ful B my lso be repired his ime, llowing dispch unil noher TLD ful occurs. Alernively, ful B my be lef in he sysem nd he ircrf my be dispched unil is ssocied dispch inervl ends ime. In he siuion jus described i my be h wih he simulneous occurrence of fuls A nd B reducion in he dispch inervl is specified. This scenrio is depiced in Figure 4, gin for he cse of MEL minennce. In his cse fuls A nd B, when occurring lone, bring bou he iniiion of inervls, which would end nd s before. However, when boh fuls re presen wihin he sysem he dispch inervl is reduced o. Thus, s ful B occurs he inervl ending 3 is iniied, rher hn he one ending. Upon reching 3 hree minennce sregies exis. These re: o Repir boh fuls, A nd B, llowing unlimied dispch of he ircrf, o Repir ful A only, llowing dispch unil, o Figure 4. Fuls Acing in Combinion MEL which poin ful B mus be repired, Repir ful B only, llowing dispch unil, which poin ful A mus be repired. When fuls combine in his mnner o reduce he dispch inervl i is possible for he ordering of he fuls o lso ply pr. I my be he cse h he dispch inervl would be reduced when eiher A ws followed by B or B ws followed by A. However, noher possibiliy is h he dispch inervl my only be reduced if A is followed by B bu no if B is followed by A. The siuions occurring bove my lso rise when PIR minennce is being used o minin he FADEC sysem. However, he reducion or oherwise of he dispch inervls cn be more complex, since he exc ime of occurrence of he fuls is no known. As n exmple of his incresed complexiy consider Figure 5, which shows he occurrence of wo fuls A nd B. Ech of hese fuls hs n ssocied inervl nd fuls re o be minined using PIR. In combinion he fuls A nd B iniie inervl. In he figure I nd I represen wo consecuive PIR inspecions. A inspecion I fuls A nd B re discovered nd ssumed o hve occurred he midpoin of he wo inspecions,. A inervl iniied his ime llows dispch of he ircrf for ime T fer inspecion I. However, if he simulneous exisence of fuls A nd B cuses he iniiion of inervl ending 3 he siuion is compliced somewh. Upon reching he minennce dedline 3 fuls A nd B, ful A lone or ful B lone could be repired. If A nd B re boh repired he ircrf my be dispched indefiniely. If jus ful A or ful B is repired he ircrf my hen be dispched unil when he remining ful mus be repired. These exmples of siuions h my rise when RAMS /05/$ IEEE Auhorized licensed use limied o: IEEE Xplore. Downloded on Ocober 3, :30 from IEEE Xplore. Resricions pply.

5 I A B I 3 Figure 5. Fuls Acing in Combinion PIR implemening TLD o FADEC sysems by no mens cover ll possibiliies. They merely highligh some of he mny cses h mus be considered when emping o model TLD. 3. MODELING TLD Before pplying TLD o ircrf i is imporn o be sure h he ircrf will sill conform o he levels of sfey desired of i. The firs sep in demonsring his conformiy is o consruc mhemicl model of he sysem in order o predic is behvior when TLD is used. In [] nd [3] wo of he commonly-used mehods of TLD nlysis re presened ful ree-bsed pproch nd Mrkov modeling pproch. These pproches re demonsred here for he simple sysem shown in Figure 6. Mone Crlo simulion is hen pplied o he sme sysem. The sysem is h given s n exmple in [] nd consiss of wo unis, U nd U. These hve filure res per hour of U = nd U = respecively. 3. Time-Weighed Averge (TWA) Ful Tree Approch In his mehod he overll verge filure re of he sysem is obined by dding he filure res of he HMC fuls, HMC, nd he uncovered fuls, UC, o imeweighed verge (TWA) of he filure res of he sysem from ech of is dispchble configurions, i.e., TWA = HMC + UC + FUFU,L +,L +,L () where TWA is he TWA filure re of he sysem. FU, nd re respecively he frcions of ime spen in he fullup, nd dispchble sysem configurions. FU,L is he LOTC re of he sysem from he FU se o LOTC.,L nd,l re he verge LOTC res wih nd fuls. If he sysem hs n dispchble configurions, le se i = represen he full-up configurion (i.e. he se wih no filed componens). Ses i =, m nd i = m +,,n will represen he nd dispchble sysem configurions respecively. If i,l is he filure re o LOTC for he i h configurion i my be clculed s follows. The filure probbiliy (of LOTC) for ech se is divided by suible ime period, such s he verge fligh ime, o Uni fils U T Sysem Fils Uni fils U Figure 6. A Dul Uni Sysem obin probbiliy per fligh hour. This is hen equed o he verge filure re o LOTC over he ime inervl []. Thus Q i,l i,l =, () where Q i,l is he filure probbiliy (of LOTC) for se i nd fl is he verge fligh ime. Define T nd T s he lenghs of ime spen in he nd dispchble sysem ses before repir, i.e. he nd dispch inervls. If i is he filure re ino filed dispchble sysem se hen he frcions of ime spen in he nd dispchble sysem configurions my be pproximed s follows (see[3]): = m T nd = T. (3) i n i m+ The frcion of ime spen in he full-up se, FU, is deermined using he fc h he frcion of ime spen in he filed ses dded o he frcion of ime in full-up se will be uniy. Therefore once he frcions of ime spen in ll he filed ses is known, FU my be clculed s follows: FU =. (4) The verge LOTC res wih nd fuls re defined in [3] s: m i n ii,l ii,l m+, L = nd = m,l n (5) i m+ The dul uni sysem shown in Figure 6 hs hree dispchble configurions, numbered ses o 3 s follows:. U nd U work (Full-up),. U is filed, U works (), 3. U works, U is filed (). If q U nd q U represen he probbiliies of filure of U nd U hen, subsiuing hese ino equion (), we ge: ququ qu qu,l =,,L =, 3,L =. (6) Subsiuing hese ino equion (3) we obin he frcions of ime spen in he nd dispchble sysem configurions: qu qu,, = T = T (7) which my hen be subsiued ino (4) o give he frcion of ime spen in he full-up se: qu qu. FU = T T (8) Subsiuing (6), (7) nd (8) ino () gives TWA for his dul uni sysem. If he filure probbiliies q Ui re given by: Ui qui = e, (9) hen TWA cn be clculed for differen vlues of T nd T, which re he dispch inervls for U nd U respecively. The fligh ime, fl, ws ssumed o be 5 hours. 3. Single Ful Se Mrkov Approch Like he ful ree pproch, his echnique hs severl RAMS /05/$ IEEE Auhorized licensed use limied o: IEEE Xplore. Downloded on Ocober 3, :30 from IEEE Xplore. Resricions pply.

6 rcive feures []. For insnce, Mrkov modeling mkes direc use of filure nd repir res, neging he need o eque he filure res o probbiliy/fligh hour, s is he cse wih he ful ree pproch. The min drwbck wih Mrkov nlysis is he poenilly huge number of sysem ses h my be presen wihin he model [4]. For his reson [3] considers filures h led o one finl filed se represening LOTC. This filed se is bsorbing, since once he sysem is in his se no rnsiion my occur o ny oher se. However, in [3] n rificil feedbck loop, simuled repir, is inroduced h kes he sysem from he finl filed se bck o he full-up se wih no componens filed. I is rgued h sedy-se soluion for he verge filure re of he sysem is required bu h wihou feedbck loop he sedy-se probbiliy of he sysem being filed will lwys pproch vlue of. In [3] his simuled repir re is se o uniy, for simpliciy. This choice is rbirry, mde since, if LOTC even occurs for n engine, he conrol mus be repired before he nex fligh kes plce. A repir inervl of one hour is deemed suible. This mehod lso differs from convenionl Mrkov echnique in h i involves only lower order ful ses, commonly only single-ful ses, lhough dul ful ses re included if imporn for he nlysis. The Mrkov model obined using his echnique is like h in Figure 7. The rnsiion res from he full-up se o he nd ful ses re nd nd he corresponding repir res re ν nd ν. In [3] nd re pproximed by he sums of he filure res of ll of he nd ype fuls respecively. ν nd ν re given by he reciprocl of he dispch inervls, T nd T. The filure res from he nd ful ses re given by,l nd,l nd he simuled repir is given by ν FB. The Mrkov model leds o sysem of 4 liner differenil equions, given by: Q = QA, (0) where A is he rnsiion re mrix, ( + ) 0 ( ν + ) ν,l 0,L A = () ν ( + ) 0 ν,l,l ν FB 0 0 ν FB nd Q = [ QFU, Q, Q, QLOTC ], () where Q FU, Q, Q nd Q LOTC re he probbiliies of he sysem being in he full-up,, nd LOTC ses respecively. A sedy se he derivives of hese probbiliies will be zero, i.e. Q A = 0, (3) nd hus sysem of 4 liner equions is obined. This sysem of equions is dependen. In order o obin sysem of independen equions one of he equions is replced by he condiion h he sum of he probbiliies of being in ech se is uniy, i.e. Q FU + Q + Q + QLOTC =. (4) Column of A is rbirrily chosen o be replced by (4), giving sysem of equions h is solved o find Q sedy se. In order o find he LOTC re for he sysem he Full-up ν ν,l LOTC,L verge rnsiion re ino he LOTC se is considered [3]. This LOTC re, obined from he single ful se Mrkov model is: Probbiliy flow ino LOTC se Mkv =, (5) Probbiliy of being in LOTC se A his poin fuls h re required o be modeled s leding direcly o he LOTC se, i.e. HMC fuls or uncovered elecronic fuls, re dded. Since hese fuls could led o LOTC from ny non-lotc se he sum of heir filure res is muliplied by he sum of he probbiliies of being in ll he ses excep he LOTC se [3]. Thus (5) becomes: ( QFU + Q + Q)( HMC + UC ) + Q,L + Q,L Mkv =.(6) QLOTC When modeling he dul uni sysem shown in Figure 6 he filure res for he single ful se Mrkov model my be clculed s hey were for he TWA model. Mrkov models of TLD ll inherenly model he MEL minennce sregy. This is becuse s soon s componen of he sysem fils he sysem will mke rnsiion o differen se nd hus he filure ime of he componen will be known. If PIR is o be used s minennce pproch he mximum llowed dispch ime clculed using Mrkov pproch is doubled in order o clcule he mximum periodic inspecion inervl. 3.3 Convenionl Mrkov Approch o Single Ful Se Model ν FB Figure 7. Single Ful Se Mrkov Digrm. Noe h he LOTC re obined from he single ful se Mrkov model, Mkv, is dependen on he vlue of he feedbck re, ν FB, nd independen of he iniil condiions of he sysem. Andrews nd Moss [4] deil mehod of finding he sympoic filure re of sysem wih bsorbing finl filed ses. This my be pplied o he Mrkov model described in he previous secion. In his cse here is no need o include simuled repir from he finl fully-filed se, i.e. he feedbck re, ν FB. Thus he rnsiion re mrix is: ( + ) 0 ( ) ν ν +,L 0,L A = (7) ν 0 ( + ) ν,l,l The reduced rnsiion re mrix, A m, is obined by runcing A by deleing enries for he bsorbing ses. Thus he finl row nd column of (7) is deleed o ge A m. The sympoic filure re, ASY, is hen given by: RAMS /05/$ IEEE Auhorized licensed use limied o: IEEE Xplore. Downloded on Ocober 3, :30 from IEEE Xplore. Resricions pply.

7 A m ASY = (8) 0 Q FU () 0 Q () 0 Q () 0,, 3, where he elemen of he reduced rnsiion re mrix in he i h row nd j h column is defined s i,j nd Q FU (0), Q (0) nd Q (0) re he iniil probbiliies of being in he FU, nd sysem ses. Therefore, noe h in his cse he soluion is dependen on he iniil condiions of he sysem.,, 3, 3.4 Mone Crlo Simulion Mehod Performing Mone Crlo simulion (MCS) involves consrucing nd using compuer model of he sysem under invesigion. This model is bsed on srucured, logicl sysem represenion nd conins se of rules h describe he response of he sysem o evens h occur o i during is lifeime [4]. These evens my be componen repirs, filures or sequences of filures or, in his cse, evens h will led o he implemenion of TLD nd is subsequen minennce dedlines. In fc, ny even my be dded o he model subjec o he complexiy of he sofwre used. MCS requires he generion of uniform se of rndom numbers. These rndom numbers re hen used o genere imes of componen filures by using he relevn filure disribuion ssocied wih ech componen. In he simulion code consruced during he course of his work he numbers genered re, in fc, pseudo-rndom. The sysem simulion my hen be run mny imes, mking noe fer ech simulion of ny prmeers of ineres, nd coninuing he simulions unil he required olernce is reched for hese prmeers. In he work presened here he prmeer of ineres h hd o be obined from he sysem ws is filure re. Ech simulion is crried ou unil suible poin in ime, wheher his be he lifeime of he sysem or he ime which he sysem fils. The lgorihm for he min module of he simulion progrm used in his work is given in Figure 8. When modeling sysem on which TLD will be pplied he correc scheduling of ll occurring evens is of gre impornce. In he MCS code used here he componen filure imes re iniilly genered nd dded o n rry h holds reference o he componen, long wih he ime h is filure will occur. A reference is lso kep of he order of ll he evens in he schedule nd his is upded s evens re removed from, or dded o, he schedule. As componen filures occur in he simulion lis of TLD crieri is checked o see if he filure of his componen will iniie, eiher on is own or in combinion wih oher componen filures, he implemenion of TLD dedline. If his is he cse he dedline iself is hen dded o he schedule long wih is ime of occurrence, which vries ccording o he dispch cegory for h priculr filure or combinion of filures. In ech simulion loop he firs even chronologiclly is removed from he schedule nd he simulion ime is dvnced o ime of his even. If his ime exceeds he,3,3 3,3 red sysem logic, componen filure nd repir disribuions nd TLD combinions from file inpu TLD opions se mx_sim_ime ol_sim_ime = 0, ol_fil = 0, filure_re = 0 while (filure re hs converged) mrker for simulion filure iniilized: flg = 0 iniil componen filures dded o empy schedule while (curren_sim_ ime mx_sim_ ime) curren_sim_ime = ime of s even in schedule remove s even from schedule if ( s even is componen filure) upde he sus of he sysem if (sysem fils due o his componen filure) simulion filed: flg = end his simulion if (TLD inervl iniied due o his filure) dd minennce dedlines o schedule else if ( s even is minennce dedline) implemen repirs ccording o minennce sregy else if ( s even is periodic inspecion) dd minennce dedline o schedule ol_sim_ime = ol_sim_ime + curren_sim_ime ol_fil = ol_fil + flg filure_re = ol_fil / ol_sim_ime check for convergence of filure re Figure 8. The Algorihm for he Min Module of he MCS. mximum lifeime of he sysem he simulion ends immediely. Oherwise, if he even is componen filure, he sus of he sysem is upded nd if he sysem fils he simulion ends. If he even is minennce dedline he pproprie repirs re crried ou o componens ccording o he minennce sregies being pplied o he sysem. In he cse of he sysem being modeled here his simply mens repiring he uni h filed o cuse he minennce dedline. Finlly, if he even is periodic inspecion, minennce dedline will be dded o he schedule he pproprie poin in ime. In he MCS code number of oher ssumpions re mde. These re: o The mximum lifeime of he sysem is se o be RAMS /05/$ IEEE Auhorized licensed use limied o: IEEE Xplore. Downloded on Ocober 3, :30 from IEEE Xplore. Resricions pply.

8 00000 fligh hours. This corresponds o lifeime of bou 37 yers for sysem used for 5 hours of fligh per dy, or 55 yers for sysem used for 0 hours per dy. o The lengh of fligh is 5 hours. This ws chosen since his ws used in he exmple in [] when fligh lengh is needed for he TWA modeling pproch. o Minennce operions nd inspecions cnno occur mid-fligh. Thus, for ny evens of his ype whose occurrence ime is clculed o occur in fligh, he ime of he even is moved o he beginning of h fligh. This is conservive pproch o his problem. I should be noed h, despie he fc h hey hve been used in order o model his priculr sysem, hese ssumpions my be esily chnged wihin he MCS code, if necessry. The code my lso be esily pplied o oher sysems. This is becuse he sysem upon which he simulion is o be implemened is supplied o he code in represenion of he form of is ful ree. Oher sysemspecific informion supplied o he code is he filure disribuion of ech componen, long wih is ssocied prmeers. The filure re of he sysem ws cully clculed fer every 000 simulions. Resuls were obined o n ccurcy of deciml plces (for filure re given in unis of number of filures per 0 6 fl. hrs.). If he filure re mched over 0 consecuive clculions convergence ws ssumed. 4. RESULTS Resuls were obined using ech of he bove modeling pproches for differing inervls of dispch for ech uni in he dul uni exmple. The dispch inervl for U () is represened s D nd he dispch inervl for U () s D. The MCS code ws used o model ll possible minennce sregies for he dul uni sysem, wih U nd U being minined using MEL nd/or PIR. The PIR minennce sregy ws simuled wih differing inspecion inervls, incremenlly from 0.5 o imes he lengh of he dispch inervl. Resuls re shown in Figure 9 for vlue of D=00 fl. hrs. s ws chosen in []. Presened in his wy i cn be seen h he resuls obined from ll of he mehods re similr. The resuls obined using he TWA pproch give n overesimion of he filure re in comprison o he oher pproches, nd ech of he MCS models gives similr resuls o he Mrkov model, which is cknowledged o be he more ccure of he wo pproches presened in [] nd [3]. Figure 0 shows he sme resuls presened s percenge difference from he resuls of he single ful se Mrkov pproch. The TWA resuls grow from 7.% bove hose of he single ful se Mrkov o 6.4% s D increses from 00 o 500 fl. hrs. All of he simulion resuls lie wihin 4% of he single ful se Mrkov resuls for D = 00 fl. hrs. As D increses from 00 o 500 fl. hrs. he simulion resuls remin close o he single ful se Mrkov resuls. However, s D rises he simulion resuls become progressively lrger in comprison o he single ful se Mrkov resuls. Indeed, when D = 500 fl. hrs. he resuls from ll of he simulions re greer hn he resuls single no. of filures per 0 6 fl. hrs. % difference D (fl. hrs.) D (fl. hrs.) fuls se model, by s much s 5.95%. A rend h is lso observed from he simulion resuls is h, for he mos pr, he sysem filure re is observed o be higher for he cse where boh unis re minined using PIR hn when one is repired using PIR nd one using MEL. The sysem filure res when boh unis re minined using MEL is he lowes obined from ll of he minennce sregies nd end o lie closes o he single ful se Mrkov resuls. These resuls re s would be expeced for his simple sysem. I should be noed h, for his simple dul uni sysem, ll of he modeling pproches br TWA give similr llowble dispch ime for he remining uni, given dispch ime for one of he unis, ssuming h filure re of 0 filures per 0 6 fl. hrs. mus no be exceeded. Noe h he resuls obined for he single ful se Mrkov model wih n rificil feedbck loop re, for his sysem, very close o hose chieved by using he convenionl Mrkov pproch o modeling he single ful se model wih iniil condiion h he sysem srs in he full-up se, i.e. boh unis work iniilly. 5. CONCLUSIONS MCS, U:PIR, U:PIR MCS, U:MEL, U:PIR MCS, U:PIR, U:MEL MCS, U:MEL, U:MEL TWA Single Ful Se Mrkov Figure 9. Comprison of Mehods D=00 fl. hrs. MCS, U:PIR, U:PIR MCS, U:MEL, U:PIR MCS, U:PIR, U:MEL MCS, U:MEL, U:MEL TWA Figure 0. Percenge Difference from Single Ful Se Mrkov D=00 fl. hrs. Ech of he mehods pplied o he dul uni sysem brings wih i boh dvnges nd disdvnges, some of he min exmples of which re shown in Tble. The MCS mehod proves o be mch for exising echniques for he simple dul uni sysem considered. However, he sligh differences in sysem filure re observed beween he MCS nd he single ful se Mrkov could grow lrger s he RAMS /05/$ IEEE Auhorized licensed use limied o: IEEE Xplore. Downloded on Ocober 3, :30 from IEEE Xplore. Resricions pply.

9 Mehod Advnges Disdvnges TWA ful ree Cler udible represenion of sysem filure logic. TLD inroduces dependencies beween he bsic evens, which re inpproprie for convenionl ful ree pproch. single ful se Mrkov/ convenionl Mrkov MCS Overcome problems ssocied wih dependencies in he sysem. Filure nd repir res re esily incorpored ino he model. Dependencies, non-consn filure res nd (ordered) combinions of fuls re esily del wih. PIR nd MEL minennce re esily modeled, including differen sregies for differen dispch cegories, e.g. -MEL, -PIR. Uses ful ree represenion of he sysem, which brings he ssocied dvnge of clriy of sysem filure logic represenion. More informion my be obined from he model hn wih oher mehods e.g. could idenify fuls wih inpproprie dispch inervls. Filure res pproximed from filure probbiliies. Only he MEL minennce sregy is modeled. Consn filure res re required. Model grows exponenilly nd becomes incresingly difficul o udi s he number of sysem ses increses (convenionl model). Involves n pproximion of he sysem, reducing he number of sysem ses nd ignoring mny combinions of fuls (single ful se only). Iniil generion of code my be ime-consuming however, he code cn be used for mny sysems. Exr CPU ime required o run mny simulions however, in comprison o he ime ken designing new ircrf nd gining cerificion his will be negligible. Tble. Advnges nd Disdvnges of he Differen Modeling Techniques. sysem being modeled becomes more complex. Preliminry work o invesige his heory, in which slighly more complex sysem is modeled, suggess he difference beween he echniques is greer hn for he dul uni sysem. REFERENCES. H. Lrsen, G. Horn, Time-Limied Dispch: An Inercive Trining nd Self-Sudy Course., Keybridge Technologies, Inc., 00.. FAA Memorndum: Policy for Time-Limied Dispch (TLD) of Engines Fied wih Full Auhoriy Digil Engine Conrols (FADEC) Sysems, June 9 00, Policy No. ANE TLD-R. 3. Guidelines for Time-Limied-Dispch (TLD) Anlysis for Elecronic Engine Conrol Sysems, SAE ARP 507, SAE Inernionl, J. D. Andrews, T. R. Moss, Relibiliy nd Risk Assessmen, nd Ediion, London, Professionl Engineering Publishing, 00. BIOGRAPHIES Drren R. Presco Deprmen of Aeronuicl nd Auomoive Engineering, Loughborough Universiy, Loughborough, Leicesershire, LE 3TU, UK e-mil: D.R.Presco@lboro.c.uk Drren Presco is currenly sudying for docore Loughborough Universiy, hving previously grdued wih firs clss honors degree in Mhemics nd gining msers degree wih disincion in Indusril Mhemicl Modeling. His curren reserch sudies re underken s pr of he Risk nd Relibiliy group in he Aeronuicl nd Auomoive Engineering deprmen of Loughborough Universiy. John D. Andrews Deprmen of Aeronuicl nd Auomoive Engineering, Loughborough Universiy, Loughborough, Leicesershire, LE 3TU, UK e-mil: J.D.Andrews@lboro.c.uk John Andrews is Professor of Sysems Relibiliy in he Deprmen of Aeronuicl nd Auomoive Engineering. He joined Loughborough Universiy in 989 hving previously gined nine yers indusril reserch experience wih Briish Gs. His curren reserch ineress concern he ssessmen of he sfey nd risk of poenilly hzrdous indusril civiies. This reserch hs been hevily suppored by indusril funding. Over recen yers grns hve been secured from BAE Sysems, MOD, Rolls-Royce, ExxonMobil nd Bechel. Professor Andrews hs numerous journl/conference publicions long wih joinly uhored book Relibiliy nd Risk Assessmen which is now in is second ediion. ACKNOWLEDGEMENT The uhors would like o cknowledge Phil Wilkinson for his involvemen nd suppor during he course of his work. RAMS /05/$ IEEE Auhorized licensed use limied o: IEEE Xplore. Downloded on Ocober 3, :30 from IEEE Xplore. Resricions pply.

e t dt e t dt = lim e t dt T (1 e T ) = 1

e t dt e t dt = lim e t dt T (1 e T ) = 1 Improper Inegrls There re wo ypes of improper inegrls - hose wih infinie limis of inegrion, nd hose wih inegrnds h pproch some poin wihin he limis of inegrion. Firs we will consider inegrls wih infinie

More information

Contraction Mapping Principle Approach to Differential Equations

Contraction Mapping Principle Approach to Differential Equations epl Journl of Science echnology 0 (009) 49-53 Conrcion pping Principle pproch o Differenil Equions Bishnu P. Dhungn Deprmen of hemics, hendr Rn Cmpus ribhuvn Universiy, Khmu epl bsrc Using n eension of

More information

The solution is often represented as a vector: 2xI + 4X2 + 2X3 + 4X4 + 2X5 = 4 2xI + 4X2 + 3X3 + 3X4 + 3X5 = 4. 3xI + 6X2 + 6X3 + 3X4 + 6X5 = 6.

The solution is often represented as a vector: 2xI + 4X2 + 2X3 + 4X4 + 2X5 = 4 2xI + 4X2 + 3X3 + 3X4 + 3X5 = 4. 3xI + 6X2 + 6X3 + 3X4 + 6X5 = 6. [~ o o :- o o ill] i 1. Mrices, Vecors, nd Guss-Jordn Eliminion 1 x y = = - z= The soluion is ofen represened s vecor: n his exmple, he process of eliminion works very smoohly. We cn elimine ll enries

More information

4.8 Improper Integrals

4.8 Improper Integrals 4.8 Improper Inegrls Well you ve mde i hrough ll he inegrion echniques. Congrs! Unforunely for us, we sill need o cover one more inegrl. They re clled Improper Inegrls. A his poin, we ve only del wih inegrls

More information

3. Renewal Limit Theorems

3. Renewal Limit Theorems Virul Lborories > 14. Renewl Processes > 1 2 3 3. Renewl Limi Theorems In he inroducion o renewl processes, we noed h he rrivl ime process nd he couning process re inverses, in sens The rrivl ime process

More information

Minimum Squared Error

Minimum Squared Error Minimum Squred Error LDF: Minimum Squred-Error Procedures Ide: conver o esier nd eer undersood prolem Percepron y i > for ll smples y i solve sysem of liner inequliies MSE procedure y i = i for ll smples

More information

Minimum Squared Error

Minimum Squared Error Minimum Squred Error LDF: Minimum Squred-Error Procedures Ide: conver o esier nd eer undersood prolem Percepron y i > 0 for ll smples y i solve sysem of liner inequliies MSE procedure y i i for ll smples

More information

September 20 Homework Solutions

September 20 Homework Solutions College of Engineering nd Compuer Science Mechnicl Engineering Deprmen Mechnicl Engineering A Seminr in Engineering Anlysis Fll 7 Number 66 Insrucor: Lrry Creo Sepember Homework Soluions Find he specrum

More information

1.0 Electrical Systems

1.0 Electrical Systems . Elecricl Sysems The ypes of dynmicl sysems we will e sudying cn e modeled in erms of lgeric equions, differenil equions, or inegrl equions. We will egin y looking fmilir mhemicl models of idel resisors,

More information

Motion. Part 2: Constant Acceleration. Acceleration. October Lab Physics. Ms. Levine 1. Acceleration. Acceleration. Units for Acceleration.

Motion. Part 2: Constant Acceleration. Acceleration. October Lab Physics. Ms. Levine 1. Acceleration. Acceleration. Units for Acceleration. Moion Accelerion Pr : Consn Accelerion Accelerion Accelerion Accelerion is he re of chnge of velociy. = v - vo = Δv Δ ccelerion = = v - vo chnge of velociy elpsed ime Accelerion is vecor, lhough in one-dimensionl

More information

1 jordan.mcd Eigenvalue-eigenvector approach to solving first order ODEs. -- Jordan normal (canonical) form. Instructor: Nam Sun Wang

1 jordan.mcd Eigenvalue-eigenvector approach to solving first order ODEs. -- Jordan normal (canonical) form. Instructor: Nam Sun Wang jordnmcd Eigenvlue-eigenvecor pproch o solving firs order ODEs -- ordn norml (cnonicl) form Insrucor: Nm Sun Wng Consider he following se of coupled firs order ODEs d d x x 5 x x d d x d d x x x 5 x x

More information

S Radio transmission and network access Exercise 1-2

S Radio transmission and network access Exercise 1-2 S-7.330 Rdio rnsmission nd nework ccess Exercise 1 - P1 In four-symbol digil sysem wih eqully probble symbols he pulses in he figure re used in rnsmission over AWGN-chnnel. s () s () s () s () 1 3 4 )

More information

A Kalman filtering simulation

A Kalman filtering simulation A Klmn filering simulion The performnce of Klmn filering hs been esed on he bsis of wo differen dynmicl models, ssuming eiher moion wih consn elociy or wih consn ccelerion. The former is epeced o beer

More information

f t f a f x dx By Lin McMullin f x dx= f b f a. 2

f t f a f x dx By Lin McMullin f x dx= f b f a. 2 Accumulion: Thoughs On () By Lin McMullin f f f d = + The gols of he AP* Clculus progrm include he semen, Sudens should undersnd he definie inegrl s he ne ccumulion of chnge. 1 The Topicl Ouline includes

More information

A Time Truncated Improved Group Sampling Plans for Rayleigh and Log - Logistic Distributions

A Time Truncated Improved Group Sampling Plans for Rayleigh and Log - Logistic Distributions ISSNOnline : 39-8753 ISSN Prin : 347-67 An ISO 397: 7 Cerified Orgnizion Vol. 5, Issue 5, My 6 A Time Trunced Improved Group Smpling Plns for Ryleigh nd og - ogisic Disribuions P.Kvipriy, A.R. Sudmni Rmswmy

More information

0 for t < 0 1 for t > 0

0 for t < 0 1 for t > 0 8.0 Sep nd del funcions Auhor: Jeremy Orloff The uni Sep Funcion We define he uni sep funcion by u() = 0 for < 0 for > 0 I is clled he uni sep funcion becuse i kes uni sep = 0. I is someimes clled he Heviside

More information

PHYSICS 1210 Exam 1 University of Wyoming 14 February points

PHYSICS 1210 Exam 1 University of Wyoming 14 February points PHYSICS 1210 Em 1 Uniersiy of Wyoming 14 Februry 2013 150 poins This es is open-noe nd closed-book. Clculors re permied bu compuers re no. No collborion, consulion, or communicion wih oher people (oher

More information

Neural assembly binding in linguistic representation

Neural assembly binding in linguistic representation Neurl ssembly binding in linguisic represenion Frnk vn der Velde & Mrc de Kmps Cogniive Psychology Uni, Universiy of Leiden, Wssenrseweg 52, 2333 AK Leiden, The Neherlnds, vdvelde@fsw.leidenuniv.nl Absrc.

More information

Chapter 2: Evaluative Feedback

Chapter 2: Evaluative Feedback Chper 2: Evluive Feedbck Evluing cions vs. insrucing by giving correc cions Pure evluive feedbck depends olly on he cion ken. Pure insrucive feedbck depends no ll on he cion ken. Supervised lerning is

More information

EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR A SECOND-ORDER ITERATIVE BOUNDARY-VALUE PROBLEM

EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR A SECOND-ORDER ITERATIVE BOUNDARY-VALUE PROBLEM Elecronic Journl of Differenil Equions, Vol. 208 (208), No. 50, pp. 6. ISSN: 072-669. URL: hp://ejde.mh.xse.edu or hp://ejde.mh.un.edu EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR A SECOND-ORDER ITERATIVE

More information

Physics 2A HW #3 Solutions

Physics 2A HW #3 Solutions Chper 3 Focus on Conceps: 3, 4, 6, 9 Problems: 9, 9, 3, 41, 66, 7, 75, 77 Phsics A HW #3 Soluions Focus On Conceps 3-3 (c) The ccelerion due o grvi is he sme for boh blls, despie he fc h he hve differen

More information

An integral having either an infinite limit of integration or an unbounded integrand is called improper. Here are two examples.

An integral having either an infinite limit of integration or an unbounded integrand is called improper. Here are two examples. Improper Inegrls To his poin we hve only considered inegrls f(x) wih he is of inegrion nd b finie nd he inegrnd f(x) bounded (nd in fc coninuous excep possibly for finiely mny jump disconinuiies) An inegrl

More information

Temperature Rise of the Earth

Temperature Rise of the Earth Avilble online www.sciencedirec.com ScienceDirec Procedi - Socil nd Behviorl Scien ce s 88 ( 2013 ) 220 224 Socil nd Behviorl Sciences Symposium, 4 h Inernionl Science, Socil Science, Engineering nd Energy

More information

ENGR 1990 Engineering Mathematics The Integral of a Function as a Function

ENGR 1990 Engineering Mathematics The Integral of a Function as a Function ENGR 1990 Engineering Mhemics The Inegrl of Funcion s Funcion Previously, we lerned how o esime he inegrl of funcion f( ) over some inervl y dding he res of finie se of rpezoids h represen he re under

More information

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 9: The High Beta Tokamak

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 9: The High Beta Tokamak .65, MHD Theory of Fusion Sysems Prof. Freidberg Lecure 9: The High e Tokmk Summry of he Properies of n Ohmic Tokmk. Advnges:. good euilibrium (smll shif) b. good sbiliy ( ) c. good confinemen ( τ nr )

More information

A LIMIT-POINT CRITERION FOR A SECOND-ORDER LINEAR DIFFERENTIAL OPERATOR IAN KNOWLES

A LIMIT-POINT CRITERION FOR A SECOND-ORDER LINEAR DIFFERENTIAL OPERATOR IAN KNOWLES A LIMIT-POINT CRITERION FOR A SECOND-ORDER LINEAR DIFFERENTIAL OPERATOR j IAN KNOWLES 1. Inroducion Consider he forml differenil operor T defined by el, (1) where he funcion q{) is rel-vlued nd loclly

More information

5.1-The Initial-Value Problems For Ordinary Differential Equations

5.1-The Initial-Value Problems For Ordinary Differential Equations 5.-The Iniil-Vlue Problems For Ordinry Differenil Equions Consider solving iniil-vlue problems for ordinry differenil equions: (*) y f, y, b, y. If we know he generl soluion y of he ordinry differenil

More information

A 1.3 m 2.5 m 2.8 m. x = m m = 8400 m. y = 4900 m 3200 m = 1700 m

A 1.3 m 2.5 m 2.8 m. x = m m = 8400 m. y = 4900 m 3200 m = 1700 m PHYS : Soluions o Chper 3 Home Work. SSM REASONING The displcemen is ecor drwn from he iniil posiion o he finl posiion. The mgniude of he displcemen is he shores disnce beween he posiions. Noe h i is onl

More information

( ) ( ) ( ) ( ) ( ) ( y )

( ) ( ) ( ) ( ) ( ) ( y ) 8. Lengh of Plne Curve The mos fmous heorem in ll of mhemics is he Pyhgoren Theorem. I s formulion s he disnce formul is used o find he lenghs of line segmens in he coordine plne. In his secion you ll

More information

A new model for limit order book dynamics

A new model for limit order book dynamics Anewmodelforlimiorderbookdynmics JeffreyR.Russell UniversiyofChicgo,GrdueSchoolofBusiness TejinKim UniversiyofChicgo,DeprmenofSisics Absrc:Thispperproposesnewmodelforlimiorderbookdynmics.Thelimiorderbookconsiss

More information

Convergence of Singular Integral Operators in Weighted Lebesgue Spaces

Convergence of Singular Integral Operators in Weighted Lebesgue Spaces EUROPEAN JOURNAL OF PURE AND APPLIED MATHEMATICS Vol. 10, No. 2, 2017, 335-347 ISSN 1307-5543 www.ejpm.com Published by New York Business Globl Convergence of Singulr Inegrl Operors in Weighed Lebesgue

More information

MATH 124 AND 125 FINAL EXAM REVIEW PACKET (Revised spring 2008)

MATH 124 AND 125 FINAL EXAM REVIEW PACKET (Revised spring 2008) MATH 14 AND 15 FINAL EXAM REVIEW PACKET (Revised spring 8) The following quesions cn be used s review for Mh 14/ 15 These quesions re no cul smples of quesions h will pper on he finl em, bu hey will provide

More information

Procedia Computer Science

Procedia Computer Science Procedi Compuer Science 00 (0) 000 000 Procedi Compuer Science www.elsevier.com/loce/procedi The Third Informion Sysems Inernionl Conference The Exisence of Polynomil Soluion of he Nonliner Dynmicl Sysems

More information

Solutions to Problems from Chapter 2

Solutions to Problems from Chapter 2 Soluions o Problems rom Chper Problem. The signls u() :5sgn(), u () :5sgn(), nd u h () :5sgn() re ploed respecively in Figures.,b,c. Noe h u h () :5sgn() :5; 8 including, bu u () :5sgn() is undeined..5

More information

MTH 146 Class 11 Notes

MTH 146 Class 11 Notes 8.- Are of Surfce of Revoluion MTH 6 Clss Noes Suppose we wish o revolve curve C round n is nd find he surfce re of he resuling solid. Suppose f( ) is nonnegive funcion wih coninuous firs derivive on he

More information

MAT 266 Calculus for Engineers II Notes on Chapter 6 Professor: John Quigg Semester: spring 2017

MAT 266 Calculus for Engineers II Notes on Chapter 6 Professor: John Quigg Semester: spring 2017 MAT 66 Clculus for Engineers II Noes on Chper 6 Professor: John Quigg Semeser: spring 7 Secion 6.: Inegrion by prs The Produc Rule is d d f()g() = f()g () + f ()g() Tking indefinie inegrls gives [f()g

More information

INTEGRALS. Exercise 1. Let f : [a, b] R be bounded, and let P and Q be partitions of [a, b]. Prove that if P Q then U(P ) U(Q) and L(P ) L(Q).

INTEGRALS. Exercise 1. Let f : [a, b] R be bounded, and let P and Q be partitions of [a, b]. Prove that if P Q then U(P ) U(Q) and L(P ) L(Q). INTEGRALS JOHN QUIGG Eercise. Le f : [, b] R be bounded, nd le P nd Q be priions of [, b]. Prove h if P Q hen U(P ) U(Q) nd L(P ) L(Q). Soluion: Le P = {,..., n }. Since Q is obined from P by dding finiely

More information

Magnetostatics Bar Magnet. Magnetostatics Oersted s Experiment

Magnetostatics Bar Magnet. Magnetostatics Oersted s Experiment Mgneosics Br Mgne As fr bck s 4500 yers go, he Chinese discovered h cerin ypes of iron ore could rc ech oher nd cerin mels. Iron filings "mp" of br mgne s field Crefully suspended slivers of his mel were

More information

REAL ANALYSIS I HOMEWORK 3. Chapter 1

REAL ANALYSIS I HOMEWORK 3. Chapter 1 REAL ANALYSIS I HOMEWORK 3 CİHAN BAHRAN The quesions re from Sein nd Shkrchi s e. Chper 1 18. Prove he following sserion: Every mesurble funcion is he limi.e. of sequence of coninuous funcions. We firs

More information

GENERALIZATION OF SOME INEQUALITIES VIA RIEMANN-LIOUVILLE FRACTIONAL CALCULUS

GENERALIZATION OF SOME INEQUALITIES VIA RIEMANN-LIOUVILLE FRACTIONAL CALCULUS - TAMKANG JOURNAL OF MATHEMATICS Volume 5, Number, 7-5, June doi:5556/jkjm555 Avilble online hp://journlsmhkueduw/ - - - GENERALIZATION OF SOME INEQUALITIES VIA RIEMANN-LIOUVILLE FRACTIONAL CALCULUS MARCELA

More information

On the Pseudo-Spectral Method of Solving Linear Ordinary Differential Equations

On the Pseudo-Spectral Method of Solving Linear Ordinary Differential Equations Journl of Mhemics nd Sisics 5 ():136-14, 9 ISS 1549-3644 9 Science Publicions On he Pseudo-Specrl Mehod of Solving Liner Ordinry Differenil Equions B.S. Ogundre Deprmen of Pure nd Applied Mhemics, Universiy

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

More information

Solutions for Nonlinear Partial Differential Equations By Tan-Cot Method

Solutions for Nonlinear Partial Differential Equations By Tan-Cot Method IOSR Journl of Mhemics (IOSR-JM) e-issn: 78-578. Volume 5, Issue 3 (Jn. - Feb. 13), PP 6-11 Soluions for Nonliner Pril Differenil Equions By Tn-Co Mehod Mhmood Jwd Abdul Rsool Abu Al-Sheer Al -Rfidin Universiy

More information

LAPLACE TRANSFORM OVERCOMING PRINCIPLE DRAWBACKS IN APPLICATION OF THE VARIATIONAL ITERATION METHOD TO FRACTIONAL HEAT EQUATIONS

LAPLACE TRANSFORM OVERCOMING PRINCIPLE DRAWBACKS IN APPLICATION OF THE VARIATIONAL ITERATION METHOD TO FRACTIONAL HEAT EQUATIONS Wu, G.-.: Lplce Trnsform Overcoming Principle Drwbcks in Applicion... THERMAL SIENE: Yer 22, Vol. 6, No. 4, pp. 257-26 257 Open forum LAPLAE TRANSFORM OVEROMING PRINIPLE DRAWBAKS IN APPLIATION OF THE VARIATIONAL

More information

Chapter Direct Method of Interpolation

Chapter Direct Method of Interpolation Chper 5. Direc Mehod of Inerpolion Afer reding his chper, you should be ble o:. pply he direc mehod of inerpolion,. sole problems using he direc mehod of inerpolion, nd. use he direc mehod inerpolns o

More information

A Simple Method to Solve Quartic Equations. Key words: Polynomials, Quartics, Equations of the Fourth Degree INTRODUCTION

A Simple Method to Solve Quartic Equations. Key words: Polynomials, Quartics, Equations of the Fourth Degree INTRODUCTION Ausrlin Journl of Bsic nd Applied Sciences, 6(6): -6, 0 ISSN 99-878 A Simple Mehod o Solve Quric Equions Amir Fhi, Poo Mobdersn, Rhim Fhi Deprmen of Elecricl Engineering, Urmi brnch, Islmic Ad Universi,

More information

Average & instantaneous velocity and acceleration Motion with constant acceleration

Average & instantaneous velocity and acceleration Motion with constant acceleration Physics 7: Lecure Reminders Discussion nd Lb secions sr meeing ne week Fill ou Pink dd/drop form if you need o swich o differen secion h is FULL. Do i TODAY. Homework Ch. : 5, 7,, 3,, nd 6 Ch.: 6,, 3 Submission

More information

Tax Audit and Vertical Externalities

Tax Audit and Vertical Externalities T Audi nd Vericl Eernliies Hidey Ko Misuyoshi Yngihr Ngoy Keizi Universiy Ngoy Universiy 1. Inroducion The vericl fiscl eernliies rise when he differen levels of governmens, such s he federl nd se governmens,

More information

A new model for solving fuzzy linear fractional programming problem with ranking function

A new model for solving fuzzy linear fractional programming problem with ranking function J. ppl. Res. Ind. Eng. Vol. 4 No. 07 89 96 Journl of pplied Reserch on Indusril Engineering www.journl-prie.com new model for solving fuzzy liner frcionl progrmming prolem wih rning funcion Spn Kumr Ds

More information

(b) 10 yr. (b) 13 m. 1.6 m s, m s m s (c) 13.1 s. 32. (a) 20.0 s (b) No, the minimum distance to stop = 1.00 km. 1.

(b) 10 yr. (b) 13 m. 1.6 m s, m s m s (c) 13.1 s. 32. (a) 20.0 s (b) No, the minimum distance to stop = 1.00 km. 1. Answers o Een Numbered Problems Chper. () 7 m s, 6 m s (b) 8 5 yr 4.. m ih 6. () 5. m s (b).5 m s (c).5 m s (d) 3.33 m s (e) 8. ().3 min (b) 64 mi..3 h. ().3 s (b) 3 m 4..8 mi wes of he flgpole 6. (b)

More information

On Source and Channel Codes for Multiple Inputs and Outputs: Does Multiple Description Meet Space Time? 1

On Source and Channel Codes for Multiple Inputs and Outputs: Does Multiple Description Meet Space Time? 1 On Source nd Chnnel Codes for Muliple Inpus nd Oupus: oes Muliple escripion Mee Spce Time? Michelle Effros Rlf Koeer 3 Andre J. Goldsmih 4 Muriel Médrd 5 ep. of Elecricl Eng., 36-93, Cliforni Insiue of

More information

Probability, Estimators, and Stationarity

Probability, Estimators, and Stationarity Chper Probbiliy, Esimors, nd Sionriy Consider signl genered by dynmicl process, R, R. Considering s funcion of ime, we re opering in he ime domin. A fundmenl wy o chrcerize he dynmics using he ime domin

More information

Systems Variables and Structural Controllability: An Inverted Pendulum Case

Systems Variables and Structural Controllability: An Inverted Pendulum Case Reserch Journl of Applied Sciences, Engineering nd echnology 6(: 46-4, 3 ISSN: 4-7459; e-issn: 4-7467 Mxwell Scienific Orgniion, 3 Submied: Jnury 5, 3 Acceped: Mrch 7, 3 Published: November, 3 Sysems Vribles

More information

Nonlinear System Modelling: How to Estimate the. Highest Significant Order

Nonlinear System Modelling: How to Estimate the. Highest Significant Order IEEE Insrumenion nd Mesuremen Technology Conference nchorge,, US, - My Nonliner Sysem Modelling: ow o Esime he ighes Significn Order Neophyos Chirs, Ceri Evns nd Dvid Rees, Michel Solomou School of Elecronics,

More information

Green s Functions and Comparison Theorems for Differential Equations on Measure Chains

Green s Functions and Comparison Theorems for Differential Equations on Measure Chains Green s Funcions nd Comprison Theorems for Differenil Equions on Mesure Chins Lynn Erbe nd Alln Peerson Deprmen of Mhemics nd Sisics, Universiy of Nebrsk-Lincoln Lincoln,NE 68588-0323 lerbe@@mh.unl.edu

More information

Mathematics 805 Final Examination Answers

Mathematics 805 Final Examination Answers . 5 poins Se he Weiersrss M-es. Mhemics 85 Finl Eminion Answers Answer: Suppose h A R, nd f n : A R. Suppose furher h f n M n for ll A, nd h Mn converges. Then f n converges uniformly on A.. 5 poins Se

More information

1. Consider a PSA initially at rest in the beginning of the left-hand end of a long ISS corridor. Assume xo = 0 on the left end of the ISS corridor.

1. Consider a PSA initially at rest in the beginning of the left-hand end of a long ISS corridor. Assume xo = 0 on the left end of the ISS corridor. In Eercise 1, use sndrd recngulr Cresin coordine sysem. Le ime be represened long he horizonl is. Assume ll ccelerions nd decelerions re consn. 1. Consider PSA iniilly res in he beginning of he lef-hnd

More information

3 Motion with constant acceleration: Linear and projectile motion

3 Motion with constant acceleration: Linear and projectile motion 3 Moion wih consn ccelerion: Liner nd projecile moion cons, In he precedin Lecure we he considered moion wih consn ccelerion lon he is: Noe h,, cn be posiie nd neie h leds o rie of behiors. Clerl similr

More information

Question Details Int Vocab 1 [ ] Question Details Int Vocab 2 [ ]

Question Details Int Vocab 1 [ ] Question Details Int Vocab 2 [ ] /3/5 Assignmen Previewer 3 Bsic: Definie Inegrls (67795) Due: Wed Apr 5 5 9: AM MDT Quesion 3 5 6 7 8 9 3 5 6 7 8 9 3 5 6 Insrucions Red ody's Noes nd Lerning Gols. Quesion Deils In Vocb [37897] The chnge

More information

Some basic notation and terminology. Deterministic Finite Automata. COMP218: Decision, Computation and Language Note 1

Some basic notation and terminology. Deterministic Finite Automata. COMP218: Decision, Computation and Language Note 1 COMP28: Decision, Compuion nd Lnguge Noe These noes re inended minly s supplemen o he lecures nd exooks; hey will e useful for reminders ou noion nd erminology. Some sic noion nd erminology An lphe is

More information

SOME USEFUL MATHEMATICS

SOME USEFUL MATHEMATICS SOME USEFU MAHEMAICS SOME USEFU MAHEMAICS I is esy o mesure n preic he behvior of n elecricl circui h conins only c volges n currens. However, mos useful elecricl signls h crry informion vry wih ime. Since

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

ON NEW INEQUALITIES OF SIMPSON S TYPE FOR FUNCTIONS WHOSE SECOND DERIVATIVES ABSOLUTE VALUES ARE CONVEX

ON NEW INEQUALITIES OF SIMPSON S TYPE FOR FUNCTIONS WHOSE SECOND DERIVATIVES ABSOLUTE VALUES ARE CONVEX Journl of Applied Mhemics, Sisics nd Informics JAMSI), 9 ), No. ON NEW INEQUALITIES OF SIMPSON S TYPE FOR FUNCTIONS WHOSE SECOND DERIVATIVES ABSOLUTE VALUES ARE CONVEX MEHMET ZEKI SARIKAYA, ERHAN. SET

More information

1. Introduction. 1 b b

1. Introduction. 1 b b Journl of Mhemicl Inequliies Volume, Number 3 (007), 45 436 SOME IMPROVEMENTS OF GRÜSS TYPE INEQUALITY N. ELEZOVIĆ, LJ. MARANGUNIĆ AND J. PEČARIĆ (communiced b A. Čižmešij) Absrc. In his pper some inequliies

More information

A comparison of modelling approaches for the time-limited dispatch (TLD) of aircraft

A comparison of modelling approaches for the time-limited dispatch (TLD) of aircraft Loughborough University Institutional Repository A comparison of modelling approaches for the time-limited dispatch (TLD) of aircraft This item was submitted to Loughborough University's Institutional

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

Introduction to LoggerPro

Introduction to LoggerPro Inroducion o LoggerPro Sr/Sop collecion Define zero Se d collecion prmeers Auoscle D Browser Open file Sensor seup window To sr d collecion, click he green Collec buon on he ool br. There is dely of second

More information

Phys 110. Answers to even numbered problems on Midterm Map

Phys 110. Answers to even numbered problems on Midterm Map Phys Answers o een numbered problems on Miderm Mp. REASONING The word per indices rio, so.35 mm per dy mens.35 mm/d, which is o be epressed s re in f/cenury. These unis differ from he gien unis in boh

More information

AJAE appendix for Is Exchange Rate Pass-Through in Pork Meat Export Prices Constrained by the Supply of Live Hogs?

AJAE appendix for Is Exchange Rate Pass-Through in Pork Meat Export Prices Constrained by the Supply of Live Hogs? AJAE ppendix for Is Exchnge Re Pss-Through in Por Me Expor Prices Consrined by he Supply of Live Hogs? Jen-Philippe Gervis Cnd Reserch Chir in Agri-indusries nd Inernionl Trde Cener for Reserch in he Economics

More information

Computer-Aided Analysis of Electronic Circuits Course Notes 3

Computer-Aided Analysis of Electronic Circuits Course Notes 3 Gheorghe Asachi Technical Universiy of Iasi Faculy of Elecronics, Telecommunicaions and Informaion Technologies Compuer-Aided Analysis of Elecronic Circuis Course Noes 3 Bachelor: Telecommunicaion Technologies

More information

Evaporation losses as a major factor in determining allowable yield from water supply reservoirs: the case of Botswana s major reservoirs

Evaporation losses as a major factor in determining allowable yield from water supply reservoirs: the case of Botswana s major reservoirs River Bsin Mngemen III 631 Evporion losses s mjor fcor in deermining llowble yield from wer supply reservoirs: he cse of Boswn s mjor reservoirs P. K. Kenbho & B. P. Prid Deprmen of Environmenl Science,

More information

Version 001 test-1 swinney (57010) 1. is constant at m/s.

Version 001 test-1 swinney (57010) 1. is constant at m/s. Version 001 es-1 swinne (57010) 1 This prin-ou should hve 20 quesions. Muliple-choice quesions m coninue on he nex column or pge find ll choices before nswering. CubeUniVec1x76 001 10.0 poins Acubeis1.4fee

More information

A Structural Approach to the Enforcement of Language and Disjunctive Constraints

A Structural Approach to the Enforcement of Language and Disjunctive Constraints A Srucurl Aroch o he Enforcemen of Lnguge nd Disjuncive Consrins Mrin V. Iordche School of Engineering nd Eng. Tech. LeTourneu Universiy Longview, TX 7607-700 Pnos J. Ansklis Dermen of Elecricl Engineering

More information

ECE Microwave Engineering

ECE Microwave Engineering EE 537-635 Microwve Engineering Adped from noes y Prof. Jeffery T. Willims Fll 8 Prof. Dvid R. Jcson Dep. of EE Noes Wveguiding Srucures Pr 7: Trnsverse Equivlen Newor (N) Wveguide Trnsmission Line Model

More information

A NUMERICAL SOLUTION OF THE URYSOHN-TYPE FREDHOLM INTEGRAL EQUATIONS

A NUMERICAL SOLUTION OF THE URYSOHN-TYPE FREDHOLM INTEGRAL EQUATIONS A NUMERICAL SOLUTION OF THE URYSOHN-TYPE FREDHOLM INTEGRAL EQUATIONS A JAFARIAN 1,, SA MEASOOMY 1,b, ALIREZA K GOLMANKHANEH 1,c, D BALEANU 2,3,4 1 Deprmen of Mhemics, Urmi Brnch, Islmic Azd Universiy,

More information

Properties of Logarithms. Solving Exponential and Logarithmic Equations. Properties of Logarithms. Properties of Logarithms. ( x)

Properties of Logarithms. Solving Exponential and Logarithmic Equations. Properties of Logarithms. Properties of Logarithms. ( x) Properies of Logrihms Solving Eponenil nd Logrihmic Equions Properies of Logrihms Produc Rule ( ) log mn = log m + log n ( ) log = log + log Properies of Logrihms Quoien Rule log m = logm logn n log7 =

More information

The order of reaction is defined as the number of atoms or molecules whose concentration change during the chemical reaction.

The order of reaction is defined as the number of atoms or molecules whose concentration change during the chemical reaction. www.hechemisryguru.com Re Lw Expression Order of Recion The order of recion is defined s he number of oms or molecules whose concenrion chnge during he chemicl recion. Or The ol number of molecules or

More information

Unified Framework for Developing Testing Effort Dependent Software Reliability Growth Models

Unified Framework for Developing Testing Effort Dependent Software Reliability Growth Models P. K. Kpur, Omr Shnwi, Anu G. Aggrwl, Rvi Kumr Unified Frmework for Developing Tesing Effor Dependen Sofwre Relibiliy Growh Models P.K. KAPUR 1, OMAR SHATNAI, ANU G. AGGARAL 1, RAVI KUMAR 1 1 Deprmen of

More information

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Swiss Federl Insiue of Pge 1 The Finie Elemen Mehod for he Anlysis of Non-Liner nd Dynmic Sysems Prof. Dr. Michel Hvbro Fber Dr. Nebojs Mojsilovic Swiss Federl Insiue of ETH Zurich, Swizerlnd Mehod of

More information

graph of unit step function t

graph of unit step function t .5 Piecewie coninuou forcing funcion...e.g. urning he forcing on nd off. The following Lplce rnform meril i ueful in yem where we urn forcing funcion on nd off, nd when we hve righ hnd ide "forcing funcion"

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

Estimating the population parameter, r, q and K based on surplus production model. Wang, Chien-Hsiung

Estimating the population parameter, r, q and K based on surplus production model. Wang, Chien-Hsiung SCTB15 Working Pper ALB 7 Esiming he populion prmeer, r, q nd K bsed on surplus producion model Wng, Chien-Hsiung Biologicl nd Fishery Division Insiue of Ocenogrphy Nionl Tiwn Universiy Tipei, Tiwn Tile:

More information

USING ITERATIVE LINEAR REGRESSION MODEL TO TIME SERIES MODELS

USING ITERATIVE LINEAR REGRESSION MODEL TO TIME SERIES MODELS Elecronic Journl of Applied Sisicl Anlysis EJASA (202), Elecron. J. App. S. Anl., Vol. 5, Issue 2, 37 50 e-issn 2070-5948, DOI 0.285/i20705948v5n2p37 202 Universià del Sleno hp://sib-ese.unile.i/index.php/ejs/index

More information

More on Magnetically C Coupled Coils and Ideal Transformers

More on Magnetically C Coupled Coils and Ideal Transformers Appenix ore on gneiclly C Couple Coils Iel Trnsformers C. Equivlen Circuis for gneiclly Couple Coils A imes, i is convenien o moel mgneiclly couple coils wih n equivlen circui h oes no involve mgneic coupling.

More information

For the reaction, R P, the is given by,

For the reaction, R P, the is given by, Dr JADU SAMUEL CHEMICAL KINETICS Inroducion Chemicl ineics is brnch of physicl chemisry, which dels wih he sudy of he re of chemicl recions nd he vrious fcors ffecing i Such sudies lso enble us o elucide

More information

M r. d 2. R t a M. Structural Mechanics Section. Exam CT5141 Theory of Elasticity Friday 31 October 2003, 9:00 12:00 hours. Problem 1 (3 points)

M r. d 2. R t a M. Structural Mechanics Section. Exam CT5141 Theory of Elasticity Friday 31 October 2003, 9:00 12:00 hours. Problem 1 (3 points) Delf Universiy of Technology Fculy of Civil Engineering nd Geosciences Srucurl echnics Secion Wrie your nme nd sudy numer he op righ-hnd of your work. Exm CT5 Theory of Elsiciy Fridy Ocoer 00, 9:00 :00

More information

Robust estimation based on the first- and third-moment restrictions of the power transformation model

Robust estimation based on the first- and third-moment restrictions of the power transformation model h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,

More information

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004 ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform

More information

3D Transformations. Computer Graphics COMP 770 (236) Spring Instructor: Brandon Lloyd 1/26/07 1

3D Transformations. Computer Graphics COMP 770 (236) Spring Instructor: Brandon Lloyd 1/26/07 1 D Trnsformions Compuer Grphics COMP 770 (6) Spring 007 Insrucor: Brndon Lloyd /6/07 Geomery Geomeric eniies, such s poins in spce, exis wihou numers. Coordines re nming scheme. The sme poin cn e descried

More information

Random Walk with Anti-Correlated Steps

Random Walk with Anti-Correlated Steps Random Walk wih Ani-Correlaed Seps John Noga Dirk Wagner 2 Absrac We conjecure he expeced value of random walks wih ani-correlaed seps o be exacly. We suppor his conjecure wih 2 plausibiliy argumens and

More information

RESPONSE UNDER A GENERAL PERIODIC FORCE. When the external force F(t) is periodic with periodτ = 2π

RESPONSE UNDER A GENERAL PERIODIC FORCE. When the external force F(t) is periodic with periodτ = 2π RESPONSE UNDER A GENERAL PERIODIC FORCE When he exernl force F() is periodic wih periodτ / ω,i cn be expnded in Fourier series F( ) o α ω α b ω () where τ F( ) ω d, τ,,,... () nd b τ F( ) ω d, τ,,... (3)

More information

T-Match: Matching Techniques For Driving Yagi-Uda Antennas: T-Match. 2a s. Z in. (Sections 9.5 & 9.7 of Balanis)

T-Match: Matching Techniques For Driving Yagi-Uda Antennas: T-Match. 2a s. Z in. (Sections 9.5 & 9.7 of Balanis) 3/0/018 _mch.doc Pge 1 of 6 T-Mch: Mching Techniques For Driving Ygi-Ud Anenns: T-Mch (Secions 9.5 & 9.7 of Blnis) l s l / l / in The T-Mch is shun-mching echnique h cn be used o feed he driven elemen

More information

15. Vector Valued Functions

15. Vector Valued Functions 1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,

More information

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 10: The High Beta Tokamak Con d and the High Flux Conserving Tokamak.

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 10: The High Beta Tokamak Con d and the High Flux Conserving Tokamak. .615, MHD Theory of Fusion Sysems Prof. Freidberg Lecure 1: The High Be Tokmk Con d nd he High Flux Conserving Tokmk Proeries of he High Tokmk 1. Evlue he MHD sfey fcor: ψ B * ( ) 1 3 ρ 1+ ν ρ ρ cosθ *

More information

14. The fundamental theorem of the calculus

14. The fundamental theorem of the calculus 4. The funmenl heorem of he clculus V 20 00 80 60 40 20 0 0 0.2 0.4 0.6 0.8 v 400 200 0 0 0.2 0.5 0.8 200 400 Figure : () Venriculr volume for subjecs wih cpciies C = 24 ml, C = 20 ml, C = 2 ml n (b) he

More information

Hermite-Hadamard-Fejér type inequalities for convex functions via fractional integrals

Hermite-Hadamard-Fejér type inequalities for convex functions via fractional integrals Sud. Univ. Beş-Bolyi Mh. 6(5, No. 3, 355 366 Hermie-Hdmrd-Fejér ype inequliies for convex funcions vi frcionl inegrls İmd İşcn Asrc. In his pper, firsly we hve eslished Hermie Hdmrd-Fejér inequliy for

More information

Inventory Control of Perishable Items in a Two-Echelon Supply Chain

Inventory Control of Perishable Items in a Two-Echelon Supply Chain Journal of Indusrial Engineering, Universiy of ehran, Special Issue,, PP. 69-77 69 Invenory Conrol of Perishable Iems in a wo-echelon Supply Chain Fariborz Jolai *, Elmira Gheisariha and Farnaz Nojavan

More information

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous

More information

Inventory Management Models with Variable Holding Cost and Salvage value

Inventory Management Models with Variable Holding Cost and Salvage value OSR Journl of Business nd Mngemen OSR-JBM e-ssn: -X p-ssn: 9-. Volume ssue Jul. - Aug. PP - www.iosrjournls.org nvenory Mngemen Models wi Vrile Holding os nd Slvge vlue R.Mon R.Venkeswrlu Memics Dep ollege

More information

ECE Microwave Engineering. Fall Prof. David R. Jackson Dept. of ECE. Notes 10. Waveguides Part 7: Transverse Equivalent Network (TEN)

ECE Microwave Engineering. Fall Prof. David R. Jackson Dept. of ECE. Notes 10. Waveguides Part 7: Transverse Equivalent Network (TEN) EE 537-635 Microwve Engineering Fll 7 Prof. Dvid R. Jcson Dep. of EE Noes Wveguides Pr 7: Trnsverse Equivlen Newor (N) Wveguide Trnsmission Line Model Our gol is o come up wih rnsmission line model for

More information