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

2D lid-driven cavity flow

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 5, 2005, 17:40
Default 2D lid-driven cavity flow
  #1
Márcio Ricardo
Guest
 
Posts: n/a
Hi, I have a problem when simulating a 2D lid-driven cavity flow at Re=100: the velocity profiles at the center lines and the center of the primary vortex match the literature, but the secondary vortices don't. For a lid velocity Ul = +1.0, the vortex located at the right lower corner is smaller than the vortex located at the left lower corner, which could mean that my formulation doesn't obey the momentum conservation. Does anybody knows why this could happen? Also, when I plot the streamlines, they seem to cross the walls when very close to these corners.

Here's my formulation:

>Fractional step for u-v coupling

>1st order - explicit time advance (dt = 1e-4)

>Re = 100

>Simplified viscous term: Dij = d(visc*(dui/dxj))/dxj

>no turbulence modeling

>Finite Element for spatial discretization (quadrilateral elements): 9 node quadratic lagrange element for velocity and 4 node linear element for pressure. The weighting functions are constant and define a control volume around each node (CVFEM). This leads to non-symmetric matrices for the velocity element.

>Lumped mass formulation for the transient term: all the terms in the transient matrix are summed and stored in each row of the matrix. The remaining columns in each row are set to zero.
  Reply With Quote

Old   July 5, 2005, 21:14
Default Re: 2D lid-driven cavity flow
  #2
Jim_Park
Guest
 
Posts: n/a
How about your nodal spacing? What is the nodal pattern? Are your quads rectangles for the test problem? In other words, how is your flow regime resolved?

I have run this problem (Re = 100) using a finite difference code. The corner vortices at 20 x 20 are much different than those for 40 x 40. Those solutions are for uniform mesh spacing. Going to 50 x 50 with the meshes stretched in the center to allow very close spacing in the corners shows even further change in the corners. For the 40 x 40 and 50 x 50, the centerline velocity profiles match the Ghia, Ghia, and Chin solutions (usually regarded as the standards) very well. I also checked the centerline velocities using CFX - with good agreement also.

These results are in a report you can download from www.ornl.gov. The report number is ORNL/TM-1999/292. The author is James E. Park.
  Reply With Quote

Old   July 6, 2005, 09:33
Default Re: 2D lid-driven cavity flow
  #3
Márcio Ricardo
Guest
 
Posts: n/a
Hi, Jim. Thanks for your help. My results agree with the report, except for the secondary vortices, and my grid has 71 x 71 points. Before that I had run a case with an unstructured grid even finer at the walls, and the same problem occured. Do you think there could be an error in the spatial discretization or in the fractional step algorithm?

TIA

Márcio Ricardo
  Reply With Quote

Old   July 6, 2005, 16:13
Default Re: 2D lid-driven cavity flow
  #4
Jim_Park
Guest
 
Posts: n/a
I'm really out of ideas.

The possibility of insufficient resolution was what your results suggested. But I'd think 71 x 71 is enough if the other aspects are OK. By that I mean code correct, time (?) step within stability bounds, etc.

I have no experience with fractional step, or with finite element, so can't help at all there.

Sorry I wasn't of any particular help.
  Reply With Quote

Old   July 7, 2005, 07:27
Default Re: 2D lid-driven cavity flow
  #5
Márcio Ricardo
Guest
 
Posts: n/a
Thank you, Jim. I'm checking the code for the advective term and the fractional step algorithm.

Thanks for your help

Márcio Ricardo
  Reply With Quote

Old   July 26, 2005, 02:51
Default Re: 2D lid-driven cavity flow
  #6
ALI RAJABPOUR
Guest
 
Posts: n/a
Would you please send me the pressure distribution in a cavity flow at low and high Re.
  Reply With Quote

Old   July 27, 2005, 00:58
Default Re: 2D lid-driven cavity flow
  #7
Zhao
Guest
 
Posts: n/a
  Reply With Quote

Old   July 27, 2005, 09:02
Default Re: 2D lid-driven cavity flow
  #8
Marcio Ricardo
Guest
 
Posts: n/a
Try the following site. There's a lot of information - including papers - about flows in different cavity shapes.

http://www.gyte.edu.tr/enerji/ercane...cavityflow.htm
  Reply With Quote

Old   July 28, 2005, 09:42
Default Re: 2D lid-driven cavity flow
  #9
Ali Rajabpour
Guest
 
Posts: n/a
Dear Marcio Ricardo,

The address that you have sent, does not have the pressure solution of cavity flow.

Thanks alot

  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
Lid Driven Cavity using Ghost Cell Method and in C++ illuminati5288 Main CFD Forum 0 August 12, 2011 22:05
Scalar transport in Lid driven cavity CFDtoy Main CFD Forum 2 February 19, 2011 09:43
2D LID driven cavity reeshav Main CFD Forum 1 April 12, 2010 01:51
Lid Driven Cavity Problem Ravi Shankar Main CFD Forum 0 May 17, 2006 08:37
InterFoam with lid driven cavity billy OpenFOAM Running, Solving & CFD 6 July 13, 2005 15:44


All times are GMT -4. The time now is 00:15.