Commit 998f8f3a authored by Henry Weller's avatar Henry Weller Committed by Andrew Heather
Browse files

thermoSingleLayer::q: Use a local "alpha" without hydrophilic/phobic adjustment

Resolves bug-report https://bugs.openfoam.org/view.php?id=2605
parent 38471767
......@@ -208,7 +208,7 @@ tmp<volScalarField> kinematicSingleLayer::pp()
void kinematicSingleLayer::correctAlpha()
{
alpha_ == pos0(delta_ - deltaSmall_);
alpha_ == pos(delta_ - deltaSmall_);
}
......
......@@ -269,21 +269,17 @@ void thermoSingleLayer::updateSubmodels()
tmp<fvScalarMatrix> thermoSingleLayer::q(volScalarField& hs) const
{
// Only apply heat transfer where the film is present
// - leads to temperature unboundedness?
// volScalarField boundedAlpha(max(alpha_, ROOTVSMALL));
// volScalarField htcst(htcs_->h()*boundedAlpha);
// volScalarField htcwt(htcw_->h()*boundedAlpha);
const volScalarField alpha(pos(delta_ - deltaSmall_));
return
(
// Heat-transfer to the primary region
- fvm::Sp(htcs_->h()/Cp_, hs)
+ htcs_->h()*(hs/Cp_ + alpha_*(TPrimary_ - T_))
+ htcs_->h()*(hs/Cp_ + alpha*(TPrimary_ - T_))
// Heat-transfer to the wall
- fvm::Sp(htcw_->h()/Cp_, hs)
+ htcw_->h()*(hs/Cp_ + alpha_*(Tw_- T_))
+ htcw_->h()*(hs/Cp_ + alpha*(Tw_- T_))
);
}
......
Supports Markdown
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