|August 30, 2013, 04:30||
Implicit source term and porousInterFoam
Join Date: Dec 2010
Posts: 275Rep Power: 7
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
+ 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()))
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?