Nanoscale Heat Transfer using Phonon Boltzmann Transport Equation
|
|
- Franklin Phelps
- 5 years ago
- Views:
Transcription
1 Excerpt frm the Prceedngs f the COMSOL Cnference 9 Bstn Nanscale Heat Transfer usng Phnn Bltzmann Transprt Equatn Sangwk Shn *, and Ajt K. Ry Ar Frce Research Labratry, Unversty f Daytn Research nsttute *3 Cllege Park, Daytn, OH , sangwk@stanfrdalumn.rg Abstract: COMSOL Multphyscs was used t slve a phnn Bltzmann transprt equatn (BTE) fr nanscale heat transprt prblems. One-dmensnal steady-state and transent BTE prblems were successfully slved based n fnte element and dscrete rdnate methds fr spatal and angular dscretzatns, respectvely, by utlzng the bult-n feature f the COMSOL, Ceffcent Frm f PDE. A senstvty study was cnducted wth varus dscretzatn refnements fr dfferent values f the Knudsen number, whch s a measure f the nanscale regme. t was fund that suffcent refnement fr angular dscretzatn s crtcal n btanng accurate slutns f the BTE. Keywrds: Nanscale heat transprt, phnn Bltzmann transprt equatn, fnte element methd, dscrete rdnates methd.. ntrductn Fr the last tw centures, many heat cnductn prblems have been mdeled successfully by a Furer dffusve equatn (FDE). The equatn can be derved usng a cnservatn law f energy and Furer s lnear apprxmatn f heat flux usng a temperature gradent. The FDE s a parablc equatn reflectng a dffusve nature f heat transprt. An underlyng assumptn s that the heat s effectvely transferred between lcalzed regns thrugh suffcent scatterng events f phnns wthn a medum. Therefre, the FDE des nt hld when the number f scatterngs s neglgble, whch culd happen when a mean free path f phnn s smlar r a larger rder f magntude than a dman sze f nterest, sgnfcant bundary scatterng ccurs at materal nterfaces causng thermal resstance, etc. Anther prblem f the FDE s that t admts an nfnte speed f the heat transprt, whch s cntradctry t the thery f relatvty. Therefre, the FDE s napprprate fr the heat transfer prblems at small tme and spatal scales. T reslve the ssue f the nfnte speed f heat carrer, hyperblc equatnns, such as a Cattane equatn (Telegraph equatn) [] r a relatvstc heat equatn [], are ften used t reflect a wave nature f the heat transprt fr cases f small speed f heat prpagatn (e.g., lw temperature, lw cnductve medum, thermal nsulatr, phase transtn, etc.) and/r hgh speed f heat flux (e.g., pc- r femtsecnd pulsed-laser heatng.) Despte sme ssues [], these hyperblc equatns can be used t cnsder the fnte speed f phnns fr a shrt tme scale. Hwever, they stll cannt be used fr small spatal scale. Fr the nanscale heat transfer analyss, ne needs equatns and methds fr small-scale smulatn n terms f bth tme and space. A mlecular dynamcs (MD) smulatn can be a useful and accurate methd n ths regard. Hwever, the MD smulatn s usually cmputatnally expensve, s t s sutable fr systems havng a few atmc layers r several thusand atms, but nt sutable fr devce-level thermal analyss. Alternatvely, a phnn Bltzmann transprt equatn (BTE), a frstrder partal dfferental equatn fr a phnn dstrbutn functn, has been used fr ths purpse [3-5]. The dstrbutn functn s a scalar quantty n the sx-dmensnal phase space (three space crdnates and three wavevectr crdnates). The phnn BTE s als called an equatn f phnn radatve transfer (EPRT) when the phnn dstrbutn functn s replaced wth a phnn ntensty functn [3]. The BTE s knwn t be dffcult t slve, and thus ften smplfed wth a relaxatn tme apprxmatn. The BTE can predct a ballstc nature f heat transfer under an assumptn that partcle-lke behavrs f phnn are much mre sgnfcant than ts wave-lke behavrs, whch make the BTE vald fr structures larger than the wavelength f phnns. The BTE can be slved analytcally fr smple gemetres [4, 6, 7], r numercally fr cmplex gemetres by usng ether determnstc methds (e.g., dscrete rdnates methd, sphercal harmncs methd,
2 Reprt Dcumentatn Page Frm Apprved OMB N Publc reprtng burden fr the cllectn f nfrmatn s estmated t average hur per respnse, ncludng the tme fr revewng nstructns, searchng exstng data surces, gatherng and mantanng the data needed, and cmpletng and revewng the cllectn f nfrmatn. Send cmments regardng ths burden estmate r any ther aspect f ths cllectn f nfrmatn, ncludng suggestns fr reducng ths burden, t Washngtn Headquarters Servces, Drectrate fr nfrmatn Operatns and Reprts, 5 Jeffersn Davs Hghway, Sute 4, Arlngtn VA -43. Respndents shuld be aware that ntwthstandng any ther prvsn f law, n persn shall be subject t a penalty fr falng t cmply wth a cllectn f nfrmatn f t des nt dsplay a currently vald OMB cntrl number.. REPORT DATE OCT 9. REPORT TYPE 3. DATES COVERED --9 t TTLE AND SUBTTLE Nanscale Heat Transfer usng Phnn Bltzmann Transprt Equatn 5a. CONTRACT NUMBER 5b. GRANT NUMBER 5c. PROGRAM ELEMENT NUMBER 6. AUTHOR(S) 5d. PROJECT NUMBER 5e. TASK NUMBER 5f. WORK UNT NUMBER 7. PERFORMNG ORGANZATON NAME(S) AND ADDRESS(ES) Ar Frce Research Labratry,3 Cllege Park,Daytn,OH, PERFORMNG ORGANZATON REPORT NUMBER 9. SPONSORNG/MONTORNG AGENCY NAME(S) AND ADDRESS(ES). SPONSOR/MONTOR S ACRONYM(S). DSTRBUTON/AVALABLTY STATEMENT Apprved fr publc release; dstrbutn unlmted 3. SUPPLEMENTARY NOTES Presented at the 9 COMSOL Cnference Bstn, 8- Oct, Bstn, MA. SPONSOR/MONTOR S REPORT NUMBER(S) 4. ABSTRACT COMSOL Multphyscs was used t slve a phnn Bltzmann transprt equatn (BTE) fr nanscale heat transprt prblems. One-dmensnal steady-state and transent BTE prblems were successfully slved based n fnte element and dscrete rdnate methds fr spatal and angular dscretzatns, respectvely by utlzng the bult-n feature f the COMSOL Ceffcent Frm f PDE. A senstvty study was cnducted wth varus dscretzatn refnements fr dfferent values f the Knudsen number, whch s a measure f the nanscale regme. t was fund that suffcent refnement fr angular dscretzatn s crtcal n btanng accurate slutns f the BTE. 5. SUBJECT TERMS 6. SECURTY CLASSFCATON OF: 7. LMTATON OF ABSTRACT a. REPORT unclassfed b. ABSTRACT unclassfed c. THS PAGE unclassfed Same as Reprt (SAR) 8. NUMBER OF PAGES 3 9a. NAME OF RESPONSBLE PERSON Standard Frm 98 (Rev. 8-98) Prescrbed by ANS Std Z39-8
3 fnte vlume methd [8], etc.) r statstcal methds (e.g., Mnte Carl smulatn [9, ]). n general, slvng the BTE s much mre effcent than the MD, and the predctns agree well wth expermental data [, ]. n the present study, we wll slve the BTE wth the relaxatn tme apprxmatn fr nedmensnal heat transprt prblems usng the COMSOL Multphyscs. We wll cmbne the fnte element (FE) capablty f the COMSOL wth the DOM t slve bth steady-state and transent heat transfer prblems. The numercal results wll be cmpared wth analytcal nes fr the steady-state prblem. A senstvty study wll be cnducted fr varus scale measures t shw the dfference between the BTE slutns wth the cnventnal FDE nes. Furthermre, we wll determne a sutable dscretzatn scheme fr spatal and angular spaces n rder t btan relable transent slutns n usng the COMSOL.. Use f COMSOL Multphyscs The BTE equatn under a relaxatn tme apprxmatn s wrtten as f f f f v f () t t scat τ, where f s the carrer dstrbutn functn, f the equlbrum dstrbutn, v the carrer grup velcty vectr, and τ the relaxatn tme. The EPRT n terms f the phnn ntensty s wrtten as v () t τ, where the phnn ntensty s expressed as ( t, v, r) v hf ( t, v, r) D( ) / 4π, where h s the reduced Planck cnstant, the phnn angular frequency, and D () the densty f state per unt vlume. An equlbrum phnn ntensty n Eq. () s wrtten as dω 4π Ω 4π. (3) Fr ne-dmensnal prblems, Eq. () becmes vcsθ (4) t x τ, where θ s a plar angle between a phnn prpagatn drectn and the glbal heat transfer drectn ( x ). A slutn prcess f the -D BTE n Eq. (4) alng wth Eq. (3) can be summarzed as fllws: Wth an ntal value f, Eq. (4) s slved fr the phnn ntensty feld, ( t, x), fr a partcular phnn prpagatn angle, θ. After btanng the slutns fr all pssble sld angles, 4 π, can be calculated usng Eq. (3). Wth ths updated value f, the slutn prcess repeats untl the cnvergence f and. Therefre, the BTE and thus the EPRT are hghly nnlnear equatns t slve. Once the phnn ntensty s slved, heat flux and an equvalent temperature can be btaned by ntegratng the phnn ntensty ver the sld angles, 4 π, such that q csθ dω and (5) T u c Ω 4π 4π dω c v c v Ω 4π, (6) where c s the vlumetrc specfc heat. One f the mst wdely used methds fr slvng the BTE numercally s the dscrete rdnates methd (DOM). The DOM s a tl t transfrm the equatn f radatve transfer nt a set f smultaneus partal dfferental equatns. Ths s based n a dscrete representatn f drectnal varatn f ntensty. A slutn t the transprt prblem can be fund by slvng the equatn f radatve transfer fr a set f dscrete drectns spannng the entre sld angle. The ntegrals ver the sld angle can be apprxmated by numercal ntegratn usng Gauss-Legender quadratures. The BTEs n Eqs. (3)-(4) were slved by the DOM usng a bult-n feature f Ceffcent Frm PDEs n COMSOL Multphyscs. The spatal dman ( x ) was dscretzed nt n elements as a nrmal FE mesh refnement, whle the angular dman ( Ω ), referred t as rdnates, was dscretzed n the drectnal crdnates by dvdng the angular space nt a fnte number, m. The angular dman was dscretzed at Gaussan quadrature pnts wthn µ csθ. Fr each plar angle ( θ m ), we slved the system f nnlnear equatns fr the phnn ntensty, ( t, x, µ ), and then updated m the equlbrum phnn ntensty, ( t, x), usng
4 fr all angles. Eqs. (3) and (5) can be wrtten n dscrete frms as ( t, x) ( t, x, µ m ) w and m (7) m q ( t, x) π ( t, x, µ ) µ w m m m m m, (8) respectvely. The weght satsfes w. The slver used fr the steady-state and transent analyses was a bult-n drect slver, UMFPACK. Default values were used fr the slver parameters except that strct tme steps were taken by the slver and a maxmum BDF rder was set t. When the relatve errr f the calculated value f the equvalent equlbrum ntensty between tw teratn steps was less than a present tlerance, we cnsdered the prblem was cnverged and then calculated the equvalent temperature and heat flux usng Eqs. (6) and (8), respectvely. 5. Results and Dscussn We have slved the BTE equatns wth the COMSOL under bth steady-state and transent states. Spatal and angular dmans were dscretzed nt n elements and m quadrature pnts, respectvely. A quadratc Lagrange element was used fr the spatal FE mesh. The ttal degree f freedm fr the -D prblem becmes m ( n ). Fr the steady-state prblems, t was bserved that the wellcnverged slutns were btaned wth a relatvely carse mesh (5 elements) whle usng 6 quadrature pnts. The senstvty f mesh and quadrature pnts fr the transent analyss wll be stated later n ths sectn. Fgure shws steady-state temperature dstrbutns va nndmensnal emssve pwer ( e * ) alng the -D dman fr varus Knudsen numbers ( Kn ). Fr cmparsn purpse, an analytcal slutn was als btaned by usng an analgy t the radatn heat transfer [5] by slvng an ntegral equatn, ξ * * e η η η η η η ( ) E( ) e ( ) E( ) d (9), where η x L s the nndmensnal crdnate, where L s the length f the dman, n E ( x) µ exp( x µ ) dµ the expnental n ntegral, ξ the ptcal thckness, whch s an m nverse f the Knudsen number, Kn Λ L, where Λ s the average mean free path f * phnns, and e ( η ) [ e ( η) J q] [ J q J q] s the emssve pwer nndmensnalzed by a dfference f bundary heat fluxes ( J and J ) at bth ends f the -D dman. n Fgure, the analytc and numercal slutns fr varus Kn (.,,, ) were drawn wth markers and lnes, respectvely, whch shws an excellent agreement f the slutns wth each ther. Whle the cnventnal Furer slutn wuld yeld the temperature gradent varyng frm at η t at η, the slutn frm the BTE yelds a reduced temperature gradent. The larger the Knudsen number, the smaller the temperature gradent, whch s attrbuted frm mre ballstc phnn transprt as cmpared t the dffusve ne. Nndm. emssve pwer, e* q Kn. Kn Kn Kn Fgure. Steady-state nndmensnalzed emssve pwer dstrbutn fr varus Knudsen numbers. We studed the effects f refnements n terms f spatal FE meshes and angular quadrature pnts n the tme-dependent transent BTE slutns. The transent tme ( t ) was nndmensnalzed by the relaxatn tme f phnns ( τ ), s that t * t τ. Fgure shws nndmensnalzed temperature dstrbutns at t *. alng the -D dman fr Kn. We used three dfferent FE meshes (5, 6 and elements) and a fxed angular dvsn wth 6 Gaussan pnts ( n gp 6 ). The three FE meshes yeld nearly dentcal temperature slutns wth a decreasng trend frm the ht sde ( η ) t the cld sde ( η ). q
5 Smth slutns were btaned except at a regn near the ht bundary. An nset rectangular area n Fgure wthn η. 3 was zmed n and repltted n Fgure 3, whch clearly shws sgnfcant wggles near the ht bundary. Therefre, the FE mesh refnement des nt help smthng ut the wggles at all. The wgglng f the slutn, called a ray effect, s knwn t be attrbuted t nsuffcent refnement n the angular dman, and thus cannt be smthed ut wth the FE mesh refnement n the spatal dman. Nndmensnal temperature, θ elements 6 elements elements Fgure. Transent nndmensnal temperature dstrbutn predcted wth 5, 6 and fnte elements and 6 Gaussan pnts. Nndmensnal temperature, θ elements 6 elements elements...3 t *. alng the -D dman fr Kn predcted wth sx dfferent Gaussan quadrature pnts ( n gp 4, 8, 6, 3, 64, 8 ). We utlzed a Matlab-scrptng feature f COMSOL t mplement the large n gp n buldng the system f equatns. A fxed spatal dvsn wth 4 FE elements was used n these calculatns. An nset rectangular area wthn η. 3 was magnfed and repltted n Fgure 5. Althugh a hghly refned FE mesh (4 elements) was used, the carse angular dvsns allevate the wgglng f the slutn (e.g., n gp 4 ). The wggles near the ht end gradually decrease wth the ncrease f n gp, and thus the ray effect. Therefre, we fund that the spatal and angular refnements are ndependent f each ther, and thus d nt recmmend t use the hghly refned spatal mesh wth the carse angular mesh. Nndmensnal temperature, θ ngp4 ngp8 ngp6 ngp3 ngp64 ngp Fgure 4. Transent nndmensnal temperature dstrbutn predcted wth 4 fnte elements and sx Gaussan quadrature pnts ( n ) fr angular gp dscretzatn. Fgure 3. A magnfed vew f Fgure wthn η.3. T elmnate the ray effect, we made further refnement n the angular dman by usng mre quadrature pnts. Fgure 4 shws nndmensnalzed temperature dstrbutns at
6 Nndmensnal temperature, θ ngp4 ngp8 ngp6 ngp3 ngp64 ngp8...3 Fgure 5. A magnfed vew f Fgure 5 wthn η.3. Fgure 6 and Fgure 7 shw the temperature and heat flux dstrbutns, respectvely, at varus tme scales wth three dfferent Knudsen numbers. The result successfully reprduced the ne reprted earler by Chen [4], but ths tme usng the FE methd usng the COMSOL. The fgures shw the heat transfer frm the ht sde t the cld sde wth the fnte speed f phnn as the shrt term (small t * ) slutns ndcate, whch cannt be predcted by the FDE. Whle cnstant temperatures were appled as bundary cndtns at bth sdes, we can bserve the temperature jump at these bundares. The hgher the Kn, the larger the jump. The reasn fr the temperature jump at these emssve bundares was well explaned n earler wrk [7]. Fr a small Knudsen number at Kn., the temperature slutn appraches that f the steady-state FDE as the tme ncreases, whle fr a large Knudsen number at Kn, the temperature slutn appraches a nearly hrzntal lne, whch ndcates the ballstc heat transfer rather than the dffusve ne, and results n a small temperature gradent and thus a large thermal cnductvty. 6. Cnclusns We have slved the phnn Bltzmann transprt equatn fr the nanscale heat transfer prblems wth the COMSOL Multphyscs. One-dmensnal steady-state and transent prblems were successfully slved usng FEM and DOM fr spatal and angular dscretzatns, respectvely, by utlzng the bult-n feature f the COMSOL, Ceffcent Frm f PDE. A senstvty study was cnducted wth varus dscretzatn refnements fr dfferent values f the Knudsen number, whch s a measure f the nanscale regme. Sgnfcant ray effects were fund wth carse angular dscretzatns. Hwever, false scatterng due t carse spatal dscretzatn was nt severe wth the COMSOL slutn. t was fund that the spatal and angular refnements were ndependent f each ther, and t was nt recmmended t use the hghly refned spatal mesh wth the carse angular mesh n slvng the BTE wth the COMSOL. The present success wth the -D smulatn encurages us t apply t fr multdmensnal cases. (a) Nndmensnal temperature, θ t*. t*. t* t* (b) Nndmensnal temperature, θ t*. t*. t* t* (c) Nndmensnal temperature, θ t*. t* t* t* Fgure 6. Temperature dstrbutns at dfferent tme scales. (a) Kn, (b) Kn and (c) Kn..
7 (a) Nndmensnal heat flux, q* t*. t*. t* t* (b) Nndmensnal heat flux, q* t*. t*. t* t* (c) Nndmensnal heat flux, q* t*. t* t* t* Fgure 7. Heat flux dstrbutns at dfferent tme scales. (a) Kn, (b) Kn and (c) Kn.. 8. References. C. Ba, et al., On Hyperblc Heat Cnductn and the Secnd Law f Thermdynamcs, Jurnal f Heat Transfer, 7(), (995).. Y. M. Al, et al., Relatvstc Heat Cnductn, nternatnal Jurnal f Heat and Mass Transfer, 48(), (5). 3. A. Majumdar, Mcrscale Heat Cnductn n Delectrc Thn Flms, Jurnal f Heat Transfer, 5(), 7-6 (993). 4. G. Chen, Ballstc-Dffusve Equatns fr Transent Heat Cnductn frm Nan t Macrscales, Jurnal f Heat Transfer, 4(), 3-38 (). 5. G. Chen, Nanscale Energy Transprt and Cnversn: A Parallel Treatment f Electrns, Mlecules, Phnns and Phtns, Oxfrd Unversty Press (5). 6. G. Chen, Nnlcal and Nnequlbrum Heat Cnductn n the Vcnty f Nanpartcles, Jurnal f Heat Transfer, 8(3), (996). 7. R. Yang, et al., Smulatn f Nanscale Multdmensnal Transent Heat Cnductn Prblems Usng Ballstc-Dffusve Equatns and Phnn Bltzmann Equatn, Jurnal f Heat Transfer, 7(3), (5). 8. J. Y. Murthy, et al., Cmputatn f Sub-Mcrn Thermal Transprt Usng an Unstructured Fnte Vlume Methd, Jurnal f Heat Transfer, 4(6), 76-8 (). 9. M.-S. Jeng, et al., Mdelng the Thermal Cnductvty and Phnn Transprt n Nanpartcle Cmpstes Usng Mnte Carl Smulatn, Jurnal f Heat Transfer, 3(4), (8).. S. Mazumder, et al., Mnte Carl Study f Phnn Transprt n Sld Thn Flms ncludng Dspersn and Plarzatn, Jurnal f Heat Transfer, 3(4), ().. G. Chen, Sze and nterface Effects n Thermal Cnductvty f Superlattces and Perdc Thn- Flm Structures, Jurnal f Heat Transfer, 9(), -9 (997). 9. Acknwledgements Ths wrk was perfrmed under U.S. Ar Frce Cntract N. FA865-5-D-55.
8 Presented at the COMSOL Cnference 9 Bstn Nanscale Heat Transfer usng Phnn Bltzmann Transprt Equatn Sangwk Shn Ar Frce Research Labratry, Materals and Manufacturng Drectrate, AFRL/RX, Wrght-Pattersn AFB, OH, Unversty f Daytn Research nsttute, Daytn, OH emal : sangwk.shn@udr.udaytn.edu Ajt K. Ry Ar Frce Research Labratry, Wrght-Pattersn AFB, OH 9 COMSOL Cnference, Bstn, MA Octber 8-, 9
9 Outlne Backgrund nfrmatn. Descrptn f phnn Bltzmann transprt equatn (BTE). Mdelng and slutn prcedure f BTE usng COMSOL. Results Steady-state and transent prblems. ssues f refnement n spatal and angular dmans. Summary and cnclusns.
10 Furer Equatn (FE) Fr last tw centures, heat cnductn has been mdeled by Furer Eq (FE). T Cnservatn f energy: ρc q t Dffusve Furer s lnear apprxmatn f heat flux: q k T L >> Λ T T T α t x Parablc equatn > Dffusve nature f heat transprt. T Heat s effectvely transferred between lcalzed regns thrugh suffcent scatterng events f phnns wthn medum. Des nt hld when number f scatterng s neglgble. e.g., mean free path ~ devce sze (chp-package level). Bundary scatterng at nterfaces causng thermal resstance. Admts nfnte speed f heat transprt > Cntradct wth thery f relatvty. L Furer Equatn cannt be used fr small tme and spatal scales. 3
11 Hyperblc Heat Cnductn Equatn (HHCE) Reslve the ssue f the Furer equatn wth the nfnte speed f heat carrer. C t T T α t Defntn f heat flux: T q τ q k T, ( τ : relaxatn tme) t Hyperblc equatn > Wave nature f heat transprt. Called as Cattane equatn r Telegraph equatn. Fnte speed f heat carrers. Ad hc apprxmatn f heat flux defntn. Vlates nd law f thermdynamcs. f heat surce vares faster than speed f sund, heat wuld appear t be mvng frm cld t ht. ( C α τ ) HHCE: culd be used fr shrt tme scale, but nt fr shrt spatal scale. 4
12 Small Scale Heat Transprt (Tme & Space) Furer Equatn cannt be used fr small tme and spatal scales. HHCE: culd be used fr shrt tme scale, but nt fr shrt spatal scale. Needs equatns and methds fr small scale smulatn n terms f bth tme and space. Mlecular dynamcs smulatn. Accurate methd. Cmputatnally expensve. Sutable fr systems havng a few atmc layers r several thusands f atms. Nt sutable fr devce-level thermal analyss. Bltzmann Transprt Equatn (BTE). Ballstc-Dffusve Equatn (BDE). Smlar t Cattane Eq. (HHCE) wth a surce term. Derved frm BTE. Gd apprxmatn f BTE wthut nternal heat surce, dsturbance, etc. 5
13 Bltzmann Transprt Equatn (BTE) BTE: als called as equatn f phnn radatve transfer (EPRT). Equatn fr phnn dstrbutn functn: f t v f Can predct ballstc nature f heat transfer. Neglects wave-lke behavrs f phnn. Vald fr structures larger than wavelength f phnns (~ RT). Slutn methds: f t scat Determnstc: dscrete rdnates methd, sphercal harmncs methd. Statstcal: Mnte Carl. Much mre effcent than MD. f f τ Agrees well wth expermental data. L << Λ x T Ballstc T L 6
14 Detals f Bltzmann Transprt Equatn (BTE) Phnn ntensty: ( t, v, r) v f ( t, v, r) D( ) / 4π BTE becmes EPRT: t v τ, 4π Ω 4π dω *Equlbrum phnn ntensty determned by Bse-Ensten statstcs Fr -D, t v csθ x τ x θ Ι L Fr each angle (θ ), slve nn-lnear equatn wth teratns fr Slvng fr Ι. Updatng Ι ο. Heat flux: q Ω D csθ d dω 4π nternal energy: u D ( ct ) f D( ) d dω 4 π Ω 4π v d dω 7
15 Mdelng & Slutn Prcedure usng COMSOL t v csθ x τ Use a bult-n feature f COMSOL, Ceffcent Frm PDEs. The spatal dman s dscretzed usng FE mesh. The angular (mmentum) dman s dscretzed usng Gaussan quadrature pnts. Fr each angle (θ ), set up the BTE wth crrespndng ceffcents (µ csθ ) and BCs (Neumann vs. Drchlet). Calculate equlbrum phnn ntensty ( ) by numercal ntegratn f ( ) usng Gaussan quadratures. Slve. Drect slver (UMFPACK). Max. BDF rder. Pstprcess and vsualze the results. 8
16 Orgnal -D BTE: Nndmensnalze wth Splt nt () and (-) drectns: Dscretze angular space at Gaussan quadrature pnts: After FE run, pstprcess: 9 Detals f Slutn Prcedure < > 4 4 * *, Drchlet BCs : ) (, ) (, η η η η π σ π σ µ η µ µ η µ T T Kn t Kn t ) cs (, θ µ τ µ x v t L Kn L x t t Λ,, * η τ ), ( ), ( gp gp gp gp n n n n w w t q w w t µ µ π η η Kn t η µ *
17 Ceffcent Frm fr BTE (6 Fnte elements, 6 Gaussan Pnts) Fle Edt Optns Draw Physcs Mesh Slve Pstprcessng Multphyscs Help fvldel Tree F t: f:: < > 3 Geml DE, Ceffcent Frm (c) ;.. POE, Ceffcent Frm (c) :.. POE, Ceffcent Frm (c3)! POE, Ceffcent Frm (c4) ; POE, Ceffcent Frm (cs) : POE, Ceffcent Frm (c6)! POE, Ceffcent Frm (c7)! POE, Ceffcent Frm (c8) : POE, Ceffcent Frm (c9)! POE, Ceffcent Frm (clo)! POE, Ceffcent Frm (el l ) ~ POE, Ceffcent Frm (c)! POE, Ceffcent Frm (c3) ' POE, Ceffcent Frm (c4) :.. POE, Ceffcent Frm (cs) :.. PDE, Ceffcent Frm (c6) Equatn V'{-cV'ul -au! y) au! ~ V'ul f Subdman selectn Ceffcents POE ceffcents Ceffcent Value/ Expressn Descrptn l Dffusn ceffcent a t Absrptn ceffcent Surce term e. l Mass ceffcent v d. t Dampng/Mass ceffcent Grup: Select by grup ~ Actve n ths dman l Cnservatve flux cnvectn ceff. ~!Kn*mu(l) Cnvectn ceffcent y l Cnservatve flux surce term OK l [ Cancel [ Apply [ Help POE, Ceffcent Frm (c) Dependent varables: ul Default element type: Lagrange- Qua Wave extensn: Off Weak cnstrants: Off.. Mesh cnssts f 5 elements. Mesh cnssts f 3 el ements. ~~~--~:_~>~~M~e~s~h_:c~nssts f 6 el ements. lllr---r Nrmal.9 v!memry: (6/)
18 Steady-State Prblem: Analytc vs. Numercal Slutns Emssve pwer ~ Temperature Gradent e * ( η) e ( η) J J q J q q * * * e e ( η ) e ( η ) Nndm. emssve pwer, e* Kn Λ L Small temperature gradent at smaller scale Kn. Kn Kn Kn Gradent f nndm. emssve pwer, e* Analytc Numercal (ngp6) Numercal (ngp8) Knudsen number, KnΛ/L
19 Steady-State Prblem: Analytc vs. Numercal Slutns Heat flux Thermal cnductvty q * q J q J q k * q * L * e Nndmensnal heat flux, q* Analytc Numercal (ngp6) Numercal (ngp8). Knudsen number, KnΛ/L k* q*l/ e* Analytc Numercal (ngp6) Numercal (ngp8). Knudsen number, KnΛ/L
20 Transent Prblem: Effect f Spatal Refnement (Mre Fnte Elements) Refne spatal (x-) drectn wth n fnte elements (n5, 6, ). Dvde angular drectn wth 6 Gaussan pnts (ngp6). Nndmensnal temperature, θ t*. 5 elements 6 elements elements Nndmensnal temperature, θ elements 6 elements elements Ray effect....3 Spatal refnement leads t a smther slutn. Hwever, t des nt slve ray effect. 3
21 Transent Prblem: Effect f Angular Refnement (Mre Gaussan Pnts) Refne spatal (x-) drectn wth 4 FE elements. Dvde angular drectn wth ngp Gaussan pnts (ngp4,8,6,3,64,8). Nndmensnal temperature, θ t*. ngp4 ngp8 ngp6 ngp3 ngp64 ngp8 Nndmensnal temperature, θ ngp4 ngp8 ngp6 ngp3 ngp64 ngp Angular refnement reslve ray effect. Spatal and angular refnements are ndependent. Hghly refned spatal mesh wth carse angular mesh allevates slutn. 4
22 Results f BTE (Temperature & Heat Flux Dstrbutns wth Tme ncrease) Nndmensnal temperature, θ Nndmensnal heat flux, q* t*. t*. t* t* t*. t*. t* t* Nndmensnal temperature, θ Nndmensnal heat flux, q* t*. t*. t* t* Kn Kn Kn t*. t*. t* t* Nndmensnal temperature, θ Nndmensnal heat flux, q* t*. t* t* t* t*. t* t* t*
23 Cmparsns f FE, HHTC, BDE vs. BTE (Temperature & Heat Flux Dstrbutns fr 3 Kn) Nndmensnal temperature, θ BTE BDE Furer Cattane Kn t* Nndmensnal temperature, θ BTE BDE Furer Cattane Kn t* Nndmensnal temperature, θ BTE BDE Furer Cattane Kn. t* Nndmensnal heat flux, q* BTE BDE Furer Cattane Nndmensnal heat flux, q* BTE BDE Furer Cattane Nndmensnal heat flux, q* BTE BDE Furer Cattane
24 Summary and Cnclusns Nanscale smulatn s cnducted fr phnn heat transfer usng Bltzmann transprt equatn. Phnns fr delectrc, thermelectrc, semcnductr materals. Electrns fr metals. Gas mlecules fr rarefed gas states. Numercal slutn f BTE has been btaned fr -D prblem (bth steady-state and transent prblems). Temperature and heat flux dstrbutns frm the nanscale smulatn yeld cmpletely dfferent results frm the slutns frm Furer and Cattane equatns. Thermal cnductvty wll be dfferent, t. 7
Transient Conduction: Spatial Effects and the Role of Analytical Solutions
Transent Cnductn: Spatal Effects and the Rle f Analytcal Slutns Slutn t the Heat Equatn fr a Plane Wall wth Symmetrcal Cnvectn Cndtns If the lumped capactance apprxmatn can nt be made, cnsderatn must be
More informationConduction Heat Transfer
Cnductn Heat Transfer Practce prblems A steel ppe f cnductvty 5 W/m-K has nsde and utsde surface temperature f C and 6 C respectvely Fnd the heat flw rate per unt ppe length and flux per unt nsde and per
More informationV. Electrostatics Lecture 27a: Diffuse charge at electrodes
V. Electrstatcs Lecture 27a: Dffuse charge at electrdes Ntes by MIT tudent We have talked abut the electrc duble structures and crrespndng mdels descrbng the n and ptental dstrbutn n the duble layer. Nw
More informationNanoscale Heat Transfer using Phonon Boltzmann Transport Equation
Excerpt fr the Prceedings f the COMSOL Cnference 29 Bstn Nanscale Heat Transfer using Phnn Bltzann Transprt Equatin Sangwk Sihn *,2 and Ajit K. Ry Air Frce Research Labratry, 2 University f Daytn Research
More informationPhysic 231 Lecture 33
Physc 231 Lecture 33 Man pnts f tday s lecture: eat and heat capacty: Q cm Phase transtns and latent heat: Q Lm ( ) eat flw Q k 2 1 t L Examples f heat cnductvty, R values fr nsulatrs Cnvectn R L / k Radatn
More informationChapter 3, Solution 1C.
COSMOS: Cmplete Onlne Slutns Manual Organzatn System Chapter 3, Slutn C. (a If the lateral surfaces f the rd are nsulated, the heat transfer surface area f the cylndrcal rd s the bttm r the tp surface
More informationAnalytical Modeling of Natural Convection in Horizontal Annuli
Analytcal Mdelng f Natural Cnvectn n Hrzntal Annul Peter Teertstra, M. Mchael Yvanvch, J. Rchard Culham Mcrelectrncs Heat Transfer Labratry Department f Mechancal Engneerng Unversty f Waterl Waterl, Ontar,
More informationChapter 7. Systems 7.1 INTRODUCTION 7.2 MATHEMATICAL MODELING OF LIQUID LEVEL SYSTEMS. Steady State Flow. A. Bazoune
Chapter 7 Flud Systems and Thermal Systems 7.1 INTODUCTION A. Bazune A flud system uses ne r mre fluds t acheve ts purpse. Dampers and shck absrbers are eamples f flud systems because they depend n the
More informationSIMULATION OF THREE PHASE THREE LEG TRANSFORMER BEHAVIOR UNDER DIFFERENT VOLTAGE SAG TYPES
SIMULATION OF THREE PHASE THREE LEG TRANSFORMER BEHAVIOR UNDER DIFFERENT VOLTAGE SAG TYPES Mhammadreza Dlatan Alreza Jallan Department f Electrcal Engneerng, Iran Unversty f scence & Technlgy (IUST) e-mal:
More information2 Analysis of the non-linear aerodynamic loads of hypersonic flow. 1 General Introduction
4 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES PRELIMINARY STUDY OF NON-LINEAR AEROELASTIC PHENOMENA IN HYPERSONIC FLOW Zhang Wewe, Ye Zhengyn, Yang Yngnan Cllege f Aernautcs, Nrthwestern Plytechncal
More informationLecture 12. Heat Exchangers. Heat Exchangers Chee 318 1
Lecture 2 Heat Exchangers Heat Exchangers Chee 38 Heat Exchangers A heat exchanger s used t exchange heat between tw fluds f dfferent temperatures whch are separated by a sld wall. Heat exchangers are
More informationNatural Convection in a Horizontal Annulus with Oscillating Inner Cylinder Using Lagrangian-Eulerian Kinematics
Natural Cnvectn n a Hrzntal Annulus wth Oscllatng Inner Cylnder Usng Lagrangan-Euleran Knematcs Esam M. Alawadh Kuwat Unversty Mechancal Engneerng Department P. O. Bx # 5969, Safat, 3060 KUWAIT Abstract
More informationA New Method for Solving Integer Linear. Programming Problems with Fuzzy Variables
Appled Mathematcal Scences, Vl. 4, 00, n. 0, 997-004 A New Methd fr Slvng Integer Lnear Prgrammng Prblems wth Fuzzy Varables P. Pandan and M. Jayalakshm Department f Mathematcs, Schl f Advanced Scences,
More informationSpring 2002 Lecture #17
1443-51 Sprng 22 Lecture #17 r. Jaehn Yu 1. Cndtns fr Equlbrum 2. Center f Gravty 3. Elastc Prpertes f Slds Yung s dulus Shear dulus ulk dulus Tday s Hmewrk Assgnment s the Hmewrk #8!!! 2 nd term eam n
More informationProblem Set 5 Solutions - McQuarrie Problems 3.20 MIT Dr. Anton Van Der Ven
Prblem Set 5 Slutns - McQuarre Prblems 3.0 MIT Dr. Antn Van Der Ven Fall Fall 003 001 Prblem 3-4 We have t derve the thermdynamc prpertes f an deal mnatmc gas frm the fllwng: = e q 3 m = e and q = V s
More informationWp/Lmin. Wn/Lmin 2.5V
UNIVERITY OF CALIFORNIA Cllege f Engneerng Department f Electrcal Engneerng and Cmputer cences Andre Vladmrescu Hmewrk #7 EEC Due Frday, Aprl 8 th, pm @ 0 Cry Prblem #.5V Wp/Lmn 0.0V Wp/Lmn n ut Wn/Lmn.5V
More information_J _J J J J J J J J _. 7 particles in the blue state; 3 particles in the red state: 720 configurations _J J J _J J J J J J J J _
Dsrder and Suppse I have 10 partcles that can be n ne f tw states ether the blue state r the red state. Hw many dfferent ways can we arrange thse partcles amng the states? All partcles n the blue state:
More informationelement k Using FEM to Solve Truss Problems
sng EM t Slve Truss Prblems A truss s an engneerng structure cmpsed straght members, a certan materal, that are tpcall pn-ned at ther ends. Such members are als called tw-rce members snce the can nl transmt
More informationMode-Frequency Analysis of Laminated Spherical Shell
Mde-Frequency Analyss f Lamnated Sphercal Shell Umut Tpal Department f Cvl Engneerng Karadenz Techncal Unversty 080, Trabzn, Turkey umut@ktu.edu.tr Sessn ENG P50-00 Abstract Ths paper deals wth mde-frequency
More informationSection 3: Detailed Solutions of Word Problems Unit 1: Solving Word Problems by Modeling with Formulas
Sectn : Detaled Slutns f Wrd Prblems Unt : Slvng Wrd Prblems by Mdelng wth Frmulas Example : The factry nvce fr a mnvan shws that the dealer pad $,5 fr the vehcle. If the stcker prce f the van s $5,, hw
More information4DVAR, according to the name, is a four-dimensional variational method.
4D-Varatnal Data Assmlatn (4D-Var) 4DVAR, accrdng t the name, s a fur-dmensnal varatnal methd. 4D-Var s actually a smple generalzatn f 3D-Var fr bservatns that are dstrbuted n tme. he equatns are the same,
More informationComparison of Building Codes and Insulation in China and Iceland
Prceedngs Wrld Gethermal Cngress 00 Bal, Indnesa, 5-9 prl 00 Cmparsn f Buldng Cdes and Insulatn n Chna and Iceland Hayan Le and Pall Valdmarssn Tanjn Gethermal esearch & Tranng Centre, Tanjn Unversty,
More informationCHAPTER 3 ANALYSIS OF KY BOOST CONVERTER
70 CHAPTER 3 ANALYSIS OF KY BOOST CONERTER 3.1 Intrductn The KY Bst Cnverter s a recent nventn made by K.I.Hwu et. al., (2007), (2009a), (2009b), (2009c), (2010) n the nn-slated DC DC cnverter segment,
More information6. ELUTRIATION OF PARTICLES FROM FLUIDIZED BEDS
6. ELUTRIATION OF PARTICLES FROM FLUIDIZED BEDS Elutratn s the prcess n whch fne partcles are carred ut f a fludzed bed due t the flud flw rate passng thrugh the bed. Typcally, fne partcles are elutrated
More information3-42. Chapter 15 Steady Heat Conduction. Heat Conduction in Cylinders and Spheres
Chapter 5 Steady Heat Cnductn Heat Cnductn n Cylnders and Spheres 3-64C When the dameter f cylnder s very small cmpared t ts length, t can be treated as an ndefntely lng cylnder. Cylndrcal rds can als
More informationAdvances in Engineering Research (AER), volume 102 Second International Conference on Mechanics, Materials and Structural Engineering (ICMMSE 2017)
Secnd Internatnal Cnference n Mechancs, Materals and Structural Engneerng (ICMMSE 2017) Materal Selectn and Analyss f Ol Flm Pressure fr the Flatng Rng Bearng f Turbcharger Lqang PENG1, 2, a*, Hupng ZHENG2,
More informationRoot Locus Properties of Adaptive Beamforming and Capon Estimation for Uniform Linear Arrays
Root Locus Propertes of Adaptve Beamformng and Capon Estmaton for Unform Lnear Arrays Allan Stenhardt Alphatech phone: 703-284-8426 emal: astenhardt@dc.alphatech.com Abstract In ths paper we explore propertes
More informationExploiting vector space properties for the global optimization of process networks
Exptng vectr space prpertes fr the gbal ptmzatn f prcess netwrks Juan ab Ruz Ignac Grssmann Enterprse Wde Optmzatn Meetng March 00 Mtvatn - The ptmzatn f prcess netwrks s ne f the mst frequent prblems
More informationMonin Obukhov Similarity and Local-Free-Convection Scaling in the Atmospheric Boundary Layer Using Matched Asymptotic Expansions
OCTOBER 08 T O N G A N D D I N G 369 Mnn Obukhv Smlarty cal-free-cnvectn Scalng n the Atmspherc Bundary ayer Usng Matched Asympttc Expansns CHENNING TONG AND MENGJIE DING Department f Mechancal Engneerng
More informationDepartment of Civil Engineering & Applied Mechanics McGill University, Montreal, Quebec Canada
Department f Cvl Engneerng & Appled Mechancs McGll Unversty, Mntreal, Quebec Canada CIVE 90 THEMODYNAMICS & HEAT TANSFE Assgnment #6 SOUTIONS. Cnsder a.-m hgh and -m-wde duble-pane wndw cnsstng f tw 3-mmthck
More informationChapter 6 : Gibbs Free Energy
Wnter 01 Chem 54: ntrductry hermdynamcs Chapter 6 : Gbbs Free Energy... 64 Defntn f G, A... 64 Mawell Relatns... 65 Gbbs Free Energy G(,) (ure substances)... 67 Gbbs Free Energy fr Mtures... 68 ΔG f deal
More informationPT326 PROCESS TRAINER
PT326 PROCESS TRAINER 1. Descrptn f the Apparatus PT 326 Prcess Traner The PT 326 Prcess Traner mdels cmmn ndustral stuatns n whch temperature cntrl s requred n the presence f transprt delays and transfer
More informationIntroduction to Electronic circuits.
Intrductn t Electrnc crcuts. Passve and Actve crcut elements. Capactrs, esstrs and Inductrs n AC crcuts. Vltage and current dvders. Vltage and current surces. Amplfers, and ther transfer characterstc.
More informationA BESTEST VALIDATION STUDY OF THE DYNAMIC GROUND-COUPLED HEAT TRANSFER MODEL USED IN ACCURATE. Dong Chen 1. PO Box 56, Highett. Vic.
A BESTEST VALIDATION STUDY OF THE DYNAMIC GROUND-COUPLED HEAT TRANSFER MODEL USED IN ACCURATE Dng Chen CSIRO Energy Transfrmed Flagshp and CSIRO Ecsystem Scences PO Bx 56, Hghett. Vc. 390, Australa ABSTRACT
More informationFeedback Principle :-
Feedback Prncple : Feedback amplfer s that n whch a part f the utput f the basc amplfer s returned back t the nput termnal and mxed up wth the nternal nput sgnal. The sub netwrks f feedback amplfer are:
More informationSection 10 Regression with Stochastic Regressors
Sectn 10 Regressn wth Stchastc Regressrs Meanng f randm regressrs Untl nw, we have assumed (aganst all reasn) that the values f x have been cntrlled by the expermenter. Ecnmsts almst never actually cntrl
More informationApproach: (Equilibrium) TD analysis, i.e., conservation eqns., state equations Issues: how to deal with
Schl f Aerspace Chemcal D: Mtvatn Prevus D Analyss cnsdered systems where cmpstn f flud was frzen fxed chemcal cmpstn Chemcally eactng Flw but there are numerus stuatns n prpulsn systems where chemcal
More informationChem 204A, Fall 2004, Mid-term (II)
Frst tw letters f yur last name Last ame Frst ame McGll ID Chem 204A, Fall 2004, Md-term (II) Read these nstructns carefully befre yu start tal me: 2 hurs 50 mnutes (6:05 PM 8:55 PM) 1. hs exam has ttal
More informationXXIX CILAMCE November 4 th to 7 th, 2008 Maceió - Brazil
XXIX CILAMCE Nvember 4 th t 7 th, 8 Maceó - Bral ELECTROMAGNETIC SCATTERING PROBLEM SOLVED BY BOTH NODAL AND GALERKIN FEM-BEM APPROACHES M. M. Afns M. O. Schreder T. A. S. Olvera marcmatas@des.cefetmg.br
More informationRegression with Stochastic Regressors
Sectn 9 Regressn wth Stchastc Regressrs Meanng f randm regressrs Untl nw, we have assumed (aganst all reasn) that the values f x have been cntrlled by the expermenter. Ecnmsts almst never actually cntrl
More informationA/2 l,k. Problem 1 STRATEGY. KNOWN Resistance of a complete spherical shell: r rk. Inner and outer radii
Prblem 1 STRATEGY KNOWN Resstance f a cmplete sphercal shell: R ( r r / (4 π r rk sphere Inner an uter ra r an r, SOLUTION Part 1: Resstance f a hemsphercal shell: T calculate the resstance f the hemsphere,
More informationConservation of Energy
Cnservatn f Energy Equpment DataStud, ruler 2 meters lng, 6 n ruler, heavy duty bench clamp at crner f lab bench, 90 cm rd clamped vertcally t bench clamp, 2 duble clamps, 40 cm rd clamped hrzntally t
More informationLucas Imperfect Information Model
Lucas Imerfect Infrmatn Mdel 93 Lucas Imerfect Infrmatn Mdel The Lucas mdel was the frst f the mdern, mcrfundatns mdels f aggregate suly and macrecnmcs It bult drectly n the Fredman-Phels analyss f the
More informationDesign of Analog Integrated Circuits
Desgn f Analg Integrated Crcuts I. Amplfers Desgn f Analg Integrated Crcuts Fall 2012, Dr. Guxng Wang 1 Oerew Basc MOS amplfer structures Cmmn-Surce Amplfer Surce Fllwer Cmmn-Gate Amplfer Desgn f Analg
More informationCTN 2/23/16. EE 247B/ME 218: Introduction to MEMS Design Lecture 11m2: Mechanics of Materials. Copyright 2016 Regents of the University of California
Vlume Change fr a Unaxal Stress Istrpc lastcty n 3D Istrpc = same n all drectns The cmplete stress-stran relatns fr an strpc elastc Stresses actng n a dfferental vlume element sld n 3D: (.e., a generalzed
More informationCIRCLE YOUR DIVISION: Div. 1 (9:30 am) Div. 2 (11:30 am) Div. 3 (2:30 pm) Prof. Ruan Prof. Naik Mr. Singh
Frst CIRCLE YOUR DIVISION: Dv. 1 (9:30 am) Dv. (11:30 am) Dv. 3 (:30 m) Prf. Ruan Prf. Na Mr. Sngh Schl f Mechancal Engneerng Purdue Unversty ME315 Heat and Mass ransfer Eam #3 Wednesday Nvember 17 010
More informationIGEE 401 Power Electronic Systems. Solution to Midterm Examination Fall 2004
Jós, G GEE 401 wer Electrnc Systems Slutn t Mdterm Examnatn Fall 2004 Specal nstructns: - Duratn: 75 mnutes. - Materal allwed: a crb sheet (duble sded 8.5 x 11), calculatr. - Attempt all questns. Make
More informationWater vapour balance in a building moisture exposure for timber structures
Jnt Wrkshp f COST Actns TU1 and E55 September 21-22 9, Ljubljana, Slvena Water vapur balance n a buldng msture expsure fr tmber structures Gerhard Fnk ETH Zurch, Swtzerland Jchen Köhler ETH Zurch, Swtzerland
More informationPHYSICS 536 Experiment 12: Applications of the Golden Rules for Negative Feedback
PHYSICS 536 Experment : Applcatns f the Glden Rules fr Negatve Feedback The purpse f ths experment s t llustrate the glden rules f negatve feedback fr a varety f crcuts. These cncepts permt yu t create
More informationThermodynamics of Materials
Thermdynamcs f Materals 14th Lecture 007. 4. 8 (Mnday) FUGACITY dg = Vd SdT dg = Vd at cnstant T Fr an deal gas dg = (RT/)d = RT dln Ths s true fr deal gases nly, but t wuld be nce t have a smlar frm fr
More information55:041 Electronic Circuits
55:04 Electrnc Crcuts Feedback & Stablty Sectns f Chapter 2. Kruger Feedback & Stablty Cnfguratn f Feedback mplfer S S S S fb Negate feedback S S S fb S S S S S β s the feedback transfer functn Implct
More informationTubular Flow with Laminar Flow (CHE 512) M.P. Dudukovic Chemical Reaction Engineering Laboratory (CREL), Washington University, St.
Tubular Flw wth Lamnar Flw (CHE 5) M.P. Dudukvc Chemcal Reactn Engneerng Labratry (CREL), Washngtn Unversty, St. Lus, MO 4. TUBULAR REACTORS WITH LAMINAR FLOW Tubular reactrs n whch hmgeneus reactns are
More informationEE 221 Practice Problems for the Final Exam
EE 1 Practce Prblems fr the Fnal Exam 1. The netwrk functn f a crcut s 1.5 H. ω 1+ j 500 Ths table recrds frequency respnse data fr ths crcut. Fll n the blanks n the table:. The netwrk functn f a crcut
More informationA Note on the Linear Programming Sensitivity. Analysis of Specification Constraints. in Blending Problems
Aled Mathematcal Scences, Vl. 2, 2008, n. 5, 241-248 A Nte n the Lnear Prgrammng Senstvty Analyss f Secfcatn Cnstrants n Blendng Prblems Umt Anc Callway Schl f Busness and Accuntancy Wae Frest Unversty,
More informationVan der Waals-coupled electronic states in incommensurate double-walled carbon nanotubes
Kahu Lu* 1, Chenha Jn* 1, Xapng Hng 1, Jhn Km 1, Alex Zettl 1,2, Enge Wang 3, Feng Wang 1,2 Van der Waals-cupled electrnc states n ncmmensurate duble-walled carbn nantubes S1. Smulated absrptn spectra
More informationBME 5742 Biosystems Modeling and Control
BME 5742 Bsystems Mdeln and Cntrl Cell Electrcal Actvty: In Mvement acrss Cell Membrane and Membrane Ptental Dr. Zv Rth (FAU) 1 References Hppensteadt-Peskn, Ch. 3 Dr. Rbert Farley s lecture ntes Inc Equlbra
More informationinhomogeneous media using the conjugate gradient
Rad Scence Vlume 34 Number 6 Pages 1339-1347 Nvember-December 1999 Slvng the vlume ntegral equatn n axsymmetrc nhmgeneus meda usng the cnjugate gradent fast Hankel transfrm methd $. Y. Chen and W. C. Chew
More informationA Proposal of Heating Load Calculation considering Stack Effect in High-rise Buildings
A Prpsal f Heatng Lad Calculatn cnsderng Stack Effect n Hgh-rse Buldngs *Dsam Sng 1) and Tae-Hyuk Kang 2) 1) Department f Archtectural Engneerng, Sungkyunkwan Unversty, 2066 Sebu-r, Jangan-gu, Suwn, 440-746,
More informationBig Data Analytics! Special Topics for Computer Science CSE CSE Mar 31
Bg Data Analytcs! Specal Tpcs fr Cmputer Scence CSE 4095-001 CSE 5095-005! Mar 31 Fe Wang Asscate Prfessr Department f Cmputer Scence and Engneerng fe_wang@ucnn.edu Intrductn t Deep Learnng Perceptrn In
More informationShell Stiffness for Diffe ent Modes
Engneerng Mem N 28 February 0 979 SUGGESTONS FOR THE DEFORMABLE SUBREFLECTOR Sebastan vn Herner Observatns wth the present expermental versn (Engneerng Dv nternal Reprt 09 July 978) have shwn that a defrmable
More informationChalcogenide Letters Vol. 11, No. 7, July 2014, p THE QUANTUM MECHANICAL STUDY OF CADMIUM SULFUR NANOPARTICLES IN BASIS OF STO s
Chalcgende Letters Vl. 11, N. 7, July 014, p. 35-364 THE QUNTUM MECHNICL STUDY OF CDMIUM SULFUR NNOPRTICLES IN BSIS OF STO s M.. RMZNOV *, F. G. PSHEV,. G. GSNOV,. MHRRMOV,. T. MHMOOD Baku State Unversty,
More informationNomenclature: number of electrons e -1. electron charge F constant number, (columbs/moles of e -1 ) atomic number g
Quanttatve Analyss f Irreversbltes Causes Vltage Drp n Fuel cell (Smulatn) Hssen Ghadaman*, Dr. Yadlah Sabh** Department f Energy Engneerng, Scence and Research Branch Azad Unversty, Islamc Republc f IRAN
More informationNUMERICAL 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 informationFaculty of Engineering
Faculty f Engneerng DEPARTMENT f ELECTRICAL AND ELECTRONIC ENGINEERING EEE 223 Crcut Thery I Instructrs: M. K. Uygurğlu E. Erdl Fnal EXAMINATION June 20, 2003 Duratn : 120 mnutes Number f Prblems: 6 Gd
More informationINTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 1, No 4, 2011
Behavr f hgh frequenc mdal respnses n nn lnear sesmc analss Dhleep. M 1, Trved. A 2, Bse.P.R 3 1 Assstant Prfessr, Department f Cvl Engneerng, SCMS Schl f Engneerng and Technlg, Karukutt, Ernakulam, Kerala
More informationDepartment of Applied Mathematics, Tsinghua University Beijing , People's Republic of China Received 17 August 1998; accepted 10 December 1998
Cmput. Methds Appl. Mech. Engrg. 79 (999) 345±360 www.elsever.cm/lcate/cma The dscrete art cal bndary cndtn n a plygnal art cal bndary fr the exterr prblem f Pssn equatn by usng the drect methd f lnes
More informationExercises H /OOA> f Wo AJoTHS l^»-l S. m^ttrt /A/ ?C,0&L6M5 INFERENCE FOR DISTRIBUTIONS OF CATEGORICAL DATA. tts^e&n tai-ns 5 2%-cas-hews^, 27%
/A/ mttrt?c,&l6m5 INFERENCE FOR DISTRIBUTIONS OF CATEGORICAL DATA Exercses, nuts! A cmpany clams that each batch f ttse&n ta-ns 5 2%-cas-hews, 27% almnds, 13% macadama nuts, and 8% brazl nuts. T test ths
More informationAnalysis The characteristic length of the junction and the Biot number are
-4 4 The temerature f a gas stream s t be measured by a thermule. The tme t taes t regster 99 erent f the ntal ΔT s t be determned. Assumtns The juntn s sheral n shae wth a dameter f D 0.00 m. The thermal
More informationMicrofacet models for reflection and refraction
Mcrfacet mdels fr reflectn and refractn Steve Marschner Crnell Unversty CS 5625 Sprng 2016 (based n presentatn fr Walter, Marschner, L, and Trrance EGSR 07) Mcrfacet scatterng mdels Rugh delectrc surface
More informationMicrofacet models for reflection and refraction
Mcrfacet mdels fr reflectn and refractn Steve Marschner Crnell Unversty CS 6630 Sprng 2012 (based n presentatn fr Walter, Marschner, L, and Trrance EGSR 07) Mcrfacet scatterng mdels Rugh delectrc surface
More informationPhysics 107 HOMEWORK ASSIGNMENT #20
Physcs 107 HOMEWORK ASSIGNMENT #0 Cutnell & Jhnsn, 7 th etn Chapter 6: Prblems 5, 7, 74, 104, 114 *5 Cncept Smulatn 6.4 prves the ptn f explrng the ray agram that apples t ths prblem. The stance between
More informationThe lattice Boltzmann method as a basis for ocean circulation modeling
Jurnal f Marne Research, 57, 503 535, 999 The lattce Bltzmann methd as a bass fr cean crculatn mdelng by Rck Salmn ABSTRACT We cnstruct a lattce Bltzmann mdel f a sngle-layer, reduced gravty cean n a square
More informationDimensional Analysis.
Dmensnal nalyss. Unts, r hw we chse t measure magntudes f quanttes such as sze, temperature, r mass are nt requred fr laws f nature t functn prperly. Snce equatns f mmentum, energy, and mass cnseratn are
More informationENGI 4421 Probability & Statistics
Lecture Ntes fr ENGI 441 Prbablty & Statstcs by Dr. G.H. Gerge Asscate Prfessr, Faculty f Engneerng and Appled Scence Seventh Edtn, reprnted 018 Sprng http://www.engr.mun.ca/~ggerge/441/ Table f Cntents
More informationA method of constructing rock-analysis diagrams a statistical basks.
130 A methd f cnstructng rck-analyss dagrams a statstcal basks. 0T~ By W. ALF~.D ll~ch).ra)so.~, ~.Se., B.Se. (Eng.), F.G.S. Lecturer n Petrlgy, Unversty Cllege, Nttngham. [Read January 18, 1921.] D R.
More informationWYSE Academic Challenge 2004 Sectional Physics Solution Set
WYSE Acadec Challenge 004 Sectnal Physcs Slutn Set. Answer: e. The axu pssble statc rctn r ths stuatn wuld be: ax µ sn µ sg (0.600)(40.0N) 4.0N. Snce yur pushng rce s less than the axu pssble rctnal rce,
More informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More information4.8 Degradation of Elastomers by Heat and/or Radiation
4.8 Degradatn f Elastmers by Heat and/r Radatn M.It Japan Atmc Energy Research Insttute, Nuclear Educatn Center 2-28-49, Hnkmagme, Bunkyu-ku Tky, 113, JAPAN Abstract Ths artcle studed sme prblems n the
More informationFinal Exam Spring 2014 SOLUTION
Appled Opts H-464/564 C 594 rtland State nverst A. La Rsa Fnal am Sprng 14 SOLTION Name There are tw questns 1%) plus an ptnal bnus questn 1%) 1. Quarter wave plates and half wave plates The fgures belw
More informationTHEORY OF HYPERBOLIC TWO-TEMPERATURE GENERALIZED THERMOELASTICITY
Materals Physcs Mechancs 40 (08) 58-7 Receed: February 3, 08 THEORY OF HYPERBOLIC TWO-TEMPERATURE GENERALIZED THERMOELASTICITY Hamdy M. Yussef,, Alaa A. El-Bary 3 Mathematcs Department, Faculty f Educatn,
More informationControl of Program Motion of Dynamic Systems Using Relative Motions
TECHNISCHE MECHANIK, Band 8, Heft 3-4, (8), - 8 Manuskrptengang: 5. Oktber 7 Cntrl f Prgram Mtn f Dynamc Systems Usng Relatve Mtns D Sanh, Dnh Van Phng, D Dang Kha, Tran Duc In ths paper, the prblem f
More informationTowards Robust Two-Level Methods For Indefinite Systems
Twards Rbust Tw-Level Methds Fr Indefnte Systems J. Fsh, Y. Qu and A. Suvrv Departments f Cvl, Mechancal and Aerspace Engneerng Rensselaer Plytechnc Insttute Try, NY 28 We present a black-bx tw-level slver
More informationCircuits Op-Amp. Interaction of Circuit Elements. Quick Check How does closing the switch affect V o and I o?
Crcuts Op-Amp ENGG1015 1 st Semester, 01 Interactn f Crcut Elements Crcut desgn s cmplcated by nteractns amng the elements. Addng an element changes vltages & currents thrughut crcut. Example: clsng a
More informationChemistry 20 Lesson 11 Electronegativity, Polarity and Shapes
Chemistry 20 Lessn 11 Electrnegativity, Plarity and Shapes In ur previus wrk we learned why atms frm cvalent bnds and hw t draw the resulting rganizatin f atms. In this lessn we will learn (a) hw the cmbinatin
More informationAnalysis of ellipsometric data obtained from curved surfaces. J. Křepelka
Analyss f ellpsmetrc data btaned frm curved surfaces J. Křepelka Jnt Labratry f Optcs f Palacky Unversty and Insttute f Physcs f Academy f Scences f the Czech Republc, 7. lstpadu 5a, 77 7 Olmuc, Czech
More informationSolving the VAR Sources Planning Problem in Multiple Load Cases Using Genetic Algorithm Based Method
Slvng the VAR Surces Plannng Prblem n Multple Lad Cases Usng Genetc Algrthm Based Methd Ch -Hsn Ln* and Shn-Yeu Ln * Department f Electrnc Engneerng Ka Yuan Unversty Kahsung, Tawan, R. O. C. e-mal: chsnln@ee.yu.edu.tw
More informationNumerical 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 informationFormal solvers of the RT equation
Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan
More informationC. Ozgur Colpan Ibrahim Dincer Feridun Hamdullahpur
THERMODYNAMIC MODELING O DIRECT INTERNAL REORMING PLANAR SOLID OXIDE UEL CELLS WITH ANODE RECIRCULATION C. Ozgur Clpan Ibrahm Dncer erdun Hamdullahpur OUTLINE INTRODUCTION TO SOCs Advantages, dsadvantages
More informationOut-of-plane orbital maneuvers using swing-bys with the Moon
Jurnal f Physcs: Cnference Seres PAPER OPEN ACCESS Out-f-plane rbtal maneuvers usng swng-bys wth the Mn Related cntent - Pwered Swng-By Maneuvers arund the Mn A F Slva, A F B A Prad and O C Wnter cte ths
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More informationChapter 4: Damping . (4-1)
6 Chapter 4: Dampng 4.: The Absrbng Atm and ts Envrnment In chapter, we saw hw the absrptn ceffcent fr a nn-nteractng atmc speces at rest s gven by the cmbnatn f the strength f the transtn, the ppulatn
More informationCHAPTER 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 informationFeature Selection for Time Series Modeling *
Jurnal f Intellgent Learnng Systems Applcatns, 013, 5, 15-164 http://dxdrg/10436/lsa01353017 Publshed Onlne August 013 (http://wwwscrprg/urnal/lsa) Feature Selectn fr Tme Seres Mdelng * Qng-Gu Wang, Xan
More informationIntegrating Certified Lengths to Strengthen Metrology Network Uncertainty
Integratng Certfed engths t Strengthen Metrlgy Netwrk Uncertanty Authrs: Jseph Calkns, PhD New Rver Knematcs je@knematcs.cm Sctt Sandwth New Rver Knematcs sctt@knematcs.cm Abstract Calbrated and traceable
More informationTHERMAL TEST LEVELS & DURATIONS
PREFERRED RELIABILITY PAGE 1 OF 7 PRACTICES PRACTICE NO. PT-TE-144 Practice: 1 Perfrm thermal dwell test n prtflight hardware ver the temperature range f +75 C/-2 C (applied at the thermal cntrl/munting
More informationA Few Basic Facts About Isothermal Mass Transfer in a Binary Mixture
Few asic Facts but Isthermal Mass Transfer in a inary Miture David Keffer Department f Chemical Engineering University f Tennessee first begun: pril 22, 2004 last updated: January 13, 2006 dkeffer@utk.edu
More informationChapter 13: Multiple Regression
Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to
More informationAir Force Research Laboratory /RVBXR 29 Randolph Road Hanscom AFB, MA
REPRT DCUMENTATN PAGE Frm Apprved MBN. 0704-0-088 The publc reprtng burden fr ths cllectn f nfrmatn s estmated t average hur per respnse, ncludng the tme fr revewng nstructns, searchng exstng data surces,
More informationWhere a licence is displayed above, please note the terms and conditions of the licence govern your use of this document.
Mdellng the seepage flw durng cassn nstallatn n a natural seabed Faramarz, Asaad; Faz, Khyar; Drar, Samr; Mehravar, Mura; Harreche, Ouahd Dcument Versn Publsher's PDF, als knwn as Versn f recrd Ctatn fr
More information