CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > CFX

Problem with a periodic boundary condition

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 7, 2012, 13:26
Default Problem with a periodic boundary condition
  #1
New Member
 
martial
Join Date: May 2012
Posts: 7
Rep Power: 13
chuck209 is on a distinguished road
Hi,

I have a problem with a periodic boundary condition and i really don't find how to solve it.

In fact I'm trying to make an iteration on my BC at each time step.For instance, my boundary condition could be:

Pstatic(t+1)=Pstatic(t)+DeltaP ,
with deltaP=constant value, and t=time.

To recuperate the Pstatique at the boundary, i use a CEL function:
Pstatique(t)=areaAve(pressure)@outlet.

The problem is that the Pstatic(t+1) isn't equal to Pstatic(t)+DeltaP, and it is always less.

That's why a try to change my CEL fundtion to:

Pstatic(t+1)=Pstatic(t)

And surprise: the average pressure is decreasing at each time step....(whereas is suppose to be constant)

I think the probleme come from the difference between the Hybrid and conservative value that are always different, but i don't find how to solve it....

is it possible to recuperate the hybrid value in CFX PRE using areaAve(pressure)@outlet?or to directly impose a conservative value?
Does somebody know how to solve this problem?

Thank a lot!!!!!
chuck209 is offline   Reply With Quote

Old   May 7, 2012, 15:18
Default
  #2
Member
 
Ahmet
Join Date: Mar 2011
Location: Ankara/ Turkey
Posts: 37
Rep Power: 15
a_cucen is on a distinguished road
Dear chuck,

Maybe you can use function to do this. You can import a chart that contains pressure changes according to time change. and you can take this function to your boundary condition.
a_cucen is offline   Reply With Quote

Old   May 7, 2012, 16:47
Default
  #3
New Member
 
martial
Join Date: May 2012
Posts: 7
Rep Power: 13
chuck209 is on a distinguished road
Hi!
Thanks for your answer.
Unfortunately, my CEL is more complicated than the one I give you.I really need to recuperate the previously value (Pressure) to find the next value.

Do you think using a routine function could help me?For this, do I need to use Fortran?
chuck209 is offline   Reply With Quote

Old   May 7, 2012, 17:00
Default
  #4
Senior Member
 
Edmund Singer P.E.
Join Date: Aug 2010
Location: Minneapolis, MN
Posts: 511
Rep Power: 20
singer1812 is on a distinguished road
Cant do that within CEL normally but depending on what you want to do, you might be able to delay the update of a variable and get what you want.

Do a search on TRANS_LOOP
singer1812 is offline   Reply With Quote

Old   May 8, 2012, 16:32
Default
  #5
New Member
 
martial
Join Date: May 2012
Posts: 7
Rep Power: 13
chuck209 is on a distinguished road
thank you. I'm trying to do it with de AV.
But for instance, if I would like to simulate a derivative term like dP/dt, may have other choice than doing:
(P(t+1)-P(t))/deltat
I have to discretize, is'nt it?
chuck209 is offline   Reply With Quote

Old   May 8, 2012, 22:38
Default
  #6
Senior Member
 
Edmund Singer P.E.
Join Date: Aug 2010
Location: Minneapolis, MN
Posts: 511
Rep Power: 20
singer1812 is on a distinguished road
are you trying to do this in POST or while the simulation is running?
singer1812 is offline   Reply With Quote

Old   May 9, 2012, 15:12
Default
  #7
New Member
 
martial
Join Date: May 2012
Posts: 7
Rep Power: 13
chuck209 is on a distinguished road
I'm trying to do it when the simulation it's running.
In fact i need to change the BC (static pressure) at the outlet and for that, i have a function:

dP/dt=f(P,axialvelocity)

And i need to use previously value for P because i'm using the following discretisation:

(P(n+1)-P(n))/deltat=f(p(n),axial velocity(n))
Do you know what i mean?
chuck209 is offline   Reply With Quote

Old   May 9, 2012, 15:21
Default
  #8
Senior Member
 
Edmund Singer P.E.
Join Date: Aug 2010
Location: Minneapolis, MN
Posts: 511
Rep Power: 20
singer1812 is on a distinguished road
Yes discretize it. Look up trans_loop for details on how to delay update of additional variable.
singer1812 is offline   Reply With Quote

Old   May 9, 2012, 16:06
Default
  #9
New Member
 
martial
Join Date: May 2012
Posts: 7
Rep Power: 13
chuck209 is on a distinguished road
Ok!
I just have a last question: how can i be sure than the delay has been correctly apply?
I mean is it possible to calculate for each step the variable with the delay and the variable once the formula compute?
chuck209 is offline   Reply With Quote

Old   May 9, 2012, 16:09
Default
  #10
Senior Member
 
Edmund Singer P.E.
Join Date: Aug 2010
Location: Minneapolis, MN
Posts: 511
Rep Power: 20
singer1812 is on a distinguished road
You should most certainly do this. Set up a monitor point for the pressure, and another monitor point for the Additional Variable that is delayed.

You should be able to see the lag in the solver manager.
singer1812 is offline   Reply With Quote

Old   May 9, 2012, 17:15
Default
  #11
New Member
 
martial
Join Date: May 2012
Posts: 7
Rep Power: 13
chuck209 is on a distinguished road
Thanks a lot!
chuck209 is offline   Reply With Quote

Reply

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
Periodic boundary condition Arif FLUENT 3 March 9, 2017 02:18
Problem with assigned inlet velocity profile as a boundary condition Ozgur_ FLUENT 5 August 25, 2015 05:58
problem about periodic boundary condition in Fluent winnawinna FLUENT 0 December 29, 2010 00:32
[Gmsh] Import problem ARC OpenFOAM Meshing & Mesh Conversion 0 February 27, 2010 11:56
Transient outlet boundary condition problem jwillie2000 CFX 1 December 7, 2009 18:07


All times are GMT -4. The time now is 22:36.