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

OpenFOAM vs Fluent for cylinder at Re%3d150

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 28, 2006, 03:02
Default Hi Frank, I needed the mean
  #21
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,779
Rep Power: 22
hjasak will become famous soon enough
Hi Frank,

I needed the mean fixed pressure boundary condition for some other stuff... so I've now implemented one. For the case I'm looking at it behaves pretty well. If you'd like to re-visit your vortex shedding simulations with the new outlet b.c. please give me a shout and I'll pass over the code.

It bothers me a little bit that OpenFOAM did not perform the way I've expected (pride is a sin, I know) :-) and this may bring some closure to the problem...

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   April 28, 2006, 03:24
Default That would be very nice. I wil
  #22
Senior Member
 
Frank Bos
Join Date: Mar 2009
Location: The Netherlands
Posts: 340
Rep Power: 10
lr103476 is on a distinguished road
That would be very nice. I will just see how the new BC behaves for my problem.

While finishing a paper on previous work, I am still working on these cylinder flows using OpenFOAM, I varied grid and time step size. My main issues at this time are:

* I did not succeed in getting a grid independent solution.
* I like to verify if the bahaviour of the solver is indeed second order accurate (in space and time).

Possibly the new BC could shed some light on these problems. As soon as I get some fruitfull results I will post them on a website.

Regards, Frank
__________________
Frank Bos
lr103476 is offline   Reply With Quote

Old   April 28, 2006, 14:42
Default Hi Frank, Can you please po
  #23
Senior Member
 
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 703
Rep Power: 13
msrinath80 is on a distinguished road
Hi Frank,

Can you please post the source code you used to obtain the frequency of the flow. I'm trying a similar problem using Openfoam and since I'm not too good in C++, I'm looking for some help. I don't know exactly how to compile and use the probe, so some instructions should be useful.

Thanks very much for your help!
msrinath80 is offline   Reply With Quote

Old   May 1, 2006, 03:40
Default Just after the runTime.write()
  #24
Senior Member
 
Frank Bos
Join Date: Mar 2009
Location: The Netherlands
Posts: 340
Rep Power: 10
lr103476 is on a distinguished road
Just after the runTime.write() call in the icoFoam.C file, I put the following header file

#include "computeForces.H"

This file includes some code to calculate the lift, drag and moment and output the values to screen. This code is taken from the liftDrag utility.

When the lift/drag are written to screen, the foamLog routine is used to write the lift/drag histories to files. Finally these files are read into matlab for post-processing, the frequencies is only one thing I am looking at.

If you still got some difficulties with this, I could send you the source code by mail, if you want.

Regards, Frank
__________________
Frank Bos
lr103476 is offline   Reply With Quote

Old   May 1, 2006, 13:35
Default Fantastic! Thanks a lot. I'll
  #25
Senior Member
 
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 703
Rep Power: 13
msrinath80 is on a distinguished road
Fantastic! Thanks a lot. I'll give it a try as soon as I can.
msrinath80 is offline   Reply With Quote

Old   May 2, 2006, 07:51
Default Hi Frank, So all I need to
  #26
Senior Member
 
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 703
Rep Power: 13
msrinath80 is on a distinguished road
Hi Frank,

So all I need to do is edit icoFoam.C in the ~/OpenFOAM/OpenFOAM-1.3/applications/solvers/incompressible/icoFoam directory and include a reference to the computeForces header file?

Do I have to recompile OpenFoam to get this to work? Also, how does the foamLog routine work? Do I have to run it after the case has converged or before?

Thanks for your help!
msrinath80 is offline   Reply With Quote

Old   May 2, 2006, 09:00
Default You have to create your own cr
  #27
Senior Member
 
Frank Bos
Join Date: Mar 2009
Location: The Netherlands
Posts: 340
Rep Power: 10
lr103476 is on a distinguished road
You have to create your own createForces.H file. I simply took some stuff from the liftDrag utility, the part where the force are being calculated.

Furthermore I would create your own solver, for example: icoMySolverFoam in your user application directory. Try not to change the original code, instead create your own solvers, utilities and libraries. As a start copy and rename icoFoam to you user's application directory. Then try the stuff I mentioned before involving createForces.H etc...

Frank
__________________
Frank Bos
lr103476 is offline   Reply With Quote

Old   July 5, 2006, 15:07
Default Does anybody know why the Cour
  #28
Senior Member
 
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 703
Rep Power: 13
msrinath80 is on a distinguished road
Does anybody know why the Courant number is still significant. Isn't Crank Nicholson supposed to be an implicit scheme which is unconditionally stable (as Frank mentioned)?
msrinath80 is offline   Reply With Quote

Old   March 21, 2008, 22:36
Default I think both stability and acc
  #29
Senior Member
 
lakeat's Avatar
 
Daniel WEI (老魏)
Join Date: Mar 2009
Location: Beijing, China
Posts: 689
Blog Entries: 9
Rep Power: 13
lakeat is on a distinguished road
Send a message via Skype™ to lakeat
I think both stability and accuracy should be noticed.
Though I am little confused why Frank said,

"The max courant number which corresponds to the chosen time step (0.02) was between 1.7 and 1.9. I also performed computations with half of the time step (0.01) such that the max Courant = 0.8 - 0.9. In that case the vortex wake was also similar."

It should be a little little different in someplaces that may have not been awared of.



Daniel
__________________
~
Daniel WEI
-------------
Boeing Research & Technology - China
Beijing, China
Email
lakeat is offline   Reply With Quote

Old   March 21, 2008, 22:58
Default What?
  #30
Senior Member
 
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 703
Rep Power: 13
msrinath80 is on a distinguished road
What?
msrinath80 is offline   Reply With Quote

Old   September 16, 2008, 09:08
Default Hello Frank and hello Foamers,
  #31
Senior Member
 
lakeat's Avatar
 
Daniel WEI (老魏)
Join Date: Mar 2009
Location: Beijing, China
Posts: 689
Blog Entries: 9
Rep Power: 13
lakeat is on a distinguished road
Send a message via Skype™ to lakeat
Hello Frank and hello Foamers,

I have trouble in the flow over circular cylinder case, could you give me some advice?

1. May I know which reference are you reading for Re=150?
2. How did you set the inlet velocity profile?
3. Are you sure that the ddtSchemes with "CrankNicholson 0.5" is good for the case? I have compared the following
a) CrankNicholson 0.5;
b) CrankNicholson 0.9;
c) Euler;
d) backward;
and I found:
1) Courant number does not change the solution significantly as you said;
2) The Cd values of CrankNicholson are higher than those of Euler and backward, I don't know which are more accurate?

Not only that, my geometry is exactly as the one in Ferziger's book, see the pic below, but my average Cd (3.415) is higher than the book's (3.228)!
What might go wrong?



Thanks...

\Daniel
__________________
~
Daniel WEI
-------------
Boeing Research & Technology - China
Beijing, China
Email
lakeat is offline   Reply With Quote

Old   September 16, 2008, 09:20
Default http://www.cfd-online.com/Open
  #32
Senior Member
 
lakeat's Avatar
 
Daniel WEI (老魏)
Join Date: Mar 2009
Location: Beijing, China
Posts: 689
Blog Entries: 9
Rep Power: 13
lakeat is on a distinguished road
Send a message via Skype™ to lakeat



Thanks...

\Daniel
__________________
~
Daniel WEI
-------------
Boeing Research & Technology - China
Beijing, China
Email
lakeat is offline   Reply With Quote

Old   September 16, 2008, 09:23
Default http://www.cfd-online.com/Open
  #33
Senior Member
 
lakeat's Avatar
 
Daniel WEI (老魏)
Join Date: Mar 2009
Location: Beijing, China
Posts: 689
Blog Entries: 9
Rep Power: 13
lakeat is on a distinguished road
Send a message via Skype™ to lakeat



Thanks...

\Daniel
__________________
~
Daniel WEI
-------------
Boeing Research & Technology - China
Beijing, China
Email
lakeat is offline   Reply With Quote

Old   September 16, 2008, 09:41
Default Sorry for that, Chrome is stil
  #34
Senior Member
 
lakeat's Avatar
 
Daniel WEI (老魏)
Join Date: Mar 2009
Location: Beijing, China
Posts: 689
Blog Entries: 9
Rep Power: 13
lakeat is on a distinguished road
Send a message via Skype™ to lakeat
Sorry for that, Chrome is still green, I have to use IE. \Daniel
__________________
~
Daniel WEI
-------------
Boeing Research & Technology - China
Beijing, China
Email
lakeat is offline   Reply With Quote

Old   October 27, 2008, 06:24
Default Hi all, I'm currently using
  #35
New Member
 
Sébastien Bocquet
Join Date: Mar 2009
Posts: 21
Rep Power: 9
seb62 is on a distinguished road
Hi all,

I'm currently using openFoam on the flow around
a cylinder in order to better understand the code
parameters and post-processing functionalities.
I see that some of you are using a lift/drag utility. That would be a very useful tool which
doesn't seem to be implemented in the version I have (1.4).
Could someone tell me where I could find an updated version of the lift/drag utility?
(and possibly a short guide or description of
what it does).

Thanks very much for your help

Seb
seb62 is offline   Reply With Quote

Old   October 27, 2008, 11:27
Default Look at this link, http://www
  #36
Senior Member
 
lakeat's Avatar
 
Daniel WEI (老魏)
Join Date: Mar 2009
Location: Beijing, China
Posts: 689
Blog Entries: 9
Rep Power: 13
lakeat is on a distinguished road
Send a message via Skype™ to lakeat
Look at this link,
http://www.cfd-online.com/OpenFOAM_D...tml?1212507143

Good luck, Seb!

\Daniel
__________________
~
Daniel WEI
-------------
Boeing Research & Technology - China
Beijing, China
Email
lakeat is offline   Reply With Quote

Old   December 16, 2008, 10:13
Default Thanks Daniel, I am testing
  #37
New Member
 
Sébastien Bocquet
Join Date: Mar 2009
Posts: 21
Rep Power: 9
seb62 is on a distinguished road
Thanks Daniel,

I am testing a similar case,flow around a cylinder at Re=100.
I'm using icoFoam, and I'm looking at the drag coefficient and vortex shedding frequency.
After checking that the results are time step independent, I varied the divScheme.
I found that the best results (providing my experimental results that I'm using as a reference are accurate enough) are obtained using
Gauss Linear, and not the upwind scheme. For the moment, I only tried these two, and I just find a bit strange that Gauss Linear gives me better results. Is this scheme centered?

Thanks for any comment you could make on this,
it will be helpful.
seb62 is offline   Reply With Quote

Old   December 16, 2008, 10:33
Default Hi Sebastien, Try also the
  #38
Senior Member
 
Frank Bos
Join Date: Mar 2009
Location: The Netherlands
Posts: 340
Rep Power: 10
lr103476 is on a distinguished road
Hi Sebastien,

Try also the vanLeerV scheme for convection. This gives nice vortex preserving results for me.

Frank
__________________
Frank Bos
lr103476 is offline   Reply With Quote

Old   December 16, 2008, 11:05
Default Hi Frank! I want to do a si
  #39
Senior Member
 
lakeat's Avatar
 
Daniel WEI (老魏)
Join Date: Mar 2009
Location: Beijing, China
Posts: 689
Blog Entries: 9
Rep Power: 13
lakeat is on a distinguished road
Send a message via Skype™ to lakeat
Hi Frank!

I want to do a simulation of FSI of a bridge deck (A very high Re), could you recommend me some papers on FSI.

{I have paid a lot of focuses on DES and DDES tech for the last several months, but I'm afraid they are still unaffordable. }

Thank you very much.

Daniel
__________________
~
Daniel WEI
-------------
Boeing Research & Technology - China
Beijing, China
Email
lakeat is offline   Reply With Quote

Old   December 16, 2008, 11:16
Default Hi Sebastien, Gauss Linear i
  #40
Member
 
Pierre Le Fur
Join Date: Mar 2009
Location: UK
Posts: 60
Rep Power: 9
pierre is on a distinguished road
Hi Sebastien,
Gauss Linear is simply central differencing. If it does not result in any numerical instability for you then why not keep it. 2nd order is all good.
Frank recommends vanLeer, which simply blends with upwind, for stability's sake, in which case limitedlinearV is not bad either...

Ciao

Pierre
pierre 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
drag coefficient of cylinder using fluent john lynam Main CFD Forum 2 February 15, 2010 10:16
RE:Change roughness of a cylinder in FLUENT!please Ken FLUENT 1 April 3, 2007 03:12
Fluent-Flow over heated cylinder Khyati Main CFD Forum 0 July 12, 2006 00:17
cylinder grid for fluent Wolfgang Main CFD Forum 1 October 21, 2002 12:11
FLUENT-FLOW AROUND A CIRCULAR CYLINDER! Paul Bennett Main CFD Forum 8 March 22, 1999 14:00


All times are GMT -4. The time now is 00:14.