View Issue Details

IDProjectCategoryView StatusLast Update
0002500OpenFOAMBugpublic2017-03-22 18:41
Reporterolle Assigned Tohenry  
PrioritynormalSeveritymajorReproducibilityalways
Status closedResolutionunable to reproduce 
Platformx86_64OSCentOSOS Version6.8
Summary0002500: Heat equation source term behaviour differ in chtMultiRegionFoam and chtMultiRegionSimpleFoam
DescriptionI have made a simple test case consising of an isolated solid "box" (there is actually two solids in the model in order to the the multi region approach to work). No fluid regions.

The solid is 1x1x1 m and has an internal heat source implemented using fvOptions and the scalarSemiImplicitSource. The heat source is 1 W/m3.
Boundaries of the box are either zeroGradient (which don't make sense in a steady state simulation with a constant source) or externalWallHeatFluxTemperature with constant heat transfer coefficient.

Running first the transient solver with zeroGradient conditions I get the expected temprature increase rate of the box, corresponding to the power input of the source.

Running then the transinet solver the heat transfer coefficient BC I also get the expected steady state temperature.

But running the same case with the steady state solver gives no temperature rise at all.
Steps To ReproduceRun the provided test case with the two solvers using the Allrun script

- alter system/box/changeDictionaryDict to change BC
- alter system/controlDict to change solver
- alter the fvSolution and fvSchemes files to fit the solver.
TagschtMultiregionFoam, fvOptions, Heat transfer

Activities

olle

2017-03-15 15:02

reporter  

sourceTestBox.tar.gz (3,908 bytes)

henry

2017-03-20 18:28

manager   ~0007953

I ran the Allrun which failed due to problems with links in the system directory:

 fvSolution -> /home/setuotma/OpenFOAM/setuotma-4.1/run/sealHeat/sourceTestBox/system/box/fvSolution

etc. I fixed these and the case ran with chtMultiRegionSimpleFoam but not chtMultiRegionFoam. Could you provide the script/setup you use for chtMultiRegionFoam?

henry

2017-03-20 18:37

manager   ~0007954

injectionRateSuSp
        {
            h (1 0);
        }

is a REALLY small injection rate of enthalpy, I would expect a REALLY small temperature rise which is what you get. If you set the "Su" value to something more reasonably you get the corresponding temperature rise.

Issue History

Date Modified Username Field Change
2017-03-15 15:02 olle New Issue
2017-03-15 15:02 olle File Added: sourceTestBox.tar.gz
2017-03-15 15:02 olle Tag Attached: chtMultiregionFoam
2017-03-15 15:02 olle Tag Attached: fvOptions
2017-03-15 15:02 olle Tag Attached: Heat transfer
2017-03-20 18:28 henry Note Added: 0007953
2017-03-20 18:37 henry Note Added: 0007954
2017-03-22 18:41 henry Assigned To => henry
2017-03-22 18:41 henry Status new => closed
2017-03-22 18:41 henry Resolution open => unable to reproduce