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

Incompressible Navier-Stokes in Stream Function - Vorticity Formulation

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 15, 2012, 10:20
Default Incompressible Navier-Stokes in Stream Function - Vorticity Formulation
  #1
New Member
 
Join Date: Aug 2012
Posts: 9
Rep Power: 4
Daco is on a distinguished road
Hello,

I am trying to solve the incompressible, steady state Navier-Stokes in stream function - vorticity formulation, for a given domain (attached pictures).

for the boundary conditions of the vorticity I use Thompson's condition (1930).

I am solving it using an iterative method I developed from the differential equations. I use central finite differences schemes (attached as well).
I solve the vorticity equation - than the stream function equation - than calculate velocities - check convergence - if not, do it again.

I use successive under relaxation for the vorticity for without it the program diverges after a few iterations.

My problem is that as I increase the Reynolds number my solution becomes extremely unstable and it diverges. let say for Re=100 I need only a relaxation factor of 0.5 for the vorticityand the program converges - I get physical results and all is well. But for Re=700 It doesn't converge even for 0.01...

I tried formulating my equation with upwind-schemes but it seemed to only worsen my stability.

Does anyone had this experience with this problem? how do I overcome it?

It seems my problem comes from the convection terms in the equation (the ones Re multiplies)


I am attaching my code that is written in matlab to who ever want to take a look.
use main.m to start it. change Re number if you need.
relax1 - relaxation factor for the vorticity
relax2 - -||- stream function
relax3 - -|| - velocities
it is set now on Re=100. and the mesh is pretty coarse but still gets nice results. sorry for not writing titles on the figures.. u can check the code to see what's in them.

Thanks,

Daniel
Attached Images
File Type: jpg pic.jpg (50.6 KB, 30 views)
File Type: jpg figure1.jpg (22.5 KB, 32 views)
Attached Files
File Type: zip N-S_incomp_steady.zip (9.2 KB, 36 views)
Daco is offline   Reply With Quote

Old   October 9, 2012, 16:22
Default
  #2
New Member
 
Thyago de Lellys
Join Date: May 2011
Posts: 1
Rep Power: 0
thylellys is on a distinguished road
Hi Daniel,

in fact .... why do you use FDM? im using another meshless method. maybe FVM will be better than FDM in your case.
thylellys is offline   Reply With Quote

Old   October 9, 2012, 16:38
Default
  #3
New Member
 
Join Date: Aug 2012
Posts: 9
Rep Power: 4
Daco is on a distinguished road
The reason is that we didn't cover FVM in our course. Only FDM

but thanks though.

If I'm here I'll update that I solved the problem using very strong relaxation (0.05). It's a banging your head against the wall solution but it worked and I was out of time so it had to do.

Thanks,

Daniel
Daco is offline   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
Coding 2d incompressible Navier Stokes Simon Main CFD Forum 15 July 20, 2012 04:24
ParaView for OF-1.6-ext Chrisi1984 OpenFOAM Installation 0 December 31, 2010 07:42
Stream function in axisymmetrical problem Daniel CD-adapco 0 September 8, 2003 11:20
Newbie:Viscoelasticity and Navier stokes equation Rajil Saraswat Main CFD Forum 2 June 9, 2003 07:21
Compressible vs. Incompressible formulations Fernando Velasco Hurtado Main CFD Forum 3 January 7, 2000 17:51


All times are GMT -4. The time now is 14:55.