ENH: Updated regionModels from internal development line
Merge request reports
Activity
Filter activity
282 305 + fvm::div(phi_, hs_) 283 306 == 284 307 - hsSp_ 285 + q(hs_) 286 + radiation_->Shs() 308 + alpha_*q(hs_) Need to check with @Sergio - this was added to enable the pyrolysis and film to act simultaneously
282 305 + fvm::div(phi_, hs_) 283 306 == 284 307 - hsSp_ 285 + q(hs_) 286 + radiation_->Shs() 308 + alpha_*q(hs_) 282 305 + fvm::div(phi_, hs_) 283 306 == 284 307 - hsSp_ 285 + q(hs_) 286 + radiation_->Shs() 308 + alpha_*q(hs_) My question was about alpha in general. alpha is clipped to be either 1 or 0 (if some height gets low?)
You add the source only where alpha=1. Where alpha=0 you do not apply the source term. Hence you're throwing away that term (which might be needed for conservation). So now you've moved an accumulation error (height low) into a conservation error.
Added 1 commit:
- 40a88a14 - ENH: Changing p_rgh BC and call to solidThermo in noPyrolysis.C
Added 196 commits:
-
40a88a14...180e69c6 - 194 commits from branch
develop
- ea642d39 - Merge remote-tracking branch 'origin/develop' into feature-regionModels
- 1837ec57 - STYLE: sigFpe: use of FatalErrorInFunction
-
40a88a14...180e69c6 - 194 commits from branch
mentioned in commit c8a9165b
Please register or sign in to reply