Complete analytical model of a loop heat pipe with a flat evaporator

Size: px
Start display at page:

Download "Complete analytical model of a loop heat pipe with a flat evaporator"

Transcription

1 Complete nlytil model of loop het pipe with flt evportor B. Siedel, V. Srtre, F. Lefèvre To ite this version: B. Siedel, V. Srtre, F. Lefèvre. Complete nlytil model of loop het pipe with flt evportor. Interntionl Journl of Therml Sienes, Elsevier, 205, 89, pp <0.06/j.ijthermlsi >. <hl > HAL Id: hl Sumitted on Mr 206 HAL is multi-disiplinry open ess rhive for the deposit nd dissemintion of sientifi reserh douments, whether they re pulished or not. The douments my ome from tehing nd reserh institutions in Frne or rod, or from puli or privte reserh enters. L rhive ouverte pluridisiplinire HAL, est destinée u dépôt et à l diffusion de douments sientifiques de niveu reherhe, puliés ou non, émnnt des étlissements d enseignement et de reherhe frnçis ou étrngers, des lortoires pulis ou privés.

2 Complete nlytil model of loop het pipe with flt evportor Benjmin Siedel, Vlérie Srtre, Frédéri Lefèvre Université de Lyon, CNRS INSA-Lyon, CETHIL UMR5008, F-6962, Villeurnne, Frne Université Lyon, F-69622, Villeurnne, Frne Astrt A stedy-stte nlytil model hs een developed to determine the thermohydruli ehviour of loop het pipe with flt evportor. Its min originlity lies in the omintion of energy lne equtions for eh omponent of the system with 2D nlytil solutions for the temperture field in the evportor. Bsed on Fourier series expnsion, het trnsfer in the wik s well s in the evportor sing re urtely modelled, enling thorough onsidertion of the prsiti het fluxes. The model is sed on the therml ontt resistne etween the wik nd the sing, the therml ondutivity of the wik nd the ommodtion oeffiient. This nlytil method offers simple solution tht n e implemented in LHP design nlysis without the need of lrge omputtionl resoures. A sensitivity nlysis hs een rried out to evlute the influene of severl prmeters on the LHP ehviour. The results show tht the min prmeters of the model re independent. Therefore, they ould e experimentlly determined using n pproprite test enh with only few temperture mesurements. The model hs een vlidted with set of experimentl dt from the literture. A good greement is found etween the theoretil nd the experimentl results. Keywords: Loop Het Pipe, Anlytil, Modeling, Sensitivity nlysis. Introdution Loop Het Pipes LHP re two-phse ooling systems le to pssively trnsport high mounts of het over distnes up to severl meters. Developed in the 970 s, these devies hve proven their reliility in mny sptil pplitions nd re now ndidtes for terrestril ooling solutions. Indeed, their speifi design offers roustness nd flexiility for wide vriety of prtil pplitions [, 2]. As onsequene, mny efforts hve een dedited to understnd their opertion in order to optimize their design. A LHP is omplex therml system inluding n evportor onneted to the het soure, ondenser to dissipte the het lod nd vpour nd liquid lines to trnsport the working fluid etween oth omponents. Coupled thermo-hydruli phenomen govern the LHP ehviour nd need to e understood to enle orret system designing. Prmeters suh s evportion effiieny, het losses to the mient nd prsiti het flux in the evportor s well s ondenstion het trnsfer n e of gret influene on the loop opertion. A lot of ppers onerning LHP omplete modelling n e found in the literture [36]. However, these numeril models often imply omplex lgorithms nd lrge omputtionl resoures whih re not lwys ville for upstrem pre-design pplitions. An nlytil model offers Corresponding uthor Emil ddress: vlerie.srtre@ins-lyon.fr Vlérie Srtre the dvntge to provide solution without n exessive omputing time nd tht n e esily implemented. Some nlytil models of LHP re found in reserh works. Mydnik et l.[7] ited y Luny et l.[8] suggested losed form solution of LHP nlytil model estlishing the energy lne in the reservoir nd the pressure lne in the whole system. Furukw[9] developed omplete nlytil model of LHP le to enhne the sizing of the system nd to study the influene of mny geometril prmeters on the loop opertion. However, the model requires the operting temperture s input prmeter. Yet in most ses, the evportor temperture is the min expeted output of LHP model. Luny et l.[8] proposed losed-form expressions of the operting temperture of the LHP, for oth the vrile nd the fixed ondutne mode. Their model is sed on energy lne equtions on eh system omponent nd on thermodynmi equtions. The therml links in the reservoir re defined s equivlent therml resistnes. Their solution is useful tool in the LHP design. However, the het trnsfer in the evportor is not urtely determined nd hs to e djusted in ordne with experimentl dt. The purpose of the present study is to present omplete nlytil model of Loop Het Pipe urtely tking into ount het nd mss trnsfer in the evportor struture. The model is developed for flt disk-shped evportor geometry. This model is sed on the nlytil study of Luny et l.[8]. However, the present study rests upon two 2-D nlytil solutions desriing the tem- Preprint sumitted to Interntionl Journl of Therml Sienes Jnury 0, 204

3 perture field in the wik nd the evportor sing. Both solutions enle to determine prsiti het losses through the wik nd the evportor ody, the sensile het given to the liquid flowing through the porous struture s well s the het dissipted y evportion t the wik-groove interfe. A similr pproh ws implemented in the se of onventionl het pipes y Lefèvre nd Lllemnd [0] nd lter extended y Lips nd Lefèvre []. These fetures, oupled with energy lnes nd thermodynmi reltionships in the rest of the LHP, give simple solution for the operting temperture. An itertive proedure is implemented to lulte the two-phse length in the ondenser. T ext Q ext,e Body Evportor/Reservoir Q in Q Qw Qev T e T we Wik T v Condenser two-phse T sink 2. Model desription T ext Q ext,r T r Q su P 2.. Anlytil Model of the LHP The therml stte of the omplete LHP n e determined using energy lne equtions nd thermodynmi reltionships. Figure presents the operting priniple of the LHP nd the links etween its omponents. The totl het lod to e dissipted y the evportor Q in is onduted through the wik or through the evportor ody so tht: Q in = Q w +Q T ext T r,in Liquid line P Thermodynmi reltionship Fluid therml resistne Condution The wik is ssumed to e fully sturted with liquid. The therml het flux Q w is trnsverslly onduted through the evportor wll, the wll-wik interfe nd then dividesup: prtq ev isevportedtthewik-grooveinterfe wheres the rest is dissipted y ondution nd onvetion with the liquid flowing through the porous struture nd with the liquid in the reservoir. Q is onduted longitudinlly through the evportor wll to the reservoir nd prt of it, Q ext,e, is given y onvetion to the mient. Both the het flux through the wik Q w nd the het flux onduted through the evportor sing Q re funtions of the reservoir, the groove, the wik nd the evportor tempertures T r, T v, T we nd T e. The sme dependene pplies for Q ev nd Q ext,e : Q w = ft r,t v,t e 2 Q ev = ṁ l h lv = ft r,t v,t we 3 Q = ft r,t e 4 Q ext,e = ft r,t e 5 Q ext,e is lso funtion of T ext, whih is given dt of our model. As result, the het lod Q in n lso e expressed s funtion of these four tempertures: Q in = ft r,t v,t we,t e 6 An nlytil expression of Q in will e derived in susetions 2.2 nd 2.3. The prt of Q w tht is not dissipted y evportion is the trnsversl prsiti het flux. Prt of this flux is onduted through the wik nd relesed to the 2 T sink T,o T v Condenser suooling Convetion Phse-hnge Contt Figure : LHP shemti nodl network reservoir wheres the rest is dissipted y onvetion due to the liquid flow inside the porousstruture. At the interfe etween the wik nd the evportor envelope, there is temperture gp T e -T we due to ontt resistne R defined s: R = S T e T we Q w 7 where S is the ontt surfe etween the wik nd the evportor ody nd T e nd T we re the tempertures on the envelope side nd on the wik side, respetively. A glol het lne on the evportor/reservoir gives the following eqution: Q in = Q ev +Q sen +Q su +Q ext,e +Q ext,r 8 where Q sen is the sensile het given to the liquid, Q su is the suooling due to the liquid entering the reservoir

4 nd Q ext,r is the het flux dissipted to the mient y the reservoir. The determintion of Q sen nd Q su leds to: Q sen = ṁ l p,l T v T r 9 Q su = ṁ l p,l T r T r,in 0 where T r,in is the temperture of the liquid oming from the ondenser nd flowing k to the reservoir. To evlute the het trnsfer given y the reservoir to the mient, it is ssumed tht its surfe is t uniform temperture equl to T r. Het trnsfer with the mient Q ext,r is then pproximted y: Q ext,r = h ext S r T r T ext where S r is the totl externl surfe of the reservoir. A single expressionof Q in n e derived y omining equtions 3, 5 nd 7 to : Q in = ft r,t v,t we,t e,ṁ l,t r,in 2 In some ses, the reservoirn e full of liquid. This phenomenon hs een extensively studied y Adoni et l.[2]. In the present model, the existene of two-phse equilirium in the reservoir is ssumed. As shown y Luny et l. [8], thermodynmi reltionship exists etween the sturtion temperture inside the grooves nd the one t the liquid-vpour interfe in the reservoir: T = T v T r = T P P v + P l ρ l g H 3 where ρ l is the liquid density nd H is the elevtion of the ondenser ompred to the evportor. The slope of the pressure-temperture sturtion urve is given y the Clusius-Clpeyron eqution: T P = T /ρ v /ρ l h lv 4 The model n lso ope with non-ondensle gses NCG, whih n e generted for vrious resons in the LHP nd umulte in the reservoir, s explined y Singh et l. [3]. We ssume tht in operting onditions, these NCG re drined to the reservoir. In order to tke into ount the overpressure generted y the NCG P NCG, eqution 3 is modified: T = T v T r = T P P v + P l ρ l g H +P NCG 5 The vpour line is onsidered diti so tht the vpour enters the ondenser with temperture equl to T v. Furthermore, the ondenstion temperture T nd the vpour temperture T v re linked with the thermodynmi reltionship 4. Sine the pressure drops in the vpour line re low, it is ssumed tht ondenstion ours t temperture T v T T v. In the prt of the ondenser where liquid suooling ours, het trnsferwiththe hetsinkndwiththe mient 3 re lulted onsidering onvetive het trnsfer h sink nd h ext, respetively: T,o = T sink +T v T sink 6 πd,i L L 2φ exp ṁ l p,l /h l +D,i /h sink D,o wheret,o isthetempertureoftheliquidttheondenser outlet, L nd L 2φ re the lengths of the ondenser nd of the two-phse region respetively. D,i nd D,o re the inner nd outer ondenser dimeters. In similr wy, the liquid line het lne is: T r,in = T ext +T,o T ext 7 πd l,i L l exp ṁ l p,l /h l +D l,i /h ext D l,o where D l,i nd D l,o re the inner nd outer dimeters of the liquid line nd L l its length. Additionlly, in the ondenser, the het exhnge with the het sink in the twophse zone is equl to the ltent het of the ondensing vpour: ṁ l h lv = πl 2φ T v T sink 8 h ond D,i + h sink D,o Equtions 3, 7, 8, 2, 5 to 7 nd 8 form set of 8 independent equtions with 8 unknowns: T we, T e, T r, T r,in, T,o, T v, ṁ l nd L 2φ. Their solution leds to the determintion of the omplete therml stte of the LHP. The vlues of Q, Q w nd Q ev equtions 2-4 n e lulted using severl methods. In the lrge mjority of LHP models from the literture, het trnsfer in the evportor is simply desried using equivlent therml resistnes sed on the geometril hrteristis nd the thermophysil properties of the evportor. However, suh method does not tke into ount dequtely the het flux in the wik nd the evportor ody sine the determintion of the therml resistnes requires 2D or 3D pproh. Thus, in the present study, n urte therml nlysis of the evportor hs een onduted, sed on n nlytil pproh Anlytil Therml Model of the Wik The first prt of the nlytil model of the evportor dels with het nd mss trnsfer inside the porous struture. As shown in Figure 2, prt of the porous wik is modelled, ordered on one side y the liquid ulk of the reservoir nd y hlf of fin nd hlf of groove on the other side. The 2-D sttionry het eqution in the wik is expressed s: 2 T w x T w y 2 = 0 9 where T w is the temperture of the porous struture nd x nd y re the xis oordintesfig. 3. A non-dimensionl temperture is defined s: T w = λ efft w T r φ 0 20

5 Reservoir vpour liquid wik groove Thus eqution 23 eomes: T wx,y = A m YosmπX 25 m=0 Comining equtions 2 nd 25 yields: mπb 2 A m Y+ 2 A m Y Y 2 = 0 26 Solving the previous differentil eqution leds to: A 0 Y = A 0 Y +A 02 if m = 0 27 wll A m Y = A m e mπby +A m2 e mπby otherwise 28 nd leds to: with X = x, Figure 2: Evportor ross-setion 2 T w X 2 + B 2 2 T w Y 2 = 0 2 Y = y nd B = 22 where nd re the lengths of the modelled region in the x nd y diretions respetively, λ eff is the wik effetive therml ondutivity nd φ 0 is n ritrry het flux. y T r 0 T we 0 T v x Figure 3: Shemti of the wik model A generl expression of the non-dimensionl temperture field expnded in 2-D Fourier series is given y: Tw X,Y = A m YosmπX+B m YsinmπX m=0 23 For X = 0 nd X =, n diti oundry ondition is onsidered sed on symmetry hypothesis: T w X = T w X=0 X = 0 24 X= 4 The oundry ondition for Y = 0 is set temperture profile orresponding to the evportor temperture T we for one side x 0 nd to the groove temperture T v for the other side x. At the juntion etween the wll nd the groove, the temperture singulrity is treted y onsidering liner vrition of the temperture etween two fitive points 0 nd. The slope of the temperture grdient nd thus the position of these points depend on the mximum evportion rte t the liquid-vpour interfe tht n e lulted using the kineti gs theory, whih gives het trnsfer oeffiient h ev [4]: h ev = ev 2 ev ρ v h 2 lv T st 2πRTst M 0.5 P st 2ρ v h lv 29 where ev is the ommodtion oeffiient. Suh n pproh ws lredy used in [5] to ope numerilly with the temperture singulrity t the triple line. In the se of the evportion of thin liquid film, the ommodtion oeffiient is defined s the rtio of the tul evportion rte to theoretil mximl phse hnge rte. A oeffiient equl to unity desries perfet evportion while lower vlue represents inomplete evportion. In the se of wter, vlues vrying from 0.0 to re suggested in the literture[6]. We ssume the following reltionship etween h ev nd the position of 0 nd : h ev = λ eff 0 30 The oundry ondition for Y = 0 is then defined s: λ eff T we T r if 0 < X 0 φ 0 λ eff Tw X,0 = T v T r +T we T v X φ 0 0 if 0 < X < λ eff T v T r φ 0 For Y =, the reservoir temperture is set: if X < 3 T wx, = 0 32

6 The non-dimensionl temperture field is then: Tw X,Y = A 0Y +A Am e mπby +A m2 e mπby osmπx with A 0 = λ eff T v T r +T we T v 0 + φ 0 A 02 = λ eff T v T r +T we T v 0 + φ 0 A m = 2 λ eff T we T v φ 0 m 2 π 2 A m2 = 2 λ eff φ 0 T we T v m 2 π 2 os mπ 0 os mπ e 2mπB 36 os mπ 0 os mπ 0 e 2mπB 37 In eqution 33, the liquid flow inside the porous struture is not tken into ount. This flow is twodimensionl. However, in the present study, we ssume D flow inside the wik, onsidering homogeneous volumetri soure inside the wik q defined s: q = 4ṁ l p,l T v T r πd 2 w 38 where ṁ l is the totl liquid mss flow rte in the wik nd p,l is the speifi het of the liquid. This ssumption respets the energy lne nd enles to derive n nlytil expression for the influene of the liquid flow inside the wik. The superposition priniple enles to dd simple model of het trnsfer in the wik with homogeneous soure to the previous developed nlysis. It hs to e noted tht 2D pproh of het nd mss trnsfer in porous wik ws presented y Co nd Fghri [7], ut the model is not entirely nlytil. The het eqution in the wik is thus: 2 T y 2 = q 39 λ eff nd the oundry onditions re: T y = 0 T y = = 2 q λ eff 40 T y = = T r 4 This leds to the non-dimensionl temperture field implied y the liquid flow inside the wik: T = q 2φ 0 Y 2 3Y This solution is dded to the previous one to give generl expression of the non-dimensionl temperture field in the wik, with onsidertion of -D onvetion: T t = T w +T 43 5 The previously desried model enles to lulte the het flux through the wik Q w nd the het dissipted y evportion Q ev : Q w = S w T t λ eff 0 y dx 44 y=0 [ 0 + = S w φ 0 A 0 3q 2φ 0 + BA m A m2 sin mπ ] 0 + Q ev = S w T t λ eff 0 + y dx 45 2 y=0 [ = S w φ A 0 3q 2φ 0 BA m A m2 sin mπ ] Anlytil Therml Model of the Evportor Body The model of the wik is suffiient if the het losses through the evportor ody re negligile. Nevertheless, the ody is usully mde of high ondutive mteril. Therefore, it is generlly neessry to estimte the het trnsferred y ondution from the evportor to the reservoir nd to the mient. A seond nlytil model is developed to desrie the evportor sing. As shown in Figure 4, the evportor wll is unwrpped nd modelled s retngulr domin. At x = 0 nd x =, n diti ondition is ssumed due to the symmetry. At y = 0, onvetive het trnsfer with the mient is tken into ount on the whole externl surfe, inluding the heting setion 0 x 0, where het flux φ in is lso pplied. The inner prt of the ody is highly influened y the reservoir temperture. The retngulr shpe hosen to represents the evportor ody does not onsider the effet of the grooves on the distortion of the temperture field etween the se plte nd the wik. We ssume, therefore, set temperture profile, with the evportor temperture T e on one side 0 x 0 nd the reservoir temperture T r on the otherside x with liner profile in-etween. The sum of two solutions is neessry to tke into ount the omplete set of oundry onditions. The first one orresponds to set temperture profile in the inner prt of the wll, in ontt with the wik nd the liquidvpour ulk in the reservoir. Convetive het losses to the mient re ssumed for the whole externl surfe of the evportor/reservoir. This hypothesis onsiders tht the eletroni omponent to ool down dissiptes het with the LHP on one side nd with the mient on the other. The seond solution dds the het input to the evportor.

7 vpour liquid wik T e φ in Q 0 h ext T ext y d T e T x,y T r 0 h T ext ext = + T r First solution 0 x y d T 2 x,y Seond solution T = φ = 0 x φ in Figure 4: Evportor sing modelling Using the sme mthemtil proedure s previously, the het eqution eomes: where: 2 T X 2 + D 2 2 T Y 2 = 0 46 X = x ; Y = y d ; D = d ; T = λ T T r φ 0 d 47 The oundry onditions for the first solution T re: T X T Y = T X=0 X = 0 symmetry 48 X= = Bi T + h ext T r T ext Y=0 φ 0 T X, = with λ T e T r φ 0 d Bi = dh ext λ if 0 < X 0 49 λ φ 0 d T e T r X if 0 0 X < 0 if < X < 50 wheres for the seond solution T2 : T2 X = T 2 X=0 X = 0 5 X= T2 φ in if 0 < X 0 Y = φ 0 Y=0 0 if X > 52 0 X, = 0 53 T 2 6 The omintion of equtions46 to53 nd the ddition of oth solutions led to the non-dimensionl temperture field in the evportor ody: TX,Y = C 0 Y +C Cm e mπdy +C m2 e mπdy osmπx with λ C 0 = Bi Bi +φ 0 d φ in T e T r 0 + +T r T ext φ 0 λ C 02 = T e T r 0 + Bi T r T ext Bi +φ 0 d φ in 56 φ 0 os mπ 0 os mπ 0 e mπd + mπd+bi mπd Bi e mπd mπ 0 +e 2mπD 57 os mπ 0 os mπ 0 e mπd + mπd Bi mπd+bi e mπd mπ 0 +e 2mπD 58 C m = 2 λ T e T r φ 0 d m 2 π 2 2 φ in φ 0 m 2 π 2 D sin C m2 = 2 λ T e T r φ 0 d m 2 π 2 +2 φ in φ 0 m 2 π 2 D sin The het dissipted to the mient in the evportor setion of the ody is: Q ext,e = 0 0 h ext πxt x,0 T ext dx 59 + h ext π 0 T x,0 T ext dx 0 = h ext π T r T ext + φ 0d λ C 02

8 +h ext π φ 0d λ mπ C m +C m2 mπ mπ0 +os mπ sin mπ The het trnsferred y therml ondution through the evportor ody to the reservoir setion is lulted y integrting the Fourier s lw t y = Fig. 4. Thus, the totl het losses through the ody is: d T Q = Q ext,e + π 0 λ 0 x dy 60 x= = Q ext,e +d 0 φ 0 πsin mπ C m e mπd C m2 e mπd 2.4. Solving proedure The solving proedure is presented in Figure 5. The set of equtions is not liner. Thus n itertive proedure is used to solve it. After initilistion of the prmeters of the model, the two-phse length in the ondenser L 2φ is set, ording to the energy lne for given het input Q in 8. L 2φ hs mjor influene on the determintion of the temperture of the liquid entering in the reservoir T r,in equtions 6-7. Then, the thermophysil properties re lulted, s well s the pressure drops in the trnsport lines. K-oeffiients n e defined to reformt the expression of Q w, Q ev, Q nd T r,in s funtions of T r, T v, T e nd T we : no Prmeter initiliztion, Q in Setting of the two-phse length in the ondenser L 2φ ording to Eq. 8 Clultion of the thermophysil properties nd the pressure drops Clultion of the K- oeffiients See Appendix A Determintion of the vpour temperture T v Eq. B. Is the energy lne in the ondenser Eq. 8 stisfied? yes Q in i + = Q in i + Q in Q in i + > Q mx yes Plot of the operting urve Figure 5: Solving lgorithm flowhrt no Q w = K T r +K 2 T v +K 3 T we using 34,36,37,44 6 Q ev = K 4 T r +K 5 T v +K 6 T we using 34,36,37,45 62 Q ext,e = K 0 T r +K T e +K 2 using 56,57,58,59 63 Q = K 7 T r +K 8 T e +K 9 using 57,58,60,63 64 T r,in = K 3 T v +K 4 using 6,7 65 The detiled expression of these oeffiients is presented in Appendix A. Equtions,7-, 5-7,6-65 re solved nd give seond-order expression whih enles to lulte the vpour temperture T v see Appendix B. This proedure is iterted until the energy lne is stisfied in the ondensereqution 8. The sme method is omputed for eh het input Q in i. Sine the het trnsfer oeffiient with the het sink is generlly muh lower thn the ondenstion het trnsfer oeffiient, the ondenstion therml resistne is negleted to simplify eqution 8. Eqution 5 is funtion of the pressure drops P v nd P l. These prmeters, due to the frition fores in 7 the vpour nd the liquid lines, depend on the fluid flow regime. They re lulted s follows: P = f 2ρD ṁ 2 L 66 whereais the ross-setionreofthe tue. Forsmooth tue wll of dimeter D nd length L, the frition ftor f is expressed y: 64/Re if Re 2000 f = if 2000 < Re < Re 0.25 if Re 950 A 67 Eqution 5 lso inludes the pressure of nonondensle gses P NCG. In order to tke into ount the NCG, it is neessry to lulte the liquid level in the reservoir tht depends on the het lod. We ssume tht the void frtion of the two-phse flow in the ondenser is 0.5 nd tht the vpour density is negligile ompred to the liquid density. The liquid level e l in the reservoir is

9 therefore expressed s: e l = S w m f εs w e w L lπdl,i 2 L 2 L 2φπD,i 2 ρ l where m f is the totl fluid hrge in the system nd e w is the wik thikness. The totl volume of NCG nd vpour in the reservoir V v is equl to: V v = S w e r e l 69 with e r the thiknessofthe reservoir. Consideringthenonondensle gses s idel gses, their prtil pressure is lulted y: P NCG = m NCGRT r M NCG V v 70 Extensive studies hve een undertken to develop models le to predit urtely the effetive therml ondutivity of porous struture. Thus, mny different orreltions n e found in the literture to determine the effetive therml ondutivity of porous mteril [224] Tle 2 presents the results otined for nikel wik of 75% porosity sturted with wter, using vrious orreltions. This thermophysil property depends not only on the ondutivity of the mterils onstituting the wik ut lso on geometril prmeters suh s the porosity, the men pore dimeter nd the pore size distriution. The results re very different ording to the hosen orreltion. This shows tht this prmeter is diffiult to evlute urtely. For the stndrd se, n effetive therml ondutivityequlto5w m K ishosen, whihorresponds to the wik properties of Singh et l.[9]. where R is the idel gs onstnt nd m NCG nd M NCG re the NCG totl mss nd the molr mss, respetively. To solve equtions 6-8, the het sink het trnsfer oeffiient h sink is determined in ordne with the ondenser design, wheres the het trnsfer oeffiient of the liquid h l is lulted ssuming lminr fully-developed flow nd onstnt Nusselt numer Nu D = 4.36 [8]. Equtions nd 7 depend on the het trnsfer oeffiient with the mient h ext. It is given y the orreltion of Churhill nd Chu [8] for free onvetion on the surfe of n isotherml ylinder: h ext = λ ir D R 6 D /Pr This orreltion is vlid for Ryleigh numers R D lower thn 0 2. The determintion of the effetive therml ondutivity of the porous struture will e disussed in the next setion. 3. Results nd disussion This setion presents sensitivity nlysis to show the influene of different prmeters on the LHP performne. The LHP geometry onsidered for this nlysis is similr to stndrd systems used for eletroni ooling pplitions. Its geometril hrteristis re sed on the experiments of Singh et l.[9] nd Choi et l.[20]. A vlidtion of the model is presented in the lst prgrph. 3.. Stndrd se The properties of this stndrd LHP, hving flt disk-shped evportor, is defined in Tle. The system is supposed to e in horizontl orienttion. The wik is mde of nikel nd the working fluid is wter. The min prmeters for the wik re the effetive therml ondutivity, the ommodtion nd the ontt therml resistne etween the wik nd the sing. 8 Correltion λ eff W m K Alexnder 5.82 Chudhry-Bndhri 4.06 Mxwell 6.96 Prllel sheme 23 Zehner-Shlunder.75 Tle 2: Effetive therml ondutivity lultion with wtersturted nikel wik, 75% porosity Informtion from the literture onerning the therml ontt etween porous struture nd solid se plte is very sre. This prmeter depends on mny geometril nd mnufturing hrteristis. Choi et l.[20] investigted new tehniques to enhne the therml ontt ondutne of evportors in LHPs. Severl hnnel designs hve een tested nd ompred. An estimtion of the ontt therml resistne inluding in prtie therml ondution in the se plte, in the wik nd n evportion resistne is given nd vlues rnging from to m 2 K W hveeen otined. Bsedon this work, onstnt ontt resistne of 0 4 m 2 K W is hosen for the stndrd se. Although the ommodtion oeffiient is key prmeter in the determintion of the LHP temperture, its vlue is very diffiult to predit. In the literture, severl studies re dedited to the determintion of this prmeter for vrious fluids[6]. However, the sttering of the results onfirms the diffiulty of desriing the evportion urtely. For wter, vlues rnging from 0.0 to hve een found. Therefore, n ommodtion oeffiient equl to 0. is hosen for the stndrd se. The tempertures t different lotions in the LHP re plotted in Figure 6 for het inputs rnging etween 0W nd 0 W. A temperture differene of severl degrees etween T e nd T we shows the impt of the therml ontt resistne etween the wik nd the sing, prtiulrly t high het inputs. Sine the mss flow rtes in the trnsport lines nd the ondenser re moderte due to the high

10 Evportor design Flt disk-shped Wik dimeter 40 mm Wik thikness 3 mm Wik effetive therml ondutivity 5W m K Vpour groove width mm Wik-wll ontt resistne 0 4 K m 2 W Aommodtion oeffiient 0. Reservoir depth 0 mm Evportor dimeter 4 mm Evportor wll mteril Stinless steel Working fluid Wter Fluid hrge 7 g Condenser nd trnsport lines i/o dimeters 2/2.4 mm Trnsport lines length 200 mm Condenser length 00 mm Het sink temperture 22 C Het trnsfer oeffiient with the het sink 2000W m 2 K Amient temperture 22 C Het trnsfer oeffiient with the mient 5W m 2 K Tle : Definition of the stndrd LHP for the sensitivity nlysis ltent het of vporiztion of wter, pressure drops in the loop re not signifint. Therefore, the sturtion tempertures in the reservoir nd in the grooves, T r nd T v, re lmost equl. The vlues of T,o nd T r,in re lso very similr euse the het trnsfer etween the liquid line nd the mient is limited. T C T e T we T v T r T,o T r,in Q in W Figure 6: Temperture of the LHP in the stndrd se A ler trnsition etween the vrile ondutne mode nd the fixed ondutne mode is oserved round 60W. For lower het inputs, the temperture of the liquid exiting the ondenser is onstnt nd equl to the het 9 sink temperture. Aove this point, T,o strts to inrese nd the shpe of ll the temperture urves in the LHP is qusi-liner. Figure 7 shows the temperture in the porous struture in the stndrd se, for totl het lod Q in of 50W. For y = 0mm, the temperture profile shows shrp disontinuity t the trnsition etween the groove nd the fin. The intense evportion ourring t x = 0.5mm uses the onvergene of the hetlines towrd this point nd shrp distortion of the lines in its neighourhood. As y inreses, the temperture profile flttens until the temperture is uniform in the x-diretion for y > 2mm. Then, the lterntion etween the fin nd the groove does not hve ny influene nd the temperture field in the wik in ontt with the liquid ulk of the reservoir is uniform. As defined in Eqution 3, the differene etween the groove temperture nd the reservoir temperture is set only y the pressure losses in the trnsport lines, the hydrostti pressure differene nd the non-ondensle gses prtil pressure. Smooth tues with reltively lrge dimeter led to redued pressure losses. Therefore, the LHP opertes with groove temperture nd reservoir temperture lmost equl. As onsequene, the onvetive ooling due to the liquid flow in the wik is extremely low nd het trnsfer in the porous struture is minly ontrolled y het ondution. Figure 8 shows the effet of non-ondensle gses on the evportor temperture. Vrious quntities of ir rngingetween ñgnd200 ñgresimultedin thereservoir. For NCG mss elow 0 ñg, orresponding to prtil pressure P NCG equl to out 400P, NCG do not hve determinnt influene on the evportor temperture. However, when the mss of NCG in the LHP is lrger, its prtil pressure is more importnt thn the sum

11 y mm C 82.2 C 83.2 C 79.7 C 80.2 C 80.7 C 8.7 C 82.7 C T r 79.7 C 80.2 C 80.7 C 80.7 C 80.2 C C 0 T we 0.5 T v x mm Figure 7: 2-D temperture field in the wik of the pressure drops in the whole loop. The influene of NCG eomes then importnt t low het input, leding to signifint inrese of the evportor temperture. As onsequene, the shpe of the hrteristi urve of the LHP is flttened. T e C No NCG µg NCG 0 µg NCG 00 µg NCG 200 µg NCG 400 µg NCG Q in W Figure 8: Influene of the NCGs on the evportor temperture These results re in ordne with Singh et l. s experimentl study [3]. Their work shows tht the net effet of the generted NCG in the LHP ws to produe n overll 0 rise in the stedy-stte operting temperture. Besides, it ws oserved tht the performne degrding effet of the NCG ws more pronouned t low het lods Sensitivity nlysis The LHP opertion depends on mny distint prmeters: the geometril design size, shpe of the system, the thermophysil properties of the working fluid nd of the mterils of the loop elements, the het trnsfer hrteristis inside nd outside the LHP nd the pressure losses in the system. Most of these prmeters re esy to determine if the LHP geometry nd the working fluid properties re known. However, some prmeters re very diffiult to determine preisely theoretilly s it hs een seen in the previous setion: the effetive therml ondutivity of the wik λ eff, the ontt resistne etween the wik nd the evportor envelope R nd the ommodtion oeffiient ev. Furthermore, the het sink het trnsfer oeffiient h sink nd the het trnsfer oeffiient with the mient h ext re lso lulted using orreltion eing inherently inurte. Thus, in the present setion, sensitivity nlysis is onduted on these prmeters to see their influene on the model. The sensitivity oeffiient of the funtion T = fx i,x j,x k,... in reltion to the prmeter x i is defined s: S i = T 72 x i xj,x k,... In order to ompre severl prmeter sensitivities, it is onvenient to define reltive sensitivities: Si T = x i S i = x i 73 x i xj,x k,... This oeffiient enles to quntify the vrition T used y reltive vrition x i /x i of the prmeter x i. The greter the solute vlue of the oeffiient, the more the funtion is sensitive to the prmeter. Unless otherwise mentioned, the sensitivity nlysis is onduted with the stndrd LHP defined in Tle. A slight vrition 5% of eh studied prmeter is pplied, the vlue of the other prmeters eing onstnt. This effet on the tempertures t different lotions of the LHP is predited y the model. Figure 9 presents the reltive sensitivity of T e, T v, T r,in nd T,o to the het trnsfer oeffiient with the mient. Sensitivities oft e nd T v resignifintforlowhet inputs nd derese with the het lod. Het trnsfer with the mient ools the loop down nd leds to lower opertionl temperture. As expeted, h ext does not hve ny influene on T r,i nd T,o in vrile ondutne mode euse their vlue is only set y the het sink. For higher het lods, the inrese of the het losses to the mient led to lower opertionl temperture. As onsequene, lrger prt of the ondenser is used to ondenste the vpour, the suooling length in the ondenser is smller nd T,o inreses. The sme effet is oserved on

12 T r,in, ut prtilly offset y the het losses from the liquid line to the mient. 0 0 dt/dh ext *h ext K T e T v T r,in T,o dt v /dx i *x i K λ eff R ev h sink Q in W Q in W Figure 0: Reltive sensitivity of T v to vrious prmeters Figure 9: Reltive sensitivity of the LHP tempertures to h ext In mny experimentl onfigurtions, n dequte therml insultion of the entire LHP enles to redue onsiderly the het losses to the mient. In the following, this prmeter is set to zero in order to etter highlight the effet of the other prmeters on the LHP. The reltive sensitivity of the vpour temperture to the other prmeters is shown in Figure 0. It is ler tht het trnsfer inside the evportor funtion of the prmeters λ eff, R nd ev governsthe opertion in vrile ondutne mode wheres t high het lods, h sink eomes the dominnt prmeter. Indeed, t low het lods, T,o equls T sink so the het trnsfer oeffiient with the het sink hs no influene on the suooling of the liquid nd on the LHP opertion in generl. In fixed ondutne mode, het trnsfer in the ondenser sets the LHP opertionl temperture nd the sensitivity of T v to h sink is liner. The prmeters ev, λ eff nd R hve moderte influene on the vpour temperture nd this influene dereses s the het lod inreses. The reltive sensitivity of T e to the sme prmeters is shown in Figure. The influene of the ommodtion oeffiient nd of the het trnsfer with the het sink is lmost the sme s for T v. However, the sensitivities to R nd λ eff hve different ehviour. While the ontt resistne hs limited effet on T v, its influene on T e is lrge nd inreses with the het input. Indeed, the differene etween T e nd T we is proportionl to the het trnsfer rte Q w nd to R. As the het input Q in inreses, Q w inreses lmost linerly nd leds to higher sensitivity of T e to the ontt resistne. The prtiulr shpe of the reltive sensitivity of T e to λ eff isdisussedinthe following; it requiresmoredetiled nlysisontheeffetofλ eff onthelhpthermlehviour. dt e /dx i *x i K λ eff R ev h sink Q in W Figure : Reltive sensitivity of T e to vrious prmeters Figure 2 shows the vrition of the vpour temperture with the effetive therml ondutivity of the wik for severl het lods. An optiml vlue of λ eff is found etween nd 2W m K whtever the het input. ThisminimumvlueofT v isonsequeneoftheevolution of the het trnsferred from the sing to the evportion zone. As the therml ondutivity of the wik inreses,

13 the het entering the wik eomes lrger t the expense of the longitudinl prsiti het flux. Figure 3 shows the distriution of het trnsfer in the evportor. The lrgest prt of the het lod Q in enters the wik Q w, wheres the rest is thermlly onduted through the evportor ody Q. When the vlue of λ eff is very low, inresing the ondutivity enhnes oth the evportion Q ev nd the trnsversl prsiti het flux Q w Q ev. When λ eff exeeds W m K, the inrese of the het flux entering the wik Q w is smller, leding to n inrese of the trnsversl prsiti het flux t the expense of Q ev, tht deresess λ eff eomes lrger. The mximum evportion het flux leds to minimum vpour temperture Figure 2. The sme onlusion hs een drwn in the numeril study of Siedel et l.[5]. Q/Q in % Q w Q ev Q w -Q ev Q Q=0W Q=30W Q=50W Q=70W Q=90W Q=0W λ W.m -.K - eff Figure 3: Distriution of het trnsfer in the evportor Q in = 50W T v C W m K in the stndrd se, the sign of the slope of the temperture urve hnges s the het flux inreses. As onsequene, the sign of the sensitivity of T e to λ eff hnges when inresing the het lod, s it is shown in Figure. This is not the se for T v, sine the minimum of the urves T v λ eff does not depend on the het flux Figure λ W.m -.K - eff Figure 2: Influene of λ eff on T v The evolution of the temperture differene T e T v with λ eff for severl het inputs is given in Figure 4. As expeted, this differene dereses when the therml ondutivity of the wik inreses. Indeed, higher vlue of λ eff leds to lower therml resistnes in the evportor Figure. Figure 5 presents the evolution of T e with λ eff, whih is onsequene of the results otined in figures 2 nd 4. An optiml effetive therml ondutivity exists, for whih T e is miniml. However, its vlue is lso dependent on the het input Q in, ontrry to the temperture of the vpour Figure 2. The optiml vlue of λ eff otined for the minimum vlue of T e inreses from 2W m K to 0W m K with the inrese of the het input from 0W to 0W. Therefore, there is not n optiml vlue of the therml ondutivity of the wik ut rnge of optiml vlues depending on the het lod. Figure 5 shows tht t given λ eff for exmple 2 T e -T v K Q=0W Q=30W Q=50W Q=70W Q=90W Q=0W λ W.m -.K - eff Figure 4: Influene of λ eff on T e T v This sensitivity nlysis shows tht the min prmeters of the model re independent. Therefore, their in-

14 T e C Q=0W Q=30W Q=50W Q=70W Q=90W Q=0W hnge of slope λ W.m -.K - eff Figure 5: Influene of λ eff on T e etween the experimentl results nd the lulted tempertures of the evportor wll nd of the vpour in the grooves. A good greement is found etween the experimentl dt nd the model. Severl unknown prmeters were identified y omprison of the predited vlues with the experimentl dt of Singh et l.: the het trnsfer oeffiient etween the ondenser wll nd the het sink, the evportor wll thikness, the ontt resistne nd the ommodtion oeffiient. As result, stright-tue equivlent ondenser is simulted with het trnsfer oeffiient h sink fixed to 3.2kW m 2 K, onsidering n outside dimeter of 2.4 mm for the tues. The ommodtion oeffiient is equl to 0.4. The vlue of R is set to 0 5 K m 2 W. The vlue of R tht enles to fit t est theresultsisverysmll, ut it hstoe notedtht T v is not relly the experimentl mesurement of the vpour temperture, ut the temperture of the tue t the exit of the evportor. Therefore, these experimentl results re not suffiient to estimte urtely the prmeters. Nevertheless, the order of mgnitude of the prmeters re lose to the one defined for the stndrd se. fluene n e differentited from eh other. Thus, the vilility of preise experimentl dt of severl representtive tempertures of the LHP for vrious het inputs my theoretilly led to preise determintion of these prmeters nd provide the model with dequte input prmeters. This sensitivity nlysis lso shows the lrge influene of the tested prmeters on the LHP opertion. Their inurte determintion n led to mjor error on the LHP opertion predition Model vlidtion The present nlytil model is ompred to n experimentl dt set from Singh et l. [9] for vlidtion purpose. These uthors studied the opertionl hrteristis of flt disk-shped evportor LHP, 30 mm in dimeter, using wter s working fluid. The 3 mm thik porous wik is mde of sintered nikel nd its therml effetive ondutivity is onsidered equl to out 5.87W m K using Alexnder s formul nd sed on Singh et l. s study[22]: T C T e experimentl T v experimentl T e model T v model Q in W ε 0.59 λl λ eff = λ l λ wm 74 Figure 6: Comprison etween the model nd dt from Singh et l. [9] where λ l nd λ wm re the therml ondutivities of the liquid nd of the wik mteril respetively nd ε is the porosity, equl to 75%. The wik is emedded in opper evportor. The vpour nd liquid lines, of internl dimeter 2mm, re50mm nd290mm longrespetively. Afinnd-tue ondenser, 50 mm long, dissiptes het y fored onvetion of ir t mient temperture i.e. 22 C. No prsiti het trnsfer through the evportor ody is tken into ount, sine n O-ring sel prevents het ondution to the reservoir. Figure 6 shows the omprison 3 4. Conlusion In this pper, omplete nlytil model of LHP hs een developed. Its originlity lies in the omintion of energy lne equtions for eh omponent of the system with nlytil solutions for the temperture field in the evportor. Bsed on Fourier series expnsions, het trnsfer in the wik s well s in the evportor sing re

15 urtely modelled. This nlytil method offers simple solution tht n e implemented in LHP design nlysis without the need of lrge omputtionl resoures. A sensitivity nlysis hs een undertken to ssess the influene of five prmeters on the loop opertion. This nlysis enles etter omprehension of the operting mehnisms of the LHP s well s omprtive study of the prmeters ffeting its tempertures. It ppers tht these prmeters n e experimentlly determined using n pproprite test enh with only few temperture mesurements. This model hs een vlidted with set of experimentl dt from the literture. A good greement hs een met etween the simultion nd the experimentl results. The model results show tht onvetion inside the wik does not ply mjor role nd n e negleted. Moreover, the temperture field in the wik is lmost uniform fr from the grooves. In ordne with experimentl dt from the literture, the presene of NCGs in the reservoir leds to n inrese of the evportor temperture. This degrding effet is more pronouned t low het lods. Nomenlture A Fourier series oeffiient ross-setionl re [m 2 ], 0, length [m] ev ommodtion oeffiient B length rtio Fourier series oeffiient length [m] C Fourier series oeffiient, 0, length [m] p speifi het [J.kg.K ] D dimeter [m] length rtio d length [m] e thikness [m] f frition ftor g grvittionl elertion [m.s 2 ] H height [m] h het trnsfer oeffiient [W.m 2.K ] h lv enthlpy of vporiztion [J.kg ] K oeffiient L length [m] M molr mss [kg.mol ] m, n Fourier series inrement m f totl fluid hrge [kg] ṁ mss flow rte [kg.s ] P pressure [P ] Q het trnsfer rte [W] q volumetri het soure [W.m 3 ] R ontt resistne [K.m 2.W ] S surfe re [m 2 ] S i solute sensitivity Si reltive sensitivity T temperture [K] T non-dimensionl temperture V volume [m 3 ] X, Y non-dimensionl oordintes x i,j,k sensitivity prmeter x, y xis oordintes [m] Greek Symols differene ε porosity λ therml ondutivity [W.m.K ] φ,φ 0 het flux [W.m 2 ] ρ density [kg.m 3 ] Susripts 2φ two-phse ir ir evportor ody ondenser, onvetive, ontt e evportor eff effetive ev evportion ext externl, mient i inner in input, inlet l liquid NCG non ondensle gs o outlet, outer r reservoir sen sensile sink het sink su suooling t totl inluding onvetion v vpour w wik we wik side of the wik-envelope interfe wm wik mteril Non Dimensionl Numers Bi Biot numer P r Prndtl numer R Ryleigh numer Re Reynolds numer Referenes [] Y. F. Mydnik, Miniture loop het pipes, in: 3 th Interntionl Het Pipe Conferene, 2004, pp [2] Y. Mydnik, Loop het pipes, Applied Therml Engineering

16 [3] T. Ky, J. Ku, T. T. Hong, M. K. Cheung, Mthemtil modeling of loop het pipes, in: 37 th AIAA Aerospe Sienes Meeting nd Exhiit, 999. [4] A. Delil, V. Bturkin, G. Gorenko, P. Gkl, V. Ruzykin, Modelling of miniture loop het pipe with flt evportor, in: Proeedings of the 32 nd Interntionl Conferene on Environmentl Systems, [5] P. Chung, An improved stedy-stte model of loop het pipes sed on experimentl nd theoretil nlyses, Ph.D. thesis, Pennsylvni Stte University [6] L. Bi, G. Lin, H. Zhng, D. Wen, Mthemtil modeling of stedy-stte opertion of loop het pipe, Applied Therml Engineering [7] Y. F. Mydnik, Y. G. Fershtter, N. N. Solodovnik, Loop het pipes: Design, investigtion, prospets of use in erospe tehnis, Teh. Rep. 9485, SAE Interntionl, Wrrendle, PA Apr [8] S. Luny, V. Srtre, J. Bonjour, Anlytil model for hrteriztion of loop het pipes, Journl of Thermophysis nd Het Trnsfer [9] M. Furukw, Model-sed method of theoretil design nlysis of loop het pipe, Journl of Thermophysis nd Het Trnsfer [0] F. Lefèvre, M. Lllemnd, Coupled therml nd hydrodynmi models of flt miro het pipes for the ooling of multiple eletroni omponents, Interntionl Journl of Het nd Mss Trnsfer [] S. Lips, F. Lefèvre, A generl nlytil model for the design of onventionl het pipes, to e pulished, Interntionl Journl of Het nd Mss Trnsfer. [2] A. Adoni, A. Amirjn, V. Jsvnth, D. Kumr, P. Dutt, Theoretil studies of hrd filling in loop het pipes, Journl of Thermophysis nd Het Trnsfer [3] R. Singh, A. Akrzdeh, M. Mohizuki, Opertionl hrteristis of the miniture loop het pipe with non-ondensle gses, Interntionl Journl of Het nd Mss Trnsfer [4] V. P. Crey, Liquid-vpor phse-hnge phenomen, Hemisphere, New York, 992. [5] B. Siedel, V. Srtre, F. Lefèvre, Numeril investigtion of the thermohydruli ehviour of omplete loop het pipe, Applied Therml Engineering [6] I. Emes, N. Mrr, H. Sir, The evportion oeffiient of wter: review, Interntionl Journl of Het nd Mss Trnsfer [7] Y. Co, A. Fghri, Anlytil solutions of flow nd het trnsfer in porous struture with prtil heting nd evportion on the upper surfe, Interntionl Journl of Het nd Mss Trnsfer [8] F. P. Inroper, D. P. DeWitt, Fundmentls of het nd mss trnsfer, John Wiley & Sons, 996. [9] R. Singh, A. Akrzdeh, M. Mohizuki, Opertionl hrteristis of miniture loop het pipe with flt evportor, Interntionl Journl of Therml Sienes [20] J. Choi, B. Sung, C. Kim, D.-A. Bor-Tsiu, Interfe engineering to enhne therml ontt ondutne of evportors in miniture loop het pipe systems, Applied Therml Engineering [2] S. Mo, P. Hu, J. Co, Z. Chen, H. Fn, F. Yu, Effetive therml ondutivity of moist porous sintered nikel mteril, Interntionl Journl of Thermophysis [22] R. Singh, A. Akrzdeh, M. Mohizuki, Effet of wik hrteristis on the therml performne of the miniture loop het pipe, Journl of Het Trnsfer [23] K. Boomsm, D. Poulikkos, On the effetive therml ondutivity of three-dimensionlly strutured fluid-sturted metl fom, Interntionl Journl of Het nd Mss Trnsfer [24] C. Li, G. P. Peterson, The effetive therml ondutivity of wire sreen, Interntionl Journl of Het nd Mss Trnsfer

17 Appendix A. Detiled expression of the K-oeffiients K = S w 0 + [ K 2 = S w 2 λ eff 0 m 2 π 2 sin λ eff 6 +ṁ l p,l πdw 2 mπ 0 + os λeff 0 + [ K 3 = S w 2 λ eff 0 m 2 π 2 sin λ ] 2 eff 0 + K 4 = S w 0 + λeff [ K 5 = S w 2 λ eff K 6 = S w [ + λ eff 0 + K 7 = h ext π 0 m 2 π 2 sin λeff 2 λ eff 0 m 2 π 2 sin ṁ l p,l mπ 0 + ṁ l p,l 6 πd 2 w mπ mπ 0 + ] Bi + mπ0 2λ 0 m 2 π 2 os [ mπ e 0 πsin mπd + πh ext πd 2 w A. mπ0 mπ e 2mπB os e 2mπB ] A.2 mπ0 os mπ0 os ] 6 ṁ l p,l os πd 2 w mπ0 os mπ mπ mπ λ e mπd +e mπdmπd+bi mπd Bi mπ e + 0 πsin mπd + πh ext λ e mπd +e +mπdmπd Bi ] mπd +Bi K 8 = h ext π Bi + mπ0 + 2λ 0 m 2 π 2 os [ mπ e 0 πsin mπd + πh ext λ e mπd +e mπdmπd +Bi os os os os +os +os mπ mπ mπ mπ mπ +os mπ e 2mπB e 2mπB A.3 A.4 mπ e 2mπB e 2mπB A.5 mπ e 2mπB e 2mπB mπ0 mπ + 0 sin mπ0 mπ + 0 sin mπ0 mπ + 0 sin mπd Bi mπ e + 0 πsin mπd + πh ext mπ0 mπ +os + 0 sin λ mπ mπ 6 A.6 A.7 A.8

18 e mπd +e +mπdmπd Bi ] mπd +Bi K 9 = h ext π 0 0 φ 0 d 0 φ in 2 λ φ 0 + 2φ in m 2 π 2 sin mπ0 [ mπ 0 πsin +e 2mπD mπ + 0 πsin +e 2mπD ] Bi + T ext e mπd πh ext λ e mπd + πh ext λ +os mπ mπ +os mπ mπ mπ0 mπ + 0 sin mπ0 mπ + 0 sin A.9 K 0 = h ext π Bi + mπ0 mπ 2λ 0 m 2 π 2 os os [ e mπd +e mπdmπd +Bi mπd Bi K = h ext π Bi + mπ0 mπ + 2λ 0 m 2 π 2 os os [ e mπd +e mπdmπd +Bi mπd Bi K 2 = h ext π 0 0 φ 0 d 0 φ in 2 λ φ 0 Bi + T ext 2φ in m 2 π 2 sin mπ0 πhext λ mπ mπ [ +e 2mπD +e 2mπD ] πl L 2φ K 3 = exp exp K 4 = exp + exp ṁ l p,l /h l D,i +/h sink D,o πl l ṁ l p,l /h l D l,i +/h ext D l,o πl l ṁ l p,l /h l D l,i +/h ext D l,o πhext λ +os mπ mπ + e mπd +e +mπdmπd Bi mπd +Bi πhext λ +os mπ mπ + e mπd +e +mπdmπd Bi mπd +Bi exp +os T ext mπ0 mπ + 0 sin πl l ṁ l p,l /h l D l,i +/h ext D l,o πl L 2φ ṁ l p,l /h l D,i +/h sink D,o mπ0 mπ + 0 sin ] mπ0 mπ + 0 sin ] T sink A.0 A. A.2 A.3 A.4 Appendix B. Seond order eqution for the determintion of T v 0 =Tv 2 [ p,l K 3 K 6 K +K 2 +K 7 K 4 +K 5 K 3 +K 8 +K 8 R /S K 6 K +K 2 K 3 K 4 +K 5 ] +T v [ p,l K 3 K 3 +K 8 K 4 T +K 3 K 4 K K 6 K 8 TR /S +K 6 K T K 7 T +K 9 Q in +h lv p,l K 4 K 6 K +K 2 +K 7 +K 8 K +K 2 K 4 +K 5 K 3 +K 8 +K 3 K 8 R /S +h lv K K +K 2 +K 7 +K 3 K 7 R /S K 3 +K 8 +K 3 K 8 R /S K 0 +h ext S r ] + p,l K 4 K 6 Q in K 9 + T K 6 K +K 7 +K K 8 R /S K 4 K 3 +K 8 +K 3 K 8 R /S 7 B.

19 +h lv T K 3 +K 8 K 4 +K 0 K 6 +K K +K 7 +K 3 K 8 K 4 +K 0 R /S K 3 K 7 K +K K 6 K 8 R /S +h lv K 3 +K 8 +K 3 K 8 R /S Q in K 2 +h ext S r T +T ext K 6 +K +K 3 K R /S Q in K 9 8

THE INFLUENCE OF MODEL RESOLUTION ON AN EXPRESSION OF THE ATMOSPHERIC BOUNDARY LAYER IN A SINGLE-COLUMN MODEL

THE INFLUENCE OF MODEL RESOLUTION ON AN EXPRESSION OF THE ATMOSPHERIC BOUNDARY LAYER IN A SINGLE-COLUMN MODEL THE INFLUENCE OF MODEL RESOLUTION ON AN EXPRESSION OF THE ATMOSPHERIC BOUNDARY LAYER IN A SINGLE-COLUMN MODEL P3.1 Kot Iwmur*, Hiroto Kitgw Jpn Meteorologil Ageny 1. INTRODUCTION Jpn Meteorologil Ageny

More information

Review Topic 14: Relationships between two numerical variables

Review Topic 14: Relationships between two numerical variables Review Topi 14: Reltionships etween two numeril vriles Multiple hoie 1. Whih of the following stterplots est demonstrtes line of est fit? A B C D E 2. The regression line eqution for the following grph

More information

(h+ ) = 0, (3.1) s = s 0, (3.2)

(h+ ) = 0, (3.1) s = s 0, (3.2) Chpter 3 Nozzle Flow Qusistedy idel gs flow in pipes For the lrge vlues of the Reynolds number typilly found in nozzles, the flow is idel. For stedy opertion with negligible body fores the energy nd momentum

More information

ANALYSIS AND MODELLING OF RAINFALL EVENTS

ANALYSIS AND MODELLING OF RAINFALL EVENTS Proeedings of the 14 th Interntionl Conferene on Environmentl Siene nd Tehnology Athens, Greee, 3-5 Septemer 215 ANALYSIS AND MODELLING OF RAINFALL EVENTS IOANNIDIS K., KARAGRIGORIOU A. nd LEKKAS D.F.

More information

Table of Content. c 1 / 5

Table of Content. c 1 / 5 Tehnil Informtion - t nd t Temperture for Controlger 03-2018 en Tble of Content Introdution....................................................................... 2 Definitions for t nd t..............................................................

More information

8 THREE PHASE A.C. CIRCUITS

8 THREE PHASE A.C. CIRCUITS 8 THREE PHSE.. IRUITS The signls in hpter 7 were sinusoidl lternting voltges nd urrents of the so-lled single se type. n emf of suh type n e esily generted y rotting single loop of ondutor (or single winding),

More information

Lecture Notes No. 10

Lecture Notes No. 10 2.6 System Identifition, Estimtion, nd Lerning Leture otes o. Mrh 3, 26 6 Model Struture of Liner ime Invrint Systems 6. Model Struture In representing dynmil system, the first step is to find n pproprite

More information

Appendix C Partial discharges. 1. Relationship Between Measured and Actual Discharge Quantities

Appendix C Partial discharges. 1. Relationship Between Measured and Actual Discharge Quantities Appendi Prtil dishrges. Reltionship Between Mesured nd Atul Dishrge Quntities A dishrging smple my e simply represented y the euilent iruit in Figure. The pplied lternting oltge V is inresed until the

More information

Generalization of 2-Corner Frequency Source Models Used in SMSIM

Generalization of 2-Corner Frequency Source Models Used in SMSIM Generliztion o 2-Corner Frequeny Soure Models Used in SMSIM Dvid M. Boore 26 Mrh 213, orreted Figure 1 nd 2 legends on 5 April 213, dditionl smll orretions on 29 My 213 Mny o the soure spetr models ville

More information

Thermodynamics. Question 1. Question 2. Question 3 3/10/2010. Practice Questions PV TR PV T R

Thermodynamics. Question 1. Question 2. Question 3 3/10/2010. Practice Questions PV TR PV T R /10/010 Question 1 1 mole of idel gs is rought to finl stte F y one of three proesses tht hve different initil sttes s shown in the figure. Wht is true for the temperture hnge etween initil nd finl sttes?

More information

Electromagnetism Notes, NYU Spring 2018

Electromagnetism Notes, NYU Spring 2018 Eletromgnetism Notes, NYU Spring 208 April 2, 208 Ation formultion of EM. Free field desription Let us first onsider the free EM field, i.e. in the bsene of ny hrges or urrents. To tret this s mehnil system

More information

Project 6: Minigoals Towards Simplifying and Rewriting Expressions

Project 6: Minigoals Towards Simplifying and Rewriting Expressions MAT 51 Wldis Projet 6: Minigols Towrds Simplifying nd Rewriting Expressions The distriutive property nd like terms You hve proly lerned in previous lsses out dding like terms ut one prolem with the wy

More information

NEW CIRCUITS OF HIGH-VOLTAGE PULSE GENERATORS WITH INDUCTIVE-CAPACITIVE ENERGY STORAGE

NEW CIRCUITS OF HIGH-VOLTAGE PULSE GENERATORS WITH INDUCTIVE-CAPACITIVE ENERGY STORAGE NEW CIRCUITS OF HIGH-VOLTAGE PULSE GENERATORS WITH INDUCTIVE-CAPACITIVE ENERGY STORAGE V.S. Gordeev, G.A. Myskov Russin Federl Nuler Center All-Russi Sientifi Reserh Institute of Experimentl Physis (RFNC-VNIIEF)

More information

Restraint of purlins for various roof systems

Restraint of purlins for various roof systems NS009 Restrint of purlins for vrious roof systems T. Vrny, M. Brhm & A. Beli ulty of ivil Engineering, zeh Tehnil University, Prh, zehi Astron Buildings S.A., Diekirh, Luxemourg, A memer of the Lind Group

More information

Learning Partially Observable Markov Models from First Passage Times

Learning Partially Observable Markov Models from First Passage Times Lerning Prtilly Oservle Mrkov s from First Pssge s Jérôme Cllut nd Pierre Dupont Europen Conferene on Mhine Lerning (ECML) 8 Septemer 7 Outline. FPT in models nd sequenes. Prtilly Oservle Mrkov s (POMMs).

More information

a) Read over steps (1)- (4) below and sketch the path of the cycle on a P V plot on the graph below. Label all appropriate points.

a) Read over steps (1)- (4) below and sketch the path of the cycle on a P V plot on the graph below. Label all appropriate points. Prole 3: Crnot Cyle of n Idel Gs In this prole, the strting pressure P nd volue of n idel gs in stte, re given he rtio R = / > of the volues of the sttes nd is given Finlly onstnt γ = 5/3 is given You

More information

1 This question is about mean bond enthalpies and their use in the calculation of enthalpy changes.

1 This question is about mean bond enthalpies and their use in the calculation of enthalpy changes. 1 This question is out men ond enthlpies nd their use in the lultion of enthlpy hnges. Define men ond enthlpy s pplied to hlorine. Explin why the enthlpy of tomistion of hlorine is extly hlf the men ond

More information

Novel Fiber-Optical Refractometric Sensor Employing Hemispherically-Shaped Detection Element

Novel Fiber-Optical Refractometric Sensor Employing Hemispherically-Shaped Detection Element Novel Fier-Optil Refrtometri Sensor Employing Hemispherilly-Shped Detetion Element SERGEI KHOTIAINTSEV, VLADIMIR SVIRID Deprtment of Eletril Engineering, Fulty of Engineering Ntionl Autonomous University

More information

PAIR OF LINEAR EQUATIONS IN TWO VARIABLES

PAIR OF LINEAR EQUATIONS IN TWO VARIABLES PAIR OF LINEAR EQUATIONS IN TWO VARIABLES. Two liner equtions in the sme two vriles re lled pir of liner equtions in two vriles. The most generl form of pir of liner equtions is x + y + 0 x + y + 0 where,,,,,,

More information

University of Sioux Falls. MAT204/205 Calculus I/II

University of Sioux Falls. MAT204/205 Calculus I/II University of Sioux Flls MAT204/205 Clulus I/II Conepts ddressed: Clulus Textook: Thoms Clulus, 11 th ed., Weir, Hss, Giordno 1. Use stndrd differentition nd integrtion tehniques. Differentition tehniques

More information

AP CALCULUS Test #6: Unit #6 Basic Integration and Applications

AP CALCULUS Test #6: Unit #6 Basic Integration and Applications AP CALCULUS Test #6: Unit #6 Bsi Integrtion nd Applitions A GRAPHING CALCULATOR IS REQUIRED FOR SOME PROBLEMS OR PARTS OF PROBLEMS IN THIS PART OF THE EXAMINATION. () The ext numeril vlue of the orret

More information

Solutions to Assignment 1

Solutions to Assignment 1 MTHE 237 Fll 2015 Solutions to Assignment 1 Problem 1 Find the order of the differentil eqution: t d3 y dt 3 +t2 y = os(t. Is the differentil eqution liner? Is the eqution homogeneous? b Repet the bove

More information

Section 4.4. Green s Theorem

Section 4.4. Green s Theorem The Clulus of Funtions of Severl Vriles Setion 4.4 Green s Theorem Green s theorem is n exmple from fmily of theorems whih onnet line integrls (nd their higher-dimensionl nlogues) with the definite integrls

More information

SECTION A STUDENT MATERIAL. Part 1. What and Why.?

SECTION A STUDENT MATERIAL. Part 1. What and Why.? SECTION A STUDENT MATERIAL Prt Wht nd Wh.? Student Mteril Prt Prolem n > 0 n > 0 Is the onverse true? Prolem If n is even then n is even. If n is even then n is even. Wht nd Wh? Eploring Pure Mths Are

More information

April 8, 2017 Math 9. Geometry. Solving vector problems. Problem. Prove that if vectors and satisfy, then.

April 8, 2017 Math 9. Geometry. Solving vector problems. Problem. Prove that if vectors and satisfy, then. pril 8, 2017 Mth 9 Geometry Solving vetor prolems Prolem Prove tht if vetors nd stisfy, then Solution 1 onsider the vetor ddition prllelogrm shown in the Figure Sine its digonls hve equl length,, the prllelogrm

More information

Dorf, R.C., Wan, Z. T- Equivalent Networks The Electrical Engineering Handbook Ed. Richard C. Dorf Boca Raton: CRC Press LLC, 2000

Dorf, R.C., Wan, Z. T- Equivalent Networks The Electrical Engineering Handbook Ed. Richard C. Dorf Boca Raton: CRC Press LLC, 2000 orf, R.C., Wn,. T- Equivlent Networks The Eletril Engineering Hndook Ed. Rihrd C. orf Bo Rton: CRC Press LLC, 000 9 T P Equivlent Networks hen Wn University of Cliforni, vis Rihrd C. orf University of

More information

Math 32B Discussion Session Week 8 Notes February 28 and March 2, f(b) f(a) = f (t)dt (1)

Math 32B Discussion Session Week 8 Notes February 28 and March 2, f(b) f(a) = f (t)dt (1) Green s Theorem Mth 3B isussion Session Week 8 Notes Februry 8 nd Mrh, 7 Very shortly fter you lerned how to integrte single-vrible funtions, you lerned the Fundmentl Theorem of lulus the wy most integrtion

More information

TOPIC: LINEAR ALGEBRA MATRICES

TOPIC: LINEAR ALGEBRA MATRICES Interntionl Blurete LECTUE NOTES for FUTHE MATHEMATICS Dr TOPIC: LINEA ALGEBA MATICES. DEFINITION OF A MATIX MATIX OPEATIONS.. THE DETEMINANT deta THE INVESE A -... SYSTEMS OF LINEA EQUATIONS. 8. THE AUGMENTED

More information

Comparing the Pre-image and Image of a Dilation

Comparing the Pre-image and Image of a Dilation hpter Summry Key Terms Postultes nd Theorems similr tringles (.1) inluded ngle (.2) inluded side (.2) geometri men (.) indiret mesurement (.6) ngle-ngle Similrity Theorem (.2) Side-Side-Side Similrity

More information

for all x in [a,b], then the area of the region bounded by the graphs of f and g and the vertical lines x = a and x = b is b [ ( ) ( )] A= f x g x dx

for all x in [a,b], then the area of the region bounded by the graphs of f and g and the vertical lines x = a and x = b is b [ ( ) ( )] A= f x g x dx Applitions of Integrtion Are of Region Between Two Curves Ojetive: Fin the re of region etween two urves using integrtion. Fin the re of region etween interseting urves using integrtion. Desrie integrtion

More information

Estimation of Global Solar Radiation in Onitsha and Calabar Using Empirical Models

Estimation of Global Solar Radiation in Onitsha and Calabar Using Empirical Models Communitions in Applied Sienes ISS 0-77 Volume, umer, 0, 5-7 Estimtion of Glol Solr dition in Onitsh nd Clr Using Empiril Models M.. nuhi, J. E. Ekpe nd G. F Ieh Deprtment of Industril Physis, Eonyi Stte

More information

On the Scale factor of the Universe and Redshift.

On the Scale factor of the Universe and Redshift. On the Sle ftor of the Universe nd Redshift. J. M. unter. john@grvity.uk.om ABSTRACT It is proposed tht there hs been longstnding misunderstnding of the reltionship between sle ftor of the universe nd

More information

Lecture 1 - Introduction and Basic Facts about PDEs

Lecture 1 - Introduction and Basic Facts about PDEs * 18.15 - Introdution to PDEs, Fll 004 Prof. Gigliol Stffilni Leture 1 - Introdution nd Bsi Fts bout PDEs The Content of the Course Definition of Prtil Differentil Eqution (PDE) Liner PDEs VVVVVVVVVVVVVVVVVVVV

More information

Matrices SCHOOL OF ENGINEERING & BUILT ENVIRONMENT. Mathematics (c) 1. Definition of a Matrix

Matrices SCHOOL OF ENGINEERING & BUILT ENVIRONMENT. Mathematics (c) 1. Definition of a Matrix tries Definition of tri mtri is regulr rry of numers enlosed inside rkets SCHOOL OF ENGINEERING & UIL ENVIRONEN Emple he following re ll mtries: ), ) 9, themtis ), d) tries Definition of tri Size of tri

More information

A Study on the Properties of Rational Triangles

A Study on the Properties of Rational Triangles Interntionl Journl of Mthemtis Reserh. ISSN 0976-5840 Volume 6, Numer (04), pp. 8-9 Interntionl Reserh Pulition House http://www.irphouse.om Study on the Properties of Rtionl Tringles M. Q. lm, M.R. Hssn

More information

Lecture Summaries for Multivariable Integral Calculus M52B

Lecture Summaries for Multivariable Integral Calculus M52B These leture summries my lso be viewed online by liking the L ion t the top right of ny leture sreen. Leture Summries for Multivrible Integrl Clulus M52B Chpter nd setion numbers refer to the 6th edition.

More information

I1 = I2 I1 = I2 + I3 I1 + I2 = I3 + I4 I 3

I1 = I2 I1 = I2 + I3 I1 + I2 = I3 + I4 I 3 2 The Prllel Circuit Electric Circuits: Figure 2- elow show ttery nd multiple resistors rrnged in prllel. Ech resistor receives portion of the current from the ttery sed on its resistnce. The split is

More information

SECOND HARMONIC GENERATION OF Bi 4 Ti 3 O 12 FILMS

SECOND HARMONIC GENERATION OF Bi 4 Ti 3 O 12 FILMS SECOND HARMONIC GENERATION OF Bi 4 Ti 3 O 12 FILMS IN-SITU PROBING OF DOMAIN POLING IN Bi 4 Ti 3 O 12 THIN FILMS BY OPTICAL SECOND HARMONIC GENERATION YANIV BARAD, VENKATRAMAN GOPALAN Mterils Reserh Lortory

More information

1.3 SCALARS AND VECTORS

1.3 SCALARS AND VECTORS Bridge Course Phy I PUC 24 1.3 SCLRS ND VECTORS Introdution: Physis is the study of nturl phenomen. The study of ny nturl phenomenon involves mesurements. For exmple, the distne etween the plnet erth nd

More information

Application of the theory of compound cores for the assessment of stress pattern in the cross section of a strengthened beam column

Application of the theory of compound cores for the assessment of stress pattern in the cross section of a strengthened beam column IOP Conferene Series: Mterils Siene nd Engineering PAPER OPEN ACCESS Applition of the theory of ompound ores for the ssessment of stress pttern in the ross setion of strengthened bem olumn To ite this

More information

New Expansion and Infinite Series

New Expansion and Infinite Series Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University

More information

THE PYTHAGOREAN THEOREM

THE PYTHAGOREAN THEOREM THE PYTHAGOREAN THEOREM The Pythgoren Theorem is one of the most well-known nd widely used theorems in mthemtis. We will first look t n informl investigtion of the Pythgoren Theorem, nd then pply this

More information

SIDESWAY MAGNIFICATION FACTORS FOR STEEL MOMENT FRAMES WITH VARIOUS TYPES OF COLUMN BASES

SIDESWAY MAGNIFICATION FACTORS FOR STEEL MOMENT FRAMES WITH VARIOUS TYPES OF COLUMN BASES Advned Steel Constrution Vol., No., pp. 7-88 () 7 SIDESWAY MAGNIFICATION FACTORS FOR STEEL MOMENT FRAMES WIT VARIOUS TYPES OF COLUMN BASES J. ent sio Assoite Professor, Deprtment of Civil nd Environmentl

More information

Numbers and indices. 1.1 Fractions. GCSE C Example 1. Handy hint. Key point

Numbers and indices. 1.1 Fractions. GCSE C Example 1. Handy hint. Key point GCSE C Emple 7 Work out 9 Give your nswer in its simplest form Numers n inies Reiprote mens invert or turn upsie own The reiprol of is 9 9 Mke sure you only invert the frtion you re iviing y 7 You multiply

More information

THE ASYMMETRY OF COASTAL WATER LEVEL RESPONSE TO LANDFALLING HURRICANES SIMULATED BY A THREE-DIMENSIONAL STORM SURGE MODEL

THE ASYMMETRY OF COASTAL WATER LEVEL RESPONSE TO LANDFALLING HURRICANES SIMULATED BY A THREE-DIMENSIONAL STORM SURGE MODEL THE ASYMMETRY OF COASTAL WATER LEVEL RESPONSE TO LANDFALLING HURRICANES SIMULATED BY A THREE-DIMENSIONAL STORM SURGE MODEL Mhun Peng *, Lin Xie nd Leonrd J. Pietrfes Deprtment of Mrine, Erth nd Atmospheri

More information

Lecture 6: Coding theory

Lecture 6: Coding theory Leture 6: Coing theory Biology 429 Crl Bergstrom Ferury 4, 2008 Soures: This leture loosely follows Cover n Thoms Chpter 5 n Yeung Chpter 3. As usul, some of the text n equtions re tken iretly from those

More information

Review of Gaussian Quadrature method

Review of Gaussian Quadrature method Review of Gussin Qudrture method Nsser M. Asi Spring 006 compiled on Sundy Decemer 1, 017 t 09:1 PM 1 The prolem To find numericl vlue for the integrl of rel vlued function of rel vrile over specific rnge

More information

Engr354: Digital Logic Circuits

Engr354: Digital Logic Circuits Engr354: Digitl Logi Ciruits Chpter 4: Logi Optimiztion Curtis Nelson Logi Optimiztion In hpter 4 you will lern out: Synthesis of logi funtions; Anlysis of logi iruits; Tehniques for deriving minimum-ost

More information

1 This diagram represents the energy change that occurs when a d electron in a transition metal ion is excited by visible light.

1 This diagram represents the energy change that occurs when a d electron in a transition metal ion is excited by visible light. 1 This igrm represents the energy hnge tht ours when eletron in trnsition metl ion is exite y visile light. Give the eqution tht reltes the energy hnge ΔE to the Plnk onstnt, h, n the frequeny, v, of the

More information

Finite Element Simulation on Frictional and Brittle Preseismic fault slip

Finite Element Simulation on Frictional and Brittle Preseismic fault slip Finite Element Simultion on Fritionl nd Brittle Preseismi fult slip Zhishen Wu (1) Yun Go (1) Yutk Murkmi (2) (1) Deprtment of Urn & Civil Engineering. Irki University, Jpn (e-mil: zswu@ip.irki..jp; goyun@hs.irki..jp,

More information

Consequently, the temperature must be the same at each point in the cross section at x. Let:

Consequently, the temperature must be the same at each point in the cross section at x. Let: HW 2 Comments: L1-3. Derive the het eqution for n inhomogeneous rod where the therml coefficients used in the derivtion of the het eqution for homogeneous rod now become functions of position x in the

More information

Forces on curved surfaces Buoyant force Stability of floating and submerged bodies

Forces on curved surfaces Buoyant force Stability of floating and submerged bodies Stti Surfe ores Stti Surfe ores 8m wter hinge? 4 m ores on plne res ores on urved surfes Buont fore Stbilit of floting nd submerged bodies ores on Plne res Two tpes of problems Horizontl surfes (pressure

More information

Activities. 4.1 Pythagoras' Theorem 4.2 Spirals 4.3 Clinometers 4.4 Radar 4.5 Posting Parcels 4.6 Interlocking Pipes 4.7 Sine Rule Notes and Solutions

Activities. 4.1 Pythagoras' Theorem 4.2 Spirals 4.3 Clinometers 4.4 Radar 4.5 Posting Parcels 4.6 Interlocking Pipes 4.7 Sine Rule Notes and Solutions MEP: Demonstrtion Projet UNIT 4: Trigonometry UNIT 4 Trigonometry tivities tivities 4. Pythgors' Theorem 4.2 Spirls 4.3 linometers 4.4 Rdr 4.5 Posting Prels 4.6 Interloking Pipes 4.7 Sine Rule Notes nd

More information

Green s Theorem. (2x e y ) da. (2x e y ) dx dy. x 2 xe y. (1 e y ) dy. y=1. = y e y. y=0. = 2 e

Green s Theorem. (2x e y ) da. (2x e y ) dx dy. x 2 xe y. (1 e y ) dy. y=1. = y e y. y=0. = 2 e Green s Theorem. Let be the boundry of the unit squre, y, oriented ounterlokwise, nd let F be the vetor field F, y e y +, 2 y. Find F d r. Solution. Let s write P, y e y + nd Q, y 2 y, so tht F P, Q. Let

More information

Learning Objectives of Module 2 (Algebra and Calculus) Notes:

Learning Objectives of Module 2 (Algebra and Calculus) Notes: 67 Lerning Ojetives of Module (Alger nd Clulus) Notes:. Lerning units re grouped under three res ( Foundtion Knowledge, Alger nd Clulus ) nd Further Lerning Unit.. Relted lerning ojetives re grouped under

More information

1 PYTHAGORAS THEOREM 1. Given a right angled triangle, the square of the hypotenuse is equal to the sum of the squares of the other two sides.

1 PYTHAGORAS THEOREM 1. Given a right angled triangle, the square of the hypotenuse is equal to the sum of the squares of the other two sides. 1 PYTHAGORAS THEOREM 1 1 Pythgors Theorem In this setion we will present geometri proof of the fmous theorem of Pythgors. Given right ngled tringle, the squre of the hypotenuse is equl to the sum of the

More information

Section 1.3 Triangles

Section 1.3 Triangles Se 1.3 Tringles 21 Setion 1.3 Tringles LELING TRINGLE The line segments tht form tringle re lled the sides of the tringle. Eh pir of sides forms n ngle, lled n interior ngle, nd eh tringle hs three interior

More information

6.5 Improper integrals

6.5 Improper integrals Eerpt from "Clulus" 3 AoPS In. www.rtofprolemsolving.om 6.5. IMPROPER INTEGRALS 6.5 Improper integrls As we ve seen, we use the definite integrl R f to ompute the re of the region under the grph of y =

More information

MATH Final Review

MATH Final Review MATH 1591 - Finl Review November 20, 2005 1 Evlution of Limits 1. the ε δ definition of limit. 2. properties of limits. 3. how to use the diret substitution to find limit. 4. how to use the dividing out

More information

Thermal energy 2 U Q W. 23 April The First Law of Thermodynamics. Or, if we want to obtain external work: The trick of using steam

Thermal energy 2 U Q W. 23 April The First Law of Thermodynamics. Or, if we want to obtain external work: The trick of using steam April 08 Therml energy Soures of het Trnsport of het How to use het The First Lw of Thermoynmis U W Or, if we wnt to otin externl work: U W 009 vrije Universiteit msterm Close yle stem power plnt The trik

More information

GM1 Consolidation Worksheet

GM1 Consolidation Worksheet Cmridge Essentils Mthemtis Core 8 GM1 Consolidtion Worksheet GM1 Consolidtion Worksheet 1 Clulte the size of eh ngle mrked y letter. Give resons for your nswers. or exmple, ngles on stright line dd up

More information

Cyclic voltammetry simulation at microelectrode arrays with COMSOL Multiphysics Ò

Cyclic voltammetry simulation at microelectrode arrays with COMSOL Multiphysics Ò J Appl Eletrohem (009) 39:9 63 DOI 0.007/s0800-009-9797- ORIGINAL PAPER Cyli voltmmetry simultion t miroeletrode rrys with COMSOL Multiphysis Ò Alessndro Lvhi Æ U. Brdi Æ C. Borri Æ S. Cporli Æ A. Fossti

More information

NON-DETERMINISTIC FSA

NON-DETERMINISTIC FSA Tw o types of non-determinism: NON-DETERMINISTIC FS () Multiple strt-sttes; strt-sttes S Q. The lnguge L(M) ={x:x tkes M from some strt-stte to some finl-stte nd ll of x is proessed}. The string x = is

More information

Linear Algebra Introduction

Linear Algebra Introduction Introdution Wht is Liner Alger out? Liner Alger is rnh of mthemtis whih emerged yers k nd ws one of the pioneer rnhes of mthemtis Though, initilly it strted with solving of the simple liner eqution x +

More information

AP Calculus BC Chapter 8: Integration Techniques, L Hopital s Rule and Improper Integrals

AP Calculus BC Chapter 8: Integration Techniques, L Hopital s Rule and Improper Integrals AP Clulus BC Chpter 8: Integrtion Tehniques, L Hopitl s Rule nd Improper Integrls 8. Bsi Integrtion Rules In this setion we will review vrious integrtion strtegies. Strtegies: I. Seprte the integrnd into

More information

Vapour compression refrigeration system

Vapour compression refrigeration system www.ookspr.om VTU NEWS VTU NOTES QUESTION PAPERS FORUMS RESULTS Vpour ompression refrigertion system Introdution In vpour ompression system, the refrigernts used re mmoni, ron dioxide, freons et. the refrigernts

More information

A Matlab/Simulink Model of a Langevin s Ultrasonic Power Transducers

A Matlab/Simulink Model of a Langevin s Ultrasonic Power Transducers Mtl/Simulink Model of Lngevin s Ultrsoni Power Trnsduers Igor Jovnović, Ugleš Jovnović nd Drgn Mnčić strt Ultrsoni sndwih trnsduer, lso known s Lngevin s trnsduer, is hlf-wve resonnt struture tht osilltes

More information

Introduction to Olympiad Inequalities

Introduction to Olympiad Inequalities Introdution to Olympid Inequlities Edutionl Studies Progrm HSSP Msshusetts Institute of Tehnology Snj Simonovikj Spring 207 Contents Wrm up nd Am-Gm inequlity 2. Elementry inequlities......................

More information

Switching properties of an arbitrarily excited nonlinear electron-wave directional coupler

Switching properties of an arbitrarily excited nonlinear electron-wave directional coupler Proeedings of the 6th WSEAS Interntionl Conferene on Miroeletronis, Nnoeletronis, Optoeletronis, Istnul, Turkey, My 7-9, 7 1 Swithing properties of n ritrrily exited nonliner eletron-wve diretionl oupler

More information

ANALYSIS OF CFD HEAT TRANSFER OF VACUUM FREEZE-DRYING SHELF

ANALYSIS OF CFD HEAT TRANSFER OF VACUUM FREEZE-DRYING SHELF HEFAT202 9 th Interntionl Conferene on Het Trnsfer, Fluid Mehnis nd Thermodynmis 6 8 July 202 Mlt ANALYSIS OF CFD HEAT TRANSFER OF VACUUM FREEZE-DRYING SHELF Hong-Ping Cheng *, Shin-min Tsi 2 * Professor,

More information

Psychrometric Applications

Psychrometric Applications Psychrometric Applictions The reminder of this presenttion centers on systems involving moist ir. A condensed wter phse my lso be present in such systems. The term moist irrefers to mixture of dry ir nd

More information

AP Calculus AB Unit 4 Assessment

AP Calculus AB Unit 4 Assessment Clss: Dte: 0-04 AP Clulus AB Unit 4 Assessment Multiple Choie Identify the hoie tht best ompletes the sttement or nswers the question. A lultor my NOT be used on this prt of the exm. (6 minutes). The slope

More information

CHAPTER 20: Second Law of Thermodynamics

CHAPTER 20: Second Law of Thermodynamics CHAER 0: Second Lw of hermodynmics Responses to Questions 3. kg of liquid iron will hve greter entropy, since it is less ordered thn solid iron nd its molecules hve more therml motion. In ddition, het

More information

MATRIX INVERSE ON CONNEX PARALLEL ARCHITECTURE

MATRIX INVERSE ON CONNEX PARALLEL ARCHITECTURE U.P.B. Si. Bull., Series C, Vol. 75, Iss. 2, ISSN 86 354 MATRIX INVERSE ON CONNEX PARALLEL ARCHITECTURE An-Mri CALFA, Gheorghe ŞTEFAN 2 Designed for emedded omputtion in system on hip design, the Connex

More information

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1. 1 [(y ) 2 + yy + y 2 ] dx,

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1. 1 [(y ) 2 + yy + y 2 ] dx, MATH3403: Green s Funtions, Integrl Equtions nd the Clulus of Vritions 1 Exmples 5 Qu.1 Show tht the extreml funtion of the funtionl I[y] = 1 0 [(y ) + yy + y ] dx, where y(0) = 0 nd y(1) = 1, is y(x)

More information

Some Aspects of Non-Orthogonal Stagnation-Point Flow towards a Stretching Surface

Some Aspects of Non-Orthogonal Stagnation-Point Flow towards a Stretching Surface Engineering, 00,, 705-709 doi:0.436/eng.00.909 Published Online September 00 (http://www.sirp.org/journl/eng) Some Aspets of Non-Orthogonl Stgntion-Point Flow towrds Strething Surfe Abstrt Mothr Rez, Andi

More information

THE ANALYSIS AND CALCULATION OF ELECTROMAGNETIC FIELD AROUND OVERHEAD POWER LINE HongWang Yang

THE ANALYSIS AND CALCULATION OF ELECTROMAGNETIC FIELD AROUND OVERHEAD POWER LINE HongWang Yang 5th Interntionl Conferene on Advned Mterils nd Computer Siene (ICAMCS 6) THE ANALYSIS AN CALCULATION OF ELECTROMAGNETIC FIEL AROUN OVERHEA POWER LINE HongWng Yng eprtment of eletril engineering, North

More information

HYDROTHERMAL ANALYSIS OF THE ABSORBER TUBES USED IN LINEAR FRESNEL REFLECTOR SOLAR THERMAL SYSTEM

HYDROTHERMAL ANALYSIS OF THE ABSORBER TUBES USED IN LINEAR FRESNEL REFLECTOR SOLAR THERMAL SYSTEM Proceedings of the 21 st Ntionl & 10 th ISHMT-ASME Het nd Mss Trnsfer Conference Decemer 27-30, 2011, IIT Mdrs, Indi Pper ID: ISHMT_IND_17_034 HYDROTHERMAL ANALYSIS OF THE ABSORBER TUBES USED IN LINEAR

More information

Part 4. Integration (with Proofs)

Part 4. Integration (with Proofs) Prt 4. Integrtion (with Proofs) 4.1 Definition Definition A prtition P of [, b] is finite set of points {x 0, x 1,..., x n } with = x 0 < x 1

More information

Linear choosability of graphs

Linear choosability of graphs Liner hoosility of grphs Louis Esperet, Mikel Montssier, André Rspud To ite this version: Louis Esperet, Mikel Montssier, André Rspud. Liner hoosility of grphs. Stefn Felsner. 2005 Europen Conferene on

More information

A Lower Bound for the Length of a Partial Transversal in a Latin Square, Revised Version

A Lower Bound for the Length of a Partial Transversal in a Latin Square, Revised Version A Lower Bound for the Length of Prtil Trnsversl in Ltin Squre, Revised Version Pooy Htmi nd Peter W. Shor Deprtment of Mthemtil Sienes, Shrif University of Tehnology, P.O.Bo 11365-9415, Tehrn, Irn Deprtment

More information

(a) A partition P of [a, b] is a finite subset of [a, b] containing a and b. If Q is another partition and P Q, then Q is a refinement of P.

(a) A partition P of [a, b] is a finite subset of [a, b] containing a and b. If Q is another partition and P Q, then Q is a refinement of P. Chpter 7: The Riemnn Integrl When the derivtive is introdued, it is not hrd to see tht the it of the differene quotient should be equl to the slope of the tngent line, or when the horizontl xis is time

More information

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies Stte spce systems nlysis (continued) Stbility A. Definitions A system is sid to be Asymptoticlly Stble (AS) when it stisfies ut () = 0, t > 0 lim xt () 0. t A system is AS if nd only if the impulse response

More information

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION XX IMEKO World Congress Metrology for Green Growth September 9,, Busn, Republic of Kore THERMAL EXPANSION COEFFICIENT OF WATER FOR OLUMETRIC CALIBRATION Nieves Medin Hed of Mss Division, CEM, Spin, mnmedin@mityc.es

More information

Section 3.6. Definite Integrals

Section 3.6. Definite Integrals The Clulus of Funtions of Severl Vribles Setion.6 efinite Integrls We will first define the definite integrl for funtion f : R R nd lter indite how the definition my be extended to funtions of three or

More information

VISIBLE AND INFRARED ABSORPTION SPECTRA OF COVERING MATERIALS FOR SOLAR COLLECTORS

VISIBLE AND INFRARED ABSORPTION SPECTRA OF COVERING MATERIALS FOR SOLAR COLLECTORS AGRICULTURAL ENGINEERING VISIBLE AND INFRARED ABSORPTION SPECTRA OF COVERING MATERIALS FOR SOLAR COLLECTORS Ltvi University of Agriulture E-mil: ilze.pelee@llu.lv Astrt Use of solr energy inreses every

More information

Calculus Cheat Sheet. Integrals Definitions. where F( x ) is an anti-derivative of f ( x ). Fundamental Theorem of Calculus. dx = f x dx g x dx

Calculus Cheat Sheet. Integrals Definitions. where F( x ) is an anti-derivative of f ( x ). Fundamental Theorem of Calculus. dx = f x dx g x dx Clulus Chet Sheet Integrls Definitions Definite Integrl: Suppose f ( ) is ontinuous Anti-Derivtive : An nti-derivtive of f ( ) on [, ]. Divide [, ] into n suintervls of is funtion, F( ), suh tht F = f.

More information

Lecture 6. CMOS Static & Dynamic Logic Gates. Static CMOS Circuit. PMOS Transistors in Series/Parallel Connection

Lecture 6. CMOS Static & Dynamic Logic Gates. Static CMOS Circuit. PMOS Transistors in Series/Parallel Connection NMOS Trnsistors in Series/Prllel onnetion Leture 6 MOS Stti & ynmi Logi Gtes Trnsistors n e thought s swith ontrolled y its gte signl NMOS swith loses when swith ontrol input is high Peter heung eprtment

More information

Parabola and Catenary Equations for Conductor Height Calculation

Parabola and Catenary Equations for Conductor Height Calculation ELECTROTEHNICĂ, ELECTRONICĂ, AUTOMATICĂ, 6 (), nr. 3 9 Prbol nd Ctenr Equtions for Condutor Height Clultion Alen HATIBOVIC Abstrt This pper presents new equtions for ondutor height lultion bsed on the

More information

Lecture 27: Diffusion of Ions: Part 2: coupled diffusion of cations and

Lecture 27: Diffusion of Ions: Part 2: coupled diffusion of cations and Leture 7: iffusion of Ions: Prt : oupled diffusion of tions nd nions s desried y Nernst-Plnk Eqution Tody s topis Continue to understnd the fundmentl kinetis prmeters of diffusion of ions within n eletrilly

More information

A Mathematical Model for Unemployment-Taking an Action without Delay

A Mathematical Model for Unemployment-Taking an Action without Delay Advnes in Dynmil Systems nd Applitions. ISSN 973-53 Volume Number (7) pp. -8 Reserh Indi Publitions http://www.ripublition.om A Mthemtil Model for Unemployment-Tking n Ation without Dely Gulbnu Pthn Diretorte

More information

Chemical Equilibrium

Chemical Equilibrium Chpter 16 Questions 5, 7, 31, 33, 35, 43, 71 Chemil Equilibrium Exmples of Equilibrium Wter n exist simultneously in the gs nd liquid phse. The vpor pressure of H O t given temperture is property ssoited

More information

Distributed Generation Placement in Unbalanced Distribution System with Seasonal Load Variation

Distributed Generation Placement in Unbalanced Distribution System with Seasonal Load Variation Distriuted Genertion Plement in Unlned Distriution System with Sesonl Lod Vrition Rvi Tej Bhimrsetti Dept. of Eletril Engg., NT Kurukshetr Kurukshetr, ndi svrtej@gmil.om Ashwni Kumr, Memer, EEE Dept. of

More information

6.3.2 Spectroscopy. N Goalby chemrevise.org 1 NO 2 CH 3. CH 3 C a. NMR spectroscopy. Different types of NMR

6.3.2 Spectroscopy. N Goalby chemrevise.org 1 NO 2 CH 3. CH 3 C a. NMR spectroscopy. Different types of NMR 6.. Spetrosopy NMR spetrosopy Different types of NMR NMR spetrosopy involves intertion of mterils with the lowenergy rdiowve region of the eletromgneti spetrum NMR spetrosopy is the sme tehnology s tht

More information

Polyphase Systems. Objectives 23.1 INTRODUCTION

Polyphase Systems. Objectives 23.1 INTRODUCTION Polyphse Systems 23 Ojetives eome fmilir with the opertion of threephse genertor nd the mgnitude nd phse reltionship etween the three phse voltges. e le to lulte the voltges nd urrents for three-phse Y-onneted

More information

Now we must transform the original model so we can use the new parameters. = S max. Recruits

Now we must transform the original model so we can use the new parameters. = S max. Recruits MODEL FOR VARIABLE RECRUITMENT (ontinue) Alterntive Prmeteriztions of the pwner-reruit Moels We n write ny moel in numerous ifferent ut equivlent forms. Uner ertin irumstnes it is onvenient to work with

More information

7.1 Integral as Net Change and 7.2 Areas in the Plane Calculus

7.1 Integral as Net Change and 7.2 Areas in the Plane Calculus 7.1 Integrl s Net Chnge nd 7. Ares in the Plne Clculus 7.1 INTEGRAL AS NET CHANGE Notecrds from 7.1: Displcement vs Totl Distnce, Integrl s Net Chnge We hve lredy seen how the position of n oject cn e

More information

QUADRATIC EQUATION. Contents

QUADRATIC EQUATION. Contents QUADRATIC EQUATION Contents Topi Pge No. Theory 0-04 Exerise - 05-09 Exerise - 09-3 Exerise - 3 4-5 Exerise - 4 6 Answer Key 7-8 Syllus Qudrti equtions with rel oeffiients, reltions etween roots nd oeffiients,

More information

A Brief Note on Quasi Static Thermal Stresses In A Thin Rectangular Plate With Internal Heat Generation

A Brief Note on Quasi Static Thermal Stresses In A Thin Rectangular Plate With Internal Heat Generation Americn Journl of Engineering Reserch (AJER) 13 Americn Journl of Engineering Reserch (AJER) e-issn : 3-847 p-issn : 3-936 Volume-, Issue-1, pp-388-393 www.jer.org Reserch Pper Open Access A Brief Note

More information

Magnetically Coupled Coil

Magnetically Coupled Coil Mgnetilly Coupled Ciruits Overview Mutul Indutne Energy in Coupled Coils Liner Trnsformers Idel Trnsformers Portlnd Stte University ECE 22 Mgnetilly Coupled Ciruits Ver..3 Mgnetilly Coupled Coil i v L

More information