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

flow in duct-fixed pressure boundary condition

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 14, 2012, 07:05
Default flow in duct-fixed pressure boundary condition
  #1
New Member
 
maryam
Join Date: Jan 2012
Posts: 11
Rep Power: 5
sepidehkavousi is on a distinguished road
hi every body
I am writing a code for flow in micro channels, but as boundary condition I do not have any information about velocity, I just know I have input pressure 150000 (pa) and output pressure atmosphere (100000 pa). I use finite volume method and simple algorithm.
In Versteeg's book it is defined az physical boundary condition bud I have some problems . I don't get convergence for my code and I think it is something about my boundary condition defining. I wonder if anyone could help me.
my Du and Dv and so for most nodes pressure corrections too small to correct velocity or pressure .
I don't know what to do about under relaxation factor, does it have such important effect on code convergence?
please help me
sepidehkavousi is offline   Reply With Quote

Old   May 14, 2012, 09:01
Default
  #2
New Member
 
sara
Join Date: May 2012
Posts: 5
Rep Power: 5
sara1987 is on a distinguished road
is it steady state?
have you done grid size independence check?
have you initialized it correctly?
can u tell me more details about the problem u wanna solve
sara1987 is offline   Reply With Quote

Old   May 14, 2012, 10:23
Default
  #3
Member
 
Join Date: Jul 2011
Posts: 59
Rep Power: 6
rmh26 is on a distinguished road
Can you assume the flow is undisturbed near the boundaries. ie just a posiuelle type flow. If so then as a reasonable approximation you can set the tangential part of the flow to zero and set the derivative of the normal part of the flow to zero. There are some better approximation you can use if the flow isn't so regular.

"Open and Traction boundary conditions for the incompressible navier stokes equations" Jie Liu JCP 228, 2009
rmh26 is offline   Reply With Quote

Old   May 14, 2012, 11:11
Default
  #4
Senior Member
 
Mihai Pruna
Join Date: Apr 2010
Location: Boston
Posts: 190
Blog Entries: 1
Rep Power: 7
mihaipruna is on a distinguished road
can you try setting gradient limiters?
__________________
Mihai Pruna's Bio
mihaipruna is offline   Reply With Quote

Old   May 18, 2012, 07:51
Default
  #5
New Member
 
maryam
Join Date: Jan 2012
Posts: 11
Rep Power: 5
sepidehkavousi is on a distinguished road
Quote:
Originally Posted by sara1987 View Post
is it steady state?
have you done grid size independence check?
have you initialized it correctly?
can u tell me more details about the problem u wanna solve

Dear sara...yes it is steady state ...at first I make my model from a paper and I used the same grid size and channel size and I give the initial condition near the real results but it does not get convergence....it is about 3d flow in a micro channel with two fixed boundary condition in inlet and outlet and no slip condition...water is flowing in the channel, I did every thing exactly written in Versteeg for fized boundary condition, but I can not find the problem
sepidehkavousi is offline   Reply With Quote

Old   May 18, 2012, 08:05
Default
  #6
New Member
 
maryam
Join Date: Jan 2012
Posts: 11
Rep Power: 5
sepidehkavousi is on a distinguished road
Quote:
Originally Posted by mihaipruna View Post
can you try setting gradient limiters?


Hi...I couldn't understand you...can you tell me more about this?
sepidehkavousi is offline   Reply With Quote

Old   May 18, 2012, 08:09
Default
  #7
New Member
 
maryam
Join Date: Jan 2012
Posts: 11
Rep Power: 5
sepidehkavousi is on a distinguished road
Quote:
Originally Posted by rmh26 View Post
Can you assume the flow is undisturbed near the boundaries. ie just a posiuelle type flow. If so then as a reasonable approximation you can set the tangential part of the flow to zero and set the derivative of the normal part of the flow to zero. There are some better approximation you can use if the flow isn't so regular.

"Open and Traction boundary conditions for the incompressible navier stokes equations" Jie Liu JCP 228, 2009

thank you my dear friend...I would study the paper and ask you If i had problem...
sepidehkavousi is offline   Reply With Quote

Old   May 18, 2012, 09:03
Default
  #8
Senior Member
 
Mihai Pruna
Join Date: Apr 2010
Location: Boston
Posts: 190
Blog Entries: 1
Rep Power: 7
mihaipruna is on a distinguished road
Quote:
Originally Posted by sepidehkavousi View Post
Hi...I couldn't understand you...can you tell me more about this?
well, when I had this problem in OpenFOAM, I used one of the following settings. I think it gives alink to the code, maybe you can implement it in your program. good luck!

http://openfoamwiki.net/index.php/Op...guide/Limiters
__________________
Mihai Pruna's Bio
mihaipruna is offline   Reply With Quote

Old   May 19, 2012, 03:58
Default
  #9
New Member
 
maryam
Join Date: Jan 2012
Posts: 11
Rep Power: 5
sepidehkavousi is on a distinguished road
Quote:
Originally Posted by mihaipruna View Post
well, when I had this problem in OpenFOAM, I used one of the following settings. I think it gives alink to the code, maybe you can implement it in your program. good luck!

http://openfoamwiki.net/index.php/Op...guide/Limiters
Hi
I am writing my code with fortran language...didn't you understand exactly what OpenFOAM did in that case?
sepidehkavousi is offline   Reply With Quote

Old   May 20, 2012, 12:11
Default
  #10
Senior Member
 
Mihai Pruna
Join Date: Apr 2010
Location: Boston
Posts: 190
Blog Entries: 1
Rep Power: 7
mihaipruna is on a distinguished road
I haven't reviewed the source code for OpenFOAM, but for your program, you should be able to calculate the difference of inflow and outflow in a cell, for any quantity/ direction, divide by cell size in that direction and if the result is larger than a certain number, repeat the iteration with a smaller time step, add relaxation factors maybe, or simply set that gradient to the maximum allowable value and proceed from there, I think this is what OpenFOAM does and although it introduces some accuracies (not much in my case, few cells have that problem), it allows the solution to converge.
__________________
Mihai Pruna's Bio
mihaipruna is offline   Reply With Quote

Reply

Tags
finite volume, inlet pressure, simple

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
Pressure boundary condition C-H Kuo Main CFD Forum 10 August 26, 2015 07:45
Wind turbine simulation Saturn CFX 34 October 16, 2014 05:27
External,incompressible flow boundary condition? John FLUENT 2 August 19, 2011 00:42
transient simulation of a rotating rectangle icesniffer CFX 1 August 8, 2009 07:25


All times are GMT -4. The time now is 11:29.