CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   Writing own solver (

nico January 16, 2006 13:02

Hi! I want to add the Darcy-
I want to add the Darcy-term in the simpleFoam solver. I tried this:

tmp<fvvectormatrix> UEqn
fvm::div(phi, U)
+ turbulence->divR(U)+fvm::SuSp(nu*G,U)


solve(UEqn() == -fvc::grad(p));

Everything works fine until i changed my Field G from a volScalarField in a volTensorField, because i have directed permeabilities. I tried nu*(G & U), but it did not give the right solution.

Hope somebody can help me.


All times are GMT -4. The time now is 13:40.