|
[Sponsors] |
Help with coding SIMPLE for Lid-Driven Cavity |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
August 12, 2013, 10:19 |
Help with coding SIMPLE for Lid-Driven Cavity
|
#1 |
New Member
Sam Gregson
Join Date: Jul 2013
Posts: 4
Rep Power: 13 |
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 lid-driven 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, 06:06 |
|
#2 |
New Member
Sam Gregson
Join Date: Jul 2013
Posts: 4
Rep Power: 13 |
UPDATE:
I am currently testing my code using Poiseuille flow between two parallel plates (following advice by HouKen in this post: http://www.cfd-online.com/Forums/mai...ng-simple.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, 12: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, 07:11 |
|
#4 |
New Member
Sam Gregson
Join Date: Jul 2013
Posts: 4
Rep Power: 13 |
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 under-relaxation though. I think I misunderstood the under-relaxation of the momentum equations. I applied under-relaxation by using the SOR method. However this under-relaxes the process of solving the system of equations rather than under-relaxing the solution itself. Is this correct?
|
|
August 19, 2013, 08:06 |
|
#5 |
Senior Member
|
yes , you are right. ap/alpha*U=sigma*Uab+Source+(1-alpha)/alpha*U*
in my blog in this web, there is a fortran code for the simple algorithm |
|
August 19, 2013, 11:22 |
|
#6 |
New Member
Sam Gregson
Join Date: Jul 2013
Posts: 4
Rep Power: 13 |
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, 23: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 21:29 |
Lid Driven Cavity using Ghost Cell Method and in C++ | illuminati5288 | Main CFD Forum | 0 | August 12, 2011 23:05 |
Traditional driven cavity | Azfarizal b Mukhtar | CFX | 0 | July 27, 2004 22:22 |
SIMPLE Cavity Flow | Kunal Jain | Main CFD Forum | 1 | August 26, 2001 12:43 |
2D driven cavity and chaos | Abhijit Joshi | Main CFD Forum | 6 | February 1, 1999 14:51 |