CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   Siemens (https://www.cfd-online.com/Forums/siemens/)
-   -   Need CAVPRO.f example (https://www.cfd-online.com/Forums/siemens/54237-need-cavpro-f-example.html)

Chris House December 6, 2004 09:31

Need CAVPRO.f example
 
Can someone give me an example of a CAVPRO.f subroutine? I am trying to use the cavitation feature of STAR-CD and need to edit the saturation pressure to be a function rather than a constant.

Chris

4xF December 15, 2004 07:34

Re: Need CAVPRO.f example
 
Just for curiosity: Why do you need to have the saturation pressure as a function? As far as I know, the saturation pressure is a function of temperature and is very little dependent on pressure variations. And the cavitation model implemented in STAR (anybody corrects me if I am wrong here) is isothermal.

Chris House December 15, 2004 14:11

Re: Need CAVPRO.f example
 
I need the saturation pressure as a function because I am working with air and oil in cavitation.

The air/oil system does not behave like a fluid/vapor system. The lower the pressure gets, the more air comes out of the solution. That's why I need to make the saturation pressure as a function.

Can you please help me out with a CAVPRO.f subroutine or the variable names for the Rayleigh model that I need to use with the FORTRAN code so it returns them to STAR correctly?

Any help you could give me would be appreciated.

Chris

4xF December 16, 2004 04:42

Re: Need CAVPRO.f example
 
The CAVPRO subroutine appears to be called on per cell basis (i.e. called for every cell). Thus, you should calculate for every cell PV and send it back to STAR. For example, you can have a variation of PV with the local pressure as:

REAL A,B PARAMETER (A=1.e-5,B=0.01)

PV = A*PR**2 - B*PR + 2370

RETURN END

Hope this helps. If you still do not understand the logic of STAR subroutines programming, I would highly recommend you to attend a training course (2days) on user subroutines. I've learnt a lot doing so.

Take care...


All times are GMT -4. The time now is 13:35.