ENH: forces FO - porosity calc updated for mesh changes

This commit is contained in:
andy
2014-04-08 11:59:41 +01:00
committed by Andrew Heather
parent 792ee881af
commit 557e7ebda8

View File

@ -922,7 +922,8 @@ void Foam::forces::calcForcesMoment()
forAllConstIter(HashTable<const porosityModel*>, models, iter) forAllConstIter(HashTable<const porosityModel*>, models, iter)
{ {
const porosityModel& pm = *iter(); // non-const access required if mesh is changing
porosityModel& pm = const_cast<porosityModel&>(*iter());
vectorField fPTot(pm.force(U, rho, mu)); vectorField fPTot(pm.force(U, rho, mu));