CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

SIMPLEC Method

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 7, 2009, 09:13
Question SIMPLEC Method
  #1
New Member
 
S.M. Hadi Meshkati
Join Date: Nov 2009
Posts: 7
Rep Power: 7
meshkati is on a distinguished road
Hey everyone,
I try to develop a CFD code by solving 3D Navier-Stocks equations in curvilinear coordinates. But I have a problem with coupling of pressure and velocity (pressure correction equation - pí) that I use SIMPLEC method for it.
In initial iterates (that no pressure gradients exists in domain) the momentum equations have no source term and then: ap=sigma(a_neighbors).
Therefore the parameter of [ap - sigma(a_neighbors)] that appears in the denominator of SIMPLEC method equals to zero and then the divides by zero error was occurred.
What can I do?
Thanks for advance.
Hadi
(Some of my code specifications are: Finite volume method, non -staggered grid, incompressible fluid, steady calculation, no free surface calculation, Rhie and Chow interpolation method, laminar or k-e turbulence model, SIMPLEC pressure correction method and the test case is straight channel with uniform velocity boundary condition in its inlet and if any another comments was needed just tell me)
meshkati is offline   Reply With Quote

Old   November 7, 2009, 21:22
Default
  #2
Senior Member
 
ztdep's Avatar
 
p ding
Join Date: Mar 2009
Posts: 287
Rep Power: 9
ztdep is on a distinguished road
Send a message via Yahoo to ztdep Send a message via Skype™ to ztdep
Ap is the value after under-relaxation. for example , Ap/0.6 Ap/0.6-Sigam (A neighbour) >0
ztdep is offline   Reply With Quote

Old   November 8, 2009, 07:03
Default
  #3
Member
 
Mohammad Reza Hadian
Join Date: Mar 2009
Location: Yazd, Iran
Posts: 50
Rep Power: 8
hadian is on a distinguished road
also, it is recommended to add Delta_F (Fe-Fw+Fn-Fs+Ft-Fb) to Ap.
hadian is offline   Reply With Quote

Old   November 14, 2009, 07:48
Default
  #4
New Member
 
S.M. Hadi Meshkati
Join Date: Nov 2009
Posts: 7
Rep Power: 7
meshkati is on a distinguished road
Quote:
Originally Posted by ztdep View Post
Ap is the value after under-relaxation. for example , Ap/0.6 Ap/0.6-Sigam (A neighbour) >0
That's right, but in SIMPLEC method we can assume all under relaxation factors equals to one (for faster convergence). Then ap equals to ap/urf.
meshkati is offline   Reply With Quote

Old   November 15, 2009, 03:47
Default
  #5
Senior Member
 
ztdep's Avatar
 
p ding
Join Date: Mar 2009
Posts: 287
Rep Power: 9
ztdep is on a distinguished road
Send a message via Yahoo to ztdep Send a message via Skype™ to ztdep
Quote:
Originally Posted by meshkati View Post
That's right, but in SIMPLEC method we can assume all under relaxation factors equals to one (for faster convergence). Then ap equals to ap/urf.
" in SIMPLEC method we can assume all under relaxation factors equals to one "
i do not agree with you. if a 1 is used. the simplec will destinated to divergence
ztdep is offline   Reply With Quote

Old   November 18, 2009, 22:58
Default
  #6
New Member
 
yu
Join Date: Nov 2009
Posts: 26
Rep Power: 7
universez is on a distinguished road
Quote:
Originally Posted by hadian View Post
also, it is recommended to add Delta_F (Fe-Fw+Fn-Fs+Ft-Fb) to Ap.
what is meaning of Delta_F (Fe-Fw+Fn-Fs+Ft-Fb) ?
universez is offline   Reply With Quote

Old   November 19, 2009, 01:48
Default
  #7
Member
 
Mohammad Reza Hadian
Join Date: Mar 2009
Location: Yazd, Iran
Posts: 50
Rep Power: 8
hadian is on a distinguished road
Delta_F=Fe-Fw+Fn-Fs+Ft-Fb and these are convective fluxes passes control volume faces. actually this is the continuity equation and should be zero but it is not zero during the iterations. have a look on the book by Versteeg and Malalaskera.
also, note that in unsteady form of equation you have a A0 which added to Ap. so that Ap-SUM(Anb) will not be zero. moreover a relaxation factor is recommended.
hadian is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Best front tracking method for surface evaporation Komkamol Main CFD Forum 0 July 17, 2009 10:23
Code for most powerfull FDV Method D.S.Nasan Main CFD Forum 6 September 4, 2008 02:08
Help about the preconditioner in gmres method Dan Gao Main CFD Forum 3 July 5, 2008 01:18
hess-smith method and fvm method yangqing FLUENT 0 March 20, 2002 20:25
Is SIMPLEC method suitable for turbomachinery? lhb Main CFD Forum 14 February 15, 2001 20:59


All times are GMT -4. The time now is 11:21.