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

karman vortex street with turbulence and heated walls

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 25, 2011, 07:59
Default karman vortex street with turbulence and heated walls
  #1
New Member
 
Jan Goebel
Join Date: Mar 2009
Location: Mannheim, Baden Wuerttemberg, Germany
Posts: 24
Rep Power: 7
jango is on a distinguished road
Hello Everybody,

To get used to OpenFOAM cases and solvers I`m playing around a bit with the karman vortex street. With the icoFoam Solver for laminar flow at low Reynolds nubmer this isn`t a problem.

But now I want to increase the Reynolds number und add a k-epsilon (near) wall treatment. The fluid is air with a temperature of 273 K and the wall temperature of the cylinder is 573 K. For this case I assumed, that rhoPimpleFoam should be the best solver. With the standard adjustments for this solver, the case runs without any convergence problems but I can see no vortexes. So I changed the runtime for the turbulent case an adjusted the material properties (thermophysical properties) to constantTransport with the according values of Pr and mu for air.

After a few iterations the solution diverges. Now I have a few questions:

Is the rhoPimpleFoam solver the best choice for this case ?
Adjusting the thermophysicalProperties according to this Problem (considering adjustments for Temperature and material properties for air) I get errormessages and convergence problems. Should I adjust the underrelaxation factors to stabilize the solution ?

Thank you a lot for your Help

Yours,

Jan
jango is offline   Reply With Quote

Old   July 25, 2011, 10:08
Default
  #2
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 773
Rep Power: 11
Bernhard is on a distinguished road
Do you understand turbulence modelling? What do you mean by 'k-epsilon wall treatment'? If you're interesting in vortex shedding then you can better use a large eddy simulations, then you're solving for at least a part of the turbulent spectrum. And why do you include temperature? This is not needed for the vortex street to develop. Start with isothermal simulations. If you're sure those are correct, then you can start with solving for non-constant densities. Good luck!
Bernhard is offline   Reply With Quote

Old   July 25, 2011, 10:50
Default
  #3
New Member
 
Jan Goebel
Join Date: Mar 2009
Location: Mannheim, Baden Wuerttemberg, Germany
Posts: 24
Rep Power: 7
jango is on a distinguished road
Hello Bernhard,

thank you for your quick response. Of course the temperature isn`t important for the vortex street. As I learnded CFD with Fluent at the university, I have various test cases based on a fluent script and I'm just trying various settings and solvers for the better understanding of OpenFOAM.

Thank you for the hint with the large Eddy- Simulation, but is rhoPimpleFoam the best choice for the turbulence modelling (and maybe heat transfer) or is simpleFOAM the better choice ?
jango is offline   Reply With Quote

Old   July 25, 2011, 11:43
Default
  #4
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 773
Rep Power: 11
Bernhard is on a distinguished road
simpleFoam is designed for steady-state problems, and your case is transient by definition, so simpleFoam is not an option.
Bernhard is offline   Reply With Quote

Old   July 28, 2011, 10:38
Default
  #5
New Member
 
Jan Goebel
Join Date: Mar 2009
Location: Mannheim, Baden Wuerttemberg, Germany
Posts: 24
Rep Power: 7
jango is on a distinguished road
Thank you for the hints, that helped me a lot. So I just successfully calculated the karman vortex street with the icoFoam solver, increased the reynolds number for the turbulent case and calculated it with the psioFoam solver. Then I changed to compressible medium (air) with the rhoPimpleFoam solver. Then I mapped the fields from this calculation adding Temperature and running the job again. But evaluating the results, I can see, that temperature is not considered in my calculation. Could it be, that my initial definitions in the files 0/T and 0/alphat are erroneous ?

My setting for T:

dimensions [0 0 0 1 0 0 0];

internalField uniform 293;
internalField2 uniform 593;

boundaryField
{
far_field
{
type symmetryPlane;
}
outlet
{
type calculated;
// inletValue uniform 293;
// value uniform 293;
}
inlet
{
type fixedValue;
value uniform 293;
}
cyl_wall
{
type fixedValue;
value $internalField2
}
frontAndBackPlanes
{
type empty;
}
}

and alphat:

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

internalField uniform 0;

boundaryField
{
far_field
{
type symmetryPlane;
}
outlet
{
type calculated;
value uniform 0;
}
inlet
{
type calculated;
value uniform 0;
}
cyl_wall
{
type alphatWallFunction;
value uniform 0;
}
frontAndBackPlanes
{
type empty;
}
}

What can I do to calculate the temperature ?
jango is offline   Reply With Quote

Old   July 28, 2011, 11:14
Default
  #6
New Member
 
charbel's Avatar
 
Charbel
Join Date: Jan 2011
Location: France
Posts: 29
Rep Power: 5
charbel is on a distinguished road
Quote:
Originally Posted by jango View Post
What can I do to calculate the temperature ?
See this link to learn how to add temperature to icoFoam solver.
__________________
--
Charbel
charbel is offline   Reply With Quote

Old   July 28, 2011, 15:29
Default
  #7
New Member
 
Jan Goebel
Join Date: Mar 2009
Location: Mannheim, Baden Wuerttemberg, Germany
Posts: 24
Rep Power: 7
jango is on a distinguished road
Thank you charbel, I think that could be the solution. But I'm still a bit confused. Taking a look at the rhoPimpleFoam solver for compressible fluids and the angledDuct tutorial case, there`s already temperature added in this tutorial case and the field must already be there, but in my case, the fluid has no temperature. As the fluid should be heated by the temperature applied to the cylinder, could it be that my initial conditions for alphat are wrong ?
jango 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
Heated Walls vitor OpenFOAM Running, Solving & CFD 5 August 16, 2010 14:47
Natural convection - Inlet boundary condition max91 CFX 1 July 29, 2008 20:28
Turbulence boundary values lego CFX 9 October 25, 2002 11:55
Turbulence boundary values lego Main CFD Forum 0 October 24, 2002 13:47


All times are GMT -4. The time now is 16:09.