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

Lid Driven Cavity

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 6, 2020, 09:07
Default Lid Driven Cavity
  #1
New Member
 
Mudit Mehta
Join Date: Mar 2020
Posts: 15
Rep Power: 2
Mudit2208 is on a distinguished road
I am working on Lid driven cavtiy problem using SIMPLE Method in python. My code is working and giving correct results for reynolds number 100. But when i run it for reynolds number 1000, it gives this:

RuntimeWarning: invalid value encountered in double_scalars:
vp[i, j] = dv[i, j] * (pp[i, j] - pp[i, j + 1])

RuntimeWarning: overflow encountered in double_scalars
i, j] * pp[i, j - 1]) + pb[i, j])) / pap[i, j]

and i have no clue why. Please help
Mudit2208 is offline   Reply With Quote

Old   May 6, 2020, 13:09
Default
  #2
New Member
 
Join Date: Feb 2017
Posts: 6
Rep Power: 5
donatello is on a distinguished road
What you are seeing is likely a divide by 0 error, probably from your iterative solver. Maybe take a look at how your residuals are computed.
donatello is offline   Reply With Quote

Old   May 6, 2020, 15:20
Default
  #3
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 5,551
Rep Power: 58
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Mudit2208 View Post
I am working on Lid driven cavtiy problem using SIMPLE Method in python. My code is working and giving correct results for reynolds number 100. But when i run it for reynolds number 1000, it gives this:

RuntimeWarning: invalid value encountered in double_scalars:
vp[i, j] = dv[i, j] * (pp[i, j] - pp[i, j + 1])

RuntimeWarning: overflow encountered in double_scalars
i, j] * pp[i, j - 1]) + pb[i, j])) / pap[i, j]

and i have no clue why. Please help



I suggest to consider increasing the grid resolution according to the Re number. But go step by step, solve Re=100, 200, 400 and check the solutions
FMDenaro is offline   Reply With Quote

Reply

Tags
cfd code, python script, simple algorithm

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
Setting up Lid driven Cavity Benchmark with 1M cells for multiple cores puneet336 OpenFOAM Running, Solving & CFD 11 April 7, 2019 01:58
Boundary Conditions at Corners, lid driven cavity t.teschner Main CFD Forum 8 August 17, 2016 09:20
Lid driven Cavity suman91 SU2 2 October 20, 2014 16:38
is there any parallel code for the famous Lid Driven Cavity flow? gholamghar Main CFD Forum 0 August 1, 2010 02:55
Lid driven cavity flow KK FLUENT 1 December 16, 2009 19:10


All times are GMT -4. The time now is 08:50.