CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   Main CFD Forum (https://www.cfd-online.com/Forums/main/)
-   -   Cavitation: implementation of Schnerr-Sauer model (https://www.cfd-online.com/Forums/main/112550-cavitation-implementation-schnerr-sauer-model.html)

 luftraudi January 30, 2013 17:06

Cavitation: implementation of Schnerr-Sauer model

Hi,

I'm working on the implementation of the Schnerr-Sauer cavitation model in our cfd code, but I had some problems with the correct implementation of the source term in the transport equation. My problems are as follows:

If I initialize my flow field with pure liquid, so the volume fraction of the vapour is zero leading to a vanishing source term (in cavitating regions where the pressure is below the vapour pressure) since the radius R become zero. Thus no cavitation occur. This is quite physical behavior since cavitation is linked to the presence of nuclei.

If I initialize my flow field with nuclei (e.g. vapour volume fraction value of 1.e-5) cavitation can occur in regions where the pressure is below the vapour pressure. Now the problem is, that the condensation source term (in regions where the pressure is above the vapour pressure) is non-vanishing leading to a reduction of the vapour volume fraction. In the worst case the vapour volume fraction become negative, which is not physical.

does anyone has experience with this topic? If so, how did you come up with this aspect? Did you implement a source term which is different from what is presented in your papers or did you modify the numerical procedure (finite volume, Euler implicit, upwind interpolation) for the transport equation for the volume fraction in order to keep the volume fraction between the physical bounds 0 and 1?

With regards,
luftraudi

 RaghavendraRohith November 21, 2013 23:52

Hi Michel

In ANSYS Fluent, as i worked in that the source term enhances phase change from fluid to vapor and vapor to fluid. The additional source term is greater and smaller than the vapour pressure depending on the pressure linked equations.

Please check it in my work.

Thanks
Rohith

 All times are GMT -4. The time now is 21:29.