NUMECA帮助文档六

1、第十二章 跨叶片截面模块12.1绪言本章针对透平机械讲述快速三维跨叶片截面模块的分析过程。这个模块是全自动完成的并且利用一些NUMECA工具。此外，附加模块FINE/Design2D这些工具联系起来，可以进行叶片重新设计，改善叶片表面压力分布，关于这些详见第13章。这个模块假设流动是轴对称的，并且流面形状和厚度也由用户提供或由参数自动生成(利用根部和顶部边界)。几何输入数据必须由用户提供：1、 流面及叶片这个流面上的截面或2、 完整的叶片轮廓及端壁本模块由网格自动生成与NS湍流方程组成。在下一节讲述这个跨叶片截面模块的界面及对用户的建议。12-4节讲述自动生成网格的理论和求解方程。12-5节讲

2、述几何数据和输出结果。12-6讲述实例。12-2跨叶片截面模块的界面在FINE/Design2D界面之下运行跨叶片截面模块，这些可以高速，简单，交互式求解。所有参数可以在用户界面中选取，并自动创建输入文件及求解。监视工具，MonitorTurbo,可以在计算中和计算后检查收敛情况及结果。它可以实时查看叶片表面压力分布的收敛过程及叶片几何形状。结果分析利用 NUMECA CFViwe 后处理工具进行，自动进入跨叶片截面模式。几何数据以ASCII输入文件列出，但是求解参数定义及边界条件在这个界面中列出。这个截面的描述由FINE/Design2D界面中的菜单创建。更详细的说明见12-5.12.2.1

3、开始新的或打开现存S1面计算在开始界面下，Project Selection窗口允许创建新工程或打开现存工程。对于创建新的跨叶片截面工程，按如下操作：1、 单击按扭Create a New Project2、 选取工程保存路径及输入文件名3、 关闭Grid File Selection窗口，Design 2D不需要输入网格文件4、 进入S1流面模块，菜单Modules/Design 2D如果要打开现存工程，在Project Selection窗口中单击Open an Existing Project按扭，并在File chooser窗口中选取一个文件。最近使用过的文件在最近工程列表中列出。如果

4、所选取的文件是以Design2D模式保存的，则FINE界面自动转到这个相应的模块，显示界面如图12.2.1-1所示。FINE/Design2D界面如同FINE/Trubo界面一样，包括菜单，工具栏，计算设置与参数区域。在菜单中同样也有一个Modules项，可以快速转到其它模块。Design2D模块的图标栏仅包括2D计算内容。界面左侧的参数列表也是与2D计算一致的。这一项的大多数内容与FINE/Turbo工程是相似的。之间的差别仅在于：l 在Flow Model页：Design2D模块不能进行非定常计算。l Boundary Conditions页的说明见12-2.3l Blade-to-bla

5、de data页的说明见12-2.2l Initial Solution页的说明见12-2.512.2.2 跨叶片截面数据Blade-to-blade data页(上图左侧倒数第三行展开)出现三个分项：ll 网格生成参数，l 反问题：仅需2D设计参数，不需要分析。详见第13章。叶片轮廓和子午流面(或端壁边界)是由一些连续的点的坐标构成的ASCII文件输入，详细说明见12-5节所有的几何数据都有一个相同的长度单位(可以任意选择)。这个长度单位设置是在Blade geometry页Blade geometry项进行选择。12-2.2.1 叶片几何定义输入数据类型有三种形式可用：-2D：仅对轴流适用

6、，且流面为圆柱面(半径为常数)。这个叶片几何模型是在这个流面上。-Q3D：The axisymmetric streamsurface is given, and the blade geometry is specified on the given surface.-3D: 这个完全的3D几何模型是由连续的叶片截面组成（最少数2个截面，次序是从根部到顶部）。a) 流面数据图12.2.2-2显示了输入参数区域。对不同的输入文件类型需要输入文件：l 2D 或 Q3D：流面定义文件l 3D ：端壁文件，两个文件，根部曲线和顶部曲线输入文件类型为3D则需要附加参数：l 几何分区或用户定义流面 Ge

7、ometrical division 或 streamtube providedl 如果是Geometrical division：展向位置为0-1之间l 如果是streamtube provided：输入两个文件(两条直线)In both cases the module automatically calculates the intersection of the 3D blade with the considered streamsurface。对3D类型，推荐展向边界网格厚度为1%In a 2D，Q3D or 3D case with streamtube provided, th

8、e blockage ratio permits to scale the streamtube thickness distribution.无论哪种输入数据类型，都必须设置 NB of points along meridional trace.这个值默认为300.b) 叶片几何数据在Blade Geometry页，依据输入数据类型，叶片几何定义需要如下参数：l 叶片几何数据文件：-2D或Q3D：吸力面与压力面文件-3D：展向截面数及每个截面的位置，每个截面的压力侧与吸力侧文件l 前缘钝角边：如果选取，则前缘钝角处理l 尾缘钝角边：如果选取，则尾缘钝角处理l 分流叶片：(我对这个没兴趣)l

9、 长度单位：用户输入数据所采用的单位，为米与所采用单位之比(如果输入的数据单位为毫米，则这个比值为0.001)。全部几何数据必须采用相同单位，即用户所选取的这个单位。lc) 叶片数据编辑选取Blade data edition将修改前缘或与尾缘数据。12-2.2.2网格生成参数如图12.2.2-5这一项必须由用户修改。系统默认值对大多数网格是适用的。可用的网格生成参数有：l H- or I网格：用户可选是否包含上游和/或下游网格。l 生成方式：全自动或半自动-对全自动模式，用户仅需指定吸力侧网格点数。其它区域的网格数按光顺需要计算出。-对半自动模式，用户指定上游或下游的网格倾斜角度，以及每个区

10、域的网格数。l 圆周方向网格数:这个数最好是8的倍数，这样有四重网格可用。l 周期性边界的类型：选取直线或曲线。默认为周期性边界线为曲线。l 在下面的三页中的参数影响了网格的质量和分布a) clustering第一个参数clustering coefficient in streamwise direction输入值范围为从0.0到1.0。如果值为1.0则网格均匀分布，如果为0.0则所有网格点集中在边缘。第二项为定义Euler网格生成还是N-S网格生成。如果选取了N-S网格生成还必须设置以下参数：l 第一层网格单元的尺寸(单位：米)l pitchwise 方向的网格数l 从叶片边缘到进口或出口

11、的clustering 类型，constant 或 decreasingb) smoothing椭圆形光顺可改善网格质量。用户可选clustering control 与orthoganallty control，并设置光顺次数Number of sweep.c) local pre-smoothing局部光顺用于前/尾缘半径的网格。12-2.2.3 反问题设计菜单这个菜单为反向设计模块。当进行流动分析时，这项菜单的内容不需要指定。对于反问题模块的详细说明见第13章。12-2.3 边界条件在Boundary conditions 这一页允许定义进出口边界及旋转速度。这个Rotational s

12、peed单位必须是RPM(revolutions per minute每分钟转数，旋转方向的正向与定义的正向相同，参看图12.5.1-12)12-2.3.1 进口边界条件以下为可用的进口边界条件：l 绝对速度方向(弧度)+绝对总量(Pa ,K)l 相对速度方向(弧度)+相对总量(Pa, K)l 绝对速度分量V(m/s)+绝对总量(Pa, K)l 相对速度分量V(m/s)+相对总量(Pa, K)l 质量流量(kg/s)+绝对速度方向(弧度)+静温(K)l 质量流量(kg/s)+相对速度方向(弧度)+静温(K)l 质量流量(kg/s)+ 绝对速度分量V(m/s)+静温(K)l 质量流量(kg/s)

13、+ 相对速度分量V(m/s)+静温(K)用户选取其中一种方式设置进口条件。12-2.3.2 出口边界条件有三种出口边界条件可用：l 出口静压(Pa)l 质量流量(kg/s)(+出口静压用于初始计算)l 超音速出口(出口静压由初始计算)零阶和一阶外推法可用，同时可用专家模式(默认为零阶)。设置质量流量为边界条件对于反问题方法则是比较稳定的。因此在分析过程中也推荐用质量流量作为边界条件。比较理想的设置是质量流量进口和压力出口，或是进口的总量与出口质量流量。12-2.4 数学模型通常Numerical Model参数的默认值是适当的。In the blade-to-blade module it i

14、s not possible to perform the calculations on a coarser mesh level as it is the case with 3D projects.12-2.5 初始求解菜单12-2.6 输出参数这个菜单允许选取一些输出量。这个模块自动创建用于CFView工程的3D与2D的输出文件。关于这个菜单命令详见12-5.2.12-2.7 控制变量在这个Control Variables 页用户需要选择流动分析或反问题设计模式。如图12.2.2-10.所有的相关参数见第15章FINE/Turbo计算。12-2.8 跨叶片截面流动分析一旦 .b2b

15、和 .b2b.inp 输入文件被正确创建，就可以开始流动分析，在Solver菜单单击Start按扭。重新开始以前的计算也同样是在Solver菜单按Start按扭，并选取一个初始求解文件。不能在粗网格中进行计算，重新计算也不可以多重网格设置。12-3 专家参数在专家模式下Control Variables 页，连续的专家参数可用：ITERZ ：迭代步数，在这个范围内不进行湍流计算。通常设为200.IATERZ : 这个转换与前一次计算有关。通常设为1(默认值)。12-4 理论这一节主要讲述网格生成工具和跨叶片截面流动求解。12-4.1 网格生成这个网格自动生成工具可创建H(周期性) 或 I(非周

16、期性) 形网格，按用户提供的参数。相反多数quasi-3D方法在轴流面上生成2D网格，这个网格相当于3D网格在展向只有一个单元。这个求解也主要是采用3D流动求解。此外12-4.2 求解器12-5 跨叶片截面模块的文件格式12-5.1 输入文件B2B运行需要三种输入文件：l 几何输入文件，叶片几何定义l 求解输入文件(.run),求解参数，流体，流动条件定义l 反问题设计输入文件允许定义反问题这个由用户创建的几何输入文件与AutoGrid输入文件的格式相同。求解和反问题设计输入文件是由用户在FINE/Design2D环境中设置参数和流动条件后由程序自动生成的。在下一节讲述几何输入文件。12-5.

17、1.1 几何输入文件叶片几何数据和流面(或端壁)是由ASCII文件中连续的点的坐标构成。三种格式的几何输入文件可用：l 2D类型：这个设置仅对完全的轴向网格，流面为圆柱(半径为常数)。The blade geometry is then given on the cascade plane(cylinder developed to generate the cascade).l Q3D类型：The axisymmetric streamsurface is given(including its thickness), and the blade geometry is specified

18、on the given surface.l 3D类型：3D叶片模型由一些叶片截面构成(最少2个，按由根部到顶部的顺序)。这个根部和顶部端壁也同样是单独列出。这个流面可由用户指出(轮缘和轮毂)，或自动计算得出(几何分区).在下一步，用户仅需指出展向位置及厚度。所有的几何数据都应用相同的长度单位(用户所选取的).a) Streamsurface Data:对2D与Q3D类型，在ASCII文件中必须指定流面和厚度。对3D类型，由两个文件分别定义根部和顶部端壁。Streamsurface 数据文件构成如下：第一行：输入文件标题第二行：坐标系统： ZRB(用于2D与Q3D类型) ZR，RZ 或XYZ(

19、用于3D类型)第三行：N 数据点的数量第四行到第3+N行：点的坐标2或3列(按进口到出口排列)第4+N行： 行标记(可选)第5+N行：Zin, Zout, Rin,Rout,指定网格进出口位置(可选)(X,Y,Z)为直角坐标(Z为旋转轴方向)，R为径向，B为流面厚度网格的进出口位置按以下方式确定：l 按流面数据(4+N行和5+N行)l 利用FINE输出参数B2BLIMb) 叶片模型数据12-5.2 输出文件12-5.2.1 几何文件12-5.2.2 变量文件12-5.2.3 数字控制文件12-6 网格验证12-6.1 VKI-LS82 叶片12-6.2 ARL SL 19 亚音速压缩机叶片第十

20、三章 Design 2D13-1 绪言本章描述了FINE/Design 2D模块进行透平反问题设计。这个模块是与第十二章跨叶处截面相关的，增加了网格生成和流动求解，2D模块包含了叶片重新设计的反问题方法。这种方法允许按叶片表面压力分布重新设计叶片型线。使用本方法需要跨叶片截面流动的许可。FINE/Design 2D单独使用一些NUMECA工具。NUMECA计划扩充FINE/Design 2D的功能，应用领域和设计方法。这个第一个版本对于透平叶处的反问题设计有局限，因此综合了跨叶片截面方法在这个界面。分析和设计这两个模块基于一个假设，即流动是轴对称的。这个几何输入是必须的，对于反问题设计则不需要

21、分析。(参看第12章)关于跨叶片截面方法的更多信息，见第12章。在下面的章节中介绍了使用反问题设计方法的注意事项。关于反问题的理论分析见第13-3节。Design 2D模块的几何数据输入和结果输出见第13-4节。最后在第13-5节以实例验证了这个模块的功能和作用。13-2 反问题设计这个反问题设计的相关描术见第十二章。全部的参数都可在界面中选择设置，自动创建文件及分析。图形控制，Monitor Turbo,可以实时监测收敛过程，如压力沿叶片表面分布，叶片生成，同样也可与压力分布相适应。这个结果可以被NUMECA CFVivew后处理。等等这一节讲述下面三个问题：l 开始一个新的或打开型线工程l

22、 计算及相关设置l 分析跨叶片截面流动13-2.1 开始这个模块运行首先要启动FINE界面，然后菜单Modules/Design 2D。启动FINE后，在Project Selection窗口，允许Create a New Project或Open an Existing Project。创建一个新跨叶片截面工程按如下操作：1、 单击Create a New Project按扭2、 选择路径，并输入文件名，确定3、 关闭Grid File Selection窗口，没有网格文件与当前工程关联。4、 转换到Design 2D模式，菜单Modules/Design 2D如果要打开一个已存在的Desi

23、gn 2D工程，单击按扭Open an Existing Project在Project Selection窗口中，并在File Chooser窗口中选取文件。如果是最近计算过的工程也可以在最近工程列表中选取。如果所选工程是在Design2D中保存的，则FINE界面自动转到这个相应的模式。Design 2D模块界面的介绍见第12章。这里仅讲述跨叶片截面分析计算，相应的操作选项页面为Inverse design.下面各部分讲述相关的设置。13-2.2创建反问题设计的输入文件反问题设计应该是从初始几何模型的流动分析开始。当跨叶片截面计算收敛后进行分析，用户可用反问题方法修改叶片。为这一目的必须准备

24、两个文件”,”，这两个文件是在分析过程中添加的。13-2.2.1 输入文件”这个文件”是在FINE/D应该选择三个可能的反问题公式中的一个：l Classical formulation: 利用压力在叶片周围的分布修改叶片的吸力侧和压力侧曲线。如果要使用这一公式，则这两条曲线必须是封闭的。l Mixed formulation: 仅对吸力侧修改，按吸力侧压力分布。这个方法也同样要求曲线是封闭的，至于压力侧则按吸力侧修改。l Loading formulation: 按从前缘到尾缘的加载分布修改叶片的曲率直线。因为叶片厚度分布是常数？，所以不要求轮廓线封闭。这个反问题设计的目标计算可按等熵马赫数

25、或压力系数(pressure coefficient Cp)进行.这个等熵马赫数的计算对于可压缩流动是限制的。允许控制前缘和尾缘。这个反问题方法允许保持一些局部不变，如前缘和尾缘区域。这些固定区域的压力侧和吸力侧的点的数量由用户指定，默认为1。无论如何，冻结前缘和尾缘区域是被推荐的。如果初始叶片流场分析中发现了尾缘的流动分离，极力推荐冻结尾缘区域。有两个松驰因子可用：l 第一个松驰因子(默认值为1.0) 两个文件”.req”和”.vel”。这个松驰因子不经常使用。l 第二个松驰因子(默认值为0.05)，对几何模型修改。这个值从0.01到0.1这间改变。高的值0.1因稳定性不推荐使用，相反，低的

26、值0.01通常不能进行反问题处理。这两个松驰因子默认值分别为1.0和0.05.最终，下列几何模型约束是可用的：l 几何模型轮廓封闭永远可用，除非：-叶片厚度是自动约束的(“loading” formulation)-叶片尾缘是钝角，并且trailing edge thickness does not need to be maintained.l Constant stagger angle 不经常使用l Constant blade chord: 这个约束允许保持叶片的弦长,这个选项不应该与constant radius 共同使用，因为这两个设置是相矛盾的。l Constant radius

27、: 这个设置推荐用于径流透平，允许保持子午面前/尾缘位置不变。这些点被转移到圆周方向。这个设置不推用于在半径尾缘。l Reference point: 选取前缘为固定参考点(默认值)，或尾缘(需选相应按扭)。13-2.2.2 输入文件 ”projectname.b2b.req”这个”projectname.b2b.req”文件如同”*.vel”和”*.load”输出文件一样，也是由B to B分析模块自动生成。 如同由初始模型流动分析的结果文件projectname.b2 b.vel 或projectname.b2b.load 一样，利用MonitorTurbo生成。利用MonitorTurb

28、o创建这个文件，应按以下次序：1、 启动 MonitorTurbo2、 选取菜单 Loading diagram3、 打开或文件4、 选取Mis或Cp分布(in case the inlet total conditions are not imposed by the boundary conditions, the inverse problem has to be formulated in terms of Cp)。这个Mis和Cp的定义见公式Eq.13-1和Eq.13-2，第13-4节。5、 激活Markers以区分离散点6、 激活Edit a curve按扭7、 选取吸力侧或压力侧曲线，把游标放在相应的区域，按<s>键选取。被激活的曲线颜色改变8、 利用鼠标中键修改这条曲线：a. Choose the lef