ENH: Create output for Su VollerPrakash

This commit is contained in:
sergio
2022-04-22 10:01:32 -07:00
parent ecd27ad4a4
commit 4f7019edc3

View File

@ -5,7 +5,7 @@
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2017 OpenCFD Ltd.
Copyright (C) 2017-2020 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -68,6 +68,16 @@ Foam::porousModels::VollerPrakash::S() const
const volScalarField& solidAlpha =
mesh_.lookupObject<volScalarField>(solidPhase_);
if (mesh_.time().outputTime() && debug)
{
volScalarField Su
(
"Su" + solidAlpha.name(),
Cu_*sqr(solidAlpha)/(pow3(1.0 - solidAlpha) + 1e-3)
);
Su.write();
}
return Cu_*sqr(solidAlpha)/(pow3(1.0 - solidAlpha) + 1e-3);
}