# Flow in cylinder

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

 December 18, 2011, 04:03 Flow in cylinder #1 New Member   SNH Join Date: Dec 2011 Posts: 8 Rep Power: 6 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.

December 18, 2011, 11:39
#2
Member

Join Date: Mar 2009
Posts: 72
Rep Power: 9
Quote:
 Originally Posted by nhosseini 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

 December 18, 2011, 12:22 #3 New Member   SNH Join Date: Dec 2011 Posts: 8 Rep Power: 6 no boundary condition, I solve mass equation in inlet for pressure

 December 18, 2011, 16:23 #4 Member     Join Date: Mar 2009 Posts: 72 Rep Power: 9 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

 December 18, 2011, 23:21 #5 New Member   SNH Join Date: Dec 2011 Posts: 8 Rep Power: 6 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

 December 19, 2011, 04:45 #6 Member     Join Date: Mar 2009 Posts: 72 Rep Power: 9 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

 December 19, 2011, 05:51 #7 New Member   SNH Join Date: Dec 2011 Posts: 8 Rep Power: 6 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?

 Tags cylinder

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post goodegg Main CFD Forum 12 January 22, 2013 12:47 vmlxb6 CFX 17 May 16, 2011 02:29 maruthamuthu_venkatraman OpenFOAM 1 November 19, 2009 14:55 cfdxue Main CFD Forum 0 November 27, 2007 00:26 Anna Main CFD Forum 9 March 24, 2006 15:32

All times are GMT -4. The time now is 05:51.