资源描述：
《化学液相沉积法制备碳_碳复合材料时的温度场研究》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
------------------------------------------------------------------------------------------------化学液相沉积法制备碳_碳复合材料时的温度场研究第42卷 第3期西 安 交 通 大 学 学 报2008年3月JOURNALOFXI′ANJIAOTONGUNIVERSITYVol.42 №3Mar.2008化学液相沉积法制备碳/碳复合材料时的温度场研究潘仡文1,杨冬1,刘红林2(1.西安交通大学动力工程多相流国家重点实验室,710049,西安所,710025,西安)摘要:、预制体对煤油,的数学模型.的预制体在汽化温度为453115K时的制备过程.给出了不同时刻的温度场分布、径向温度随,.将计算所得的电加热功率与制备过程中实际的电加热功率进行了比较,二者较为接近,验证了所建数学模型的可靠性.关键词:化学液相沉积法;碳/碳复合材料;温度场中图分类号:TQ028 文献标志码:A 文章编号:0253Ο987X(2008)03Ο0309Ο04ResearchonTemperatureFieldofKerosenewhenPreparingCarbon/CarbonCompositesbyChemicalLiquidDepositionPANYiwen1,YANGDong1,LIUHonglin2(1.StateKeyLaboratoryofMultiphaseFlowinPowerEngineering,Xi′anJiaotongUniversity,Xi′an710049,China;2.The43rdInstituteofthe4thAcademyofCASE,Xi′an710025,China)—————————————————————————————————————— ------------------------------------------------------------------------------------------------Abstract:Amodeloftemperaturefieldandphasechangeinthecarbonfiberpreformwasestab2lishedbasedonheatandmasstransferofgasandfluidphasesflow,heatconductionbetweenthepreformandkerosene,andheatandmasstransferinthecarbonfiberperform.Anumericalin2vestigationonpreparationatthevaporizationtemperature453115Kforthepreformwithadiame2terof70mmusinginductionheatingwasperformed.Thechangeofthetemperaturefieldwithtimeinthepreformwasobtained.Thetemperaturedistributionatdifferenttime,thevariationofradicaltemperaturewithtime,andvaporratesalongtheradicaldirectioninthephasechangeprocesswereindicated.Thecomputingresultsareinagreementwiththedataoftheactualpowerinputofelectricheating.Keywords:chemicalliquiddeposition;carbon/carboncomposites;temperaturefield 碳/碳复合材料广泛应用于航空、航天、生物及多种民用领域.目前,制备碳/碳复合材料的方法大体分为气相法和液相法两大类.气相法主要有化学气相沉积(CVD)和化学气相渗透(CVI).气相法出现较早,无论是在实验还是数值模拟方面都较为成熟[1Ο4].化学液相沉积法是在液态碳氢化合物中进行的,在预制体的周围不断形成高密度的液Ο气相包裹层,这种包裹层的密度比传统气相沉积法的纯气相层高出50～—————————————————————————————————————— ------------------------------------------------------------------------------------------------100倍,从而预制体的致密化速度比气相沉积法也高出50～100倍,大大缩短了制备时间,降低了制备成本.基于上述优点,本文对碳纤维预制体骨架浸在煤油中发生池沸腾换热时的骨架温度场进行了数值计算,计算时考虑了气液两相间的传热与传质以及相变的产生.收稿日期:2007Ο06Ο18. 作者简介:潘仡文(1982-),男,硕士生;杨冬(联系人),男,副教授. 基金项目:国家科技支撑计划资助项目(2006BAA03B02Ο03).310西 安 交 通 大 学 学 报第42卷 1 数学模型图1所示为化学液相沉积法(CLD)制备碳/碳复合材料的示意图,石墨棒被交变磁场的磁力线所切割,产生涡流,涡流回路使石墨棒发热,从而加热预制体.采用液态煤油作为前躯体,将碳纤维预制体浸泡其中,在高温作用下含碳前躯体分解,裂解产物中的氢离开体系而碳沉积于预制件上.随着反应进行,整个预制体最终达到预期的致密程度.,较为恒定的高温,始致密化,.由于热梯度的原因,(芯部)沉积,沉积带逐渐由内向外推进,直至预制体内外热梯度消失,整个预制体达到致密化.本文在计算中作了如下假设.(1)因为预制体紧密包裹着加热棒,接触热阻可以忽略,从而近似认为模型中与石墨加热棒接触的表面具有和石墨加热棒相同的温度.(2)边界上的换热为池沸腾换热.—————————————————————————————————————— ------------------------------------------------------------------------------------------------(3)预制体的孔隙率不在模型中体现,而是在后面的计算中考虑进去.(4)煤油的热裂解是一个极为复杂的化学反应过程,且煤油裂解时消耗的反应热与整个研究体系所消耗的总热量相比可以忽略不计,故在计算过程中忽略了煤油裂解时消耗的反应热.,并且假.实际问题中液相煤油与气相煤油为互相贯通的连续体,与混合模型吻合较好,因此选择混合模型来模拟气液两相问题.对于混合模型,其连续方程为(ρ)?ρ+???(mmvm)=mtn(1)式中:vm=nk=1-1αρs;kkvk/ρm为质量平均速度,kg?∑ρm=k=1—————————————————————————————————————— ------------------------------------------------------------------------------------------------αρ为混合密度,kg?m∑kk-3;αk为第k相的浓度;m??为气穴或质量源的质量传递,kg.能量方程为(αρkkEk)+???∑(αkvk(ρkEk+p))=tk∑=1k=1???(keff??T)+SE式中:keff为有效热传导率,W?m-2?K-1.动量方程为(ρ)vm+???(ρ??p+ρmvmvm)=-mg+F+t???[μm(??vm+??v)]+??Tmnk=1nn(2)αρvv∑kk—————————————————————————————————————— ------------------------------------------------------------------------------------------------rkrk(3)式中:n为相数;F为体积力,N;μm为混合黏性,kg?(m?s)-1;vrk为第k相(k=2)的飘移速度m?s-1;p为体系所受压强,Pa;g为重力加速度,图1 碳/碳复合材料制备示意图m?s-2.11112 多孔介质控制方程 求解计算软件中有特基于上述假设与实际情况,建立了用于数值模拟的模型,大圆柱中间的空心小圆柱为石墨加热棒,大圆柱为预制体骨架,大圆柱外围被液体煤油所包裹.对模型进行网格划分,划分示意图如图2所示,每个网格是一个计算单元.模型划分采用结构化网格,网格为三维的六面体网格,网格数为41800.111 控制方程组1.1.1 两相流方程 混合模型是一种简化的多相为本文研究设计的多孔介质模型,其动量方程是在式(1)的基础上加上一个动量源项33Si=-—————————————————————————————————————— ------------------------------------------------------------------------------------------------j=1∑νj+Dijμj=1∑Cijννmagj2(4)式中:第1项为黏性损失项;第2项为惯性损失项;Si是第i个(x、y或z方向)动量方程中的源项;Dij是以1/α为对角单元的对角矩阵,α为描述多孔介质渗透性的参数;Cij是以C2为对角单元的对角矩 第3期 潘仡文,等:化学液相沉积法制备碳/碳复合材料时的温度场研究311阵,C2是惯性阻力因子. 能量方程为(γρ)ρ(ν(ρfEf+(1-γsEs)+???fEf+p))=t??keff??T-hJ∑ii—————————————————————————————————————— ------------------------------------------------------------------------------------------------i113 初始条件τ?ν)+(+Sfh(5)整个计算为非稳态过程,操作压强为121159kPa,反应初始温度为室温.石墨棒初始温度为室温,用自定义函数控制其升温过程,使其在之后的10min内以指数规律上升到1273115K,然后保持温度不变.大圆柱的外侧面设置为压力出口,上下端,.,,孔,终止时的孔隙率为10%.用自定义函数来控制相变过程以及孔隙率的变化.ρ式中:f为流体密度;ρs为固体密度;p为孔隙内的压强;τ为曲折因子;ν为流体速度;Ji为i方向的摩尔通量密度;hi为i方向的溶解度系数;Ef为流体总能;Es为固体介质总能;γ参数;Sfh为流体焓的源项.112 边界条件压力为011,293115K时的密度ρm-3.在本计算研究的煤油20=81816?温度和密度范围内,煤油的密度随温度的变化呈线ρ=ρ(t-20),β为平均密度系数,β=性关系:20-β0176kg?m-3?℃-1[5].本文计算中,煤油在达到沸2 计算结果与分析—————————————————————————————————————— ------------------------------------------------------------------------------------------------所建模型的轴向长度为0115m,径向长度为0107m,石墨棒的半径为0103m.通过计算发现,随着时间的推移,整个预制体骨架的温度场也随之发生变化.选取4个比较有代表性的时间来显示整个骨架的温度分布,如图3～图6所示.点以后,始终处于饱和状态.压力为011MPa时,其饱和温度为453115K,液相密度ρ=697kg?m-3.煤油是多种化合物的复杂混合物,计算煤油的平均相对分子质量的公式为M=a+bt+ct2(6)式中:系数a、b、c与特性因子K的关系见表1.表1 参数a、b、c与特性因子K的关系Kab10.0560123010155701240111059012401115630122512106901180c01008000100090010010001001150100140 K是由平均沸点和相对密度所决定的,取K=111362[6],—————————————————————————————————————— ------------------------------------------------------------------------------------------------经插值计算得M=139106.煤油的导热系数为λ=011008(1-0100054)/d204式中:d20时煤油和水的密度之比.计算得到4是20℃煤油的导热系数为011231W?(m?K)-1.Cooper提出的以相对压力为基础的换热系数图3 50s时骨架的温度分布计算公式为[7]0112-014343lnRp(-014343lnpr)-0155M-015q0167h=55pr(7)式中:pr=p/pcrit=112/3018=0103896,pcrit为液体煤油的临界压力;M=139106;加热面的表面粗糙度Rp=112×10-6图4 300s时骨架的温度分布从图5中可以看出,反应进行到600s时轴向温度分布不均匀,因为反应中液相煤油汽化时会带走一部分热量,同时形成气液两相流动.气相煤油一部分被排走,而一部分液相煤油补充进来,形成回;热流密度q=1320179kW?m-2.根据式(7)计算可得碳纤维表面沸腾换热系数—————————————————————————————————————— ------------------------------------------------------------------------------------------------h=706219W?m-2?K-1.312西 安 交 通 大 学 学 报第42卷 图5 600s积的乘积,计算得到的电加热功率为28191kW,与实验中的实际电加热功率2910kW相比,误差较小,表现出较好的一致性.3 结 论(1)通过模拟计算,得到了不同时刻的温度场分布、径向温度随时间的变化规律,以及相变过程含气率沿径向的分布.从径向温度随时间的变化规律可以看出,随着反应的进行,形成了内高外低的温度梯度.研究结果为优化工艺过程、改善产品质量、降低生产成本提供了一定的参考价值.(2)根据以上所得的计算结果,不难看出由于和加热棒比较接近的内层温度高,故内层的煤油首先由液相转变为气相,然后发生裂解反应,生成的碳沉积下来.沉积过程由里向外进行,最终生成致密而结实的碳/碳复合材料.参考文献:[1] 侯向辉,李贺军,李克智,等.CVI制备碳/碳复合材料图6 900s时骨架的温度分布流.另外3个时刻的回流不明显,汽化所带走的热量很少,故轴向温度分布较图5所显示的时刻均匀.4个时刻沿着半径方向的温度分布如图7所示,气相体积分数φg的分布如图8所示.在内热源的作用下,整个预制体的温度从里向外逐渐降低,—————————————————————————————————————— ------------------------------------------------------------------------------------------------形成了一个内高外低的温度梯度.石墨棒升温至1273115K后保持温度不变,成为一个恒温内热源.在其作用下,温度梯度逐渐减小.同时,处于沸腾状态的液相煤油逐渐发生相变,由于温度梯度的作用,会得到如图8所示的预制体的气相体积分数的分布情况.可以明显看出,在反应进行到900s时,预制体内的液相煤油已基本转化为气相煤油.致密化行为模拟研究[J].兵器材料科学与工程,1999,22(2):28Ο33.HOUXianghui,LIHejun,LIKezhi,etal.Numericalstudyondensificationinpreparingcarbon/carboncom2positesbyCVI[J].OrdnanceMaterialScienceandEn2gineering,1999,22(2):28Ο33.[2] 姜开宇,李贺军,侯向辉,等.轴对称C/C复合材料件等温CVI过程的数值模拟研究[J].西北工业大学学报,2000,18(4):665Ο668.JIANGKaiyu,LIHejun,HOUXianghui,etal.Nu2mericalstudyonaxisymmetriccarbon/carboncompos2itesbyconstanttemperatureCVI[J].JournalofNorthwesternPolytechnicalUniversity,2000,18(4):665Ο668.图7 不同时刻沿径向的温度分布(下转)384西 安 交 通 大 学 学 报第42卷 —————————————————————————————————————— ------------------------------------------------------------------------------------------------参考文献:[1] DOMENK,KUDOA,OHNISHIT.MechanismofphotocatalyticdecompositionofwaterintoH2andO2overNiOΟSrTiO3[J].JCatal,1986,102(1):92Ο98.[2] KATOH,ASAKURAK,KUDOA.HighlyefficientwatersplittingintoH2andO2overlanthanum2dopedNaTaO3photocatalystswithhighcrystallinityandsur2facenanostructure[J].JAmChemSoc,2003,125(10):3082Ο3089.[3] ZOUZhigang,YEJinhua,SAYAMAK,etal.splittingofwaterundervisiblewithoxidesemiconductor[J,414(6864):[4] HITOKIG,TAT,KONDOJ,etal.Anoxynitride,TaON,asanefficientwateroxidationpho2[9] YOUNHC,BARALS,FENDLERJ.Dihexadecylphosphate,vesicle2stabilizedandinsitugeneratedmixedCdSandZnSsemiconductorparticles:prepara2tionandutilizationforphotosensitizedchargesepara2tionandhydrogengeneration[J].JPhysChem,1988,92(22):6320Ο6327.[10]XINGChanjuan,ZHANG,YANWei,etal.Band2solutionofCd1-xZnxSbywatersplit2IntEnergy,2006,31(14):2018Ο[11]TSUJII,KATOH,KOBAYASHIH,etal.Photo2—————————————————————————————————————— ------------------------------------------------------------------------------------------------catalyticH2evolutionreactionfromaqueoussolutionsoverbandstructure2controlled(AgIn)xZn2(1-x)S2solidsolutionphotocatalystswithvisible2lightresponseandtheirsurfacenanostructures[J].JAmChemSoc,2004,126(41):13406Ο13413.[12]TSUJII,KATOH,KOBAYASHIH,etal.Photo2catalyticH2evolutionundervisible2lightIrradiationoverband2structure2controlled(CuIn)xZn2(1-7329.[13]TSUJII,KATOH,KUDOA.Visible2light2inducedH2evolutionfromanaqueoussolutioncontainingsul2fideandsulfiteoveraZnS2uInS22AgInS2solid2solutionphotocatalyst[J].AngewChem,2005,44(23):3565Ο3568.[14]MAEDAK,TERAMURAK,LUDaling,etal.Pho2tocatalystreleasinghydrogenfromwater[J].Nature,2006,440(7082):295.x)λ≤500nm)tocatalystundervisiblelightirradiation([J].ChemCommun,2002(16):1698Ο1699.—————————————————————————————————————— ------------------------------------------------------------------------------------------------[5] NIISHIROR,KATOH,KUDOA.Nickelandeithertantalumorniobium2codopedTiO2andSrTiO3.photo2catalystswithvisible2lightresponseforH2orO2evo2lutionfromaqueoussolutions[J].PhysChemChemPhys,2005,7(10):2241Ο2245.[6] KUDOA,SEKIZAWAM.PhotocatalyticH2evolu2tionundervisiblelightirradiationonZn1-xCuxSsolidsolution[J].CatalLett,1999,58(4):241Ο243.[7] KUDOA,SEKIZAWAM.PhotocatalyticH2evolu2tionundervisiblelightirradiationonNi2dopedZnS[J].ChemCommun,2000(15):1371Ο1372.[8] REBERJF,RUSEKM.Phtotchemicalhydrogenpro2ductionwithplatinizedsuspensionsofcadmiumsulfideandcadmiumzincsulfidemodifiedbysilversulfide[J].JPhysChem,1986,90(5):824Ο834.S2solidsolutions[J].JPhysChem:B,2005,109(15):7323Ο(编辑 荆树蓉 赵大良)(上接)[3] VAIDYARAMANS,LACKEYWJ,AGRAWALP—————————————————————————————————————— ------------------------------------------------------------------------------------------------K,etal.1Dmodelforforcedflow2thermalgradientchemicalvaporinfiltrationprocessforcarbon/carboncomposites[J].Carbon,1996,34(9):1123Ο1133.[4] GOLECKII,MORRISRC,NARASIMHAND.Rapiddensificationofporouscarbon/carboncompositesbythermalgradientCVI[J].AppPhysLett,1995,66(18):2334Ο2336.[5] 胡志宏.火箭发动机特种煤油传热特性的试验研究[D].西安:西安交通大学能源与动力工程学院,1996.[6] 梁文杰.石油化学[M].东营:中国石油大学(华东)出版社,1996:96Ο104.[7] 杨世铭,陶文铨.传热学[M].4版.北京:高等教育出版社,1998:221Ο227.(编辑 荆树蓉 赵大良)——————————————————————————————————————
