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

Problem with BICGSTAB

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 28, 2011, 01:39
Smile Problem with BICGSTAB
  #1
New Member
 
Amir
Join Date: Oct 2010
Posts: 16
Rep Power: 6
Amir_Ghasemi is on a distinguished road
Hello all,

I have written a BICGSTAB solver for Poisson equation in SPH (Smoothed particle hydrodynamics) method. As all the boundaries are Neumann kind, the resulted matrix is singular, to remove singularity, pressure of one point in the domain is fixed. The solution converge satisfactory but results show inhomogeneity(sharp peak) near the fixed point which is nonphysical.

Any idea that why the issue arise is appreciated.

I have one idea by myself but I am not sure how correct it is. There is some compatibility condition which should be satisfied between Neumann BCs and source term of the Poisson equation. Such constraint dictate that sum of the source terms should equal zero (in case of Neumann pressure BCs and impermeable boundaries). Now when we fix one point in the domain, what will happen to this compatibility condition?, should it be satisfied any more and if yes, how we can do this?

Thank you so much.
Amir_Ghasemi is offline   Reply With Quote

Old   July 31, 2011, 10:42
Default
  #2
New Member
 
Amir
Join Date: Oct 2010
Posts: 16
Rep Power: 6
Amir_Ghasemi is on a distinguished road
Hello all,

I have just found what the problem was. May matrix was singular with constant vector in null space and its determinate was zero and it had emerged from discretizing PPE with Neumann boundary. One should be aware that these are not enough to decide to remove one equation (or fixing the pressure of one point). you can do this, if the system of equation has infinite solutions. Ax = b has infinite solution if det(A) = 0 and b is in the range of A, if these condition exist you can remove one equation and be sure that the result of new system also satisfy the original one and in may system this was not the case.

Ghasemi

Last edited by Amir_Ghasemi; August 4, 2011 at 01:28.
Amir_Ghasemi is offline   Reply With Quote

Old   October 13, 2011, 09:48
Default
  #3
New Member
 
Tom
Join Date: Mar 2011
Posts: 12
Rep Power: 6
TheBoyce is on a distinguished road
Hello Ghasemi,

I have the same problem as you initially described in the first post. How did you solve this and obtain a solution without the non-physical pressure gradients?
TheBoyce is offline   Reply With Quote

Old   October 14, 2011, 13:16
Default
  #4
New Member
 
Amir
Join Date: Oct 2010
Posts: 16
Rep Power: 6
Amir_Ghasemi is on a distinguished road
Hello TheBoyce,

As I have said, I think that before fixing the pressure of one node you should become sure that your matrices possess infinite solution, if this is not the case you obtain non physical distribution around the fixed point. I was dealing with a Poisson eq. in SPH method and I changed my pressure boundary condition to get a symmetric matrix with constant vector in it's null space. for this kind of matrix, sum of your source terms should be zero. I have done this with subtracting each source term from the average.

Again I emphasize that zero determination does not always lead to infinite solution, no solution is also possible. in first case you can fix a point, in second case, fixing a point is not allowed.
Amir_Ghasemi is offline   Reply With Quote

Reply

Tags
bicgstab, neuman boundary, sph

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
UDF compiling problem Wouter Fluent UDF and Scheme Programming 6 June 6, 2012 04:43
Gambit - meshing over airfoil wrapping (?) problem JFDC FLUENT 1 July 11, 2011 05:59
natural convection problem for a CHT problem Se-Hee CFX 2 June 10, 2007 06:29
Adiabatic and Rotating wall (Convection problem) ParodDav CFX 5 April 29, 2007 19:13
Is this problem well posed? Thomas P. Abraham Main CFD Forum 5 September 8, 1999 14:52


All times are GMT -4. The time now is 23:45.