Code_Saturne User's Forum

Hello, Generally the mass flow rate is computed as \dot{m}=\rho \cdot \vec{V} \cdot \vec{n} \cdot A . In the documentation is mention that the surface vector of the faces are oriented outwards. Also, I think that a simple example like a fluid flow in a tube with an inlet and outlet can answer this q...

Hello, I'm not a Mac user...but I have searched about this problem and several Mac users report this PyQt/Qt problem using dialog box in different softwares. My first recommendation is try to install CS 4.0 and check if works. May be the developers can give you same additional information about this...

Hello,

I think that this error come from Qt when managed the dialog box. Is this error only produced with the velocity initialization box?

Could you give more information about what OS are you using? Did you try to compile another version of CS (e.g. 4.0).

Regards,

Luciano

Hello,

I attach the case, the listing and the source files. I noted that you changed the Cdest and Cprod parameters. I hope that you find it useful and can solve this problem

Regards,

Luciano

Hello, 1) You can define multiple zones to be initialized with different values (e.g. Different velocities). This is can be done in the GUI under the Volume Condition Tab -> Definition of volume regions. Here you can add multiple zones and specify a selection criterion for the cells (e.g. coordinate...

Hello, The cavitation number and vapor pressure are absolute pressures. So, if you want to use a P_{ref}=P_0>0 , the outlet gauge pressure P_g for a given cavitation number is given by, P_{abs}=0.5 \rho V^2+P_{v}=P_g+P_{ref} , wherewith P_g=0.5 \rho V^2+P_{v}-P_{ref} . In the figure that you attach ...

Hello, As Benoît mention the problem was in the pressure. In your case if you set a P_{ref}=0 and the outlet pressure according to the cavitation number \sigma you will get the cavitation of the fluid at the throat. I have run you case with a dt=5e-6 [s] and the results are reasonable according to t...

Hello, Inside of the boundary faces loop, you can access to the coordinates of the centers of the boundary faces thought the pointer cdgfbo ! Boundary faces of group '3' call getfbr('3', nlelt, lstelt) do ilelt = 1, nlelt ifac = lstelt(ilelt) ! y-coordinate Face center ycoord = cdgfbo(2, ifac) enddo...

Hello Amir, I was looking the cavitation.f90 file and the theory documentation, because I never work with cavitation. In the theory documentation of the CS 4.0.0 the vaporization source terms are m^+ = -\dfrac{C_{prod} \rho_l \min \left( p-p_V,0 \right)}{0.5\rho_lu_\infty^2t_\infty}, m^- = -\dfrac{C...

Hello amir,

Could you also post the case and the mesh files?

Regards,

Luciano