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

why doesnt my code work for small Re No??

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 1, 2000, 02:13
Default why doesnt my code work for small Re No??
  #1
lostsoul
Guest
 
Posts: n/a
hi all. i have written a 2d ns solver using the projection method for pressure in the primitive variables. i can run my simulation around a square for Re>100 and it shows the karman vortices. however for Re<100 the code gives all garbage! could somebody enlighten me on this? thanks.. lost soul.

  Reply With Quote

Old   December 1, 2000, 05:04
Default Re: why doesnt my code work for small Re No??
  #2
Fabian
Guest
 
Posts: n/a
Hi Lostsoul (hope this is not your state),

your question is the well known Stokean paradoxon of creeping flow. In a limited flowfield there is no solution for Re->0 as the flowfield has a singularity at infinit distance from the cylinder. See for instance: Kaplun, Lagerstrom "Asymtotic expansions of Navier Stokes solutions for small reynolds numbers" J. Math. Mech. Vol. 6, 585-593. And if you know German: Gersten, Herwig "Stroemungsmechanik", Vieweg Verlag gives a brief introduction to this problem.
  Reply With Quote

Old   December 1, 2000, 07:22
Default Re: why doesnt my code work for small Re No??
  #3
a less-lost soul
Guest
 
Posts: n/a
hi fabian! thanks a lot for the feedback.i am surely getting back my lost confidence well.. umm. i shall surely refer the paper that u have mentioned. btw. do u have any idea how i can make my code work for small Re No? i want to simulate the laminar flow around the square where are no karman vortices. my code starts garbaging below 60!! thanks a lot for all. less-lost soul
  Reply With Quote

Old   December 1, 2000, 10:30
Default Re: why doesnt my code work for small Re No??
  #4
Fabian
Guest
 
Posts: n/a
Hi less-lost soul,

I am not actually an expert in this topic, but for my knoledge you should try to put your outlet boundary as far as possible from the square (in theory at infinite distance). Good luck. But it may be better looking into the theory about creeping flow.
  Reply With Quote

Old   December 1, 2000, 13:04
Default Re: why doesnt my code work for small Re No??
  #5
Tony Chen
Guest
 
Posts: n/a
I guess there are probably two ways to deal with low Re in the projection method:

(1) in the 1st step of projecting the velocity, solving the viscous terms implicitly, i.e. mv the viscous terms to the LHS;

(2) Or just simply reduce the time increment.

Let me know if these work. Good luck.

Tony
  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
SIMPLE (or SIMPLER) fortran code Francesco Di Maio Main CFD Forum 5 March 6, 2014 01:54
From 2D compressible code to 3D code David Liu Main CFD Forum 22 June 26, 2012 17:59
Fast Polygon Triangulation with source code. Liang, Wu Main CFD Forum 0 November 2, 2005 20:42
user friendly cfd code waqar Main CFD Forum 19 August 18, 2000 16:31
own Code vs. commercial code Bernhard Mueck Main CFD Forum 10 February 16, 2000 11:07


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