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

velocity increases extremely

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 30, 2010, 08:55
Default velocity increases extremely
  #1
Member
 
Axel Söhngen
Join Date: Jan 2010
Location: Germany, Trier
Posts: 31
Rep Power: 16
chelvistero is on a distinguished road
Hello,
I'm trying to create an external flow over a cylinder. I made 2D mesh with blockMesh. First I try the case with a velocity of 1 m/s and Re=150. I'm using icoFoam but when I start the calculation the velocity increses with every timestep but it has to be constant.
Whats the mistake?? Here are p and U dictionaries:

p:

dimensions [0 2 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
frontAndBack
{
type empty;
}

inlet
{
type zeroGradient;
}

outlet
{
type fixedValue;
value uniform 0;
}

lowerWall
{
type symmetryPlane;
}

upperWall
{
type symmetryPlane;
}

zylinder
{
type zeroGradient;
}
}


U:

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
inlet
{
type fixedValue;
value uniform (1 0 0);
}

outlet
{
type zeroGradient;
}

frontAndBack
{
type empty;
}

lowerWall
{
type symmetryPlane;
}

upperWall
{
type symmetryPlane;
}

zylinder
{
type fixedValue;
value uniform (0 0 0);
}



Thanks in advance
chelvistero is offline   Reply With Quote

Old   January 30, 2010, 11:32
Default
  #2
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 20
sega is on a distinguished road
Are you doing a transient or steady-state simulation?
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   January 31, 2010, 07:36
Default
  #3
Member
 
Axel Söhngen
Join Date: Jan 2010
Location: Germany, Trier
Posts: 31
Rep Power: 16
chelvistero is on a distinguished road
Hello sega,
I'm doing a transient simulation!
chelvistero is offline   Reply With Quote

Old   January 31, 2010, 07:54
Default
  #4
Senior Member
 
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 18
idrama is on a distinguished road
How does it themselves felt? I mean, do ya have, for instance, a CFL-drop or do you see it by doing pros-processing?

Cheers
idrama is offline   Reply With Quote

Old   January 31, 2010, 08:20
Default
  #5
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 20
sega is on a distinguished road
I'm not sure, but have you tried a different boundary condition for pressure on the outlet? Prescribing a pressure may be not always the right choice.
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   January 31, 2010, 08:27
Default
  #6
Member
 
Axel Söhngen
Join Date: Jan 2010
Location: Germany, Trier
Posts: 31
Rep Power: 16
chelvistero is on a distinguished road
Hello Claus,
the cfl increases extremely, and in paraview is shown that the velocity increases every timestep too. I try to simulate the karman vortex street with a Re of 100.
chelvistero is offline   Reply With Quote

Old   January 31, 2010, 08:53
Default
  #7
Member
 
Axel Söhngen
Join Date: Jan 2010
Location: Germany, Trier
Posts: 31
Rep Power: 16
chelvistero is on a distinguished road
No I haven't tried another boundary for pressure at outlet! What else can I choose for the pressure at outlet? And is the pressure at the cylinder (uniform 0) the right choice?
Thank you for your reply
chelvistero is offline   Reply With Quote

Old   January 31, 2010, 09:12
Default
  #8
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 20
sega is on a distinguished road
Quote:
Originally Posted by chelvistero View Post
And is the pressure at the cylinder (uniform 0) the right choice?
I don't think so. Think of the Bernoulli equation.
I would choose zeroGradient!
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   January 31, 2010, 09:25
Default
  #9
Member
 
Axel Söhngen
Join Date: Jan 2010
Location: Germany, Trier
Posts: 31
Rep Power: 16
chelvistero is on a distinguished road
Ok that's right. But the courant number and velocity still explode.
chelvistero is offline   Reply With Quote

Old   January 31, 2010, 17:21
Default
  #10
Senior Member
 
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 18
idrama is on a distinguished road
Choosing a zeroGradient is the worse thing you can do, because of the pressure term in the Navier-Stokes-Equation, i.e. you have to specified the pressure at least in one point.

Have you checked your mesh? A blow up is also related to singularities which can occur by bifurcation, so try a finer mesh or check if it fit the physical behavior of your scenario.

Cheers
idrama is offline   Reply With Quote

Reply


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
ATTENTION! Reliability problems in CFX 5.7 Joseph CFX 14 April 20, 2010 16:45
Velocity in Porous medium : HELP! HELP! HELP! Kali Sanjay Phoenics 0 November 6, 2006 07:10
Variables Definition in CFX Solver 5.6 R P CFX 2 October 26, 2004 03:13
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 10:11
what the result is negatif pressure at inlet chong chee nan FLUENT 0 December 29, 2001 06:13


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