@李东岳 谢谢李老师。这个问题我当时解决了，忘记回论坛答复了。我的处理办法是，在并行的时候，分节点node_n把所有壁面信息如壁面摩擦速度、壁面摩擦速度传递给node0，由node0统一处理把全部壁面信息传递给所有分节点node_n，然后分节点分别计算壁距离、y+。思路是这样，写在这里供大家参考吧

最近想用下MK低雷诺数湍流模型，但是阻尼函数里面要用到y+，我只知道近壁面第一层网格的y+的求解方法：
$y+=\frac{y}{\upsilon /\sqrt{\tau _{w}/\rho }}$
请问，流体域内部任意网格处的y+也要用这个公式求么？当地y+的摩擦速度怎么取呢？取距离当地网格最近的壁面处的摩擦速度么？
如下图，在整个边界层域的y+都表达出来了，整个域中的摩擦速度用的是当地的还是壁面的呢？

MK，NH，CH此类低雷诺数湍流模型中，其中一个阻尼系数fu是y+的函数: $$f_{\mu }=f\left (\mathrm{Re,}y+ \right )$$而这个阻尼系数是控制方程一部分，要在整个求解域求解，,照这样说y+在整个求解域中都应该是有意义的，不应该仅仅局限于近壁面。但是y+的计算公式中要用到近壁面摩擦速度:$$y+=y{u}_{\tau }/\nu $$，请问远离近壁面的网格单元中的y+的计算应该如何计算呢？是否是用当地离壁距离与法向壁面处的磨擦速度之积除以当地粘度呢？