Implicit source term and porousInterFoam
 

Hi all,

I want to add an additional implicit source term in the momentum equation (using interFoam). The term is simply a darcy term (=mu/K * U) to account for additional resistance. I would have done it in the following way: in UEqn.H add

fvVectorMatrix UEqn
(
fvm::ddt(rho, U)
+ fvm::div(rhoPhi, U)
- fvm::laplacian(muEff, U)
- (fvc::grad(U) & fvc::grad(muEff))
//- fvc::div(muEff*(fvc::interpolate(dev(fvc::grad(U)) ) & mesh.Sf()))
+fvm::Sp(mu/K, U)
);

with K defined as volVectorField and constant for each cell.

Then i noticed this porousInterFoam solver. Is it doing the same thing? The additional resistance term is treated implicit or explicit?
Is there any tutorial for this solver?


Best
andrea