-
Notifications
You must be signed in to change notification settings - Fork 40
Description
Dear authors:
Hello and thank you for providing such an excellent open source toolkit! I learned a lot from this kit.
I have a question about the code in the toolkit (using OpenFoam version V8), and there is a piece of code in the UEqn.H file
fvVectorMatrix UEqn
(
(1./eps)(fvm::ddt(rho, U) + fvm::div(rhoPhiByEps, U) + MRF.DDt(rho, U))
+ (1./eps)(turbulence->divDevTau(rho, U)) // Why is it multiplied by 1/eps?
+ fvm::Sp(Drag,U)
==
fvOptions(rho, U)
);
This code is a finite volume discrete assembly of the macroscopic momentum equation multiplied by a system of linear algebraic equations, which I think is the corresponding equation (refer to the formula (21) in the paper “Multiphase flow modeling in multiscale porous media:An open-source micro-continuum approach”).
One of the things I am confused about is why the viscous term
(1./eps)*(turbulence->divDevTau(rho, U))
I am looking forward to your reply and best wishes : )