void thermoSingleLayer::evolveRegion()
{
if (debug)
{
tabAdd();
Info<<tab.c_str()<< "thermoSingleLayer::evolveRegion()" << endl;
}
//kvm, if updateSubmodels occurs after solveContinuity, then vaporization and separation don't get called until after deltaRho_ has been calculated, resulting in errors
updateSubmodels();
// Solve continuity for deltaRho_
solveContinuity();
for (int oCorr=0; oCorr<nOuterCorr_; oCorr++)
{
// Explicit pressure source contribution
tmp<volScalarField> tpu(this->pu());
// Implicit pressure source coefficient
tmp<volScalarField> tpp(this->pp());
// Solve for momentum for U_
tmp<fvVectorMatrix> UEqn = solveMomentum(tpu(), tpp());
// Solve energy for hs_ - also updates thermo
solveEnergy();
// Film thickness correction loop
for (int corr=1; corr<=nCorr_; corr++)
{
// Solve thickness for delta_
solveThickness(tpu(), tpp(), UEqn());
}
}
T_ == T(hs_);
#include "diagnostics.H"
// Update deltaRho_ with new delta_
deltaRho_ == delta_*rho_;
// Update film wall and surface velocities
updateSurfaceVelocities();
// Update film wall and surface temperatures
// only need to call once?
// updateSurfaceTemperatures();
// Reset source terms for next time integration
resetPrimaryRegionSourceTerms();
if (debug)
{
Info<<tab.c_str()<< "leaving thermoSingleLayer::evolveRegion()" << endl;
tabSubtract();
}
}
void kinematicSingleLayer::evolveRegion()
{
if (debug)
{
InfoInFunction << endl;
}
// Update film coverage indicator
correctAlpha();
// Update film wall and surface velocities
updateSurfaceVelocities();
// Update sub-models to provide updated source contributions
updateSubmodels();
// Solve continuity for deltaRho_
solveContinuity();
// Implicit pressure source coefficient - constant
tmp<volScalarField> tpp(this->pp());
for (int oCorr=1; oCorr<=nOuterCorr_; oCorr++)
{
// Explicit pressure source contribution - varies with delta_
tmp<volScalarField> tpu(this->pu());
// Solve for momentum for U_
tmp<fvVectorMatrix> UEqn = solveMomentum(tpu(), tpp());
// Film thickness correction loop
for (int corr=1; corr<=nCorr_; corr++)
{
// Solve thickness for delta_
solveThickness(tpu(), tpp(), UEqn());
}
}
// Update deltaRho_ with new delta_
deltaRho_ == delta_*rho_;
// Reset source terms for next time integration
resetPrimaryRegionSourceTerms();
}