Transcription of Heat transfer in OpenFOAM
1 heat transfer in OpenFOAMDr. Johann Turnow,silentdynamics GmbH2018-02-22c Copyright silentdynamics GmbH2018-02-221/81 ContentsClassification | Solver | Conduction | Convection | CHT | Radiation |ClassificationSolverConductionConvectio nCHTR adiationc Copyright silentdynamics GmbH2018-02-222/81 Main heat transfer typesClassification | Solver | Conduction | Convection | CHT | Radiation |1. heat conduction2. heat convection3. heat radiationIsteady, transientIcompressible, Copyright silentdynamics GmbH2018-02-223/81 Solver overviewClassification | Solver | Conduction | Convection | CHT | Radiation |Overview of OpenFOAM solvers for heat transfer analysisIlaplacianFoam:Transient, incompressible, thermal diffusion according to Fourier slawIscalarTransportFoam:Steady-state, incompressible, laminar, passive scalar for a given velocity fieldIbuoyantBoussinesqSimpleFoam.
2 Steady-state, thermal, natural convection, incompressible,Boussinesq s approximationIbuoyantBoussinesqPimpleFoa m:Transient, thermal, natural convection, incompressible, Boussinesq sapproximationc Copyright silentdynamics GmbH2018-02-224/81 Solver overviewClassification | Solver | Conduction | Convection | CHT | Radiation |Overview of OpenFOAM solvers for heat transfer analysisIbuoyantSimpleFoam:Steady-state, natural convection, compressible (sub-sonic), includingradiationIbuoyantPimpleFoam:tra nsient, natural convection, compressible(sub-sonic), includingradiationIrhoSimpleFoam:Steady- state, thermal, compressible(sub-sonic)IrhoSimplecFoam:S teady-state, thermal, compressible(sub-sonic) -Pressure underrelaxiation =1 IrhoPimpleFoam:Transient, thermal, compressible(sub-sonic)c Copyright silentdynamics GmbH2018-02-225/81 Solver overviewClassification | Solver | Conduction | Convection | CHT | Radiation |Overview of OpenFOAM solvers for heat transfer analysisIchtMultiRegionFoam.
3 Transient, compressible, conjugate heat transfer between solid andfluidIchtMultiRegionSimpleFoam:Steady -state, compressible, conjugate heat transfer between solidand fluidIthermoFoam:Transient, evolves the thermophysical properties for a frozen velocityfieldc Copyright silentdynamics GmbH2018-02-226/81 heat conductionClassification | Solver | Conduction | Convection | CHT | Radiation |Basic solver: laplacianFoamISimple heat conduction equation according to Fourier s law T t= cp 2T x2(1)ITake a look at the solverIcd $FOAM_SOLVERS or solIcd basic/laplacianIgedit (fvm::ddt(T) - fvm::laplacian(DT, T));c Copyright silentdynamics GmbH2018-02-227/81 heat conductionClassification | Solver | Conduction | Convection | CHT | Radiation |Basic solver: laplacianFoamIDefine the heat diffusivityDT:Igedit constant/transportProperties//DT = heat diffusivityDT DT [ 0 2 -1 0 0 0 0 ] ; //air//DT DT [ 0 2 -1 0 0 0 0 ].
4 //water//DT DT [ 0 2 -1 0 0 0 0 ] ; //aluc Copyright silentdynamics GmbH2018-02-228/81 heat conductionClassification | Solver | Conduction | Convection | CHT | Radiation |Example coffee cupIUsinglaplacianFoamto simulation usual problemsILet s try to analyze the temperature distribution in our coffee cupIQuestion: Can you touch the cup without any pain?c Copyright silentdynamics GmbH2018-02-229/81 heat conductionClassification | Solver | Conduction | Convection | CHT | Radiation |Example coffee cupISetting the boundary conditionsIgedit 0/TinternalField uniform 273;boundaryField{sideWalls{type zeroGradient; //adiabatic}coffee{type fixedValue; // fixed Temperature uniform 373.}}
5 }}c Copyright silentdynamics GmbH2018-02-2210/81 heat conductionClassification | Solver | Conduction | Convection | CHT | Radiation |Example coffee cupISetting the boundary conditionsIgedit 0/TinternalField uniform 273;boundaryField{sideWalls{type zeroGradient; //adiabatic}coffee{type fixedGradient; //fixed heat flux 10000;value uniform 373;}}c Copyright silentdynamics GmbH2018-02-2211/81 heat conductionClassification | Solver | Conduction | Convection | CHT | Radiation |Example coffee cupIDefine the heat diffusivityDTforalu:Igedit constant/transportProperties//DT = heat diffusivity//DT DT [ 0 2 -1 0 0 0 0 ] ; //air//DT DT [ 0 2 -1 0 0 0 0 ] ; //waterDT DT [ 0 2 -1 0 0 0 0 ].
6 //aluIdecomposeParIfoamJob -parallel laplacianFoamItail -f logc Copyright silentdynamics GmbH2018-02-2212/81 heat conductionClassification | Solver | Conduction | Convection | CHT | Radiation |Example coffee cupITake a look at the temperature after for ouralucupIThealugives pretty hot fingers after ,c Copyright silentdynamics GmbH2018-02-2213/81 heat conductionClassification | Solver | Conduction | Convection | CHT | Radiation |Example coffee cupIComparison to a usual porcelain cupITheporcelaincup gives us cool fingers fingers after (.)
7 C Copyright silentdynamics GmbH2018-02-2214/81 heat conductionClassification | Solver | Conduction | Convection | CHT | Radiation |OutcomeILaplacian solver gives a fairly good overview for simple heatconduction problemsIAlways the first choice for simple heat conduction solutionsIFirst step: Think about which results you expectIImportant to avoid nonphysical solutions .. :-)IAlways take a look at the residualsIAlways remember that the mesh resolution influences the results incase of heat transfer dramatically!IA Priori: Which boundary conditions should be applied?
8 IBe careful with the constant heat flux boundary conditionc Copyright silentdynamics GmbH2018-02-2215/81 heat convectionClassification | Solver | Conduction | Convection | CHT | Radiation |Which solvers can we use?IscalarTransportFoamfor laminar, unsteady/steady flowsIbuoyantBoussinesqSimpleFoam:Steady -state, thermal, natural convection, incompressible,Boussinesq s approximationIbuoyantBoussinesqPimpleFoa m:Transient, thermal, natural convection, incompressible, Boussinesq sapproximation Set the gravitation to Zero for simple passive scalar flowsc Copyright silentdynamics GmbH2018-02-2216/81 heat convectionClassification | Solver | Conduction | Convection | CHT | Radiation |Wich equation is solved?
9 VolScalarField alphaEff("alphaEff", turbulence->nu()/Pr+ alphat);fvScalarMatrix TEqn(fvm::div(phi, T)- fvm::laplacian(alphaEff, T)==radiation->ST(rhoCpRef, T)+ fvOptions(T));c Copyright silentdynamics GmbH2018-02-2217/81 heat convectionClassification | Solver | Conduction | Convection | CHT | Radiation |Let s take a look at our cup!IQuestion: How much is the coffee cooled down when you hold thecup in the cold wind of 0 C and a wind speed of Copyright silentdynamics GmbH2018-02-2218/81 heat convectionClassification | Solver | Conduction | Convection | CHT | Radiation |Let s take a look at our cup!
10 IOnly the fluid is treated firstc Copyright silentdynamics GmbH2018-02-2219/81 heat convectionClassification | Solver | Conduction | Convection | CHT | Radiation |Let s take a look at our cup!IDo we need turbulence?Re=U L =1m/s 10 06m/s2=16666(2)IYes weed need model kOmegaSST (wallbounded)Igedit constant/RASP ropertiessimulationType RAS;RAS{RASM odel kOmegaSST;turbulence on;printCoeffs on;}c Copyright silentdynamics GmbH2018-02-2220/81 heat convectionClassification | Solver | Conduction | Convection | CHT | Radiation |Let s take a look at our cup!