Skip to content
Snippets Groups Projects

Wrong flux with externalWallHeatFluxTemperature

    • Open Issue created by Rémi Bonin @rbonin

    Summary

    When using the boundary condition externalWallHeatFluxTemperature in coefficient mode and with thicknessLayers, kappaLayers and emissivity different from 0, the computed flux is wrong.

    Steps to reproduce

    1. Mesh a solid square
    2. Left wall and right wall both using externalWallHeatFluxTemperature in coefficient mode with the same Ta, same h and emissivity at 0.9 for both.
    3. On the right wall, add one thicknessLayers at 0.1 and kappaLayers at 1
    4. Top and bottom walls using zeroGradient
    5. Run the case with chtMultiRegionSimpleFoam (laplacianFoam does not accept externalWallHeatFluxTemperature, so you have to manipulate the mesh files a little bit, or you can set up a fluid simulation instead with buoyantSimpleFoam, the same behavior appears).
    6. Check the heat flux at both boundaries

    Example case

    case_report.zip

    What is the current bug behaviour?

    The heat flux converges to around 2000 W

    What is the expected correct behavior?

    The heat flux should converge towards 0 W

    Relevant logs and/or images

    Environment information

    • OpenFOAM version : v2212
    • Operating system : Windows 10
    • Hardware info : 28 cores / 128 GB
    • Compiler : -

    Possible fixes

    Edited

      Child items

      0
      Show labels

      No child items are currently assigned. Use child items to break down this issue into smaller parts.

        Activity

        Please register or sign in to reply