[BUAA LSM]

When solving cavitation, the multiphase flow model uses the vof model, and the source term is added to the phase fraction equation and the pressure Poisson equation, as shown in the figure. However, it is strange that when the source term is not added to the pressure Poisson equation, the pressure field is convergent, but the phase fraction field is wrong. After the source term is added, when cavitation occurs, the source term is not 0 and the pressure field diverges. I'm pretty sure the openfoam code has source term in it.
I used the code written by myself, the pressure and speed coupling used the SIMPLE algorithm, the interface capture used the CICSAM method, the processing of source terms used explicit expression and implicit expression, and the results were divergent.
I tried a lot of methods and I couldn't solve the divergence problem, the cavitation model was the Schnerr and Sauer model