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

No-Slip Boundary Conditions for Pressure Equation in SIMPLE Algorithm

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

Like Tree1Likes
  • 1 Post By FMDenaro

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 4, 2020, 12:55
Default No-Slip Boundary Conditions for Pressure Equation in SIMPLE Algorithm
  #1
New Member
 
Rafael March
Join Date: Mar 2020
Posts: 5
Rep Power: 6
rafaelmarch3 is on a distinguished road
Hello all,

I have implemented the SIMPLE algorithm to solve steady-state Stokes Equations in 3D structured grid.

My Geometry is a cube. My boundary conditions are: 1) Fixed pressure at right and left boundary faces (call it x- and x+) and no-slip elsewhere.

I have a question regarding the consistent formulation of no-slip boundary conditions.

For the momentum equations, it is trivial to implement no-slip conditions: simply set ux and uy=0 at the y+, y-, z+ and z- boundary faces.

However, what's the correct boundary condition for pressure in these boundary faces? Right now I'm simply assuming dP/dn = 0, which means a zero pressure gradient in at these boundary sides. Is this the correct assumption? I have seen some references e.g. [1] suggesting to project the momentum equation in the normal direction to get a boundary condition for pressure. What do you think?

REFERENCES:
[1] ON PRESSURE BOUNDARY CONDITIONS FOR THE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS, GRESHO AND SANI, INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, VOL. 7, 1 1 1 1 - 1 145 (1987)
rafaelmarch3 is offline   Reply With Quote

Old   November 5, 2020, 05:44
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,760
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by rafaelmarch3 View Post
Hello all,

I have implemented the SIMPLE algorithm to solve steady-state Stokes Equations in 3D structured grid.

My Geometry is a cube. My boundary conditions are: 1) Fixed pressure at right and left boundary faces (call it x- and x+) and no-slip elsewhere.

I have a question regarding the consistent formulation of no-slip boundary conditions.

For the momentum equations, it is trivial to implement no-slip conditions: simply set ux and uy=0 at the y+, y-, z+ and z- boundary faces.

However, what's the correct boundary condition for pressure in these boundary faces? Right now I'm simply assuming dP/dn = 0, which means a zero pressure gradient in at these boundary sides. Is this the correct assumption? I have seen some references e.g. [1] suggesting to project the momentum equation in the normal direction to get a boundary condition for pressure. What do you think?

REFERENCES:
[1] ON PRESSURE BOUNDARY CONDITIONS FOR THE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS, GRESHO AND SANI, INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, VOL. 7, 1 1 1 1 - 1 145 (1987)



The BCs for the pressure equation have nothing to do with the no-slip condition. They are derived from the mass conservation, for this reason the Neuman condition for the pressure is prescribed. The condition dp/dn=0 can be applied provided that the source term is congruently modified, this way the compatibility relation will be fulfilled and a solution exists (apart a constant).
rafaelmarch3 likes this.
FMDenaro is offline   Reply With Quote

Old   November 5, 2020, 06:24
Default
  #3
New Member
 
Rafael March
Join Date: Mar 2020
Posts: 5
Rep Power: 6
rafaelmarch3 is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
The BCs for the pressure equation have nothing to do with the no-slip condition. They are derived from the mass conservation, for this reason the Neuman condition for the pressure is prescribed. The condition dp/dn=0 can be applied provided that the source term is congruently modified, this way the compatibility relation will be fulfilled and a solution exists (apart a constant).
Hi Filippo, many thanks for your reply.

Yes, when I mentioned the no-slip condition, I wanted to make sure the compatibility relation is being fulfilled. Right now I don't have a source term in the pressure Equation. So I assume my set of BCs are not consistent. I notice that my code converges or diverges depending on the initial condition, which made me suspicious on my choice of BCs for the momentum and pressure equations.

I would implement the consistent conditions for pressure like in the equations below. Note that I have a Stokes-Brinkman system, meaning that I have a term proportional to the velocity that vanishes at the boundary (no-slip). I would apply these boundary conditions for pressure using the velocity values in the previous iteration. I guess my concern is that this condition does not have any information on the wall shear (tau_x = dux/dy). So I'm not entirely sure this is the proper way.

https://imgur.com/ktaJoOY

Could you please point me to a reference on how to find the proper compatible boundary conditions for the pressure equation to have a well-posed system?

Thank you,
Rafael.
rafaelmarch3 is offline   Reply With Quote

Old   November 5, 2020, 07:03
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,760
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by rafaelmarch3 View Post
Hi Filippo, many thanks for your reply.

Yes, when I mentioned the no-slip condition, I wanted to make sure the compatibility relation is being fulfilled. Right now I don't have a source term in the pressure Equation. So I assume my set of BCs are not consistent. I notice that my code converges or diverges depending on the initial condition, which made me suspicious on my choice of BCs for the momentum and pressure equations.

I would implement the consistent conditions for pressure like in the equations below. Note that I have a Stokes-Brinkman system, meaning that I have a term proportional to the velocity that vanishes at the boundary (no-slip). I would apply these boundary conditions for pressure using the velocity values in the previous iteration. I guess my concern is that this condition does not have any information on the wall shear (tau_x = dux/dy). So I'm not entirely sure this is the proper way.

https://imgur.com/ktaJoOY

Could you please point me to a reference on how to find the proper compatible boundary conditions for the pressure equation to have a well-posed system?

Thank you,
Rafael.

Here You can find many details about the derivativon of the pressure equation and the compatibility relation
https://www.researchgate.net/publica...ary_conditions
FMDenaro is offline   Reply With Quote

Reply

Tags
boundary conditions, pressure poisson, simple algorithm

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Question about adaptive timestepping Guille1811 CFX 25 November 12, 2017 18:38
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 07:27
An error has occurred in cfx5solve: volo87 CFX 5 June 14, 2013 18:44
RPM in Wind Turbine Pankaj CFX 9 November 23, 2009 05:05
Convective Heat Transfer - Heat Exchanger Mark CFX 6 November 15, 2004 16:55


All times are GMT -4. The time now is 15:10.