第七章 网格设计 Grid Design 网格是流场计算的基础 It is the basis of calculating flow field 7-1 几何方法构筑叶栅通道网格

1. 第七章 网格设计 • Grid Design • 网格是流场计算的基础 • It is the basis of calculating flow field • 7-1 几何方法构筑叶栅通道网格 • Geometry method to construct grids • 两种方法： Two Methods: 1.物理空间构筑曲线网格，变换到计算空间的正交网格，计算在 正交网格中进行。To construct curve meshes in physic space，then transfer to orthogonal grid to compute in computational space • 2.直接在物理平面内构筑网格并求解，此方法比较容易变化网格密度以适应参数梯度。Directly construct grids in physic plane and solve the equations in physic plane， this method makes it easy to change density of grids.

2. 网格边界分别平行于求解域边界或与边界相适应。网格边界分别平行于求解域边界或与边界相适应。 • The grid boundary parallel to the boundary of flow field.

3. 通道内X方向间隔相等 ∆X=const Keep constant grid gap inside channel • 通道进口前 In front of channel • 通道出口后 Behind the channel • Y方向网格间距常数（Euler方程计算） The grid gaps in Y direction kept constant • 几何关系如下 • Geometry relationship 按指数规律伸展∆X ∆X changes in power law

4. 变换后的方程 Transferred equation • 变换关系导数Transform derivatives

5. 变换后在 平面内的控制方程为 The equation in transformed plane 二.局部加密的叶栅通道网格 The mesh in cascade tunnel with local refining • 翼型或叶栅头部附近加密 Refining at leading edge of airfoil and cascade • 叶栅通道出现激波时，在激波位置加密 Refining when shock appears in cascade channel

6. 壁面附近边界层 Refining in boundary layer of the wall • 例：example • 通道内X方向均匀 Uniform grids in X direction of cascade tunnel 通道前与第一个网格宽度 ,其余采用 • Keep the width of first mesh as same as that of in channel mesh on others grids

7. i=2，3，4,I • 下游类似 The same way used for down-stream meshes • Y方向网格划分Grids in y direction 即y方向按指数规律变化 为ACG边（下） 为上边界 That is power law used in y direction

8. 一般要求第一个间隔<0.015，可调节J数来保证 The first grid gap is 0.015,it can be ensure by adjusting the total number of grid J • 在计算平面内 ，各个相邻网格点的间距 • in computational plane ,the grid gap is used equal

9. §7-2 几何方法构造c形网格 • Construction of C-Type mesh using geometric method • 对外流问题（如NACA-0012翼型绕流）通常可构筑C形网格 • For out flow question , usually construct c type grids y • 将前缘至后缘流向划分50个点，按2次曲线（抛物线分布） • Discrete chord wise into 50 points , using 2nd parabolic curve L x

10. 将坐标向左平移0.00125，则NACA-0012的翼向y坐标 • Moving coordinates 0.00125 in –direction, then the y coordinates can be • determined by • 用两组抛物线构筑c型网格 • Construct the mesh using two serial of parabolic curves

11. §7-3 保角变化法 Conformal transformation 一、基本概念 basic concept （一）在(x,y)平面内两族网线若存在 on (x,y) plane, the serial carves as follows 则两族网正交 then the two curves are orthorgal

12. (二)：用复数代表(x,y)平面 Define the (x,y) plane with a complex variety • 两族曲线用复函数表示为 • two curves can be denotes with a complex function • (三)正交坐标系 称为 • The orthogonality coordinate expressed as • 两平面的关系 the relationship between two coordinates can be expressed as

13. 在 平面内两线族 in plane the curves are • 彼此关系可写成 The relationship between physical plane and computational plane is follows

14. 在z平面和 平面内，点与点相互对应 in physic plane z and computational plane the points responds each other • 两平面内的正交曲线相互对应 the orthorgal curves in two planes are corresponding • 相应的交点也相互对应 the intersect points are corresponding

15. 二：保角变换方法生成网格 • Grid generation with conformal transformation (isogonality) • 一、例一:复平面内正交直线族 • Ex1, the orthogonality curves in complex domain • 代表与坐标轴平行的正交网格的两个函数 Two function denote the grids parallel with coordinates • 另一复平面 ，两平面之间的关系 • In complex domain , the relationship between two domains are

16. 两平面间的网格族是 • The grids in two domains are 双曲线 Hyperbolic curves 双曲线 Hyperbolic curves 取0，0.1，0.2，0.3……，1.0时对应的网格 Where is specified as 0，0.1，0.2，0.3……，1.0,the grids are

17. （二）例二：复平面内正交线族的函数 Ex2, orthogoral curves in complex plane

18. 变换关系： • Transformation relationship

19. 网格线族为 • The grids are • 取0，0.1，0.2，0.3……，1.0时对应的网格线如下图 • When are given as 0，0.1，0.2，0.3… the grids in (x,y) and plane is follow

20. 利用保角变换可以生成任意封闭曲体外部的网格利用保角变换可以生成任意封闭曲体外部的网格 • The grids around an arbitrary close body • 通常用于变换的方程可以写成 • The relation used for conformal transformation can be written as

21. 7-4 生成网格的微分方程The partial differential equation method for mesh generation 适/贴体网格生成可以用一个微分方程边值问题实现 Body fitted grids can be gained using the solution of a partial difference EQ TTM方法（Thompson, Thames and Martin）方法: 利用求解微分方程方法生成空间区域网格的方法 TTM method is the method to generate grids using PDE 要求 requirements 边界形状已知，可以正确描述 Boundary as known, can be described, correct

22. 网格密度可调，有伸展性（可变网格密度） the density of grids can be adjusted and extended • 附面层：尺度为 ~ Boundary layer :size • 分离流：尺度为弦长C Separated flow :size C • 尾迹区：尺度为 ，弦长C Wake flow :size C • 无粘区（势流）：尺度为 Inviscous flow:size

23. 正交性：可以减少差分方程所依赖的网格点数 Orthorgality :may decrease dependent grid of FDE • 保形性：物理平面与计算平面方程相同 Conformality: the EQ in physic plane and computational pla are same (一)TTM方法建立物理平面与计算平面的关系 TTM method is used to found the relationship between physic plane and computational plane • 可用微分方程表示上述变换 The transformation can be expressed as follow

24. 正确表述外边界 • Define the boundary correctly

25. 在计算平面内，为了差分方程简便，一般用等间隔网格在计算平面内，为了差分方程简便，一般用等间隔网格 • To be convenient for FD，the grids in computational domain is equal gap • 在 平面内构造正交网格，再在 平面内生成 相应的点坐标，即求逆变换方程的解 Normally, construct the orthogonal grids in plane first, and then create the corresponding points in plane, that is to find the inverse transformation • 变换关系

26. 导数关系：

27. 导数关系：

28. 且有 and then • 同理可得 can be obtained in the same way

29. 最后可得逆变换方程 Finally the inverse transformation equation obtained • 求解此方程组可得离散表达形式 To solve above equation and obtain transformation relationship in discreted form

30. 在 平面 和 代表 和 代表 denote line and In plane, , • 在 平面对应点 • 边界点给出后就可以得到内点上的各点

31. 正交性和保形性（采用Laplace变换） Orthogonality and confor • 由坐标变换关系可知 From coordinate transformation • 采用Laplace变换，则 Using Laplace transformation, then

32. 反变换为 • inverse transformation is • 没有混合导数项，具有保形性和正交性 No mix derivatives existed

33. 采用正交，保形变换后，全速势方程的形式不变采用正交，保形变换后，全速势方程的形式不变 Using orthogonal, conformal transformation , the fully potential equation unchanged. • 拉普拉斯和泊松方程的变换关系式 Laplace vs. poisson transformation • 拉普拉斯变换 • Laplace transformation

34. Laplace变换是单值变换， 变换前后的点一一对应 The Laplace transformation is a unique value transformation • 物理平面是单连域，则计算平面也是 The single connected field physical plane is still a single connected field in computational plane • 例：机翼（翼型） • Wing as an example • 极值原理：参变量的极大、极小值必定在边界上 • Extremum principle: The maximum and minimum is on boudary

35. ＞ (a,b) b 2 3 4 1 0 a

36. （二）泊松方程 • poisson equation • 可以控制网格密度 the grid density can be controlled 解为 The solution is u u u 0 x x x 0 0 1 P=2>0向下加密densed to lower surface P=0均匀uniform P=-2<0向上加密densed to up surface

37. 用于网格生成的poisson方程 • Poisson equation using for mesh generation • 反变换方程 • The inverse transformation

38. P（ξ，η）和Q（ξ，η）可以调整网格密度 • P（ξ，η）and Q（ξ，η）may adjust the grid density • 其中 表示要求网格向 点靠拢， 是调整量。n，m表示要向 靠拢的网格数量 • Where denotes the grids will approach to ， are adjust parameters.n，m denotes the number grids to be closed to point.

39. 另一种网格加密办法是由Thomas和Middle coff 提出的 • Other method for refining grids is gained by Thomas and Middle coff

40. 7-5 代数法和混合法 一、代数法 采用几何剖分方法，利用代数运算生成计算区域网格，无需求解微分方程。 • 等比网格： • 等差网格： • 指数网格：

41. 二、混合法 先利用方法生成一个方向（平面内的代数网格），再利用TTM方法生成另一个方向（平面）内的网格，最后将所生成的网格联接成一个整体网格 对收扩喷管：横截面用TTM方法，而轴向横截面用代数法。

42. 7-6 动网格设计 非定常流中，边界运动的处理方法。（每一时间步需要重新生成网格）。 • 考虑计算流场边界和网格边界的运动，采用积分分解式的质量守恒律。 ：流体流动速度矢 ：求解域周边界移动速度矢

43. 将求解域划分成J个网格，则空间离散后的方程为将求解域划分成J个网格，则空间离散后的方程为 • 网格体积变化 • 每个网格来说 ——几何守恒律。描述几何性质，具有积分形式的守恒律表达式。 • 在物理空间求解时，需联立上述系数方程

44. 以密度场为例 • ---平均密度 • 在计算空间求解（比较方便） • 计算空间的网格不随时间变化，可以用积分形式的守恒方程也可以用微分形式的守恒方程 • 计算空间的守恒方程是经过变换的方程，如连方程

45. 坐标系之间的面积体积关系 • 三维Jacobian行列式 • 流体质量 • 物理空间 • 计算空间 • 两者之间的关系

46. 物理空间微分形式的流体动力学方程为 • 计算空间微分形式的动力学方程为 • 每一时间步都要计算Jacobian行列式J的值 • 计算求解式和边界分别为

47. 物理域和网格边界分别为 • 边界移动速度 • 体积变化守恒律 • （运动速度）

48. 即 • 在计算平面内网格不随时间变化，故 • 方程左侧空间导数

49. 故 由于

50. 同理 • 由三维空间坐标变换量关系式得