在OpenFOAM中如何设置变通量边界条件(variably heat flux boundary condition)

[Jhhu]

算例为土柱冻结算例，起始土柱内部温度为6.7℃，土柱顶部施以温度为-6℃的循环液流(the circulating fluid)，使土柱从上到下逐步冻结，下图为示意图，

温度边界条件：
变热通量上边界条件/a variable heat flux upper boundary condition：
qh=hc(TCoolant-TTop), hc=28 W m-2 K-1
其中：TCoolant是循环液流的温度值(-6 ℃)，而TTop是土柱顶部温度，我理解的是每一次时间步长计算后都会更新，所以顶部通量是一个动态更新的过程。
1、之前有看过关于fixed heat flux boundary的帖子，采用fixedGradient边界条件代替，但没找到关于 variably heat flux boundary condition
2、 采用swak4Foam？
希望能给我一些建议，共同探讨下。非常感谢！