#259 Bug in bubbleFoam solver: UEqns.H, line 46

[Alexander Vakhrushev]

In bubbleFoam: UEqns.H:46

For the phase beta instead of

  (scalar(1) + Cvm*rhob*alpha/rhob)*

should be

  (scalar(1) + Cvm*rhoa*alpha/rhob)*'

see RuschePhD2002, p.105, Tab. 3.1

[Alexander Vakhrushev]

Line 64 as well:

 + alpha/rhob*(liftCoeff + Cvm*rhob*DDtUa)

should be

 + alpha/rhob*(liftCoeff + Cvm*rhoa*DDtUa)

[Alexander Vakhrushev]

Here comes the patch