adding the viscous force to the UEqn.H

This commit is contained in:
Behrad Esgandari
2023-04-26 13:39:33 +02:00
committed by GitHub
parent e4de205ef5
commit 3e3de6295b

View File

@ -4,6 +4,7 @@ tmp<fvVectorMatrix> tUEqn
( (
fvm::ddt(voidfraction,U) - fvm::Sp(fvc::ddt(voidfraction),U) fvm::ddt(voidfraction,U) - fvm::Sp(fvc::ddt(voidfraction),U)
+ fvm::div(phi,U) - fvm::Sp(fvc::div(phi),U) + fvm::div(phi,U) - fvm::Sp(fvc::div(phi),U)
+ particleCloud.divVoidfractionTau(U, voidfraction)
- fOther/rho - fOther/rho
== ==
fvOptions(U) fvOptions(U)