FLUENT中壁面函数和近壁面模型
湍流流动的近壁处理详解
壁面对湍流有明显影响。
在很靠近壁面的地方，粘性阻尼减少了切向速度脉动，壁面也阻止了法向的速度脉动。
离开壁面稍微远点的地方，由于平均速度梯度的增加，湍动能产生迅速变大，因而湍流增强。
因此近壁的处理明显影响数值模拟的结果，因为壁面是涡量和湍流的主要来源。
实验研究表明，近壁区域可以分为三层，最近壁面的地方被称为粘性底层，流动是层流状态，分子粘性对于动量、热量和质量输运起到决定作用。
外区域成为完全湍流层，湍流起决定作用。
在完全湍流与层流底层之间底区域为混合区域（Blending region），该区域内分子粘性与湍流都起着相当的作用。
近壁区域划分见图4－1。
图4－1，边界层结构第一节，壁面函数与近壁模型近壁处理方法有两类：第一类是不求解层流底层和混合区，采用半经验公式（壁面函数）来求解层流底层与完全湍流之间的区域。
采用壁面函数的方法可以避免改进模型就可以直接模拟壁面存在对湍流的影响。
第二类是改进湍流模型，粘性影响的近壁区域，包括层流底层都可以求解。
对于多数高雷诺数流动问题，采用壁面函数的方法可以节约计算资源。
这是因为在近壁区域，求解的变量变化梯度较大，改进模型的方法计算量比较大。
由于可以减少计算量并具有一定的精度，壁面函数得到了比较多的应用。
对于许多的工程实际流动问题，采用壁面函数处理近壁区域是很好的选择。
如果我们研究的问题是低雷诺数的流动问题，那么采用壁面函数方法处理近壁区域就不合适了，而且壁面函数处理的前提假设条件也不满足。
这就需要一个合适的模型，可以一直求解到壁面。
FLUENT提供了壁面函数和近壁模型两种方法，以便供用户根据自己的计算问题选择。
4.1.1壁面函数FLUENT 提供的壁面函数包括:1，标准壁面函数；2，非平衡壁面函数两类。
标准壁面函数是采用Launder and Spalding [L93]的近壁处理方法。
该方法在很多工程实际流动中有较好的模拟效果。
4.1.1.1 标准壁面函数根据平均速度壁面法则，有：**1ln()U Ey k = 4－1其中，1/41/2*/p pw U C k U μτρ≡，1/41/2*p pC k y y μρμ≡，并且k ＝0.42，是V on Karman 常数；E ＝9.81，是实验常数；p U 是P 点的流体平均速度；p k 是P 点的湍动能；p y 是P 点到壁面的距离；μ是流体的动力粘性系数。
「胡言」壁面函数与近壁面处理
「胡言」壁面函数与近壁面处理Fluent提供了众多的近壁面处理方法，这里简单扒一扒。
1y 的基本概念在临近壁面位置，法向速度非常大的梯度。
在非常小的壁面法向距离内，速度从相对较大的值下降到与壁面速度相同。
因此对于该区域内流场的计算，通常采用两种方式：（1）利用壁面函数法；（2）加密网格，利用壁面模型法。
对于这两类方法的选取，可以通过对于这两类方法的选取，可以通过y 来体现。
如图所示为近壁面位置无量纲速度分布情况。
图中横坐标所表示为无量纲壁面距离y ，纵坐标为无量纲速度u 。
其中：从图中看出，近壁面区域可分为3个区间：•粘性子层（Viscous sublayer region）：黏性子层是一个紧贴壁面的极薄层，在该区域中，粘性力在动量、质量、能量交换过程中起主导地位，湍流剪切应力在该区域可以忽略不计。
该区域中的流动可看成层流流动，在平行于壁面方向上的速度分量沿壁面法向方向呈线性分布。
在该区域中，y <5。
•对数律层（log law region）：对数律层位于近壁区域的最外层，在该区域中，湍流剪切应力占主导地位，粘性力几乎可以忽略，流动为充分发展湍流状态，流速近似成对数分布。
该区域中，y >60。
•过渡层（Buffer layer region）：过渡层位于粘性子层与对数律层之间，该区域中粘性力与湍流剪切应力相当，流动状态极其复杂，难以用模型进行表达。
由于过渡层非常薄，工程中常将其合并到对数律层进行处理。
过渡层中y 在5~60之间。
对于近壁区域求解，主要集中在粘性子层的求解上，主要有两种方式：1、利用壁面模型直接求解粘性子层若想要求解粘性子层，则需要保证y 值小于1（建议接近1）。
由于y 直接影响第一层网格节点位置，因此对于求解粘性子层的情况，需要非常细密的网格，通常要求有10~20层边界层网格。
2、利用壁面函数近似处理。
对于湍流模型，需要选择低雷诺数湍流模型（如k-omega模型）。
FLUENT中文学习资料-建模与数值分析
4.1.1 壁面函数 FLUENT 提供的壁面函数包括:1，标准壁面函数；2，非平衡壁面函数两类。标准壁面函 数是采用 Launder and Spalding [L93]的近壁处理方法。该方法在很多工程实际流动中有较好的 模拟效果。 4.1.1.1 标准壁面函数 根据平均速度壁面法则，有：
U* = 1 ln( Ey* ) k
Gk ≈ τ w
τw ∂U =τw 1/ 4 1/ 2 ∂y kρC µ k p y p
4－5
耗散率不需要求解输运方程，直接用如下公式计算：
εp =
3/ 4 3/ 2 Cµ kp
ky p
4－6
以上所介绍的标准壁面函数是 FLUENT 程序的默认设置。标准壁面函数包含了定常剪切 和局部平衡假设条件，如果壁面有很强的压力梯度，并且很强的非平衡性，则我们可以选择非 平衡壁面函数方法。 4.1.1.2 非平衡壁面函数
* ( y* < yT )
4－3
* ( y* > yT )
其中 P 的计算采用下列公式[L93]
Pr Prt 4－4 − 1 Prt Pr & ′′ 壁面热流； Tp 近邻 其中， k f 是流体导热系数； ρ 是流体密度； c p 是流体定压比热； q π /4 A P= sin(π / 4) k
47
PDF 文件以 "FinePrint pdfFactory Pro" 试用版创建
衡壁面函数则把壁面函数方法推广到有压力梯度和非平衡的流动过程中。但是，如果流动情况 偏离了壁面函数的理想条件，则壁面函数就不合适了。如：高粘度流体流过狭窄的通道，壁面 由渗透的流动，大压力梯度并导致边界层分离的流动，由强体积力的流动，近壁区域三维性很 强的流动问题。 如果要成功解决上述问题， 必须采用改进模型的方法来模拟近壁流动。 FLUENT 提供了双层区模型（Two-Layer Zonal Model ） 。 4.1.1.3 双层区模型
Fluent物理模型概述
Fluent物理模型概述Fluent为各种不可压缩和可压缩、层流和湍流流体流动问题提供了全面的模拟能力。
可以进行稳态或瞬态分析。
在Fluent中，大量传输现象的数学模型(如传热和化学反应)与复杂几何模型的能力相结合。
Fluent应用实例包括:工艺设备层流非牛顿流;叶轮机械与汽车发动机部件的共轭传热；电站锅炉中煤粉燃烧的分析；外部空气动力学;通过压缩机、泵和风扇的流量;以及气泡塔和流化床中的多相流。
为了模拟工业设备和过程中的流体流动和相关的运输现象，本教程提供了各种有用的特性。
包括多孔介质、集总参数(风扇和热交换器)、流向周期性流动和传热、涡流和移动参考系模型。
模型的移动参照系系包括对单个或多个参照系建模的能力。
此外，还提供了一种时间精确的滑动网格方法，用于叶轮机械应用中的多级建模，例如，计算时间平均流场的混合平面模型。
Fluent中另一组非常有用的模型是一组自由面和多相流模型。
这些可用于分析气-液、气-固、液-固和气-液-固流动。
针对这类问题，Fluent提供了(VOF)、混合模型、欧拉模型以及离散相模型(DPM)。
DPM对分散相(粒子、液滴或气泡)进行拉格朗日轨迹计算，包括与连续相耦合。
多相流的例子包括明渠流、喷雾、沉降、分离和空化。
在Fluent模型中，鲁棒性和准确性是湍流模型至关重要的组成部分。
所提供的湍流模型具有广泛的适用性，而且还包括其他物理现象的影响，如浮力和压缩性。
通过使用壁面函数和分区处理模型来求解近壁区域。
各种传热模式可以模拟，包括自然对流、强迫对流、混合对流、多孔介质等。
辐射模型和一些子模型都是可以使用的，还可以计算燃烧。
Fluent的一个特别的优点是它能够使用多种模型来模拟燃烧现象，包括涡流耗散模型和概率密度函数模型。
还有许多其他模型对于反应流应用非常有用，包括煤和液滴燃烧、表面反应和污染物形成模型。
总之，fluent提供了丰富的模型让你来模拟你所感兴趣的问题。
对于所有流动，Fluent求解质量和动量守恒方程。
fluent基本物理模型介绍
基本物理模型本章介绍了FLUENT所提供的基本物理模型以及相关的定义和使用。
基本物理模型概述FLUENT提供了从不可压到可压、层流、湍流等很大范围模拟能力。
在FLUENT中,输运现象的数学模型与所模拟的几何图形的复杂情况是结合在一起的。
FLUENT应用的例子包括层流非牛顿流的模拟，涡轮机和汽车引擎的湍流热传导，锅炉内煤炭粉碎机的燃烧，可压射流，空气动力外流，以及固体火箭发动机的可压化学反应流。
为了与工业应用相结合，FLUENT提供了很多有用的功能。
如多孔介质，块参数（风扇和热交换），周期性流动和热传导，涡流，以及移动坐标系模型。
移动参考系模型可以模拟单一或者多个参考系。
FLUENT还提供了时间精度滑动网格方法以及计算时间平均流动流场的混合平面模型，滑动网格方法在模拟涡轮机多重过程中很有用。
FLUENT中另一个很有用的模型是离散相模型，这个模型何以用于分析喷雾和粒子流。
，多项流模型可以用于预测射流的破散以及大坝塌陷之后流体的运动，气穴现象，沉淀和分离。
湍流模型是FLUENT中很重要的一部分，湍流会影响到其它的物理现象如浮力和可压缩性。
湍流模型提供了很大的应用范围，而不需要对特定的应用做出适当的调节，而且它涵括了其它物理现象的影响，如浮力和可压缩性。
通过使用扩展壁面函数和区域模型，它可以对近壁面的精度问题有很好的考虑。
各种热传导模式可以被模拟，其中包括具有或不具有其它复杂性如变化热传导的，多孔介质的自然的、受迫的以及混合的对流。
模拟相应介质的辐射模型及子模型的设定通常可以将燃烧的复杂性考虑进来。
FLUENT一个最强大的功能就是它可以通过耗散模型或者和概率密度函数模型来模拟燃烧现象。
对于燃烧应用十分有用的其它模型也可以在FLUENT中使用，其中包括碳和液滴的燃烧以及污染形成模型。
连续性和动量方程对于所有的流动，FLUENT都是解质量和动量守恒方程。
对于包括热传导或可压性的流动，需要解能量守恒的附加方程。
FLUENT壁面函数的选择
FLUENT壁面函数的选择壁面函数问题1、无论是标准k—ε模型、RNGk—ε模型，还是Realizable k—ε模型,都是针对充分发展的湍流才有效的，也就是说，这些模型均是高Re数的湍流模型。
它们只能用于求解处于湍流核心区的流动。
而壁面函数是对近壁区的半经验描述，是对某些湍流模型的补充（近壁区对整体流动影响较大和低雷诺数Re的情况），通过壁面函数法和低Re数k—ε模型与标准k—ε模型和RNGk—ε模型配合，成功解决整个整个管道的流动计算问题。
2、在壁面区，流动情况变化很大。
解决这个问题目前有两个途径：一、是不对粘性影响比较明显的区域（粘性底层和过渡层）进行求解，而是用一组半经验的公式（即壁面函数）将壁面上的物理量与湍流核心区内的相应物理量联系起来。
这就是壁面函数法。
在划分网格的时候，不需要在壁面区加密，只需要把第一个节点布置在对数律成立的区域内，即配置在湍流充分发展区域。
如果要用到壁面函数的话，在define---modle--viscous面板里有near wall treatment一项。
可以选择标准壁面函数、不平衡壁面函数等。
二、是采用低Re数的k—ε模型来求解粘性底层和过渡层，此时需要在壁面区划分比较细密的网格，越靠近壁面，网格越细。
当局部湍流的Re数小于150时，就应该使用低Re数的k—ε模型。
总结：相对于低Re数的k—ε模型，壁面函数法计算效率高，工程实用性强。
但当流动分离过大或近壁面流动处于高压之下时，不是很理想。
在划分网格的时候，需要在壁面的位置设置边界层网格，原因也是如此。
为什么要用壁面函数？？就是因为，k-epsilon模型中，k的boundary condition已知，在壁面上为零，而epsilon的boundary condition 在壁面上为一未知的非零量，如此如何来解两方程模型？？？所以，我们就需要壁面函数来确定至少第一内节点上的值，当然也包括壁面上的值。
实际上就是把epsilon方程的boundary condition放到了流体内部。
FLUENT中壁面函数-vs-近壁面模型
FLUENT中壁面函数 vs 近壁面模型在数值模拟中，如何有效处理固体壁面附近的流场一直是一个比较棘手的问题.一个稍复杂一点算例，简单更换一下壁面处理方法对计算结果都有较显著的影响，在缺少实验数据验证和流场涉及多种流动形态时,如何选择行之有效和经济合理的算法是一个艰难的考验，一般需要仔细考察流场与算法机理之间的契合度。
边界层分为层流边界层和湍流边界层，层流边界层为最靠近壁面或者层流流动时的边界层，对于一般湍流流动，两种边界层都有。
按参数分布规律划分时，边界层分为内区和外区，内区分为:粘性底层,Laminar sublayer(y+<5，Amano的三层模型),粘性起主导作用，在粘性支层中与壁面平行的速度与离开壁面的距离成线性关系（陶文铨,《数值传热学》）;过渡层,Buffer region（5〈y+<30），湍流作用与粘性作用共同作用;对数律层，Log—law region（30〈y+），湍流起主导作用，无量纲速度与温度分布服从对数分布律；外区：惯性力主导，上限取决于雷诺数图1 边界层结构（引自中科大Fluent讲稿）FLUENT中有两种方法处理近壁面区域：A.壁面函数法。
不求解粘性影响内部区域（粘性子层及过渡层），使用一种称之为“wallfunction”的半经验方法去计算壁面与充分发展湍流区域之间的粘性影响区域.采用壁面函数法,省去了为壁面的存在而修改湍流模型.Fluent中的standard wallfunctions， scalable wall functions，Non—Equilibrium wall functions和Enhanced wall treatment都属于壁面函数法的模型。
B。
近壁模型法。
修改湍流模型以使其能够求解近壁粘性影响区域，包括粘性底层.此处使用的方法即近壁模型.（近壁模型不需要使用壁面函数，如一些低雷诺数模型，K—W 湍流模型是一种典型的近壁湍流模型）。
Fluent学习笔记之壁面网格与yplus的选择
壁面网格与y+的选择（参考官方手册ANSYS_Fluent_Theory_Guide19.2）大量的实验表明，近壁区域可以大致细分为三层。
在最内层，称为“粘性子层”，流动几乎是层流的，(分子)粘性在动量和热量或质量传递中起主导作用。
在被称为“湍流充分发展层”的外层，湍流起着主要作用。
最后，在粘性子层和湍流充分发展层之间有一个过渡区域，叫做“缓冲层”或“混合区”，在这里分子粘度和湍流的影响是同等重要的其中y+定义为y+≝ρuτy/μ其中uτ定义为uτ=√τw/ρ“壁面函数法”与“近壁模型法”传统的近壁区域建模方法有两种。
在第一种方法中，不解决粘度影响的内部区域(粘性子层和缓冲层)。
相反，被称为“壁面函数”的半经验公式被用来连接壁面和全紊流区域之间的粘滞区域。
使用壁面函数避免了需要修改湍流模型来考虑壁面的存在。
在另一种方法中，对湍流模型进行了修改，使粘滞区通过网格一直到壁面(包括粘滞子层)进行求解。
为了便于讨论，我们将其称为“近壁建模”方法。
这两种方法如图4.14所示。
图4.14: ANSYS Fluent中的近壁面处理所有壁面函数(scalable wall function除外)的主要缺点是在壁面法向网格细化的情况下，计算结果会恶化。
Y+小于15会逐渐导致壁面剪应力和壁面传热的误差增大。
但这已经是几年前的工业标准，ANSYS Fluent已经能够提供更先进的壁面方程，允许一致的网格细化没有恶化的结果。
对所有基于ω方程的湍流模型都采用这种与y+无关的公式。
对于基于ε-方程的模型，Menteri-Lechner和增强壁面处理(Enhanced Wall Treatment, EWT)具有相同的作用。
对y+不敏感的壁面处理也是Spalart-Allmaras模型的默认值，它允许你不用考虑近壁面的y+值而运行此模型。
只有当边界层的整体分辨率足够时，才能得到高质量的边界层数值结果。
这个要求实际上比实现特定的y+值更重要。
matlab涡流方程
matlab涡流方程01—湍流模型standard, RNG, and realizable k - ε模型，这三种模型的形式都很相似，都有k和ε的输运方程。
这些模型的主要区别如下：湍流粘度的计算方法控制K和湍流扩散的湍流普朗特数方程中的生成项和消耗项standard k - ε模型是一个半经验模型，模型方程的推导依赖于现象和经验。
standard模型是基于湍流动能K及其耗散率ε的输运方程的模型。
K的模型传输方程是从精确方程推导出来的，而ε的模型传输方程是通过物理推理得到的，与数学上的精确方程相似性很小。
RNG k-ε模型是使用一种称为重整化群理论的统计方法推导出来的。
它在形式上类似于standard k-ε模型,但包括以下改进: RNG模型在其ε方程中增加了一项，提高了高速流动的准确性。
RNG模型考虑了涡流对湍流的影响，提高了旋涡流动的精度。
RNG理论提供了湍流普朗特数的解析公式，而standard k-ε使用用户指定的常数值。
虽然标准k-ε模型是高雷诺数模型,RNG微分公式理论提供了一种从分析中获得有效的粘度,考虑低雷诺数的影响，然而，这个特性的有效使用取决于对近壁区域的适当处理。
这些特性使得在更广泛的流动中RNG k-ε模型比standard k-ε更准确和可靠的。
Realizable k-ε模型和standard k-ε模型主要有两点不同：Realizable k-ε模型包含另一种湍流粘度公式。
修正后的耗散率ε输运方程是从准确的涡度波动的输运方程推导而来。
02—Fluent中的壁面函数和近壁模型近壁区域可大致细分为三层。
在最内层，称为“粘性底层”，流动几乎是层流的，粘度在动量和传热传质中起主导作用。
外层被称为完全湍流层，湍流起着主要作用。
在粘性底层和完全湍流层之间存在一个过渡区域，分子粘度和湍流的影响同样重要。
壁面流动 fluent学习笔记
近壁区流动及fluent求解对策一、边界层1、边界层（boundary layer）是高雷诺数绕流中紧贴物面的粘性力不可忽略的流动薄层，又称流动边界层、附面层。
在这部分区域中，沿着固壁面切向速度由固壁处的0速度发展到接近来流的速度，一般定义为在边界处的流速达到来流流速的99%。
在这部分区域中，由于厚度很小，故速度急剧变化，速度梯度很大，流体的粘性效应也主要体现在这一区域中。
2、边界层有层流、湍流、混合流3、边界层分离边界层流动从物体表面脱离的现象。
二维边界层分离有两种情况，一是发生在光滑物面上，另一是发生在物面有尖角或其他外形中断或不连续处。
光滑物面上发生分离的原因在于，边界层内的流体因克服粘性阻力而不断损失动量，当遇到下游压力变大（即存在逆压梯度）时，更需要将动能转变为压力能，以便克服前方压力而运动，当物面法向速度梯度在某位置上小到零时，表示一部分流体速度已为零，成为“死水”，边界层流动无法沿物面发展，只能从物面脱离，该位置称为分离点。
分离后的边界层在下游形成较大的旋涡区；但也可能在下游某处又回附到物面上，形成局部回流区或气泡。
尖点处发生边界层分离的原因在于附近的外流流速很大，压强很小，因而向下游必有很大的逆压梯度，在其作用下，边界层即从尖点处发生分离。
三维边界层的分离比较复杂，是正在深入研究的课题。
边界层分离导致绕流物体压差阻力增大、飞机机翼升力减小、流体机械效率降低、螺旋桨性能下降等，一般希望避免或尽量推迟分离的发生；但有时也可利用分离，如小展弦比尖前缘机翼的前缘分离涡可导致很强的涡升力。
二、flunet求解壁面区流动对于有固体壁面的充分发展的湍流流动，沿壁面法线的不同距离上可将流动划分为壁面区和核心区。
核心区是完全湍流区，壁面区又分为：粘性底层、过渡层、对数律层。
粘性底层中粘性力占主导，对数律层粘性力影响不明显，流动处于充分发展的湍流状态。
在flunet中，无论是标准k—ε模型、RNGk—ε模型，还是Realizable k—ε模型,都是针对充分发展的湍流才有效的，也就是说，这些模型均是高Re数的湍流模型。
CFD中的壁面函数是怎么回事？
CFD中的壁面函数是怎么回事？近壁区域流体变量具有较大的梯度，因此对近壁区域的精细模拟基本上可以决定壁面流（wall-bounded flow）的结果是否成功。
大部分高雷诺数湍流模型，如kepsilon湍流模型、雷诺应力模型、LES模型，仅仅适用于离开壁面一定距离的湍流区域。
在高雷诺数区域，层流粘度相对于湍流粘度可以忽略不计。
然而在壁面区域附近，比如粘性之层中，雷诺数很低。
因此在这一部分区域需要考虑层流粘度的影响。
Spalart-Allmaras和komega湍流模型在壁面网格足够细化的情况下，可以适用于全雷诺数范围的流动。
对于kEpsilon模型，其中适用于粘性支层的kepsilon模型称之为低Re-kepsilon模型。
低雷诺数kepsilon湍流模型需要在粘性支层以及缓冲层内布置20个左右网格节点。
三维情况下，低雷诺数kepsilon湍流模型需要耗费大量的计算资源在近壁区。
如果采用传统的高Re数ke模型来计算，对于壁面附近区域，需要采用壁面函数法。
壁面函数的精髓在于，在那些流体变量梯度较大的近壁区域，不需要进行求解。
然而，壁面函数不适用于求解区域各处的雷诺数均较低的情况。
在使用壁面函数的时候，壁面的第一层网格需要布置在湍流边界层内，同时需要布置在足够旺盛的湍流区内。
例如log区是很好的选择。
如下图所示。
上图中，我们的网格节点Up便布置在log区。
这样便节省了粘性支层过密网格带来的资源消耗。
同时，我们还需要确保网格节点不能布置的太远。
那么在CFD计算前，如何能粗略的获取第一层网格高度的值呢？1. 首先我们需要管壁的摩擦因数，其大体可以这样计算：•内流：Cf=0.079Re^-0.25•外流：Cf=0.058Re^-0.22. 在获取Cf之后，计算壁面剪切力τ_w=0.5 Cf ρ U^23. 计算摩擦速度u_τ= (τ_w/ρ)^0.54. 第一层网格高度为h=(y+ μ)/(ρ u_τ)上面只是简单几何的情况下，第一层网格高度的计算公式。
大涡模拟壁面函数wener
大涡模拟壁面函数wener大涡模拟（LES）是一种高保真度流动模拟方法，适用于计算高雷诺数下的湍流流动，常用于工程应用中的气体和液体流动。
而壁面函数也是LES模拟中一个重要的问题，因为在现实流动中，壁面附近的细节变化对整个流场的影响非常重要。
这篇文章将介绍LES中涉及到的壁面函数wener，以及其特点和应用。
1、壁面函数概述壁面函数是一种通过数学公式来描述流体流经固体壁面附近流动的方法。
在LES计算中，由于涡模拟的栅格大小要比实际物理尺度小得多，因此需要壁面函数来描述栅格内的流动，提高模拟精度。
壁面函数通常包括平均速度和摩擦应力模型，能够模拟近壁区域的速度和摩擦应力分布。
2、wener壁面函数wener壁面函数是一种常用的LES壁面函数，由Popovac和Hanjalic （WEN model)在1989年提出，并在1999年被Weller、Tabor、Jasak等人进行了修正，称为WEN-TAB。
wener壁面函数可以分为两部分：内部和外部。
内部可以用logarithmic函数表示，外部可以用指数函数表示。
wener壁面函数适用于气体和液体流动，可以通过修正来适应不同的流动条件。
3、 wener壁面函数的特点wener壁面函数的主要特点可以总结为以下几点：（1）适用范围广：wener壁面函数适用于各种流动条件，包括气体和液体流动。
（2）高精度：wener壁面函数使用logarithmic和指数函数来表示壁面内外的流动，能够高精度地描述流场中的速度和摩擦应力分布。
（3）易于实现：wener壁面函数的公式比较简单，容易实现，并且可以通过修正来适应不同的流动条件。
（4）计算效率高：wener壁面函数的计算效率高，对模拟结果的影响也比较小，能够提高模拟速度和精度。
4、wener壁面函数在工程应用中的应用wener壁面函数在工程应用中广泛使用，例如汽车、航空、机械等领域的气流模拟，以及海洋、河流等领域的水流模拟。
fluent湍流模型
k–ω 湍流模型
k–ω 湍流模型得到广泛特点：
模型方程不包括在壁面上没有定义的项，例如不需要壁面函数可以 在壁面积分
对于有压力梯度的大范围边界层流动是精确稳定的
FLUENT 提供k–ω 模型下的两个子模型 标准k–ω (SKW) 模型
在航天和涡轮机械领域得到最广泛的应用 几个k–ω子模型选项：压缩效果，转錑，剪切流修正. 剪切应力输运k–ω (SSTKW) 模型(Menter, 1994) SST k–ω 模型使用混合函数从壁面附近的标准k–ω 模型逐渐过渡到边
解决低雷诺数下的differential viscosity模型 由解析方法得到的 Prandtl / Schmidt数的代数公式 旋流修正 对更复杂的剪切流来说比SKE 表现更好，比如剪切流、旋涡和分离 流
© 2006 ANSYS, Inc. All rights reserved.
在FLUENT中可用的湍流模型
基于RANS的模 型
1－方程模型 Spalart-Allmaras 2－方程模型
标准 k–ε RNG k–ε realizable k–ε 标准 k–ω SST k–ω 雷诺德应力模型 分离涡模拟 大涡模拟
增加 每个计算迭代步
消耗
© 2006 ANSYS, Inc. All rights reserved.
Rij uiuj
(Reynolds 应力张量)
Reynolds 应力是由附加的平均过程引起的，因此为了封闭控制方程 组，必须对Reynolds应力建模
© 2006 ANSYS, Inc. All rights reserved.
6-7
ANSYS, Inc. Proprietary
2.1.5 壁面和近壁区网格处理原则[共2页]
![2.1.5 壁面和近壁区网格处理原则[共2页]](https://img.taocdn.com/s3/m/7ddf76fae2bd960591c6770f.png)
32 精通CFD工程仿真与案例实战—
—FLUENT GAMBIT ICEM CFD Tecplot
图2-11 翼形的非结构三角形网格 图2-12 非结构四面体网格
图2-13 三角形与四边形混合网格 图2-14 多面体网格 2.1.5 壁面和近壁区网格处理原则
通常使用的湍流模型（如k-ε模型）是针对充分发展的湍流才有效的，它们只在高Re （雷诺数）的湍流模拟中运用。
但是，近壁区附近Re 较低，湍流发展并不充分，湍流的脉动影响不如分子粘性的影响大，该区域不能使用高Re 数的湍流模型，必须采用特殊的处理方法，而且在网格划分上也要进行特殊的处理。
本节介绍的壁面函数法可以与k-ε模型等配合，解决近壁区下的流动计算问题。
1．壁面边界层与壁面函数
对于有壁面的流动，当主流为充分发展湍流时，根据离壁面法线的距离的不同，可将流动划
分为壁面区（或称内区、近壁区）和核心区（或称外区）。
核心区是完全湍流区，为充分发展的湍流。
在壁面区，由于有壁面的影响，流动与核心区不同。
壁面区可分为3个子层，如图2-15所示。
● 粘性底层。
● 过渡层。
● 对数律层。
粘性底层是一个紧贴壁面的极薄层，在动量、热量
和质量的交换过程中粘性力起主要作用，而湍流切应力可以忽略，因此流动几乎可以看成层流流动，且在平行图2-15 近壁区3个子层的划分 l。
fluent壁面函数
fluent壁面函数
Fluent壁面函数是3D流体力学软件Fluent的一种材料流体模型的计算方法，用于解决复杂流动的数学模型。
Fluent壁面函数可以模拟从一个表面发出的流量，以及它如何随着它在壁面上移动而改变。
Fluent壁面函数使用壁面平衡来计算壁面流体的性质。
它们通过简单的流体力学原理来描述壁面上流体的行为。
在壁面平衡中，以特定的壁面上的单位面积为基础，在该面积上求和所有流体运动控制件，以z轴为准确地计算壁面的准确运动。
Fluent壁面函数可用于描述复杂的表面形状，如管道、箱子、板条和船等。
可以模拟在特定表面上的流动，以及表面上流体的改变情况。
它们还可以用来预测整体流体行为，如流量和速度。
Fluent壁面函数可以帮助研究者更好地理解流体是如何在某一表面上传播的，从而更好地设计流体系统。
它们还可以用来确定流体系统中存在的压力、温度和流速。
- 1 -。
fluent 对流换热边界
fluent 对流换热边界Fluent对流换热边界流体力学领域中，流动与换热问题一直是研究的热点之一。
在工程实践中，准确预测流体流动和换热行为对于设计高效的换热设备和系统至关重要。
Fluent作为一种流体力学分析软件，具备强大的模拟和计算能力，能够对流换热边界进行准确的模拟和分析。
在Fluent中，对流换热边界是指通过流体流动进行热量传递的边界条件。
换热过程是指热量从热源传递到冷源的过程，而对流则是指流体通过流动方式进行热量传递的过程。
对流换热边界的准确模拟可以帮助工程师理解和优化热交换设备的性能。
对流换热边界的模拟需要考虑多种因素，如流体的性质、流动的速度、壁面的温度等。
Fluent提供了多种模型和方法来模拟不同类型的对流换热边界。
其中，最常用的模型是壁面函数模型和湍流模型。
壁面函数模型是一种简化模型，适用于壁面边界层很薄的情况。
在这种模型中，通过壁面函数来描述壁面附近的流体流动和换热行为。
这种模型计算速度快，适用于大规模的工程计算。
湍流模型则是考虑了流动中的湍流效应。
湍流是流体流动中的一种不规则运动状态，会对流动和换热行为产生重要影响。
Fluent提供了多种湍流模型，如k-ε模型、k-ω模型和RNG k-ε模型等。
这些模型基于不同的假设和方程，适用于不同的流动情况。
在进行对流换热边界的模拟时，需要设置边界条件和模型参数。
边界条件包括壁面温度、壁面热流和壁面热传导系数等。
这些边界条件可以根据实际情况进行设定，以准确模拟实际的换热过程。
除了边界条件，模型参数的选择也对模拟结果的准确性和稳定性有重要影响。
模型参数的选择需要考虑流动特性、壁面特性和计算资源等因素。
在Fluent中，可以通过网格独立性分析和参数敏感性分析来优化模型参数的选择，以获得更准确的模拟结果。
使用Fluent进行对流换热边界的模拟可以帮助工程师更好地理解和优化换热设备的性能。
通过模拟分析，可以获得壁面温度分布、热传导率、热流分布等重要参数，为工程设计和优化提供依据。
7 Fluent软件的使用(6)
涡粘模型
• 在动量扩散方面的作用，涡粘性类似于分子粘性。 • 涡粘性不是流体材料的属性，而是湍流的特性。 • 涡粘模型在CFD的湍流模拟中应用最为广泛。
• 涡粘模型的局限性:
– 各向同性假设，真实的许多流动中具有明显的各向异性特征 （流线弯曲严重的流动、强旋转流动等等）。
– 没有考虑流动旋转对于雷诺应力的影响。
outer layer
inner layer
近壁区网格要求
• 标准的壁面函数:
– k–ε模型
– RSM模型
• 非平衡的壁面函数:
– k–ε模型
– RSM模型
• 加强的壁面函数
– k–ε模型
– RSM模型 – k–ω模型（默认均采用）
近壁网格单元尺寸的估算
• 对于平板，湍流摩擦力系数的幂指数定律为
RANS 方程和封闭问题
• 物理量的时间平均定义如下 （所有湍流都是非定常的，但非定 常并不都是湍流引起的）
• 物理量的瞬态值定义为平均值与脉动值之和
• 对N-S方程进行时间平均，可以得到 RANS 方程。
雷诺应力, Rij
雷诺应力张量
• Rij 是二阶对称张量，该项来源于动量方程中对流加速 度项的时间平均。
标准的 k–ε模型
• 湍能耗散率ε 是指将动能通过粘性转化为内能。
• 机理上，湍动能将通过粘性效应耗散。产生和耗散的 不平衡将会引起湍动能的增加或者减少。
标准的 k–ε模型
• 标准模型中ε 方程完全是半经验模拟的
• 耗散率是湍动能和湍流特征长度的函数
• 涡粘性系数可以表示为
标准的 k–ε模型
• SKE 是在工业界应用最广的湍流模型。
– 通过log分布律可以获得速度剖面和壁面剪切应力，用来设置近壁网格 单元中剪切应力的边界值。 – 平衡假设被用来设置k， ε或ω的物面边界条件。 – 非平衡的壁面函数方法用来改善下列流动的计算结果：大的逆压梯度、 分离、再附和驻点。 – 在能量方程和组分方程中也建立了类似的log分布律。 – 优点：采用壁面函数后，可以在壁面附近采用粗网格，减少计算成本。
胡言｜Fluent湍流模型II
胡言｜Fluent湍流模型II本文描述RANS系列湍流模型对网格的需求。
注：本文内容译自Fluent UserGuide.1RANS模型壁面处理建议对所有可用的模型(Spalart-Allmaras、ε-equation和ω-equation)使用y 非敏感壁面处理。
它提供了对Y 值最低敏感度的一致的壁面剪切应力和壁面传热预测。
当使用壁面函数（Wall Functions）时，应避免在近壁面区域使用过于细密的网格，一般建议在整个计算区域中保证y >30。
通常不推荐在计算过程中使用壁面函数，因为壁面函数不允许对近壁面网格进行系统的细化。
对于低雷诺数到中雷诺数(Re~1e4 - 1e6)的流动，使用壁面函数会引起较大的计算误差。
此时若使用壁面函数，可使用Scalable Wall Functions（Scalable壁面函数）。
该模型可以避免网格限制，允许在精细网格上使用。
2RANS模型的网格需求网格生成对模型的精度有很大的影响。
在生成高质量CFD计算网格时，需要考虑很多因素。
从湍流模型的角度来看，沿剪切层法向方向至少存在不低于10层计算网格。
若网格过于粗糙，则模型无法达到其标定的性能。
特别是对于位置未知的自由剪切流动，划分计算网格更加麻烦。
然而应该清楚，网格分辨率不足会导致模型性能下降。
对于有界壁面流动，强烈建议沿壁面法向采用结构网格（边界层网格）。
结构网格部分应覆盖整个边界层并超出边界层厚度，避免限制边界层的增长。
对于壁面边界层的高级湍流模型，如Spalart-Allmaras模型和SST模型，只有当边界层内包含至少10个或更多的结构网格(六面体或棱柱网格)时，才能为其模型提供改进的结果。
此外，应确保棱柱层完全覆盖壁面边界层。
注意，这些不是这些模型的具体要求，而是壁面边界层模拟的一般要求。
ε-based和ω-based模型提供的y 不敏感壁面处理方案，使模型对壁面网格的y 值相对不敏感。
一般来说，更重要的是要保证边界层被足够的网格填充，然后才能达到一定的y 标准。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
FLUENT中壁面函数和近壁面模型
技术邻作者：Jessica_4643
文章所包含相关领域及技术点：壁面函数、近壁面模型、fluent
在数值模拟中，如何有效处理固体壁面附近的流场一直是一个比较棘手的问题。
一个稍复杂一点算例，简单更换一下壁面处理方法对计算结果都有较显著的影响，在缺少实验数据验证和流场涉及多种流动形态时，如何选择行之有效和经济合理的算法是一个艰难的考验，一般需要仔细考察流场与算法机理之间的契合度。
边界层分为层流边界层和湍流边界层，层流边界层为最靠近壁面或者层流流动时的边界层，对于一般湍流流动，两种边界层都有。
按参数分布规律划分时，边界层分为内区和外区，内区分为：
粘性底层，Laminar sublayer(y+<5，Amano的三层模型)，粘性起主导作用，在粘性支层中与壁面平行的速度与离开壁面的距离成线性关系（陶文铨，《数值传热学》）；
过渡层，Buffer region(5<y+<30)，湍流作用与粘性作用共同作用；
对数律层，Log-law region（30<y+），湍流起主导作用，无量纲速度与温度分布服从对数分布律；
外区：惯性力主导，上限取决于雷诺数
图1边界层结构（引自中科大Fluent讲稿）
FLUENT中有两种方法处理近壁面区域：
A.壁面函数法。
不求解粘性影响内部区域（粘性子层及过渡层），使用一种称之为“wall function”的半经验方法去计算壁面与充分发展湍流区域之间的粘性影响区域。
采用壁面函数法，省去了为壁面的存在而修改湍流模型。
Fluent中的standard wa ll functions，scalable wall functions,Non-Equilibrium wall functions和En hanced wall treatment都属于壁面函数法的模型。
B.近壁模型法。
修改湍流模型以使其能够求解近壁粘性影响区域，包括粘性底层。
此处使用的方法即近壁模型。
（近壁模型不需要使用壁面函数，如一些低雷诺数模型，K-W湍流模型是一种典型的近壁湍流模型）。
所有壁面函数（除scalable壁面函数外）的最主要缺点在于：沿壁面法向细化网格时，会导致数值结果恶化。
当y+小于15时，将会在壁面剪切力及热传递方面逐渐导致产生无界错误。
然而这是若干年前的工业标准，如今ANSYS FLUENT采取了措施提供了更高级的壁面格式，以允许网格细化而不产生结果恶化。
而y+无关的格式是默认的基于w方程的低湍流模型，其采用网格求解的方式计算近壁面粘性区域。
对于基于epsil on方程的模型，增强壁面函数（EWT）提供了相同的功能。
这一选项同样是SA模型所默认的，该选项允许用户使其模型与近壁面y+求解无关。
只有当所有的边界层求解都达到要求了才可能获得高质量的壁面边界层数值计算结果。
这一要求比单纯的几个Y+值达到要求更重要。
使用近壁模型法时，覆盖边界层的最小网格数量在10层左右，最好能达到20层。
还有一点需要注意的是，提高边界层求解常常可以取得稳健的数值计算结果，因为只需要细化壁面法向方向网格。
对于非结构网格，建议划分10~20层棱柱层网格以提高壁面边界层的预测精度。
棱柱层厚度应当被设计为保证有15层或更多网格节点。
另外，棱柱层大于边界层厚度是必要的，否则棱柱层会限制边界层的增长。
这可以在获得计算结果后，通过查看边界层中心的最大湍流粘度，该值提供了边界层的厚度（最大值的两倍位置即边界层的边）。
一些建议：（1）对于epsilon方程，使用enhanced壁面函数。
（2）若壁面函数有助于epsilon方程，则可以使用scalable壁面函数。
（3）对于基于w方程的模型，使用默认的增强壁面函数。
（4）SA模型，使用增强壁面处理。
文章后续内容可到技术邻平台继续观看，进入技术邻后在页面最上面搜索文章的原标题即可观看。