
[Sponsors] 
March 6, 2001, 02:09 
SMAC

#1 
Guest
Posts: n/a

Hi,
While using SMAC, only velocity B.C need be specified. However i am observing that after convergence, the freestrem pressure which i specify once, only at the begining, is not the same. Here I am using the free stream vel as the B.C. What could be the problem? Also are there any special considerations for the ficticious boundary cells? (the simulation is for the unconfined flow over a cylinder) Any advice would be appreciated. Does anyone know where an electronic copy of papers relating to SMAC etc can be found. Thank you. 

March 6, 2001, 10:09 
Re: SMAC

#2 
Guest
Posts: n/a

Teeda,
Sounds as if you're getting a good solution. For an incompressible flow, only the pressure gradient is important. This is not due to the algorithm you have chosen (SMAC or otherwise). So, during the iterative solution of the pressure (continuity) equation used in SMAC, the pressure level will drift somewhat. However, if the iteration converges, the pressure gradient is correct and the continuity equation will be satisfied. If you need to fix the pressure at some point(s) on the boundary (trying to match an experiment, for example), that is permitted. But the velocity boundary condition at that point will be overridden. 

March 8, 2001, 13:17 
Re: SMAC

#3 
Guest
Posts: n/a

Dear Jim
Thanx, your reply was helpful. I have some more questions though, 1. during the convergence of the simulation, the maximum cell div did not drop below 0.002...it kept oscillating about this value. However the results at this stage were fairly accurate..why isnt the max div dropping further? 2. in your last reply u said the vel b.c will be over ridden...could u please elaborate Thanx again 

March 8, 2001, 13:50 
Re: SMAC

#4 
Guest
Posts: n/a

I'll take a shot at answering.
First I've used the original SMAC code as distributed by the T3 group at Los Alamos. More recently the SOLAVOF code, same source. These codes [in fact most of those developed by T3 in the late 60's, 70's and early 80's] are remarkably similar in the strategy used to implement boundary conditions. A ring of 'image' cells or 'ghost' cells is used and appropriate values are set in or on the appropriate edges of those cells to impose the desired boundary condition. For velocity, it's a specified normal velocity (including noflow, the zero value), AND the tangential velocity (no slip or free slip  zero shear). When the pressure is specified at one (or more) boundary location, it is set in the first 'real' cell and that cell is excluded from updating in the momentum and pressure equation loops. So the velocity values set in the image cells just outside the cell where the pressure is set are never referenced and thus ignored. The pressure condition replaces the velocity condition at that point in the solution domain. As to your first question regarding convergence, I can very little. This is always a black art 

March 8, 2001, 14:28 
Re: SMAC

#5 
Guest
Posts: n/a

Dear Jim
Thanx for the reply. I hope u wont mind if I bother u with another question. my problem is the simulation of the flow over a cyl. for one edge of the domain, the bc is no slip and for the other the b.c is free stream vel. The grid is staggered. My ques is as follows: while implementing noslip, say, consider an image cell. This would contain 2 tangential vel components and two vert. vel. comp. The no slip boundary lies along one vert comp, therefore this is set to zero. As for the tangential components, they r combined with the corresponding tang comp in the real cell just above it such that their avg is zero. What about the 2nd vert. comp. of vel. in the image cell. Is there any condition that should be imposed on it? Do u have an electronic copy of the SMAC paper? Thanks. 

March 8, 2001, 18:15 
Re: SMAC

#6 
Guest
Posts: n/a

Sorry, I've not dealt with irregular boundary situation you describe.
Also, I don't have a copy of the SMAC code. Left that at work when I retired. 

March 9, 2001, 13:45 
Re: SMAC

#7 
Guest
Posts: n/a

teeda,
Concerning your first question about the limited convergence of the iterations: one possibility is that you are encountering roundoff errors, especially if the code is being run with 32bit precision. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Periodic pressure B.C. on collocated SMAC code  waku2005  Main CFD Forum  0  August 18, 2010 23:03 
Looking for the original SMAC article  Italo  Main CFD Forum  0  September 2, 2009 13:26 
About Poisson equation's B.C. in SMAC algorithm  snakejoe  Main CFD Forum  0  April 7, 2009 05:58 
what is difference of SMAC, project,fraction step  qunwuhe@hotmail.com  Main CFD Forum  0  November 1, 2007 02:15 
SMAC  Niru  Main CFD Forum  0  May 23, 2003 01:03 