Tecnologia e Inovação, Lisboa, Portugal. ABB Corporate Research Center, Wallstadter Str. 59, Ladenburg, Germany,

Size: px
Start display at page:

Download "Tecnologia e Inovação, Lisboa, Portugal. ABB Corporate Research Center, Wallstadter Str. 59, Ladenburg, Germany,"

Transcription

1 A New Connuous-Tme Schedulng Fomulaon fo Connuous Plans unde Vaable Eleccy Cos Pedo M. Caso * Io Hajunkosk and Ignaco E. Gossmann Depaameno de Modelação e Smulação de Pocessos Insuo Naconal de Engenhaa Tecnologa e Inovação Lsboa Pougal ABB Copoae eseach Cene Wallsade S Ladenbug Gemany Depamen of Chemcal Engneeng Canege Mellon Unvesy Psbugh Pennsylvana Absac Ths wok addesses he schedulng of connuous plans subjec o enegy consans elaed o me-dependen eleccy pcng and avalably. Dscee and connuous-me fomulaons ae pesened ha can addess hese ssues ogehe wh mulple nemedae due daes. Boh fomulaons ely on he esouce-ask newok pocess epesenaon. The compuaonally pefomance s compaed fo he objecve of oal eleccy mnmzaon wh he esuls favong he dscee-me model due o he moe naual way of handlng such a wde vaey of dscee evens. In pacula can successfully handle poblems of ndusal sze. Neveheless he new connuous-me model s a majo beakhough snce s he fs model of s ype ha s able o effecvely ncopoae me-vaable uly pofles. When compaed o a smple manual schedulng pocedue he poposed schedulng appoaches can lead o poenal eleccy savngs aound 20% by swchng poducon fom peods of hgh o low eleccy cos. * To whom coespondence should be addessed. Tel.: Fax: E- mal: pedo.caso@ne.p. 1

2 1. Movaon Enepses ae cuenly unde pessue o poduce a he lowes possble cos whn connuously changng economc consans 1. To acheve hs goal hey mus acvely look a he bes opeang pacces and opmze hese boh globally and locally. Whn hs oveall goal schedulng plays an mpoan pa. Ths wok s movaed by a eal ndusal poblem ha we canno dsclose fo confdenaly easons. I nvolves he fnal sage of a mulpoduc plan wh eleccy nensve paallel equpmen uns whee schedulng nvolves decdng when each un has o poduce a cean poduc. Mos of he mes hs s made manually by he opeao accodng o heusc ules. The poducs ae hen sen o soage uns whee hey ae soed unl dspachng akes place. Meeng cusome demands on me s val and fo hs eason n some plans no ohe facos ae aken no accoun besdes yng o keep he soage uns full n ode o be able o fulfll he odes. Plan schedulng s dffcul due o he followng facos: lage combnaoal sze asng fom he numbe of equpmen uns poducs and soage uns; vaous opeang and conacual consans; lbealzed eleccy make wh nonanspaen bllng pacces. Due o he nheen complexy he opeao schedulng choces may be fa fom he opmal ones. The mos challengng aspec of plan schedulng s undoubedly he ncopoaon of enegy consans elaed o eleccy pcng and avalably. We consde n hs pape he case whee n he plannng sage conacs ae ageed beween he eleccy supple and he plan whch ofen specfes maxmum levels of powe usage. If eleccy consumpon exceeds hs heshold he plan ncus n sff penales wheeas undepoducon coss he same as planned poducon. Even fo nomal poducon eleccy cos vaes sgnfcanly houghou he day and hs mus be aken no consdeaon whn he modelng famewok. Anohe mpoan aspec concens meeng he sales foecass whch epesen he mnmum amoun of poducs ha has o be avalable n he soage uns dung he me hozon. These ae ypcally consdeed o occu a he end of each day. Accounng fo evens ha occu a pedeemned pons of me fo nsance a a specfc hou of a gven day beng ehe a change n eleccy cos level o he occuence of a demand pon may 2

3 be elavely easy o exemely complex dependng on he ype of me epesenaon employed. Ths s saghfowad wh a dscee-me appoach wheneve a suffcenly fne me gd can be used. In such cases he peseleced me of some of he gd me pons wll mach exacly hose of he pons of change and so he consans ae ease o enfoce. In conas wh connuous-me he absolue me of all even pons s deemned by he solve and hus s much hade o elae he evens wh he pons of change. Ths pape pesens a new connuous-me fomulaon ha effecvely handles me dependen cos paamees and dscee demand pons. Incopoaon of he fome aspec whn a connuousme fomulaon has no been epoed befoe o he bes of ou knowledge wheeas he consans used o model he lae aspec ae concepually smla o hose used by Maavelas and Gossmann 2. The poposed appoach bulds on he geneal mulpupose fomulaon of Caso e al. 3 whch can addess poblems nvolvng bach and connuous asks effcenly. Indeed fo bach plans a sudy 4 has found o be he bes sngle me gd fomulaon. Neveheless he consans ha ae gven ae only suable fo connuous plans meely because bach asks unlke connuous asks canno be dvded n as many mes as equed. In ohe wods whle he oupu maeal fom a bach ask s poduced enely a s end afe a specfed peod of me he one fom a connuous ask s connuously beng poduced. Thus he execuon of a sngle nsance of a connuous ask s equvalen o he execuon (a he same pocessng ae) of mulple nsances poducng he same oal amoun of maeal n sequence dung he same me peod. The new fomulaon s also moe dealed wh espec o he esouce balances. In ode o goously accoun fo nvenoy consans balances need o be consdeed boh a he sa and end of me slos so ha he effec of esouce consumpon and geneaon occung due o connuous asks s consdeed n full as explaned n Schllng & Paneldes 5. The new fomulaon s bul on a unfed famewok fo pocess epesenaon he esouce-task Newok (TN) of Paneldes 6. Ths means ha he model vaables and consans ae wen s ems of absac enes lke esouces asks and even pons so has a much wde scope ha he sngle sage ndusal case sudy used fo llusave puposes. 3

4 2. Poblem Defnon Dung he fnal sage of he pocess an nemedae maeal s ansfomed no one of dffeen fnal poducs chaacezed by chemcal composon and pacle sze dsbuon hough he use of eleccy. These ae hen sen o soage uns whee hey wa unl cusome dspach akes place. Ths pocess s llusaed n Fgue 1. I emphaszes ha a pacula soage un may be suable fo jus a subse of he poducs. In fac he poduc allocaed o a soage un nomally neve changes. Snce s saghfowad o ncopoae such consans n he upcomng models s pefeable o consde full connecvy beween uns o ncease he flexbly. In ohe wods we wll be assumng ha evey soage un can handle all poducs (shaed soage) bu only one a a me. Machne 1 Poduc 1 Soage 1 Cusome 1 Inemedae Machne 2 Poduc 2 Soage 2 Cusome 2 Machne M Poduc P Soage S Cusome 3 Fnal pocessng sage Fgue 1. Fnal pocessng sage of ndusal case sudy Typcal plan schedules ae esablshed ove one week so hs wll be he me hozon (H=168 h) assumed fo he emande of he pape. Naually s saghfowad o consde ohe values. Le M epesen he se of machnes P he se of poducs and S he soage uns. Machnes ae chaacezed by: () powe equemens fo each poduc pw pm [MW]; () pocessng aes ρ pm [on/h]. Le DY and H be ses whose elemens ae he days of he week and he hous of he day especvely. Each poduc may have mulple demands ove he week and can occu a any hou of he day d pdyh [on]. Fnally soage uns have known maxmum capaces cap s [on]. 4

5 The enegy conac sgned by he plan and eleccy povde esablshes a cean pcng polcy. Eleccy cos s ypcally lowe dung he ngh and hghe dung he day. Smlaly o he poduc demands s geneally assumed ha he cos can change evey hou of evey day 19. The maxmum level of oal powe consumpon s n un specfed hough he paamee pwx hdy [MW]. Fo llusaon puposes he enegy polcy gven n Duae e al. s used. I consss of hee enegy levels E wh pces c e of and [ /kwh]. The weekly dsbuon s gven n Fgue 2. MO F SA SU 0.25 Eleccy cos ( /kwh) Hou Fgue 2. Eleccy cos polcy whn a wokng week. The objecve wll be o mnmze he oal enegy cos subjec o consans on esouce avalably ha ncludes pocessng and soage uns and ules (eleccy). 3. esouce-task Newok epesenaon of he Pocess The models poposed n secons 5 and 6 ae bul on he esouce-task Newok o make hem as geneal as possble. The nex sep s hus o denfy he se of pocess esouces and asks whch wll be denfed as a membe of one o moe subypes n ode o bee undesand he sucual paamees geneaon pocess and he dffeen ems n he model consans. Fo ha pupose we use he gudelnes gven n Caso e al. 8 o defne new subypes such ha he specfcs of he ndusal poblem a hand can be ncopoaed. The fnal esul s he TN gven n Fgue 3. Noe ha s helpful o geneae he complee dawng even hough s no absoluely necessay. 5

6 EL S 1 M 1 Soe_P 1 _S 1 Pocess_P 1 _M 1 ae P 1M1... pw P 1M1 pw P 1MM Pocess_P 1 _M M ae P 1MM P 1 _M Send_o_soage ae=gndng ae Hold_n_soage Duaon=1 me neval Soe_P 1 _S S Send_o_soage ae=gndng ae Hold_n_soage Duaon=1 me neval P 1 _S 1 P 1 _S S emove_p 1 _S 1 Insananeous emove_p 1 _S S Insananeous P 1 M M M... S S... Soe_P P _S S Pocess_P P _M M ae P PMM Send_o_soage ae=gndng ae... pw P PMM pw P PM1 P P _M Hold_n_soage Duaon=1 me neval P P _S S emove_p P _S S Insananeous Soe_P P _S 1 Pocess_P P _M 1 ae P PM1 Send_o_soage ae=gndng ae P P Hold_n_soage Duaon=1 me neval P P _S 1 emove_p P _S 1 Insananeous Equpmen esouces nvolved n mng consans ( TC ) Equpmen esouces no nvolved n mng consans esouces poduced o consumed connuously ( CT ) Fnal poduc esoces ( FP ) Uly esoces ( UT ) Connuous asks (I c ) Soage asks (I s ) Tansfe asks (I ) Connuous neacon Dscee neacon Fgue 3. esouce-task Newok epesenaon of fnal pocessng sage esouces Idenfcaon In an TN esouces ae epesened as ccles. Gong fom lef o gh n Fgue 3 we have he maeal npu o he schedulng poblem ha wll be called aw-maeal and named M (he sole membe of se M ) despe he fac ha s n ealy an nemedae pocess maeal. I s connuously consumed whle a poduc s beng poduced. Whle aleady n he sae of fnal poduc sll has o go hough soage and wa unl dspachng. I s hus equed o denfy poduc locaon whch can be n hee dffeen places: () mmedaely afe he machnes LM ; () conaned n soage LS ; o () nsde he fnal anspoaon vessel FP. In he fs suaon hee s no need o dffeenae he machne whee he poduc s poduced so a oal of P esouces ae 6

7 geneaed fom P 1 _M o P P _M whee he M sands fo afe he machnes. In conas n he second possble place we have o dsngush n whch soage un (o goup of uns wokng ogehe) he poduc s locaed snce hey can have dffeen capaces. The esul s he geneaon of addonal P S esouces. esouces connuously poduced and/o consumed ae he elemens of se CT = M LM LS. These ae coloed n dak gey o faclae denfcaon. The fnal locaon gves se o he fnal poduc esouces P 1 o P P whch ae flled n black. And hs ends he maeal esouces. Movng on o he equpmen esouces se EQ we have pocessng equpmen soage equpmen and evenually anspoaon equpmen. They ae howeve handled dffeenly so subse TC s needed. Ths ncludes solely he fome ype whch ae geneally he mos mpoan esouces n a schedulng poblem. They wll have a song mpac on he fnal fom of he connuous-me gd due o he pacpaon n he mng consans. In hs specfc case he machnes ae he membes of TC and ae flled n shadowed whe unlke soage ubs whch have a lgh gey backgound. Fnally hee ae he ules UT he subjec of hs pape. Connuous-me fomulaons 9-12 have so fa consdeed a maxmum uly avalably ha emans consan houghou he me hozon. Fo seam and coolng wae we ae dealng wh flowae values ypcally gven n [on/h] whle fo eleccy we deal wh powe [MW]. Howeve as dscussed eale (secon 2) we wll geneally assume a uly avalably pofle whee n he lm he specfed maxmum avalably can change evey hou of he day dung he ene week Tasks Idenfcaon Tasks ae epesened as ecangles n an TN and mus be chaacezed fully n ems of s esouces unlke n a Sae-Task Newok (STN) epesenaon Moe specfcally he STN epesenaon of hs pocess would show jus P pocessng asks snce hee s no need o explcly say n whch un he poduc s gong o be execued. In conas n he TN of Fgue 3 he pocessng asks mus also be dsaggegaed no all he possble machnes.

8 In Fgue 3 hee ae hee ypes of asks. The pocessng asks flled n whe ae pefomed connuously a a known ae se I c. The execuon of a pocessng ask nvolves fou esouces he ones lnked by aows o he asks. In ode o dffeenae beween connuous and dscee neacon sold and dashed lnes ae especvely used. Some lnes wll have aows on boh decons denong consumpon of he esouce a he sa of he ask and poducon a s end. As an example Pocess_P 1 _M 1 connuously consumes M whle poducng P 1 _M. I also consumes powe fom eleccy EL and machne M 1 a s sa only o egeneae M 1 a s end. Soage asks come nex (I s n lgh gey). They ae moe complex han hose descbed n Caso e al. 8 snce now he soage uns ae shaed ahe han dedcaed o a sngle maeal esouce. These soage asks ae hybd conssng of wo pas: (a) Send_o_soage connuous; (b) Hold_n_soage bach. Any pa wll gge he consumpon of he coespondng desnaon soage un (e.g. S 1 ) such ha s made unavalable fo ohe poducs dung ha me. Boh can be acve smulaneously n ode fo a poduc o ene a paally flled soage un wh he same qualy. The connuous pa consumes a poduc-afe-machnes esouce (e.g. P 1 M) whle s beng poduced by Pocess_P 1 _M 1 Pocess_P 1 _M M and poduces a poduc-a-soage esouce (e.g. P 1 _S 1 ). The bach pa can be vewed as empoaly hdng he maeal esouce (e.g. P 1 _S 1 ) fo one me neval. The las subse of asks ansfe he poducs fom soage o cusomes I (dak gey). They ae assumed nsananeous meanng ha hey las much shoe when compaed o he pocessng asks. They could have also been easly defned as bach o connuous. We ae also assumng ha no equpmen esouce s nvolved whch s agan no n any way lmng Fuhe equemens fo Connuous-me Model Ideally he TN epesenaon of he pocess should be ndependen on he model ha s employed o solve he schedulng poblem. Howeve hs s no always he case. In he poblem a hand he eleccy cos paamee s me dependen and we need o make sue ha he coec paamee s consdeed. Whle n dscee-me hs s no paculaly dffcul snce he enegy level can be decly lnked o a subse of he me nevals n connuous-me he lnk s a he ask level. 8

9 In ohe wods we need o dsngush whehe he ask s execued a e 1 e 2 o e E wh subse I e c ndcang he pocessng asks pocessed a enegy level e. Fuhe dsaggegaon of he pocessng asks s hus equed (Fgue 4) leadng o a oal P M E of such asks. Neveheless he exac same esouces ae nvolved so Fgue 3 s sll vey much useful. Pocess_P 1 _M 1 _E 1 ae P 1M1 Pocess_P 1 _M 1 ae P 1M1... Pocess_P 1 _M 1 _E E ae P 1M1 Fgue 4. Due o me dependen eleccy coss pocessng asks need o be fuhe dsaggegaed fo he connuous-me model. 4. TN-based Model Enes Ths secon focus on he model enes ha ae used egadless of he me epesenaon employed by he esouce-task Newok based model. In he case of he paamees hs does no necessaly mean ha hey have he same doman o ake he same values. Addonal vaables fo he connuous-me model ae explaned followng he descpon of he undelyng me gd n secon Vaables Tasks ae geneally chaacezed by exen vaables one bnay N and one connuous se ξ. The fome denfy he sa of ask a me pon whle he lae nomally gve he amoun handled by he ask. We ae mplcly assumng ha he asks ae ehe nsananeous (sa and end a ) o connuous whch can be made o las a sngle me neval whou loss of genealy 3. In pacce hs means ha a few consecuve nsances of he ask wll have o be pefomed o mee lage demands. Bach asks ha las one me neval ndependenly of s duaon ae also ncluded. In fac he soage asks descbed n secon 3.2 ae hybd conssng of a connuous and a bach * pa ha can occu smulaneously so anohe se of connuous exen vaables ξ s equed fo 9

10 he full chaacezaon. I s also mpoan o hghlgh ha no equpmen esouce s nvolved n he ansfe asks. Thus bnay vaables ae no necessay. The beakhough of he TN 6 comes fom he unfed eamen of esouces. The mahemacal fomulaons keep ack of esouce avalably ove me hough he excess esouce vaables. The wod excess s cucal wh non-zeo values ndcang ha hee s sll some amoun of esouce avalable a even pon besde he amoun ha s allocaed o bach asks ha wll sa o nsananeous asks ha wll be execued a. Equpmen esouces ae mosly eaed as ndvduals so havng =0 ndcaes ha hee s one ask sang a ha uses un. If one wans o foce such equpmen o be dle dung hen mus be avalable n excess and he consan o use s =1. In he pesence of connuous asks anohe se of excess esouce vaables 5 end s somemes equed. These gve he excess amoun of esouce mmedaely befoe he end of neval. In hs pacula poblem hee wee wo easons fo he use: () he need o occupy a soage un gh fom he sa of a pocessng ask n ode o ensue ha hee s one avalable; () he need o know he amoun n soage o guaanee ha he maxmum capacy s no exceeded. Evens ggeed by he sa o end of asks wll affec he wo ypes of excess esouce vaables dffeenly (see secon 5.3). esouce avalably a he begnnng of he me hozon 0 s ofen known fo all esouces and so s nomally defned as a paamee. The value fo equpmen uns s gven n Eq. (1). Neveheless he scope of he mahemacal fomulaon s wde f such an eny s a vaable. Ths s paculaly useful f one wans o solve he smulaneous desgn and schedulng poblem 8 whee one of he goals s o selec he mos appopae uns fo he plan. In hs pacula case and even hough he aw-maeal equemens can easly be deemned fom he poduc demands ha ae fxed he nal avalably of he aw-maeal esouce M wll be a vaable n he model. I wll be he only npu of aw-maeal o he sysem ( s saghfowad o consde ohewse) and so wll be equal o he oal amoun equed ove he full me hozon. Fo all ohe esouces wll be se o zeo (Eq. 2). 10

11 EQ 0 = 1 (1) EQ M 0 = 0 (2) 4.2. Sucual Paamees The TN epesenaon of he pocess s bough no he model by he sucual paamees. These can be slghly dffeen dependng on whehe we ae usng a dscee o connuous-me fomulaon 14. Howeve n poblems wh asks lasng a mos a sngle me neval he dsceeme fomulaon can use he exac same paamees as fo he connuous-me 1 model. Sucual paamees gve ehe he oal esouce consumpon/poducon o he popoon elave o he amoun handled by he ask. Thee ae fve ses of sucual paamees coespondng o dffeen ypes of neacon wh he exen vaables. Dscee neacons occu ehe a he sa o end of he ask whle a connuous neacon akes place houghou he execuon of he ask. Paamees µ and μ ae he dscee neacons assocaed wh vaables N especvely fo he sa and end of he ask. They ae used wheneve he amoun of esouce consumed/poduced s ndependen on he amoun handled by he ask as fo equpmen uns. Fo maeal esouces one nomally eles on he dscee neacons assocaed wh he connuous exen vaables ξ.e. paamees ν and ν. Fnally paamees * λ hold he connuous neacons assocaed wh he connuous exen vaables ehe ξ o ξ. Whle hee can be many paamees n oal he lage majoy wll be equal o zeo and hose ha ae no wll mosly have a value of ehe 1 o -1. Afe some pacce devng hem fom he TN becomes a elavely easy ask. As an example fo he case sudy a hand we have: µ M1 =-1; μ 1; µ EL =-pw P1M1 ; λ M =-1; λ P1_M =1 - =Pocess_P 1 _M 1 (3) M1 = µ S1 =-1; μ 1; ν P1_S1 =-1; ν 1; λ P1_M =-1; λ P1_S1 =1 =Soe_P 1 _S 1 (4) S1 = P1 _ S1 = ν P1_S1 =-1; ν 1 =emove_p 1 _S 1 (5) P1 = 11

12 4.3. Ohe Paamees Whle sucual paamees lnk asks wh esouces hee ae also paamees ha ae ehe elaed o asks o o esouces. These have aleady been gven fo he pocess enes such as machnes and poducs so now we jus need o make he coespondence o he absac enes of he TN befoe assgnng he appopae values. An algohm was devsed fo hs pupose ha smply geneaes he asks and esouces fom he se of poducs soage uns machnes and eleccy levels whle defnng funcons (f x ) o keep ack of he coespondence. These may also depend on whehe he poblem s o be solved wh a connuous o dscee-me fomulaon (c o d supescp especvely). I c = ; I e c = =1; p=1 m=1 e=1 I c =I c { }; I ec =I ec { } ρ max =ρ pm =1 e=e1 e=e? YES m=m? YES p=p? YES END NO NO NO m=m1 p=p1 Fgue 5. Pa of he algohm esponsble fo geneaon of ses I c c and I e. Fgue 5 llusaes he pa of he algohm esponsble fo geneang he connuous asks (I c ) and he elemens of se I e c. As an example fo P=2 M=1 and E=3 we ge I c ={ 1 6 } and 12

13 c I e 1 ={ 1 4 } c c I e2 ={ 2 5 } and I e1 ={ 3 6 }. Along he way he asks maxmum pocessng aes ae made equal o he pocessng aes of he coespondng poduc-machne pa ρ pm. Ths can be genecally epesened by eq 6. max ρ c = f ρ p m) I (6) 1( max ρ The algohm also geneaes he dffeen demand pons dtd and he me peods ptp e fo eleccy levels e ogehe wh he fxed fx d and sang and endng mes lb ep and ub ep [h]. Boh he pofle of cos and maxmum powe avalably nfluence he numbe of me peods. Takng as an example he pofle of Fgue 2 and assumng consan maxmum powe avalably hee ae a oal of 8 1 and 10 peods fo he cheap medum and expensve cos levels especvely. Fo he dscee-me fomulaon he eleccy cos levels c e can be easly assocaed o specfc me nevals. I was assumed ha he cos can change evey hou so hee s no poblem f 1 s an nege mulple of he duaon of evey me neval δ. Fo he ohe cases a weghed mean can be used even hough he oal cos wll no be exac wheneve hee ae paally execued asks (see lae on secon 6.1). Equaon () emphaszes ha he cos fo neval ce s a funcon of δ. ce d = f ( c δ ) T () 2 e T end max mn Excess esouce vaables may be subjec o gven uppe ( ) and lowe bounds ( end mn ). The lae ae se o zeo whle one should also do he same fo he fome wheneve possble. In hs way soluon degeneacy wll be educed by ensung ha asks ae execued sequenally whou wang peods n beween. Fo nsance dong hs fo poducs locaed mmedaely afe he machnes wll foce sendng o soage he exac same amoun ha s pocessed. On he ohe hand poducs should be allowed o exs n soage vessels ( LS ) up o an amoun equal o he maxmum capacy Eq. (8). Neveheless an excess of such esouces s only empoay ehe because hey wll be ansfeed o he clen o because hey need o be hdden empoaly so ha ohe poducs do no occupy he same soage vessel. Ths s accomplshed hough Eq. (9). We ensue ha he sysem eceves he mnmum amoun of aw-maeal f hee s max none avalable a he las even pon (Eq. 10). 13

14 end max CT M 0 f3( cap ) \ T T = s LS (8) max CT FP M = 0 \ T (9) = 0 max M T (10) The sysem can exchange maeal wh s suoundngs a any me. We know exacly when a cean quany of a esouce becomes avalable o o s emoved fom he sysem. In a dsceeme model hee s a dec lnk beween eal me and he me pons of he gd so paamees ou Π and Π n can easly be deved see Eqs. (11-12). In hs case hey wll be used o povde he machnes wh he maxmum level of powe consumpon and o mee he poduc demands. When usng connuous-me fomulaons howeve he coespondence beween he dscee npus/oupus and he even pons of he me gd wll be accomplshed hough he use of wo addonal ses of bnay vaables. Thus he me ndex of he npu paamees needs o be changed o me peod p of enegy level e whle ha of he oupu paamees s changed o demand pon d (Eqs ). d UT Π = f pwx ) T (11) n 4 ( h dy T ou d FP Π = f d p dy h ) T 1 (12) 5 ( n c UT Π p e = f4 ( pwxh dy) ptp e E (13) ou c FP Π d = f5 ( d p dy h ) d TD (14) 5. New Connuous-Tme Fomulaon (CT) The new connuous-me fomulaon fo sho-em schedulng of connuous mulpoduc plans unde vaable uly avalably coss/pofles and mulple nemedae due daes s gven below. The man focus wll be on ssues elaed o me modelng and he esouce balances ove me. The model enes ae summazed n he Nomenclaue secon wh all connuous vaables beng nonnegave. 14

15 5.1. Tme epesenaon The poposed connuous-me fomulaon uses a sngle me gd o keep ack of evens akng place see Fgue 6. I uses T even pons also named global me pons 15 ha can be placed anywhee beween he ogn and end (H) of he me hozon. Tadonally asks sang a even pon have been assumed o sa a he absolue me deemned fo ha even pon T. Fo bach asks wh duaon shoe han T 1 -T s assumed ha he oupu esouces say n he coespondng equpmen un fo he emanng me. Wh connuous asks s assumed ha hey las exacly he duaon of he me neval whle beng pocessed a a ae lowe han he pedefned maxmum ae. In ealy howeve he asks have oal feedom o sa anywhee povded ha hey end befoe T 1. Gménez e al. 16 hghlgh ha he numbe of even pons needed o epesen a soluon can be educed f asks ae no equed o sa (end) exacly a a me pon whle pesenng a new connuous-me fomulaon fo he sho-em schedulng of mulpupose bach plans. The cuen case sudy s a pefec example fo he advanages of such an appoach as wll be explaned nex. Slo 1 Slo 2 Slo T-2 Slo T T-2 T-1 1 T 0 H Fgue 6. Connuous-me epesenaon. Consde a smple example nvolvng he execuon of a sngle pocessng ask a slo/neval see Fgue. The neval boundaes T and T 1 can esul fom wo consecuve demand pons n cases whee he duaon of he ask s lowe han he duaon of he lowes cos peod locaed nsde he neval (shown n geen). Vaable Ts wll gve he sang me of he ask execued dung neval (o he eales sang me amongs all asks execued). Noe ha no fuhe even pons ae equed. Hgh cos Low cos Pocessngask Medum cos T ub e2p-1 ub e3p Ts ub e1p T 1 lb e 3p lb e 1p lb e 2p Fgue. Connuous asks execued a neval do no necessaly sa a even pon. 15

16 In he moe geneal case he amoun o poduce wll eque an exenson lage han he duaon of he cheapes me peod. In such case he pocessng ask wll be dvded by he solve no as many mes as he numbe of dffeen enegy cos levels whch may nclude consecuve me peods of he same enegy level. Each such dvson wll eque an addonal even pon as can be seen n Fgue 8. Medum Hgh cos Low cos Medum cos P. Task Pocessngask Pocessngask T ub e2p-1 T 1 ub e3p T 2 ub e1p T 3 Ts lb e 3p Ts 1 lb e 1p Ts 2 lb e 2p Fgue 8. Mulple ask nsances mgh be equed o mee he demands one fo evey acve enegy cos level Tmng Consans The fundamenal mng consan fo a sngle me gd fomulaon 3 saes ha he dffeence n me beween wo consecuve even pons mus be geae han he duaon of he ask akng place. Equaon (15) s wen fo all equpmen esouces nvolved n he mng consans. The summaon s used o educe he negaly gap 3 and mplcly assumes ha hee can only be one ask execued n such equpmen un a a cean me. Ths s ensued by he nal esouce avalably and he excess esouce balances (see Eqs. (1) and (2)). μ ξ TC T 1 T T T max c ρ I (15) A smla consan can be used o guaanee ha he asks ae fully execued whn me neval. Equaon (16) ogehe wh Eq. (1) also sasfy Eq. (15). Howeve compuaonal sudes have shown ha s bee o keep Eq. (15). μ ξ TC T 1 Ts T T max c ρ I (16) Ts T T T (1) We now defne a new bnay vaable Y pe whch akes he value of one f dung neval asks ae pocessed whn me peod p of enegy level e. Smlaly we defne bnay vaable ou Y d o denfy whehe o no even pon coesponds o demand pon d. I can be assumed whou loss 16

17 of genealy ha hee s: () one acve me peod of an enegy level dung ; () one even pon assocaed o due dae d eqs p e ee ptpe Y = 1 T T (18) d T 1 ou Y = 1 d TD (19) I me neval s locaed whn me peod p of enegy level e hen he sang me of asks mus be geae han he me peod lowe bound Eq. (20). Lkewse hey mus end befoe s uppe bounday (see Fgue 8). Equaon (21) s a bg-m consan ha s only acve f hee s a connuous ask beng execued ha belongs o enegy level e. Ths consan s he one dffeenang beween enegy levels fo equvalen asks (efe o Fgue 4) so ha he pope eleccy cos s accouned fo n he objecve funcon see Eq. (33). Ts ee ptp lbe py p e T T e (20) Ts μ ξ ube py p e H ( 1 μ N ) max c ρ c I e ptpe I e TC T T e E (21) If even pon coesponds o demand pon d hen he me values mus mach as expessed n Eqs. (22)-(23). ou T fxdy d T 1 (22) dtd ou ou T fxdy d H 1 Y d ) T 1 (23) dtd ( dtd Fnally hee ae geneal bounds on he mng vaables. Accodng o Fgue 6 T [0H] so Eqs. (-25) apply. Equaon (25) enfoces ha no ask sas a he las even pon. Snce he model assumes ha asks can sa afe T he me of he fs even pon can be se o zeo whou loss of genealy (Eq. 26). T H T () Ts H T T (25) T 1 = 0 (26) 1

18 5.3. Excess esouce Balances esouce avalably ove he me gd s managed by he excess esouce balances. These ae mulpeod maeal balance expessons n whch he excess amoun a even pon s equal o ha a he pevous even pon (-1) adjused by he amouns dsceely o connuously poduced/consumed by all asks sang o endng a. Fgue 9 llusaes he conbuons o he values of he vaables and suppos he explanaon of Eqs. (2-28). In Eq. (2) he nal esouce avalably (fs em on he HS) s only o be consdeed a he fs even pon =1. Fom Eqs. (1-2) we know ha wll be dffeen fom zeo only fo a subse of he esouces 1 fo EQ and a value o be deemned by he solve 0 fo M. Fo he emanng pons he value fom he pevous even pon s used nsead. Fo connuous esouces CT he ones nvolved n Eq. (28) eceves he conbuon of end 1 whch gves he excess amoun mmedaely befoe he end of neval -1 (even pon ). Fo he ohe esouces bu he ules one should use he value a he begnnng of he pevous neval -1. Ths hd em on he HS s no used fo UT snce we wan o peven esouce avalably o popagae fom one even pon o he nex. Ths s because hee wll be npus fom he exeo a evey even pon (sxh em on he HS). EQ M UT CT CT Π n -1 npu maeal -pw I I s I c -1 I c oupu oupu npu oupu npu EQ -1 I c I s 1 EQ CT -Π ou npu FP end Fgue 9. Evens affecng he value of he excess esouce vaables. 18

19 19 T Y Y v N v N TD d ou d ou d T E e TP p e p n e p I I T end FP UT e UT CT CT Π Π = = ) ( ) ( ) ( 1 1 ) ( ξ μ ξ μ (2) ) ( * T T v CT I I end s c = ξ λ ξ ξ λ (28) In Eq. (2) he fouh and ffh ems deal wh mos of he dscee neacons of asks wh esouces. The excess amoun a even pon deceases by asks sang a and nceases by asks sang a -1 (endng a ) o (f hey ae nsananeous I ). The dscee neacon a he end of soage asks (I s ) affecs excess vaable end ahe han 1 see las em on he HS of Eq. (28). The pupose s o ge he amoun n soage fo evey u so ha he maxmum capacy consan can be enfoced. Vaable end also deals wh he connuous neacons of connuous (I c ) and soage asks. Fnally he las em n Eq. (2) emoves poduc esouces fom he sysem wheneve he me of even pon s equal o a demand pon d. The excess esouce vaables should le whn gven uppe and lowe bounds. The geneal consans ae gven n Eqs. (29-30) whle he values fo hs case sudy ae dscussed n secon 4.3 and gven n Eqs. (8-10). T max mn (29) max mn T T end end end (30) 5.4. Ohe Consans The emanng se of consans ensues ha he values of he connuous exen vaables ae se o zeo wheneve he ask s no execued.e. when he bnay exen vaable s equal o zeo. The uppe bound n Eq. (31) U can be gven by he poduc of he me hozon mes he sum of he maxmum pocessng ae on evey machne see Eq. (32). * T T I I N U s c I s ξ ξ (31) = M m m p P p H U max ρ (32)

20 5.5. Objecve Funcon The mahemacal fomulaon s compleed wh he objecve funcon. Equaon (33) maxmzes he oal eleccy cos [k ]. Ths s gven by he sum ove all eleccy levels e me nevals and asks of he poduc of eleccy cos c e [ /kwh] powe consumpon [MW] and duaon of he ask [h]. mn ee Ic UT e T T ξ ce ( μ ) (33) max ρ 6. Dscee-Tme fomulaon (DT) Exenal neacons wh he sysem a specfc pons n me such as dffeen eleccy coss and due daes ae handled much moe easly by a dscee-me fomulaon. Ths of couse povded ha he locaon of some of he gd s me pons s whn accepable accuacy of he exenal evens mng. Afe ecallng a few mpoan aspecs elaed o he epesenaon of me he model consans ae gven. These ae analogous o he non-mng consans of he connuousme fomulaon whch s no supsng gven he fac ha he dscee-me fomulaon s n fac a sub-model of he moe geneal connuous me fomulaon Tme epesenaon The sngle me gd employed by he dscee-me fomulaon has equal lengh me nevals of duaon δ see Fgue 10. As a consequence he absolue me of all even pons s known a po. Fxed lengh me nevals have mpoan consequences fo boh bach and connuous asks. Fo he fome he duaon mus ofen be ounded o a mulple of he neval lengh leadng o he consdeaon of an appoxmaed veson of he eal poblem. Fo connuous asks he exen can be made o vay whn [0 δ ρ max ] o ensue ha all possble demands can be me. Howeve hee may be nevals whee he ask s beng pefomed below s maxmum pocessng ae whch n vey consaned poblems may compomse opmaly o even feasbly snce no ohe ask can ake advanage of he emanng capacy of he un. Ths can be llusaed wh a smple example. 20

21 Slo1 Slo 2 Slo 3 Slo T-2 Slo T-1 δ T-2 T-1 1 T 0 H Fgue 10. Dscee-me epesenaon. Demand 20 Demand 35 Demand h Idle Subopmal 5 Idle Opmal Fgue 11. Lmaon of dscee-me fomulaon when handlng connuous asks. In Fgue 11 can be seen ha fo δ=1 h we need o occupy he pocessng un fo he fs egh hous of he me hozon o mee he blue and ed poducs demands. The poblem s ha he las one s a medum-cos neval. Even hough he objecve funcon accouns only fo he half an hou ha he un s acve opmaly s compomsed snce he demands of boh poducs can be me enely n he low-cos peod (fs seven hous). Neveheless such soluon canno be obaned snce he excess esouce balances allow poducon of a mos one poduc on each me neval. Howeve can be acheved wh δ=0.5 h (o wh a connuous-me fomulaon) fo whch he nevals lengh wll mach exacly he me equed o poduce 5 of poduc. Thus he ed poduc can sa o be poduced mmedaely afe he blue ends n ode o fnsh befoe he medum-cos me peod sas. Moe mpoanly would sll be possble o poduc anohe 10 on befoe he due dae n case of hghe demands Model Consans The excess esouce balances ae gven by Eqs. (28) and (34). The lae dffes fom Eq. (2) n he ems nvolvng paamees Π n and Π ou. Now he doman s exacly equal o ha of he equaon (ecall he las paagaph of secon 4.3) so hee s no need o employ addonal bnay vaables. 21

22 22 T v N v N ou T n I I T end FP UT UT CT CT Π Π = = ) ( 1 1 ) ( ξ μ ξ μ (34) Snce hee ae no mng consans he capacy consans need o ensue ha he maxmum amoun of maeal handled by connuous asks does no exceed he maxmum pocessng ae mes he duaon of he me neval. Fo soage asks Eq. (35) s dencal o Eq. (31). ) ( max * T T I I N U s c I I I c s s ρ δ ξ ξ (35) 6.3. Objecve Funcon The dffeences n compason o he objecve funcon of he connuous-me fomulaon gven n Eq. (33) esul fom he eleccy cos paamees whch have a dffeen doman. Because of hs Eq. (36) has one less summaon. c UT I T T ce max ) ( mn ρ ξ μ (36). Compuaonal Sudes The pefomance of he models s llusaed hough he soluon of 10 es cases wh daa geneaed andomly fom a eal ndusal poblem. In all of hem he poduc demands lead o plan opeaon below he maxmum capacy. Due o he dffeences n eleccy cos among he enegy levels assgnng poducon o he lowe levels wll have he bgges mpac on he oal cos. As we appoach full plan capacy all he enegy levels become acve and he solve swches focus o fndng he bes combnaon of poduc-machne assgnmen accodng o he pocessng aes and powe needs whch ae poduc dependen. Wh he excepon of EX5 he machnes ae occuped 41 o 56 % of he me. Decly elaed o he complexy of he poblem s he numbe of poducs machnes and soage uns so he poblems wll be manly chaacezed by (P M S). The mahemacal fomulaons gve se o mxed-nege lnea pogammng (MILP) poblems. They wee mplemened n GAMS 22.8 ogehe wh he algohm ha conves he poblem daa no he TN foma. All poblems wee solved by CPLEX 11.1 wh defaul opons up o a elave

23 opmaly oleance=10-6 unless ohewse saed. The hadwae conssed on a lapop wh an Inel Coe2 Duo T9300 pocesso unnng a 2.5 GHz 4 GB of AM and unnng Wndows Vsa Enepse..1. Illusave Example To llusae he capables of he new connuous-me fomulaon we sa wh a vey smple example nvolvng wo poducs one machne and one soage un (211). The machne powe equemens s n hs case poduc ndependen and equal o pw M1 = 5 [MW] and hee ae no consans on maxmum powe consumpon. The pocessng aes ae ρ P1M1 = 80 ρ P2M1 =0 [on/h] soage capacy s cap S1 =1200 [on] and he poduc demands can be found n Table 1. Table 1. Poduc Demands a he End of he Day fo Examples EX1-EX3 [on] MO ( h) TU (48 h) WE (2 h) TH (96 h) F (120 h) SA (144 h) SU (168 h) P P The opmal soluon fo EX1 coesponds o a oal eleccy cos of The esulng schedule and soage pofles ae shown n Fgue 12 whee he numbes nsde he boxes ndcae he ask lengh [h]. I eques a oal of T=11 even pons whch s he mnmum numbe ha ensues feasbly. No mpovemen n cos s obseved fo T=12. Such behavo s no common n connuous-me fomulaons 3411 whee ypcally he objecve funcon nceases wh T a leas a couple of mes. Theefoe s an ndcaon ha he poblem s hghly consaned whch s a dec consequence of: () lage vaably of he eleccy cos leadng o a oal of 35 me peods wh consan pce; () he asks no beng allowed o coss me peods. Thus some demands can only be me by execung mulple nsances of a pacula ask. Moe specfcally he Monday (MO) and Sauday (SA) demands of P1 need wo asks whle he Fday (F) demand of P2 eques hee. In Fgue 12 8 ou of 10 nsances ae execued n low-cos peods (geen egon) wh he emanng beng n medum-cos levels (yellow egon). Ideally all should be execued n he geen egon bu ha would only be possble fo opeaon well below he maxmum capacy of he plan. Powe shoages may foce asks o swch o a hghe cos level. In ode o show ha he model can 23

24 cope wh consans of hs ype EX2 was solved. I s uses he same daa of EX1 bu now he maxmum powe consumpon dung he fs seven hous of Tuesday and Thusday (pwx 0-6TU ; pwx 0-6TH ) s se o zeo. The opmal soluon fo EX2 s shown n Fgue 13 and coesponds o a oal cos of When compaed o EX1 anslaes no a 15.8% ncease n cos whch s vey sgnfcan consdeng such small changes n daa. Such a schedule s obvously a feasble soluon o EX1 whch s n many aspecs smla o he pevous one. Hence also seves o llusae he mpac ha advanced schedulng ools lke he models poposed n hs pape can have on plan pofably. We wll eun o hs dscusson n secon.3. M S S1 Tme (h) P1 200 P2 0 (on) Tme (h) Fgue 12. Opmal soluon fo EX1. M S S1 Tme (h) P1 200 P2 0 (on) Tme (h) Fgue 13. Opmal soluon fo EX2.

25 M S S S1 P1 Tme (h) 1000 P (on) Tme (h) S (on) Tme (h) P1 P2 Fgue 14. Opmal soluon fo EX3. The amoun of avalable soage s ofen lmng. The pevous examples feaued a shaed soage un ha could only handle one poduc a a me. Ths made mpossble o mee demands of dffeen poducs n he same day whch explans he foma of such daa. Moe mpoanly poducon of p could only sa afe he oal amoun poduced of p had been dspached. In ode o allow fo moe flexbly EX3 consdes wo dencal soage uns and no maxmum powe consans. Wh he addonal soage un he oal cos can be educed o (2.5%). The opmal soluon s gven n Fgue 14. I bascally moves all he poducon of P2 (n ed) o low-cos peods whle moe P1 (blue) ges deced o medum-cos peods. Whle one mgh be emped o explan hs behavo based on he dffeen pocessng aes hee exss a degeneae soluon ha feaues P2 n boh low- and medum-cos peods. On Monday he poducon of P1 has nceased fom 1000 o 1089 on so slo S2 wll eman paally flled unl he end of Wednesday. S2 also holds P2 on Fday even hough could all be allocaed o S1. 25

26 .2. Compuaonal Pefomance The compuaonal esuls of he 10 examples ae gven n Table 2. The example poblems can be dvded no wo levels of complexy. EX1-5 can be addessed by he connuous-me model (CT) whle fo he ohes one has o ely on he dscee-me fomulaon (DT). Ths s he fs majo esul: he connuous-me fomulaon s lmed o small poblems. Even fo EX5 whch feaues hee poducs wo machnes and wo soage uns we had o specfy low poduc demands (leadng o a 33% capacy) o ensue ha mos asks could f no a sngle peod of consan eleccy cos. In hs way a good soluon could sll be found wh a elavely small numbe of even pons. Howeve ook aleady moe han 1 hou o pove opmaly fo T=11 ( 26911) see Table 2 a value ha s sll 0.5% above he one fom DT. Fo lage poblems CT s nacable due o he followng facs: () he compuaonal effo s songly dependen on T wh expeence ellng us ha we ge ypcally a one ode of magnude ncease fo a sngle ncease n T; () s no saghfowad o fnd a numbe ha ensues feasbly; () may be even dffcul fo he solve o fnd ou ha he value of T s nsuffcen.e. fndng ha he poblem s nfeasble. The las examples have moe o less he complexy one expecs o fnd n a eal suaon. The encouagng esul s ha he dscee-me fomulaon wh 1-hou nevals can fnd vey good soluons o he poblem n sho compuaonal me. In fac fom he esuls n Table 2 one can see ha whn up o 5 mnues he opmal/bes soluons euned ae whn a elave opmaly gap of 0.8 %. Thus DT povdes an effcen decson-makng ool fo hs poblem. The gaps could be educed mos of he mes by up o one hou of compuaonal me essenally due o he fndng of bee soluons snce he value of he elaxed poblem emaned vually unchanged. Ths s an ndcaon of a hgh degee of degeneacy (moe on hs n secon.4). The esuls fom Table 2 also show ha he dscee-me fomulaon s sgnfcanly ghe han he connuous-me one. In EX1-4 and EX he negaly gap s even zeo and wh he excepon of EX8 dd no exceed 0.2%. The subsanally lage negaly gap fo he lae (3.3%) s pobably due o he use of hashe consans on maxmum powe avalably ha pevened he 26

27 machnes o opeae smulaneously n a sgnfcan pa of he week. Neveheless CT sll exhbed a easonable negaly gap (maxmum=23% fo EX2 mnmum=4.8% fo EX5 aveage=11.6%). Fnally was supsng o fnd ou ha he numbe of bnay vaables of DT was of he same ode of magnude as ha of CT n spe of he lage dffeence beween he T numbes used n he wo me gds. ecall ha T s one of he ndces of he bnay exen vaables N. Thee ae wo easons fo hs: (a) hee ae hee mes as many asks of he connuous ype (I c ) fo CT due o he enegy levels dsaggegaon see secon 3.2.1; bu moe mpoanly (b) CT employs wo ou addonal ses of bnay vaables Y and Y pe wh he lae feaung hee ndces and beng he d majo conbuo wheneve hee s a lage vaably on eleccy cos. Takng EX4 (T=12) as an example such vaables ae esponsble fo 69% of he oal numbe. Table 2. Compuaonal esuls (P Model T bnay sngle consans MIP MIP CPU nodes M S) vaables vaables ( ) ( ) s EX1 (211) DT CT CT EX2 (211) DT CT CT EX3 (212) DT CT CT EX4 (212) DT CT CT EX5 (322) DT a 1332 CT CT CT EX6 (323) DT b EX (334) DT EX8 (335) DT c EX9 (434) DT d 400 EX10 (534) DT e 500 a Inegaly gap a 300s (IG300) [%]=0.06; Inegaly gap a me of emnaon (IGT) [%]=0.04; b IG300=IGT=0.02; c IG300=0.61 IGT=0.18; d IG300=IGT=0.06; e IG300=0.8 IGT=

28 .3. Poenal Cos Savngs One mpoan queson s wha ae he poenal savngs ha one may ge by usng such a compehensve opmzaon schedulng ool? In ohe wods how mpoan s o accoun fo he eal eleccy cos pofle when devng he schedule? Gven: () he many consans n ems of soage avalably demand sasfacon and maxmum powe avalably; () he need o povde a schedule fas no only fo nomal opeaon bu moe mpoanly as a esponse o unpedced evens such as machne beakdown o ugen odes ha need o be sasfed; I s fa o assume ha a feasble soluon o he poblem povdes a vald bass fo compason. A so called blnd schedule can easly be deved by he dscee-me fomulaon unde he assumpon of consan eleccy pce houghou he me hozon. Unde such scenao and accodng o Eq. (36) he model wll smply y o poduce he poducs n he fase machnes. If n ems of pocessng ae we ge he same poduc ank on evey machne hen he schedule wll end o a sees poducon ahe han a paallel one whch as we can see n Fgue 15 s he pefeed choce unde a vaable cos pofle. Neveheless he nemedae poduc demands paally avod hs. Also he effec becomes less mpoan fo hghe plan capaces. Keepng hs n mnd we sck o he plan. Afe obanng he blnd schedule he eal cos paamees ae used o compue he ue oal cos. Table 3. Poenal Cos Savngs as a Funcon of he Plan Capacy. EX1 EX2 EX3 EX4 EX5 EX6 EX EX8 EX9 EX10 Capacy (%) Savngs (%) Table 3 lss he esuls obaned. Even afe consdeng ha an expeenced schedule by means of he heuscs n use a he plan could do a sgnfcanly bee job han he blnd schedule fo he easons menoned above hese ae vey sgnfcan cos savngs. The end s ha as he plan appoaches full capacy he poenal savngs dmnsh whch s no supsng gven ha hee ae fewe degees of feedom. I s elevan o hghlgh he esul fo EX3 when compaed o ha fo EX2. As dscussed n secon.1 hey have he same daa bu EX3 has one moe soage un whch consdeably nceases he plan flexbly. Ths makes much ease o geneae a feasble 28

29 schedule bu f one s no caeful one may make wong decsons concenng he mpac on cos (ncease n poenal savngs fom 23 o 39%). M M S S S1 P1 Tme (h) 1000 P2 P (on) Tme (h) S2 P P2 P (on) Tme (h) Fgue 15. Bes soluon found by connuous-me fomulaon fo EX5..4. Fnal emaks To end he analyss we dscuss he ssue of soluon degeneacy and s mpac on he fom of he schedule. Fo a schedule o be mplemened n pacce should have as few evens as possble snce hee ae always facos ha ae no consdeed n he model. In hs case sudy he ssue of sequence dependen changeoves s no elevan. Howeve ha does no mean ha we ae allowed o change poducs on a gven machne evey hou! Theefoe we wsh o fs complee he poducon of a poduc and only hen swch o anohe one. The ealy s hus much close o he concep of evens used by me gd based connuous-me fomulaons o of pecedence used by sequencng based models 15. The eave pocedue used by CT of nceasng he numbe of even pons one by one mplcly ensues smplcy. In conas n DT a much lage numbe of evens ae allowed so hee wll ypcally be many moe ask nsances beng execued o acheve a pacula amoun of poduc. Thus hee wll be a lage numbe of soluons ha coespond o he same value of he objecve funcon wh he lage majoy of hem beng undesable due o many changeoves see 29

30 Fgue 16. To ovecome such soluons a pos-pocessng pocedue s equed ha emoves he supefluous bu hs s beyond he scope of hs pape Unaccepable OK Fgue 16. Hgh soluon degeneacy may lead o an unaccepable numbe of changeoves n soluons fom he dscee-me fomulaon. 8. Conclusons Ths pape has focused on he modelng of dscee evens ha occu a pedeemned pons n me wh a connuous-me schedulng fomulaon. These ncluded mulple nemedae due daes uly avalably and vaable eleccy coss. A concepually new geneal model has been poposed ha eles on he esouce-task Newok fo pocess epesenaon. I s neveheless lmed o connuous asks. The valdy of he appoach has been demonsaed on a few es cases adaped fom a eal ndusal poblem. In such pocess soage uns ae used fo fnal poducs and may ac as shaed soage uns. To effcenly model hs ssue a novel hybd bach/connuous ask has been poposed. Despe he majo modelng beakhough he esuls have shown ha only poblems of small sze can be handled effecvely. Ths uned ou aenon o dscee-me models whee dscee evens can be handled n a moe naual and saghfowad way. The only doub was o whehe a suffcenly fne me gd could be used o epesen he poblem daa accuaely. The goal was o ensue a weekly schedule ha accouned fo houly changes n eleccy cos and hs has been successfully accomplshed. Poblems of ndusal sgnfcance wee ackled down o opmaly gaps below 1% wh paccal compuaonal mes of 5 mnues. Ths s mosly due o he chaacescally low negaly gap of TN dscee-me fomulaons whch wee found o be sgnfcanly lowe han he connuous-me counepas. The las pa of he pape has hghlghed he mpoance of akng vaable eleccy coss no consdeaon when devng he schedule. Sae-of-he-a schedulng fomulaons have he poenal o acheve majo savngs when compaed o pocedues ha ae mosly focused on 30

31 feasbly. Whle accuae values ae obvously dependen on poblem daa paculaly on he dffeen cos levels ageed wh he eleccy povde and on he schedulng pacce a he plan esuls have shown poenal cos savngs aound 20%. Clealy values of hs ode of magnude povde enough movaon fo he ncopoaon of he dscee-me mahemacal fomulaon pesened n hs pape as a cenal elemen of he decson makng ool used by ndusy. Acknowledgmens The auhos gaefully acknowledge fnancal suppo fom Fundação Luso-Amecana and he Cene fo Advanced Pocess Decson-makng a Canege Mellon Unvesy. Nomenclaue Ses/Indces DY/dy=days of he week E/e=Eleccy cos levels H/h=hous of he day I/=asks I c =connuous asks c I e =connuous asks execued dung enegy level e I s =soage asks I =nsananeous ansfe asks M/m=machnes P/p=poducs /=esouces CT =esouces connuously poduced/consumed EQ =equpmen esouces FP =fnal poduc esouces LM =poducs a a locaon afe he machnes LS =poducs a a locaon nsde he soage uns M =aw-maeal esouces TC =equpmen esouces nvolved n mng consans UT =uly esouces S/s=soage uns T/=even pons TD/d=demand pons TP e /p=me peods of eleccy level e Paamees c e =eleccy cos a level e [ /kwh] cap s =capacy of soage un s [on] ce =eleccy cos dung me neval [ /kwh] d pdyh =demand of poduc p a he end of hou h of day dy [on] H= me hozon [h] lb ep =sang me of me peod p of eleccy level e [h] 31

32 pw pm =powe equemen fo poduc p n machne m [MW] pwx hdy =maxmum powe consumpon dung hou h of day dy [MW] max =uppe bound on avalably of esouce a even pons end max = uppe bound on avalably of esouce dung nevals mn = lowe bound on avalably of esouce a even pons end mn = lowe bound on avalably of esouce dung nevals fx d =absolue me of demand pon d [h] ub ep =endng me of me peod p of eleccy level e [h] δ=duaon of evey neval on dscee-me gd [h] λ =connuous neacon of esouce dung execuon of ask µ =dscee neacon of esouce wh ask a s sa acng on bnay exen vaables μ =dscee neacon of esouce wh ask a s end acng on bnay exen vaables ν = dscee neacon of esouce wh ask a s sa acng on connuous exen vaables ν = dscee neacon of esouce wh ask a s end acng on connuous exen vaables Π Π n ou =amoun eceved no he sysem of esouce a even pon =amoun emoved fom he sysem of esouce a even pon max ρ =maxmum pocessng ae of ask [on/h] ρ pm =pocessng ae of poduc p n machne m [on/h] Vaables N =execuon of ask dung neval (bnay exen vaables) =excess amoun of esouce a even pon = nal avalably of esouce (can be a paamee fo some esouces) 0 end =excess amoun of esouce mmedaely befoe he end of neval T =absolue me of even pon [h] Ts =sang me of asks execued dung neval [h] Y pe =bnay vaable denfyng f dung neval asks ae execued whn peod p of level e ou Y d =bnay vaable denfyng f even pon coesponds o demand pon d ξ =amoun handled by ask a even pon/dung neval (connuous exen vaables) [on] ξ =amoun connuously sen o soage by ask dung neval [on] * 32

33 efeences (1) Casagnol D.; Gallesey E. Techncal epo CH-D ABB Copoae eseach Däwl Swzeland (2) Maavelas C.; Gossmann I. A Geneal Connuous Sae Task Newok Fomulaon fo Sho Tem Schedulng of Mulpupose Bach Plans wh Due Daes. Poceedngs ESCAPE-13 (Eds. A. Kaslawsk and I. Tuunen) 2003 pp (3) Caso P.; Babosa-Póvoa A.; Maos H.; Novas A. Smple Connuous-me Fomulaon fo Sho-Tem Schedulng of Bach and Connuous Pocesses. Ind. Eng. Chem. es (4) Shak M.; Janak S.; Floudas C. Connuous-Tme Models fo Sho-Tem Schedulng of Mulpupose Bach Plans: A Compaave Sudy. Ind. Eng. Chem. es (5) Schllng G.; Paneldes C. Opmal Peodc Schedulng of Mulpupose Plans. Compu. Chem. Eng (6) Paneldes C.C. Unfed Famewoks fo he Opmal Pocess Plannng and Schedulng. In Poceedngs of he Second Confeence on Foundaons of Compue Aded Opeaons; Cache Publcaons: New Yok 1994; pp 253. () Duae B.; Sanos L.; Maano J. Opmal szng schedulng and shf polcy of he gndng secon of a ceamc le plan. Compues Opeaons eseach (8) Caso P.; Babosa-Póvoa A.; Novas A. Smulaneous Desgn and Schedulng of Mulpupose Plans Usng esouce Task Newok Based Connuous-Tme Fomulaons. Ind. Eng. Chem. es (9) Maavelas C.T.; Gossmann I.E. New Geneal Connuous-Tme Sae-Task Newok Fomulaon fo Sho-Tem Schedulng of Mulpupose Bach Plans. Ind. Eng. Chem. es (10) Janak S.L.; Ln X.; Floudas C.A. Enhanced Connuous-Tme Un-Specfc Even-Based Fomulaon fo Sho-Tem Schedulng of Mulpupose Bach Pocesses: esouce Consans and Mxed Soage Polces. Ind. Eng. Chem. es (11) Shak M.; Floudas C. Novel Unfed Modelng Appoach fo Sho-Tem Schedulng. Submed o Ind. Eng. Chem. es. 33

5-1. We apply Newton s second law (specifically, Eq. 5-2). F = ma = ma sin 20.0 = 1.0 kg 2.00 m/s sin 20.0 = 0.684N. ( ) ( )

5-1. We apply Newton s second law (specifically, Eq. 5-2). F = ma = ma sin 20.0 = 1.0 kg 2.00 m/s sin 20.0 = 0.684N. ( ) ( ) 5-1. We apply Newon s second law (specfcally, Eq. 5-). (a) We fnd he componen of he foce s ( ) ( ) F = ma = ma cos 0.0 = 1.00kg.00m/s cos 0.0 = 1.88N. (b) The y componen of he foce s ( ) ( ) F = ma = ma

More information

1 Constant Real Rate C 1

1 Constant Real Rate C 1 Consan Real Rae. Real Rae of Inees Suppose you ae equally happy wh uns of he consumpon good oday o 5 uns of he consumpon good n peod s me. C 5 Tha means you ll be pepaed o gve up uns oday n eun fo 5 uns

More information

Handling Fuzzy Constraints in Flow Shop Problem

Handling Fuzzy Constraints in Flow Shop Problem Handlng Fuzzy Consans n Flow Shop Poblem Xueyan Song and Sanja Peovc School of Compue Scence & IT, Unvesy of Nongham, UK E-mal: {s sp}@cs.no.ac.uk Absac In hs pape, we pesen an appoach o deal wh fuzzy

More information

Name of the Student:

Name of the Student: Engneeng Mahemacs 05 SUBJEC NAME : Pobably & Random Pocess SUBJEC CODE : MA645 MAERIAL NAME : Fomula Maeal MAERIAL CODE : JM08AM007 REGULAION : R03 UPDAED ON : Febuay 05 (Scan he above QR code fo he dec

More information

FIRMS IN THE TWO-PERIOD FRAMEWORK (CONTINUED)

FIRMS IN THE TWO-PERIOD FRAMEWORK (CONTINUED) FIRMS IN THE TWO-ERIO FRAMEWORK (CONTINUE) OCTOBER 26, 2 Model Sucue BASICS Tmelne of evens Sa of economc plannng hozon End of economc plannng hozon Noaon : capal used fo poducon n peod (decded upon n

More information

Modern Energy Functional for Nuclei and Nuclear Matter. By: Alberto Hinojosa, Texas A&M University REU Cyclotron 2008 Mentor: Dr.

Modern Energy Functional for Nuclei and Nuclear Matter. By: Alberto Hinojosa, Texas A&M University REU Cyclotron 2008 Mentor: Dr. Moden Enegy Funconal fo Nucle and Nuclea Mae By: lbeo noosa Teas &M Unvesy REU Cycloon 008 Meno: D. Shalom Shlomo Oulne. Inoducon.. The many-body poblem and he aee-fock mehod. 3. Skyme neacon. 4. aee-fock

More information

Integer Programming Models for Decision Making of. Order Entry Stage in Make to Order Companies 1. INTRODUCTION

Integer Programming Models for Decision Making of. Order Entry Stage in Make to Order Companies 1. INTRODUCTION Inege Pogammng Models fo Decson Makng of Ode Eny Sage n Make o Ode Companes Mahendawah ER, Rully Soelaman and Rzal Safan Depamen of Infomaon Sysems Depamen of Infomacs Engneeng Insu eknolog Sepuluh Nopembe,

More information

( ) ( )) ' j, k. These restrictions in turn imply a corresponding set of sample moment conditions:

( ) ( )) ' j, k. These restrictions in turn imply a corresponding set of sample moment conditions: esng he Random Walk Hypohess If changes n a sees P ae uncoelaed, hen he followng escons hold: va + va ( cov, 0 k 0 whee P P. k hese escons n un mply a coespondng se of sample momen condons: g µ + µ (,,

More information

p E p E d ( ) , we have: [ ] [ ] [ ] Using the law of iterated expectations, we have:

p E p E d ( ) , we have: [ ] [ ] [ ] Using the law of iterated expectations, we have: Poblem Se #3 Soluons Couse 4.454 Maco IV TA: Todd Gomley, gomley@m.edu sbued: Novembe 23, 2004 Ths poblem se does no need o be uned n Queson #: Sock Pces, vdends and Bubbles Assume you ae n an economy

More information

s = rθ Chapter 10: Rotation 10.1: What is physics?

s = rθ Chapter 10: Rotation 10.1: What is physics? Chape : oaon Angula poson, velocy, acceleaon Consan angula acceleaon Angula and lnea quanes oaonal knec enegy oaonal nea Toque Newon s nd law o oaon Wok and oaonal knec enegy.: Wha s physcs? In pevous

More information

CHAPTER 10: LINEAR DISCRIMINATION

CHAPTER 10: LINEAR DISCRIMINATION HAPER : LINEAR DISRIMINAION Dscmnan-based lassfcaon 3 In classfcaon h K classes ( k ) We defned dsmnan funcon g () = K hen gven an es eample e chose (pedced) s class label as f g () as he mamum among g

More information

to Assess Climate Change Mitigation International Energy Workshop, Paris, June 2013

to Assess Climate Change Mitigation International Energy Workshop, Paris, June 2013 Decomposng he Global TIAM-Maco Maco Model o Assess Clmae Change Mgaon Inenaonal Enegy Wokshop Pas June 2013 Socaes Kypeos (PSI) & An Lehla (VTT) 2 Pesenaon Oulne The global ETSAP-TIAM PE model and he Maco

More information

Lecture 5. Plane Wave Reflection and Transmission

Lecture 5. Plane Wave Reflection and Transmission Lecue 5 Plane Wave Reflecon and Tansmsson Incden wave: 1z E ( z) xˆ E (0) e 1 H ( z) yˆ E (0) e 1 Nomal Incdence (Revew) z 1 (,, ) E H S y (,, ) 1 1 1 Refleced wave: 1z E ( z) xˆ E E (0) e S H 1 1z H (

More information

University of California, Davis Date: June xx, PRELIMINARY EXAMINATION FOR THE Ph.D. DEGREE ANSWER KEY

University of California, Davis Date: June xx, PRELIMINARY EXAMINATION FOR THE Ph.D. DEGREE ANSWER KEY Unvesy of Calfona, Davs Dae: June xx, 009 Depamen of Economcs Tme: 5 hous Mcoeconomcs Readng Tme: 0 mnues PRELIMIARY EXAMIATIO FOR THE Ph.D. DEGREE Pa I ASWER KEY Ia) Thee ae goods. Good s lesue, measued

More information

I-POLYA PROCESS AND APPLICATIONS Leda D. Minkova

I-POLYA PROCESS AND APPLICATIONS Leda D. Minkova The XIII Inenaonal Confeence Appled Sochasc Models and Daa Analyss (ASMDA-009) Jne 30-Jly 3, 009, Vlns, LITHUANIA ISBN 978-9955-8-463-5 L Sakalaskas, C Skadas and E K Zavadskas (Eds): ASMDA-009 Seleced

More information

Hierarchical Production Planning in Make to Order System Based on Work Load Control Method

Hierarchical Production Planning in Make to Order System Based on Work Load Control Method Unvesal Jounal of Indusal and Busness Managemen 3(): -20, 205 DOI: 0.389/ujbm.205.0300 hp://www.hpub.og Heachcal Poducon Plannng n Make o Ode Sysem Based on Wok Load Conol Mehod Ehsan Faah,*, Maha Khodadad

More information

Chapter 3: Vectors and Two-Dimensional Motion

Chapter 3: Vectors and Two-Dimensional Motion Chape 3: Vecos and Two-Dmensonal Moon Vecos: magnude and decon Negae o a eco: eese s decon Mulplng o ddng a eco b a scala Vecos n he same decon (eaed lke numbes) Geneal Veco Addon: Tangle mehod o addon

More information

Chapter Finite Difference Method for Ordinary Differential Equations

Chapter Finite Difference Method for Ordinary Differential Equations Chape 8.7 Fne Dffeence Mehod fo Odnay Dffeenal Eqaons Afe eadng hs chape, yo shold be able o. Undesand wha he fne dffeence mehod s and how o se o solve poblems. Wha s he fne dffeence mehod? The fne dffeence

More information

CptS 570 Machine Learning School of EECS Washington State University. CptS Machine Learning 1

CptS 570 Machine Learning School of EECS Washington State University. CptS Machine Learning 1 ps 57 Machne Leann School of EES Washnon Sae Unves ps 57 - Machne Leann Assume nsances of classes ae lneal sepaable Esmae paamees of lnea dscmnan If ( - -) > hen + Else - ps 57 - Machne Leann lassfcaon

More information

) from i = 0, instead of i = 1, we have =

) from i = 0, instead of i = 1, we have = Chape 3: Adjusmen Coss n he abou Make I Movaonal Quesons and Execses: Execse 3 (p 6): Illusae he devaon of equaon (35) of he exbook Soluon: The neempoal magnal poduc of labou s epesened by (3) = = E λ

More information

Field due to a collection of N discrete point charges: r is in the direction from

Field due to a collection of N discrete point charges: r is in the direction from Physcs 46 Fomula Shee Exam Coulomb s Law qq Felec = k ˆ (Fo example, f F s he elecc foce ha q exes on q, hen ˆ s a un veco n he decon fom q o q.) Elecc Feld elaed o he elecc foce by: Felec = qe (elecc

More information

A BRANCH-AND-PRICE METHOD FOR THE VEHICLE ROUTING PROBLEM WITH CROSS-DOCKING AND TIME WINDOWS

A BRANCH-AND-PRICE METHOD FOR THE VEHICLE ROUTING PROBLEM WITH CROSS-DOCKING AND TIME WINDOWS A BRANCH-AN-PRCE EHO FOR HE VEHCLE ROUNG PROBLE WH CROSS-OCKNG AN E WNOWS Rodolfo ondo ABSRAC: One mpoan faco n supply chan managemen s o effcenly conol he supply chan flows. ue o s mpoance, many companes

More information

Outline. GW approximation. Electrons in solids. The Green Function. Total energy---well solved Single particle excitation---under developing

Outline. GW approximation. Electrons in solids. The Green Function. Total energy---well solved Single particle excitation---under developing Peenaon fo Theoecal Condened Mae Phyc n TU Beln Geen-Funcon and GW appoxmaon Xnzheng L Theoy Depamen FHI May.8h 2005 Elecon n old Oulne Toal enegy---well olved Sngle pacle excaon---unde developng The Geen

More information

Graduate Macroeconomics 2 Problem set 5. - Solutions

Graduate Macroeconomics 2 Problem set 5. - Solutions Graduae Macroeconomcs 2 Problem se. - Soluons Queson 1 To answer hs queson we need he frms frs order condons and he equaon ha deermnes he number of frms n equlbrum. The frms frs order condons are: F K

More information

CHAPTER 3 DETECTION TECHNIQUES FOR MIMO SYSTEMS

CHAPTER 3 DETECTION TECHNIQUES FOR MIMO SYSTEMS 4 CAPTER 3 DETECTION TECNIQUES FOR MIMO SYSTEMS 3. INTRODUCTION The man challenge n he paccal ealzaon of MIMO weless sysems les n he effcen mplemenaon of he deeco whch needs o sepaae he spaally mulplexed

More information

Solution in semi infinite diffusion couples (error function analysis)

Solution in semi infinite diffusion couples (error function analysis) Soluon n sem nfne dffuson couples (error funcon analyss) Le us consder now he sem nfne dffuson couple of wo blocks wh concenraon of and I means ha, n a A- bnary sysem, s bondng beween wo blocks made of

More information

ajanuary't I11 F or,'.

ajanuary't I11 F or,'. ',f,". ; q - c. ^. L.+T,..LJ.\ ; - ~,.,.,.,,,E k }."...,'s Y l.+ : '. " = /.. :4.,Y., _.,,. "-.. - '// ' 7< s k," ;< - " fn 07 265.-.-,... - ma/ \/ e 3 p~~f v-acecu ean d a e.eng nee ng sn ~yoo y namcs

More information

A hybrid method to find cumulative distribution function of completion time of GERT networks

A hybrid method to find cumulative distribution function of completion time of GERT networks Jounal of Indusal Engneeng Inenaonal Sepembe 2005, Vol., No., - 9 Islamc Azad Uvesy, Tehan Souh Banch A hybd mehod o fnd cumulave dsbuon funcon of compleon me of GERT newos S. S. Hashemn * Depamen of Indusal

More information

When to Treat Prostate Cancer Patients Based on their PSA Dynamics

When to Treat Prostate Cancer Patients Based on their PSA Dynamics When o Tea Posae Cance Paens Based on he PSA Dynamcs CLARA day on opeaons eseach n cance eamen & opeaons managemen Novembe 7 00 Mael S. Lave PhD Man L. Pueman PhD Sco Tyldesley M.D. Wllam J. Mos M.D CIHR

More information

( ) () we define the interaction representation by the unitary transformation () = ()

( ) () we define the interaction representation by the unitary transformation () = () Hgher Order Perurbaon Theory Mchael Fowler 3/7/6 The neracon Represenaon Recall ha n he frs par of hs course sequence, we dscussed he chrödnger and Hesenberg represenaons of quanum mechancs here n he chrödnger

More information

Real-coded Quantum Evolutionary Algorithm for Global Numerical Optimization with Continuous Variables

Real-coded Quantum Evolutionary Algorithm for Global Numerical Optimization with Continuous Variables Chnese Jounal of Eleconcs Vol.20, No.3, July 2011 Real-coded Quanum Evoluonay Algohm fo Global Numecal Opmzaon wh Connuous Vaables GAO Hu 1 and ZHANG Ru 2 (1.School of Taffc and Tanspoaon, Souhwes Jaoong

More information

The Unique Solution of Stochastic Differential Equations. Dietrich Ryter. Midartweg 3 CH-4500 Solothurn Switzerland

The Unique Solution of Stochastic Differential Equations. Dietrich Ryter. Midartweg 3 CH-4500 Solothurn Switzerland The Unque Soluon of Sochasc Dffeenal Equaons Dech Rye RyeDM@gawne.ch Mdaweg 3 CH-4500 Solohun Swzeland Phone +4132 621 13 07 Tme evesal n sysems whou an exenal df sngles ou he an-iô negal. Key wods: Sochasc

More information

A Methodology for Detecting the Change of Customer Behavior based on Association Rule Mining

A Methodology for Detecting the Change of Customer Behavior based on Association Rule Mining A Mehodology fo Deecng he Change of Cusome Behavo based on Assocaon Rule Mnng Hee Seo Song, Soung He Km KAIST Gaduae School of Managemen Jae Kyeong Km KyungHee Unvesy Absac Undesandng and adapng o changes

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 5 Lecure 9 Hamlonan Equaons of Moon (Chaper 8) Wha We Dd Las Tme Consruced Hamlonan formalsm H ( q, p, ) = q p L( q, q, ) H p = q H q = p H = L Equvalen o Lagrangan formalsm Smpler, bu

More information

Accelerated Sequen.al Probability Ra.o Test (SPRT) for Ongoing Reliability Tes.ng (ORT)

Accelerated Sequen.al Probability Ra.o Test (SPRT) for Ongoing Reliability Tes.ng (ORT) cceleaed Sequen.al Pobably Ra.o Tes (SPRT) fo Ongong Relably Tes.ng (ORT) Mlena Kasch Rayheon, IDS Copygh 25 Rayheon Company. ll ghs eseved. Cusome Success Is Ou Msson s a egseed adema of Rayheon Company

More information

Variability Aware Network Utility Maximization

Variability Aware Network Utility Maximization aably Awae Newok ly Maxmzaon nay Joseph and Gusavo de ecana Depamen of Eleccal and Compue Engneeng, he nvesy of exas a Ausn axv:378v3 [cssy] 3 Ap 0 Absac Newok ly Maxmzaon NM povdes he key concepual famewok

More information

N 1. Time points are determined by the

N 1. Time points are determined by the upplemena Mehods Geneaon of scan sgnals In hs secon we descbe n deal how scan sgnals fo 3D scannng wee geneaed. can geneaon was done n hee seps: Fs, he dve sgnal fo he peo-focusng elemen was geneaed o

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 5 Lecure 9 Hamlonan Equaons of Moon (Chaper 8) Wha We Dd Las Tme Consruced Hamlonan formalsm Hqp (,,) = qp Lqq (,,) H p = q H q = p H L = Equvalen o Lagrangan formalsm Smpler, bu wce as

More information

SCIENCE CHINA Technological Sciences

SCIENCE CHINA Technological Sciences SIENE HINA Technologcal Scences Acle Apl 4 Vol.57 No.4: 84 8 do:.7/s43-3-5448- The andom walkng mehod fo he seady lnea convecondffuson equaon wh axsymmec dsc bounday HEN Ka, SONG MengXuan & ZHANG Xng *

More information

Solution of Non-homogeneous bulk arrival Two-node Tandem Queuing Model using Intervention Poisson distribution

Solution of Non-homogeneous bulk arrival Two-node Tandem Queuing Model using Intervention Poisson distribution Volume-03 Issue-09 Sepembe-08 ISSN: 455-3085 (Onlne) RESEARCH REVIEW Inenaonal Jounal of Muldscplnay www.jounals.com [UGC Lsed Jounal] Soluon of Non-homogeneous bulk aval Two-node Tandem Queung Model usng

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 5 Lecure 0 Canoncal Transformaons (Chaper 9) Wha We Dd Las Tme Hamlon s Prncple n he Hamlonan formalsm Dervaon was smple δi δ Addonal end-pon consrans pq H( q, p, ) d 0 δ q ( ) δq ( ) δ

More information

Molecular dynamics modeling of thermal and mechanical properties

Molecular dynamics modeling of thermal and mechanical properties Molecula dynamcs modelng of hemal and mechancal popees Alejando Sachan School of Maeals Engneeng Pudue Unvesy sachan@pudue.edu Maeals a molecula scales Molecula maeals Ceamcs Meals Maeals popees chas Maeals

More information

Simulation of Non-normal Autocorrelated Variables

Simulation of Non-normal Autocorrelated Variables Jounal of Moden Appled Sascal Mehods Volume 5 Issue Acle 5 --005 Smulaon of Non-nomal Auocoelaed Vaables HT Holgesson Jönöpng Inenaonal Busness School Sweden homasholgesson@bshse Follow hs and addonal

More information

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS R&RATA # Vol.) 8, March FURTHER AALYSIS OF COFIDECE ITERVALS FOR LARGE CLIET/SERVER COMPUTER ETWORKS Vyacheslav Abramov School of Mahemacal Scences, Monash Unversy, Buldng 8, Level 4, Clayon Campus, Wellngon

More information

Answers to Tutorial Questions

Answers to Tutorial Questions Inoducoy Mahs ouse Answes.doc Answes o Tuoal Quesons Enjoy wokng hough he examples. If hee ae any moe quesons, please don hesae o conac me. Bes of luck fo he exam and beyond, I hope you won need. Tuoal

More information

Physics 201 Lecture 15

Physics 201 Lecture 15 Phscs 0 Lecue 5 l Goals Lecue 5 v Elo consevaon of oenu n D & D v Inouce oenu an Iulse Coens on oenu Consevaon l oe geneal han consevaon of echancal eneg l oenu Consevaon occus n sses wh no ne eenal foces

More information

An introduction to Support Vector Machine

An introduction to Support Vector Machine An nroducon o Suppor Vecor Machne 報告者 : 黃立德 References: Smon Haykn, "Neural Neworks: a comprehensve foundaon, second edon, 999, Chaper 2,6 Nello Chrsann, John Shawe-Tayer, An Inroducon o Suppor Vecor Machnes,

More information

Variants of Pegasos. December 11, 2009

Variants of Pegasos. December 11, 2009 Inroducon Varans of Pegasos SooWoong Ryu bshboy@sanford.edu December, 009 Youngsoo Cho yc344@sanford.edu Developng a new SVM algorhm s ongong research opc. Among many exng SVM algorhms, we wll focus on

More information

Set of square-integrable function 2 L : function space F

Set of square-integrable function 2 L : function space F Set of squae-ntegable functon L : functon space F Motvaton: In ou pevous dscussons we have seen that fo fee patcles wave equatons (Helmholt o Schödnge) can be expessed n tems of egenvalue equatons. H E,

More information

I-Hsuan Hong Hsi-Mei Hsu Yi-Mu Wu Chun-Shao Yeh

I-Hsuan Hong Hsi-Mei Hsu Yi-Mu Wu Chun-Shao Yeh Poceedngs of he 8 Wne Smulaon Confeence S. J. ason, R. R. Hll, L. önch, O. Rose, T. Jeffeson, J. W. Fowle eds. PRICING DECISION AND LEAD TIE SETTING IN A DUOPOL SEICONDUCTOR INDUSTR I-Hsuan Hong Hs-e Hsu

More information

APPROXIMATIONS FOR AND CONVEXITY OF PROBABILISTICALLY CONSTRAINED PROBLEMS WITH RANDOM RIGHT-HAND SIDES

APPROXIMATIONS FOR AND CONVEXITY OF PROBABILISTICALLY CONSTRAINED PROBLEMS WITH RANDOM RIGHT-HAND SIDES R U C O R R E S E A R C H R E P O R APPROXIMAIONS FOR AND CONVEXIY OF PROBABILISICALLY CONSRAINED PROBLEMS WIH RANDOM RIGH-HAND SIDES M.A. Lejeune a A. PREKOPA b RRR 7-005, JUNE 005 RUCOR Ruges Cene fo

More information

Reflection and Refraction

Reflection and Refraction Chape 1 Reflecon and Refacon We ae now neesed n eplong wha happens when a plane wave avelng n one medum encounes an neface (bounday) wh anohe medum. Undesandng hs phenomenon allows us o undesand hngs lke:

More information

If there are k binding constraints at x then re-label these constraints so that they are the first k constraints.

If there are k binding constraints at x then re-label these constraints so that they are the first k constraints. Mathematcal Foundatons -1- Constaned Optmzaton Constaned Optmzaton Ma{ f ( ) X} whee X {, h ( ), 1,, m} Necessay condtons fo to be a soluton to ths mamzaton poblem Mathematcally, f ag Ma{ f ( ) X}, then

More information

Lecture 11 SVM cont

Lecture 11 SVM cont Lecure SVM con. 0 008 Wha we have done so far We have esalshed ha we wan o fnd a lnear decson oundary whose margn s he larges We know how o measure he margn of a lnear decson oundary Tha s: he mnmum geomerc

More information

Efficient Bayesian Network Learning for System Optimization in Reliability Engineering

Efficient Bayesian Network Learning for System Optimization in Reliability Engineering Qualy Technology & Quanave Managemen Vol. 9, No., pp. 97-, 202 QTQM ICAQM 202 Effcen Bayesan Newok Leanng fo Sysem Opmzaon n Relably Engneeng A. Gube and I. Ben-Gal Depamen of Indusal Engneeng, Faculy

More information

A VISCOPLASTIC MODEL OF ASYMMETRICAL COLD ROLLING

A VISCOPLASTIC MODEL OF ASYMMETRICAL COLD ROLLING SISOM 4, BUCHAEST, - May A VISCOPLASTIC MODEL OF ASYMMETICAL COLD OLLING odca IOAN Spu Hae Unvesy Buchaes, odcaoan7@homal.com Absac: In hs pape s gven a soluon of asymmecal sp ollng poblem usng a Bngham

More information

Inventory Policy Implications of On-Line Customer Purchase Behavior

Inventory Policy Implications of On-Line Customer Purchase Behavior Invenoy Polcy Implcaons of On-Lne Cusome Puchase Behavo Hay Goenevel Smon School of Busness Unvesy of Rochese Rochese, NY 4627 goenevel@smon.ochese.edu Panab Majumde Fuqua School of Busness Duke Unvesy

More information

Monetary policy and models

Monetary policy and models Moneay polcy and odels Kes Næss and Kes Haae Moka Noges Bank Moneay Polcy Unvesy of Copenhagen, 8 May 8 Consue pces and oney supply Annual pecenage gowh. -yea ovng aveage Gowh n oney supply Inflaon - 9

More information

ScienceDirect. Behavior of Integral Curves of the Quasilinear Second Order Differential Equations. Alma Omerspahic *

ScienceDirect. Behavior of Integral Curves of the Quasilinear Second Order Differential Equations. Alma Omerspahic * Avalable onlne a wwwscencedeccom ScenceDec oceda Engneeng 69 4 85 86 4h DAAAM Inenaonal Smposum on Inellgen Manufacung and Auomaon Behavo of Inegal Cuves of he uaslnea Second Ode Dffeenal Equaons Alma

More information

MCTDH Approach to Strong Field Dynamics

MCTDH Approach to Strong Field Dynamics MCTDH ppoach o Song Feld Dynamcs Suen Sukasyan Thomas Babec and Msha Ivanov Unvesy o Oawa Canada Impeal College ondon UK KITP Sana Babaa. May 8 009 Movaon Song eld dynamcs Role o elecon coelaon Tunnel

More information

The balanced budget multiplier and labour intensity in home production

The balanced budget multiplier and labour intensity in home production Inenaonal Jounal of Economc Behavo and Oganzaon 205; 3(2-): 23-30 Publshed onlne Febuay 26, 205 (hp://www.scencepublshnggoup.com/j/jebo) do: 0.648/j.jebo.s.20503020.5 ISSN: 2328-7608 (Pn); ISSN: 2328-766

More information

EEL 6266 Power System Operation and Control. Chapter 5 Unit Commitment

EEL 6266 Power System Operation and Control. Chapter 5 Unit Commitment EEL 6266 Power Sysem Operaon and Conrol Chaper 5 Un Commmen Dynamc programmng chef advanage over enumeraon schemes s he reducon n he dmensonaly of he problem n a src prory order scheme, here are only N

More information

Modeling Background from Compressed Video

Modeling Background from Compressed Video Modelng acgound fom Compessed Vdeo Weqang Wang Daong Chen Wen Gao Je Yang School of Compue Scence Canege Mellon Unvesy Psbugh 53 USA Insue of Compung echnology &Gaduae School Chnese Academy of Scences

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Lnear Response Theory: The connecon beween QFT and expermens 3.1. Basc conceps and deas Q: ow do we measure he conducvy of a meal? A: we frs nroduce a weak elecrc feld E, and hen measure

More information

Tight results for Next Fit and Worst Fit with resource augmentation

Tight results for Next Fit and Worst Fit with resource augmentation Tgh resuls for Nex F and Wors F wh resource augmenaon Joan Boyar Leah Epsen Asaf Levn Asrac I s well known ha he wo smple algorhms for he classc n packng prolem, NF and WF oh have an approxmaon rao of

More information

Suppose we have observed values t 1, t 2, t n of a random variable T.

Suppose we have observed values t 1, t 2, t n of a random variable T. Sppose we have obseved vales, 2, of a adom vaable T. The dsbo of T s ow o belog o a cea ype (e.g., expoeal, omal, ec.) b he veco θ ( θ, θ2, θp ) of ow paamees assocaed wh s ow (whee p s he mbe of ow paamees).

More information

An Exact Resource Allocation Model with Hard and Soft Resource Constraints

An Exact Resource Allocation Model with Hard and Soft Resource Constraints An Exac Resouce Allocaon Model wh Had and Sof Resouce Consans Sxeenh Annual Confeence of POMS, Chcago, IL, Apl 29 - May 2, 2005. Feenc Kuzslcz (kuzslc@kk.pe.hu) nvesy of Pécs, Depamen of Busness Infomacs,

More information

Reactive Methods to Solve the Berth AllocationProblem with Stochastic Arrival and Handling Times

Reactive Methods to Solve the Berth AllocationProblem with Stochastic Arrival and Handling Times Reacve Mehods o Solve he Berh AllocaonProblem wh Sochasc Arrval and Handlng Tmes Nsh Umang* Mchel Berlare* * TRANSP-OR, Ecole Polyechnque Fédérale de Lausanne Frs Workshop on Large Scale Opmzaon November

More information

Physics 2A Chapter 11 - Universal Gravitation Fall 2017

Physics 2A Chapter 11 - Universal Gravitation Fall 2017 Physcs A Chapte - Unvesal Gavtaton Fall 07 hese notes ae ve pages. A quck summay: he text boxes n the notes contan the esults that wll compse the toolbox o Chapte. hee ae thee sectons: the law o gavtaton,

More information

, on the power of the transmitter P t fed to it, and on the distance R between the antenna and the observation point as. r r t

, on the power of the transmitter P t fed to it, and on the distance R between the antenna and the observation point as. r r t Lecue 6: Fiis Tansmission Equaion and Rada Range Equaion (Fiis equaion. Maximum ange of a wieless link. Rada coss secion. Rada equaion. Maximum ange of a ada. 1. Fiis ansmission equaion Fiis ansmission

More information

calculating electromagnetic

calculating electromagnetic Theoeal mehods fo alulang eleomagne felds fom lghnng dshage ajeev Thoapplll oyal Insue of Tehnology KTH Sweden ajeev.thoapplll@ee.kh.se Oulne Despon of he poblem Thee dffeen mehods fo feld alulaons - Dpole

More information

An axisymmetric incompressible lattice BGK model for simulation of the pulsatile ow in a circular pipe

An axisymmetric incompressible lattice BGK model for simulation of the pulsatile ow in a circular pipe INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS In. J. Nume. Meh. Fluds 005; 49:99 116 Publshed onlne 3 June 005 n Wley IneScence www.nescence.wley.com). DOI: 10.100/d.997 An axsymmec ncompessble

More information

Chapter Fifiteen. Surfaces Revisited

Chapter Fifiteen. Surfaces Revisited Chapte Ffteen ufaces Revsted 15.1 Vecto Descpton of ufaces We look now at the vey specal case of functons : D R 3, whee D R s a nce subset of the plane. We suppose s a nce functon. As the pont ( s, t)

More information

Heuristic Unit Commitment of Wind Farms Integrated in Power System Consideration with Demand Response

Heuristic Unit Commitment of Wind Farms Integrated in Power System Consideration with Demand Response Avalable onlne a www.scnze.com Scnze Jounal of Engneeng, ol 3, Issue 1, (2017): 40-47 DOI: 10.21634/SJE.3.1.4047 ISSN 2415-105X Heusc Un Commmen of Wnd Fams Inegaed n Powe Sysem Consdeaon wh Demand Response

More information

On One Analytic Method of. Constructing Program Controls

On One Analytic Method of. Constructing Program Controls Appled Mahemacal Scences, Vol. 9, 05, no. 8, 409-407 HIKARI Ld, www.m-hkar.com hp://dx.do.org/0.988/ams.05.54349 On One Analyc Mehod of Consrucng Program Conrols A. N. Kvko, S. V. Chsyakov and Yu. E. Balyna

More information

TRANSIENTS. Lecture 5 ELEC-E8409 High Voltage Engineering

TRANSIENTS. Lecture 5 ELEC-E8409 High Voltage Engineering TRANSIENTS Lece 5 ELECE8409 Hgh Volage Engneeng TRANSIENT VOLTAGES A ansen even s a sholved oscllaon (sgnfcanly fase han opeang feqency) n a sysem cased by a sdden change of volage, cen o load. Tansen

More information

Midterm Exam. Thursday, April hour, 15 minutes

Midterm Exam. Thursday, April hour, 15 minutes Economcs of Grow, ECO560 San Francsco Sae Unvers Mcael Bar Sprng 04 Mderm Exam Tursda, prl 0 our, 5 mnues ame: Insrucons. Ts s closed boo, closed noes exam.. o calculaors of an nd are allowed. 3. Sow all

More information

L4:4. motion from the accelerometer. to recover the simple flutter. Later, we will work out how. readings L4:3

L4:4. motion from the accelerometer. to recover the simple flutter. Later, we will work out how. readings L4:3 elave moon L4:1 To appl Newon's laws we need measuemens made fom a 'fed,' neal efeence fame (unacceleaed, non-oang) n man applcaons, measuemens ae made moe smpl fom movng efeence fames We hen need a wa

More information

A. Thicknesses and Densities

A. Thicknesses and Densities 10 Lab0 The Eath s Shells A. Thcknesses and Denstes Any theoy of the nteo of the Eath must be consstent wth the fact that ts aggegate densty s 5.5 g/cm (ecall we calculated ths densty last tme). In othe

More information

UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 2017 EXAMINATION

UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 2017 EXAMINATION INTERNATIONAL TRADE T. J. KEHOE UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 27 EXAMINATION Please answer wo of he hree quesons. You can consul class noes, workng papers, and arcles whle you are workng on he

More information

CS 188: Artificial Intelligence Fall Probabilistic Models

CS 188: Artificial Intelligence Fall Probabilistic Models CS 188: Aificial Inelligence Fall 2007 Lecue 15: Bayes Nes 10/18/2007 Dan Klein UC Bekeley Pobabilisic Models A pobabilisic model is a join disibuion ove a se of vaiables Given a join disibuion, we can

More information

Chapters 2 Kinematics. Position, Distance, Displacement

Chapters 2 Kinematics. Position, Distance, Displacement Chapers Knemacs Poson, Dsance, Dsplacemen Mechancs: Knemacs and Dynamcs. Knemacs deals wh moon, bu s no concerned wh he cause o moon. Dynamcs deals wh he relaonshp beween orce and moon. The word dsplacemen

More information

Computer Propagation Analysis Tools

Computer Propagation Analysis Tools Compue Popagaion Analysis Tools. Compue Popagaion Analysis Tools Inoducion By now you ae pobably geing he idea ha pedicing eceived signal sengh is a eally impoan as in he design of a wieless communicaion

More information

Numerical Study of Large-area Anti-Resonant Reflecting Optical Waveguide (ARROW) Vertical-Cavity Semiconductor Optical Amplifiers (VCSOAs)

Numerical Study of Large-area Anti-Resonant Reflecting Optical Waveguide (ARROW) Vertical-Cavity Semiconductor Optical Amplifiers (VCSOAs) USOD 005 uecal Sudy of Lage-aea An-Reonan Reflecng Opcal Wavegude (ARROW Vecal-Cavy Seconduco Opcal Aplfe (VCSOA anhu Chen Su Fung Yu School of Eleccal and Eleconc Engneeng Conen Inoducon Vecal Cavy Seconduco

More information

8 Baire Category Theorem and Uniform Boundedness

8 Baire Category Theorem and Uniform Boundedness 8 Bae Categoy Theoem and Unfom Boundedness Pncple 8.1 Bae s Categoy Theoem Valdty of many esults n analyss depends on the completeness popety. Ths popety addesses the nadequacy of the system of atonal

More information

THERMODYNAMICS 1. The First Law and Other Basic Concepts (part 2)

THERMODYNAMICS 1. The First Law and Other Basic Concepts (part 2) Company LOGO THERMODYNAMICS The Frs Law and Oher Basc Conceps (par ) Deparmen of Chemcal Engneerng, Semarang Sae Unversy Dhon Harano S.T., M.T., M.Sc. Have you ever cooked? Equlbrum Equlbrum (con.) Equlbrum

More information

Garrettsville, Ohio Public Water System Lead (Pb) Components

Garrettsville, Ohio Public Water System Lead (Pb) Components Gaesvlle, Oho Publc Wae Sysem Lead (Pb) Componens PWSID# OH6701412 Pepaed Febuay 2017 To comply wh Secon 6109.21 of he Oho Revsed Code, enaced n Sepembe 2016, he Vllage of Gaesvlle n Poage Couny, Oho has

More information

WORK POWER AND ENERGY Consevaive foce a) A foce is said o be consevaive if he wok done by i is independen of pah followed by he body b) Wok done by a consevaive foce fo a closed pah is zeo c) Wok done

More information

Lecture 17: Kinetics of Phase Growth in a Two-component System:

Lecture 17: Kinetics of Phase Growth in a Two-component System: Lecue 17: Kineics of Phase Gowh in a Two-componen Sysem: descipion of diffusion flux acoss he α/ ineface Today s opics Majo asks of oday s Lecue: how o deive he diffusion flux of aoms. Once an incipien

More information

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!) i+1,q - [(! ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL The frs hng o es n wo-way ANOVA: Is here neracon? "No neracon" means: The man effecs model would f. Ths n urn means: In he neracon plo (wh A on he horzonal

More information

Lecture 18: The Laplace Transform (See Sections and 14.7 in Boas)

Lecture 18: The Laplace Transform (See Sections and 14.7 in Boas) Lecure 8: The Lalace Transform (See Secons 88- and 47 n Boas) Recall ha our bg-cure goal s he analyss of he dfferenal equaon, ax bx cx F, where we emloy varous exansons for he drvng funcon F deendng on

More information

Sections 3.1 and 3.4 Exponential Functions (Growth and Decay)

Sections 3.1 and 3.4 Exponential Functions (Growth and Decay) Secions 3.1 and 3.4 Eponenial Funcions (Gowh and Decay) Chape 3. Secions 1 and 4 Page 1 of 5 Wha Would You Rahe Have... $1million, o double you money evey day fo 31 days saing wih 1cen? Day Cens Day Cens

More information

CS 268: Packet Scheduling

CS 268: Packet Scheduling Pace Schedulng Decde when and wha pace o send on oupu ln - Usually mplemened a oupu nerface CS 68: Pace Schedulng flow Ion Soca March 9, 004 Classfer flow flow n Buffer managemen Scheduler soca@cs.bereley.edu

More information

Contact, information, consultations

Contact, information, consultations ontact, nfomaton, consultatons hemsty A Bldg; oom 07 phone: 058-347-769 cellula: 664 66 97 E-mal: wojtek_c@pg.gda.pl Offce hous: Fday, 9-0 a.m. A quote of the week (o camel of the week): hee s no expedence

More information

Dynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005

Dynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005 Dynamc Team Decson Theory EECS 558 Proec Shruvandana Sharma and Davd Shuman December 0, 005 Oulne Inroducon o Team Decson Theory Decomposon of he Dynamc Team Decson Problem Equvalence of Sac and Dynamc

More information

Lecture-V Stochastic Processes and the Basic Term-Structure Equation 1 Stochastic Processes Any variable whose value changes over time in an uncertain

Lecture-V Stochastic Processes and the Basic Term-Structure Equation 1 Stochastic Processes Any variable whose value changes over time in an uncertain Lecue-V Sochasic Pocesses and he Basic Tem-Sucue Equaion 1 Sochasic Pocesses Any vaiable whose value changes ove ime in an unceain way is called a Sochasic Pocess. Sochasic Pocesses can be classied as

More information

Course Outline. 1. MATLAB tutorial 2. Motion of systems that can be idealized as particles

Course Outline. 1. MATLAB tutorial 2. Motion of systems that can be idealized as particles Couse Oulne. MATLAB uoal. Moon of syses ha can be dealzed as pacles Descpon of oon, coodnae syses; Newon s laws; Calculang foces equed o nduce pescbed oon; Deng and solng equaons of oon 3. Conseaon laws

More information

Energy in Closed Systems

Energy in Closed Systems Enegy n Closed Systems Anamta Palt palt.anamta@gmal.com Abstact The wtng ndcates a beakdown of the classcal laws. We consde consevaton of enegy wth a many body system n elaton to the nvese squae law and

More information

TSS = SST + SSE An orthogonal partition of the total SS

TSS = SST + SSE An orthogonal partition of the total SS ANOVA: Topc 4. Orhogonal conrass [ST&D p. 183] H 0 : µ 1 = µ =... = µ H 1 : The mean of a leas one reamen group s dfferen To es hs hypohess, a basc ANOVA allocaes he varaon among reamen means (SST) equally

More information

A Tour of Modeling Techniques

A Tour of Modeling Techniques A Tour of Modelng Technques John Hooker Carnege Mellon Unversy EWO Semnar February 8 Slde Oulne Med neger lnear (MILP) modelng Dsuncve modelng Knapsack modelng Consran programmng models Inegraed Models

More information