Merkle's model (cavitation module) in Code_Saturne v5.0

Questions and remarks about code_saturne usage
Forum rules
Please read the forum usage recommendations before posting.
Post Reply
So_nuc
Posts: 45
Joined: Wed Sep 05, 2018 12:16 pm

Merkle's model (cavitation module) in Code_Saturne v5.0

Post by So_nuc »

Hello,

I have noticed that the source terms from the Merkle's model seem to be computed differently in the sources of the code. The Theory Guide v5.0.0 (page 145) indeed indicates the factor alpha * (1- alpha) instead of alpha for m- or 1- alpha for m+ according to the original model [1]. Moreover a recent PhD Thesis [2] indicates "By default, the model of Code_Saturne with cavitation module is a variation of the Merkle's model" and [3] seems to give the beginning of an explanation in Section 3.3 but it is unclear.

[1] Ding Li and Charles L. Merkle, A UNIFIED FRAMEWORK FOR INCOMPRESSIBLE AND COMPRESSIBLE FLUID FLOWS (2006)

[2] Christophe Leclercq, Simulation numérique du chargement mécanique en paroi généré par les écoulements cavitants, pour application à l’usure par cavitation des pompes centrifuges (2017)

[3] R. CHEBLI a,b , O. COUTIER -D ELGOSAH a , B. AUDEBERT b, Development of algorithm based on the fractional time
step for the simulation of cavitation.

Thank you in advance.

Regards,

So_nuc
Last edited by So_nuc on Tue Nov 05, 2019 2:44 pm, edited 1 time in total.
Yvan Fournier
Posts: 4070
Joined: Mon Feb 20, 2012 3:25 pm

Re: Merkle's model (cavitation module) in Code_Saturne v5.0

Post by Yvan Fournier »

Hello,

I just forwarded a link to you post to colleagues more familiar with this part of the code, so I hope they can provide the answers soon.

Best regards,

Yvan
BLM
Posts: 11
Joined: Thu Oct 11, 2012 9:14 am

Re: Merkle's model (cavitation module) in Code_Saturne v5.0

Post by BLM »

Hello,

You are right, the cavitation source term model implemented in Code_Saturne differs from the original Merkle's model and the reason should be explain in ref [3] : it is because the alpha(1-alpha) factor makes the min/max principle (alpha must be between 0 and 1) is easier to satisfy with this formulation. I think Cdest=50 and Cprod=10000 constants were calibrated in order to take into account this modification in the source term.

I don't know if you have already modify some code_saturne source code but if you want to use the original Merkle model, I think you may just have to modifify the "cavitation_compute_source_term" routine in the cavitation.f90 file.

Best regards,
Benoît
Post Reply