《FLuent-换热器的相变模拟计算》由会员分享,可在线阅读,更多相关《FLuent-换热器的相变模拟计算(28页珍藏版)》请在金锄头文库上搜索。
1、 UGM 2001New Initiatives at Fluent Inc.Phase Change in Heat ExchangersBrian Bell, Fluent Inc.UGM 20011Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001MotivationDemonstrate the use of Fluent to model phase change in heat exchangersProcesses of interestuCondensationuEvapo
2、rationuBoiling Illustrate how to model one such process through use of a detailed exampleShell-and-tube condenserProvide motivation for users to begin developing models of their own devices2Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001OutlineProblem DescriptionShell-
3、and-tube condenseruPure vapor condensationuNon-condensable gasesModeling ApproachPorous mediumHeat and mass transfer modelingModel ImplementationUser-Defined Functions and User-Defined MemoryResultsSteam condenser with non-condensable gasesCommercial chiller condenser3Company Confidential Copyright
4、2001 Fluent Inc. All rights reserved. UGM 2001Description of ProblemShell-and-tube condenser 4Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001Goals of CFD ModelingCondenser performance characterized by heat and mass transfer rateCFD allows evaluation of factors affecti
5、ng heat and mass transfer in condenseruTube bundle configurationtTube arrangementtNumber of passestLocation of inlet portstBafflesuPressure dropuVelocity fielduNon-condensablestLocation and configuration of purge systemResults allow identification of potential design improvements5Company Confidentia
6、l Copyright 2001 Fluent Inc. All rights reserved. UGM 2001Film Condensation ProcessDriving potential for condensation is the temperature difference between vapor and cooling waterDriving potential variation caused by Pressure dropRise of cooling water temperatureNon-condensablesTPH2OPairairCondensa
7、te layerTube wallCooling water6Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001CFD Modeling TheoryPorous medium approachTube bundle treated as porous mediumEnables computationally efficient modeling of entire condenserComparison with detailed modeling approachuIn 2-D, O
8、(100)-O(1000) control volumes per tube versus more than one tube per control volumeHeat and mass transfer modelsCondensation rate calculationuCondensation rate determined from local flow field and cooling water temperatureuLiquid film flow rate tracked in bundle from top to bottomuCooling water temp
9、erature tracked from inlet to outlet7Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001Porous Medium ApproachRepresentation of tube bundle as porous mediumPorosity is only required parameterPorosity defined as ratio of fluid volume to total volume PduExample: staggered tu
10、be bundle with equilateral triangular layoutPorosity, b, expressed as:8Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001Transport EquationsGeneric transport equation for porous medium approachconvectiondiffusiondistributed resistanceEqn.continuity1x-mom.uy-mom.vspeciesw
11、w Distributed resistance takes form of source terms that model details of the flow that are not resolved by the gridPorosity in convection and diffusion terms not modeled in FluentDistributed resistance terms most significant in tube bundle region9Company Confidential Copyright 2001 Fluent Inc. All
12、rights reserved. UGM 2001Evaluation of Modeling ApproachAdvantagesComputationally efficientuDoes an alternate, tractable approach exist?Approach demonstrated to give meaningful data by several authorsDisadvantagesLoss of some flow details due to averagingCan be overcome by detailed modeling of smal
13、l regions of condenser10Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001Heat Transfer ProcessFilm condensation on horizontal tubeCooling WaterTube WallCondensate FilmLiquid-vapor InterfaceRefrigerantVaporLatent heat released at liquid-vapor interface transferred to cool
14、ing water11Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001Heat Transfer ModelHeat transfer is modeled by coupling of thermal resistance network with CFD codeTcwTt,iTt,oTiRcwRtubeRcondCooling WaterCFD code provides interface temperature, Ti Cooling water and tube therma
15、l resistances are generally well-knownFilm heat transfer coefficient is required for Rcond12Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001Film Heat Transfer CoefficientCritical component of heat transfer modelObtain from experimentOr obtain from literatureSteam conden
16、sation on smooth tubesFigure courtesy of Kansas State University,Professor Steve Eckles, and Duane L. Randall13Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001Modeling AssumptionsEffect of liquid on flow field is neglectedApproach can also be implemented in Eulerian-Eul
17、erian multiphase frameworkuSatisfactory model for liquid phase representation not currently availableuPublished results of this type of model do not appear to show significant advantageVapor is assumed to be saturatedNo superheatingVapor temperature determined from pressure field calculated by CFD c
18、ode14Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001Implementation of Model with UDFsUDFs are required for:Source terms required by porous medium approachuCondensation rateuPressure drop in porous regionRepresentation of tube bundleuPorosityuCondensate film flow rate a
19、ccountinguCooling water temperature calculation with multiple tube passes15Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001Cooling water temperature calculation for each segmentuEvery iteration, condensation rate is summed over each segmentuInlet cooling water temperatu
20、re = outlet temperature from previous segmentuSegment outlet cooling water temperature calculated by energy balance.uLog-mean temperature for each segment calculated based on vapor temperature and cooling water inlet and outlet temperaturesTube Bundle RepresentationBundle consists of N passes and M
21、segmentsEach segment defined as unique cell zoneExample:2 Pass bundleN = 2, M = 416Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001Tube Bundle Grid StructureStructured, cartesian grid used in tube bundleEach control volume has unique i,j,k indexi=1j=1k=1i=1j=2k=1i=1j=3k
22、=1i=1j=3k=2i=1j=2k=2i=1j=1k=2i =1j=1k=3i =1j=2k=3i=1j=3k=3Grid structure created with UDFsuGrid generator, solver do NOT utilize structureUsed to track condensate film flow rate17Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001Source TermsAlgorithm for source term in co
23、ntinuity equationObtain pressure, velocity and species mass fraction (if necessary) from current solution valuesObtain film Reynolds number and cooling water temperature from User-Defined MemoryCalculate heat flux based on current value of solution variables Translate heat flux into volumetric mass
24、source termUnder-relax source termuSi+1 = Si + a (So Si)uRequired for solution stability. Alpha typically 0.01 0.10uValue of source term from previous iteration, So, stored in User-Defined MemorySource term in momentum equations Calculated using empirical correlations with tube bundle porosity and c
25、urrent velocity18Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001Define_On_Demand FunctionsDefine_On_Demand functions executed once per iterationUpdate condensate film mass flow rateUpdate cooling water temperatureuAssume uniform temperature for each bundle segmentNew v
26、alues stored in User-defined memoryAutomatic Define_On_Demand execution possibleExample:19Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001Solution AlgorithmInitialize Solution: Assign porosity, tube bundle orientationUpdate cooling water temperature and liquid condensat
27、e mass flow rateCalculate source termsSolve flow equationsYesNoSolution Converged?Stop20Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001ExamplesSteam condensation with non-condensable gasesMcAllister Condenserfrom: Bush et al., 1990, Proc. Int. Symp. On Condensers and C
28、ondensation21Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001McAllister Condenser GeometryBoundary conditions and model inputsShell Dimensions 1.02 m X 1.22 m X 0.78 mCooling water flow directionInlet temperature: 17.8 CInlet velocity: 1.19 m/sTube BundleSingle pass, 4
29、segmentsOuter Diameter: .0254 mInner Diameter: .0242 mPitch: .0349 mPorosity: 0.52PurgeMass flow rate: .011 kg/sInletPressure: 27670 PaAir mass fraction: 0.0012222Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001Condenser Grid15,000 Control VolumesSimple geometry allows
30、structured grid throughout domainGrid profile in x-z plane23Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001ResultsCondensation RateInlet mass flow rateCFD: 2.124 kg/sExp.: 2.032 kg/sError: 4.5%Cooling water temperature contoursVolumetric condensation rate contours24Com
31、pany Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001McAllister Condenser Flow FieldVelocity MagnitudeMax: 34.4 m/sMin: 0.02 m/sPressureMax: 27,663 PaMin: 27,530 PaAir Mass FractionMax: .534Min: .00122Condensation Rate *Max: 6.1 kg/smMin: 0.0 kg/sm* Minimum condensation rate in
32、 tube bundle is 0.18 kg/sm25Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001Effect of Air on Condensation RateVolumetric condensation rate contours without airVolumetric condensation rate contours with airComparison of modeling results with and without non-condensable g
33、ases26Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001Effect of Inundationon Condensation RateTracking condensate film flow rate from upper tubes to lower tubes allows use of inundation correction factorContours of condensate film mass flow rateContours of volumetric co
34、ndensation rate27Company Confidential Copyright 2001 Fluent Inc. All rights reserved. UGM 2001ConclusionsModeling phase change processes in heat exchangers is possible in Fluent 5 through the use of UDFs This approach is well-suited for falling film condensation or evaporation processesHeat and mas
35、s transfer models must be provided by the userCFD modeling of McAllister condenserDescription of model development process intended to serve as a reference for users who wish to develop similar modelsAccurate results demonstrate the potential of this approach28Company Confidential Copyright 2001 Fluent Inc. All rights reserved.