自定义边界条件的设置问题

[摇滚小鱼1007]

本人OF和C++小白一枚。求助一个边界条件设置问题。

方程是一个泥沙含沙量的标量输运方程，在壁面需要设置一个泥沙被冲起的边界条件，给了一个法向梯度边界

P表示冲刷量，是一个以近壁面摩阻流速值U_tau作为因变量的公式。P=f(U_tau)
水流在壁面用的是spalding壁面函数。U_tau的求解在spalding壁面函数类中以一个protected函数的方式给出。

我想请教一下大佬们：
（1）这样的边界条件需要自定义库来做吗？难不难实现？
（2）如果想引入spalding类中的保护型函数来给出U_tau，能否通过继承来实现？或者有没有其他好的方法？
（3）还有就是，有没有类似边界的，已有的程序或者其他编程资料可以参考？

[李东岳]

不难实现。是一个非常简单的边界条件。难度是一个OKS课的教学内容。我给你信心，你自己试试吧...很难手把手把这个写出来给你 参考下固定梯度边界条件，然后把你的$-P/V$包含进去就好

[Samuel-Tu]

@摇滚小鱼1007 对于第2个问题，定梯度的边界条件和spalding边界有继承关系吗。。我有把类似计算utau的代码放入自定义的代码中，自己计算utau来用。
另外你的这个边界条件中v代表什么，C是泥沙在水中的浓度吗?这个边界是怎么推导出来的啊，有无论文参考啊，像是梯度等于泥沙冲刷量除以切向速度，为什么梯度是怎么算的，我一直想不通。我是积雪冲刷问题，和这个很像，但是我发现论文里面都没提冲刷面的边界条件的计算方法，我自己算雪颗粒输运方程，发现雪颗粒浓度一直算不对，可能就是这个边界没设置对。

[李东岳]

@Samuel-Tu 在 自定义边界条件的设置问题 中说：

雪颗粒输运方程

这个考虑颗粒粒径了么？
颗粒传输速度给的是流体速度还是自行封闭的？

[Samuel-Tu]

@东岳 雪颗粒速度是流体速度，相当于单向耦合。。

[Samuel-Tu]

@东岳 方程里也没有考虑粒径。。

[李东岳]

https://www.cfd-china.com/topic/1324 想到了这个贴

[Samuel-Tu]

@东岳 我做的就是这个。。

[摇滚小鱼1007]

@Samuel-Tu sorry很久没上来看了。这个边界就是简单的泥沙冲刷通量边界，v是粘度。P是冲刷量，在我们泥沙界，这个量是比较难确定的，一般是前人老爷子做试验给出的经验公式。这个量到底用什么公式，里面的参数怎么取，对结果都会造成决定性影响。我不清楚你的雪颗粒，是否指的是浓度？