View Issue Details
| ID | Project | Category | View Status | Date Submitted | Last Update |
|---|---|---|---|---|---|
| 0001751 | OpenFOAM | Bug | public | 2015-06-17 23:44 | 2015-06-18 11:48 |
| Reporter | dhora | Assigned To | henry | ||
| Priority | normal | Severity | minor | Reproducibility | always |
| Status | resolved | Resolution | fixed | ||
| Product Version | dev | ||||
| Summary | 0001751: div(phi, symmTensor) with vanLeer and QUICK | ||||
| Description | Not 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 Reproduce | Use the attached files and solver for the scalarTransportFoam tutorial. | ||||
| Tags | No tags attached. | ||||
|
|
|
|
|
Thanks for the detailed bug-report and test code and case. Resolved by commit 0575911f947d2ab2bb31e315361728d0b712ceb6 in OpenFOAM-dev |
| 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 |
