From 3e3de6295b4971a2ddc6ce9dc6bdd561ba85c488 Mon Sep 17 00:00:00 2001 From: Behrad Esgandari <93246630+behradesg@users.noreply.github.com> Date: Wed, 26 Apr 2023 13:39:33 +0200 Subject: [PATCH] adding the viscous force to the UEqn.H --- applications/solvers/cfdemSolverPimple/UEqn.H | 1 + 1 file changed, 1 insertion(+) diff --git a/applications/solvers/cfdemSolverPimple/UEqn.H b/applications/solvers/cfdemSolverPimple/UEqn.H index 022516ac..d72c8c23 100644 --- a/applications/solvers/cfdemSolverPimple/UEqn.H +++ b/applications/solvers/cfdemSolverPimple/UEqn.H @@ -4,6 +4,7 @@ tmp tUEqn ( fvm::ddt(voidfraction,U) - fvm::Sp(fvc::ddt(voidfraction),U) + fvm::div(phi,U) - fvm::Sp(fvc::div(phi),U) + + particleCloud.divVoidfractionTau(U, voidfraction) - fOther/rho == fvOptions(U)