Role of divergence
Hi,
I have finished my CFD calculation and to check it's accuracy, I analized the divergence in each cell. It is in order of 1.0E4 ~ 1.0E5 for the "internal" cells (where velocity and pressure fields are determined entirely by N.S. equation). However, it is in order of 1.0E2 at the worstgrid(near the inlet). Please note that one of the neignbouring cells at this location is the Ghostcell (at which the velocity field is estimated after satisfying the massconservation as mentioned in MAC by Harlow and Welch,...and then the pressure estimation is done using normal component of the N.S. equation). I don't know how much does this worst grid spoil the solution and how to improve the situation. Also, can a same criterion on divergence (e.g. Divergence < 1.0E4 (say)) be used irrespective of grid size to check accouracy of the solution....i.e. if we get same divergence = 9.4E5 for two cells having different grid sizesn a1 x b1 x c1 and a2 x b2 x c2 respectively, can we say that accuracy is almost same at both the cells. Here, "Divergence*Volume = Mass accumulation" in the cell. Thanks in anticipation, chandra. 
Re: Role of divergence
The error will be slightly larger at boundaries because you are enforcing mass flow internally and then updating boundary velocities.
Mass flow conservation and grid resolution effects are two different issues. Whatever the grid size, mass flow should be conserved. Since solution to poisson equation is expensive (every time step), convergence is a tradeoff between accuracy and computational time and you may have to do some tests before you can decide the right numbers for your application. Another check for mass flow rate would be to actually integrate velocity over normal surface area at various cross sections. 
Re: Role of divergence
The convergence criteria may be stricter for very fine grids (as small as 1e07).

Re: Role of divergence
Thanks a lot for the suggestions...however, the variation of convergence criterion based on size of the grid seems to be a little confusing....if we see physical meaning of divergence, it is nothing but 'rate of mass accumulation per unit volume'...and so the selected divergencecriterion should be independent of the grid size...please suggest me if there are other things also which I am not able to point out as far as choosing of this divergencecriterion is concerned..
Also, as I already have mentioned in my previous posting, the inletgrids having neighbouringcells a Ghostcell show very bad performance, could any of you please suggest me about the extent to which this bad cell can spoil the solution. I think answer of it is a matter of experience...and I am very new in it:) thanks in advance, chandra 
Re: Role of divergence
chandra,
Satisfying divergence does depend a lot on the order of accuracy of numerical scheme. For example, using a second order finite difference scheme you might be able to get divergence as low as 1E6. Using spectral methods divergence can be as low as 1E16. Also, in your code, you control divergence through solution to Poisson equation (is that right?). Hence you should use trial/error changing convergence criteria for pressure solution to see how if affects divergence. Regarding higher divergence values at boundaries, i think that should not be a problem, since you do not actually solve NavierStokes equations at the boundary. So as long as your solution looks ok, I wouldn't worry too much about it. 
Re: Role of divergence
heyy, thanks a lot but the divergence near the boundary is really worrying me...
Divergence in the cell adjacent to the inlet having ghostcell a neighbour is having divergence in order of 100 times more than internal grids. I dont know how to improve the situation. I think refining the grid in this problematic region may do it but I am not sure. Please suggest me if there is any other way... Thanks in anticipation... Chandra 
Re: Role of divergence
I'm Sorry I do not know how to alleviate problem related to divergence near boundary. There might be papers that discuss this in detail.
Good luck. 
Re: Role of divergence
Heyy man, thanks a lot for yr suggestions...I will search if there is any such kind of paper dealing with the issue..
thanks...chandra 
All times are GMT 4. The time now is 02:00. 