diff --git a/src/lagrangian/parcel/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C b/src/lagrangian/parcel/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C index daa306ae10..4f86af6a7f 100644 --- a/src/lagrangian/parcel/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C +++ b/src/lagrangian/parcel/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C @@ -72,6 +72,14 @@ Foam::LiquidEvaporation::LiquidEvaporation : PhaseChangeModel(dict, owner, typeName), liquids_(owner.thermo().liquids()), + condensation_ + ( + this->coeffDict().template lookupOrDefault + ( + "condensation", + false + ) + ), activeLiquids_(this->coeffDict().lookup("activeLiquids")), liqToCarrierMap_(activeLiquids_.size(), -1), liqToLiqMap_(activeLiquids_.size(), -1) @@ -200,7 +208,13 @@ void Foam::LiquidEvaporation::calculate const scalar Cinf = Xc[gid]*pc/(RR*Ts); // molar flux of vapour [kmol/m2/s] - const scalar Ni = max(kc*(Cs - Cinf), 0.0); + scalar Ni = kc*(Cs - Cinf); + + // limit if not permitting condensation + if (!condensation_) + { + Ni = max(Ni, 0); + } // mass transfer [kg] dMassPC[lid] += Ni*pi*sqr(d)*liquids_.properties()[lid].W()*dt; diff --git a/src/lagrangian/parcel/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.H b/src/lagrangian/parcel/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.H index 1b4218a154..09c80615a2 100644 --- a/src/lagrangian/parcel/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.H +++ b/src/lagrangian/parcel/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.H @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org - \\ / A nd | Copyright (C) 2011-2020 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2021 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -56,6 +56,9 @@ protected: //- Global liquid properties data const liquidMixtureProperties& liquids_; + //- Whether or not to permit condensation (default false) + Switch condensation_; + //- List of active liquid names List activeLiquids_;