exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

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

Browse source
This commit is contained in:
Dominik Christ 2015-04-27 14:46:52 +01:00
commit fb637562c9
2 changed files with 3 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

1
ListOfContributors
View file

@ -74,3 +74,4 @@
Zeljko Tukovic Zeljko Tukovic
Niklas Wikstrom Niklas Wikstrom
Vanja Skuric Vanja Skuric
Alexander Vakhrushev

4
applications/solvers/multiphase/bubbleFoam/UEqns.H
View file

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