CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Bugs

Possible bug in fvMatrix residual member

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   April 14, 2011, 03:46
Default Possible bug in fvMatrix residual member
  #1
Member
 
Luca M.
Join Date: Mar 2009
Location: Luzern, Switzerland
Posts: 59
Rep Power: 8
luca is on a distinguished road
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
luca is offline   Reply With Quote

Old   June 6, 2011, 04:35
Default
  #2
Member
 
Luca M.
Join Date: Mar 2009
Location: Luzern, Switzerland
Posts: 59
Rep Power: 8
luca is on a distinguished road
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
luca is offline   Reply With Quote

Old   June 6, 2011, 05:12
Default
  #3
Senior Member
 
Laurence R. McGlashan
Join Date: Mar 2009
Posts: 370
Rep Power: 14
l_r_mcglashan will become famous soon enough
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
l_r_mcglashan is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
transsonic nozzle with rhoSimpleFoam Unseen OpenFOAM Running, Solving & CFD 7 April 16, 2014 03:38
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33
Extrusion with OpenFoam problem No. Iterations 0 Lord Kelvin OpenFOAM 6 April 12, 2011 11:24
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


All times are GMT -4. The time now is 11:35.