|
[Sponsors] |
November 6, 2002, 18:58 |
SIMPLER boundary condition
|
#1 |
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. |
|
November 6, 2002, 22:20 |
Re: SIMPLER boundary condition
|
#2 |
Guest
Posts: n/a
|
The header of your posting appears to have the answer - check outlet BC.
|
|
November 6, 2002, 23:43 |
Re: SIMPLER boundary condition
|
#3 |
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!
|
|
November 7, 2002, 03:39 |
Re: SIMPLER boundary condition
|
#4 |
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. |
|
November 7, 2002, 13:23 |
Re: SIMPLER boundary condition
|
#5 |
Guest
Posts: n/a
|
The flow is laminar, Re=400. Thanks for the help.
|
|
November 8, 2002, 00:24 |
Re: SIMPLER boundary condition
|
#6 |
Guest
Posts: n/a
|
Could you mail me your code for reference? Thx!
|
|
November 8, 2002, 02:43 |
Re: SIMPLER boundary condition
|
#7 |
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 |
|
November 8, 2002, 12:06 |
Re: SIMPLER boundary condition
|
#8 |
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 |
|
November 8, 2002, 13:56 |
Re: SIMPLER boundary condition
|
#9 |
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. |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Boundary Conditions | Thomas P. Abraham | Main CFD Forum | 20 | July 7, 2013 05:05 |
inlet velocity boundary condition | murali | CFX | 5 | August 3, 2012 08:56 |
External Radiation Boundary Condition for Grid Interface | CFD XUE | FLUENT | 0 | July 9, 2010 02:53 |
problem with boundary condition??? | smn | CFX | 5 | November 24, 2009 06:37 |
RPM in Wind Turbine | Pankaj | CFX | 9 | November 23, 2009 04:05 |