CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Writing own solver (http://www.cfd-online.com/Forums/openfoam-solving/60348-writing-own-solver.html)

nico January 16, 2006 13:02

Hi! I want to add the Darcy-
 
Hi!
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)
);

UEqn().relax();

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.
Thanks

Nico


All times are GMT -4. The time now is 18:44.