View Issue Details

IDProjectCategoryView StatusLast Update
0001751OpenFOAM[All Projects] Bugpublic2015-06-18 11:48
ReporterdhoraAssigned Tohenry 
PrioritynormalSeverityminorReproducibilityalways
Status resolvedResolutionfixed 
Product Versiondev 
Fixed in Version 
Summary0001751: div(phi, symmTensor) with vanLeer and QUICK
DescriptionNot sure whether this is really a bug but the behaviour seems to be a bit strange at first sight. I didn't investigate further. The problem can be reproduced with a slightly modified version of scalarTransportFoam for tensors and symmTensors and the corresponding tutorial.

The symmTensor results start to differ from the scalar and tensor values in the second time step. It seems to converge to identical values at the end but at the beginning there is a significant discrepancy. I observed this behaviour with the vanLeer and QUICK schemes. Linear, upwind or SFCD for example are not affected.

Output from modified solver:

Time = 0.0002

Value Cell 25
(0.00091745 0 0 0 0.00091745 0 0 0 0.00091745)
(0.00091745 0 0 0.00091745 0 0.00091745)
0.00091745
Value Cell 0
(0.000912834 0 0 0 0.000912834 0 0 0 0.000912834)
(0.000912834 0 0 0.000912834 0 0.000912834)
0.000912834
Value Cell 1
(0.00218276 0 0 0 0.00218276 0 0 0 0.00218276)
(0.00218276 0 0 0.00218276 0 0.00218276)
0.00218276
Value Cell 2
(0.000400171 0 0 0 0.000400171 0 0 0 0.000400171)
(0.000400171 0 0 0.000400171 0 0.000400171)
0.000400171
Value Cell 3
(0.000923626 0 0 0 0.000923626 0 0 0 0.000923626)
(0.000923626 0 0 0.000923626 0 0.000923626)
0.000923626
Laplacian(DT, X)
(5.11345 0 0 0 5.11345 0 0 0 5.11345)
(5.11345 0 0 5.11345 0 5.11345)
5.11345
div(phi, X)
(-5.85312 0 0 0 -5.85312 0 0 0 -5.85312) // TENSOR IDENTICAL TO SCALAR
(-4.91988 0 0 -4.91988 0 -4.91988) // SYMMTENSOR DIFFERENT
-5.85312 // TENSOR IDENTICAL TO SCALAR
DILUPBiCG: Solving for Txx, Initial residual = 0.173109, Final residual = 3.59718e-07, No Iterations 10
DILUPBiCG: Solving for Txy, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Tyx, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Tyy, Initial residual = 0.173109, Final residual = 3.59718e-07, No Iterations 10
Steps To ReproduceUse the attached files and solver for the scalarTransportFoam tutorial.
TagsNo tags attached.

Activities

dhora

2015-06-17 23:44

reporter  

tensorTransportFoam.tar.gz (2,659 bytes)

henry

2015-06-18 11:48

manager   ~0004965

Thanks for the detailed bug-report and test code and case.

Resolved by commit 0575911f947d2ab2bb31e315361728d0b712ceb6 in OpenFOAM-dev

Issue History

Date Modified Username Field Change
2015-06-17 23:44 dhora New Issue
2015-06-17 23:44 dhora File Added: tensorTransportFoam.tar.gz
2015-06-18 11:48 henry Note Added: 0004965
2015-06-18 11:48 henry Status new => resolved
2015-06-18 11:48 henry Resolution open => fixed
2015-06-18 11:48 henry Assigned To => henry