如何实现壁函数.docVIP

================================================================================================== ? 要准确求解壁面处的流动，需要很细的网格，用壁面函数就是为了避开这一点采用的近似处理。壁面函数在很多书和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朋友参加讨论，更希望提供详细的步骤和有关壁面函数法的子程序! ? ================================================================================================== 我的理解：由于k-e方程要求高雷诺数，所以壁面第一点应布置在粘性支层外，粘性支层外一定范围内速度分布呈对数分布（这是流体理论的研究结果），而壁面函数主要处理的是湍流黏度，k，e，处理这些要用到这个粘性支层厚度，和速度和切应力。如果知道厚度了，就可以根据对数分布求出速度，然后计算其他的。所以壁面函数就是要先求出粘性子层厚度Y+，然后求U+（不要求u和u*，u*是为了无量纲用的，用以简化推导和计算，事实上后边用的都是u+，y+，知道u+和y+就可以干求其他）。壁面函数要求第一点布置在湍流旺盛区（就是确定y+，有推导出来的表达式），而对数分布的成立也是有范围的，所以壁面处网格的划分才是最关键的，一般通过试算搞出来。这就个大概思想，主要的理论你还得看看陶的书，peric的书在怎么推讲得更细致一些。希望大家补充。 ================================================================================================== ?? 壁面函数：???在划分网格时，把第一个内节点布置在对数律成立的范围内，即配置到旺盛湍流区域，11..5~30 y+200~400。???流场计算完后，查看 ：DisplayContours… Contours of /Turbulence…/Wall Yplus???如果y+ 的值大于该范围，应该加密该区域网格，重新计算，再查看y+ ，如果仍不在其范围，继续加密网格。壁面网格加密可采用自适应网格：???AdaptY+/Y*…，Options选项，只选Refine ；Type选Y+；点击Mark，再点击Adapt；及完成网格加密。非平衡壁面函数（Non-Equilibrium Wall Function）主要应用于以下情况：????涉及分离、再附着、冲击等受压力梯度影响的远离平衡的复杂流动Enhanced Wall Treatment要求y+ 4~5。 ================================================================================================== 1．基本思想 对于湍流充分发展的核心流动区域使用标准的K-epsilon模型或其改进模型求解；对壁面分子粘性影响明显的区域，直接用半经验公式将壁面上的物理量与湍流核心区内的求解量联系起来，而不对壁面区内流动求解。也就是将求解的第一个内节点布置在近壁区域的对数律成立的区域里，即湍流充分发展的区域，其内不再配置任何节点，如下图所示：  2.???????第一个内节点动量方程中u?和能量方程中T?与壁面函数值间的关联 在湍流充分发展的对数律层，无量纲速度和温度服从对数律分布。流体力学理论所得到的速度表达式为  =================================================================================================不清楚上式中的那些量是已知的，那些量是待求的。 =================================================================================================