From c47c287d45f05fef7d050f9b1c044744ce065684 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Mon, 12 May 2014 23:15:27 +0100
Subject: [PATCH] DarcyForchheimer: add support for solvers which provide nu
rather than mu
---
.../DarcyForchheimer/DarcyForchheimer.C | 16 ++++++++++++++--
1 file changed, 14 insertions(+), 2 deletions(-)
diff --git a/src/finiteVolume/cfdTools/general/porosityModel/DarcyForchheimer/DarcyForchheimer.C b/src/finiteVolume/cfdTools/general/porosityModel/DarcyForchheimer/DarcyForchheimer.C
index 1ce0a83871d..54df8c5badf 100644
--- a/src/finiteVolume/cfdTools/general/porosityModel/DarcyForchheimer/DarcyForchheimer.C
+++ b/src/finiteVolume/cfdTools/general/porosityModel/DarcyForchheimer/DarcyForchheimer.C
@@ -203,9 +203,21 @@ void Foam::porosityModels::DarcyForchheimer::correct
if (UEqn.dimensions() == dimForce)
{
const volScalarField& rho = mesh_.lookupObject<volScalarField>(rhoName);
- const volScalarField& mu = mesh_.lookupObject<volScalarField>(muName);
- apply(Udiag, Usource, V, rho, mu, U);
+ if (mesh_.foundObject<volScalarField>(muName))
+ {
+ const volScalarField& mu =
+ mesh_.lookupObject<volScalarField>(muName);
+
+ apply(Udiag, Usource, V, rho, mu, U);
+ }
+ else
+ {
+ const volScalarField& nu =
+ mesh_.lookupObject<volScalarField>(nuName);
+
+ apply(Udiag, Usource, V, rho, rho*nu, U);
+ }
}
else
{
--
GitLab