CFX边界条件及CCL语言ppt课件.ppt

© 2005 ANSYS CHINA, Inc.3-1CFX 10.0边界条件边界条件/CFX/CFX表达式语言表达式语言讲座 3© 2005 ANSYS CHINA, Inc.3-2CFX 10.0l在计算域的所有外部边界上都需要l尽量采用意义明确的名字〔你不需要接受系统的默认名字）l选择边界所在的域 (当计算中有多个域的时候适用).边界条件© 2005 ANSYS CHINA, Inc.3-3CFX 10.0l边界类型l进口，出口，开放式，墙，对称面l位置l选择相应的二维几何体或组合l坐标系l如果有多于一个，请选择合适的l坐标类型l只在旋转坐标系中可选这个选项允许用户指定基于静止或旋转坐标系的某些参数基本设定基本设定© 2005 ANSYS CHINA, Inc.3-4CFX 10.0l可选项取决于边界类型及域的设定l例如：超音速的可选性取决于域设定里面的热传递模型l在右边例子中，只模拟了湍流 (i.e. 没有传热或多组分/多相流模拟）边界条件详细内容边界条件详细内容© 2005 ANSYS CHINA, Inc.3-5CFX 10.0l进口l流体自进口流入计算域: 用白色箭头表示l出口l流体自出口流出计算域: 用黄色箭头表示l开放式边界l根据当地条件，流体可以流入或流出计算域：用蓝色双向箭头表示. 设定过程和进口类似，也需要设定流体方向和压力l墙l用绿色的八角体表示. 非滑移/自由滑移, 热传递性质和墙壁的粗糙度等特性都可以设定l对称面l流体在该面一侧的特性是另一面的镜像. 当对称流动存在时可以用来减少计算网格的数量.边界类型边界类型© 2005 ANSYS CHINA, Inc.3-6CFX 10.0预设边界条件•可以通过对已知数据文件的插值来设定边界条件.•我们经常使用以前的模拟结果或已知的实验数据作为当前模拟的边界条件•CFX-Pre 可以用插分功能来根据读入的数据生成CEL 表达式. •当我们采用“Profile〞方法设定边界条件时，这些数据是自动生成的.© 2005 ANSYS CHINA, Inc.3-7CFX 10.0l创建一个边界条件文件。

这个可以用CFX-Post当中的输出功能实现从CFX-Post中产生的边界条件文件带颜色的格中的信息是需要的，CFX-Pre会自动读入施加预设边界条件的步骤施加预设边界条件的步骤© 2005 ANSYS CHINA, Inc.3-8CFX 10.0l在CFX-Pre中读入边界条件文件这个可以通过 Tools>Initialize Profile Data 并选择相应的文件来实现我们也可以载入多个边界条件文件，每一个文件也可以应用在多个地方施加边界条件文件的步骤施加边界条件文件的步骤© 2005 ANSYS CHINA, Inc.3-9CFX 10.0l在边界条件处勾选：Use Profile Data l从下拉菜单中选择合适的文件，然后点击 Generate Values. 施加边界条件文件的步骤施加边界条件文件的步骤© 2005 ANSYS CHINA, Inc.3-10CFX 10.0l边界条件的详细录入界面会有相应改变，以适应文件输入只有当你点击了Apply以后，这些改变才会生效施加边界条件文件的步骤施加边界条件文件的步骤© 2005 ANSYS CHINA, Inc.3-11CFX 10.0l边界条件文件〔以及别的边界条件〕可以通过在CFX-Pre的边界条件编辑器中选择绘图选项来实现可视化。

用户可以根据边界条件文件创建边界的等值线或矢量图l每次求解器启动时，边界条件文件将被读入也就是说边界条件文件可以在两次模拟当中进行修改，而不必回到CFX-Pre当中施加边界条件文件的步骤施加边界条件文件的步骤© 2005 ANSYS CHINA, Inc.3-12CFX 10.0l在多处地方施加边界条件文件是可能的:l对于有相同法向量 X, Y 或 Z的边界, 可以把边界条件输出成2D 数据文件(比如，对两个有相同法向Z的边界来说，就输出X和Y方向的数据). 第一个边界上的数据可以不经修改用于第二个边界上l如果两个边界没有相同的法向，那么可以直接修改你的边界条件文件使第一个边界上的条件适合第二个边界在多处地方施加边界条件文件在多处地方施加边界条件文件© 2005 ANSYS CHINA, Inc.3-13CFX 10.0l有些变量需要一个前缀来表明材料的名字 (例如 air.vf 就表示空气的体积百分比)l非标准变量名和用户定义的变量名标准变量名标准变量名© 2005 ANSYS CHINA, Inc.3-14CFX 10.0l对一维的离散数据来说，数据间的拓扑关系可以通过对原始数据的坐标排序来确定。

排序后的原始数据就可以应用线形插分数据通过这种方式来分类，所以是否提前特殊规定并不重要l对于二维和三维的离散数据来说，“点云图〞的方法被用来实现插分整个过程包括对被插分点最近的三个原始数据点的快速定位，以及根据它们离被插分点距离远近的一个反向加权的平均过程 如果某一个原始数据点正好落在该插分点的位置上，那么我们就用这个原始数据点的值作为该插分点的值l在求解过程当中，根据离散和数值积分过程的具体要求，求解器需要各个不同位置上的值，比如在积分点上，节点上和各个面的中心点上 在所有的情况下，这个需要的位置都会被确定，原始的数据就会被插分到该位置数据插分方法数据插分方法数据插分方法© 2005 ANSYS CHINA, Inc.3-15CFX 10.0l当边界条件文件被CFX-Solver读入后，如果文件大小不超过16K，它会被写到.out文件的Profile Data 部分16K是一个默认值，这个默认值是可以改变的所有的文件里的数据都会被写入结果文件，并可以在DOS环境下用cfx5dfile命令分离出来 从结果文件中分离出边界条件文件从结果文件中分离出边界条件文件© 2005 ANSYS CHINA, Inc.3-16CFX 10.0l在当前结果文件中找到相应的数据文件，我们需要键入:lcfx5dfile file_001.res –list-profile-files (这个命令会列出结果文件中所有的边界条件文件，每行列一个) l如果要找到某个特定的边界条件文件，键入:lcfx5dfile file_001.res –read-profile-file .csvl或者:l cfx5dfile file_001.res –extract-profile-file .csvl 从结果文件中分离出边界条件文件从结果文件中分离出边界条件文件© 2005 ANSYS CHINA, Inc.3-17CFX 10.0lCEL - CFX Expression LanguagelCEL 是一种说明性的编译语言，可以让用户不借助外挂Fortran子程序而对模拟进行强化l - 可以涉及 CFX 内部变量l - are evaluated by the CFX Solver and CFX-PostCFX Expression LanguageCFX Expression LanguageCFX Expression Language© 2005 ANSYS CHINA, Inc.3-18CFX 10.0lCEL是说明性的l描述一个表达式的名字和定义〔也可以加入注释等）l语句的格式必须符合提前设定好的语法。

这个语法和大多数的编程语言中的数学表达式是非常类似的l语句必须包含下列:l一个数字〔如果需要，也可以带单位）l一个或多个表示常数，系统变量，用户自定义变量，函数或其它CEL表达式，用 +, -, * , / 和 ^ 分隔也可以用括号来进行分组以上表达式的语法规则和传统数学规则是一样的以上表达式的语法规则和传统数学规则是一样的CEL StatementsCEL StatementsCEL Statements© 2005 ANSYS CHINA, Inc.3-19CFX 10.0l变量和表达式应区分大小写l用括号来规定运算的优先级l对加减运算，表达式的单位必须统一l1.0 [mm] + 0.45 [yds] (valid)l2.5 [s m^-1] - (3.0 [m s^-1])^-1 (valid)l1.0 [mm] + 0.3 [kg] (invalid)表达式规则表达式规则© 2005 ANSYS CHINA, Inc.3-20CFX 10.0l分数幂和小数幂是允许的la^(1/2) (valid)l1.0^0.5 (valid) l符合运算的表达式单位并不需明确写出---靠计算结果确定l(a [kg m^-3] * b [m s^-1]) has units of [kg m^-2 s^-1]表达式规则表达式规则© 2005 ANSYS CHINA, Inc.3-21CFX 10.0l系统中定义了很多系统变量以供调用l这些系统变量可以被用在任意表达式l系统变量中已经包含了单位l变量列表可以:l - 通过在Pre里点击 System Variables 或在Expression Editor里点击Functions 按钮来列出l - 在CFX-Post里通过观察所有可选变量的方法列出系统变量系统变量© 2005 ANSYS CHINA, Inc.3-22CFX 10.0xDirection 1 in Reference Coordinate FrameyDirection 2 in Reference Coordinate FramezDirection 3 in Reference Coordinate FramerRadial spatial location, r = (x^2+y^2)^0.5thetaAngle, arctan(y/x)tTimeuVelocity in the x coordinate directionvVelocity in the y coordinate directionwVelocity in the z coordinate directionp(absolute) PressurekeTurbulent kinetic energyedTurbulent eddy dissipationTTemperaturesstrnrShear strain ratedensityDenstiyrNoDimNon-dimensional radius (rotating frame only)viscosityDynamic ViscosityCpSpecific Heat Capacity at Constant PressurecondThermal ConductivityenthalpySpecific EnthalpybetaThermal ExpansivityspeedofsoundLocal speed of sound in fluidsubdomainSub-domain variable (1.0 in Sub-domain, 0.0 elsewhere)mean diameterMean DiameterdeneffEffective DensityAV nameAdditional Variable namemfMass Fraction系统变量系统变量© 2005 ANSYS CHINA, Inc.3-23CFX 10.0l根据模拟的类型和创建的表达式的形式，CFX-5 系统变量和用户定义的表达式会变为可选或不可选l在某些情况下，系统变量会成为逻辑上不可选l时间 (t) 在稳态模拟时不可选l温度 (T) 在系统没有传热问题的时候不可选l在其它情况下，某些系统变量还问因为所选的物理模型的原因变成不可选系统变量系统变量© 2005 ANSYS CHINA, Inc.3-24CFX 10.0l一些数学上的函数和操作符在CEL里也可以运用l用户也可以通过FORTRAN语言创建自己的函数Function Operand’s Dimensions [x]Operand’s ValuesResult’s Dimensionssin(x) Angle Any Dimensionlesscos(x) Angle Any Dimensionlesstan(x) ** Angle Any Dimensionlessasin(x) Dimensionless -1 £ x £ 1 Angleacos(x) Dimensionless -1 £ x £ 1 Angleatan(x) Dimensionless Any Angleexp(x) Dimensionless Any Dimensionlessloge(x) Dimensionless 0 < x Dimensionlesslog10(x) Dimensionless 0 < x Dimensionlessabs(x) Any Any [x]sqrt(x) Any 0 £ x [x]^0.5min(x,y) *** Any Any [x]max(x,y) *** Any Any [x]step(x) * Dimensionless Any Dimensionless*step(x) is 0 for negative x, 1 for positive x and 0.5 for x=0.** note that tan(x) is undefined for np/2 where n=1, 3, 5 .. .*** both x and y must have the same dimensions.内置函数内置函数内置函数/ / /常数常数常数© 2005 ANSYS CHINA, Inc.3-25CFX 10.0内置运算符/函数和常数l在CEL里一些常用的数学常数也可以被用来创建新的表达式:le Constant: 2.7182818lg Acceleration due to gravity: 9.806 [m s^-2]lpi Constant: 3.1415927lR Universal Gas Constant: 8314.5 [m^2 s^-2 K^-1]l用户可以定义自己的1D线性或3D云图式的插分函数l在输入的数据点之间线性插分，然后输出一个单个的值l输入的单位和输出的单位由用户决定l如果用户需要的函数在CEL里没有，或者需要使用某些特定的变量〔比如梯度），用户可以通过链接到一个FORTRAN函数库来自定义函数© 2005 ANSYS CHINA, Inc.3-26CFX 10.0•CFX-Pre•从文件读入•在Expression Editor里创建•在需要的地方直接输入•修改求解器的CCL文件•编辑 .DEF 或 .RES 文件•通过DOS控制行修改•后处理中的表达式可以通过这种方式定义，让求解器使用定义表达式的方法定义表达式的方法© 2005 ANSYS CHINA, Inc.3-27CFX 10.0数学表达式粘性剪切流体的粘度: 这儿  代表剪切应力率CEL 等效表达式粘度写成温度的函数K = 10.0 [kg m^-1 s^-0.5]n = 1.5ViscT = K * sstrnr ^ (n-1)orViscT = K*(min(UpperS,(sstrnr+LowerS))^(n-1))这儿 sstrnr 是一个系统函数，表示剪切应力率CELCEL的第二种表达式包含了对应力率的上限和下限，这主要的第二种表达式包含了对应力率的上限和下限，这主要是为了使计算结果在物理上合理。

是为了使计算结果在物理上合理CEL CEL 例子例子: : 可变粘度可变粘度© 2005 ANSYS CHINA, Inc.3-28CFX 10.0在 CFX-Preexpression editor 有一个calculate的功能，可以用来测试表达式，或者画出一维的结果有些值可能需要手工输入，因为系统变量没有值在求解器中当需要值的时候，该表达式就会被检验初始值: 在计算开始的时候随时间变化的边界条件: 在每一个时间步的开始流体特性: 在两个时间步长之间的内循环里累积量: 在每一个时间步长的开始检验表达式检验表达式© 2005 ANSYS CHINA, Inc.3-29CFX 10.0l表达式编辑器 l在CFX中开发和管理表达式的一个有效的工具l从很多面板上以及主菜单的Create选项上都可以进入表达式编辑器表达式编辑器表达式编辑器 © 2005 ANSYS CHINA, Inc.3-30CFX 10.0l下面例子演示了如何把粘度设定为温度的函数. 粘度和温度之间的关系是：CEL CEL CEL 例子例子例子© 2005 ANSYS CHINA, Inc.3-31CFX 10.0l或者，你可以用任何的文字编辑器编写出含有上述函数式的文本文件，然后用表达式编辑器读入l例子: lC1 = 10. [ K^-1 ] # constantlC2 = 1. [ kg/ms ] # constantlvis = C2*exp(-C1*T) # viscosityCEL CEL CEL 例子例子例子© 2005 ANSYS CHINA, Inc.3-32CFX 10.0CEL CEL 例子例子 ( (边界条件边界条件) )l下面例子演示了如何用速度的U和V分量给一个旋转的强设上角速度:© 2005 ANSYS CHINA, Inc.3-33CFX 10.0l在指定位置允许对非逻辑合成变量进行计算.l例子:l计算在等值面上Cp 的面积平均:lareaAve(Cp)iso1l经过某一个位置的某种流体的流量:loil.massFlow()slice1l注意:l“” - 必须作用于某一个几何位置.l流体中的某相或某个组分可以通过下面方式选取 [.][.]对于高级用户来说，在求解器中对合成变量进行操作也是可能的。

这些变量可以被输入到表达式中，也可以在Solver Manager中被观测〔详见文献中关于Output Control的描述）合成变量合成变量© 2005 ANSYS CHINA, Inc.3-34CFX 10.0l系统中有一些已经定义好了的合成变量l0 或 1 取决于具体的函数---参看帮助文档lIf 1 argument:l可以是Post当中的一个表达式; 在求解器当中只允许以变量形式存在l返回值得单位取决于宗数的单位合成变量合成变量© 2005 ANSYS CHINA, Inc.3-35CFX 10.0EXPRESSIONS: ReqT = 350 [K] TempOut = areaAve(T)outlet TCFilter = step(TempOut/1[K]-ReqT/1[K]) TCTemp = 400[K]*TCFilter+285[K]*(1-TCFilter) TCFlow = 10[m/s]*TCFilter+2[m/s]*(1-TCFilter)END[…]BOUNDARY : TempControl Boundary Type = INLET Location = TempControl Coord Frame = Coord BOUNDARY CONDITIONS : MASS AND MOMENTUM : Option = Normal Speed Normal Speed = TCFlow END HEAT TRANSFER : Option = Static Temperature Static Temperature = TCTemp END ENDEND合成变量合成变量© 2005 ANSYS CHINA, Inc.3-36CFX 10.0在CFX-Post 中可以用如下方式定义表达式:用CFX-Post 中的Expression Editor直接在当地输入表达式直接在commandline (Line input mode)或 Command Editor处输入从CCL文件中读入从session 文件中读入从state 文件中读入在在CFX-PostCFX-Post当中定义表达式当中定义表达式© 2005 ANSYS CHINA, Inc.3-37CFX 10.0lCFX-Post 还提供一系列额外的函数:l[.][_[_]]([])larea面积大小 (projected to axis optional)lareaAve面积加权平均lareaInt面积加权积分lave代数平均lcount计算数目lforce某个面上特定方向受力lforceNorm指定方向的法向力大小llength曲线长度llengthAve长度加权平均llengthInt长度加权积分lmassFlow总的质量流量lmassFlowAve质量加权平均lmassFlowInt质量加权积分lmaxVal最大值lminVal最小值CFX-Post CFX-Post CFX-Post 中的函数中的函数中的函数© 2005 ANSYS CHINA, Inc.3-38CFX 10.0l接上页:l[.][_[_]]([])lprobe某点上的值lsum求总和ltorque在指定轴下对于某个轴的力矩lvolume某个三维体的体积lvolumeAve体积加权平均lvolumeInt体积加权积分CFX-Post CFX-Post CFX-Post 函数函数函数© 2005 ANSYS CHINA, Inc.3-39CFX 10.0数学表达式压力系数:CEL 等效压力系数Pref = massFlowAve(P)inflowDynH = massFlowAve(.5*Density*V^2)inflowCpress = (P - Pref)/DynHorCpress = (P - massFlowAve(P)inflow) / massFlowAve(.5*Density*V^2)inflowCEL CEL 例子例子: : 压力系数压力系数。