CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   Another way to solve Stokes (

lam June 25, 2007 09:57

Hi all, I want to use an ot
Hi all,

I want to use an other way to solve a steady Stokes equation for incompressible flow:

laplacian(nu,U)=Grad(p) (2)


div(U)=epsi*p, where epsi is very small
i put this expression in (2), and I have:

Laplacian(nu,U)=grad(div(U)/epsi) (2')

But my problem is I have not managed to discretize "grad(div(U)/epsi)"

I've tried :

tmp<fvvectormatrix> UEqn (fvm::laplacian(nuEff(), U));

But it tells me someting like in compilation:
error: no matching function for call to 'div(Foam::volVectorField&)

I think it a intressing way to test, even if it could be expensive for the memory.

Hope taht anybody will give me some advice to make a good discretisation for (2')

Thanks in advance,


lam June 25, 2007 12:15

nobody has an idea?
nobody has an idea?

lam June 26, 2007 05:54

I have tried this: tmp
I have tried this:

tmp<fvvectormatrix> UEqn (turbulence->divR(U));
where epsi is a dimensioned scalar [0 -2 1 0 0 0 0]

But my solution doesn't obey the continuity!

Can somebody tell me where I am wrong?

All times are GMT -4. The time now is 10:37.