strange behaviour of rhoPisoFoam, circular cylinder
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?
For completeness, I use this schemes setup:
default cellMDLimited Gauss linear 1;
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;
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;
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:
log of the calculation:
Then, when the calculation go crazy
Cl and Cd versus time:
turbulent kinetic energy:
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
Have you tried reducing your Courant number or using another time-integration scheme such as Euler or boundedBackward ?
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?
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
http://www.cfd-online.com/Forums/ope...calar-les.html, 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...
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?
|All times are GMT -4. The time now is 09:01.|