Turbulent Prandtl Number as function of Density
 

Hi everybody!

Thank you for your support.
I have problem to use and expression with Density variable on the Turbulent Prandtl Number during a STEADY and TRANSIENT simulation.

CFX-Solver crashes immediately with this error:

----------------------------------
Error in subroutine GETCORE :
There is circularity in recursive calls to GETVAR. A variable depends upon
itself! Oper = What = Static Enthalpy Where = Default Domain When = TIME-0
Action = RETURN. Stopping
GETVAR originally called by subroutine DEF_DIFTRBC

+--------------------------------------------------------------------+
| Writing crash recovery file |
+--------------------------------------------------------------------+

+--------------------------------------------------------------------+
| ERROR #001100279 has occurred in subroutine ErrAction.
| Message:
| Stopped in routine GETCORE
| |
+--------------------------------------------------------------------+

I tried to create a new variable with Density and to set this with the "special function" of " Update Loop = TRANS_LOOP" (during TRANS simulations) but it does not work the same.

I think it's necessary to create a Junction Box Routine (to call every time-step) to save a additional variable of density at the preceding iteration...

Can anyone help me to create a Junction Box routine?
Ansys CFX manual was not useful for me to understand well how to create a FORTRAN routines...

Thank you for your time.
Best regards,
Gaspare