zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

BUG: interfaceHeatResistance: remove max condensation limiter

Browse Source
This commit is contained in:
sergio
2022-05-25 14:45:17 -07:00
committed by Kutalmis Bercin
parent 8f8612fd0c
commit 0f0a5c5c0f
1 changed files with 2 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files

14
applications/solvers/multiphase/interCondensatingEvaporatingFoam/temperaturePhaseChangeTwoPhaseMixtures/interfaceHeatResistance/interfaceHeatResistance.C
View File

@ -5,7 +5,7 @@
\\ / A nd | www.openfoam.com \\ / A nd | www.openfoam.com
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
Copyright (C) 2020-2021 OpenCFD Ltd. Copyright (C) 2020-2022 OpenCFD Ltd.
Copyright (C) 2020 Henning Scheufler Copyright (C) 2020 Henning Scheufler
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -265,23 +265,13 @@ correct()
const dimensionedScalar& TSat = thermo.TSat(); const dimensionedScalar& TSat = thermo.TSat();
const dimensionedScalar T0(dimTemperature, Zero); const dimensionedScalar T0(dimTemperature, Zero);
dimensionedScalar L = mixture_.Hf2() - mixture_.Hf1(); const dimensionedScalar L(mag(mixture_.Hf2() - mixture_.Hf1()));
// interface heat resistance // interface heat resistance
mDotc_ = interfaceArea_*R_*max(TSat - T, T0)/L; mDotc_ = interfaceArea_*R_*max(TSat - T, T0)/L;
mDote_ = interfaceArea_*R_*max(T - TSat, T0)/L; mDote_ = interfaceArea_*R_*max(T - TSat, T0)/L;
// Limiting max condensation
forAll(mDotc_, celli)
{
scalar rhobyDt = mixture_.rho1().value()/mesh_.time().deltaTValue();
scalar maxEvap = mixture_.alpha1()[celli]*rhobyDt; // positive
scalar maxCond = -mixture_.alpha2()[celli]*rhobyDt; // negative
mDotc_[celli] = min(max(mDotc_[celli], maxCond), maxEvap);
}
// Calculate the spread sources // Calculate the spread sources
dimensionedScalar D dimensionedScalar D
( (
"D", "D",