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

SIMPLER boundary condition

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 6, 2002, 19:58
Default SIMPLER boundary condition
  #1
Ramesh
Guest
 
Posts: n/a
Hi

I've written a SIPLER code for a 2D flow in a pipe. The transient part is very important for my problem. I have notice the radial velocity for the incopresible flow is nealy zero far from inlet however, its value increases as it approaches the outlet(Say from 10^-18 to 10^-14). Also I can not reach steady state solution, I get close but the term dv/dt does not become zero. Your help on this matter is greatly appreciated.
  Reply With Quote

Old   November 6, 2002, 23:20
Default Re: SIMPLER boundary condition
  #2
mukhopadhyay
Guest
 
Posts: n/a
The header of your posting appears to have the answer - check outlet BC.
  Reply With Quote

Old   November 7, 2002, 00:43
Default Re: SIMPLER boundary condition
  #3
Ramesh
Guest
 
Posts: n/a
Well, according to Patankar it is not necessary to impose and velocity boundary condition at the outlet. Knowing the outlet pressure should be enough. However, I observed the strange behavoir of radial velocity and am not able to reach the steady state!
  Reply With Quote

Old   November 7, 2002, 04:39
Default Re: SIMPLER boundary condition
  #4
STN
Guest
 
Posts: n/a
Could you post the Reynolds number?

Is it in the range of Turbulent flow.

If yes which type of turbulent model, do you use?

These information may indicate your problem much more clear to all CFD man in this forum.
  Reply With Quote

Old   November 7, 2002, 14:23
Default Re: SIMPLER boundary condition
  #5
Ramesh
Guest
 
Posts: n/a
The flow is laminar, Re=400. Thanks for the help.
  Reply With Quote

Old   November 8, 2002, 01:24
Default Re: SIMPLER boundary condition
  #6
seaspace
Guest
 
Posts: n/a
Could you mail me your code for reference? Thx!
  Reply With Quote

Old   November 8, 2002, 03:43
Default Re: SIMPLER boundary condition
  #7
Anthony Wachs
Guest
 
Posts: n/a
Hi Ramesh,

If I read well, you are talking in your mail about radial velocity whose magnitude increases from 10^-18 to 10^-14 ?

My opinion is that you do not have to care about such velocity magnitude of 10^-14 or 10^-18 (or anyway any value less than 10^-8 !!), this is definitively zero for your computer since this is less than the usual machine precision.

What I am trying to say is that 10^-14 for a computer is zero for you. Usually, if you are performing calculations in single precision, than velocity value less than 10^-8 have nearly no meaning (for double precision, the limite is 10^-16) and could be considered to be zero.

I think your code is working very well, there is probably nothing that need to be modified.

You are likely to be a bit too optimistic or to expect too much from your code (I'm just kidding, obviously)

I hope it could help,

Anthony
  Reply With Quote

Old   November 8, 2002, 13:06
Default Re: SIMPLER boundary condition
  #8
Ambrose
Guest
 
Posts: n/a
Well, actually I have very similar problems as Ramesh does. Mine is steady, incompressible, Poiseulle-Couette flow with Re=100. Uniform flow in and fully developed flow out with no information on Pressure; pressure at exit is not given. The problem is that the axial velocity profile at exit is decreasing as the length of the pipe is getting longer; 10% error. One of the things that I suspect is the mass is decreasing from the inlet to exit. Also, have to keep relaxation factor for velocities very small like less than 0.1. Surprisingly, results are really good with 0.01. However, I have no idea why. I looked at my codes very carefully again and again, but everything looks fine. I appreciate any comments.

Ambrose
  Reply With Quote

Old   November 8, 2002, 14:56
Default Re: SIMPLER boundary condition
  #9
Ramesh
Guest
 
Posts: n/a
Hi Anthony,

Thanks for your comments, If that is the case I still have problem to reach steady state solution. As I mentioned the term dvz/dt (z is the axial direction) never becomes zero (again it ossilats from 10^-13 to 1-^-15). Thanks for your help.
  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
Boundary Conditions Thomas P. Abraham Main CFD Forum 20 July 7, 2013 06:05
inlet velocity boundary condition murali CFX 5 August 3, 2012 09:56
External Radiation Boundary Condition for Grid Interface CFD XUE FLUENT 0 July 9, 2010 03:53
problem with boundary condition??? smn CFX 5 November 24, 2009 07:37
RPM in Wind Turbine Pankaj CFX 9 November 23, 2009 05:05


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