CFD中文网

是的，都是绝对坐标系的

算例为土柱冻结算例，起始土柱内部温度为6.7℃，土柱顶部施以温度为-6℃的循环液流(the circulating fluid)，使土柱从上到下逐步冻结，下图为示意图，
捕获.JPG
温度边界条件：
变热通量上边界条件/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？
希望能给我一些建议，共同探讨下。非常感谢！

@东岳 第一次发帖求助，就能得到东岳老师的回复，真的非常感谢！
（1）是根据laplacianFOAM稍微改了下的求解器，但是如果是探针probe的话是可以运算出结果的。那cellSource不运行function的话，我需要怎么改下求解器吗
（2）还想请教下，算例边界条件设置的问题，左图是压力，右图是温度。
左图压力上下边界条件均为zero flux，右侧问题上下右边界条件为zero conductive flux。对于zero flux和zero conductive flux，边界条件是不是应设置为symmetryPlane，还是zeroGradient。我计算出的结果差别很大，不知道问题出在哪里，所以想确定下是不是我边界条件设置的不对。
case.JPG