
[Sponsors] 
Help with coding SIMPLE for LidDriven Cavity 

LinkBack  Thread Tools  Search this Thread  Display Modes 
August 12, 2013, 09:19 
Help with coding SIMPLE for LidDriven Cavity

#1 
New Member
Sam Gregson
Join Date: Jul 2013
Posts: 4
Rep Power: 11 
Hi all,
I am new to CFD but very interested in the field. I have been trying to write a very simple code to solve the liddriven cavity problem with a staggered grid using a steady SIMPLE method and have encountered some problems. Currently, solving the momentum equations gives reasonable results for the first iteration and the pressure correction does correct for divergence after updating the velocities. However as the SIMPLE iterations progress the solution becomes unrealistic and I notice that the pressure corrections increases wildly with every iteration along with the divergence of each cell. It is difficult to ask a specific question (sorry) as I don't know where the problem lies but maybe people could give me pointers on the implementation or advice on how to check various parts of the code? Here is some info on what I have done:
I would really appreciate any help. Thank you in advance, Sam Here are the x velocities and pressures after 1 iteration (including correcting velocities and ghost cells are included): Here they are after 2: 

August 16, 2013, 05:06 

#2 
New Member
Sam Gregson
Join Date: Jul 2013
Posts: 4
Rep Power: 11 
UPDATE:
I am currently testing my code using Poiseuille flow between two parallel plates (following advice by HouKen in this post: http://www.cfdonline.com/Forums/mai...ngsimple.html). I am reasonably convinced that my momentum equations are in order since they do not change the initial conditions when the precise solution is given as initial conditions. Furthermore when setting arbitrary initial conditions (apart from the inflow velocity) the momentum equations give good answers. This result also gives me confidence that boundary conditions are in order. As for the pressure equation, when given a velocity field it seems to produce a solution that has a better mass balance (for each cell), which is good news. And therefore I am reasonably happy with the pressure correction equations. So, what's the problem???? Well, when I iterate in the outer loop between the momentum equations and the pressure equation the solution just gives nonsense results, the mass imbalance grows (between outer iterations, but diminishes within an iteration) and the momentum equations finally just diverge. What is going on here??? I just don't seem to be getting anywhere. Any tips? I would be very grateful! 

August 18, 2013, 11:04 

#3 
Senior Member

Since you are writing a new code, please make sure your matrix coefficient are right calculated.
Since your procedure diverge, one reason is the underelaxtion is too high. and the other is your matrix coeffcient is wrong. good luck to you. 

August 19, 2013, 06:11 

#4 
New Member
Sam Gregson
Join Date: Jul 2013
Posts: 4
Rep Power: 11 
Thank you for your reply P Ding. Yes I spent a long time staring at the coefficients but can't see anything wrong with them. I think you are right about the underrelaxation though. I think I misunderstood the underrelaxation of the momentum equations. I applied underrelaxation by using the SOR method. However this underrelaxes the process of solving the system of equations rather than underrelaxing the solution itself. Is this correct?


August 19, 2013, 07:06 

#5 
Senior Member

yes , you are right. ap/alpha*U=sigma*Uab+Source+(1alpha)/alpha*U*
in my blog in this web, there is a fortran code for the simple algorithm 

August 19, 2013, 10:22 

#6 
New Member
Sam Gregson
Join Date: Jul 2013
Posts: 4
Rep Power: 11 
Great, thanks. I'm getting somewhere now.
I can't seem to find your blog, do you have a link? I was wondering whether the coefficients for momentum updated in the inner iterations or only in the outer iterations? It seems that my code is slightly more stable when they are updated in the inner iterations but I don't know what the norm is and whether this has any disadvantages? 

September 12, 2013, 22:48 

#7  
Senior Member

Quote:
yes, outer iteration resolve the nonlinearity , it update the coefficients. and inner iteration only solve the linear equation with a constant coefficients 

Tags 
boundary condition, code, lid driven cavity, simple method, staggered grid 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Boundary Condition for Driven Cavity problem  dado  Main CFD Forum  6  August 15, 2015 20:29 
Lid Driven Cavity using Ghost Cell Method and in C++  illuminati5288  Main CFD Forum  0  August 12, 2011 22:05 
Traditional driven cavity  Azfarizal b Mukhtar  CFX  0  July 27, 2004 21:22 
SIMPLE Cavity Flow  Kunal Jain  Main CFD Forum  1  August 26, 2001 11:43 
2D driven cavity and chaos  Abhijit Joshi  Main CFD Forum  6  February 1, 1999 13:51 