diff --git a/src/phaseSystemModels/multiphaseInter/phasesSystem/MassTransferPhaseSystem/MassTransferPhaseSystem.C b/src/phaseSystemModels/multiphaseInter/phasesSystem/MassTransferPhaseSystem/MassTransferPhaseSystem.C index 5de5b50e70bae998f98d724bb9a5ceeaf62e8867..559fee3e63915cb5a225bf687a6b77a475a96843 100644 --- a/src/phaseSystemModels/multiphaseInter/phasesSystem/MassTransferPhaseSystem/MassTransferPhaseSystem.C +++ b/src/phaseSystemModels/multiphaseInter/phasesSystem/MassTransferPhaseSystem/MassTransferPhaseSystem.C @@ -121,7 +121,7 @@ Foam::MassTransferPhaseSystem<BasePhaseSystem>::calculateL const word species(speciesName.substr(0, speciesName.find('.'))); - L += pos(dmdtNetki)*interfacePtr->L(species, T); + L -= pos(dmdtNetki)*interfacePtr->L(species, T); } return tL; @@ -240,7 +240,7 @@ Foam::MassTransferPhaseSystem<BasePhaseSystem>::heatTransfer if (KSp.valid()) { - Sp -= KSp.ref(); + Sp += KSp.ref(); } tmp<volScalarField> KSu = @@ -248,13 +248,13 @@ Foam::MassTransferPhaseSystem<BasePhaseSystem>::heatTransfer if (KSu.valid()) { - Su -= KSu.ref(); + Su += KSu.ref(); } // If linearization is not provided used full explicit if (!KSp.valid() && !KSu.valid()) { - Su -= *dmdt_[keyik]; + Su += *dmdt_[keyik]; } } diff --git a/src/phaseSystemModels/multiphaseInter/phasesSystem/massTransferModels/Lee/Lee.C b/src/phaseSystemModels/multiphaseInter/phasesSystem/massTransferModels/Lee/Lee.C index b8313ea33e3eb9a0f16778257d9d45fe95e4725a..d0cde4679e2a3b2276052e3d248ebc12a5a1b79d 100644 --- a/src/phaseSystemModels/multiphaseInter/phasesSystem/massTransferModels/Lee/Lee.C +++ b/src/phaseSystemModels/multiphaseInter/phasesSystem/massTransferModels/Lee/Lee.C @@ -78,7 +78,7 @@ Foam::meltingEvaporationModels::Lee<Thermo, OtherThermo>::Kexp { return ( - coeff*pos(Tactivate_ - refValue) + coeff*pos(Tactivate_ - refValue) ); } } diff --git a/src/phaseSystemModels/multiphaseInter/phasesSystem/massTransferModels/interfaceHeatResistance/interfaceHeatResistance.C b/src/phaseSystemModels/multiphaseInter/phasesSystem/massTransferModels/interfaceHeatResistance/interfaceHeatResistance.C index ca76c43a51e0ed96e173105a34d29f48f1ab6d30..0dfd7e8b117b3af5044a07a8b66ace78ca4e79c5 100644 --- a/src/phaseSystemModels/multiphaseInter/phasesSystem/massTransferModels/interfaceHeatResistance/interfaceHeatResistance.C +++ b/src/phaseSystemModels/multiphaseInter/phasesSystem/massTransferModels/interfaceHeatResistance/interfaceHeatResistance.C @@ -6,7 +6,7 @@ \\/ M anipulation | ------------------------------------------------------------------------------- Copyright (C) 2020 Henning Scheufler - Copyright (C) 2020-2021 OpenCFD Ltd. + Copyright (C) 2020-2022 OpenCFD Ltd. ------------------------------------------------------------------------------- License This file is part of OpenFOAM. @@ -223,7 +223,7 @@ Foam::meltingEvaporationModels::interfaceHeatResistance<Thermo, OtherThermo> auto tempOpen = fullSpeciesName.find('.'); const word speciesName(fullSpeciesName.substr(0, tempOpen)); - tmp<volScalarField> L = this->L(speciesName, T); + tmp<volScalarField> L = mag(this->L(speciesName, T)); htc_ = R_/L(); diff --git a/src/phaseSystemModels/multiphaseInter/phasesSystem/massTransferModels/kineticGasEvaporation/kineticGasEvaporation.C b/src/phaseSystemModels/multiphaseInter/phasesSystem/massTransferModels/kineticGasEvaporation/kineticGasEvaporation.C index 7b89464004e46d45aea0ff3d08e634820d9a683f..a0a45133992337acff8f38736be8d4ce94b059d0 100644 --- a/src/phaseSystemModels/multiphaseInter/phasesSystem/massTransferModels/kineticGasEvaporation/kineticGasEvaporation.C +++ b/src/phaseSystemModels/multiphaseInter/phasesSystem/massTransferModels/kineticGasEvaporation/kineticGasEvaporation.C @@ -5,7 +5,7 @@ \\ / A nd | www.openfoam.com \\/ M anipulation | ------------------------------------------------------------------------------- - Copyright (C) 2017-2021 OpenCFD Ltd. + Copyright (C) 2017-2022 OpenCFD Ltd. ------------------------------------------------------------------------------- License This file is part of OpenFOAM. @@ -206,7 +206,7 @@ Foam::meltingEvaporationModels::kineticGasEvaporation<Thermo, OtherThermo> ); word speciesName = IOobject::member(this->transferSpecie()); - tmp<volScalarField> L = this->L(speciesName, T); + tmp<volScalarField> L = mag(this->L(speciesName, T)); updateInterface(T); diff --git a/tutorials/multiphase/icoReactingMultiphaseInterFoam/evaporationMultiComponent/constant/thermophysicalProperties.liquid b/tutorials/multiphase/icoReactingMultiphaseInterFoam/evaporationMultiComponent/constant/thermophysicalProperties.liquid index 8f35a391441839d5ca32df492593cf558a686076..a12acc2ec7857d9525e5e47b2295590438fe2ac5 100644 --- a/tutorials/multiphase/icoReactingMultiphaseInterFoam/evaporationMultiComponent/constant/thermophysicalProperties.liquid +++ b/tutorials/multiphase/icoReactingMultiphaseInterFoam/evaporationMultiComponent/constant/thermophysicalProperties.liquid @@ -39,7 +39,7 @@ mixture thermodynamics { Cp 4181; - Hf 2.45e6;//-1.5833e7;//deltaHv 2.45e6; //[J/Kg] + Hf -2.45e6; //[J/Kg] } transport { diff --git a/tutorials/multiphase/icoReactingMultiphaseInterFoam/mixerVesselAMI2D/constant/thermophysicalProperties.liquid b/tutorials/multiphase/icoReactingMultiphaseInterFoam/mixerVesselAMI2D/constant/thermophysicalProperties.liquid index 7c58fc3403ba1d0a3cdd6e29a19516d8cf5ef77e..104d664c5c1d270da3fb60890a184afc4c893cf7 100644 --- a/tutorials/multiphase/icoReactingMultiphaseInterFoam/mixerVesselAMI2D/constant/thermophysicalProperties.liquid +++ b/tutorials/multiphase/icoReactingMultiphaseInterFoam/mixerVesselAMI2D/constant/thermophysicalProperties.liquid @@ -40,7 +40,7 @@ mixture thermodynamics { Cp 4181; - Hf 2.45e6; //[J/Kg] + Hf -2.45e6; //[J/Kg] } transport { diff --git a/tutorials/multiphase/icoReactingMultiphaseInterFoam/poolEvaporation/constant/thermophysicalProperties.liquid b/tutorials/multiphase/icoReactingMultiphaseInterFoam/poolEvaporation/constant/thermophysicalProperties.liquid index f1d85ecdca54d1900a9f51da2a08f07a2e2e3cec..ea44fd872f33aaafae29871315effdad50461cdb 100644 --- a/tutorials/multiphase/icoReactingMultiphaseInterFoam/poolEvaporation/constant/thermophysicalProperties.liquid +++ b/tutorials/multiphase/icoReactingMultiphaseInterFoam/poolEvaporation/constant/thermophysicalProperties.liquid @@ -39,7 +39,7 @@ mixture thermodynamics { Cp 4216; - Hf 2.45e6; + Hf -2.45e6; } transport { diff --git a/tutorials/multiphase/icoReactingMultiphaseInterFoam/solidMelting2D/constant/thermophysicalProperties.liquid b/tutorials/multiphase/icoReactingMultiphaseInterFoam/solidMelting2D/constant/thermophysicalProperties.liquid index 6f12cd4b91b4bd071ad130244ecc635d7d810936..2de7a9a9bab3640974018077fe24965b788b5940 100644 --- a/tutorials/multiphase/icoReactingMultiphaseInterFoam/solidMelting2D/constant/thermophysicalProperties.liquid +++ b/tutorials/multiphase/icoReactingMultiphaseInterFoam/solidMelting2D/constant/thermophysicalProperties.liquid @@ -41,7 +41,7 @@ mixture thermodynamics { Cp 381; - Hf 0; + Hf 8.016e4; } transport { diff --git a/tutorials/multiphase/icoReactingMultiphaseInterFoam/solidMelting2D/constant/thermophysicalProperties.solid b/tutorials/multiphase/icoReactingMultiphaseInterFoam/solidMelting2D/constant/thermophysicalProperties.solid index 590032365eb7a0220b10924404b79134004f3627..2c6e5a5e93bfe1ad0e943fe0e7c463e3a10e5fd0 100644 --- a/tutorials/multiphase/icoReactingMultiphaseInterFoam/solidMelting2D/constant/thermophysicalProperties.solid +++ b/tutorials/multiphase/icoReactingMultiphaseInterFoam/solidMelting2D/constant/thermophysicalProperties.solid @@ -39,7 +39,7 @@ mixture thermodynamics { Cp 381; - Hf 8.016e4;//[J/Kg]; + Hf 0; //[J/Kg]; } transport { @@ -48,4 +48,4 @@ mixture } -// ************************************************************************* // \ No newline at end of file +// ************************************************************************* //