CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   Explicit Solver and pressure correction (

hjasak October 5, 2006 00:39

You can't help it: in order to
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.


eugene October 5, 2006 07:16

I suggest you try and reformul
I suggest you try and reformulate your convection term into implicit part with an explicit correction.

Ustar = c*U_implicit + f(U_explicit)
phi_star = fvc::surfaceInterpolate(Ustar) <- can remain fully explicit


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.

All times are GMT -4. The time now is 21:32.