Hi @cloner0110,
thanks a lot for having shared this code in a recent version of OF!
This is a question more than a bug:
why rhoPhi is not updated using limitedH*rho1 + (1.0 - limitedH)*rho2;? From the code it looks like rhoPhi for the momentum equation will be the one computed as rhoPhi = tphiAlpha()(rho1 - rho2) + phirho2; in alphaEq.H
Cheers
Hi @cloner0110,
thanks a lot for having shared this code in a recent version of OF!
This is a question more than a bug:
why rhoPhi is not updated using
limitedH*rho1 + (1.0 - limitedH)*rho2;? From the code it looks like rhoPhi for the momentum equation will be the one computed as rhoPhi = tphiAlpha()(rho1 - rho2) + phirho2; in alphaEq.HCheers