# SMAC

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 LinkBack Thread Tools Display Modes
 March 6, 2001, 02:09 SMAC #1 Tee Guest   Posts: n/a Sponsored Links 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.
 Sponsored Links

 March 6, 2001, 10:09 Re: SMAC #2 Jim Park 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 teeda 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 Jim Park Guest   Posts: n/a I'll take a shot at answering. First I've used the original SMAC code as distributed by the T-3 group at Los Alamos. More recently the SOLA-VOF code, same source. These codes [in fact most of those developed by T-3 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 no-flow, 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 teeda 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 no-slip, 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 Jim Park 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 Dean Guest   Posts: n/a teeda, Concerning your first question about the limited convergence of the iterations: one possibility is that you are encountering round-off errors, especially if the code is being run with 32-bit precision.

 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 waku2005 Main CFD Forum 0 August 18, 2010 23:03 Italo Main CFD Forum 0 September 2, 2009 13:26 snakejoe Main CFD Forum 0 April 7, 2009 05:58 qunwuhe@hotmail.com Main CFD Forum 0 November 1, 2007 02:15 Niru Main CFD Forum 0 May 23, 2003 01:03

 Sponsored Links

All times are GMT -4. The time now is 22:54.

 Contact Us - CFD Online - Top