Where is the solve of UEqn in twoPhaseEulerFoam?
 

Hi Foamers,
In the manual "Derivation, Modelling and Solution of the Conditionally Averaged Two-Phase Flow Equations" by Weller et al. they said that momentum equations are constructed and solved. But in the solver code of twoPhaseEulerFoam, I just see UaEqn.relax(), not solve(). So could anyone explain to me how about this issue? Thank you!:confused:

I think the momentum equation is not predicted and is only corrected in pressure loop which is conservative. (see pEqn.H)

The momentum equation is defined (excluding some terms, which are transferred to the pressure equation), and the velocity is predicted with a single Jacobi step.

You could predict the velocity with a direct solution, including the transferred term by means of the reconstruct operator, however this might slow convergence down.

could you please explain how to use a reconstruct operator?



I am using twoPhaseEulerFoam(5.0) and have read some literatures about the solving process. As you said, the solver only predict the velocity rather than solve it... how could I get the value of velocity?



thanks! It would be a lot of help.