CFD Online Discussion Forums - forced convection on a cylinder

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)

- - forced convection on a cylinder (https://www.cfd-online.com/Forums/openfoam/92593-forced-convection-cylinder.html)

forced convection on a cylinder

Hello,

I want to simulate a crosswise to flow to an heated cylinder and compare my results with a analytic solution. The cylinder is heated with a fixed value of 1000 W/m^2 and the air has at the inlet a temperature of 293 Kelvin. When I run the case, it will converge, but the results are nonsense. The velocity profile looks good, but the temperature on the cylinder are nearly the same as the surrounding fluid. In my analytic solution the temperature of the cylinder is about 457 Kelvin.

I've uploadete the case here: http://dl.dropbox.com/u/1445847/forced_convection.zip

Thanks for any help

greetings

can anyone give me a hint?

Hey Maximilian,

Which solver are you using? Where do you have find the specification for:

Code:

cylinder                                                                                                                                                                    

        {                                                                                                                                                                           

        type    compressible::turbulentHeatFluxTemperature;                                                                                                                         

        heatSource      flux;                                                                                                                                                       

        q       uniform 1000;                                                                                                                                                       

        K       basicThermo;                                                                                                                                                        

        KName   K;                                                                                                                                                                  

        value   uniform 460;                                                                                                                                                        

        }

I've been searching for something like that but not havent found anything neither in the tutorials nor in the UserGuide

Hi,
I got the specification for this Bc from the source Code:

Code:

$FOAM_SRC//turbulenceModels/compressible/turbulenceModel/derivedFvPatchFields/turbulentHeatFluxTemperature/turbulentHeatFluxTemperatureFvPatchScalarField.H

I'm using the solver rhoSimpleFoam.

regards Maximilian