浅话边界条件与初始条件
边界条件
在说边界条件之前，先谈谈初值问题和边值问题。
初值和边值问题：
对一般的微分方程，求其定解，必须引入条件，这个条件大概分两类---初始条件和边界条件，如果方程要求未知量y(x)及其导数y′(x)在自变量的同一点x=x0取给定的值，即y(x0 )=y0，y′(x0)= y0′，则这种条件就称为初始条件，由方程和初始条件构成的问题就称为初值问题;
而在许多实际问题中，往往要求微分方程的解在在某个给定的区间a ≤ x ≤b的端点满足一定的条件，如y(a) = A , y(b) = B 则给出的在端点(边界点)的值的条件，称为边界条件，微分方程和边界条件构成数学模型就称为边值问题。
三类边界条件:
边值问题中的边界条件的形式多种多样，在端点处大体上可以写成这样的形式，Ay+By'=C，若B=0，A≠0,则称为第一类边界条件或狄里克莱(Dirichlet)条件;B≠0,A=0，称为第二类边界条件或诺依曼(Neumann)条件;A≠0,B≠0，则称为第三类边界条件或洛平 (Robin)条件。
总体来说，
第一类边界条件：给出未知函数在边界上的数值;
第二类边界条件：给出未知函数在边界外法线的方向导数;
第三类边界条件：给出未知函数在边界上的函数值和外法向导数的线性组合。
对应于comsol，只有两种边界条件：
Dirichlet boundary(第一类边界条件)—在端点，待求变量的值被指定。
Neumann boundary(第二类边界条件)—待求变量边界外法线的方向导数被指定。
再补充点初始条件：
初始条件，是指过程发生的初始状态，也就是未知函数及其对时间的各阶偏导数在初始时刻t=0的值.在有限元中，好多初始条件要预先给定的。不同的场方程对应不同的初始条件。
总之，为了确定泛定方程的解，就必须提供足够的初始条件和边界条件.
边界条件与初始条件是控制方程有确定解的前提。边界条件是在求解区域的边界上所求解的变量或其导数随时间和地点的变化规律。对于任何问题，都需要给定边界条件。初始条件是所研究对象在过程开始时刻各个求解变量的空间分布情况，对于瞬态问题，必须给定初始条件，稳态问题，则不用给定。
对于边界条件与初始条件的处理，直接影响计算结果的精度。
在CFD模拟时，基本边界条件有：
1流动进口边界
包括速度进口边界，压力进口边界，质量进口边界(可压流动)。
在使用流动进口边界时，需要涉及到某些流动参数，如绝对压力，湍动能及耗散率，这些参数要做特殊考虑。关于参考压力，在流场数值计算中，压力总是按相对值表示的，实际求解的压力并不是绝对值，而是相对于进口压力而言的。在有些情况下，可以通过设定进口压力为0，求解其他点的压力。还有时，为了减小数字截断误差，往往故意抬高或降低参考压力场的值，可使其余各处的计算压力场与整体数值计算的量级相吻合。
2流动出口边界
一般选在离几何扰动足够远的地方来施加。在这样的位置，流动是充分发展的，沿流动方向没有变化。该边界只有在进入计算域的流动是以进口边界条件给定时才使用，而且在只有一个出口的计算域中使用。
3给定压力边界
在流动分布的详细信息未知，但边界的压力值已知的情况下，使用恒压边界条件。应用该边界条件的典型问题包括：物体外部绕流，自由表面流，自然通风及燃烧等浮力驱动流和有多个出口的内部流动。
4壁面边界
对于湍流计算，在近壁面区，要用壁面函数来做处理。壁面函数法的前提是：流动平行于壁面而且只在垂直壁面的方向变化，流动方向不存在压力梯度，壁面处不存在化学作用，雷诺数足够大。若这些条件之一不满足，则使用壁面函数法的预测精度会大大降低，甚至完全不可用。
5对称边界
在对称边界上，垂直边界的速度取为0，而其他物理量的值在该边界内外是相等的，即计算域外紧邻边界的节点的值等于对应的计算域内紧邻边界的节点的值。
6周期性(循环)边界
流出循环边界出口的所有流动变量的通量等于进入循环边界的对应变量的通量。
使用边界条件要注意的地方：
①选择合理的边界条件，边界条件的4种组合方式：只有壁面，壁面、进口和至少一个出口，壁面、进口和至少一个恒压边界，壁面和恒压边界。
②为得到准确的结果，出口边界位置的选取要位于最后一个屏障物后10倍于屏障高度或更远的位置。
③近壁面网格。要保证壁面函数法有效，就需要使离壁面最近的一内节点位于湍流的对数律层中，即y+必须大于11.63，这就相当于给最靠近壁面的网格到壁面的距离△yp设定了一个下限。
初始条件
在瞬态问题中，除了要在计算开始之前初始化相关的数据外，不需要其他特殊处理。给定初始条件时要注意的是：要针对所有计算变量，给定整个计算域内各单元的初始条件;初始条件一定是物理上合理的，要靠经验或实测结果。
来源：元计算官网
Fluent技巧 边界条件 定义边界条件概述 边界条件包括流动变量和热变量在边界处的值。它是FLUENT分析得很关键的一部分，设定边界条件必须小心谨慎。 边界条件的分类：进出口边界条件：压力、速度、质量进口、进风口、进气扇、压力出口、压力远场边界条件、质量出口、通风口、排气扇；壁面、repeating, and pole boundaries:壁面，对称，周期，轴；内部单元区域：流体、固体(多孔是一种流动区域类型) ；内部表面边界：风扇、散热器、多孔跳跃、壁面、内部。(内部表面边界条件定义在单元表面，这意味着它们没有有限厚度，并提供了流场性质的每一步的变化。这些边界条件用来补充描述排气扇、细孔薄膜以及散热器的物理模型。内部表面区域的内部类型不需要你输入任何东西。) 下面一节将详细介绍上面所叙述边界条件，并详细介绍了它们的设定方法以及设定的具体合适条件。周期性边界条件在本章中介绍，模拟完全发展的周期性流动将在周期性流动和热传导一章中介绍。 使用边界条件面板 边界条件(Figure 1)对于特定边界允许你改变边界条件区域类型，并且打开其他的面板以设定每一区域的边界条件参数 菜单：Define/Boundary Conditions... Figure 1: 边界条件面板 改变边界区域类型 设定任何边界条件之前，必须检查所有边界区域的区域类型，如有必要就作适当的修改。比方说：如果你的网格是压力入口，但是你想要使用速度入口，你就要把压力入口改为速度入口之后再设定。 改变类型的步骤如下：: 1.在区域下拉列表中选定所要修改的区域 2.在类型列表中选择正确的区域类型 3.当问题提示菜单出现时，点击确认 确认改变之后，区域类型将会改变，名字也将自动改变 (如果初始名字时缺省的请参阅边界条件区域名字一节),设定区域边界条件的面板也将自动打开。 ！注意：这个方法不能用于改变周期性类型，因为该边界类型已经存在了附加限制。创建边界条件一节解释了如何创建和分开周期性区域。需要注意的是，只能在图一中每一个类别中改变边界类型(注意：双边区域表面是分离的不同单元区域.) Figure 1: 区域类型的分类列表 设定边界条件 在FLUENT中，边界条件和区域有关而与个别表面或者单元无关。如果要结合具有相同边界条件的两个或更多区域请参阅合并区域一节。 设定每一特定区域的边界条件，请遵循下面的步骤： 1.在边界条件区域的下拉列表中选择区域。 2. 点击Set...按钮。或者，1.在区域下拉列表中选择区域。 2.在类型列表中点击所要选择的类型。或者在区域列表中双击所需区域.，选择边界条件区域将会打开，并且你可以指定适当的边界条件
湍流边界条件设置 在流场的入口、出口和远场边界上，用户需要定义流场的湍流参数。在FLUENT 中可以使用的湍流模型有很多种。在使用各种湍流模型时，哪些变量需要设定，哪些不需要设定以及如何给定这些变量的具体数值，都是经常困扰用户的问题。本小节只讨论在边界上设置均匀湍流参数的方法，湍流参数在边界上不是均匀分布的情况可以用型函数和UDF（用户自定义函数）来定义，具体方法请参见相关章节的叙述。 在 大多数情况下，湍流是在入口后面一段距离经过转捩形成的，因此在边界上设置均匀湍流条件是一种可以接受的选择。特别是在不知道湍流参量的分布规律时，在边 界上采用均匀湍流条件可以简化模型的设置。在设置边界条件时，首先应该定性地对流动进行分析，以便边界条件的设置不违背物理规律。违背物理规律的参数设置 往往导致错误的计算结果，甚至使计算发散而无法进行下去。 在Turbulence Specification Method （湍流定义方法）下拉列表中，可以简单地用一个常数来定义湍流参数，即通过给定湍流强度、湍流粘度比、水力直径或湍流特征长在边界上的值来定义流场边界上的湍流。下面具体讨论这些湍流参数的含义，以保证在设置模型时不出现违背流动规律的错误设置： （1）湍流强度（Turbulence Intensity） 湍流强度I的定义为：I=Sqrt(u’*u’+v’*v’+w’*w’)/u_avg (8-1) 上式中u',v' 和w' 是速度脉动量，u_avg是平均速度。 湍流强度小于1％时，可以认为湍流强度是比较低的，而在湍流强度大于10％时，则可以认为湍流强度是比较高的。在来流为层流时，湍流强度可以用绕流物体的几何特征粗略地估算出来。比如在模拟风洞试验的计算中，自由流的湍流强度可以用风洞的特征长度估计出来。
2011-8-30蓝色流体|流体专业论坛专注流体 - Pow… 标题: [fluent相关]湍流边界条件参数的设置 作者: ifluid 时间: 2009-4-14 15:02 标题: 湍流边界条件参数的设置 在流场的入口、出口和远场边界上，用户需要定义流场的湍流参数。在FLUENT 中可以使用的湍流模型 有很多种。在使用各种湍流模型时，哪些变量需要设定，哪些不需要设定以及如何给定这些变量的具 体数值，都是经常困扰用户的问题。本小节只讨论在边界上设置均匀湍流参数的方法，湍流参数在边 界上不是均匀分布的情况可以用型函数和UDF（用户自定义函数）来定义，具体方法请参见相关章节的 叙述。 在大多数情况下，湍流是在入口后面一段距离经过转捩形成的，因此在边界上设置均匀湍流条件是一种可以接受的选择。特别是在不知道湍流参量的分布规律时，在边界上采用均匀湍流条件可以简 化模型的设置。在设置边界条件时，首先应该定性地对流动进行分析，以便边界条件的设置不违背物 理规律。违背物理规律的参数设置往往导致错误的计算结果，甚至使计算发散而无法进行下去。在 Turbulence Specification Method （湍流定义方法）下拉列表中，可以简单地用一个常数来定义湍 流参数，即通过给定湍流强度、湍流粘度比、水力直径或湍流特征长在边界上的值来定义流场边界上 的湍流。下面具体讨论这些湍流参数的含义，以保证在设置模型时不出现违背流动规律的错误设置： （1）湍流强度（Turbulence Intensity） 湍流强度I的定义为： I=Sqrt(u’*u’+v’*v’+w’*w’)/u_avg 上式中u',v' 和w' 是速度脉动量，u_av g是平均速度。 湍流强度小于1％时，可以认为湍流强度是比较低的，而在湍流强度大于10％时，则可以认为湍流强 度是比较高的。在来流为层流时，湍流强度可以用绕流物体的几何特征粗略地估算出来。比如在模拟 风洞试验的计算中，自由流的湍流强度可以用风洞的特征长度估计出来。在现代的低湍流度风洞中， 自由流的湍流强度通常低于0.05%。 内流问题进口处的湍流强度取决于上游流动状态。如果上游是没有充分发展的未受扰流动，则进口处可以使用低湍流强度。如果上游是充分发展的湍流，则进口处湍流强度可以达到几个百分点。如 果管道中的流动是充分发展的湍流，则湍流强度可以用公式(8-2)计算得到，这个公式是从管流经验公 式得到的： I=u’/u_avg=0.16*Re_DH^-0.125 其中Re_DH是Hy draulic Diameter（水力直径）的意思，即式(8-2)中的雷诺数是以水力直径为特 征长度求出的。 （2）湍流的长度尺度与水力直径 湍流能量主要集中在大涡结构中，而湍流长度尺度l则是与大涡结构相关的物理量。在充分发展的管流中，因为漩涡尺度不可能大于管道直径，所以l 是受到管道尺寸制约的几何量。湍流长度尺度l 与管道物理尺寸L关系可以表示为： l = 0.07L 式中的比例因子0.07是充分发展管流中混合长的最大值，而L则是管道直径。在管道截面不是圆形 时，L可以取为管道的水力直径。
浅话边界条件与初始条件 边界条件 在说边界条件之前，先谈谈初值问题和边值问题。 初值和边值问题： 对一般的微分方程，求其定解，必须引入条件，这个条件大概分两类---初始条件和边界条件，如果方程要求未知量y(x)及其导数y′(x)在自变量的同一点x=x0取给定的值，即y(x0 )=y0，y′(x0)= y0′，则这种条件就称为初始条件，由方程和初始条件构成的问题就称为初值问题; 而在许多实际问题中，往往要求微分方程的解在在某个给定的区间a ≤ x≤b 的端点满足一定的条件，如y(a) = A , y(b) = B则给出的在端点(边界点)的值的条件，称为边界条件，微分方程和边界条件构成数学模型就称为边值问题。 三类边界条件: 边值问题中的边界条件的形式多种多样，在端点处大体上可以写成这样的形式，Ay+By'=C，若B=0，A≠0,则称为第一类边界条件或狄里克莱(Dirichlet)条件;B≠0,A=0，称为第二类边界条件或诺依曼(Neumann)条件;A≠0,B≠0，则称为第三类边界条件或洛平(Robin)条件。 总体来说， 第一类边界条件：给出未知函数在边界上的数值; 第二类边界条件：给出未知函数在边界外法线的方向导数; 第三类边界条件：给出未知函数在边界上的函数值和外法向导数的线性组合。 对应于comsol，只有两种边界条件： Dirichlet boundary(第一类边界条件)—在端点，待求变量的值被指定。
Neumann boundary(第二类边界条件)—待求变量边界外法线的方向导数被指定。 再补充点初始条件： 初始条件，是指过程发生的初始状态，也就是未知函数及其对时间的各阶偏导数在初始时刻t=0的值.在有限元中，好多初始条件要预先给定的。不同的场方程对应不同的初始条件。 总之，为了确定泛定方程的解，就必须提供足够的初始条件和边界条件.边界条件与初始条件是控制方程有确定解的前提。边界条件是在求解区域的边界上所求解的变量或其导数随时间和地点的变化规律。对于任何问题，都需要给定边界条件。初始条件是所研究对象在过程开始时刻各个求解变量的空间分布情况，对于瞬态问题，必须给定初始条件，稳态问题，则不用给定。 对于边界条件与初始条件的处理，直接影响计算结果的精度。 在CFD模拟时，基本边界条件有： 1流动进口边界 包括速度进口边界，压力进口边界，质量进口边界(可压流动)。 在使用流动进口边界时，需要涉及到某些流动参数，如绝对压力，湍动能及耗散率，这些参数要做特殊考虑。关于参考压力，在流场数值计算中，压力总是按相对值表示的，实际求解的压力并不是绝对值，而是相对于进口压力而言的。 在有些情况下，可以通过设定进口压力为0，求解其他点的压力。还有时，为了减小数字截断误差，往往故意抬高或降低参考压力场的值，可使其余各处的计算压力场与整体数值计算的量级相吻合。 2流动出口边界 一般选在离几何扰动足够远的地方来施加。在这样的位置，流动是充分发展的，沿流动方向没有变化。该边界只有在进入计算域的流动是以进口边界条件给定时才使用，而且在只有一个出口的计算域中使用。
定解条件和定解问题 含有未知函数的偏导数的方程叫偏微分方程，常微分方程可以看成是特殊的偏微分方程。方程的分数是1的称为方程式，个数多于1的叫做方程组。方程（组）中出现的未知函数的最高阶偏导数的阶数称为方程（组）的阶数。如果方程（组）中的项关于未知函数及其各阶偏导数的整体来讲是线性的，就称方程（组）为线性的，否则就称为非线性的。非线性又分为半线性、拟线性和完全非线性。 一、定解条件 给定一个常微分方程，有通解和特解的概念。通解只要求满足方程，即满足某种物理定律，而不能完全确定一个物理状态。特解除了要求满足方程还要满足给定的外加（特殊）条件。对偏微分方程也是如此，换句话说，只有偏微分方程还不足以确定一个物理量随空间和时间的变化规律，因为在特定情况下这个物理量还与它的初始状态和它在边界受到的约束有关。描述初始时刻的物理状态和边界的约束情况，在数学上分别称为初始条件（或初值条件）和边界条件（或边值条件），他们统称为定解条件。 初始条件：能够用来说明某一具体物理现象初始状态的条件，即描述物理过程初始状态的数学条件。 边界条件：能够用来说明某一具体物理现象边界上的约束情况的条件，即描述物理过程边界状态的数学条件。 定解条件：初始条件和边界条件的统称。 非稳态问题：定解条件包括初始条件和边界条件。
稳态问题：定解条件为边界条件。 1、弦振动方程 ( 2(,),0,0tt xx u a u f x t x l t -=<<>) 初始条件是指初始时刻（0t =）弦的位移和速度。若以()x ?， ()x ψ分别表示弦上任意点x 的初始位移和初始速度，则初始条件为： 边界条件是指弦在两端点的约束情况，一般有三种类型。 （1）第一类边界条件（狄利克雷(Dirichlet ）边界条件）:已知端点()x a a o a l ===或处弦的位移是()a g t ，则边界条件为： (0,)(0,)u t g t = 或 (,)(,)u l t g l t = 当0()0()0l g t g t ≡≡或时，表示在该点处弦是固定的。 （2）第二类边界条件（诺伊曼（Neumann ）边界条件）：已知端点0x x l ==或处弦所受的垂直于弦线的外力0()g t 或()l g t ，则边界条件为： 0(0,)()x Tu t g t -= 或 (,)()x l Tu l x g t = 当00()0l g g t ≡≡或时，表示弦在端点0x x l ==或处自由滑动。 （3）第三类边界条件（混合边界条件或罗宾（Robin ）边界条件：已知端点处弦的位移和所受的垂直于弦线的外力的和： 000(0,)(0,)g (t),0,x Tu t k u t k -+=> 或 (,)(,)(),0x l l l Tu l t k u l t g t k +=>， (,0)(),0(,0)(), t u x x x l u x x ?ψ=?<
FLUENT边界条件（2）—湍流设置 （fluent教材—fluent入门与进阶教程于勇第九章） Fluent：湍流指定方法（Turbulence Specification Method） 2009-09-16 20:50 使用Fluent时，对于velocity inlet边界，涉及到湍流指定方法（Turbulence Specification Method），其中一项是Intensity and Hydraulic Diameter （强度和水利直径），本文对其进行论述。 其下参数共两项， （1）是Turbulence Intensity，确定方法如下： I=0.16/Re_DH^0.125 (1) 其中Re_DH是Hydraulic Diameter（水力直径）的意思，即式(1)中的雷诺数是以水力直径为特征长度求出的。 雷诺数 Re_DH=u×DH/υ(2) u为流速，DH为水利直径，υ为运动粘度。 水利直径见（2）。 （2）水利直径 水力直径是水力半径的二倍，水力半径是总流过流断面面积与湿周之比。 水力半径 R＝A/X (3) 其中，A为截面积（管子的截面积）＝流量/流速 X为湿周(字面理解水流过各种形状管子外圈湿一周的周长) 例如：方形管的水利半径 R=ab/2(a+b) 水利直径 DH=2×R (4) 举例如下： 如果水流速度u=10m/s，圆形管路直径2cm，水的运动粘度为1×10-6 m2/s。 则 DH=2×3.14*r^2/(2*3.14*r)=2*3.14*0.01^2/(3.14*0.02)=0.01 r为圆管半径 Re_DH=u×DH/υ=10*0.02/10e-6=20000 I=0.16/Re_DH^0.125=0.16/20000^0.125=0.0463971424017634≈5%
边界条件 定义边界条件概述 边界条件包括流动变量和热变量在边界处的值。它是FLUENT分析得很关键的一部分，设定边界条件必须小心谨慎。 边界条件的分类：进出口边界条件：压力、速度、质量进口、进风口、进气扇、压力出口、压力远场边界条件、质量出口、通风口、排气扇；壁面、repeating, and pole boundaries:壁面，对称，周期，轴；内部单元区域：流体、固体(多孔是一种流动区域类型) ；内部表面边界：风扇、散热器、多孔跳跃、壁面、内部。(内部表面边界条件定义在单元表面，这意味着它们没有有限厚度，并提供了流场性质的每一步的变化。这些边界条件用来补充描述排气扇、细孔薄膜以及散热器的物理模型。内部表面区域的内部类型不需要你输入任何东西。) 下面一节将详细介绍上面所叙述边界条件，并详细介绍了它们的设定方法以及设定的具体合适条件。周期性边界条件在本章中介绍，模拟完全发展的周期性流动将在周期性流动和热传导一章中介绍。 使用边界条件面板 边界条件(Figure 1)对于特定边界允许你改变边界条件区域类型，并且打开其他的面板以设定每一区域的边界条件参数 菜单：Define/Boundary Conditions... Figure 1: 边界条件面板 改变边界区域类型 设定任何边界条件之前，必须检查所有边界区域的区域类型，如有必要就作适当的修改。比方说：如果你的网格是压力入口，但是你想要使用速度入口，你就要把压力入口改为速度入口之后再设定。 改变类型的步骤如下：: 1.在区域下拉列表中选定所要修改的区域
2.在类型列表中选择正确的区域类型 3.当问题提示菜单出现时，点击确认 确认改变之后，区域类型将会改变，名字也将自动改变(如果初始名字时缺省的请参阅边界条件区域名字一节),设定区域边界条件的面板也将自动打开。 ！注意：这个方法不能用于改变周期性类型，因为该边界类型已经存在了附加限制。创建边界条件一节解释了如何创建和分开周期性区域。需要注意的是，只能在图一中每一个类别中改变边界类型(注意：双边区域表面是分离的不同单元区域.) Figure 1: 区域类型的分类列表 设定边界条件 在FLUENT中，边界条件和区域有关而与个别表面或者单元无关。如果要结合具有相同边界条件的两个或更多区域请参阅合并区域一节。 设定每一特定区域的边界条件，请遵循下面的步骤： 1.在边界条件区域的下拉列表中选择区域。 2. 点击Set...按钮。或者，1.在区域下拉列表中选择区域。 2.在类型列表中点击所要选择的类型。或者在区域列表中双击所需区域.，选择边界条件区域将会打开，并且你可以指定适当的边界条件 在图像显示方面选择边界区域 在边界条件中不论你合适需要选择区域，你都能用鼠标在图形窗口选择适当的区域。如果你是第一次设定问题这一功能尤其有用，如果你有两个或者更多的具有相同类型的区域而且你想要确定区域的标号（也就是画出哪一区域是哪个）这一功能也很有用。要使用该功能请按下述步骤做： 1.用网格显示面板显示网格。 2.用鼠标指针（默认是鼠标右键——参阅控制鼠标键函数以改变鼠标键的功能）在图形窗口中点击边界区域。在图形显示中选择的区域将会自动被选入在边界条件面板中的区域列表中，它的名字和编号也会自动在控制窗口中显示改变边界条件名字 每一边界的名字是它的类型加标号数（比如pressure-inlet-7）。在某些情况下你可能想要对边界区域分配更多的描述名。如果你有两个压力入口区域，比方说，你可能想重名名它们
1 湍流边界条件设置 2 在流场的入口、出口和远场边界上，用户需要定义流场的湍流参数。在 3 FLUENT 中可以使用的湍流模型有很多种。在使用各种湍流模型时，哪些变量需4 要设定，哪些不需要设定以及如何给定这些变量的具体数值，都是经常困扰用5 户的问题。本小节只讨论在边界上设置均匀湍流参数的方法，湍流参数在边界6 上不是均匀分布的情况可以用型函数和UDF（用户自定义函数）来定义，具体方7 法请参见相关章节的叙述。 8 在 9 大多数情况下，湍流是在入口后面一段距离经过转捩形成的，因此在边界10 上设置均匀湍流条件是一种可以接受的选择。特别是在不知道湍流参量的分布11 规律时，在边 12 界上采用均匀湍流条件可以简化模型的设置。在设置边界条件时，首先应13 该定性地对流动进行分析，以便边界条件的设置不违背物理规律。违背物理规14 律的参数设置 15 往往导致错误的计算结果，甚至使计算发散而无法进行下去。 16 在Turbulence Specification Method （湍流定义方法）下拉列表中，可17 以简单地用一个常数来定义湍流参数，即通过给定湍流强度、湍流粘度比、水18 力直径或湍流特征长在边界上的值来定义流场边界上的湍流。下面具体讨论这19 些湍流参数的含义，以保证在设置模型时不出现违背流动规律的错误设置：20 （1）湍流强度（Turbulence Intensity） 21 湍流强度I的定义为：22 I=Sqrt(u’*u’+v’*v’+w’*w’)/u_avg
24 (8-1) 25 上式中u',v' 和w' 是速度脉动量，u_avg是平均速度。 26 湍流强度小于1％时，可以认为湍流强度是比较低的，而在湍流强度大于27 10％时，则可以认为湍流强度是比较高的。在来流为层流时，湍流强度可以用28 绕流物体的几何特征粗略地估算出来。比如在模拟风洞试验的计算中，自由流29 的湍流强度可以用风洞的特征长度估计出来。在现代的低湍流度风洞中，自由30 流的湍流强度通常低于0.05%。 31 内流问题进口处的湍流强度取决于上游流动状态。如果上游是没有充分发32 展的未受扰流动，则进口处可以使用低湍流强度。如果上游是充分发展的湍流，33 则进口处湍流强度可以达到几个百分点。如果管道中的流动是充分发展的湍流，34 则湍流强度可以用公式(8-2)计算得到，这个公式是从管流经验公式得到的：35 I=u’/u_avg=0.16*Re_DH^-0.125 36 (8-2) 37 其中Re_DH是Hydraulic Diameter（水力直径）的意思，即式(8-2)中的雷38 诺数是以水力直径为特征长度求出的。 39 （2）湍流的长度尺度与水力直径 40 湍流能量主要集中在大涡结构中，而湍流长度尺度l则是与大涡结构相关41 的物理量。在充分发展的管流中，因为漩涡尺度不可能大于管道直径，所以l 是42 受到管道尺寸制约的几何量。湍流长度尺度l 与管道物理尺寸L关系可以表示43 为： 44 45 l =
1、湍流强度 定义：速度波动的均方根与平均速度的比值 小于1%为低湍流强度，高于10%为高湍流强度。 计算公式： I=0.16*(re)^(-1/8) 式中：I—湍流强度，re—雷诺数 2、湍流尺度及水力直径 湍流尺度（turbulence length)：a physical quantity related to the size of the large eddies that contain the energy in turbulent flows。 通常计算方式： l=0.07L L为特征尺度，可认为是水力直径，因数0.07是基于充分发展的湍流管流中的混合长度的最大值。 湍流参数的选取： （1）充分发展的内部流动，选取湍流强度(intensity)和水力直径(hydraulic diameter) （2）导流叶片流动、穿孔板等流动，选取强度(intensity)和长度尺度(length scale)。 （3）四周为壁面引起湍流边界层的流动，选取强度（intensity)和长度尺度(length scale)，使用边界层厚度，特征长度等于0.4倍边界层，输入此值到turbulence length scale中。 3、湍动能（Kinetic energy) 湍流模型中最常见的物理量（k）。利用湍流强度估算湍动能： k=3/2*(u*I)^2 其中：u—平均速度，I—湍流强度 4、湍流耗散率(turbulent disspipation rate)
湍流耗散率即传说中的ε。通常利用k和湍流尺度l估算ε计算公式为： cu通常取0.09，k为湍动能，l为湍流尺度 5、比耗散率ω 计算公式为： ω=k^0.5/(l*c^0.25) 式中：k为湍动能，l为湍流尺度，c为经验常数，常取0.09
在入口、出口或远场边界流入流域的流动，FLUENT 需要指定输运标量的值。本节描述了对于特定模型需要哪些量，并且该如何指定它们。也为确定流入边界值最为合适的方法提供了指导方针。 使用轮廓指定湍流参量 在入口处要准确的描述边界层和完全发展的湍流流动，你应该通过实验数据和经验公式创建边界轮廓文件来完美的设定湍流量。如果你有轮廓的分析描述而不是数据点，你也可以用这个分析描述来创建边界轮廓文件，或者创建用户自定义函数来提供入口边界的信息。一旦你创建了轮廓函数，你就可以使用如下的方法： ● Spalart-Allmaras 模型：在湍流指定方法下拉菜单中指定湍流粘性比，并在在湍流粘性 比之后的下拉菜单中选择适当的轮廓名。通过将m_t/m 和密度与分子粘性的适当结合， FLUENT 为修改后的湍流粘性计算边界值。 ● k-e 模型：在湍流指定方法下拉菜单中选择K 和Epsilon 并在湍动能（Turb. Kinetic Energy ）和湍流扩散速度（Turb. Dissipation Rate ）之后的下拉菜单中选择适当的轮廓名。 ● 雷诺应力模型：在湍流指定方法下拉菜单中选择K 和Epsilon 并在湍动能（Turb. Kinetic Energy ）和湍流扩散速度（Turb. Dissipation Rate ）之后的下拉菜单中选择适当的轮廓名。在湍流指定方法下拉菜单中选择雷诺应力部分，并在每一个单独的雷诺应力部分之后的下拉菜单中选择适当的轮廓名。 湍流量的统一说明 在某些情况下流动流入开始时，将边界处的所有湍流量指定为统一值是适当的。比如说，在进入管道的流体，远场边界，甚至完全发展的管流中，湍流量的精确轮廓是未知的。 在大多数湍流流动中，湍流的更高层次产生于边界层而不是流动边界进入流域的地方，因此这就导致了计算结果对流入边界值相对来说不敏感。然而必须注意的是要保证边界值不是非物理边界。非物理边界会导致你的解不准确或者不收敛。对于外部流来说这一特点尤其突出，如果自由流的有效粘性系数具有非物理性的大值，边界层就会找不到了。 你可以在使用轮廓指定湍流量一节中描述的湍流指定方法，来输入同一数值取代轮廓。你也可以选择用更为方便的量来指定湍流量，如湍流强度，湍流粘性比，水力直径以及湍流特征尺度，下面将会对这些内容作一详细叙述。 湍流强度I 定义为相对于平均速度u_avg 的脉动速度u^'的均方根。 小于或等于1%的湍流强度通常被认为低强度湍流，大于10%被认为是高强度湍流。从外界，测量数据的入口边界，你可以很好的估计湍流强度。例如：如果你模拟风洞试验，自由流的湍流强度通常可以从风洞指标中得到。在现代低湍流风洞中自由流湍流强度通常低到0.05%。. 对于内部流动，入口的湍流强度完全依赖于上游流动的历史，如果上游流动没有完全发展或者没有被扰动，你就可以使用低湍流强度。如果流动完全发展，湍流强度可能就达到了百分之几。完全发展的管流的核心的湍流强度可以用下面的经验公式计算： ()81Re 16.0-?'≡H D avg u u I
浅谈边界条件 对有限元计算，无论是ansys，abaqus，msc还是comsol等，归结为一句话就是解微分方程。而解方程要有定解，就一定要引入条件，这些附加条件称为定解条件。定解条件的形式很多，只讨论最常见的两种——初始条件和边界条件。 在说边界条件之前，先谈谈初值问题和边值问题。 初值和边值问题： 对一般的微分方程，求其定解，必须引入条件，这个条件大概分两类---初始条件和边界条件，如果方程要求未知量y(x)及其导数y′(x)在自变量的同一点x=x0取给定的值，即y(x0 )=y0，y′(x0)= y0′，则这种条件就称为初始条件，由方程和初始条件构成的问题就称为初值问题； 而在许多实际问题中，往往要求微分方程的解在在某个给定的区间a ≤ x ≤b的端点满足一定的条件，如y(a) = A , y(b) = B 则给出的在端点(边界点)的值的条件，称为边界条件，微分方程和边界条件构成数学模型就称为边值问题。 三类边界条件: 边值问题中的边界条件的形式多种多样，在端点处大体上可以写成这样的形式，Ay+By'=C，若B=0，A≠0,则称为第一类边界条件或狄里克莱(Dirichlet)条件；B≠0,A=0，称为第二类边界条件或诺依曼(Neumann)条件；A≠0,B≠0，则称为第三类边界条件或洛平(Robin)条件。 总体来说， 第一类边界条件：给出未知函数在边界上的数值； 第二类边界条件：给出未知函数在边界外法线的方向导数； 第三类边界条件：给出未知函数在边界上的函数值和外法向导数的线性组合。 对应于comsol，只有两种边界条件： Dirichlet boundary（第一类边界条件）—在端点，待求变量的值被指定。 Neumann boundary（第二类边界条件）—待求变量边界外法线的方向导数被指定。 再补充点初始条件： 初始条件，是指过程发生的初始状态，也就是未知函数及其对时间的各阶偏导数在初始时刻t＝0的值．在有限元中，好多初始条件要预先给定的。不同的场方程对应不同的初始条件。 总之，为了确定泛定方程的解，就必须提供足够的初始条件和边界条件！
第五章，边界条件 5-1, FLUENT 程序边界条件种类 进口 出口 壁面 orifice (interior) orifice_plate and orifice_plate-shadow 流体 Example: Face and Cell zones associated with Pipe Flow through orifice plate FLUENT 的边界条件包括： 1， 流动进、出口边界条件 2， 壁面，轴对称和周期性边界 3， Internal cell zones: fluid, solid (porous is a type of fluid zone ) 4， Internal face boundaryies: fan, radiator, porous jump, wall, interior 5－2，流动进口、出口边界条件 FLUENT 提供了10种类型的流动进、出口条件，它们分别是：
一般形式： 可压缩流动： 压力进口 质量进口 压力出口 压力远场 不可压缩流动： 特殊进出口条件： 速度进口 进口通分，出口通风 自由流出 吸气风扇，排气风扇 1， 速度进口：给出进口速度及需要计算的所有标量值 2， 压力进口：给出进口的总压和其它需要计算的标量进口值 3， 质量流进口：主要用于可压缩流动，给出进口的质量流量。对于不可压缩流动，没有必要 给出该边界条件，因为密度是常数，我们可以用速度进口条件。 4， 压力出口：给定流动出口的静压。对于有回流的出口，该边界条件比outflow 边界条件更 容易收敛。 5， 压力远场：该边界条件只对可压缩流动适合。 6， outflow ： 该边界条件用以模拟在求解问题之前，无法知道出口速度或者压力；出口流动 符合完全发展条件，出口处，除了压力之外，其它参量梯度为零。该边界条件不适合可压缩流动。 7， inlet vent ：进口风扇条件需要给定一个损失系数，流动方向和环境总压和总温。 8， intake fan ：进口风扇条件需要给定压降，流动方向和环境总压和总温。 9， out let vent ：排出风扇给定损失系数和环境静压和静温。 10, exhaust fan.：排除风扇给定压降，环境静压。 5－3 压力进口边界条件 压力进口边界条件通常用于给出流体进口的压力和流动的其它标量参数，对计算可压和不可压问题都适合。压力进口边界条件通常用于不知道进口流率或流动速度时候的流动，这类流动在工程中常见，如浮力驱动的流动问题。压力进口条件还可以用于处理外部或者非受限流动的自由边界。 压力边界条件需要表压输入。 5－1 operating gauge absolute p p p +=Operating pressure 输入： Define-operating conditions
在流场的入口、出口和远场边界上，用户需要定义流场的湍流参数。在FLUENT 中可以使用的湍流模型有很多种。在使用各种湍流模型时，哪些变量需要设定，哪些不需要设定以及如何给定这些变量的具体数值，都是经常困扰用户的问题。本小节只讨论在边界上设置均匀湍流参数的方法，湍流参数在边界上不是均匀分布的情况可以用型函数和UDF（用户自定 义函数）来定义，具体方法请参见相关章节的叙述。 在大多数情况下，湍流是在入口后面一段距离经过转捩形成的，因此在边界上设置均匀湍流条件是一种可以接受的选择。特别是在不知道湍流参量的分布规律时，在边界上采用均匀湍流条件可以简化模型的设置。在设置边界条件时，首先应该定性地对流动进行分析，以便边界条件的设置不违背物理规律。违背物理规律的参数设置往往导致错误的计算结果，甚至使计算发散而无法进行下去。 在Turbulence Specification Method （湍流定义方法）下拉列表中，可以简 单地用一个常数来定义湍流参数，即通过给定湍流强度、湍流粘度比、水力直径或湍流特征长在边界上的值来定义流场边界上的湍流。下面具体讨论这些湍流参数的含义，以保证在设置模型时不出现违背流动规律的错误设置： （1）湍流强度（Turbulence Intensity） 湍流强度I的定义为：I=Sq rt(u’*u’+v’*v’+w’*w’)/u_avg (8-1) 上式中u',v' 和w' 是速度脉动量，u_avg是平均速度。 湍流强度小于1％时，可以认为湍流强度是比较低的， 而在湍流强度大于10％时，则可以认为湍流强度是比较高的。 在来流为层流时，湍流强度可以用绕流物体的几何特征粗略地估算出来。比如在模拟风洞试验的计算中，自由流的湍流强度可以用风洞的特征长度估计出来。在现代的低湍流度风洞中，自由流的湍流强度通常低于0.05%。 内流问题进口处的湍流强度取决于上游流动状态。如果上游是没有充分发展的未受扰流动，则进口处可以使用低湍流强度。如果上游是充分发展的湍流，则进口处湍流强度可以达到几个百分点。 如果管道中的流动是充分发展的湍流，则湍流强度可以用公式(8-2)计算得到，这个公式是从管流经验公式得到的： I=u’/u_avg=0.16*Re_DH^-0.125 (8-2) 其中Re_DH是Hydraulic Diameter（水力直径）的意思，即式(8-2)中的雷诺数是以水力直径为特征长度求出的。 （2）湍流的长度尺度与水力直径 湍流能量主要集中在大涡结构中，而湍流长度尺度l则是与大涡结构相关的物理量。在充分发展的管流中，因为漩涡尺度不可能大于管道直径，所以l 是受到 管道尺寸制约的几何量。湍流长度尺度l 与管道物理尺寸L关系可以表示为： l = 0.07L (8-3) 式中的比例因子0.07 是充分发展管流中混合长的最大值，而L则是管道直径。在管道截面不是圆形时，L可以取为管道的水力直径。 湍流的特征长取决于对湍流发展具有决定性影响的几何尺度。在上面的讨论中，管道直径是决定湍流发展过程的唯一长度量。如果在流动中还存在其他对流动影响更大的物体，比如在管道中存在一个障碍物，而障碍物对湍流的发生和发展过程起着重要的干扰作用。在这种情况下，湍流特征长就应该取为障碍物的特征长度。
定义边界条件概述 边界条件包括流动变量和热变量在边界处的值。它是FLUENT分析得很关键的一部分，设定边界条件必须小心谨慎。 边界条件的分类：进出口边界条件：压力、速度、质量进口、进风口、进气扇、压力出口、压力远场边界条件、质量出口、通风口、排气扇；壁面、repeating, and pole boundaries:壁面，对称，周期，轴；内部单元区域：流体、固体(多孔是一种流动区域类型) ；内部表面边界：风扇、散热器、多孔跳跃、壁面、内部。(内部表面边界条件定义在单元表面，这意味着它们没有有限厚度，并提供了流场性质的每一步的变化。这些边界条件用来补充描述排气扇、细孔薄膜以及散热器的物理模型。内部表面区域的内部类型不需要你输入任何东西。) 下面一节将详细介绍上面所叙述边界条件，并详细介绍了它们的设定方法以及设定的具体合适条件。周期性边界条件在本章中介绍，模拟完全发展的周期性流动将在周期性流动和热传导一章中介绍。 使用边界条件面板 边界条件(Figure 1)对于特定边界允许你改变边界条件区域类型，并且打开其他的面板以设定每一区域的边界条件参数 菜单：Define/Boundary Conditions... Figure 1: 边界条件面板 改变边界区域类型 设定任何边界条件之前，必须检查所有边界区域的区域类型，如有必要就作适当的修改。比方说：如果你的网格是压力入口，但是你想要使用速度入口，你就要把压力入口改为速度入口之后再设定。 改变类型的步骤如下：: 1.在区域下拉列表中选定所要修改的区域 2.在类型列表中选择正确的区域类型 3.当问题提示菜单出现时，点击确认 确认改变之后，区域类型将会改变，名字也将自动改变(如果初始名字时缺省的请参阅边界条件区域名字一节),设定区域边界条件的面板也将自动打开。 ！注意：这个方法不能用于改变周期性类型，因为该边界类型已经存在了附加限制。创建边界条件一节解释了如何创建和分开周期性区域。需要注意的是，只能在图一中每一个类别中改变边界类型(注意：双边区域表面是分离的不同单元区域.) Figure 1: 区域类型的分类列表 设定边界条件 在FLUENT中，边界条件和区域有关而与个别表面或者单元无关。如果要结合具有相同边界条件的两个或更多区域请参阅合并区域一节。 设定每一特定区域的边界条件，请遵循下面的步骤： 1.在边界条件区域的下拉列表中选择区域。 2. 点击Set...按钮。或者，1.在区域下拉列表中选择区域。 2.在类型列表中点击所要选择的类型。或者在区域列表中双击所需区域.，选择边界条件区域将会打开，并且你可以指定适当的边界条件
在入口、出口或远场边界流入流域的流动，FLUENT 需要指定输运标量的值。本节描述了对于特定模型需要哪些量，并且该如何指定它们。也为确定流入边界值最为合适的方法提供了指导方针。 使用轮廓指定湍流参量 在入口处要准确的描述边界层和完全发展的湍流流动，你应该通过实验数据和经验公式创建边界轮廓文件来完美的设定湍流量。如果你有轮廓的分析描述而不是数据点，你也可以用这个分析描述来创建边界轮廓文件，或者创建用户自定义函数来提供入口边界的信息。一旦你创建了轮廓函数，你就可以使用如下的方法： ● Spalart-Allmaras 模型：在湍流指定方法下拉菜单中指定湍流粘性比，并在在湍流粘性 比之后的下拉菜单中选择适当的轮廓名。通过将m_t/m 和密度与分子粘性的适当结合， FLUENT 为修改后的湍流粘性计算边界值。 ● k-e 模型：在湍流指定方法下拉菜单中选择K 和Epsilon 并在湍动能（Turb. Kinetic Energy ）和湍流扩散速度（Turb. Dissipation Rate ）之后的下拉菜单中选择适当的轮廓名。 ● 雷诺应力模型：在湍流指定方法下拉菜单中选择K 和Epsilon 并在湍动能（Turb. Kinetic Energy ）和湍流扩散速度（Turb. Dissipation Rate ）之后的下拉菜单中选择适当的轮廓名。在湍流指定方法下拉菜单中选择雷诺应力部分，并在每一个单独的雷诺应力部分之后的下拉菜单中选择适当的轮廓名。 湍流量的统一说明 在某些情况下流动流入开始时，将边界处的所有湍流量指定为统一值是适当的。比如说，在进入管道的流体，远场边界，甚至完全发展的管流中，湍流量的精确轮廓是未知的。 在大多数湍流流动中，湍流的更高层次产生于边界层而不是流动边界进入流域的地方，因此这就导致了计算结果对流入边界值相对来说不敏感。然而必须注意的是要保证边界值不是非物理边界。非物理边界会导致你的解不准确或者不收敛。对于外部流来说这一特点尤其突出，如果自由流的有效粘性系数具有非物理性的大值，边界层就会找不到了。 你可以在使用轮廓指定湍流量一节中描述的湍流指定方法，来输入同一数值取代轮廓。你也可以选择用更为方便的量来指定湍流量，如湍流强度，湍流粘性比，水力直径以及湍流特征尺度，下面将会对这些内容作一详细叙述。 湍流强度I 定义为相对于平均速度u_avg 的脉动速度u^'的均方根。 小于或等于1%的湍流强度通常被认为低强度湍流，大于10%被认为是高强度湍流。从外界，测量数据的入口边界，你可以很好的估计湍流强度。例如：如果你模拟风洞试验，自由流的湍流强度通常可以从风洞指标中得到。在现代低湍流风洞中自由流湍流强度通常低到0.05%。. 对于内部流动，入口的湍流强度完全依赖于上游流动的历史，如果上游流动没有完全发展或者没有被扰动，你就可以使用低湍流强度。如果流动完全发展，湍流强度可能就达到了百分之几。完全发展的管流的核心的湍流强度可以用下面的经验公式计算： ()81Re 16.0-?'≡H D u u I
CFD-边界条件 定义边界条件概述 边界条件包括流动变量和热变量在边界处的值。它是FLUENT分析得很关键的一部分，设定边界条件必须小心谨慎。 边界条件的分类：进出口边界条件：压力、速度、质量进口、进风口、进气扇、压力出口、压力远场边界条件、质量出口、通风口、排气扇；壁面、repeating, and pole boundaries:壁面，对称，周期，轴；内部单元区域：流体、固体(多孔是一种流动区域类型) ；内部表面边界：风扇、散热器、多孔跳跃、壁面、内部。(内部表面边界条件定义在单元表面，这意味着它们没有有限厚度，并提供了流场性质的每一步的变化。这些边界条件用来补充描述排气扇、细孔薄膜以及散热器的物理模型。内部表面区域的内部类型不需要你输入任何东西。) 下面一节将详细介绍上面所叙述边界条件，并详细介绍了它们的设定方法以及设定的具体合适条件。周期性边界条件在本章中介绍，模拟完全发展的周期性流动将在周期性流动和热传导一章中介绍。 使用边界条件面板 边界条件(Figure 1)对于特定边界允许你改变边界条件区域类型，并且打开其他的面板以设定每一区域的边界条件参数 菜单：Define/Boundary Conditions... Figure 1: 边界条件面板 改变边界区域类型 设定任何边界条件之前，必须检查所有边界区域的区域类型，如有必要就作适当的修改。比方说：如果你的网格是压力入口，但是你想要使用速度入口，你就要把压力入口改为速度入口之后再设定。 改变类型的步骤如下：: 1.在区域下拉列表中选定所要修改的区域 2.在类型列表中选择正确的区域类型 3.当问题提示菜单出现时，点击确认 确认改变之后，区域类型将会改变，名字也将自动改变 (如果初始名字时缺省的请参阅边界条件区域名字一节),设定区域边界条件的面板也将自动打开。 ！注意：这个方法不能用于改变周期性类型，因为该边界类型已经存在了附加限制。创建边界条件一节解释了如何创建和分开周期性区域。需要注意的是，只能在图一中每一个类别中改变边界类型(注意：双边区域表面是分离的不同单元区域.) Figure 1: 区域类型的分类列表 设定边界条件 在FLUENT中，边界条件和区域有关而与个别表面或者单元无关。如果要结合具有相同边界条件的两个或更多区域请参阅合并区域一节。 设定每一特定区域的边界条件，请遵循下面的步骤：