CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Zigzagging in the pressure

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 16, 2010, 19:39
Default Zigzagging in the pressure
  #1
New Member
 
Francisco Cortés
Join Date: Sep 2010
Location: Santiago, Chile
Posts: 7
Rep Power: 6
fcortes is on a distinguished road
Hello everyone,
I'm using groovyBC to make a pulsatile velocity with a parabolic profile at the inlet of a tube. The solver I'm using is pisoFoam; turbulence model is a RAS kOmega; fluid is Newtonian.

The problem is that when I plot the pressure over time (anywhere, but especially near the inlet), it oscillates heavily from one time step to the next (zigzagging). In spite of this, the pressure does follow a certain tendency related to the velocity. The pressure at the outlet is set on zero and it's zeroGradient at the other patches.

Inlet //a circular patch in the XY plane
{
type groovyBC;
value (0 0 0);
valueExpression "(time()<0.2) ? fourierSeries1*para : fourierSeries2*para";
variables "yp=pts().y;minY=min(yp);maxY=max(yp);para=(-(maxY-pos().y)*(pos().y-minY)-(maxY-pos().x)*(pos().x-minY))/(0.5*pow(maxY-minY,2))*normal();";
}

The function fourierSeries1 stars on cero and is used as a transition to connect to fourierSeries2, which is a periodic pulse.

This phenomenon also occurs when using timeVaryingFixedValue (for a time series equivalent to the one used with groovyBC but a flat profile) but it does not occur when using a fixedValue (for an average value of the time series) nor with its equivalent on timeVaryingFixedValue. In other words, the phenomenon occurs when the velocity is time dependent.

If someone knows what's wrong (or needs more info to do so) please let me know.

Francisco
fcortes is offline   Reply With Quote

Old   November 19, 2010, 15:46
Default
  #2
New Member
 
Francisco Cortés
Join Date: Sep 2010
Location: Santiago, Chile
Posts: 7
Rep Power: 6
fcortes is on a distinguished road
Finally I solved it. I was using CrankNicholson 1 as time scheme and I changed it for Euler.
fcortes is offline   Reply With Quote

Old   January 25, 2012, 11:28
Default
  #3
Member
 
fisch
Join Date: Feb 2010
Posts: 97
Rep Power: 7
fisch is on a distinguished road
Hello,

is it not possible to calculate with groovyBC and 2nd order time integration?
fisch is offline   Reply With Quote

Reply

Tags
parabolic inlet, time dependent velocity, zigzagging

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
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 06:27
Pressure instabilities with interDyMFoam for the floatingObject case nbadano OpenFOAM Running, Solving & CFD 14 March 8, 2011 09:19
Using pressure outlet BC Shreyansh FLUENT 0 November 5, 2010 15:57
Help ASAP! pressure inlet & outlet engahmed FLUENT 0 June 13, 2010 15:33
Setting pressure and velocity in inlet Asghari FLUENT 5 September 22, 2006 13:23


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