BUGFIX: Corrected coefficient for the virtual mass force of the beta phase in bubbleFoam. Author: Alexander Vakhrushev. Merge: Dominik Christ.

This commit is contained in:
Dominik Christ 2015-04-27 10:35:03 +01:00
commit 11ab9eebd8
2 changed files with 3 additions and 2 deletions

View file

@ -81,3 +81,4 @@ Contents:
Zeljko Tukovic
Niklas Wikstrom
Vanja Skuric
Alexander Vakhrushev

View file

@ -43,7 +43,7 @@ fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime);
UbEqn =
(
(scalar(1) + Cvm*rhob*alpha/rhob)*
(scalar(1) + Cvm*rhoa*alpha/rhob)*
(
fvm::ddt(Ub)
+ fvm::div(phib, Ub, "div(phib,Ub)")
@ -61,7 +61,7 @@ fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime);
// g // Buoyancy term transfered to p-equation
- fvm::Sp(alpha/rhob*dragCoef, Ub)
//+ alpha/rhob*dragCoef*Ua // Explicit drag transfered to p-equation
+ alpha/rhob*(liftCoeff + Cvm*rhob*DDtUa)
+ alpha/rhob*(liftCoeff + Cvm*rhoa*DDtUa)
);
UbEqn.relax();