Elmer Discussion Forum

Hello,

I don't want to impose a Dirichlet Temp BC when the surface velocity is out of the free surface (v.n >0) so I have this in my top surface BC:

Temp = Variable Coordinate 3
Real MATC "min(273.15 273.15-5.3-0.006*tx)"

Temp Condition = Variable Velocity 1 , Velocity 2 , Velocity 3, Normal Vector 1, Normal Vector 2, Normal Vector 3
Real MATC "Cond(tx)"


and this for Cond:

$ function Cond(tx) {\
dotprod=tx(0)*tx(3)+tx(1)*tx(4)+tx(2)*tx(5);\
_Cond=-dotprod;\
}


However, this does not seem to work. Actually the code stops without any warnings or error message.

Anybody has tried this before?

Thanks

Denis

Denis,
it is odd that you do not get any error message. What I could think of is, that you perhaps dash into an issue if you evaluate this condition before the variables you need (Velocity and Normal Vector) have been assigned with a value, hence leaving the Temp Cond with Zero.

A in the Simulation section might help in that case.

Regards,

Thomas

Hi Thomas,
Thanks for the info. It's running now. Will check results.
Denis