songtianlun/openfoam
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

ENH: atmPlantCanopyUSource: improve implicit behaviour (fixes #2343)

Browse Source
This commit is contained in:
Nima Samkhaniani 2022-03-30 11:37:49 +01:00 committed by Andrew Heather
parent 123fe09fda
commit 583fc4fb0d
1 changed files with 5 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
src/atmosphericModels/fvOptions/atmPlantCanopyUSource/atmPlantCanopyUSource.C
View File

@ -6,7 +6,7 @@
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
Copyright (C) 2020 ENERCON GmbH Copyright (C) 2020 ENERCON GmbH
Copyright (C) 2020-2021 OpenCFD Ltd. Copyright (C) 2020-2022 OpenCFD Ltd.
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
This file is part of OpenFOAM. This file is part of OpenFOAM.
@ -28,6 +28,7 @@ License
#include "atmPlantCanopyUSource.H" #include "atmPlantCanopyUSource.H"
#include "addToRunTimeSelectionTable.H" #include "addToRunTimeSelectionTable.H"
#include "fvmSup.H"
// * * * * * * * * * * * * * Static Member Functions * * * * * * * * * * * * // // * * * * * * * * * * * * * Static Member Functions * * * * * * * * * * * * //
@ -99,7 +100,7 @@ void Foam::fv::atmPlantCanopyUSource::addSup
if (V_ > VSMALL) if (V_ > VSMALL)
{ {
// (SP:Eq. 42) // (SP:Eq. 42)
eqn -= (plantCd_*leafAreaDensity_*mag(U))*U; eqn -= fvm::Sp(plantCd_*leafAreaDensity_*mag(U), U);
} }
} }
@ -115,7 +116,7 @@ void Foam::fv::atmPlantCanopyUSource::addSup
if (V_ > VSMALL) if (V_ > VSMALL)
{ {
eqn -= rho*(plantCd_*leafAreaDensity_*mag(U))*U; eqn -= fvm::Sp(rho*plantCd_*leafAreaDensity_*mag(U), U);
} }
} }
@ -132,7 +133,7 @@ void Foam::fv::atmPlantCanopyUSource::addSup
if (V_ > VSMALL) if (V_ > VSMALL)
{ {
eqn -= alpha*rho*(plantCd_*leafAreaDensity_*mag(U))*U; eqn -= fvm::Sp(alpha*rho*plantCd_*leafAreaDensity_*mag(U), U);
} }
} }