Ok, so then if I understand correctly:



  • At a given time step the program first solves the equations of motion for the mixture phase according to the multiphase model (mixture model in this case) along with the turbulence.
  • Next, the cavitation model separately solves for the vapor volume fraction by accounting for evaporation and condensation mass transfer with the Rayleigh-Plesset equation for a specified nucleation site density. The pressure field from the multiphase model is an “input” to this calculation, and the radius of the bubbles is allowed to change according to the Rayleigh-Plesset equation.
  • The vapor volue fraction data is then passed from the cavitation model back to the multi-phase model, which uses a fixed bubble diameter along with the vapor fraction from the cavitation model to determine a number of bubbles. The drag on these bubbles is then calculated and the motion of the secondary phase relative to the mixture phase is obtained.



Is my understanding correct? If so, then there are two separate considerations of bubbles, one within the cavitation model determined by the nucleation site density where the bubble size is allowed to change, and one within the multiphase model where the size of the bubble is not allowed to change and is only used for calculation of the drag between phases. Is this also correct?

Thank you,