Skip to content

turbulentHeatFluxTemperature cannot be used anymore with buoyantBoussinesqSimpleFoam

Description : In OpenFOAM 2.4.0 it was possible to use the boundary condition turbulentHeatFluxTemperature with an incompressible solver like buoyantBoussinesqSimpleFoam. But in OpenFOAM v3.0+ I face the error:

--> FOAM FATAL IO ERROR: Unknown patchField type turbulentHeatFluxTemperature for patch type wall

Steps to reproduce:

Take the hotRoom tutorial: tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/hotRoom

Change the patch floor in the 0/T.org file to:

floor { type turbulentHeatFluxTemperature; gradient uniform 0; heatSource power; q uniform 64; alphaEff alphaEff; value uniform 300; }

Launch the Allrun

SIMPLE: convergence criteria field p_rgh tolerance 0.01 field U tolerance 0.0001 field T tolerance 0.01 field "(k|epsilon|omega)" tolerance 0.001

Reading thermophysical properties

Reading field T

--> FOAM FATAL IO ERROR: Unknown patchField type turbulentHeatFluxTemperature for patch type wall

Valid patchField types are :

104 ( advective alphatJayatillekeWallFunction atmBoundaryLayerInletEpsilon atmBoundaryLayerInletK calculated codedFixedValue codedMixed compressible::alphatJayatillekeWallFunction compressible::alphatWallFunction compressible::thermalBaffle1D compressible::thermalBaffle1D compressible::turbulentHeatFluxTemperature

Additional information :

If we repeat the same procedure in OF2.4.0 : tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/hotRoom

  1. change the patch floor 0/T.org:

floor { type turbulentHeatFluxTemperature; gradient uniform 0; heatSource power; q uniform 64; alphaEff alphaEff; value uniform 300; }

  1. add the line rhoCp0 1173; in transportProperties

  2. launch the Allrun

SIMPLE: convergence criteria field p_rgh tolerance 0.01 field U tolerance 0.0001 field T tolerance 0.01 field "(k|epsilon|omega)" tolerance 0.001

Starting time loop

Time = 1

DILUPBiCG: Solving for Ux, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for T, Initial residual = 1, Final residual = 0.0475703, No Iterations 1 DICPCG: Solving for p_rgh, Initial residual = 1, Final residual = 0.00653368, No Iterations 7 time step continuity errors : sum local = 5.65209e-09, global = -4.97342e-25, cumulative = -4.97342e-25 DILUPBiCG: Solving for epsilon, Initial residual = 0.0458815, Final residual = 0.00188758, No Iterations 1 DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 0.0721881, No Iterations 1 ExecutionTime = 0.05 s ClockTime = 0 s

Time = 2

DILUPBiCG: Solving for Ux, Initial residual = 0.651198, Final residual = 0.00701648, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.558322, Final residual = 0.00593882, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.651198, Final residual = 0.00701648, No Iterations 1 DILUPBiCG: Solving for T, Initial residual = 0.431799, Final residual = 0.0275141, No Iterations 1 DICPCG: Solving for p_rgh, Initial residual = 0.86867, Final residual = 0.00688476, No Iterations 27 time step continuity errors : sum local = 3.05697e-07, global = -6.73656e-24, cumulative = -7.2339e-24 DILUPBiCG: Solving for epsilon, Initial residual = 0.115159, Final residual = 0.00602713, No Iterations 1 DILUPBiCG: Solving for k, Initial residual = 0.624686, Final residual = 0.0490608, No Iterations 1 ExecutionTime = 0.06 s ClockTime = 0 s

OK no error.

why in OF240 the turbulentHeatFluxTemperature is present in incompressible https://github.com/OpenFOAM/OpenFOAM-2.4.x/blob/2b147f41daf9ca07d0fb4c6b0576dc3d10a435f3/src/turbulenceModels/incompressible/turbulenceModel/derivedFvPatchFields/turbulentHeatFluxTemperature/turbulentHeatFluxTemperatureFvPatchScalarField.H

but not in OFv3.0+? https://develop.openfoam.com/search?utf8=%E2%9C%93&search=turbulentHeatFluxTemperature&group_id=&project_id=5&search_code=true&repository_ref=master