《边界条件CFX表达式语言ppt课件.ppt》由会员分享，可在线阅读，更多相关《边界条件CFX表达式语言ppt课件.ppt（39页珍藏版）》请在三一办公上搜索。
1、边界条件/CFX表达式语言,讲座 3,在计算域的所有外部边界上都需要尽量采用意义明确的名字（你不需要接受系统的默认名字）选择边界所在的域 (当计算中有多个域的时候适用).,边界条件,边界类型进口，出口，开放式，墙，对称面位置选择相应的二维几何体或组合坐标系如果有多于一个，请选择合适的坐标类型只在旋转坐标系中可选。这个选项允许用户指定基于静止或旋转坐标系的某些参数,基本设定,可选项取决于边界类型及域的设定例如：超音速的可选性取决于域设定里面的热传递模型在右边例子中，只模拟了湍流 (i.e. 没有传热或多组分/多相流模拟）,边界条件详细内容,进口流体自进口流入计算域: 用白色箭头表示出口流体自出口
2、流出计算域: 用黄色箭头表示开放式边界根据当地条件，流体可以流入或流出计算域：用蓝色双向箭头表示. 设定过程和进口类似，也需要设定流体方向和压力墙用绿色的八角体表示. 非滑移/自由滑移, 热传递性质和墙壁的粗糙度等特性都可以设定对称面流体在该面一侧的特性是另一面的镜像. 当对称流动存在时可以用来减少计算网格的数量.,边界类型,预设边界条件,可以通过对已知数据文件的插值来设定边界条件.我们经常使用以前的模拟结果或已知的实验数据作为当前模拟的边界条件CFX-Pre 可以用插分功能来根据读入的数据生成CEL 表达式. 当我们采用“Profile”方法设定边界条件时，这些数据是自动生成的.,创建一个边
3、界条件文件。这个可以用CFX-Post当中的输出功能实现。,从CFX-Post中产生的边界条件文件。带颜色的格中的信息是需要的，CFX-Pre会自动读入。,施加预设边界条件的步骤,在CFX-Pre中读入边界条件文件。这个可以通过 ToolsInitialize Profile Data 并选择相应的文件来实现。我们也可以载入多个边界条件文件，每一个文件也可以应用在多个地方。,施加边界条件文件的步骤,在边界条件处勾选：Use Profile Data 从下拉菜单中选择合适的文件，然后点击 Generate Values.,施加边界条件文件的步骤,边界条件的详细录入界面会有相应改变，以适应文件输入
4、。只有当你点击了Apply以后，这些改变才会生效。,施加边界条件文件的步骤,边界条件文件（以及别的边界条件）可以通过在CFX-Pre的边界条件编辑器中选择绘图选项来实现可视化。用户可以根据边界条件文件创建边界的等值线或矢量图。每次求解器启动时，边界条件文件将被读入。也就是说边界条件文件可以在两次模拟当中进行修改，而不必回到CFX-Pre当中。,施加边界条件文件的步骤,在多处地方施加边界条件文件是可能的:对于有相同法向量 X, Y 或 Z的边界, 可以把边界条件输出成2D 数据文件(比如，对两个有相同法向Z的边界来说，就输出X和Y方向的数据). 第一个边界上的数据可以不经修改用于第二个边界上。如
5、果两个边界没有相同的法向，那么可以直接修改你的边界条件文件使第一个边界上的条件适合第二个边界。,在多处地方施加边界条件文件,有些变量需要一个前缀来表明材料的名字 (例如 air.vf 就表示空气的体积百分比)非标准变量名和用户定义的变量名,标准变量名,对一维的离散数据来说，数据间的拓扑关系可以通过对原始数据的坐标排序来确定。排序后的原始数据就可以应用线形插分。数据通过这种方式来分类，所以是否提前特殊规定并不重要。对于二维和三维的离散数据来说，“点云图”的方法被用来实现插分。整个过程包括对被插分点最近的三个原始数据点的快速定位，以及根据它们离被插分点距离远近的一个反向加权的平均过程。 如果某一个
6、原始数据点正好落在该插分点的位置上，那么我们就用这个原始数据点的值作为该插分点的值。在求解过程当中，根据离散和数值积分过程的具体要求，求解器需要各个不同位置上的值，比如在积分点上，节点上和各个面的中心点上。 在所有的情况下，这个需要的位置都会被确定，原始的数据就会被插分到该位置。,数据插分方法,当边界条件文件被CFX-Solver读入后，如果文件大小不超过16K，它会被写到.out文件的Profile Data 部分。16K是一个默认值，这个默认值是可以改变的。所有的文件里的数据都会被写入结果文件，并可以在DOS环境下用cfx5dfile命令分离出来,从结果文件中分离出边界条件文件,在当前结果
7、文件中找到相应的数据文件，我们需要键入:cfx5dfile file_001.res list-profile-files (这个命令会列出结果文件中所有的边界条件文件，每行列一个) 如果要找到某个特定的边界条件文件，键入:cfx5dfile file_001.res read-profile-file .csv或者: cfx5dfile file_001.res extract-profile-file .csv,从结果文件中分离出边界条件文件,CEL - CFX Expression LanguageCEL 是一种说明性的编译语言，可以让用户不借助外挂Fortran子程序而对模拟进行强化
8、- 可以涉及 CFX 内部变量 - are evaluated by the CFX Solver and CFX-Post,CFX Expression Language,CEL是说明性的描述一个表达式的名字和定义（也可以加入注释等）语句的格式必须符合提前设定好的语法。这个语法和大多数的编程语言中的数学表达式是非常类似的语句必须包含下列:一个数字（如果需要，也可以带单位）一个或多个表示常数，系统变量，用户自定义变量，函数或其它CEL表达式，用 +, -, * , / 和 分隔。也可以用括号来进行分组。,以上表达式的语法规则和传统数学规则是一样的,CEL Statements,变量和表达式应区
9、分大小写用括号来规定运算的优先级对加减运算，表达式的单位必须统一1.0 mm + 0.45 yds (valid)2.5 s m-1 - (3.0 m s-1)-1 (valid)1.0 mm + 0.3 kg (invalid),表达式规则,分数幂和小数幂是允许的a(1/2) (valid)1.00.5 (valid) 符合运算的表达式单位并不需明确写出-靠计算结果确定(a kg m-3 * b m s-1) has units of kg m-2 s-1,表达式规则,系统中定义了很多系统变量以供调用这些系统变量可以被用在任意表达式系统变量中已经包含了单位变量列表可以: - 通过在Pre里点
10、击 System Variables 或在Expression Editor里点击Functions 按钮来列出 - 在CFX-Post里通过观察所有可选变量的方法列出,系统变量,xDirection 1 in Reference Coordinate FrameyDirection 2 in Reference Coordinate FramezDirection 3 in Reference Coordinate FramerRadial spatial location, r = (x2+y2)0.5thetaAngle, arctan(y/x)tTimeuVelocity in the
11、 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 Viscosity
12、CpSpecific 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 nam
13、emfMass Fraction,系统变量,根据模拟的类型和创建的表达式的形式，CFX-5 系统变量和用户定义的表达式会变为可选或不可选在某些情况下，系统变量会成为逻辑上不可选时间 (t) 在稳态模拟时不可选温度 (T) 在系统没有传热问题的时候不可选在其它情况下，某些系统变量还问因为所选的物理模型的原因变成不可选,系统变量,一些数学上的函数和操作符在CEL里也可以运用用户也可以通过FORTRAN语言创建自己的函数,Function Operands Dimensions xOperands ValuesResults Dimensionssin(x) Angle Any Dimensionl
14、esscos(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 xsqr
15、t(x) Any 0 x x0.5min(x,y) * Any Any xmax(x,y) * Any Any xstep(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.,内置函数/常数,内置运算符/函数和常数,在CEL里一些常用的数学常
16、数也可以被用来创建新的表达式:e Constant: 2.7182818g Acceleration due to gravity: 9.806 m s-2pi Constant: 3.1415927R Universal Gas Constant: 8314.5 m2 s-2 K-1用户可以定义自己的1D线性或3D云图式的插分函数在输入的数据点之间线性插分，然后输出一个单个的值输入的单位和输出的单位由用户决定如果用户需要的函数在CEL里没有，或者需要使用某些特定的变量（比如梯度），用户可以通过链接到一个FORTRAN函数库来自定义函数,CFX-Pre从文件读入在Expression Edit
17、or里创建在需要的地方直接输入修改求解器的CCL文件编辑 .DEF 或 .RES 文件通过DOS控制行修改后处理中的表达式可以通过这种方式定义，让求解器使用,定义表达式的方法,数学表达式粘性剪切流体的粘度:这儿 代表剪切应力率,CEL 等效表达式粘度写成温度的函数K = 10.0 kg m-1 s-0.5n = 1.5ViscT = K * sstrnr (n-1)orViscT = K*(min(UpperS,(sstrnr+LowerS)(n-1)这儿 sstrnr 是一个系统函数，表示剪切应力率,CEL的第二种表达式包含了对应力率的上限和下限，这主要是为了使计算结果在物理上合理。,CEL
18、 例子: 可变粘度,在 CFX-Preexpression editor 有一个calculate的功能，可以用来测试表达式，或者画出一维的结果有些值可能需要手工输入，因为系统变量没有值。在求解器中当需要值的时候，该表达式就会被检验初始值: 在计算开始的时候随时间变化的边界条件: 在每一个时间步的开始流体特性: 在两个时间步长之间的内循环里累积量: 在每一个时间步长的开始,检验表达式,表达式编辑器 在CFX中开发和管理表达式的一个有效的工具从很多面板上以及主菜单的Create选项上都可以进入表达式编辑器,表达式编辑器,下面例子演示了如何把粘度设定为温度的函数. 粘度和温度之间的关系是：,CEL
19、 例子,或者，你可以用任何的文字编辑器编写出含有上述函数式的文本文件，然后用表达式编辑器读入例子: C1 = 10. K-1 # constantC2 = 1. kg/ms # constantvis = C2*exp(-C1*T) # viscosity,CEL 例子,CEL 例子 (边界条件),下面例子演示了如何用速度的U和V分量给一个旋转的强设上角速度:,在指定位置允许对非逻辑合成变量进行计算.例子:计算在等值面上Cp 的面积平均:areaAve(Cp)iso1经过某一个位置的某种流体的流量:oil.massFlow()slice1注意:“” - 必须作用于某一个几何位置.流体中的某相或
20、某个组分可以通过下面方式选取 .,对于高级用户来说，在求解器中对合成变量进行操作也是可能的。这些变量可以被输入到表达式中，也可以在Solver Manager中被观测（详见文献中关于Output Control的描述）,合成变量,系统中有一些已经定义好了的合成变量0 或 1 取决于具体的函数-参看帮助文档If 1 argument:可以是Post当中的一个表达式; 在求解器当中只允许以变量形式存在返回值得单位取决于宗数的单位,合成变量,EXPRESSIONS: ReqT = 350 K TempOut = areaAve(T)outlet TCFilter = step(TempOut/1K-
21、ReqT/1K) TCTemp = 400K*TCFilter+285K*(1-TCFilter) TCFlow = 10m/s*TCFilter+2m/s*(1-TCFilter)ENDBOUNDARY : TempControl Boundary Type = INLET Location = TempControl Coord Frame = Coord BOUNDARY CONDITIONS : MASS AND MOMENTUM : Option = Normal Speed Normal Speed = TCFlow END HEAT TRANSFER : Option = Sta
22、tic Temperature Static Temperature = TCTemp END ENDEND,合成变量,在CFX-Post 中可以用如下方式定义表达式:用CFX-Post 中的Expression Editor直接在当地输入表达式直接在commandline (Line input mode)或 Command Editor处输入从CCL文件中读入从session 文件中读入从state 文件中读入,在CFX-Post当中定义表达式,CFX-Post 还提供一系列额外的函数:._()area面积大小 (projected to axis optional)areaAve面积加权
23、平均areaInt面积加权积分ave代数平均count计算数目force某个面上特定方向受力forceNorm指定方向的法向力大小length曲线长度lengthAve长度加权平均lengthInt长度加权积分massFlow总的质量流量massFlowAve质量加权平均massFlowInt质量加权积分maxVal最大值minVal最小值,CFX-Post 中的函数,接上页:._()probe某点上的值sum求总和torque在指定轴下对于某个轴的力矩volume某个三维体的体积volumeAve体积加权平均volumeInt体积加权积分,CFX-Post 函数,数学表达式压力系数:,CEL 等效压力系数Pref = massFlowAve(P)inflowDynH = massFlowAve(.5*Density*V2)inflowCpress = (P - Pref)/DynHorCpress = (P - massFlowAve(P)inflow) / massFlowAve(.5*Density*V2)inflow,CEL 例子: 压力系数,
