|
[Sponsors] |
October 5, 2006, 00:39 |
You can't help it: in order to
|
#1 |
Senior Member
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,905
Rep Power: 33 |
You can't help it: in order to satisfy mass conservation, you have to solve the pressure equation implicitly. Whatever else you do will give you checkerboarding.
As for the convection term in the momentum equation, thing about what it means in terms of cost to make the convection implicit. Also, using upwind differencing introduces a lot of discretisation error and will kill your LES. Hrv
__________________
Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk |
|
October 5, 2006, 07:16 |
I suggest you try and reformul
|
#2 |
Senior Member
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 21 |
I suggest you try and reformulate your convection term into implicit part with an explicit correction.
i.e. Ustar = c*U_implicit + f(U_explicit) phi_star = fvc::surfaceInterpolate(Ustar) <- can remain fully explicit thus: fvm::div(c*phi_star, U_implicit) + fvc::div(phi_star, f(U_explicit) (I offer no guarantee that this analysis is correct!) Of course this would mean redoing your gridFilter function to output the components of c and f(U_explicit), instead of just doing an explicit filter operation. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Explicit Solver | wessels | OpenFOAM Running, Solving & CFD | 1 | September 26, 2005 05:39 |
AMG SOLVER: PRESSURE CORRECTION | alex | FLUENT | 0 | August 12, 2005 13:17 |
Explicit Pressure Based Scheme | Apurva | Main CFD Forum | 4 | June 3, 2004 20:01 |
pressure driven flow by pressure correction method | justentered | Main CFD Forum | 0 | December 29, 2003 23:52 |
explicit CFD solver and multigrid | John Buckley | Main CFD Forum | 5 | March 8, 2000 05:00 |