# Problem with BICGSTAB

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

 July 28, 2011, 01:39 Problem with BICGSTAB #1 New Member   Amir Join Date: Oct 2010 Posts: 16 Rep Power: 6 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.

 July 31, 2011, 10:42 #2 New Member   Amir Join Date: Oct 2010 Posts: 16 Rep Power: 6 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.

 October 13, 2011, 09:48 #3 New Member   Tom Join Date: Mar 2011 Posts: 12 Rep Power: 6 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?

 October 14, 2011, 13:16 #4 New Member   Amir Join Date: Oct 2010 Posts: 16 Rep Power: 6 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.

 Tags bicgstab, neuman boundary, sph

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Wouter Fluent UDF and Scheme Programming 6 June 6, 2012 04:43 JFDC FLUENT 1 July 11, 2011 05:59 Se-Hee CFX 2 June 10, 2007 06:29 ParodDav CFX 5 April 29, 2007 19:13 Thomas P. Abraham Main CFD Forum 5 September 8, 1999 14:52

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