Commit e5b64f77 authored by sergio's avatar sergio Committed by Andrew Heather
Browse files

ENH: Calculating L from saturation states.

parent 43f904ba
......@@ -639,18 +639,19 @@ void alphatWallBoilingWallFunctionFvPatchScalarField::updateCoeffs()
const fvPatchScalarField& hew =
liquid.thermo().he().boundaryField()[patchi];
const scalarField hw
const scalarField hwLiqSat
(
liquid.thermo().he().member() == "e"
? hew.patchInternalField() + pw/rhow.patchInternalField()
: hew.patchInternalField()
? liquid.thermo().he(pw, Tsatc, patchi)
+ pw/rhow.patchInternalField()
: liquid.thermo().he(pw, Tsatc, patchi)
);
const scalarField L
(
vapor.thermo().he().member() == "e"
? vapor.thermo().he(pw, Tsatc, patchi) + pw/rhoVaporw - hw
: vapor.thermo().he(pw, Tsatc, patchi) - hw
? vapor.thermo().he(pw, Tsatc, patchi) + pw/rhoVaporw - hwLiqSat
: vapor.thermo().he(pw, Tsatc, patchi) - hwLiqSat
);
// Liquid phase fraction at the wall
......
......@@ -99,7 +99,7 @@ Foam::wallBoilingModels::CHFModels::Jeschar::MHF
Kmhf_*rhoVapor*L
*(
pow(sigma/(mag(g.value())*(rhoLiq - rhoVapor)), 0.25)
* sqrt(mag(g.value())*(rhoLiq - rhoVapor)/(rhoLiq + rhoVapor))
* sqrt(mag(g.value())*(rhoLiq - rhoVapor)/(rhoLiq + rhoVapor + VSMALL))
);
}
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment