MODELING OF PLATE HEAT EXCHANGERS WITH GENERALIZED CONFIGURATIONS

Size: px
Start display at page:

Download "MODELING OF PLATE HEAT EXCHANGERS WITH GENERALIZED CONFIGURATIONS"

Transcription

1 XV OGRESSO BRASLERO DE EGEHARA MEÂA 6t BRAZLA OGRESS OF MEHAAL EGEERG MODELG OF LATE HEAT EXHAGERS WTH GEERALZED OFGURATOS Jorge Andrey Wlelms Gut Department of emal Engneerng - Unversty of São aulo Av. rof. Luano Gualberto, trav. 3, 380 / São aulo - S / / BRAZL andrey@lsp.pq.ep.usp.br José Maurío nto Department of emal Engneerng - Unversty of São aulo Av. rof. Luano Gualberto, trav. 3, 380 / São aulo - S / / BRAZL jompnto@usp.br Abstrat. A matematal model s developed for te smulaton of gasketed plate eat exangers, operatng n steady state, wt a general onfguraton. Te man purposes of ts model are to study te nfluene of te onfguraton on te equpment performane and to furter develop a metod for optmzng te exanger onfguraton. Te onfguraton s defned by sx parameters, w are as follows: number of annels, number of passes at ea sde, flud loatons, feed onneton loatons and type of annel-flow. Te matematal model s developed troug an assemblng algortm, sne t s not possble to represent te model expltly as a funton of te sx parameters. Te resultng system s omposed of ordnary dfferental equatons of te boundary value type, w s solved by te fnte dfferene metod, usng te software groms (roess Systems Enterprse, 00). Te man smulaton results are te followng: temperature profles n all annels, termal effetveness, dstrbuton of te overall eat transfer oeffent along te exanger and flud pressure drops. Examples of eat exangers wt approxmately 50 plates yeld algebra models wt tousands of equatons, w are solved wtn mnutes. Moreover, te assumpton of nvarant overall eat transfer oeffent s verfed. Keywords: plate eat exanger, matematal modelng, fnte dfferene metod, eat exanger onfguraton. ntroduton For beng ompat, easy to lean, effent and very flexble, te gasketed plate eat exanger (HE) s wdely employed n te emal, food and parmaeutal proess ndustres. Te HE onssts of a pak of gasketed orrugated metal plates, pressed togeter n a frame (see Fg. ). Te gaskets on te orners of te plates form a seres of parallel flow annels, were te fluds flow alternately and exange eat troug te tn metal plates. Te gasket desgn and te losed ports of te plates determne te flud flow dstrbuton, w an be parallel, seres or any of ter varous possble ombnatons. Te number of plates, flow dstrbuton, type of gaskets and te flud feed loatons araterze te exanger onfguraton. Fgure. Te plate eat exanger assemblage, examples of flow dstrbutons and plate man dmensons. Te smplfed termal modelng of a HE n steady state yelds a lnear system of frst order ordnary dfferental equatons, omprsng te energy balane for ea annel and te requred boundary ondtons. Te man assumptons are as follows: plug-flow nsde te annels, onstant overall eat transfer oeffent along te exanger, unform dstrbuton of flow n te annels, no eat loss and no eat exange n te flow dreton. Ts bas model was presented by MKllop & Dunkley (960) for 3 dfferent onfguratons and also by Masubu & to (977), were te

2 roeedngs of OBEM 00, Energy and Termal Systems, Vol.4, 377 dynam responses of some usual onfguratons were studed. n bot works, te Runge-Kutta-Gll ntegraton metod was used to solve te system of equatons. Te ntegraton s non-trval beause te boundary ondtons are defned at dfferent extremes of te annel. Approxmate solutons were developed by Settar & Venart (97) n polynomal form, and by Zalesk & Klepaka (99) n exponental form. Bot metods lead to good approxmatons of te exat soluton, but tey may not be relable wen tere s a large dfferene between flud eat apates. Kandlkar & Sa (989b) developed a metod to alulate an approxmated termal effetveness for large exangers, were te effets of te end plates and of te anges of passes an be negleted. n ts ase, te exanger s dvded nto a group of smpler exangers tat are nteronneted, wt known effetveness. Te analytal soluton of te system of equatons n matrx form was studed by Zalesk & Jarzebsk (973) and Zalesk (984) for exangers wt seres and parallel arrangements. Ts soluton metod may lead to numeral problems on te alulaton of egenvalues and egenvetors, and t s not reommended for large szed exangers. Kandlkar & Sa (989a) and Georgads et al. (998) used te fnte dfferene metod for te smulaton of HEs. Kandlkar & Sa (989a) smulated and ompared several onfguratons. t was verfed tat ger effetveness s aeved wen te exanger s symmetral, wt te same numbers of passes for bot streams, beause te annels tat are next to te anges of passes as well as te end annels ave a lower effetveness. However, wen te fluds ave very dfferent flow rates or eat apates, a non-symmetral onfguraton must be used. n su ases, tere s no rgorous desgn metod to selet te best onfguraton, w s made by omparson among te usual onfguratons from termal effetveness and pressure drop vewponts. Georgads et al. (998) presented a detaled modelng of a HE used for mlk pasteurzaton tat ouples te dynam termal model wt te proten-foulng model. Tree dfferent onfguratons were ompared and te reduton of te overall eat transfer oeffent, aused by te proten adeson on te plates, was studed. Te model was solved wt te fnte dfferene metod, mplemented n te software groms (roess System Enterprse, 00). To te autors knowledge tere s no rgorous desgn metod for HEs n te open lterature, as tere are for te sell-and-tube exangers (Taborek, 983). Sa & Foke (988) ave presented a detaled step-by-step desgn proedure for ratng and szng a HE, w s owever restrted to parallel flow arrangements. n all of tose works, te overall eat transfer oeffent was onsdered nvarable along te exanger beause ts assumpton brngs a great smplfaton for te soluton of te system of equatons, lnearzng te dfferental equatons of te annel flud temperatures. However, tere may be a onsderable varaton of te overall oeffent for some ases, su as te seres flow arrangement wt equal flow rates of te fluds (Buonopane et al, 963). Te am of ts work s to present a HE modelng framework tat s sutable for any onfguraton. Te purpose of developng su model s to study te nfluene of te onfguraton on te exanger performane and to furter develop an optmzaton metod for onfguraton seleton. Te frst step for te modelng onerns te parameterzaton of te dfferent onfguratons, followng te work of gnott & Tamborenea (988). Te varaton of te overall eat transfer oeffent along te exanger s also studed n ts work, wt respet to te assumpton of a onstant value. omenlature a gener model parameter A effetve eat transfer area (m ) b annel average tkness (m) p flud spef eat at onstant pressure (J/kg.K) D e equvalent dameter of annel (m) D port dameter (m) E exanger termal effetveness (%) f Fannng frton fator g gravtatonal aeleraton (g 9,8 m/s ) G annel mass veloty (kg/m.s) G port mass veloty (kg/m.s) onvetve eat transfer oeffent (W/m.K) k flud termal ondutvty (W/m.K) k plate termal ondutvty (W/m.K) L effetve plate lengt, measured between ports (m) M dagonal matrx defned wt Eq. (0) number of annels per pass number of annels u usselt number number of passes r randtl number R flud foulng fator (m.k/w) Re Reynolds number s annel flow dreton parameter (s + or ) T temperature (K) U overall eat transfer oeffent (W/m.K) w effetve plate wdt, measured between gaskets (m) W flud mass flow rate (kg/s) x oordnate, tangental to annel flud flow (m) Y f bnary parameter for type of annel-flow Y bnary parameter for ot flud loaton α dmensonless varable defned wt Eq. (9) β evron orrugaton nlnaton angle (degrees) flud pressure drop (a) ρ flud densty (kg/m 3 ) ε maxmum devaton for alulated termal effetveness ε tkness of metal plate (m) φ parameter for feed onneton relatve loaton Φ enlargement fator of plate η dmensonless oordnate, tangental to annel flud flow µ flud vsosty (a.s) θ dmensonless flud temperature Subsrpts n max mn out w old flud ot flud gener element flud nlet maxmum value mnmum value flud outlet at te plate wall Supersrpts sde of exanger sde of exanger

3 roeedngs of OBEM 00, Energy and Termal Systems, Vol.4, 378. onfguraton araterzaton Te onfguraton of a HE s defned by te nformaton tat allows te detalng of te equpment assemblage, nludng te onnetons on te fxed and moveable overs, te losed and open ports n ea plate and te type and poston of ea gasket. To araterze su onfguratons, sx dstnt parameters are used:,,, φ, Y and Y f. Tese parameters are defned wt detals n Tab. (). Table. araterzaton of te sx onfguraton parameters. : umber of annels arameter desrpton Te spae omprsed between two plates s a annel, and te HE an be represented by a row of annels, numbered from to. Te odd-numbered annels belong to sde, and te even-numbered ones belong to sde (as an analogy to te tube and sell sdes n a sell-and-tube exanger). and are te numbers of annels n ea sde. f s even, bot sdes ave te same number of annels, oterwse sde as one more annel. Allowable values:, 3, 4, 5 and : umber of passes at sdes and A pass s a set of annels were te stream s splt and dstrbuted. For a regular onfguraton, ea sde of te HE s splt nto passes wt te same number of annels per pass ( and ). asses wt dfferent numbers of annels are unusual (Kakaç & Lu, 998). Te relatonsp between, and s gven n Tab. (). Allowable values: from te fatorzaton of and respetvely φ : Feed onneton relatve loaton Te feed onneton of sde s arbtrarly set n annel at η0. Te relatve poston of te feed of sde s gven by te parameter φ, as sown n te dagram (gnott & Tamborenea, 988). Te dmensonless lengt η s not assoated wt te top and bottom of te HE, neter annel s assoated wt te fxed over. Te onfguraton an be freely rotated or mrrored. Allowable values:,, 3 and 4 Y : Hot flud loaton Ts bnary parameter assgns te fluds to te exanger sdes: llustratve fgure - f Y : te ot flud s at sde, and te old flud at sde - f Y 0 : te old flud s at sde, and te ot flud at sde Y f : Type of flow n annels Y f s a bnary parameter tat defnes te type of flow nsde te annels. As sown n te dagram, te flow an be stragt or rossed, dependng on te gasket type. Te rossed flow avods te formaton of stagnaton areas, but te stragt flow type s easer to assemble. t s not possble to use bot types togeter. - f Y f : te flow s rossed n all annels - f Y f 0 : te flow s stragt n all annels Te sx parameters an represent any regular onfguraton. For a fxed number of annels, te fve remanng parameters ave a known set of allowable values. Te ombnaton gves a fnte number of possble regular onfguratons for ea value of, as presented n Fg. (). Te dsperse pattern s due to te varaton of te number of nteger fators of and for ea value of. A onfguraton example s sown n Fg. (3) for llustraton. t represents an egt-plate HE, were te ot flud n sde (Y ) makes passes ( ) and te old flud on sde makes 3 passes ( 3). n ts example, te nlet of

4 roeedngs of OBEM 00, Energy and Termal Systems, Vol.4, 379 sde s loated next to te nlet of sde (φ) and te type of annel-flow s rossed flow (Y f ). Te parameter Y f for type of annel-flow s mostly useful for te exanger pysal onstruton and t may not be neessary for te smulaton sne ts nfluene over te onvetve oeffents and frton fators s usually unknown. Table. Relatonsp between numbers of annels and passes. Man equatons +.. f s even:.. + f s odd: +.. 7,000 umber of possble regular onfguratons 6,000 5,000 4,000 3,000,000, umber of annels - Fgure. umber of possble regular onfguratons as a funton of te number of annels. Fgure 3. Example of onfguraton for a HE wt egt plates. 3. Equvalent onfguratons For a gven value of number of annels and a fxed type of flow, te exstene of equvalent onfguratons (tat ave te same termal effetveness and pressure drops) s possble. Te dentfaton of te equvalent onfguratons s mportant to avod unneessary smulatons. Te equvalene ours due to 3 remarks, as follows: A) Aordng to te property of flow reversblty (gnott & Tamborenea, 988), te nverson of te flud flow dreton n bot sdes does not alter te effetveness of te HE. A) Wen tere s a sngle pass n a sde, te flow dreton s te same n all annels, regardless f te feed s loated at te frst annel or at te last one. A3) Smply nvertng te dreton of η, or numberng te annels n reverse order, may yeld to a new set of onfguraton parameters. A metodology to detet equvalent onfguratons s presented n Tab. (3). For ea set,, and Y f tere are groups of te parameter φ tat result n equvalent onfguratons. n te ase of an even-numbered, tere may be equvaleny between Y 0 and Y, beause sdes and ave te same number of annels and terefore an support te same passes. An example of equvaleny between 4 dfferent onfguratons s sown n Fg. (4).

5 Table 3. dentfaton of equvalent onfguratons for gven values of and Y f. roeedngs of OBEM 00, Energy and Termal Systems, Vol.4, 380 (, ) Groups of equvalent values of φ Reduton n te number of smulatons (, ) ; (, odd) ; (odd, ) {, 3} ; {, 4} 50 % odd (, even) ; (even, ) {,, 3, 4} 75 % (odd, odd) ; (even, even) {} ; {} ; {3} ; {4} 0 % (odd, even) ; (even, odd) {, } ; {3, 4} 50 % (, ) ; (, odd) ; (odd, ) {, 3,, 3} ; {, 4,, 4} 75 % (, even) {, 4,, 4} ; {, 3,, 3} 75 % even (even, ) {, 3,, 3} ; {, 4,, 4} 75 % (odd, odd) ; (even, even) {, } ; {, } ; {3, 3} ; {4, 4} 50 % (odd, even) ; (even, odd) {, } ; {, } ; {3, 3} ; {4, 4} 50 % ote: wen s even, denotes Y and denotes Y 0. Wen s odd, all equvalent onfguratons ave te same value for Y. Fgure 4. Example of four equvalent onfguratons wt 6 and Y f. 4. Modelng of te HE Te followng assumptons were made n order to derve te matematal model: B) Steady state operaton. B) o eat loss to surroundngs. B3) o eat exange on te dreton of flow. B4) lug-flow nsde te annels. B5) Unform dstrbuton of flow troug te annels of a pass. B6) erfet mxture of flud n te end of a pass. B7) Fluds wt ewtonan beavor. B8) o pase anges. Te flud nsde a annel exanges eat wt te negbor annels troug te tn metal plates, as sown n Fg. (5). Te effetve eat exange area s A Φ.w. L, were Φ s te area enlargement fator, w aounts for te orrugaton wrnkles. Te lengt of te pat x and te flud temperature T(x) an be onverted nto dmensonless form as sown n Eqs () and (). x η( x ) 0 η () L T T,n θ (T ) 0 θ () T T,n,n

6 roeedngs of OBEM 00, Energy and Termal Systems, Vol.4, 38 Applyng te energy balane to te ontrol volume sown n Fg. (5) (MKllop & Dunkley, 960) t s possble to derve te dfferental equatons for te annel temperature (Eqs 3, 4 and 5). Te overall eat transfer oeffent U between annels and + s gven by Eq. (6) as a funton of te flud onvetve eat transfer oeffent, te plate termal ondutvty k and te foulng fator R for ot and old streams. Te mass flow rate nsde annel s alulated by Eq. (7) aordng to assumpton B5, were sde() refers to te sde tat ontans annel. t s known tat ts assumpton may not old for passes wt a large number of annels (Bassouny & Martn, 984), but tere s no rgorous modelng to ts beavor. Te flow rates n ea sde, W and W, are assoated to te ot and old flud flow rates troug te parameter Y. Te dreton of te flow n annel s gven by te varable s. f te flow follows te dreton of η, ten s +, oterwse s. For te frst annel: dθ s.a U. θ dη W.p [ ( θ )] (3) For annel : ( < < ) dθ s.a [ U.( θ θ ) + U.( θ+ θ )] (4) dη W.p For te last annel: dθ s.a dη W.p [ U.( θ θ )] (5) U + + ε + k + R + R, ( ) (6) Fgure 5. ontrol volume for dervaton of energy balane nsde an upward flow annel. sde( ) W W, sde( ) sde() {, } (7) Te overall eat transfer oeffent s a funton of te flud temperature n te annels (Eq. 6), w depends on η. Sne ts makes U also a funton of η, te soluton of te system of dfferental Eqs (3), (4) and (5) s not smple. However, f te flud pysal propertes are assumed onstant, Eq. (4) an be smplfed to Eq. (8), were te oeffents α and α are gven by Eq. (9) as a funton of te onstant overall eat transfer oeffent U. onsequently, te system of equatons s redued to a lnear system of ordnary dfferental equatons, w s represented n matrx form n Eq. (0), were M s a trdagonal matrx and θ s te vetor of annel dmensonless temperatures θ (η). Te analytal soluton of Eq. (0) was studed by Zalesk & Jarzebsk (973). dθ dη s. α sde( ).( θ. θ + θ ) +, sde() {, } (8) sde( ) sde( ) A.U. α, sde( ) sde( ) sde() {, } (9) W.p dθ M. θ dη (0) Te neessary boundary ondtons for te smplfed model, defned by Eq. (8), and for te rgorous model, defned by Eqs (3) troug (6), are presented n Tab. (4). Every annel requres a boundary ondton equaton for ts nlet temperature. Te nlet poston of annel s gven by te value of s. f s + te nlet s loated at poston η 0, oterwse te nlet s loated at η. Te number and struture of te requred boundary equatons are funtons of te onfguraton parameters. For te performane evaluaton of te eat exanger, te termal effetveness and te pressure drop alulatons are needed. One te system of equatons s solved, te termal effetveness E an be alulated by Eq. () usng te average value of te flud spef eats, were W and W are related to W sde() by te parameter Y. W.p E mn W.( θ ) (.p,w.p ) mn( W.p,W.p ),out W.p. θ,out ()

7 Table 4. Types of termal boundary ondtons for te HE annels. roeedngs of OBEM 00, Energy and Termal Systems, Vol.4, 38 Boundary ondton Flud entrane: te temperature at te entrane of te frst pass s te same as te stream nlet temperature. anges of pass: tere s a perfet mxture of te flud leavng te annels of a pass, before enterng te next one. Flud ext: te stream outlet temperature results from a perfet mxture of te flud leavng te last pass. θ η) Equaton Form θ ( θ, frst pass flud,n ( η ). θ j ( η ) j prevous pass θ flud,out. θ j ( η ) j last pass, urrent pass Te flud pressure drops at sdes and, and, an be alulated by Eq. () tat reles on assumptons B7 and B8 (Sa & Foke, 988; Kakaç & Lu, 998). Te frst term n te rgt-and sde evaluates te frton loss nsde te annels, were G denotes te annel mass veloty (Eq. 3a). Te seond term represents te pressure drop for port flow, were G s te port mass veloty (Eq. 3b). Te last term s te pressure varaton due to an elevaton ange. Sne te gravtatonal aeleraton dreton s not assoated to te vertal dmenson η and tere s no nformaton on te pump loaton, ts term s always onsdered. Terefore te pressure drop may be overestmated..( L + D ). f..g G +, ρ.g. +.D e. ρ ρ ( L D ) for sdes and () W G,.b.w G 4W. for sdes and (3) π.d Te neessary onsttutve equatons for te alulaton of onvetve oeffents and frton fator, aordng to assumpton B7, are presented n Tab. (5), as well as te neessary dmensonless numbers and te annel equvalent dameter D e. Usual values for te empral parameters a to a 6 of te onsttutve equatons an be found on Saunders (988) and Sa & Foke (988). Table 5. Equatons for onvetve oeffents and frton fator alulaton, for sdes and. onsttutve Equatons Dmensonless umbers Equvalent Dameter u a.re a.r a3 a f a + Re. 5 4 a6 µ µ w 0, 7.De u, k. r k µ D e Re G. µ 4.b.w.b D e.( b + w. Φ) Φ 5. Smulaton of te HE Te matematal modelng of a HE, for te alulaton of ts termal effetveness and flud pressure drops, was presented n te prevous seton. However, t s not possble to derve a model tat s expltly a funton of te onfguraton parameters, espeally beause of te bnary varable s and of te requred boundary ondtons (te equatons n Tab. 4). To overome ts lmtaton, te modelng was developed n te form of an assemblng algortm. Gven te onfguraton parameters, ts algortm gudes te onstruton of te omplete matematal model of te HE and ts smulaton at steady state operaton. n ts work, te soluton of te system of equatons s arred out by te software groms (roess Systems Enterprse, 00), usng te seond order entered fnte dfferene metod. n ts metod, all varables dependng on te plate lengt η are dsretzed wt respet to ts dmenson. Several semes were tested and, wt 0 dsretzaton ntervals wtn η [0,], an exellent approxmaton of te analytal result was aeved. Te smplfed model tat reles on te assumpton of onstant overall eat transfer oeffent and te rgorous model developed by relaxng ts assumpton were tested and ompared. Te man struture of te assemblng algortm for te smplfed model s presented n Fg. (6). Te algortm as 4 steps as follows:

8 roeedngs of OBEM 00, Energy and Termal Systems, Vol.4, 383 Fgure 6. Man struture of te model assemblng algortm. ) Te requred data for te HE (L, w, b, D, ε, Φ, k ), for te ot and old fluds (T n, W, ρ, µ,, k, R, a, a 6 ) and te onfguraton parameters (,,, φ, Y, Y f ) are read. ) Te parameter Y assgns all flud data to sdes and. 3) umbers of annels per pass and are alulated dependng on te value of (see Tab. ). 4) ressure drops are alulated for bot sdes of te HE wt Eqs (), (3) and te equatons n Tab. (5). 5) oeffents α and α are obtaned by Eq. (9), usng te onstant overall eat transfer oeffent U, w s obtaned from Eq. (6) usng onstant eat transfer oeffents for sdes and, and. 6) Te values of s (, ) are determned. Tese depends on te onfguratons parameters,, and φ. 7) Te assemblng of te system of equaton starts wt te dfferental equatons for te HE nner annels (Eq. 8). 8) Dfferental equatons for te temperature n te frst and last annels are nluded n te system. 9) Boundary ondtons for te annels n sde are generated n te format sown n Tab. (4). Te flud pat nsde te exanger, from nlet to outlet, needs to be followed n order to determne te onnetons between annels. 0) Boundary ondtons for sde are also generated, but ea value of te parameter φ requres a spef treatment beause ts parameter determnes te flow dreton nsde annels. ) Equatons of termal effetveness for bot sdes (E and E ) are nluded n te dmensonless form of Eqs (4a) and (4b), w are obtaned usng Eqs (9) and (). E.max α α,. θ n θ out, α E.max α α,. θ n θ out (4) α ) Te resultng system of dfferental equatons, defned by te dfferental equatons on te temperature n ea annel, te boundary ondtons equatons and te effetveness equatons, s solved by numeral or analytal metods. 3) f a numeral soluton metod s used, te overall energy onservaton an be verfed by te effetveness values E and E. f E E > ε, were ε s te maxmum allowable devaton, te parameters of te numeral metod and/or te nput data sould be revsed. 4) Te man smulaton results, su as te pressure drops and outlet temperatures for sdes and (w are assgned to te ot and old streams usng Y ), as well as te HE termal effetveness, are obtaned. A omputer program was developed to run steps troug of te assemblng algortm n Fg. (6). After readng te data, te program makes all neessary alulatons, generates a report and reates te formatted nput fle for groms. Ts proedure made te smulaton of dfferent onfguratons mu smpler and effetve. Te dervaton of te assemblng algortm for te rgorous model s stragtforward. Step 5 sould be removed and all te neessary equatons for te alulaton of te annel overall eat transfer oeffent U (Eq. 6), su as dmensonless numbers, onsttutve equatons and te flud pysal propertes dependene upon annel temperature, sould be nserted n steps 7 and 8. Equaton (4) sould also be wrtten n te form of Eq. () for te effetveness alulaton n step. Tese modfatons nrease te sze and omplexty of te system of equatons and make te soluton more dffult beause of te larger number of dsretzed varables. 6. Smulaton Example An example of a HE applaton s presented to sow te smulaton results. A ot stream of toluene exanges eat wt a old stream of benzene n a medum-szed HE wt 5 plates (see Fg. 7). Bot rgorous and smplfed models were smulated and te man results are presented n Tab. (6) and Fgs (8), (9) and (0). For te smplfed model, te flud average temperatures were estmated targetng a termal effetveness of 90%.

9 roeedngs of OBEM 00, Energy and Termal Systems, Vol.4, 384 Table 6. Man smulaton results for te example Varable Rgorous Model Smplfed Model Devaton (%) HE termal effetveness (%) Benzene outlet temp. ( ) Toluene outlet temp. ( ) Benzene pressure drop (ka) Toluene pressure drop (ka) umber of varables after dsretzaton 4, U tme on workstaton (s) ,700,650 U 3,600 U7 U U6 U (W/m.K),550 U, U 9,500 U8 U Fgure 7. aratersts of te HE used for smulaton example,400,450 Rgorous model Smplfed model η Fgure 8. Dstrbuton of te overall eat transfer oeffent along te exanger 75.0 Sde - Benzene 75.0 Sde - Toluene nd ass st ass T ( o ) T ( o ) st ass nd ass rd ass η Fgure 9. Temperature profles for sde annels (rgorous model) η Fgure 0. Temperature profles for sde annels (rgorous model) Te temperature dstrbuton along te exanger, obtaned by te rgorous model, s sown n Fgs (8) and (9), were te stream passes and annel numbers are ndated. Te temperature varaton n te outer annels and 4 are lower tan n te oter annels n te same pass beause eat s exanged wt only one negbor annel. Te dstrbuton of te overall eat transfer oeffent U along te exanger, obtaned by bot models, s presented n Fg. (8). Te results from te rgorous model smulaton sowed tat te oeffent U vares from,464 to,665 W/m.K, wle te smplfed model was solved wt an average value of,56 W/m.K. Despte ts sgnfant

10 dfferene, te man smulaton results obtaned by bot models are very lose (Tab. 6), wt a devaton of only 0.35% for te exanger effetveness. Several examples ave sown tat, even wt a remarkable dfferene n te dstrbuton of te overall eat transfer oeffent, tere s lttle dfferene between te man smulaton results obtaned from bot rgorous and smplfed models. Te obtaned devatons were under.5%, wt respet to E. Te number of varables after te dsretzaton by te fnte dfferene metod s approxmately 9 tmes larger for te rgorous model. To make te soluton of large exangers relable, te blok deomposton opton of groms s used, were te system of equatons s deomposed nto smaller systems. Even wt ts opton, te soluton of large exangers (e.g. 50 plates) usng te rgorous model requres 3 to 5 mn n a DE-UX workstaton. For te ase of smaller exangers, as te one presented n ts paper, te soluton s aeved n a matter of seonds (Tab. 6). 7. onlusons Te onfguraton of a plate eat exanger (HE) was araterzed by a set of sx parameters and a metodology to detet equvalent onfguratons was presented. Based on ts parameterzaton, a detaled model for te smulaton of a HE n steady state wt a general onfguraton was developed n algortm form. Te developed assemblng algortm made te smulaton and omparson of dfferent onfguratons more flexble. An mportant feature of te proposed algortm s tat t may be oupled to any proedure to solve te system of dfferental and algebra equatons. Te assumpton of onstant overall eat transfer oeffent along te exanger, often used for te matematal modelng, was tested and sowed lttle nfluene over te man smulaton results for eat exange (termal effetveness and outlet temperatures). Te presented algortm s an mportant tool for te study of te nfluene of te onfguraton over te exanger performane, and an be furter used to develop an optmzaton metod for seletng te pleat eat exanger onfguraton. A resear n ts subjet s urrently n development. 8. Aknowledgments Te autors would lke to tank te fnanal support from FAES (grants 98/5808- and 00/3635-4). 9. Referenes Bassouny, M.K.; Martn, H., Flow Dstrbuton and ressure Drop n late Heat Exangers -, U-Type Arrangement, emal Engng Sene, v.39, n.4, pp , 984. Buonopane, R.A.; Troupe, R.A.; Morgan, J.., Heat Transfer Desgn Metod for late Heat Exangers, emal Engng rogress, v.57, n.7, pp.57-6, July 963. Georgads, M..; Rotsten, G.E.; Maetto, S., Modelng and Smulaton of omplex late Heat Exanger Arrangements under Mlk Foulng omputers & em. Engng, v. sups, Mar/Aprl, pp.s33-s338, 998. Kakaç, S.; Lu, H., Heat Exangers: Seleton, Ratng and Termal Desgn, R ress, ew York, 998. Kandlkar, S.G.; Sa, R.K., Multpass late Heat Exangers - Effetveness-TU Results and Gudelnes for Seletng ass Arrangements, ASME Journal of Heat Transfer, v., pp , 989a. Kandlkar, S.G.; Sa, R.K., Asymptot Effetveness-TU Formulas for Multpass late Heat Exangers, ASME Journal of Heat Transfer, v., pp.34-3, 989b. Masubu, M.; to, A., Dynam Analyss of a late Heat Exanger System, Bulletn of te JSME, v.0, n.4, pp , Aprl 977. MKllop, A.A.; Dunkley, W.L., late Heat Exangers: Heat Transfer, ndustral and Engng emstry, v.5, n.9, pp , Sept gnott, A.; Tamborenea,.., Termal Effetveness of Multpass late Exangers, nt. Journal of Heat and Mass Transfer, v.3, n.0, pp , Ot roess Systems Enterprse Ltd, groms ntrodutory User Gude, Release.0, London, 00. Saunders, E.A.D., Heat Exangers: Seleton, Desgn & onstruton, Longman S.&T., ew York, 988. Settar, A.;Venart, J.E.S., Approxmate Metod for te Soluton to te Equatons for arallel and Mxed-Flow Mult- annel Heat Exangers, nt. Journal of Heat Mass Transfer, v.5, pp.89-89, 97. Sa, R.K.; Foke, W.W., late Heat Exangers and ter Desgn Teory, n: Sa, R.K.; Subbarao, E..; Maselkar, R.A., Edtors, Heat Transfer Equpment Desgn, pp.7-54, Hemspere.., ew York, 988. Taborek, J., Sell-and-tube Heat Exangers: Sngle-pase Flow, n: Hewtt, G.F., Edtor, Handbook of Heat Exanger Desgn, s.3.3, Begell House, ew York, 99. Zalesk, T., A General Matematal-Model of arallel-flow, Multannel Heat-Exangers and Analyss of ts ropertes, emal Engng Sene, v.39, n.7/8, pp.5-60, 984. Zalesk, T.; Jarzebsk, A.B., Remarks on Some ropertes of Equaton of Heat-Transfer n Multannel Exangers, nt. Journal of Heat Mass Transfer, v.6, n.8, pp , 973. Zalesk, T.; Klepaka, K., Approxmate Metods of Solvng Equatons for late Heat-Exangers, nt. Journal of Heat Mass Transfer, v.35, pp.5-30, 99.

Gravity Drainage Prior to Cake Filtration

Gravity Drainage Prior to Cake Filtration 1 Gravty Dranage Pror to ake Fltraton Sott A. Wells and Gregory K. Savage Department of vl Engneerng Portland State Unversty Portland, Oregon 97207-0751 Voe (503) 725-4276 Fax (503) 725-4298 ttp://www.e.pdx.edu/~wellss

More information

Scientific Research of the Institute of Mathematics and Computer Science 1(11) 2012, 23-30

Scientific Research of the Institute of Mathematics and Computer Science 1(11) 2012, 23-30 Please te ts artle as: Grażyna Kałuża, Te numeral soluton of te transent eat onduton roblem usng te latte Boltzmann metod, Sentf Resear of te Insttute of Matemats and Comuter Sene,, Volume, Issue, ages

More information

Bidimensional Analysis of a Thermoelectric Module using Finite Element Techniques

Bidimensional Analysis of a Thermoelectric Module using Finite Element Techniques Bdmensonal Analyss of a Termoeletr Module usng Fnte Element Tenques *Antono Arenas, Jorge Vázquez, Rafael Palaos Unversdad Pontfa Comllas Esuela Téna Superor de ngenería *Departamento de Fludos y Calor,

More information

Stanford University CS254: Computational Complexity Notes 7 Luca Trevisan January 29, Notes for Lecture 7

Stanford University CS254: Computational Complexity Notes 7 Luca Trevisan January 29, Notes for Lecture 7 Stanford Unversty CS54: Computatonal Complexty Notes 7 Luca Trevsan January 9, 014 Notes for Lecture 7 1 Approxmate Countng wt an N oracle We complete te proof of te followng result: Teorem 1 For every

More information

Module 7: Solved Problems

Module 7: Solved Problems Mdule 7: Slved Prblems 1 A tn-walled nentr tube eat exanger f 019-m lengt s t be used t eat denzed water frm 40 t 60 at a flw rate f 5 kg/s te denzed water flws trug te nner tube f 30-mm dameter wle t

More information

Stress analysis by local integral equations

Stress analysis by local integral equations Boundary Elements and Oter Mes Reduton Metods XXIX 3 Stress analyss by loal ntegral euatons V. Sladek 1, J. Sladek 1 & C. Zang 2 1 Insttute of Construton and Arteture, Slovak Aademy of Senes, Bratslava,

More information

DYNAMIC ANALYSIS OF SEMI-RIGID FRAMES

DYNAMIC ANALYSIS OF SEMI-RIGID FRAMES Matematal and Computatonal pplatons, Vol., No., pp. -8, 5. ssoaton for Sentf esear DYNMIC NYSIS OF SMI-IGID FMS l Ugur Ozturk and Hkmet H. Catal Department of Cvl ngneerng, Dokuz ylul Unversty, 5, Izmr,

More information

Physics 41 Chapter 22 HW Serway 7 th Edition

Physics 41 Chapter 22 HW Serway 7 th Edition yss 41 apter H Serway 7 t Edton oneptual uestons: 1,, 8, 1 roblems: 9, 1, 0,, 7, 9, 48, 54, 55 oneptual uestons: 1,, 8, 1 1 Frst, te effeny of te automoble engne annot exeed te arnot effeny: t s lmted

More information

Multivariate Ratio Estimator of the Population Total under Stratified Random Sampling

Multivariate Ratio Estimator of the Population Total under Stratified Random Sampling Open Journal of Statstcs, 0,, 300-304 ttp://dx.do.org/0.436/ojs.0.3036 Publsed Onlne July 0 (ttp://www.scrp.org/journal/ojs) Multvarate Rato Estmator of te Populaton Total under Stratfed Random Samplng

More information

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850) hermal-fluds I Chapter 18 ransent heat conducton Dr. Prmal Fernando prmal@eng.fsu.edu Ph: (850) 410-6323 1 ransent heat conducton In general, he temperature of a body vares wth tme as well as poston. In

More information

Machine Learning: and 15781, 2003 Assignment 4

Machine Learning: and 15781, 2003 Assignment 4 ahne Learnng: 070 and 578, 003 Assgnment 4. VC Dmenson 30 onts Consder the spae of nstane X orrespondng to all ponts n the D x, plane. Gve the VC dmenson of the followng hpothess spaes. No explanaton requred.

More information

Influence of Gravity on the Performance Index of Microchannel Heat Exchangers-Experimental Investigations

Influence of Gravity on the Performance Index of Microchannel Heat Exchangers-Experimental Investigations Proeedngs of the World Congress on Engneerng 011 Vol III WCE 011, July 6-8, 011, London, U.K. Influene of Gravty on the Performane Index of Mrohannel Heat Exhangers-Expermental Investgatons Thanhtrung

More information

UNIVERSITY OF BOLTON RAK ACADEMIC CENTRE BENG(HONS) MECHANICAL ENGINEERING SEMESTER TWO EXAMINATION 2017/2018 FINITE ELEMENT AND DIFFERENCE SOLUTIONS

UNIVERSITY OF BOLTON RAK ACADEMIC CENTRE BENG(HONS) MECHANICAL ENGINEERING SEMESTER TWO EXAMINATION 2017/2018 FINITE ELEMENT AND DIFFERENCE SOLUTIONS OCD0 UNIVERSITY OF BOLTON RAK ACADEMIC CENTRE BENG(HONS) MECHANICAL ENGINEERING SEMESTER TWO EXAMINATION 07/08 FINITE ELEMENT AND DIFFERENCE SOLUTIONS MODULE NO. AME6006 Date: Wednesda 0 Ma 08 Tme: 0:00

More information

JSM Survey Research Methods Section. Is it MAR or NMAR? Michail Sverchkov

JSM Survey Research Methods Section. Is it MAR or NMAR? Michail Sverchkov JSM 2013 - Survey Researh Methods Seton Is t MAR or NMAR? Mhal Sverhkov Bureau of Labor Statsts 2 Massahusetts Avenue, NE, Sute 1950, Washngton, DC. 20212, Sverhkov.Mhael@bls.gov Abstrat Most methods that

More information

HEAT PERFORMANCE ANALYSIS OF MULTI-FLUID HEAT EXCHANGERS

HEAT PERFORMANCE ANALYSIS OF MULTI-FLUID HEAT EXCHANGERS HEA PERFORANCE ANALYSIS OF ULI-FLUID HEA EXCHANGERS A hess Submtted to the Faulty Of Purdue Unversty By Sang Bae Par Advsor: Dr. S. S. Krshnan /3/2003 HEA EXCHANGERS Deves that provde flow of thermal energy

More information

Multigrid Methods and Applications in CFD

Multigrid Methods and Applications in CFD Multgrd Metods and Applcatons n CFD Mcael Wurst 0 May 009 Contents Introducton Typcal desgn of CFD solvers 3 Basc metods and ter propertes for solvng lnear systems of equatons 4 Geometrc Multgrd 3 5 Algebrac

More information

ALGEBRAIC SCHUR COMPLEMENT APPROACH FOR A NON LINEAR 2D ADVECTION DIFFUSION EQUATION

ALGEBRAIC SCHUR COMPLEMENT APPROACH FOR A NON LINEAR 2D ADVECTION DIFFUSION EQUATION st Annual Internatonal Interdsplnary Conferene AIIC 03 4-6 Aprl Azores Portugal - Proeedngs- ALGEBRAIC SCHUR COMPLEMENT APPROACH FOR A NON LINEAR D ADVECTION DIFFUSION EQUATION Hassan Belhad Professor

More information

Controller Design for Networked Control Systems in Multiple-packet Transmission with Random Delays

Controller Design for Networked Control Systems in Multiple-packet Transmission with Random Delays Appled Mehans and Materals Onlne: 03-0- ISSN: 66-748, Vols. 78-80, pp 60-604 do:0.408/www.sentf.net/amm.78-80.60 03 rans eh Publatons, Swtzerland H Controller Desgn for Networed Control Systems n Multple-paet

More information

FAULT DETECTION AND IDENTIFICATION BASED ON FULLY-DECOUPLED PARITY EQUATION

FAULT DETECTION AND IDENTIFICATION BASED ON FULLY-DECOUPLED PARITY EQUATION Control 4, Unversty of Bath, UK, September 4 FAUL DEECION AND IDENIFICAION BASED ON FULLY-DECOUPLED PARIY EQUAION C. W. Chan, Hua Song, and Hong-Yue Zhang he Unversty of Hong Kong, Hong Kong, Chna, Emal:

More information

The Finite Element Method: A Short Introduction

The Finite Element Method: A Short Introduction Te Fnte Element Metod: A Sort ntroducton Wat s FEM? Te Fnte Element Metod (FEM) ntroduced by engneers n late 50 s and 60 s s a numercal tecnque for solvng problems wc are descrbed by Ordnary Dfferental

More information

Numerical Simulation of One-Dimensional Wave Equation by Non-Polynomial Quintic Spline

Numerical Simulation of One-Dimensional Wave Equation by Non-Polynomial Quintic Spline IOSR Journal of Matematcs (IOSR-JM) e-issn: 78-578, p-issn: 319-765X. Volume 14, Issue 6 Ver. I (Nov - Dec 018), PP 6-30 www.osrournals.org Numercal Smulaton of One-Dmensonal Wave Equaton by Non-Polynomal

More information

Tenth Order Compact Finite Difference Method for Solving Singularly Perturbed 1D Reaction - Diffusion Equations

Tenth Order Compact Finite Difference Method for Solving Singularly Perturbed 1D Reaction - Diffusion Equations Internatonal Journal of Engneerng & Appled Senes (IJEAS) Vol.8, Issue (0)5- Tent Order Compat Fnte Dfferene Metod for Solvng Sngularl Perturbed D Reaton - Dffuson Equatons Faska Wondmu Galu, Gemes Fle

More information

3D Numerical Analysis for Impedance Calculation and High Performance Consideration of Linear Induction Motor for Rail-guided Transportation

3D Numerical Analysis for Impedance Calculation and High Performance Consideration of Linear Induction Motor for Rail-guided Transportation ADVANCED ELECTROMAGNETICS SYMPOSIUM, AES 13, 19 MARCH 13, SHARJAH UNITED ARAB EMIRATES 3D Numeral Analss for Impedane Calulaton and Hgh Performane Consderaton of Lnear Induton Motor for Ral-guded Transportaton

More information

On Adaptive Control of Simulated Moving Bed Plants. Plants Using Comsol s Simulink Interface. Speaker: Marco Fütterer

On Adaptive Control of Simulated Moving Bed Plants. Plants Using Comsol s Simulink Interface. Speaker: Marco Fütterer daptve Smulated Movng ed Plants Usng Comsol s Smulnk Interfae Speaker: Maro Fütterer Insttut für utomatserungstehnk Otto-von-Guerke Unverstät Unverstätsplatz, D-39106 Magdeburg Germany e-mal: maro.fuetterer@ovgu.de

More information

Earlier Lecture. This gas tube is called as Pulse Tube and this phenomenon is called as Pulse Tube action.

Earlier Lecture. This gas tube is called as Pulse Tube and this phenomenon is called as Pulse Tube action. 31 1 Earlier Leture In te earlier leture, we ave seen a Pulse Tube (PT) ryoooler in wi te meanial displaer is removed and an osillating gas flow in te tin walled tube produes ooling. Tis gas tube is alled

More information

Computers and Geotechnics

Computers and Geotechnics Computers and Geotens 36 (2009) 1298 1307 Contents lsts avalable at SeneDret Computers and Geotens journal omepage: www.elsever.om/loate/ompgeo A system relablty approa for evaluatng stablty of rok wedges

More information

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(

More information

Adaptive Multilayer Neural Network Control of Blood Pressure

Adaptive Multilayer Neural Network Control of Blood Pressure Proeedng of st Internatonal Symposum on Instrument Sene and Tenology. ISIST 99. P4-45. 999. (ord format fle: ISIST99.do) Adaptve Multlayer eural etwork ontrol of Blood Pressure Fe Juntao, Zang bo Department

More information

Main components of the above cycle are: 1) Boiler (steam generator) heat exchanger 2) Turbine generates work 3) Condenser heat exchanger 4) Pump

Main components of the above cycle are: 1) Boiler (steam generator) heat exchanger 2) Turbine generates work 3) Condenser heat exchanger 4) Pump Introducton to Terodynacs, Lecture -5 Pro. G. Cccarell (0 Applcaton o Control olue Energy Analyss Most terodynac devces consst o a seres o coponents operatng n a cycle, e.g., stea power plant Man coponents

More information

The calculation of ternary vapor-liquid system equilibrium by using P-R equation of state

The calculation of ternary vapor-liquid system equilibrium by using P-R equation of state The alulaton of ternary vapor-lqud syste equlbru by usng P-R equaton of state Y Lu, Janzhong Yn *, Rune Lu, Wenhua Sh and We We Shool of Cheal Engneerng, Dalan Unversty of Tehnology, Dalan 11601, P.R.Chna

More information

2 Finite difference basics

2 Finite difference basics Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

GEL 446: Applied Environmental Geology

GEL 446: Applied Environmental Geology GE 446: ppled Envronmental Geology Watershed Delneaton and Geomorphology Watershed Geomorphology Watersheds are fundamental geospatal unts that provde a physal and oneptual framewor wdely used by sentsts,

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance Calculation for Conductors of Arbitrary Shape CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors

More information

Improving the Performance of Fading Channel Simulators Using New Parameterization Method

Improving the Performance of Fading Channel Simulators Using New Parameterization Method Internatonal Journal of Eletrons and Eletral Engneerng Vol. 4, No. 5, Otober 06 Improvng the Performane of Fadng Channel Smulators Usng New Parameterzaton Method Omar Alzoub and Moheldn Wanakh Department

More information

Characteristic Analysis of Exponential Compact Higher Order Schemes for Convection-Diffusion Equations

Characteristic Analysis of Exponential Compact Higher Order Schemes for Convection-Diffusion Equations Ameran Journal of Computatonal Matemats 0 9- do:0./ajm.0.00 Publsed Onlne June 0 (ttp://www.srp.org/journal/ajm) 9 Caraterst Analss of Eponental Compat Hger Order emes for Conveton-Dffuson Equatons Abstrat

More information

Solving Singularly Perturbed Differential Difference Equations via Fitted Method

Solving Singularly Perturbed Differential Difference Equations via Fitted Method Avalable at ttp://pvamu.edu/aam Appl. Appl. Mat. ISSN: 193-9466 Vol. 8, Issue 1 (June 013), pp. 318-33 Applcatons and Appled Matematcs: An Internatonal Journal (AAM) Solvng Sngularly Perturbed Dfferental

More information

Queueing Networks II Network Performance

Queueing Networks II Network Performance Queueng Networks II Network Performance Davd Tpper Assocate Professor Graduate Telecommuncatons and Networkng Program Unversty of Pttsburgh Sldes 6 Networks of Queues Many communcaton systems must be modeled

More information

DOAEstimationforCoherentSourcesinBeamspace UsingSpatialSmoothing

DOAEstimationforCoherentSourcesinBeamspace UsingSpatialSmoothing DOAEstmatonorCoherentSouresneamspae UsngSpatalSmoothng YnYang,ChunruWan,ChaoSun,QngWang ShooloEletralandEletronEngneerng NanangehnologalUnverst,Sngapore,639798 InsttuteoAoustEngneerng NorthwesternPoltehnalUnverst,X

More information

Phase Transition in Collective Motion

Phase Transition in Collective Motion Phase Transton n Colletve Moton Hefe Hu May 4, 2008 Abstrat There has been a hgh nterest n studyng the olletve behavor of organsms n reent years. When the densty of lvng systems s nreased, a phase transton

More information

ORDINARY DIFFERENTIAL EQUATIONS EULER S METHOD

ORDINARY DIFFERENTIAL EQUATIONS EULER S METHOD Numercal Analss or Engneers German Jordanan Unverst ORDINARY DIFFERENTIAL EQUATIONS We wll eplore several metods o solvng rst order ordnar derental equatons (ODEs and we wll sow ow tese metods can be appled

More information

Verification of a Chemical Nonequilibrium Flows Solver Using the Method of Manufactured Solutions

Verification of a Chemical Nonequilibrium Flows Solver Using the Method of Manufactured Solutions Avalable onlne at www.senedret.om SeneDret Proeda Engneerng 00 (014) 000 000 www.elsever.om/loate/proeda APISAT014, 014 Asa-Paf Internatonal Symposum on Aerospae Tenology, APISAT014 Verfaton of a Cemal

More information

technische universiteit eindhoven Analysis of one product /one location inventory control models prof.dr. A.G. de Kok 1

technische universiteit eindhoven Analysis of one product /one location inventory control models prof.dr. A.G. de Kok 1 TU/e tehnshe unverstet endhoven Analyss of one produt /one loaton nventory ontrol models prof.dr. A.G. de Kok Aknowledgements: I would lke to thank Leonard Fortun for translatng ths ourse materal nto Englsh

More information

Opportunities in Analytical Approaches to Spray Drying of Solid Dosage Forms

Opportunities in Analytical Approaches to Spray Drying of Solid Dosage Forms Opportuntes n Analytal Approahes to Spray Dryng of Sold Dosage Forms Dr. Renhard Vehrng Assoate Professor and George Ford Char n Materals Engneerng Unversty of Alberta, Department of Mehanal Engneerng

More information

Example

Example Chapter Example.- ------------------------------------------------------------------------------ sold slab of 5.5 wt% agar gel at 78 o K s.6 mm thk and ontans a unform onentraton of urea of. kmol/m 3.

More information

Turbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH

Turbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH Turbulence classfcaton of load data by the frequency and severty of wnd gusts Introducton Oscar Moñux, DEWI GmbH Kevn Blebler, DEWI GmbH Durng the wnd turbne developng process, one of the most mportant

More information

Dr. M. Perumal Professor & Head Department of Hydrology Indian Institute of Technology Roorkee INDIA Co-authors: Dr. B. Sahoo & Dr. C.M.

Dr. M. Perumal Professor & Head Department of Hydrology Indian Institute of Technology Roorkee INDIA Co-authors: Dr. B. Sahoo & Dr. C.M. Dr.. Perumal Professor & Head Department of Hdrolog Indan Insttute of Tehnolog Roorkee INDIA o-authors: Dr. B. Sahoo & Dr... Rao Dr. Dr... Perumal, Professor & & Head, Dept. Dept. of of Hdrolog, I.I.T.

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

Principles of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer

Principles of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer Prncples of Food and Boprocess Engneerng (FS 31) Solutons to Example Problems on Heat Transfer 1. We start wth Fourer s law of heat conducton: Q = k A ( T/ x) Rearrangng, we get: Q/A = k ( T/ x) Here,

More information

11/19/2013. PHY 113 C General Physics I 11 AM 12:15 PM MWF Olin 101

11/19/2013. PHY 113 C General Physics I 11 AM 12:15 PM MWF Olin 101 PHY 113 C General Pyss I 11 AM 12:15 PM MWF Oln 101 Plan or Leture 23: Capter 22: Heat engnes 1. ermodynam yles; work and eat eeny 2. Carnot yle 3. Otto yle; desel yle 4. Bre omments on entropy 11/19/2013

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

Analytical calculation of adiabatic processes in real gases

Analytical calculation of adiabatic processes in real gases Journal of Physs: Conferene Seres PAPER OPEN ACCESS Analytal alulaton of adabat roesses n real gases o te ths artle: I B Amarskaja et al 016 J. Phys.: Conf. Ser. 754 11003 Related ontent - Shortuts to

More information

Finite Element Analysis of the Stability of Tunnel Surrounding Rock with Weak Rock Layer

Finite Element Analysis of the Stability of Tunnel Surrounding Rock with Weak Rock Layer Vol., No. 2 Modern Appled Sene Fnte Element Analyss of the Stablty of Tunnel Surroundng Rok wth Weak Rok Layer Yangsong Zhang Nanjng Unversty of Sene and Tehnology, Nanjng 294, Chna Tel: 86-25-84-577 E-mal:

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

Voltammetry. Bulk electrolysis: relatively large electrodes (on the order of cm 2 ) Voltammetry:

Voltammetry. Bulk electrolysis: relatively large electrodes (on the order of cm 2 ) Voltammetry: Voltammetry varety of eletroanalytal methods rely on the applaton of a potental funton to an eletrode wth the measurement of the resultng urrent n the ell. In ontrast wth bul eletrolyss methods, the objetve

More information

The Similar Structure Method for Solving Boundary Value Problems of a Three Region Composite Bessel Equation

The Similar Structure Method for Solving Boundary Value Problems of a Three Region Composite Bessel Equation The Smlar Struture Method for Solvng Boundary Value Problems of a Three Regon Composte Bessel Equaton Mngmng Kong,Xaou Dong Center for Rado Admnstraton & Tehnology Development, Xhua Unversty, Chengdu 69,

More information

The optimal delay of the second test is therefore approximately 210 hours earlier than =2.

The optimal delay of the second test is therefore approximately 210 hours earlier than =2. THE IEC 61508 FORMULAS 223 The optmal delay of the second test s therefore approxmately 210 hours earler than =2. 8.4 The IEC 61508 Formulas IEC 61508-6 provdes approxmaton formulas for the PF for smple

More information

Module 3: Element Properties Lecture 1: Natural Coordinates

Module 3: Element Properties Lecture 1: Natural Coordinates Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers

More information

A Theorem of Mass Being Derived From Electrical Standing Waves (As Applied to Jean Louis Naudin's Test)

A Theorem of Mass Being Derived From Electrical Standing Waves (As Applied to Jean Louis Naudin's Test) A Theorem of Mass Beng Derved From Eletral Standng Waves (As Appled to Jean Lous Naudn's Test) - by - Jerry E Bayles Aprl 4, 000 Ths paper formalzes a onept presented n my book, "Eletrogravtaton As A Unfed

More information

425. Calculation of stresses in the coating of a vibrating beam

425. Calculation of stresses in the coating of a vibrating beam 45. CALCULAION OF SRESSES IN HE COAING OF A VIBRAING BEAM. 45. Calulaton of stresses n the oatng of a vbratng beam M. Ragulsks,a, V. Kravčenken,b, K. Plkauskas,, R. Maskelunas,a, L. Zubavčus,b, P. Paškevčus,d

More information

Lecture 26 Finite Differences and Boundary Value Problems

Lecture 26 Finite Differences and Boundary Value Problems 4//3 Leture 6 Fnte erenes and Boundar Value Problems Numeral derentaton A nte derene s an appromaton o a dervatve - eample erved rom Talor seres 3 O! Negletng all terms ger tan rst order O O Tat s te orward

More information

Direct Methods for Solving Macromolecular Structures Ed. by S. Fortier Kluwer Academic Publishes, The Netherlands, 1998, pp

Direct Methods for Solving Macromolecular Structures Ed. by S. Fortier Kluwer Academic Publishes, The Netherlands, 1998, pp Drect Metods for Solvng Macromolecular Structures Ed. by S. Forter Kluwer Academc Publses, Te Neterlands, 998, pp. 79-85. SAYRE EQUATION, TANGENT FORMULA AND SAYTAN FAN HAI-FU Insttute of Pyscs, Cnese

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI

More information

Chapter 7. Ab initio Theory

Chapter 7. Ab initio Theory Chapter 7. Ab nto Theory Ab nto: from the frst prnples. I. Roothaan-Hall approah Assumng Born-Oppenhemer approxmaton, the eletron Hamltonan: Hˆ Tˆ e + V en + V ee Z r a a a + > r Wavefunton Slater determnant:

More information

Title: Practical reliability-based calculation tool for the post-fire assessment of concrete beams

Title: Practical reliability-based calculation tool for the post-fire assessment of concrete beams Cover page Ttle: Pratal relablty-based alulaton tool for the post-fre assessment of onrete beams Authors (sample names): Ruben Van Cole Robby Caspeele Lu Taerwe PAPER DEADLINE: **JANUARY 15, 015** PAPER

More information

COMP4630: λ-calculus

COMP4630: λ-calculus COMP4630: λ-calculus 4. Standardsaton Mcael Norrs Mcael.Norrs@ncta.com.au Canberra Researc Lab., NICTA Semester 2, 2015 Last Tme Confluence Te property tat dvergent evaluatons can rejon one anoter Proof

More information

DEMO #8 - GAUSSIAN ELIMINATION USING MATHEMATICA. 1. Matrices in Mathematica

DEMO #8 - GAUSSIAN ELIMINATION USING MATHEMATICA. 1. Matrices in Mathematica demo8.nb 1 DEMO #8 - GAUSSIAN ELIMINATION USING MATHEMATICA Obectves: - defne matrces n Mathematca - format the output of matrces - appl lnear algebra to solve a real problem - Use Mathematca to perform

More information

On Pfaff s solution of the Pfaff problem

On Pfaff s solution of the Pfaff problem Zur Pfaff scen Lösung des Pfaff scen Probles Mat. Ann. 7 (880) 53-530. On Pfaff s soluton of te Pfaff proble By A. MAYER n Lepzg Translated by D. H. Delpenc Te way tat Pfaff adopted for te ntegraton of

More information

TR/95 February Splines G. H. BEHFOROOZ* & N. PAPAMICHAEL

TR/95 February Splines G. H. BEHFOROOZ* & N. PAPAMICHAEL TR/9 February 980 End Condtons for Interpolatory Quntc Splnes by G. H. BEHFOROOZ* & N. PAPAMICHAEL *Present address: Dept of Matematcs Unversty of Tabrz Tabrz Iran. W9609 A B S T R A C T Accurate end condtons

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017 U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that

More information

Key words: path synthesis, joint clearances, Lagrange s equation, Differential evaluation (DE), optimization.

Key words: path synthesis, joint clearances, Lagrange s equation, Differential evaluation (DE), optimization. Rub Mshra, T.K.Naskar, Sanb Ahara / nternatonal Journal of Engneerng Researh and Applatons (JERA) SSN: 8-96 www.era.om Vol., ssue, Januar -Februar 0, pp.9-99 Snthess of oupler urve of a Four Bar Lnkage

More information

Development of a computer model for long-throated flumes based on manning equation and different side slopes of trapezoidal channels

Development of a computer model for long-throated flumes based on manning equation and different side slopes of trapezoidal channels Iran Agrultural Researh (08) 37() Development of a omputer model for long-throated flumes based on mannng equaton and dfferent sde slopes of trapezodal hannels M. Mahbod * and Sh. Zand-Parsa Department

More information

Solution for singularly perturbed problems via cubic spline in tension

Solution for singularly perturbed problems via cubic spline in tension ISSN 76-769 England UK Journal of Informaton and Computng Scence Vol. No. 06 pp.6-69 Soluton for sngularly perturbed problems va cubc splne n tenson K. Aruna A. S. V. Rav Kant Flud Dynamcs Dvson Scool

More information

MULTICRITERION OPTIMIZATION OF LAMINATE STACKING SEQUENCE FOR MAXIMUM FAILURE MARGINS

MULTICRITERION OPTIMIZATION OF LAMINATE STACKING SEQUENCE FOR MAXIMUM FAILURE MARGINS MLTICRITERION OPTIMIZATION OF LAMINATE STACKING SEENCE FOR MAXIMM FAILRE MARGINS Petr Kere and Juhan Kos Shool of Engneerng, Natonal nversty of ruguay J. Herrera y Ressg 565, Montevdeo, ruguay Appled Mehans,

More information

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth

More information

MATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018

MATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018 MATH 5630: Dscrete Tme-Space Model Hung Phan, UMass Lowell March, 08 Newton s Law of Coolng Consder the coolng of a well strred coffee so that the temperature does not depend on space Newton s law of collng

More information

Chemical Engineering Department University of Washington

Chemical Engineering Department University of Washington Chemcal Engneerng Department Unversty of Washngton ChemE 60 - Exam I July 4, 003 - Mass Flow Rate of Steam Through a Turbne (5 onts) Steam enters a turbne at 70 o C and.8 Ma and leaves at 00 ka wth a qualty

More information

ECE 6340 Intermediate EM Waves. Fall Prof. David R. Jackson Dept. of ECE. Notes 3

ECE 6340 Intermediate EM Waves. Fall Prof. David R. Jackson Dept. of ECE. Notes 3 C 634 Intermedate M Waves Fall 216 Prof. Davd R. akson Dept. of C Notes 3 1 Types of Current ρ v Note: The free-harge densty ρ v refers to those harge arrers (ether postve or negatve) that are free to

More information

PHYSICS 212 MIDTERM II 19 February 2003

PHYSICS 212 MIDTERM II 19 February 2003 PHYSICS 1 MIDERM II 19 Feruary 003 Exam s losed ook, losed notes. Use only your formula sheet. Wrte all work and answers n exam ooklets. he aks of pages wll not e graded unless you so request on the front

More information

APPROXIMATE OPTIMAL CONTROL OF LINEAR TIME-DELAY SYSTEMS VIA HAAR WAVELETS

APPROXIMATE OPTIMAL CONTROL OF LINEAR TIME-DELAY SYSTEMS VIA HAAR WAVELETS Journal o Engneerng Sene and ehnology Vol., No. (6) 486-498 Shool o Engneerng, aylor s Unversty APPROIAE OPIAL CONROL OF LINEAR IE-DELAY SYSES VIA HAAR WAVELES AKBAR H. BORZABADI*, SOLAYAN ASADI Shool

More information

Clustering. CS4780/5780 Machine Learning Fall Thorsten Joachims Cornell University

Clustering. CS4780/5780 Machine Learning Fall Thorsten Joachims Cornell University Clusterng CS4780/5780 Mahne Learnng Fall 2012 Thorsten Joahms Cornell Unversty Readng: Mannng/Raghavan/Shuetze, Chapters 16 (not 16.3) and 17 (http://nlp.stanford.edu/ir-book/) Outlne Supervsed vs. Unsupervsed

More information

Prediction of the reliability of genomic breeding values for crossbred performance

Prediction of the reliability of genomic breeding values for crossbred performance Vandenplas et al. Genet Sel Evol 217 49:43 DOI 1.1186/s12711-17-318-1 Genets Seleton Evoluton RESERCH RTICLE Open ess Predton of the relablty of genom breedng values for rossbred performane Jéréme Vandenplas

More information

Natural Convection Experiment Measurements from a Vertical Surface

Natural Convection Experiment Measurements from a Vertical Surface OBJECTIVE Natural Convetion Experiment Measurements from a Vertial Surfae 1. To demonstrate te basi priniples of natural onvetion eat transfer inluding determination of te onvetive eat transfer oeffiient.

More information

COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN

COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN Int. J. Chem. Sc.: (4), 04, 645654 ISSN 097768X www.sadgurupublcatons.com COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN R. GOVINDARASU a, R. PARTHIBAN a and P. K. BHABA b* a Department

More information

modeling of equilibrium and dynamic multi-component adsorption in a two-layered fixed bed for purification of hydrogen from methane reforming products

modeling of equilibrium and dynamic multi-component adsorption in a two-layered fixed bed for purification of hydrogen from methane reforming products modelng of equlbrum and dynamc mult-component adsorpton n a two-layered fxed bed for purfcaton of hydrogen from methane reformng products Mohammad A. Ebrahm, Mahmood R. G. Arsalan, Shohreh Fatem * Laboratory

More information

Lecture Notes on Linear Regression

Lecture Notes on Linear Regression Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

An Application of the IMSC on a Non-linear Flexible Structure: Numerical Analysis and Experimental Validation

An Application of the IMSC on a Non-linear Flexible Structure: Numerical Analysis and Experimental Validation , July 6-8, 2011, London, U.K. An Applaton of the IMSC on a Non-lnear Flexble Struture: Numeral Analyss and Expermental Valdaton G. Caulan, F. Resta, F. Rpamont Abstrat he ndependent modal ontrol to suppress

More information

Comparison of Newton Raphson and Hard Darcy methods for gravity main nonlinear water network

Comparison of Newton Raphson and Hard Darcy methods for gravity main nonlinear water network IOSR Journal of Mecancal and Cvl Engneerng (IOSRJMCE) eissn:,pissn: 0X, Volume, Issue Ver. III (May. June. 0), PP 0 www.osrjournals.org Comparson of Newton Rapson and Hard Darcy metods for gravty man nonlnear

More information

CME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13

CME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13 CME 30: NUMERICAL LINEAR ALGEBRA FALL 005/06 LECTURE 13 GENE H GOLUB 1 Iteratve Methods Very large problems (naturally sparse, from applcatons): teratve methods Structured matrces (even sometmes dense,

More information

A New Recursive Method for Solving State Equations Using Taylor Series

A New Recursive Method for Solving State Equations Using Taylor Series I J E E E C Internatonal Journal of Electrcal, Electroncs ISSN No. (Onlne) : 77-66 and Computer Engneerng 1(): -7(01) Specal Edton for Best Papers of Mcael Faraday IET Inda Summt-01, MFIIS-1 A New Recursve

More information

Adaptive Kernel Estimation of the Conditional Quantiles

Adaptive Kernel Estimation of the Conditional Quantiles Internatonal Journal of Statstcs and Probablty; Vol. 5, No. ; 206 ISSN 927-7032 E-ISSN 927-7040 Publsed by Canadan Center of Scence and Educaton Adaptve Kernel Estmaton of te Condtonal Quantles Rad B.

More information

P REVIEW NOTES

P REVIEW NOTES P34 - REIEW NOTES Capter 1 Energy n Termal Pyss termal equlbrum & relaxaton tme temperature & termometry: fxed ponts, absolute temperature sale P = nrt deal gas law: ( ) ( T ) ( / n) C( T ) ( ) + / n vral

More information

THERMODYNAMICS Lecture 15: Heat exchangers

THERMODYNAMICS Lecture 15: Heat exchangers HERMODYNAMICS Leture 5: Heat exangers Pierwsza strona Introdution to Heat Exangers Wat Are Heat Exangers? Heat exangers are units designed to transfer eat from a ot flowing stream to a old flowing stream

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Not-for-Publication Appendix to Optimal Asymptotic Least Aquares Estimation in a Singular Set-up

Not-for-Publication Appendix to Optimal Asymptotic Least Aquares Estimation in a Singular Set-up Not-for-Publcaton Aendx to Otmal Asymtotc Least Aquares Estmaton n a Sngular Set-u Antono Dez de los Ros Bank of Canada dezbankofcanada.ca December 214 A Proof of Proostons A.1 Proof of Prooston 1 Ts roof

More information

Exercise 10: Theory of mass transfer coefficient at boundary

Exercise 10: Theory of mass transfer coefficient at boundary Partle Tehnology Laboratory Prof. Sotrs E. Pratsns Sonneggstrasse, ML F, ETH Zentrum Tel.: +--6 5 http://www.ptl.ethz.h 5-97- U Stoffaustaush HS 7 Exerse : Theory of mass transfer oeffent at boundary Chapter,

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information