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

strange behaviour of rhoPisoFoam, circular cylinder

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

Like Tree1Likes
  • 1 Post By ivan_cozza

Reply
 
LinkBack Thread Tools Display Modes
Old   May 25, 2010, 04:35
Default strange behaviour of rhoPisoFoam, circular cylinder
  #1
Senior Member
 
Ivan Flaminio Cozza
Join Date: Mar 2009
Location: Torino, Piemonte, Italia
Posts: 207
Rep Power: 9
ivan_cozza is on a distinguished road
Send a message via MSN to ivan_cozza
Hi Foamers,
I'm doing some compressible CFD on a 2D circular cylinder, Mach 0.12, Re 1.4x10^5, k-omega SST turbulence.
First I ran a wall-modeled simulations, y+ around 30, wall functions, (mutSpalartAllmarasWallFunction for mut, as suggested in other threads), and the run converged. I get an error of about 20% on St with respect 3D LES and experiments found in litterature, so I switched to low-Re modeling.
I remeshed up to y+ max of 1.2, and I set al turbulent variables at 1x10^-12 at wall, excepting for omega (omegaWallFunction should work also for low-Re in the 1.6.x version). I mapped from the previous run and I started the simulation, but now I can't obtain a stable result, the calculation blows up.

I did something wrong in the set-up? Did someone experienced similar problems?

Thanks, Ivan
ivan_cozza is offline   Reply With Quote

Old   May 25, 2010, 04:37
Default
  #2
Senior Member
 
Ivan Flaminio Cozza
Join Date: Mar 2009
Location: Torino, Piemonte, Italia
Posts: 207
Rep Power: 9
ivan_cozza is on a distinguished road
Send a message via MSN to ivan_cozza
For completeness, I use this schemes setup:


ddtSchemes
{
default backward;
}

gradSchemes
{
default cellMDLimited Gauss linear 1;
}

divSchemes
{
div(U,p) Gauss linearUpwind cellLimited Gauss linear 1;
div(phi,U) Gauss linearUpwind cellLimited Gauss linear 1;
div(phiU,p) Gauss linearUpwind cellLimited Gauss linear 1;
div((muEff*dev2(grad(U).T()))) Gauss linear;
div(phi,h) Gauss linearUpwind cellLimited Gauss linear 1;
div(phi,omega) Gauss linearUpwind cellLimited Gauss linear 1;
div(phi,k) Gauss linearUpwind cellLimited Gauss linear 1;
div(phid,p) Gauss linearUpwind cellLimited Gauss linear 1;
}

laplacianSchemes
{
laplacian(muEff,U) Gauss linear corrected;
laplacian(alphaEff,h) Gauss linear corrected;
laplacian((rho*rAU),p) Gauss linear corrected;
laplacian(DomegaEff,omega) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(1,p) Gauss linear corrected;
laplacian((rho*(1|A(U))),p) Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p ;
}
mm.abdollahzadeh likes this.
ivan_cozza is offline   Reply With Quote

Old   May 28, 2010, 05:21
Default
  #3
Senior Member
 
Ivan Flaminio Cozza
Join Date: Mar 2009
Location: Torino, Piemonte, Italia
Posts: 207
Rep Power: 9
ivan_cozza is on a distinguished road
Send a message via MSN to ivan_cozza
The story becomes more intricate:

I tried to do a longer run with the wall-modeled case, that up to 0.5 sec goes well. After a certain number of timesteps, it goes crazy like the low-Re model! I post some pictures of the problem.

First, when everything was ok:

Cl and Cd versus time of my cylinder




turbulent kinetic energy:



Omega:



log of the calculation:



Then, when the calculation go crazy

Cl and Cd versus time:



turbulent kinetic energy:



Omega:



log of the wrong calculation:



It seems that the dissipation of the turbulence model go crazy, destroying all the k in the simulation. The stranger thing is that I did not change anything between the two calculations, I just let the run go on with the same setup.
I have exactly the same problem with the low-Re mesh, the only difference is that this phenomenon appears in a fewer number of timesteps.

Please OpenFOAM gurus, give me some hints!

Have a nice day, Ivan
ivan_cozza is offline   Reply With Quote

Old   May 28, 2010, 08:46
Default
  #4
New Member
 
David Huckaby
Join Date: Jul 2009
Posts: 21
Rep Power: 7
dhuckaby is on a distinguished road
Ivan,
Have you tried reducing your Courant number or using another time-integration scheme such as Euler or boundedBackward ?

Dave
dhuckaby is offline   Reply With Quote

Old   May 28, 2010, 08:50
Default
  #5
Senior Member
 
Ivan Flaminio Cozza
Join Date: Mar 2009
Location: Torino, Piemonte, Italia
Posts: 207
Rep Power: 9
ivan_cozza is on a distinguished road
Send a message via MSN to ivan_cozza
Quote:
Originally Posted by dhuckaby View Post
Ivan,
Have you tried reducing your Courant number or using another time-integration scheme such as Euler or boundedBackward ?

Dave
Dave,
no I didn't, but my Max Courant in the whole calculation is below 0.9...
I have to try with boundedBackward... what's the difference between it and backward? Is more diffusive?
ivan_cozza is offline   Reply With Quote

Old   May 28, 2010, 09:30
Default
  #6
New Member
 
David Huckaby
Join Date: Jul 2009
Posts: 21
Rep Power: 7
dhuckaby is on a distinguished road
Ivan,

The general guidelines from various posts on the message board has been Co < 0.5 for stability Co < 0.2 for accurracy. I have also found that some simulations require a fixed time-step for stability as opposed to a Courant number bound.

The thread below briefly mentions the "boundedBackward" scheme
Schemes for scalar in LES, I would assume it is locally more diffusive.
dhuckaby is offline   Reply With Quote

Old   May 28, 2010, 09:57
Default
  #7
Senior Member
 
Ivan Flaminio Cozza
Join Date: Mar 2009
Location: Torino, Piemonte, Italia
Posts: 207
Rep Power: 9
ivan_cozza is on a distinguished road
Send a message via MSN to ivan_cozza
Quote:
Originally Posted by dhuckaby View Post
Ivan,

The general guidelines from various posts on the message board has been Co < 0.5 for stability Co < 0.2 for accurracy. I have also found that some simulations require a fixed time-step for stability as opposed to a Courant number bound.

The thread below briefly mentions the "boundedBackward" scheme
Schemes for scalar in LES, I would assume it is locally more diffusive.

Mmm... I'm not so experienced in unsteady simulations in OF, I used more frequently steady state, but Co < 0.5 for stability seems to me quite a severe limitation. But, I will try to limit my Co to less than 0.5...
ivan_cozza is offline   Reply With Quote

Old   December 19, 2010, 08:20
Default
  #8
Senior Member
 
aerothermal's Avatar
 
Guilherme da Silva
Join Date: Aug 2010
Location: Sao Paulo - Brazil
Posts: 105
Rep Power: 6
aerothermal is on a distinguished road
Hello Cozza,

have you managed to converge?
what was your fvSchemes?

I am trying to use the linearUpwind in OF16-ext and OF171 and they are giving me errors saying those schemes are not acceptable. Did they change the name in new versions?

Regards,

Guilherme
aerothermal 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
benchmark: flow over a circular cylinder goodegg Main CFD Forum 12 January 22, 2013 12:47
flow around a cylinder pXYZ Main CFD Forum 14 July 25, 2011 10:05
Solver and geometry choose for drag coefficient calculation around circular cylinder at large Re lin OpenFOAM Running, Solving & CFD 3 April 16, 2009 10:50
flow over a cylinder urgent! kevin FLUENT 7 June 8, 2006 01:19
Turbulent steady flow around a circular cylinder Mirek Kabacinski FLUENT 0 July 23, 2003 18:40


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