壁面函数的求解步骤

我是大神

大家好：我对壁面函数的求解步骤不是很理解，我已经推出来了所有壁面函数的方程。接下来怎么求解呢。

有个人这么说的：要准确求解壁面处的流动， 需要很细的网格， 用壁面函数就是为了避开这一点 采用的近似处理。壁面函数在很多书和PAPER里都提到过， 但不同模型和不同的人相差很远，而且没有完整的步骤。 我在编程中用到高雷诺数两方程模型，碰到了壁面函数的问题： 1)由初始的速度U，按对数律计算U+； 2)由U+计算出Y+； 3)判断Y+>11.5，第一内点P位于旺盛湍流区，符合对数律，求P点U，K，E以及壁面W点的U，K，E 4)若Y+<11.5，第一内点P位于粘性支层，按U+=Y+计算。 这是我的理解，但更详细的细节，我还没弄清？比如P、W点的U、V、K、E的具体计算表达式。 以上谈到的是规则域的壁面函数法处理，对于贴体坐标转换的壁面函数法处理起来更复杂，因为与壁面平行的速度才满足对数律。 希望CFD朋友参加讨论，更希望提供详细的步骤和有关壁面函数法的子程序!


但是他第一步 初始的速度u是怎么得出来的？ 比如一个长方形，左边进口右边出口，上下是壁面，第一个节点的速度不能初始吧？ 并且即使可以初始速度 壁面函数的uplus和k也有关 k都不知道怎么能出来uplus呢？  谢谢

zdong_hn

我想这 个也想了很久，对于壁面函数，关键在摩擦速度的确定，当由对数律或线性率确定第一层网格节点速度时，实际给出的是混合边界条件。可以用迭代的方式求得摩擦速度后确定第一层网格节点速度。
对于湍流参数，应该好办吧，因为速度分布有了。
不知道真实情况是不是这样子的。

且听风吟R

同样是这个问题