|
[Sponsors] |
April 14, 2011, 03:46 |
Possible bug in fvMatrix residual member
|
#1 |
Member
Luca M.
Join Date: Mar 2009
Location: Luzern, Switzerland
Posts: 59
Rep Power: 17 |
Hello,
supposing that the fvMatrix.residual() member give us the residual of the matrix in the form R = Ax - B, I've noticed a strange behavior. In some way .residual() member applied to fvMatrix<vector> modifies the extra diagonal elements of the matrix: the operator .H() returns different results for fvMatrix<vector> before and after calling the fvMatrix.residual() member. Attached you find a simple code to reproduce the "phenomena" , you can run it under tutorial directory of the scalarTransportFoam solver Tested on 1.7.x version yesterday update scalarTransportFoam.C Thanks in advance |
|
June 6, 2011, 04:35 |
|
#2 |
Member
Luca M.
Join Date: Mar 2009
Location: Luzern, Switzerland
Posts: 59
Rep Power: 17 |
Hello
I was wondering why after two months nobody of the developers reply to this post.. Anyway for the guys that need residual() functions works properly I figure out the error: Change in file fvMatrixSolve.C the line 192 with: tmp<Field<Type> > tres(new Field<Type>(source_)); in this way there is no more modification of the source term (that was the problem). Regards Luca P.S. you have to patch only the 1.7 official release |
|
June 6, 2011, 05:12 |
|
#3 |
Senior Member
Laurence R. McGlashan
Join Date: Mar 2009
Posts: 370
Rep Power: 23 |
If you think there's a bug then post it here (as mentioned in the sticky thread in this forum):
http://www.openfoam.com/mantisbt
__________________
Laurence R. McGlashan :: Website |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
transsonic nozzle with rhoSimpleFoam | Unseen | OpenFOAM Running, Solving & CFD | 8 | July 1, 2022 06:54 |
Extrusion with OpenFoam problem No. Iterations 0 | Lord Kelvin | OpenFOAM Running, Solving & CFD | 8 | March 28, 2016 11:08 |
How to write k and epsilon before the abnormal end | xiuying | OpenFOAM Running, Solving & CFD | 8 | August 27, 2013 15:33 |
Differences between serial and parallel runs | carsten | OpenFOAM Bugs | 11 | September 12, 2008 11:16 |
Unknown error | sivakumar | OpenFOAM Pre-Processing | 9 | September 9, 2008 12:53 |