Inserting beta_grad for vtk-output in InOut.cpp.
Eliminating setRHS_PBC and adding unCoupledBetaPBC(). It was discovered the reason why the rising bubble was incorrect. There was an error when multiplying the inverse matrix invMLumped by the term of M beta. In fact, the matrix should be invMrhoLumped. With this modification, the profile was recovered and the code could be validated. However, the value of beta is variable concerning the average pressure calculation. The 2D version got excelent accordance with beta = 0.97. Now, the term M betaFlow was placed into unCoupledBetaPBC(), which will be the method to be used for beta, but without PBC. Method's name will be changed later, however.
Commiting version with the beta validated. mainRisingBubblePBC like this is OK to recover the base profile. The best value for beta is 0.9625. However, some way to optimize the calculation of beta should be found.
Included different mains for risingBubble tests. The PBC implementation is validated for these cases. However, the pressure is still with convergence problems in the two last cases, viz. beta and PBC, and I'm simulating during more time to check if they'll break. Not sure about what's generating the convergence problem. The Poiseuille mains were tested but they are not good. Check...