View Issue Details

Jump to NotesJump to History
IDProjectCategoryView StatusDate SubmittedLast Update
0001682OpenFOAMBugpublic2015-05-01 04:482017-06-30 10:39
Reporteruser1131Assigned Tochris  
PrioritynormalSeveritymajorReproducibilityalways
Status closedResolutionfixed 
PlatformLinuxOSMintOS Version17
Summary0001682: turbulentTemperatureCoupledBaffleMixed resistance
DescriptionI found that using chtMultiRegionFoam with a constant volumetric heat source and the contact resistance results in erroneous answers as the mesh is refined.

I have added more details in the cfd-online forum:

http://www.cfd-online.com/Forums/openfoam-solving/152444-potential-bug-turbulenttemperaturecoupledbafflemixed.html
TagsNo tags attached.

Activities

henry

2015-05-01 10:11

manager   ~0004710

Have you tested OpenFOAM-2.3.x?

See bug-report http://www.openfoam.org/mantisbt/view.php?id=1447

user1131

2015-05-05 17:30

  ~0004711

I just tested OpenFOAM-2.3.x. I am received with the same problem. The solver can be found here:

https://jessejohns@bitbucket.org/jessejohns/scrapemultiregionfoam

henry

2015-05-05 18:10

manager   ~0004715

Can you reproduce the problem with one of our standard solvers? If so could you provide a simple test-case which demonstrates the issue?

user1131

2015-05-05 18:24

  ~0004716

Last edited: 2015-05-05 18:26

 When defining a heat source on the boundaries, I do not get the same problem. Here, I simply used chtMultiRegionFoam with a volumetric heat source. I think this behavior can be mimicked with fvOptions and chtMultiRegionFoam; however, I could not get this working properly hence my added source term in the solver provided.

Some tests cases are provided.

user1131

2015-05-05 18:24

 

cubeContact_resistence_1900.tar.gz (130,351 bytes)

user1131

2015-05-05 18:25

 

cubeContact_resistence_9700.tar.gz (512,228 bytes)

user1131

2015-05-05 18:26

 

cubeContact_resistence_32000.tar.gz (1,534,677 bytes)

henry

2015-05-05 18:55

manager   ~0004718

I am a bit confused, are you able to reproduce the problem with one of our standard solvers or does the issue only exist with your modified solver?

user1131

2015-05-05 19:47

  ~0004720

Last edited: 2015-05-05 19:52

 None of the standard solvers allow me to apply a volumetric heat source term. I tried using fvOptions, but I failed to get it working.

I apologize, I'll be more thorough:

If I apply a heat flux to the boundary at x=0, using:

compressible::turbulentHeatFluxTemperature;

with the chtMultiRegionSolver and the contact resistance at the interface then I do not get the error of increasing temperature as the mesh is refined. This is, inherently, not the problem I am trying to solve. It seems to me that this kind of definition enforces the heat flux at the interface due to conservation.

In the volumetric heat source case, the heat flux out of the system may not be calculated correctly. Without the contact resistance definition, the solver, scrapeMutliRegionFoam, matches analytical results very well. So, I do think that the heat flux, prior to the contact resistance, is computed correctly. I am working on learning the controlDict functions currently to verify this claim.

In summary:

chtMultiRegionFoam, WITHOUT volumetric heat source, but with a surface heat flux and WITH contact resistance works.

chtMutliRegionFOAM, WITH volumetric heat source (using scrapeMultiRegionFoam) WITHOUT contact resistance works.

chtMutliRegionFOAM, WITH volumetric heat source (using scrapeMultiRegionFoam) WITH contact resistance produces significant errors as the mesh is refined.

user1131

2015-05-09 21:57

  ~0004736

I can't really track down the cause of this problem. The surface flux grows with mesh count on the surface for cases with contact resistance; cases without performed as expected and match analytical results.

chris

2017-06-30 10:39

manager   ~0008292

No feedback from the reporter.

Issue History

Date Modified Username Field Change
2015-05-01 04:48 user1131 New Issue
2015-05-01 10:11 henry Note Added: 0004710
2015-05-05 17:30 user1131 Note Added: 0004711
2015-05-05 18:10 henry Note Added: 0004715
2015-05-05 18:24 user1131 Note Added: 0004716
2015-05-05 18:25 user1131 File Added: cubeContact_resistence_1900.tar.gz
2015-05-05 18:25 user1131 File Added: cubeContact_resistence_9700.tar.gz
2015-05-05 18:26 user1131 Note Edited: 0004716
2015-05-05 18:26 user1131 File Added: cubeContact_resistence_32000.tar.gz
2015-05-05 18:55 henry Note Added: 0004718
2015-05-05 19:47 user1131 Note Added: 0004720
2015-05-05 19:48 user1131 Note Edited: 0004720
2015-05-05 19:52 user1131 Note Edited: 0004720
2015-05-09 21:57 user1131 Note Added: 0004736
2017-06-30 10:39 chris Assigned To => chris
2017-06-30 10:39 chris Status new => closed
2017-06-30 10:39 chris Resolution open => fixed
2017-06-30 10:39 chris Note Added: 0008292