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

Periodic Boundary Condition

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

Like Tree2Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   October 1, 2012, 11:50
Default Periodic Boundary Condition
  #1
C.C
Member
 
CC
Join Date: Jun 2011
Posts: 52
Rep Power: 6
C.C is on a distinguished road
Hi all,
I have a problem... I want simulate a flow in a pipe considering 2D axissymetric geometry with periodic boundary condition, but the pressure across the pipe length have a strange behavior with negative values. If I consider the inlet velocity boundary condition and outlet pressure boundary condition, the pressure profile is correct, but I need an entrance region and with the periodic boundary condition I can reduce the length...
Please, someone can help me... What is wrong? Is the periodic boundary condition advisable for my problem?
Thanks
C.C is offline   Reply With Quote

Old   October 2, 2012, 04:15
Default
  #2
Member
 
Hassan Nemati
Join Date: Apr 2012
Posts: 30
Rep Power: 5
hnemati is on a distinguished road
do you use pressure force in streamwise direction?
hnemati is offline   Reply With Quote

Old   October 2, 2012, 05:47
Default
  #3
C.C
Member
 
CC
Join Date: Jun 2011
Posts: 52
Rep Power: 6
C.C is on a distinguished road
I specified the mass flow rate and I want the value of pressure drop to compare with experimental data.
C.C is offline   Reply With Quote

Old   October 2, 2012, 05:48
Default
  #4
Senior Member
 
flotus1's Avatar
 
Alex
Join Date: Jun 2012
Location: Germany
Posts: 1,098
Rep Power: 19
flotus1 will become famous soon enoughflotus1 will become famous soon enough
Applying a periodic boundary condition in streamwise direction by simply copying the values from the outlet to the inlet results in the situation you observe.
Since there is a pressure drop along the channel, the pressure level gets lower every iteration. You could check this by running two simulations with different numbers of iterations.
flotus1 is offline   Reply With Quote

Old   October 2, 2012, 05:59
Default
  #5
Member
 
Hassan Nemati
Join Date: Apr 2012
Posts: 30
Rep Power: 5
hnemati is on a distinguished road
Yes I know, but for priodic simulation the important part is dp which should be included in stramwise momentum equation to force the flow. for example for a channel this pressure force is :dpdx = 12./Re ,where Re is based on height of channel.
hnemati is offline   Reply With Quote

Old   October 2, 2012, 06:08
Default
  #6
Senior Member
 
flotus1's Avatar
 
Alex
Join Date: Jun 2012
Location: Germany
Posts: 1,098
Rep Power: 19
flotus1 will become famous soon enoughflotus1 will become famous soon enough
So then what is the problem?
in an incompressible flow, you can compare the pressure drop (i.e. pressure derivative) to experimental data, independent of the pressure level.

Of course it would be more correct to rescale the pressure at the periodic interface, but the derivative remains unchanged by this procedure.
flotus1 is offline   Reply With Quote

Old   October 3, 2012, 05:17
Default
  #7
Senior Member
 
Join Date: Aug 2011
Posts: 251
Rep Power: 6
leflix is on a distinguished road
Quote:
Originally Posted by C.C View Post
Is the periodic boundary condition advisable for my problem?
Thanks
Yes for sure, but I guess the problem is how you implement these periodic boundary conditions.
first you should apply Phy(NI)=Phy(1) and Phy(0)=Phy(NI-1)
(BC are on node 0 and NI)
Phy stands for u,v, P.
If your solver is explicit it is easy, if your solver is implicit then you have to
deal with this in your matrix.
After you have to include in the x momentum the pressure drop to maintain the flow. It is related to the mass flow rate and in the case of a pipe flow or a channel flow there is an analytical expression.
leflix is offline   Reply With Quote

Old   October 3, 2012, 05:49
Default
  #8
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,599
Rep Power: 22
FMDenaro will become famous soon enough
Quote:
Originally Posted by leflix View Post
Yes for sure, but I guess the problem is how you implement these periodic boundary conditions.
first you should apply Phy(NI)=Phy(1) and Phy(0)=Phy(NI-1)
(BC are on node 0 and NI)
Phy stands for u,v, P.
If your solver is explicit it is easy, if your solver is implicit then you have to
deal with this in your matrix.
After you have to include in the x momentum the pressure drop to maintain the flow. It is related to the mass flow rate and in the case of a pipe flow or a channel flow there is an analytical expression.
I agree, if you use the non-dimensional form of the equation, where the Re_tau is the relevant non-dimensional number, the non-dimensional pressure gradient the force the flow in the stream-wise direction is simply -1
FMDenaro is offline   Reply With Quote

Old   October 3, 2012, 06:30
Default
  #9
C.C
Member
 
CC
Join Date: Jun 2011
Posts: 52
Rep Power: 6
C.C is on a distinguished road
Quote:
Originally Posted by leflix View Post
Yes for sure, but I guess the problem is how you implement these periodic boundary conditions.
first you should apply Phy(NI)=Phy(1) and Phy(0)=Phy(NI-1)
(BC are on node 0 and NI)
Phy stands for u,v, P.
If your solver is explicit it is easy, if your solver is implicit then you have to
deal with this in your matrix.
After you have to include in the x momentum the pressure drop to maintain the flow. It is related to the mass flow rate and in the case of a pipe flow or a channel flow there is an analytical expression.
I'm using Ansys Fluent to simulate the pipe flow and in the pressure profile the pressure has negative values
C.C is offline   Reply With Quote

Old   October 3, 2012, 06:41
Default
  #10
Senior Member
 
Join Date: Aug 2011
Posts: 251
Rep Power: 6
leflix is on a distinguished road
Quote:
Originally Posted by C.C View Post
I'm using Ansys Fluent to simulate the pipe flow and in the pressure profile the pressure has negative values

But if you use Fluent , periodic BC are already implemented and you have nothing to do, just click on "Periodic Boundary Condition"
If it doesn't work it is because you certainly did something wrong previously...
leflix is offline   Reply With Quote

Old   October 3, 2012, 06:42
Default
  #11
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,599
Rep Power: 22
FMDenaro will become famous soon enough
Quote:
Originally Posted by C.C View Post
I'm using Ansys Fluent to simulate the pipe flow and in the pressure profile the pressure has negative values
negative values of the pressure means nothing useful... you must assess that youy have a basic (constant) negative pressure gradient with a superimposed periodic pressure fluctuation
FMDenaro is offline   Reply With Quote

Old   October 7, 2012, 13:51
Default 2D Channel Flow eith Periodic boundaries in the streamwise directions
  #12
Member
 
Prasanth P
Join Date: May 2009
Posts: 30
Rep Power: 8
prasanthnitt is on a distinguished road
Hi All,

I have written a 2D Incompressible NS solver for solving a channel flow with periodic boundary conditions. I am using a fractional step method(Kim and Moin 1985 JCP). I have non dimensionalized the equations. The rey_no is set to 100.

Channel Dimensions:
2 units: in the stream-wise and in the normal direction
Time integration scheme:
AB-2

A constant pressure gradient of value -2.0/rey_no is specified.
Staggered grid.
Periodic bcn in the stream-wise direction:
u(ni,.)=u(1,.)
u(0,.)=u(ni-1,.)

V(ni+1,.)=v(1,.)
v(0,.)=v(ni,.)

p(ni+1,.)=p(1,.)
p(0,.)=p(ni,.)

u and v represent the horizontal velocity components.
p represents pressure.

Boundary conditions in the normal direction:
u(:,0)=-u(:,1) //No Slip
u(:,nj+1)=-u(:,nj) //No penetration

v(:,0)=0.0
v(:,nj)=0.0

dp/dn =0 at y=0 and y=2

Initial condition:
I have tried several:
u(:, j )=1.0
u(:,j)=parabola with centerline velocity=1

So, the analytical solution for this problem is a parabola with centerline velocity=1.
But in my case the centerline velocity is greater than 1. I don't know why this is happening.
Help of any kind would be appreciated. I would be more than happy to provide any extra information.

Cheers
Prasanth
prasanthnitt is offline   Reply With Quote

Old   October 12, 2012, 12:48
Default
  #13
Senior Member
 
Join Date: Aug 2011
Posts: 251
Rep Power: 6
leflix is on a distinguished road
Quote:
Originally Posted by prasanthnitt View Post

So, the analytical solution for this problem is a parabola with centerline velocity=1.
But in my case the centerline velocity is greater than 1. I don't know why this is happening.

Prasanth
Hi Prasanth,

No in the case of a channel flow the dimensionless velocity on the centerline is 1.5
Is it what you obtained?
leflix is offline   Reply With Quote

Old   October 13, 2012, 01:29
Default
  #14
Member
 
Prasanth P
Join Date: May 2009
Posts: 30
Rep Power: 8
prasanthnitt is on a distinguished road
[QUOTE

No in the case of a channel flow the dimensionless velocity on the centerline is 1.5
Is it what you obtained?[/QUOTE]

Its 1.5 if you normalize with the average velocity. I am normalizing with the centre-line velocity.
The issue has been fixed. I applied the pressure gradient in the velocity update step (from u* to u n+1) and it worked.

Thanks.
prasanthnitt is offline   Reply With Quote

Old   October 14, 2012, 11:41
Default
  #15
Senior Member
 
Join Date: Aug 2011
Posts: 251
Rep Power: 6
leflix is on a distinguished road
Quote:
Originally Posted by prasanthnitt View Post
[QUOTE

I am normalizing with the centre-line velocity.
The issue has been fixed. I applied the pressure gradient in the velocity update step (from u* to u n+1) and it worked.

Thanks.
ok it's fine!!
leflix is offline   Reply With Quote

Old   November 17, 2012, 10:28
Default
  #16
Member
 
Join Date: Sep 2011
Posts: 38
Rep Power: 5
Kamu is on a distinguished road
Hello all
I am using a periodic boundary condition pipe flow with an insert. I read in FLUENT user guide that the pressure drop has two components, a linear varying component and a periodic component. How do you get the periodic component? Is there a way of post processing the pressure drop that includes the two components???
Kamu is offline   Reply With Quote

Old   November 18, 2012, 11:21
Default
  #17
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 531
Blog Entries: 14
Rep Power: 17
sbaffini will become famous soon enough
There is no reason to post process the linearly varying part, as ALL the information you need about it is in the periodic condition panel, that is gradient and direction. The periodic component is the only one you get as solution and it is the only one available in the post processing.

If you want, you can define a custom field function with the following variable:

p + dp/dx*(x-x0) + dp/dy*(y-y0) + dp/dz*(z-z0)

and you get them both
sbaffini is offline   Reply With Quote

Old   November 19, 2012, 02:52
Default
  #18
Member
 
Join Date: Sep 2011
Posts: 38
Rep Power: 5
Kamu is on a distinguished road
Quote:
Originally Posted by sbaffini View Post
There is no reason to post process the linearly varying part, as ALL the information you need about it is in the periodic condition panel, that is gradient and direction. The periodic component is the only one you get as solution and it is the only one available in the post processing.

If you want, you can define a custom field function with the following variable:

p + dp/dx*(x-x0) + dp/dy*(y-y0) + dp/dz*(z-z0)

and you get them both
Thanks Sbaffini for your reply. When i put in the custom field function i will get the pressure at a given point. My dilemma is that when you define periodic conditions the inlet and outlet are combined and the pressure I get is at inlet. to get the pressure drop, dp = Pin-Pout! How do I get Pout?
Kamu is offline   Reply With Quote

Old   November 19, 2012, 05:50
Default
  #19
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 531
Blog Entries: 14
Rep Power: 17
sbaffini will become famous soon enough
In the custom field function above, the pressure is the one effectively solved for by Fluent; in this case it gives you p_in = p_out (at least for the straight pipe case).

If you performed a periodic computation by fixing the mass flow, then dp/dx is available in the periodic b.c. panel (it also has a direction, which you needed to fix... hopefully, along the pipe axis). So, to know the presure jump in your case you just need to compute dp/dx * L where L is the length of the pipe
C.C and Kamu like this.
sbaffini is offline   Reply With Quote

Old   November 19, 2012, 05:58
Default Periodic pressure drop
  #20
Member
 
Join Date: Sep 2011
Posts: 38
Rep Power: 5
Kamu is on a distinguished road
Quote:
Originally Posted by sbaffini View Post
In the custom field function above, the pressure is the one effectively solved for by Fluent; in this case it gives you p_in = p_out (at least for the straight pipe case).

If you performed a periodic computation by fixing the mass flow, then dp/dx is available in the periodic b.c. panel (it also has a direction, which you needed to fix... hopefully, along the pipe axis). So, to know the presure jump in your case you just need to compute dp/dx * L where L is the length of the pipe
Dear Paolo,
Thanks very much for the help. This is what i had settled for i would get dp/dx and then multiply by the length of the pipe. I realized that when i compute the friction factor it does not agree with that from experimental data. And when I use the linear gradient in the periodic panel it gives me an accurate answer! Could I be getting dp/dx at a wrong surface?
Kamu is offline   Reply With Quote

Reply

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
periodic boundary condition (translational) fluent 6.3 haihek FLUENT 3 January 22, 2014 05:10
Problem with a periodic boundary condition chuck209 CFX 10 May 9, 2012 16:15
periodic boundary condition in LES martor FLUENT 2 April 10, 2012 19:17
External Radiation Boundary Condition for Grid Interface CFD XUE FLUENT 0 July 9, 2010 02:53
translational periodic boundary condition Rola Afify FLUENT 2 September 12, 2006 08:39


All times are GMT -4. The time now is 21:43.