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

Lid Driven Cavity - Finite Differences BCs?

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 1, 2016, 13:49
Default Lid Driven Cavity - Finite Differences BCs?
  #1
New Member
 
Join Date: Dec 2015
Posts: 4
Rep Power: 10
aem0512 is on a distinguished road
Hello all,

I thought I'd set up the steady-state LDC problem correctly using finite differences, but upon further review...I'm pretty sure there's something wrong here.

I'm using finite differences with a staggered grid to discretize the [0,1] x [0,1] square. I've got dirichlet BC's all the way around except on the lid where u = 1, v = 0.

I've prescribed "fictitious" boundary nodes around the boundary in cases where u or v doesn't fall on the boundary. I've attached a figure which shows the values I have given to the fictitious boundary nodes. The idea is that the average of the inside node and outside node would satisfy the BCs.

The problem is that the vortex seems to grow wildly with the Reynold's number. The figures I'm seeing elsewhere don't do this. I've attached to figures that show what I'm getting. One with Re = 500 and the other with Re = 1000.

Any ideas? I'd really appreciate some input. I've been beating my head against a wall for too long with this already.

Thanks!
Attached Images
File Type: png BCs2.png (41.1 KB, 24 views)
File Type: jpg re1000.jpg (44.1 KB, 27 views)
File Type: jpg re500.jpg (44.6 KB, 22 views)
aem0512 is offline   Reply With Quote

Old   May 1, 2016, 18:39
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,838
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
I need some further infos..

1) you are using the u,v,p formulation, the node for pressure is the green circle, u is the red and v the blue, right?
2) you are using linear extrapolation for the velocity values in the ghost cells, right? therefore the values for u on the top row are uex=-uint+2*uwall
3) the BC.s enters only by the diffusion of u, what if you try the case Re=1?
4) what about the threshold for the steady state convergence?
FMDenaro is offline   Reply With Quote

Old   May 1, 2016, 20:18
Default
  #3
New Member
 
Join Date: Dec 2015
Posts: 4
Rep Power: 10
aem0512 is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
I need some further infos..

1) you are using the u,v,p formulation, the node for pressure is the green circle, u is the red and v the blue, right?
2) you are using linear extrapolation for the velocity values in the ghost cells, right? therefore the values for u on the top row are uex=-uint+2*uwall
3) the BC.s enters only by the diffusion of u, what if you try the case Re=1?
4) what about the threshold for the steady state convergence?
1) yes, that's correct

2) yes, I believe so. For the top fictitious boundary nodes I have, uex = 2 - uint so that when I average (uex + uint)/2 = 1 = u for lid (top wall). For all other walls I have uex = - uint or vex = -vint as appropriate.

3) I have attached a figure showing the case Re=1. The problem only occurs as I increase the Reynolds number. But, it shouldn't take over the whole bottom portion of my figure, right?

4) The tolerance for my linearization (Picard Iteration) convergence is set to 1e-7.

Thanks so much for your input! I've really got to figure this one out, pronto.
Attached Images
File Type: jpg re1.jpg (32.1 KB, 14 views)
aem0512 is offline   Reply With Quote

Old   May 1, 2016, 22:50
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,838
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Re=1 seems correct, i suggest:
- check the Div v=0 constraint in each cell, for staggered grid and second order discretization it must be zero at machine accuracy
- check the solution for refined grids, at higher Re, the non linear terms require a good grid
- check again the tolerance for the steady state, it Could be necessary to decrease the value
FMDenaro is offline   Reply With Quote

Reply

Tags
boundary conditions, finite differences, lid driven cavity, navier stokes

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Lid driven Cavity suman91 SU2 3 May 24, 2022 05:23
lid driven cavity varying results yasmil OpenFOAM Running, Solving & CFD 2 October 6, 2016 21:42
finite volume method - lid driven cavity ainal Main CFD Forum 0 May 15, 2011 11:06
is there any parallel code for the famous Lid Driven Cavity flow? gholamghar Main CFD Forum 0 August 1, 2010 01:55
Lid driven cavity flow KK FLUENT 1 December 16, 2009 18:10


All times are GMT -4. The time now is 06:19.