33 5 Vol.33,No JournalofHebeiUniversityofScienceandTechnology Oct.2012 : (2012) /,, ( 河北科技大学机械工程学院, 河北石家庄 ) : 利用计算流体

Size: px

Start display at page:

Download "33 5 Vol.33,No JournalofHebeiUniversityofScienceandTechnology Oct.2012 : (2012) /,, ( 河北科技大学机械工程学院, 河北石家庄 ) : 利用计算流体"

过壑夏
5 years ago
Views:

1 33 5 Vol.33,No JournalofHebeiUniversityofScienceandTechnology Oct.2012 : (2012) /,, ( 河北科技大学机械工程学院, 河北石家庄 ) : 利用计算流体力学软件 FLUENT 的凝固 / 熔化模型, 对双层壁圆筒内填充的铝硅合金相变材料的熔化 / 凝固过程进行了数值模拟, 得到了在第 3 类边界条件下圆筒内相变材料在凝固 / 熔化过程中的动态温度场分布相界面移动规律及凝固 / 熔化时间等, 对蓄热装置的设计及实验研究提供了重要的参考价值 : 相变材料 ; 熔化 ; 凝固 ;FLUENT :TK02;TK124 :A Numericalsimulationonmeltingandsolidification processofaluminum-siliconaloy CUIHai-ting,ZHANG Gai,JIANGJing-zhi (ColegeofMechanicalEngineering,HebeiUniversityofScienceandTechnology,ShijiazhuangHebei050018,China) Abstract: Byusingthesolidification /meltingmodeloffluentsoftware,themeltingandsolidificationprocessofthealumi- num-siliconaloyinacylinderheatstoragewithdoublewalweresimulatednumericaly.consideringthethirdkindofboundary condition,thedynamictemperaturefielddistributionandthepositionofthephaseinterfacechangemovingwithtimeareob- tained.andthemeltingandsolidificationdurationtimeisalsocalculated.theresultshaveimportantreferencevalueforthe numericalsimulationofthephasechangeproblemandthedesignofphasechangeenergystoragedevice. Keywords: phasechangematerial ;melting;solidification;fluent [1],,, ;, 45%~55%;,, [2-4] FLUENT, /, , H=600mm, Φ=180mm, : ; : ; : : ( ); (ZH ); ( A) : (1964-),,,,

2 Δ Δ Δ , 20, 10mm;, 4mm, 1,, :1), ;2) ; 3) 2 [5-6] : 20 图 1 储热装置示意图 Fig.1 Schematicofheatstoragedevice 20 3 PCM( ), 1 表 1 铝硅合金相变材料物性参数 Tab.1 PhysicalpropertiesofAl-Sialoy PCM /(kg m -3 ) /(kj kg -1 K -1 ) /(W m -1 K -1 ) / /(kj kg -1 ) [7] FLUENT / β= ( ρ H) t + ( ρ 珗 vh)= (κ H =h+δh, T h =h ref+ c p dt, T ref ΔH =βl, β=0ift < T, β=1ift > T, T -T T -T T)+S, ift < T < T : ρ ;t ;H ;h ref ( );h ;ΔH ;L ; β ;c p ;κ ;T 3 FLUENT 3.1 FLUENT GAMBIT,FLUENT GAMBIT Quad, GAMBIT FLUENT, 2D ( ), Solidification & Melting [8] 3 80 W/ (m 2 K) SIMPLE

$3 TemperaturechangeonPCMareainmeltingprocess Fig.4 LiquidfractionchangeonPCMareainmeltingprocess 5 6 35000s 5,,,, 6,, 35000s,PCM, 7 8 421000s, PCM 8, PCM 1.$

3 5 / 455, PCM PCM patch, 298 K, 923K 10~ 30s 图 2 Fig.2 蓄热装置网格划分 Meshofheatstoragedevice 3 PCM, 4 PCM PCM 17000s PCM, PCM ;, PCM 1.0, PCM ;, 图 3 熔化过程 PCM 区平均温度曲线图 4 熔化过程 PCM 区液相率曲线 Fig.3 TemperaturechangeonPCMareainmeltingprocess Fig.4 LiquidfractionchangeonPCMareainmeltingprocess s 5,,,, 6,, 35000s,PCM, s, PCM 8, PCM 1.0, s s, 243min 7, PCM,,PCM,PCM

$456 2012 图 5 35000s 时温度分布图 Fig.5 Temperaturedistributionat35000s 图 6 35000s 时液相率分布图 Fig.6 Liquidfractiondistributionat35000s 图 7 421000s 时温度分布图 Fig.$

4 图 s 时温度分布图 Fig.5 Temperaturedistributionat35000s 图 s 时液相率分布图 Fig.6 Liquidfractiondistributionat35000s 图 s 时温度分布图 Fig.7 Temperaturedistributionat421000s 图 s 时液相率分布图 Fig.8 Liquidfractiondistributionat421000s 4.2, 923K (673K), patch PCM 923K 673K 9 PCM, 10 PCM 9 10,PCM, PCM PCM 0,,, s PCM, 12, 5000s,PCM

$10 LiquidfractionchangeonPCMarea insolidificationprocess 图 11 5000s 时温度分布图 Fig.$

5 5 / 457 图 9 凝固过程 PCM 区平均温度曲线图 Fig.9 TemperaturechangeonPCMareainsolidificationprocess 图 10 凝固过程 PCM 区液相率曲线图 Fig.10 LiquidfractionchangeonPCMarea insolidificationprocess 图 s 时温度分布图 Fig.11 Temperaturedistributionat5000s 图 s 时液相率分布图 Fig.12 Liquidfractiondistributionat5000s s, PCM 14 图 s 时温度分布图 Fig.13 Temperaturedistributionat10000s 图 s 时液相率分布图 Fig.14 Liquidfractiondistributionat10000s

6 s 9000s 117min PCM, 5 FLUENT GAMBIT, FLUENT / FLUENT, : [1],,,. [J]. (SolarEnergy),2009(6): ca),2006,27(1): [3],,,.Al-Si [J]. (ChineseJournalofMaterialsResearch),2006,20(2): [4],,,. [J]. (MaterialsResearchandApplication),2012, 4(4): [5],,,. (JournalofHebeiUniversityofScience andtechnology),2010,31(2): [6],,. FLUENT [J]. (Low TemperatureandSpecialtyGases),2006,24(2): [7],,.FLUENT [M]. :,2007. [8],,.FLUENT [M]. :,2004. gy),2011,40(1): [10],,,. [J]. (JournalofYanshan University), 2012,36(2): [11],,,. / [J]. (HebeiJournalofIndustrial ScienceandTechnology),2005,22(5): [2],,,. Al-Si [J]. (ActaEnergiaeSolarisSini- [9],,. (JournalofHebeiUniversityofTechnolo- [12],,. [J]. (HebeiJournalofIndustrialScienceandTech- nology),2011,28(1): 1-4. 櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴 ( 383 ) : [1] WANG Hai-yan.Positiveperiodicsolutionsforfunctionaldiferentialequations[J].JDiferentialEquations,2004,202(26): [2] CHUJ,TORRESPJ,ZHANG M.Periodicsolutionsofsecondordernon-autonomoussingulardynamicalsystems[J].JDiferential Equations,2007,239(1): tions,2005,218(1): [4] DONGShi-jie,GE Wei-gao.Positivesolutionsforquasilinearsecondorderdiferentialequation[J].ApplicableAnalysis,2005,84(12): [5],. p-laplacian (JournalofHebeiUniversityofScience andtechnology),2010,31(5): [6],,,. (JournalofHebeiUniversityofScience andtechnology),2012,33(3): [3] YUJ,GUO E.Multiplicityresultsforperiodicsolutionstodelaydiferentialequationsviacriticalpointtheory[J].JDiferentialEqua- [7],,. (JournalofHebeiUniversityofSci- enceandtechnology),2012,33(4): [8] MA Ru-yun,MA Hui-li.Positivesolutionsofnonlineardiscreteperiodicboundaryvalueproblems[J].ComputMathAppl,2010,59(1): [9] ATICIF M,GUSEINOV GS.Positiveperiodicsolutionsfornonlineardiferenceequationswithperiodiccoeficients[J].JMathAnal Appl,1999,232: [10] RABINOWITZP H.Someglobalresultsfornonlineareigenvalueproblems[J].JFunctAnal,1971,7:

<313031A4C9BEC7C160BA5DB3E62831303130383135A457BAF4A4BDA769AAA9292E584C53>

<313031A4C9BEC7C160BA5DB3E62831303130383135A457BAF4A4BDA769AAA9292E584C53> 機械三甲 01 811001 王振祥國立高雄應用科技大學模具工程系甄選入學嘉義縣縣立水上國中機械三甲 02 811002 王紹誠弘光科技大學生物醫學工程系登記分發嘉義縣縣立水上國中機械三甲 03 811003 江彥廷中臺科技大學牙體技術暨