CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Implicit source term and porousInterFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   August 30, 2013, 04:30
Default Implicit source term and porousInterFoam
Senior Member
Andrea Ferrari
Join Date: Dec 2010
Posts: 314
Rep Power: 10
Andrea_85 is on a distinguished road
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?

Andrea_85 is offline   Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

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