# 2D lid-driven cavity flow

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

 July 5, 2005, 17:40 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.

 July 5, 2005, 21:14 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.

 July 6, 2005, 09:33 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

 July 6, 2005, 16:13 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.

 July 7, 2005, 07:27 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

 July 26, 2005, 02:51 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.

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

 July 27, 2005, 09:02 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

 July 28, 2005, 09:42 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

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post illuminati5288 Main CFD Forum 0 August 12, 2011 22:05 CFDtoy Main CFD Forum 2 February 19, 2011 10:43 reeshav Main CFD Forum 1 April 12, 2010 01:51 Ravi Shankar Main CFD Forum 0 May 17, 2006 08:37 billy OpenFOAM Running, Solving & CFD 6 July 13, 2005 15:44

All times are GMT -4. The time now is 22:13.