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

pressure boundary conditions for collocated grid for Navier-Stokes

Register Blogs Community New Posts Updated Threads Search

Like Tree9Likes

 
 
LinkBack Thread Tools Search this Thread Display Modes
Prev Previous Post   Next Post Next
Old   May 9, 2017, 15:38
Default pressure boundary conditions for collocated grid for Navier-Stokes
  #1
Senior Member
 
Selig
Join Date: Jul 2016
Posts: 213
Rep Power: 10
selig5576 is on a distinguished road
Assumptions:
1. Uniform equally spaced grid
2. 2D

When using the Chorin Projection method to solve the incompressible NS on a collocated grid we end up having to solve the following elliptic PDE

\frac{\delta^{2} P}{(\delta x)^{2}} + \frac{\delta^{2} P}{(\delta y)^{2}} = 
\frac{1}{dt} \left(\frac{uf_{ij}^{*} - uf_{i-1,j}^{*}}{dx} + \frac{vf_{i,j}^{*} - vf_{i,j-1}^{*}}{dx} \right)

where we have the face velocities

uf_{i,j}^{*} = \frac{u^{*}_{i+1,j} + u^{*}_{i,j}}{2},
vf_{i,j}^{*} = \frac{v^{*}_{i,j+1} + v^{*}_{i,j}}{2}

Our boundary conditions on all four sides are
\frac{\partial P}{\partial n} = 0

When solving the pressure Poisson equation on the 4 boundaries is the procedure the same as if we were using the staggered grid? I've read some literature on it and I still don't quite understand if the pressure boundary conditions are treated differently. In my code I have

Code:
P(1,:) = P(2,:)
P(nx,:) = P(nx-1,:)
P(:,1) = P(:,2)
P(:,ny) = P(:,ny-1)
%Corners
P(1,1) = 1.0 % degree of freedom
P(nx,ny) = P(nx-1,ny-1) top right
P(nx,1) = P(nx-1,2) bottom right
P(1,ny) = P(2,ny-1) top left
Which gives me good results for the lid driven cavity, however according to http://www3.nd.edu/~gtryggva/CFD-Cou...re-15-2017.pdf
the pressure boundaries are treated differently. I was curious if my understanding of how to implement the pressure BCs on a collocated is correct. In my code I am using standard central differencing with extra treatment for the handling the ghost points.
selig5576 is offline   Reply With Quote

 


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
Centrifugal fan j0hnny CFX 13 October 1, 2019 13:55
Radiation in semi-transparent media with surface-to-surface model? mpeppels CFX 11 August 22, 2019 07:30
My radial inflow turbine Abo Anas CFX 27 May 11, 2018 01:44
Wrong multiphase flow at rotating interface Sanyo CFX 14 February 7, 2017 17:19
Velocity vector in impeller passage ngoc_tran_bao CFX 24 May 3, 2016 21:16


All times are GMT -4. The time now is 12:13.