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

Flow in cylinder

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 18, 2011, 03:03
Question Flow in cylinder
  #1
New Member
 
SNH
Join Date: Dec 2011
Posts: 6
Rep Power: 4
nhosseini is on a distinguished road
Hi every body!

I want to solve a transient, incompressible inviscid flow in a cylinder, which its initial condition is zero velocity everywhere and at delta t=1 left boundary condition is V boundary= Vb and for right boundary condition we have pressure= ambient pressure
it seems that all nodes velocities should reach Vb after 1 time step but I have problem while using conservation methods, it takes some more steps for nodes to reach Vb.
I would be so thankfull if somebody could help me.
nhosseini is offline   Reply With Quote

Old   December 18, 2011, 10:39
Default
  #2
Member
 
Ford Prefect's Avatar
 
Join Date: Mar 2009
Posts: 69
Rep Power: 7
Ford Prefect is on a distinguished road
Quote:
Originally Posted by nhosseini View Post
Hi every body!

I want to solve a transient, incompressible inviscid flow in a cylinder, which its initial condition is zero velocity everywhere and at delta t=1 left boundary condition is V boundary= Vb and for right boundary condition we have pressure= ambient pressure
it seems that all nodes velocities should reach Vb after 1 time step but I have problem while using conservation methods, it takes some more steps for nodes to reach Vb.
I would be so thankfull if somebody could help me.
What boundary condition do you have for pressure at the velocity inlet?
__________________
"Trying is the first step to failure." - Homer Simpson
Ford Prefect is offline   Reply With Quote

Old   December 18, 2011, 11:22
Default
  #3
New Member
 
SNH
Join Date: Dec 2011
Posts: 6
Rep Power: 4
nhosseini is on a distinguished road
no boundary condition, I solve mass equation in inlet for pressure
nhosseini is offline   Reply With Quote

Old   December 18, 2011, 15:23
Default
  #4
Member
 
Ford Prefect's Avatar
 
Join Date: Mar 2009
Posts: 69
Rep Power: 7
Ford Prefect is on a distinguished road
What method do you use to couple velocity and pressure? Do you use some limited compressibility method or a fully incompressible method where div U = 0 ?

Is the behavior independent of time-step and Reynolds number (i.e. if you change dt or Re, will you reach Vb in the entire domain in fewer time-steps)?
__________________
"Trying is the first step to failure." - Homer Simpson
Ford Prefect is offline   Reply With Quote

Old   December 18, 2011, 22:21
Default
  #5
New Member
 
SNH
Join Date: Dec 2011
Posts: 6
Rep Power: 4
nhosseini is on a distinguished road
I use a pressure based method, and yes its behavior is depends on delta t, by increasing time step it reaches the Vb in fewer steps. and by omitting transient term whole field velocity reaches to Vb in one step which is desired, thank you
nhosseini is offline   Reply With Quote

Old   December 19, 2011, 03:45
Default
  #6
Member
 
Ford Prefect's Avatar
 
Join Date: Mar 2009
Posts: 69
Rep Power: 7
Ford Prefect is on a distinguished road
Ok. I assume that you are solving some kind of Poisson equation for pressure and then you need a boundary condition for pressure at the inlet.

You state that you use continuity equation for this. Could you please show what you mean here?

Anyways, some methods specify a homogeneous Neumann condition for pressure at boundaries. However it might be better to use the normal momentum equation instead. By doing so you will compensate for grid size and time-step and the solution will be independent of the two. It will also converge within one time-step, given that your linear solver converges sufficiently.
__________________
"Trying is the first step to failure." - Homer Simpson
Ford Prefect is offline   Reply With Quote

Old   December 19, 2011, 04:51
Default
  #7
New Member
 
SNH
Join Date: Dec 2011
Posts: 6
Rep Power: 4
nhosseini is on a distinguished road
first of all thanks for your quick reply, as I am really confused with it.

I solve du/dx=0 for mass and I use finite element control volume method, so I use u* and u* is pressure based and is

Mass: U*e - Uinlet = 0

u*e =( Ui+Ui+1)/2 + (Pi-Pi+1)/(2*ro*u*e)
U inlet = Vb

and
Momentum: U inlet = Vb

you mean that I need one more condition for pressure?




nhosseini is offline   Reply With Quote

Reply

Tags
cylinder

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
benchmark: flow over a circular cylinder goodegg Main CFD Forum 12 January 22, 2013 11:47
1 Was FSI for a flow over a cylinder vmlxb6 CFX 17 May 16, 2011 02:29
Moving a Cylinder in cross flow after solving Dynamics equations of motion maruthamuthu_venkatraman OpenFOAM 1 November 19, 2009 13:55
Flow over a flat plate & Flow over a cylinder cfdxue Main CFD Forum 0 November 26, 2007 23:26
Flow over a cylinder Anna Main CFD Forum 9 March 24, 2006 14:32


All times are GMT -4. The time now is 07:08.