Muroran-IT Academic Resources Archive Title Author(s) Citation Issue Date URL Rights Type 蓄熱槽中の蓄熱要因レンガの表面形状と配列に対する熱 伝達および流れ特性の研究 アムバリータ, ヒムサール; 鈴木, 淳; 岸浪, 紘機; 台丸谷, 政志; 松坂, 覚 寒地技術論文・報告集, 21, pp.537-542, 2005 2005 http://hdl.handle.net/10258/1614 著作権は著作者に帰属。ダウンロードを含むこの著作物 の利用は、著作権法の私的使用及び引用の範囲に限り認 められます。その範囲外の利用については、著作者の承 諾が必要です。 Journal Article See also Muroran-IT Academic Resources Archive Copyright Policy Muroran Institute of Technology 第21回寒地技術シンポジウム 2005年12月5,8,7日 CTCO5−I−084 蓄熱槽中の蓄熱要因レンガの表面形状と酉日 列に対する熱伝達および流れ特性の研究 ヒムサール・アムバリータ (室蘭工業大学) 〃 鈴 木 淳 〃 岸浪紘機 台丸谷政志 ( 〃 松 坂 覚 ⑧ ! ‘ 〃 ANumericaIStudyonFlowandHeatTransferCharacteristics 骨 CRTC forBrick'sA『TangementandSurfaceFominaRockBedHeat StorageSystem HAmbarita(MuroranlnstituteofTechnology) J、Swn'kl(MuroranInstituteofTechnology) K、Kishin潟mi(MuroranlnstituteofTechnology) M 、 D a i m a r u y a ( M u r o r a n l n s t i t u t e o f T e c h n o l o g y ) S、Matu咳負い(MuroranInstituteofTechnology) Abstract Thispaperattemptstostudynumerica]lyheattrans他randnuidmowcharacteristiofbrbrick(s)elementinarock b e d h e a t s t o r a g e t a n k s y s t e m o i l h e e m p h a s i s o f t h i s s t u d y i s t o m a k e c l e a r t h e e f f e c t o f t h e b r i c k s u r f a c e s a n d a r r a n g e m e n t o n t h e t o t a l p e r f b r m a n c e o f t h e h e a t s t o r a g e t a n k 、 A p l 8 j n s u r f a c e b r i c k w i t h t w o t y p e s d i f f e r e n t p o s i t i o n s i . e 、 i n c l i n e a n d paranelwiththestreamflowwerestudied,Thebrickwasheatedbyhotair50oCandvelocity0.1m/shominitialcondition 20oC,Lammar,transientandtwo-dimensionalgovermngequationsaresolvedbyusingSIMPLEalgorithm、Thecalculation resultsshowthatthepresenceofthedimplesonbricksur色cesdidnotenhanoetheheattransfbrratefromtheairmtothe brick、Theinc]inepositionofthebricktowardupstreamenhancestheheattransibrratesigniEcantly・Basedonthese r e s u l t s a f U n s c a l e h e a t s t o r a g e s y s t e m s c o n s i s t o f a s e t s o f p l a i n b r i c k s w e r e a l s o s t u d i e d 、 A s e ] 叩 e c t e d , t h e h e a t s t o r a g e tankwithinclmedpositionsofbricksrevealthebetterperfbrmanceofstoringheatfromthehotajrcomparedtothe p a r a l l e l p o s i t i o n . Keywords:bricks,heatstorage,plainsurface/discontinuoussurfaceofbrick,inclinedarrangement 1.Introduction Solarenergyhasrecejvedbigconsiderationsln recentyearsdueto]imitationoffossilsourceenergy、 Althoughthesolarenergyisabundantbuttherearemany problemstobeagednstwhenusmgsolarenergyNoticeable example,lowqua]ityofheatandstronglydependon weather・Basedonthoseproblemsithasbeenrealizedthat difficulttodoconverslonsolarenergytothehighleve] energysuchasmechanicalenergyorelectric8jenergy、 FromefEciencyandeffbctivenesspointofview, usmgsolarenergyasaheatisbetterthantransformitto theotherenergyfbrmTherefbre,usethesolarenergyasa heatsourceforroomheatinginthewinterseasonlsa ofasetofbricksasaheatstoragetankelement, Theobjectiveofthisresearchistostudytllenuid Howandheattrans色rcharacteristicoftheheatstorage tankbyusmgasetofthebIicks・Theemphasisofthis researchistomakeclearthee鮭ctofthebrickpositions andbricksurfacesonthefluidnowandheattransfer charactelistics、Theresultscanbeexpectedfromthis studyistosupplyanecessaryinformationofbricksheat storagetankdesignandoptimization. 2.Problemdehmtion皿dmode]jng lnordertoacquiretheseobjectivestwommnparts ofnumericalcalculationswerecarriedout・Thefirst promisingapp]ication、However』thisapplicationstin calculationlsabasicstudyontheinnuenceofabrick needsanmnOvation.Aheatstoragetankisneeded becausethedjfferentoftheusingtimeandtheav8jlable timeofthesolarenergy,AmaJorpartofstoringenergy problemliesinthenuidnowcharacteristicandheat p o s i t i o n a n d i t s s u r f a c e s i n d j y i d u a l l y , S e c o n d i s c a l c u l a t i o n foraheatstoragetanksystemwhichcontamsasetof bricks、Thebrickssur色cesandpositionsinthesecond transfer,Powerupperibrmanceoftheheatstoragetankis amosttol〕icinthisarea・Thereareseveraltechmquesm ordertostoretheheati.e,usmgPCMmaterial,porous medium,andsetsofbrick、Inthjspaper,wefbcusonusing Inthefirstcalculation,abrickisheatedbyhotair 500Cwithveloc此y0.1m/sh,ominitialtemperature200C・ Ahorizontalbrickwiththreediffもrentsurfacesl‘e・plain surface,discontinuoussurfaceswithsquaredimples,and calculationarebasedontheresultofthefirstcalculation・ AMBARITAHimsar,MuroranlnstituteofTechnologybMuroranshiMizumotocho27・lO50−0071JAPANTel、0143−46.5302 −537− 第21回寒地技術シンポジウム(2005) triangledimPleswerestudied・Thentheplamblickin Combinationofthein−mneandstaggeredwith inc]inedPosition(450)towardtotheupstreamwerealso thehorizontalandinc]inedpositionsofthebrickswere studied・PositionsofthebricksaredepictedinFigurel・ considered、Thesecombinationsweredjvidedintofour rIiffbrentmodelsandnamedasfollows・ ② ひ。。=0.1"1/s 鵬 ,│ Model-1:Heatstoragewithin-Iinehorizontalbricks Model-2:HeatStoragewithstaggeredhorizontalbricks Model-3:HeatStoragewithin-linemc]inedbrioks ハーl︸杵goヨコTl上 T=50oC ↓↓↓↓↓↓↓ トー300-う( Model-4:HeatStoragewithstaggeredinchnedbricks. A]lofthemodelsaredepictedinFigure2,Figure 3,Figure4,andFigure5. 一L=1000mm ( a ) 灘 D ← 3 0 0 う 1 −> T ‐ −> 即=U鯛 許 T=7m 則= 鍵‘ →……、…昼………蕊…J≦ Al〈産1000mm>,B Figure4Model-8:In-hneinclmedbricksarrangement ( b ) F i g u r e ’ P o s i t i o n o f t h e b r i c k ( a ) h o r i z o n t a l p O s i t i o n a n d ( b ) i n c l i n e d p o s i t i o n Inthesecondpartofcalculation,aheatstorage tankthatcontainsasetofplainbrioksisheatedbyhotair 500CwithvelocityO,1m/shominitialtemperature200C・ ThereasonofusingtheplEdnsurfacesofthebrickswj1lbe explainednextintheresultpartofthisPaper. 0 伽鋼職蝿 一癖 癖癖癖 蝿一 睡睡怖 n帥 蝿緬緬麺哩哩 施唖 廼﹄ 堅 = 0 , " = 0 , v = o 綱蝿蝿鯛一 一綱 綱蝿 蝿蝿 蝿 罪一斗 一癖 癖蝿 蝿蝿 蝿l引 研蝿蝿蝿緬一 Figure5Model-4Staggeredinchnedbricksarrangement 獄側鵬職|剛剛脇蝿、棚、 │蕊職||鯛綱|蝿剛I剛、M1、剛I’ 3.Mathematicalmodel Themowisassumedtobetransient,laminar,and 0 帥一欽 凡a−a p加一欺 冴一 一釦 諏 釦m 両m 耐m 両両 _ 2 8 0 m m − * , 。 。 」 ÷ 3 . J 緬一蝿 ヨーミ 銅蝿 一蝿 蝿射 蝿蝿趣 鯛鯛唖 T=、i, 占さ一コ弓 壷9−ミュ 弓菅 ハ出さ−コ弓.m ー’一︲|︲|︲|︲一︲−1〃一・〃 〃 = U 、 , AliWaIls 蕊 、=U1月 T=ヱパ two-dimensional・Thecompressibi]ity,radiationheat exchange,buoyancyfbrce,anddissipationsareneg1igible, Anofthethermalpropertiesareoonstant、Thegovernmg e q u a t i o n s a r e ・ Continuityequation L26onun- 1180mm 2L+聖=0 ( 1 ) aray FiRure2Modell:In-Ijnehorizontalblicksarrangement 1 調一が1唖 飾7座が︾ 乱 窯池 聖馴I哩 538 ’一一n arrangement 1|βk|碑︾ Figure8Model-2:Staggeredhorizontalbricks 脚一妙趣釘﹄砂姉 │劇雑鼎錨I−1抑日:捧紳:= 1 a十m 帥一欲平諏一欧可 帥斗団蝿訂﹃釘恥 蝿癖銅一騨蝿識 ︾蝿蝿鋼蝉鏑一 瞬唖蝿一蝿蝿銅 ︾蝿麺鍛鯛蝿一 r=m“ 癖癖麺一蝿蝿蝿 ィ=U伽 XandYmomentumequations ( 2 ) ( 8 ) ( 4 ) 第21回寒地技術シンポジウム(2005) r""……’二,-0‘祭-. ( 5 ) Theinlet:〃=U…v=0,T=Zo ( 6 ) Figure6Thebluecolorindicatesthelowertemperature andtheredcolorhighertemperature、Theincrementln thefigureislO%ofthehighestandlowesttemperature, Thefigureshowsthatattheleftsideofthebrick (upstream)theisotherm]inesarecloselyduetothehigher T h o ・ u t l o 纏 器 − 0 器 0 器 … d , = 恥 ( 7 ) heattransfbrrate・Basedonthedensityoftheisotherm linesinsidethebrick,thefiguresshowthattheheat Brickregionsz4=,ノー0 ( 8 ) transferratethroughtheleftandrightsides(upstream anddownstream)ofthebrickarehigherthanthatofthe topandthebottomsides、n1eisotherm]inesaresmooth duetoplainsurface・Thereisnosignificantdjfferentof Totalheatstoredbybricksiscalculatedby: isotherm1inesfbrmofthedisContinuoussurfacebricksand g‘。,=Zpか虚心雲c‘(zルー囚") theplainbricks(theisotherm]inesfordiscontinuous ( 9 ) bricksarenotshown). 融串●N●.・き●・聯・り Therateoftheheattransfbrbetweenthebrickandtheair # i s c a l c u l a t e d b y : ,富z噺…+z鋤劃…!” q 1鑑謹 癖韓『;#諜緯撚羅;;M蕊?『蕊蕊織罰 21.8. Equation(10)isusedwhentheparticularnodeoftl'e 蝉#溌燕識溌憩 4 2 . 3 7 4r=lh brickshasbothverticalandhorizontalcontactsurfacesto 識識韓韓韓:識騨i鍔剥 2 4 ‘ 7 9 イ 4 , 1 5 7 4/=2h theair・Ifthenodehasonlyonecontactsurfacetotheair onlyonepartoftheequation(10)isapplied,verticalor h o r i z o n t a l o n l y ・ Thermalproperties,kqp,P,respectively t h e r m a l c o n d u c t i v i t y ( W / m ・ ' 0 , s p e c i h c h e a t ( k J / k 9 . K ) , a n d densityOKg/m3)ofthematenalsarelistedasfollows・The 燃溌寵鍵溌鍵錘蕊謹選騨織鐘鍛錨脹蝿翻蝿禰癖閥零関顕瞬詞..."..:溺・鑑銅 2 7 . 6 9 4 5 . 7 8 弱 , 縮 4 7 . 8 4 bricksA=0.1,c宏0.67,P=1984,andfbrairA=0.026, 0p=1.006,P=1.165.Theinterfaceconductivitiesduetonon 4r=8h4j=6h Figure61sothermlmesinsidethebrickforhorizonta] plainbrick uniformconductivitiesofthematerialmsidethe computationaldomainwerehandledbyharmonicmean conductivities. 4.Solution皿ethods I ︾卿 層謹鱗 評鐸 } O frominletareais≦10-3).Basedonthisprocedurethe ﹄癖 ︵つ茎︶で①﹄○︾叫︺再⑳工 00 000 000 000 5 3 35 2 25 115 A]lofthegovemmgequationsarediscretized basedoncontrolvolumeapproachonstaggeredgrid system、Inordertoavoidthephysicallyunrealisticresult duetotlmestepofthetransientproblemthefUⅡyimphcit schemeisadopted、TohandUetheconvective-djffusio、 problem,thepowerlawscheme1sused、IIhesetsof discreatizedljnearequationsaresolvedbyusmg line-by-linemethodwhichiscombinedwithThomas algorithm.’Ibcouplethepressuredistributionsand velocitiestheSIMPLEalgorithmisemployed・Iterations processwi]lbestoppedifthecontinuityequationis satisfied(totalmassresidualcomparedtototEdmass Inordertomakecleartheeffectofthebrick surface,totalheatabsorbedbyplgdnsurfacebrickand discontmuousbrick,squareandtriangledimples,are showninFigure7. 謬 鐸 』二−@ 騨軍 髭鱗 Ⅱ 2 3 4 5 6 7 time(h) FORTRANcodeshavebeendeveloped '一一Piain-←squaredimple--triangeIdimple Figure7Totalheatstoredbyplainbriokand・discontinuous brick 5.Resultm-,rldiscussion Calculationsforallcaseswithtimestep10 second,lminute,andlOminutefortota16.5hourswere caTTiedout. ltcanbeseenfromthegraphthatthebriokwith plainsurfacesstoresthehighestheatcomparedtothe discontinuousbricks、Itisknownthattheelementofthe 丘.Z・邸姥“of坊e6画cksuz造ces Theisothermlinesinsidethebrickfbr4r=1h, 2h,3h,and6hforhorizontalplainbrickareshownm brickadjacenttotheflowingairstoresbiggerheat, comparedtotheelementinsidethebricklcausedbythis elementhasthehighertemperatureduetohigherheat 539 第21回寒地技術シンポジウム(2005) transferrate・Sincetheplainbrickhasnoreducedelement thetotalheatstoredbythisbrickisthehighest、InthG discontinuousblicksbecauseoftheelementofthebIicks adjacenttothenowingairisreducedduetopresenceofthe dimples,thetotalheatstoredbythisbncksarelowerthan theplainbrick、Ontheotherwaycanbesaidthattotal heatstoredintheplainbrickisbiggercomparedtothe discontinuousbricksbecauseofthetotalmassofthebrick discontinuousbricksarelowerthantotalmassofthebrick mcase-1. rateisenhanced. 400 l 悪 j 葱 1 3 5 C 篭 軸蹴、鋸 二畳。記 2 】 0 ︵つ茎︶で。﹄。]吻型pu工 ■ も L コ カ inc]inedpositionswerecarriedoutandresultispresented inFigure9andFigurelO・ Figure9showthatfbralmostallsidesofthe brickisothermlinesareclosely.Itdoesmeanthattheheat transfbrrateshomthesesidesarebigger・Thisresult 』readye沖ected・Smcethepositionofthebrickisinc]ined towardstotheupstream,thesurfacesclosetothe acceleratednowislarger・Theacceleratednowresultsthe biggertemperaturegradientconsequentlyheattransfbr 曽蕊I頚霧 堂iミ蕊│塞副陸 3 4 05 005 00 0 0 3 20 211 ︵三︶且囚﹄ロ﹄帥匡図曽ロエ 5 02 52 01 51 050 33 k = 齢 a2’厘鈴“of坊eムr2mkpos必z,on Inordertomakecleareffectofthebrickposition tothestreamflow,calculationabrickwithhorizontaland 一陶三閏一屋 崎 塁 壁 ・ 瞳 聴 5 冠 ; fimQ(h) 垂 帯 ∼ $ 〆 ウ 〆 0 蝉、孟幽 ▲垂 今 勾 ぽ ゆ ← す 芸斗… 心〆▼ 守・ 5 : 瞳 1一一円ain−鰯一Squa”dlmpla-←TriangeldimpI‘’ ロ Figure8IIbtalheattransfbrratefi,omthe8drintothebrick 2 3 4 5 6 7 time(h) 一一Horizontal−←Inclined Thetotalheattransferratefromairintothe brickforplajnbrickanddiscontinuousbricksareshownm Figure8.Forallbricks,ase]叩ected,therateoftheheat transferisdecreasmgwithtimeelapse、Thegraphshows thatforpl8dnbrickthetotaltransfもrheatrateisthe FigurelOTotalheatstoredbybrickforhorizontalblick andinc]inedbrick b i g g e s t 、 B a s e d o n t h o s e a n a l y s i s c i t e d a b o v e 1 i t i s c l e a r t h a t FigurelOshowsthattotalheatstoredbybrick employplainsurfacebxickisbetter. fbrhorizontalandinc]medbrick・Sincetheheattransfer 1 総 l I i鱗 i 鍵 1 l識蕊識識騨蕊識鞘識舞穀蕊識識殿1 2 1 . 8 1 5 t .:ふ懲戟鍵鶏溌1 | f 錘 24.76 5 0 α=2h α=lh ratefbrinc]inedbrickisbiggerthanhorizontalbrickthe totalheatstoredforinclinedbrickisbigger、Basedonthis discussionisbettertoemploytheplainbrickwithinclined p o s i t i o n . 丘8.4ppZzba如刀”坊eheaオs”渡g巴aj7ちけeノ刀 Calculationfbrheatstoragetankconsistof24 plecesofplainbnckswithfburdjffbrentarrangements werecarxiedout・Streamlinesandtemperature distributionsinsidetheheatstorageat4r=lh6minute forallmodelsareshowninFigurellandFigurel2・Itcan beseenclearlythatinclinedandstaggeredpositionhas stronglyinmuencetotheflowcharacteristicmsidetheheat storage,Inadditiontheflowcharacteristicinfluencesthe I I l 鱗 l I i N heattrans生rratefromtheajrintothebricks・The influencecanbeseeninFigurel6.Thefigureshowsthat f b r h o r i z o n t a J i n l i n e p o s i t i o n o f t h e b r i c k s ( m o d e l l ) o u t p u t airtGmperatureishighercomparedtotheothers, I f 錘 2 7 . 9 ? 群騨韓蝉簿難癖鱈ii鍵#1 5 1 0 Ni燕鶴熱戦1鱗懲騨蝉蝿:識'蝿鶴騨:割 3 6 . 2 5 5 0 α = 8 h α = 6 h Figure91sotherm]jnesinsidetheinc]inedbrick 540 第21回寒地技術シンポジウム(2005) !》鱗:鍵灘繍鱗鱗議 砿型型 画 雲 霞 議 回= =三雲 率謹悪鱈溌 灘 蕊蕊蕊蕊 :灘驚蕊蕊 蛎認穂癖 5 0 2 0 Modell Modell 輿 一 一霧 i 灘 i 蕊I 一 1 画 :娯尋鐸諏罪浬オ lE蛍 5 0 2 0 Model2 Model2 恥班 5 0 2 0 Model8 Model8 霞§ 霊 譲灘 蕊 識 灘繍 驚灘灘蕊 、IIIII& 蕊擁盤溌嬢雛 琴海蕊踊騨I │醗理鍵鍵 5 0 2 0 Model4 鱗 Model4 F i g u r e l 2 r l l e m p e r a t u r e d i s t r i b u t i o n a t 4 r = l h 6 m i n u t e F i g u r e l l S t r e a m l i n e s a t α = l h 6 m i n u t e f b r a ] l m o d e l s foranmodels 541 第21回寒地技術シンポジウム(2005) Totalheatstoredforallmodelsaredepictedin Figurel3,Thefigureshowsthatthemodel4revealsthe bestperfbrmancetostoreheatfromthe8dr‘ 5F 一一model−1−←modcI-2ヨーmodeI-3-う←model-4 4,5 一 4 一蓋戸 一 宅 ー 一 〆、 一一塾 庫 一 一 〆〆 窒乏 一 や 一判彦一 onlyO、5N7m2becausethemletvelocityissman. 0.6 ;一学 一一Moclell-←Model−2=←ModeI-3-今←MDdel-4 古 =露 0.51 2 4 6 3 0 8 1 0 1 2 time(h) 唾。﹄ロU﹄. 0 0 0隠 4 NE、芝 z 多 ≦ 硯5 一生一 / 兎一〆 / /〆 ℃⑪﹄B叫甘①二一日。﹂ 32 52 15 塁 3 . 5 一沢 一×一 Pressuredropisanessentialprobleminthe heatstorageapplication・Inthisnumericalstudythe pressuredropmsidetheheatstoragefbranmodelalso calculatedandtheresultsarepresentedinFigurel5、The graphshowsthatpressuredropfbrmodel-4isthehigher comparedtotheothermodels、Itistruethatpressuredrop formodel-4ismorethan500%comparedtomodel-1and model-2、ButthevalueofthepressuredroPfbrmodel-4is :"i一一F. 一 Figurel3'Ibtalheatstoredfor8Jlmodels ‘;&=‐=-− 両 0 Efficiencyoftheheatstoragelsacomparisonof つ む 訂 一 日 2 6 8 1 0 1 2 time(h) theheatstoredbyallbrickswithtotalheatreleasedby WarInalr. Figurel5PressuredroPfbrallmodG1s gsr”eα妙6ric肺 7 7 ' ル ー o γ e / g a s e ム y a j r ( 1 1 ) Efficiencyfbra]lmodelsisdepictedinFigurel4、Inthe graphx-a漣sistimeinsecondandlogaxithmicsc8Jeisused inthisa五s、Thegraphshowsthatfbrallmodelsatthe timelessthanlO3secondtheefficiencyisincreasing、It doesmeantherateofheatstorageinstoringtheheatis positive,Afterlhourtheratefbrancasesarenegative、 Fromthisgraphitcanbeseenthatmodel"4reveB1sthe bestefiiciencyandstmcanstoretheheatmorethan50% until2.104second(5hours)ofstoriI1gtime,Thejudgment alsocanbemadebasedonthegraphibrthisparticular typeanddimensionitisnotrecommendedtooperatethis heatstoragetankmorethan5hoursbecausetheefficiency 6.Conclusions Theconclusionscanbemadeareasfonows: 1.Thepresenceofthedjmplesonthesurfacesofthe bricksisreducingthetotalheatstoredinsidethe b r i c k 、 2.IIhedimPlesalsoreducingtherateofheattransfbr fromtheairintothebncko 3、Inclinepositionofthebrickrevealsthebestpositionin ordertoenhancetheheajjtransfbronthebncks、 4.Apphcationtotheheatstoragetankcombinationof inclinedandstaggeredPositionoftheblicksreveals thebestperfbrmance, 5.Fortheparticularproblemthatcitedabovethestoring durationshouldnotbemorethan5hours. wmbelessthan50%. 一弓虫一一主 発 〆 』 r j r A〆 l‘S、.V:Patankar,NumericalHeatIIi,ansfbrandmuid Flow,Hemisphere,Wもlshington,DC,1980. 2A.Bejan,OonvectionHeatnansfbr,JohnWiley& 一一 ︵鋲︶語。仁②一。躯画 06 05 04 03 02 01 00 7 7.Reibrences , 坐 4 国 h A − :≦空 Sons,Inc.,NewYbrk,1995, 3.J・SuzukiIS,Matsusaka,H・Ambarita,T、YUta,M 1−←modol-1-甚一modeI-2言←model-3--mod。1-41 SatouandK、Kishinami,睡甲erlmentalStudyon ThermalPerfbrmanceofRockBedHeatStorage 1.00E十011.00E+021.00E+031.00E+041.00E+O5 tjme<sBcond) System,2005SymposlumonEnvironmental Engineermg,MuroranHo]dKaidoJaPan2005,pp Figurel4StorageE値ciencyfbrallmodels 552−555. 542
© Copyright 2024 ExpyDoc