@星星星星晴 老师您好,因为我写的是新的两个标量方程,类似于求解脉动速度的k与epsilon,尝试的时候我先采用的就是类似于k和epsilon的边界条件,用的也是k与epsilon的壁面函数,我想先算一下,有必要的时候我再添加自己的壁面函数。我的计算是用的simpleFoam,所以时间步用的是一秒一步,slover我编译了一个我自认为相匹配的求解器,读取了一下由kt和epsilont求解更新的湍流热扩散系数。数学公式的话我是采用的文献上的写的代码,我觉得老师说的除以0以及负数开根的情况很有可能,我再仔细阅读一下文献。