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

Circular Cylinder at low Reynold number mesh and non orthogonal corrections

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 13, 2008, 05:10
Default Hello everybody, I try to s
  #1
New Member
 
Christophe Kassiotis
Join Date: Mar 2009
Location: Paris
Posts: 17
Rep Power: 17
kassiotis is on a distinguished road
Hello everybody,

I try to solve the benchmark proposed ans studied by Schäfer and Turek in this paper. I have found some linked topic here and here, but they don't really answer my questions.

For the moment, I consider the first case proposed in the benchmark (steady simulation, parabolic input velocity with max 0,3m/s leading to a Reynold number around 20).

You will find here some pictures of typical meshes I use:





And the "annoying" results:



As one can see, the velocity field is not really physical as the flow don't want to enter the fine mesh surrounding the cylinder.

I use icoFoam with steadyState time integration scheme. The initial condition are zero velocity and pressure fields (according to the description of the benchmark). I'm doing one time step (Is the steady solution independant of the size of the time step ?). For PISO algorithm I tried many solution with high number of orthogonal and nonOrthogonal correctors (let's say respectively 300 and 20)...leading to more time expensive computation but not really better results. I also checked the mesh with the checkMesh utility:

Mesh non-orthogonality Max: 44.01408683 average: 11.32517179
Non-orthogonality check OK.


My question are the following:
1 - is it absolutely necessary to compute first a solution with potential foam (as proposed in one of the thread)?...I think it will not really respect the benchmark proposed I my case.
2 - Do you see any mistake in my procedure?
3 - Is they a way to link the non-orthogonality results obtained by checkMesh to the number of orthogonal and nonOrthogonal correctors.
4 - Other question : in the paper by Schäfer and Turek they speak about a recirculation area that as to be measured. What is this? Is they an easy way to automatically measure this in OpenFoam (not seeing the results in paraFoam)?

Thank you for any help or comments.
kassiotis is offline   Reply With Quote

Old   May 13, 2008, 06:19
Default 1) No, it is no need to comput
  #2
Senior Member
 
mkraposhin's Avatar
 
Matvey Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 355
Rep Power: 21
mkraposhin is on a distinguished road
1) No, it is no need to compute potential solution - this may only improve solution convergence
2) YOU SHOULD NOT USE icoFoam with steady-state time integration scheme, because icoFoam uses TRANSIENT prec-vel coupling algo, PISO
3) Your mesh thikness in Y-direction is too small, mesh thikness should be about 50-100 diamterers of cylinder (character length).

what is your BC's on upper and lower walls?

opt. number of non orthogonal correctors is 1-3

4) I can't donload papers by your references
mkraposhin is offline   Reply With Quote

Old   May 13, 2008, 06:47
Default the solver icoFoam isn't stati
  #3
New Member
 
Mark Michael
Join Date: Mar 2009
Location: Rostock, Germany
Posts: 5
Rep Power: 17
mark_michael is on a distinguished road
the solver icoFoam isn't stationary ! this solver needs an time integration scheme, because it solves the isothermal incompressible Navier-Stokes-equation ! so it also need time-steps !
Get some Books about "Fluidmechanics" and "wake flow around a circular cylinder" there you find many specifications about "recirculation area" !
Can you send me the paper per email ? I can't open the link ! astadfm@gmx.de
mark_michael is offline   Reply With Quote

Old   May 13, 2008, 07:00
Default Thank you for the hints. 1
  #4
New Member
 
Christophe Kassiotis
Join Date: Mar 2009
Location: Paris
Posts: 17
Rep Power: 17
kassiotis is on a distinguished road
Thank you for the hints.

1 - Ok, It was what I was thinking about potential solution.
2 - Not use icoFoam means use something like simpleFoam ?
3 - The mesh thikness, is, I admit, really not enough to avoid influence of boundary conditions on the flow around the cylinder, especially with the no-slip boundary condition I impose...but this is how the problem is defined in the paper (hope the link will work now).
kassiotis is offline   Reply With Quote

Old   May 14, 2008, 08:03
Default 2 - yes, use simpleFoam with t
  #5
Senior Member
 
mkraposhin's Avatar
 
Matvey Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 355
Rep Power: 21
mkraposhin is on a distinguished road
2 - yes, use simpleFoam with turbulence "off" and
turbulenceModel "laminar" (file constant/turbulenceProperties)

3 - for freestream boundary i'm using pressureInletOutletVelocity for U and outletInlet for p. And, of course, if you want to compare your results with those in paper, you must use BC's from paper

4 - the link work
mkraposhin is offline   Reply With Quote

Old   May 14, 2008, 12:09
Default Hello again, I try for the
  #6
New Member
 
Christophe Kassiotis
Join Date: Mar 2009
Location: Paris
Posts: 17
Rep Power: 17
kassiotis is on a distinguished road
Hello again,

I try for the moment two strategies in parallel, with, for both, initial state computed by potential foam solver.

with simpleFoam (steadyState for time integration, output specified through the use of inletOutlet, laminar flow without any turbulence as a model):

The results are the following





They are not so bad. Unfortunately in notice two problems :
  • The computation blows up after the first time step (does it really make sense to make more than one time step, as steadyState is for static case).
  • The drag coefficient obtain is around 82...far from the 5.5 awaited.

with icoFoam:

I run the computation with an Euler time integration scheme. This seems to behave as awaited. Unfortunately, the Drag coefficient exhibits the following shape when drawing respect to time (the value at the end is not so bad, as one awaited around 5.5).



I think this is linked to the computation of pressure, no ? When one plot the pressure respect to time, one see a kind of "slow" diffusion.

So my question are the following :
  • Can one rely on the computation at the first time step (for both icoFoam and simpleFoam) ?
  • Is the behavior shown here normal ?
  • How to avoid the computation divergence (if possible) for simpleFoam.

Thank you again for your help and advices !
kassiotis is offline   Reply With Quote

Old   May 15, 2008, 07:07
Default Hello everybody, In fact, i
  #7
New Member
 
Christophe Kassiotis
Join Date: Mar 2009
Location: Paris
Posts: 17
Rep Power: 17
kassiotis is on a distinguished road
Hello everybody,

In fact, it seems that, according to results from F. Bos computing the drag coefficient for a circular cylinder that the big overestimate at the beginning with icoFoam is a normal behavior before convergence

In the following, you will find different numerical experiments:
  • Time integration schemes: Euler, backward and Crank-Nicholson. Euler seems to give the better results.
  • Grad schemes: one can't notice significant differences between Gauss linear and least squares.
  • Time step: the results given are not really significant as for time step 0.001 (res. 0.005) one save is made each 10 (res. 2) times steps.
  • Initialization: use of potentialFoam really smooth the solution for the first time steps.









So, I think that I will use icoFoam, let it run a certain time and take the final value - which is close to the one expected at first glance - for the drag coefficient.

I remain interested by any comments and/or remarks.
kassiotis 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
reynold number on inlet Fayyaz Hussain FLUENT 2 June 18, 2007 13:24
LOW REYNOLD NUMBER Limonghi Main CFD Forum 1 March 13, 2007 12:37
urgent:low reynold number K-epsilon model prasanna FLUENT 3 January 1, 2007 14:20
low turbulent reynold number C-CARE CFX 4 March 12, 2004 07:36
Reynold number??? Jame D.J. CFX 5 July 30, 2002 09:27


All times are GMT -4. The time now is 18:33.