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

Problems obtaining the correct contours for Lid driven cavity

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 3, 2016, 12:56
Default Problems obtaining the correct contours for Lid driven cavity
  #1
DSS
New Member
 
Join Date: Feb 2016
Posts: 6
Rep Power: 10
DSS is on a distinguished road
Hi

I am solving the problem of lid driven cavity in C++. I have used the Finite volume discretization and I have used the Semi Explicit method for solving the problem. However, I am obtaining results that do not match with that of the ones already existing.

The results I plotted are about symmetric about the central vertical line where in that should not be the case. I tried running the code for a finer mesh and higher reynolds number only to encounter the same problem.

I checked my boundary conditions and they are as follows.

u, v = 0 at all the walls except at the top where u = 1
dp/dn = 0 at all the walls
dp'/dn = 0 at all the walls
p' = pressure correction


I have attached the U and V velocity contours and also the pressure contours that I have obtained.

Can someone help me figure out the problem? Any help will be appreciated
Attached Images
File Type: jpg pressure_smaple.jpg (24.8 KB, 12 views)
File Type: jpg u_velocity_smaple.jpg (28.9 KB, 15 views)
File Type: jpg v_velocity_smaple.jpg (30.8 KB, 11 views)
DSS is offline   Reply With Quote

Old   February 3, 2016, 12:58
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
what about the Re number? what main differences you found from the Ghia & Ghia benchmark?
FMDenaro is offline   Reply With Quote

Old   February 3, 2016, 13:28
Default
  #3
DSS
New Member
 
Join Date: Feb 2016
Posts: 6
Rep Power: 10
DSS is on a distinguished road
I carried out this simulation at reynolds number 100. I tried for various reynolds number but still got similar plots for all the cases.

I have attached the correct solution of the problem which is also the same as that given by Ghia.

I am unable to figure out the problem causing the u velocity contours to be symmetric about the central vertical line. As compared to the existing result, the nature is same only that the contours are not inclined towards the corner in my solution.
Attached Images
File Type: jpg images.jpg (11.8 KB, 8 views)
File Type: png V&tc_liddriven_fluentu.png (10.4 KB, 11 views)
DSS is offline   Reply With Quote

Old   February 3, 2016, 13:31
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
if your solution is at Re=100, I see that the contour pattern is too specular as respect with the centerline x=L/2. That seems the convective does not enter into the solution. Are you sure you are gettin a steady state solution?
FMDenaro is offline   Reply With Quote

Old   February 3, 2016, 13:59
Default
  #5
DSS
New Member
 
Join Date: Feb 2016
Posts: 6
Rep Power: 10
DSS is on a distinguished road
Yes, this is the steady state solution that I have obtained. I checked my convection flux equation a lot of times and couldn't find any error in that.

I also compared the magnitudes of the convection and diffusion fluxes and found that the diffusion flux was comparatively larger than the convection flux. But I can't figure out any solution to this since the error could be anywhere (even in my pressure correction equations or somewhere else)

I cross-checked all my equations and formulas and couldn't get the correct results yet.

I think I am missing out on something but I can't figure out what
DSS is offline   Reply With Quote

Old   February 3, 2016, 14:06
Default
  #6
DSS
New Member
 
Join Date: Feb 2016
Posts: 6
Rep Power: 10
DSS is on a distinguished road
Yes, this is the steady state solution that I have obtained. I checked my convection flux equation a lot of times and couldn't find any error in that.

I also compared the magnitudes of the convection and diffusion fluxes and found that the diffusion flux was comparatively larger than the convection flux. But I can't figure out any solution to this since the error could be anywhere (even in my pressure correction equations or somewhere else)

I cross-checked all my equations and formulas and couldn't get the correct results yet.

I think I am missing out on something but I can't figure out what
DSS is offline   Reply With Quote

Old   February 3, 2016, 14:10
Default
  #7
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by DSS View Post
Yes, this is the steady state solution that I have obtained. I checked my convection flux equation a lot of times and couldn't find any error in that.

I also compared the magnitudes of the convection and diffusion fluxes and found that the diffusion flux was comparatively larger than the convection flux. But I can't figure out any solution to this since the error could be anywhere (even in my pressure correction equations or somewhere else)

I cross-checked all my equations and formulas and couldn't get the correct results yet.

I think I am missing out on something but I can't figure out what
Using non-dimensional equations, the convective flux should be O(1), compared to the O(10^-2) of the diffusive one, irrespective of the pressare term.
FMDenaro is offline   Reply With Quote

Old   February 3, 2016, 14:24
Default
  #8
DSS
New Member
 
Join Date: Feb 2016
Posts: 6
Rep Power: 10
DSS is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Using non-dimensional equations, the convective flux should be O(1), compared to the O(10^-2) of the diffusive one, irrespective of the pressare term.
That is exactly what I have not been able to solve. My convective flux is O(10^-3) of the diffusive one. Hence, the convective flux doesn't contribute much to the equations.
DSS is offline   Reply With Quote

Old   February 10, 2016, 14:55
Default
  #9
DSS
New Member
 
Join Date: Feb 2016
Posts: 6
Rep Power: 10
DSS is on a distinguished road
I'd really appreciate if anyone could help me with this problem. I want to know what could be the possible errors so that I can work on them.
DSS is offline   Reply With Quote

Old   February 10, 2016, 18:06
Default
  #10
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
As I previously wrote, the symmetry could be due to an error in the discretization of the convective terms
FMDenaro is offline   Reply With Quote

Reply

Tags
boundary condition, boundary condition p, finite volume method, lid driven cavity in c, lid- driven cavity


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
[ICEM] Problems with coedge curves and surfaces tommymoose ANSYS Meshing & Geometry 6 December 1, 2020 12:12
driven cavity code amir Main CFD Forum 9 April 17, 2014 12:19
Lid driven cavity ganesh Main CFD Forum 0 August 5, 2009 21:33
Traditional driven cavity Azfarizal b Mukhtar CFX 0 July 27, 2004 22:22
2D Driven square cavity problem rvndr Main CFD Forum 6 February 25, 2004 11:35


All times are GMT -4. The time now is 16:28.