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

Need help to solve N-S eqns with F.V.M

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 25, 2008, 08:27
Default Need help to solve N-S eqns with F.V.M
  #1
vasilis
Guest
 
Posts: n/a
Hi,

I am trying to solve the Navier-Stokes eqns in a 3D duct, with the aid of Finite Volume Method. At the inlet, we know the velocity in the x-direction, and at the outlet we assume zero derivative of the velocity with respect to x.

When I descritized the eqns, I considered that there are cv nodes not only in the center of each volume cell, but also on the boundaries. These extra cv nodes helped me to calculate the gradients. For those cv nodes, I impose zero velocity in the x, y and z direction at the walls, zero velocity in the y and z direction at the inlet. Also, at those nodes the pressure is calculated by doing an extrapolation from the adjacent nodes.

Also, because the flow is incompressible, I set the value of the pressure to a certain value somewhere in the domain. The gradients at the cv nodes are calculated with central finite difference, whereas the gradients at the faces are calculated by doing a linear interpolation between the two adjacent cells.

I am solving the eqns following a Newton-Raphson scheme, and I observe that the solution does not converge. The pressure and the velocity obtain some really large values, i.e. 10^(13)

Does anybody have any idea what am I doing wrong?

Thank you for your time

  Reply With Quote

Old   September 25, 2008, 08:45
Default Re: Need help to solve N-S eqns with F.V.M
  #2
vasilis
Guest
 
Posts: n/a
An update, If I increase the number of cells then Newton-Raphson diverges faster, i.e. even at the first iteration the pressure has value of the order 10^(5).

  Reply With Quote

Old   September 26, 2008, 04:52
Default Re: Need help to solve N-S eqns with F.V.M
  #3
vasilis
Guest
 
Posts: n/a
I have tried many different things: I tried to solve the poisson eqn for the pressure, I tried setting the pressure to zero at the inlet or the outlet, but nothing. The pressure keeps acquiring large values.

Does anybody know how can I bound the pressure? What boundary conditions should I use for the pressure?

Thank you
  Reply With Quote

Old   September 26, 2008, 14:16
Default Re: Need help to solve N-S eqns with F.V.M
  #4
andy2o
Guest
 
Posts: n/a
This book by Patankar covers many of the basic causes of non-convergence and is a great guide:

http://www.cfd-online.com/Books/show_book.php?book_id=8

1) How are you coupling the momentum and pressure equations?

2) Have you tried for example the SIMPLE method - it's not quite clear from your posts?

3) Have you tested your linear solvers, etc on some simpler equations? Try a really small system of equations (<10 unknowns even) and compare with results from Matlab or Octave

4) Once you've tested your linear solvers, move on to try a simple PDE. For example, try using your mesh and solver to solve the heat equation. This help help show that there are no code bugs.

Good luck, Andy
  Reply With Quote

Old   October 7, 2008, 02:56
Default Re: Need help to solve N-S eqns with F.V.M
  #5
vasilis
Guest
 
Posts: n/a
Hi,

I am using a collocated mesh. The old code was implementing the SIMPLE method. After some reading I realized that is not possible to solve the eqns on a collocated grid without using the SIMPLE (or a similar) method. If I do not want to use the SIMPLE method, then I need to use a staggered grid.

With respect to your questions about the solver, I have tested it and it works.
  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
Solve Flow or VOF simultaneously ? Ramsey FLUENT 1 February 16, 2011 14:16
Linearized NS euqations: how to solve them?(problem with Matrix operations..) matteoL OpenFOAM Running, Solving & CFD 0 November 18, 2009 07:58
How to model the NR eqns in a domain with empty space Vasilis Main CFD Forum 1 April 14, 2009 04:35
Help to solve 2 eqns simultaneously in MPI Zonexo Main CFD Forum 2 June 1, 2008 20:21
How to solve another continuum and momentum eqn? west_wing FLUENT 0 August 25, 2003 10:00


All times are GMT -4. The time now is 17:31.