# Bpundary condition of SIMPLER

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 October 7, 2002, 03:44 Bpundary condition of SIMPLER #1 jack Guest   Posts: n/a Dear ALL: I am using the SIMPLER algorithm to firstly calculate the 2D duct flow. Giving the following BC: 1) Entry of a duct: imposed velocity, pressure unknown. 2) Exit of a duct: velocity unknown (0 neumann boundary condition: pu/pn=0 ), pressure imposed and equal to 0. Both pressure and pressure correction equations were solved. The pressure correction results are only used to correct velocities. The pressure was diectly solved by the pressure equation. The result for Stokes flow (no inertia term) is correct. However, with the Reynolds number increases (low only 100), the solution become worse. The calculation still converged, but result is not correct. The solved pressure at inlet is not constant and distorted. Should I give a pressure BC at inlet? How to deal with pressure at inlet? Thanks in advance. Jack

 October 8, 2002, 08:27 Re: Bpundary condition of SIMPLER #2 Seb Perron Guest   Posts: n/a Some simple rules when dealing with incompressible flows (then always subsonic...): -whenever you impose the normal velocity, the pressure has to be unknown (free), hence if the velocity is imposed at the inlet pressure must not. -whenever the normal velocity is unknown (such as an outlet), the pressure must be imposed. Your problem might not be related to the simpler algorithm, -do you use an upwind scheme to stabilize your scheme? -do you use a staggered grid or re-interpolation to make sure the space of approximation for pressure and velocity are different? -your relaxation parameters might not be correct. regards.

 October 8, 2002, 21:00 Re: Bpundary condition of SIMPLER #3 Halim Choi Guest   Posts: n/a

 October 8, 2002, 21:29 Re: Bpundary condition of SIMPLER #4 Halim Choi Guest   Posts: n/a See the following references for treatment of pressure boundary conditions in the SIMPLE algorithm. (1) S. R. Mathur and J. Y. Murthy, Pressure boundary coditions for incompressible flow using unstructured meshes, Numerical Heat Transfer, Part B, Vol. 32, pp. 283-298, 1997, (2) K. M. Kelkar and D. Choudhury, Numerical method for the prediction of incompressible flow and heat transfer in domains with specified pressure boundary conditions, Numerical Heat Transfer, Part B. vol. 38, pp 15-36, 2000. There are two ways to treat the outlet boundary conditions (1) to specify the pressure at outlet : read above refernces. I think it is not correct to specify the outlet velocities using the Neuman conditions(du/dn=0.) when you use the pressure boundary condition. (2) to correct the outlet velocities to satisfy the overall mass conservation : for example the flow in the channel, during the iteration process, you can calculate alpa=(sum of inlet mass flow rate)/ ( sum of outlet mass flow rate) and then multiply alpa for the oulet velocities such as u(imax,j)=alpa*u(imax,j), and the pressure at the outlet can be obtained from linear extrapolation of that in the adjacent interior nodes. I hope this helps, Halim Choi

 October 9, 2002, 04:33 Re: Bpundary condition of SIMPLER #5 jack Guest   Posts: n/a Thank you very much for your help. If SIMPLE algorithm, the problem should be Ok because the pressure is obtained by pressure correction equation. Pressure can be calculated refert to a reference node where p is set to zero. But, This is SIMPLER, pressure is obtained by pressure equation. Yes, an upwind scheme is used to stabilize scheme. I use co-located mesh and interpolation like Chow and Chie to make sure the space approximation for pressure and velocity coupling. The relaxation parameters use implicit way in the solution of the algebraic equations. The relaxation parameter is 0.5 for velocities, and 1 for pressure. regards. Jack

 October 9, 2002, 04:36 Re: Bpundary condition of SIMPLER #6 jack Guest   Posts: n/a Thanks for Halim Choi's help for this peoblem. I will try to read the mentioned paers. Best regards, Jack

 October 10, 2002, 20:24 Re: Bpundary condition of SIMPLER #7 Ramadan Bassiouny Guest   Posts: n/a I would agree with the points Seb Perron recorded in his message. However, I would expect that you have to check the value of Peclet number (convective/diffusion). You might try the power law scheme. For more information about treating the outflow boundary conditions read the two pages (102-103) and pages (129-130) for pressure boundary conditions in Patanker'a book titled" Numerical heat transfer and fluid flow". Best of luck

 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 fard Main CFD Forum 0 September 26, 2011 10:16 bearcharge Main CFD Forum 0 May 14, 2010 11:32 shyamprasad Main CFD Forum 0 March 17, 2009 00:40 Ambrose Main CFD Forum 3 November 18, 2002 10:27 Ramesh Main CFD Forum 8 November 8, 2002 14:56

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

 Contact Us - CFD Online - Top