
[Sponsors] 
July 28, 2011, 01:39 
Problem with BICGSTAB

#1 
New Member
Amir
Join Date: Oct 2010
Posts: 16
Rep Power: 7 
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: 7 
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: 7 
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 nonphysical pressure gradients? 

October 14, 2011, 13:16 

#4 
New Member
Amir
Join Date: Oct 2010
Posts: 16
Rep Power: 7 
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  


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  SeHee  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 