|
[Sponsors] |
September 9, 2013, 04:58 |
Higher order Schemes are unstable, Why?
|
#1 |
Senior Member
|
Hi,
I used Gauss Linear for pressure gradient, the solution is stable, but the results are bit far from the expected, so switched Gauss Linear pressure gradient to leastSquares, the simulation is not stable, to make it stable I have given some limitation to the laplacianSchemes (0.5) for other than pressure, now the solution is stable but the results are not acceptable, it far away from the experimental results. gradSchemes { default Gauss linear; grad(p) leastSquares; grad(omega) Gauss linear; grad(epsilon) Gauss linear; grad(U) Gauss linear; grad(k) Gauss linear; } divSchemes { default none; div(phi,U) Gauss linearUpwindV grad(U); div(phi,k) Gauss upwind; div(phi,omega) Gauss upwind; div(phi,epsilon) Gauss upwind ; div((nuEff*dev(T(grad(U))))) Gauss linear; } laplacianSchemes { default none; laplacian(nuEff,U) Gauss linear limited 0.5; laplacian((1|A(U)),p) Gauss linear limited 1; laplacian(1,p) Gauss linear limited 1; laplacian(DkEff,k) Gauss linear limited 0.5; laplacian(DepsilonEff,epsilon) Gauss linear limited 0.5; laplacian(DomegaEff,omega) Gauss linear limited 0.5; } interpolationSchemes { default linear; interpolate(U) linear; } snGradSchemes { default limited 0.5; } is there any way to solve this problem? any suggestions? Thanks, Sivakumar |
|
September 10, 2013, 09:32 |
|
#2 |
Senior Member
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,905
Rep Power: 33 |
You need to read some CFD books: you are using unbounded convection schemes on bounded variables: this is why the simulations blow up for you.
Hrv
__________________
Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk |
|
September 10, 2013, 11:00 |
|
#3 |
Senior Member
|
Dear Prof (Hrv),
Thanks for your reply and suggestion. Thanks, Sivakumar |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Higher order formulation question | peterciaran | Main CFD Forum | 1 | July 27, 2012 13:47 |
2nd order boundary conditions for 2nd order discretization? | quarkz | Main CFD Forum | 30 | December 26, 2011 07:12 |
divergence with higher order scheme | shekharc | Main CFD Forum | 1 | July 23, 2009 13:53 |
4th order schemes in channelOodles | maka | OpenFOAM Bugs | 9 | January 19, 2009 11:58 |
Higher order FVM | Sergey Smirnov | Main CFD Forum | 10 | April 15, 2000 01:49 |