已阅读5页,还剩3页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
PergamonComputers&FluidsVol.24,No.1,pp.55-62,1995Copyright01995ElsevierScienceLtd0045-7930(94)00020-4PrintedinGreatBritain.Allrightsreserved0045-7930/95$9.50+0.00PETROV-GALERKINFINITEELEMENTANALYSISFORADVANCINGFLOWFRONTINREACTIONINJECTIONMOLDINGNITINR.ANTURKARFordResearchLaboratory,FordMotorCompany,P.O.Box2053,MD3198,Dearborn,MI48121-2053,U.S.A.(Received4August1993;inrevisedform4May1994)Abstract-Anumericalschemeforcomputingtheadvancementofaflowfrontandrelatedvelocity,pressure,confersionandtemperaturedistributionsduringmoldfillinginreactioninjectionmolding(RIM)isdescribedinthiswork.IntheRIMprocess,theconvectivetermintheenergyequationisdominant.Therefore,thenumericalschemehasincorporatedaPetrov-Galerkinfiniteelementmethodtosuppressspuriousoscillationsandtoimproveaccuracyofthecalculations.Theotherfeatureofthenumericalschemeisthattheflowfrontlocationsarecomputedsimultaneouslywithprimaryvariablesbyusingasurfaceparameterizationtechnique.Thenumericalresultscomparewellwiththereportedexperimentaldata.ImprovedaccuracyobtainedbythisnumericalschemeintheflowfrontregionisexpectedtoassistinthepredictionsofthefiberorientationsandthebubblegrowthinRIM,whicharedeterminedprimarilybytheflowfrontregion.I.INTRODUCTIONReactioninjectionmolding(RIM)isawidelyusedprocesstomanufactureexteriorfasciasintheautomobileindustry.Inthisprocess,aprepolymerizedisocyanateandapolyol/aminemixturearemixedtogether,andinjectedintoamold,wherepolymerizationoccurs.Afountainfloweffectintheadvancingflowfrontregionduringthemold-fillingstageplaysanimportantroleindeterminingtheresidencetimeofthefluidelementsandincontrollingthefiberorientationsinthefinalproduct11.Anaccurates:imulationofthisflowfront,however,posesachallengingproblem.Evolvingflowdomainwithadvancingflowfrontrequiresupdatingofthenumericalgridsandpredictionofthemovingboundaryateverytimestep.Lowthermalconductivityofthematerial,highflowratesintheRIMprocess,andhighlyexothermicrapidreactionsresultinconvection-dominatedenergytransportequation,whichneedsaspecialnumericaltreatment.Besides,movingcontactlinesnearthewallsneedsuitableboundaryconditionsthatdonotintroducenumericalinstability.AnumericalschemethatincorporatesallthesecomplexfeaturesoftheRIMprocessisrequiredforaccuratepredictionsneartheflowfrontregion.Previousstudieseitherhavemadesimplifyingassumptionsregardingtheflowfrontregion2&l,orhavenotcomparedtheirresultswiththeexperiments5,6.Inthispaper,wedescribeanumericalschemeindetail,whichwilladdresstheabove-mentionedcomplexities,and.presentthereleventresultsthathighlightthenumericalscheme(refertoourearlierwork7forthedetaileddiscussionofthegoverningequationsandadditionalresults).Noaprioriassumptionsaremadeinthenumericalschemeregardingtheshapeofthenewfrontorthevelocitydistributionintheflowdomain.Afree-surfaceparameterizationtechniqueisused,inwhichtheshapeoftheflowfrontiscalculatedsimultaneouslywithotherfieldvariables,suchaspressure,velocitiesandconversion,byincorporatingkinematicboundaryconditionatthesurfaceoftheflowfrontasoneofthegoverningequations.AconventionalGalerkinfinite-elementtechniqueisnotoriousforitsnumericalinstabilityinconvection-dominatedtransportproblems8.Theresultingspuriousoscillationscanbeusuallyeliminatedbymeshrefinement.However,fortransientproblemdescribedhere,meshrefinementisanimpracticalandexpensivealternative.Theotheralternativesincludevariousupwindingschemes9-121,amethodofcharacteristics6,13,141,andaGalerkin/least-squarestechnique151.Althoughthe“conservative”methods,suchasmethodsofcharacteristicsandGalerkin/least-squarestechniquesaremoreaccurate,asimplePetrov-Galerkinupwindingmethodiseasierto5556NITINR.ANTURKARimplementandcosteffective,particularlyforatransientprobleminvestigatedinthiswork.Therefore,suchaschemeisimplementedherefollowingAdornatoandBrown9tosuppressnumericalinstabilitywithoutresolvingtoextremelyrefinedmeshes.ThegoverningequationsarepresentedbrieflyinSection2,andthenumericalmethodisdescribedindetailinSection3.ThetypicalresultsofthemoldfillingstageoftheRIMprocessinatwo-dimensionalrectangularplaquearepresentedinSection4.Theresultsarealsocomparedwiththereportedexperimentaldata2,andwiththenumericalresultsobtainedbyusingconventionalGalerkinfiniteelementmethod.2.GOVERNINGEQUATIONSThelumpedkineticrateexpressionforpolymerizationreactionsinRIMis16,171:ri=-A,exp(-E,/RT)Cr,(1)where,Ciistheisocyanateconcentration,Tthetemperature,Rthegas-lawconstant,mtheorderofthereaction,E,theactivationenergyofthereaction,andA,therateconstant.Theviscositydependsontheconversionandtemperature,andisexpressedintheformofCastroMacoskoviscosityfunction2,(X,T)=rl(X)-II(T)=A,exp()(iBXi,(2)whereXistheisocyanateconversion,X,thegelconversion,andA,E,AandBaretheconstants.Forconstantthermalpropertiesanddensityofthereactivemixture,andfornegligiblemoleculardiffusion,thedimensionlessgoverningequationsare,continuityequation:v.v=o;(3)conservationofmomentumequation:Re$+v.Vv=-pV.I+v:(rcj);Gz7,$+v-VX=Dak.(l-X)“;molebalanceequation:(4)(5)conservationofenergyequation:Gzg+vVT=V*T+Brrc(j:Vv)+Darc,(l-X)m;L.1(6)where,visthevelocityvector,qtherate-of-straintensor,tthetime,pthepressure,andk,isthedimensionlessrateconstant,definedasexp(-E,/R)(l/T-l/T,).TheequationsaremadedimensionlessusingtheaveragevelocityV,halfofthethicknessofthemoldH,andthetemperatureT,andtheviscosityqO(=r(X=0,T=T,)attheinletofthemold.AllthedimensionlessgroupsandtheirdefinitionsarelistedinTable1.Theboundaryconditionsintermsofdimensionlessvariablesare1.atthewalls:v,=0(no-slip),T=T,;2.atthemid-plane:aTjay=0,&Jay=0,V,=0;3.attheinlet:v=fullydevelopedflow,T=1,X=X,;4.atthecontactline:n*(-PI+2)=0(full-slip):5.attheflowfront:n.(-PI+2)=0(forcebalance),n.(v-ah/at)=0(kinematiccondition);Table1.Dimensionlessgroupsingoverningequations,whereAH,istheheatofreaction,AT,theadiabatictemperaturerise,andC,theinitialconcentrationofisocyanateGZGraetznumberVHpC,lkReReynoldsnumberHVlrloKviscosityratio41%BrBrinkmannnumbertoV=lkT,DaDamkohlernumber(AH,H*C$/kT,)A,exp(-E,/RT)TadbadiabatictemperatureriseAT,IT,Flowfrontadvancementinreactioninjectionmolding57wherea.,andvYarethecomponentsofthevelocityvectorv,IItheunitnormalvector,rtheextrastresstensor,hthelocationvectoroftheflowfrontandTwal,thedimensionlesstemperatureatthemoldwall.Thedetailsofincorporatingtheboundaryconditionsinthenumericalanalysisareexplainedinthenextsection.3.NUMERICALANALYSISInthefiniteelementformulationtheunknownvelocities,temperatureandconversionareexpandedintermsofthebiquadraticbasisfunctions4,thepressureintermsofthebilinearbasisfunctionsll/iandtheflowfrontshapehintermsofthequadraticbasisfunctions:(7)wherelandqarethecoordinatesinisoparametrictransformation,definedasi=1i=lintheisoparametricdomain(-14+1,-1qL(5),where,PeisthelocalelementPecletnumber(=VA/D),Atheelementsizeandc,(s)isthecubicpolynomial(=(5/8)5(-l)(t+1).Theindexi=1correspondstothevertexnodes,andi=258NITINR.ANTIJRKARcorrespondstothecentroidnodesintheelement.Thestandardone-dimensionalconvec-tiondiffusionproblemhasexactsolutionatthenodesif25,9c(Pe)=2tanh(Pe/2)l+(3/Pe)coth(Pe/4)-(X/Pe)-coth(Pe/4),(1la)c2=(16/Pe)-4coth(Pe/4).(1lb)Inatwo-dimensionalproblem,thetensorialproductofequations(10)and(11)providesthefunctioncintheweightingfunctionsdescribedinequation(9).ThelocalPecletnumberiscomputedforeachthree-nodegroupbasedontheaveragevelocitiesattherelevantboundariesinthetwo-dimensionalelement9.Therearesixsuchgroups(threeinthex-direction,andthreeinthey-direction)andthus,thereare12upwindingparametersE.ThecalculationsofthePecletnumberinvolvelineardistances,whichessentiallyneglectthecurvilinearsidesoftheelements.However,itisagoodapproximationsinceflowfrontisnotseverelydeformedinourproblem.Thediffusivitiesarel/GzfortheenergyequationandisK/Rforthemomentumequation.ThePetrov-Galerkinweightedresidualequationsare,-R:=(V.v)$dl=O,s-RL=IvReg+v.VvWfdV(12)+y-PI+(K+)VWidV-ssn.-pI+(lcf)WdS=O,(13)s-Brrc(j:Vv)-Dak,(l-X)”WdV1+sVT.VWdV-s(n.VT)WdS=O,(15)VS-RI=sn.(v-ah/&)4(+=1)dS=0.(16)swhere,VistheflowdomainandStheflowboundary.Theboundarytermsappearintheenergyandmomentumequationsbecausedivergencetheoremisappliedtothehigher-orderterms.TheresidualsR,R,R,R,andR,correspondtothevariablesp,v,X,Tandh,respectively.ThePetrov-Galerkinweightingfunctionsareusedonlyformomentumandenergyequationsduetothepresenceofconvectiontermsintheseequations.Beforeintegratingtheaboveequationsusinganine-pointGaussianquadrature,theequationsaremappedintheisoparametricdomain(referto26fordetails)andtheboundaryconditionsareapplied.TheessentialboundaryconditionsforvandTatthewalls;forv,TandXattheinletofthemold;andforv
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 上海市宝山中学2023-2024学年高三第四次模拟考试化学试卷含解析
- 内蒙古一机集团第一中学2023-2024学年高考冲刺模拟化学试题含解析
- 教师评职称个人的工作总结
- 七年级数学下学期期末达标测试题课件-(新版)浙教版
- 幼儿园读书笔记15篇
- 穿越中国走进青海海西
- 《寄东鲁二稚子》原文翻译及赏析
- 刍议家园合作新思路
- 健康在你手上指甲透露身体小秘密
- 母亲母亲与母亲相处的日子随笔
- 2024版N1叉车司机模拟测试题库含答案 全套
- 行业现状报告组织结构优化与运营效率提升
- 胸腺恶性肿瘤课件查房
- 古建筑木结构防火保护技术培训
- 厂房空调调试方案(大型厂房空调配套方案)
- 通信行业深度洞察报告
- 临床营养诊疗流程的建立与发展
- 包场协议书正规范本(通用版)
- 学校窗帘施工方案
- 2023年新劳动法全文
- 2023年全国各地中考试题分类汇编:文言文阅读(94份)
评论
0/150
提交评论