硕士论文-汽车空调系统与发动机冷却系统的耦合分析

上海交通大学硕士学位论文汽车空调系统与发动机冷却系统的耦合分析姓名：祁照岗申请学位级别：硕士专业：制冷与低温工程指导教师：陈江平20041201COP(1114)Theairconditioningsystemandenginecoolingsystemareinter-actedinvehicleoperation.Thesystemperformancecanbeimprovedbycoupledanalysis.Numericalmodelsaredevelopedforairconditioningsystemandenginecoolingsysteminthispaper.Thepressuredropanalysisofairsideofheatexchangers(HX)wasstudieddetailly,andthedatabasewassetuptoimplementthenumericalmodels.Parametersincludingairdistributionflowresistancevehiclebodyandsolarradiationweremodeled.ExperimentsarecarriedoutforairconditioningsystemandenginecoolingsystemperformancesintheEnvironmentalSimulationLab.Coupledcalculationalresultspresentthemodeldevelopedinthispapercanwellsimulatethesystemperformanceandcompartmentcool-downcurve.Numericalandexperimentaldatahaveagoodagreement.Themodeldevelopedinthispaperishelpfulforthedesignoftheautomobileairconditioningsystemandenginecoolingsystem.TheresultsofrelativitystudyshowthatairconditioningsystemcoolingcapacityandCOPwillretainsteadyandkeeppassengerscompartmentcomfortableunderawidevehiclespeedrange.Buttheheatoutputofradiatorchangesirregularly.Underdifferentnormalcondenserairoutlettemperature(notallowedtobe11to14abovetheambienttemperature),heatoutputofradiatorchangeslittlebutthecoolantandoiltemperaturecanberetaininatemperaturerangebetweenwhichrangeenginecouldoperateverywell.Keywords:airconditioningsystem,coolingsystem,coupledsimulation,experimentalanalysis,relativitystudyiiiairsideA−mrefsideA−mRAmAm2gam/scJ/kgKdmadmgm/s2Hmhm0hJ/kgkmLmMkg/kmolmkgdryairmkg/smkg/sn1/scriticalpbarp∆PaqW/m2QwaRmsmtKvm/saVm3hVm3Vm3/sxx/vapormmxx/waterdryairmmα(w/m2K)βρkg/m2sρkg/m3σN/mkw/m2KηPasVηisentηλW/mKiv12brcondirdirodirwIindindoindwmormotro/oflowpanr[1]CFC-12CFC-121974RowlandMoin1987CFC-12ODPGWPCFCsHCFCs2030CFCsHFC-134aCFC-121992HFC-134aCFC-12HFC-134aCFC-12DuPontBateman[2]HFC-134aCFC-12CFC-12HFC-134a6%4%HFC-134aCFC-12175Kpa8NippondensoHirata[3]CFC-12HFC-134aHFC-134aHFC-134aCFC-12HFC-134aCFC-12HFC-134aCFC-127NissanAkioYamada[4]HFC-134aCFC-12HFC-134a1020[5]HFC-134aHFC-134aCalsonicRefkiEl-BouriniHFC-134aCFC123HFC-134aCFC-12WangTouber[6]WangTouberNyersStoyan[7]106X.Jia[8][9,10][11,12]LeeKyuHyunWonJong-Pil[13][14]4[15]CFC-12HFC-134a[38][16]Hamilton[17]RamanAli[18]5[19]SANTANA[20],HFC-134a[40][58]SaizJabardo[21]20%10%LeeG.H.YooJ.Y.[22][23]20003000K6Fig.1-1Enginecoolingsystem(1-1)1-1[40]7[41]CAD[42],[43][44][45]MATLAB/SIMULINK8Fig.1-2coolingsystemandairconditioningsystem(1-2)9[23]CFD(ComputationalFluidDynamics)()()CAE(Computer-AidedEngineering)CFDCFDCFDCFDCFD12103411CCOTCCTXV2.1.1.1[24]Lp∆rp∆ap∆sp∆2.1.1.2rp∆0()Lrrdppdldl∆=∫(2-1)112(1)1/7xFrxFrρβρ′=≤′′−+Fr2()[1(1)][1(1)(1)]2rrdpxxKdldρρρξρρρ′′=+−−−−′′′′′(2-2)12Re2300ZPξ1/2.512log()3.7ReZPdκξξ=−+(2-3)Re(1(1/))ZPdxρηηη=′′′′−−(2-4)Re2300ZP64/ReZPξ=ReZP1121/7FrFrβ+22rdpxdldρξρ⎛⎞=Φ⎜⎟′′⎝⎠(2-5)Φ2.1.1.301(01)11()apρρρ−∆=−′′′(2-6)ap∆''''22211124[(1)(1)]amxwxwxwxwpdπ−+−−−∆=(2-7)ws0.02050.016()()wdswρρρη−′′′==′′′′(2-8)2.1.1.40sp∆=132.1.1.5CFD[46-56]CFD(2-1)CFD2-D[57]2-1Fig.2-12-DcomputationaldomainFp:finpitch,Fd:theflowdepthofthefin,Lp:louverpitch,:louverangle.FpmmFdmmLpmm11.4191.12721.121.81.1272.1.1.6Fig.2-2PressuredropcomparisonbetweenCFDresultsandtheexperimentaldata(StyleI)14Fig.2-3PressuredropcomparisonbetweenCFDresultsandtheexperimentaldata(StyleII)2-22-3CFD10%30%1200Re≤Lp400Re≈Lp2-22-3CFD2.1.2.12-415Fig.2-4Controlvolumeforcalculation[25]()airasurfaceairsidedQttdAα−=−(2-9).()kSsurfaceairsidedmxxdAσ−=−(2-10)σxt.Ssurfacexsurfacet.ktairdmmdx=−(2-11).tairzm∗.(1)*tairzm−.(1)tairzm−∗surfacet,,,insurfacedtdxdtt**0.()()0irefrigsidesurfacerefrigsideccpDwaterwatersurfaceptairpAttdxdtcthctcmcα−−−++−+=(2-12).()pDaairsidesurfacetairppcAdtdxttmccα−⎛⎞=−+−⎜⎟⎜⎟⎝⎠(2-13)..()aairsidesSsurfacetairpAdxxxmcαη−⎛⎞=−−⎜⎟⎜⎟⎝⎠(2-14)16pc(1)ppApVcxcxc=−+indt.()irefrigsiderefrigerantsurfacerefrigerantrefrigerantprefrigerantAdtttmcα−=−(2-15)1(1)RSRSairsideAAηη−=−−(2-16).01()lRSSSsurfacexxdLxxη=−−∫(2-17)xt.SsurfacextFinL[25]LL0L1()cosh[(1)]()cosh()FinLsurfaceabLttttab−=−−(2-18)b2.1.2.2aα[26]20.331RePrcaaaacdλαη=⋅⋅⋅⋅(2-19):ReafacemdAη=Prpcηλ=⋅c1c2Reo[26]η[27]η171(1)RRairsideAAηη−=−−(2-20)[27]r0Rη00tanh()/Rararηϕϕ=(2-21)2aFinasαλ=aαRη2.1.2.3[28][28]iα333iKBααα=+(2-22)2-5BαKαBαKαFig.2-5CourseofheattransfercoefficientBαKαoverthetubelength2.1.2.3.1x18Kα1.50.60.010.352.20.010.70.6720.5[((1)1.9(1)())((18(1)()))]GKFlFlxxxxxααρραραρ−−−′′=−+−++−′′′′(2-23)3.5/1500ρρ′′′≤≤FlαGα[28],/FlGNudαλ=(2-24)Nu331/3(3.661.077RePr/)Nudz=+(2-25)Re2300()()()2/3/8Re1000Pr112.7/8Pr1Nuξξ−=+−(2-26)(Re1.64)210(1.82log)ξ−−=(2-27)Re2300Pr/pcηλ=[28]xKα0.40.010.372.20.010.70.6720.5[((1)1.2(1)())((18(1)()))]GKFlFlxxxxxααρραραρ−−−′′=−+−++−′′′′(2-28)3.5/1500ρρ′′′≤≤Kα2.1.2.3.2[28]Bα0()()(.)()0()nBFdwxpqCFFFFqραα+=(2-29)190.270.36FCM=(2-30)*0.760.80.110pn=−⋅R134a(2-31)*p*/criticalppp=0.453.7(*)71.72.816*(3.4)*1*pFppp=++−(2-32)0.4()0(/)dFdd=d00.01m(2-33)0.133()0(/)waaFRR=Ra01m(2-34)0q0αR134a202000/qWm=203500/WmKα=[28]Bα25%[28]OBq2()sFlOBKrTqrhhσαρ=′′′′′−(2-35)Krr0.3e-6mF