|
[Sponsors] |
calculation of kappa in buoyantBoussinesqPimpleFoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
November 6, 2012, 08:52 |
calculation of kappa in buoyantBoussinesqPimpleFoam
|
#1 |
New Member
Join Date: Nov 2010
Posts: 9
Rep Power: 15 |
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; kappat.correctBoundaryConditions(); 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; kappat.correctBoundaryConditions(); 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. |
|
November 6, 2012, 09:53 |
|
#2 |
New Member
Join Date: Nov 2010
Posts: 9
Rep Power: 15 |
ok, i might have been wrong all the time. So, kinematic thermal conductivity is not the same as normal thermal conductivity coefficient...
|
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Transfer output data as input for second calculation | mannobot | Main CFD Forum | 3 | December 7, 2018 09:37 |
Questions on Boundary Layer Spring Method for 3D Calculation | didiean | FLUENT | 15 | December 18, 2014 21:09 |
MRF and Heat transfer calculation | Susan YU | FLUENT | 0 | June 2, 2010 08:46 |
Warning 097- | AB | Siemens | 6 | November 15, 2004 04:41 |
Heat Flux Calculation under REPORTS | Alberto Schroth | FLUENT | 0 | May 16, 2000 08:19 |