CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Bugs

calculation of kappa in buoyantBoussinesqPimpleFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   November 6, 2012, 09:52
Default calculation of kappa in buoyantBoussinesqPimpleFoam
New Member
Join Date: Nov 2010
Posts: 9
Rep Power: 9
yosuu is on a distinguished road
Hello Foamers!

I have discovered a serious bug (in my opinion) in the buoyantBoussinesq solvers (version 2.1.1). I checked how the kappa is implemented, and i have my doubts that it is correctly done.

In the original code is:

kappat = turbulence->nut()/Prt;

volScalarField kappaEff("kappaEff", turbulence->nu()/Pr + kappat);

but i think that there should be a Cp multiplying the kappa terms, so the final code

kappat = Cp*turbulence->nut()/Prt;

volScalarField kappaEff("kappaEff", Cp*turbulence->nu()/Pr + kappat);

because from definition of Prt=nut/alphat=nut*Cp*rho/kappa. Since the equations are incompressible, that ends up as Prt=nut*Cp/kappa.

Note: to make it work, Cp must be added in readTransportProperties.H from buoyantBoussinesqSimpleFoam.

What do you think about it? When I run laminar cases, it takes forever to see the temperature field to change from the initial condition when using the original code.
yosuu is offline   Reply With Quote

Old   November 6, 2012, 10:53
New Member
Join Date: Nov 2010
Posts: 9
Rep Power: 9
yosuu is on a distinguished road
ok, i might have been wrong all the time. So, kinematic thermal conductivity is not the same as normal thermal conductivity coefficient...
yosuu is offline   Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Questions on Boundary Layer Spring Method for 3D Calculation didiean FLUENT 15 December 18, 2014 22:09
MRF and Heat transfer calculation Susan YU FLUENT 0 June 2, 2010 08:46
Transfer output data as input for second calculation mannobot Main CFD Forum 1 June 2, 2010 03:55
Warning 097- AB Siemens 6 November 15, 2004 05:41
Heat Flux Calculation under REPORTS Alberto Schroth FLUENT 0 May 16, 2000 08:19

All times are GMT -4. The time now is 14:00.