How to set boundary condition like q=k*exp(T)? DEFINE_HEAT_FLUX or global iteration?
 

Hi, all!
The problem is following.
I have rather complicated model of heat transfer.
So I need boundary condition for heat flux like q=k*exp(T) (really - some exps).
Usual conditions allow me only linear connection between temperature in cell and temperature on face and q=kT^4 - for emission.
How can I use my own law for heat flux?
Macro DEFINE_HEAT_FLUX contains only constant coefficients, but may be I don't know smth about it.
Another way - to fix initial temperature field, then calculate boundary heat flux and solve with this flux. But I should repeat it manually - how can I automate this action?Thanks a lot!