热物理过程的数值模拟

《热物理过程的数值模拟》由会员分享，可在线阅读，更多相关《热物理过程的数值模拟（25页珍藏版）》请在装配图网上搜索。

1、四、非线笥问题迭代式解法的收敛性每一层次上满足迭代法求解的收敛条件+相邻次间代数方程的系数变化不太大(亦即未知量 的变化不太大一多数情形下非线性问题迭代式解法是可以收敛的)。使相邻两层次间未知量变化不太大的措施：1、欠松弛迭代常用逐次欠弛线迭法(SLUR)： 一组临时系数下逐线迭代求解+对所得的解 施以欠松弛，再用欠松弛后的解去计算新的系数，常数，以进入下一层次的迭代。实施：常把欠松弛处理纳入迭代过程，而不是在一个层次迭代完成后再行欠松弛。Za bt bt (n+1) = t (n) + W (nn t (n)ppappaa(p)t(n+i)=Za b +1 b + b + (1 W)pt(n

2、)W pn nW pa t(n+i) = Y a bt b + bp pn na = a W, b = b + (1 W)(a W)t(n)，用交替方向线迭代法求解这一方程，就实现了 SLUR p pp p的迭代求解。为一般化起见，上式中t b上没有标以迭代层次的符号J, GS时不相同)。n2、采用拟非稳态法 前面已指出，稳态问题的迭代解法与非稳态问题的步进法十分相似。对于非线性稳态问题，从代数方程的一组临时系数进入到另一组临时系数亦好象非稳态问题前进了一个时间层，非稳态 问题的物理特性：系数热惯性越大(ao = pcAv/ At T)，温度变化越慢，仿此，对稳态非线性 p问题，可在离散方程中

3、加入拟非稳态项，以减小未知量托两个层次间的变化，即由(Za b S AV)t(n+1) = Za bt b + b n (Za b S AV + a。)t(n+1)=Za bt b + b + aot(n)nppn nnpp pn np pZa bt b + b + aot (n)t (n+1) = n np ppZ a b S A V + a onpp一直进行到t ,t b收敛，虚拟时间步At的大小通过计算实践确定。pn3、采用 Jacobi 点迭代法 中止迭代的判据(该层次迭代)除前述变化率判据外，还可以规定迭代的轮数，例如规定进行4-6次ADI线迭代就结束该层次上的计算。此时，用收敛速度

4、低的丁迭代也就起到了欠松弛的 作用。五、迭代法的收敛速度1、收敛速度 对给定的代数方程组(包括是临时系数的情形)，采用不同的迭代方法求解时，使一定的初始 误差缩小成Q倍所需要的迭代轮数K是不相的。 1因为 G 是对称的，所以|G|(GtG) = l：P(G2)= p(g)2 所以/|e(0)|2 2e(o) | =P (G) k |e(o)|2 2 2 (P (G) k, l a l a /1 P (G) = _l a /(I P (G)a, p (G)均 l a / R or R l a / knn即所需迭代轮数与收敛速度成反比，收敛速度又与谱半径成反比，收敛速度愈快，迭代轮数 愈少。注意：

5、不同的迭代方法每进行一轮迭代所需的运算次数不同，最终所需的计算时间的多少取 决于迭代轮数及每一轮迭代所需的时间。2、收敛性的定性分析为什么不同的迭代方法的收敛速度不同，亦即为达到满足一定精度要求所需的迭代轮数不yxAttTB/ tRB同？以二维常物性、无内热源、稳态导热问题来进行讨论。dx2dy21B、C a t = a t + a t + a t + a tp p E E ww N N S S迭代法需要假定一个初场，例如假定一个均场， 从微分方程为看，均匀是其一个解，但却不是所研究 问题的解，为什么？因为它虽然满足内部节点上的离 散方程，却不满足与边界有关的节点的离散方程（图 中红点），即不

6、满足边界条件。所以迭代法的实质是要 通过迭代，尽快建立起与边界条件相适应的e变量场, 关键：必须使B、C的影响迅速传入计算区域内部，以改进节点e变量值，尽快与B、C相应，B、C 的影响传入愈快，逼近真解就愈快，收敛就越快！B、C的影响传入计算区域内部的快慢与哪些因素有关？（1）与迭代方法有关J 迭代：节点温度的更新均用上轮迭代所得的“旧”值来计算，所以完成一轮迭代后， B、 C的影响只能传入与边界相邻的一批节点上，即仅可传入一个网格，且扫描方向与收敛快慢无关。 要在以后各轮迭代中，B、C的影响才由这些节点逐步向内渗透，所以收敛慢。GS 迭代：假设从左向右扫描，则每做完一轮迭代，左边界和下边界的

7、影响传遍全区域，而 右边界的影响只能传入一个网格，且收敛速度受迭代扫描方向的影响。线迭代：GS线迭代。自左一右扫描，完成一轮迭代不仅左边界的影响逐步传入，而且在每 一列的直接求解中，上、下边界的影响全部传入到该列的各节点上，即一轮迭代使左、上、下边 界影响传入全区域，但右边界影响仍仅传入一个网格。ADI： 轮迭代包括一次逐行、一次逐列的扫描；所以在每一轮迭代后所有边界的影响均传 入计算区域内部，从而加快了收敛速度。收敛速度的比较，正方形区域，IB、C，Laplace方程五点格式，均匀网格步长为h。2）与边界条件的性质有关迭代方法点迭代线迭代Jacobih2/2h2Gouss-Seidelh22

8、h2SOR2h2屮2 hIB、C规定了边界节点的温度，影响直接传入计算区域内部； 九dt3B、C规定了环境温度及定向点位置，二t -1，对边界温度的限定程度比IB、C时l dxg弱，所以对内部的影响也较弱；或将tf视为外部温度，其对计算区域内部的影响被外部换热热阻 削弱，而IB、C可视为a fg或外部热阻T 0的极限情况，故3B、C的影响比IB、C时弱。2B、C仅规定了壁面的钭率，壁温完全不确定，对内部节点温度值的改进提供的信息最少， 收敛最慢。可见，为了提高代数方程迭代解法的收敛速度，应力求使边界条件的影响迅速传入计算区域 内部，措施： 增加迭代解法中直接解法的成分，从点迭代一线迭代一ADI

9、； 适当选择扫描的始边，多以IB、C或3B、C的边界为始边，少以2B、C （尤其是绝热边 界）为扫描始边。5-4 不规则区域的处理网格生成技术如何对不规则区域进行有效的处理，以便于进行传热与流动过程的数值模拟，是近年来计算 传热研究中的一个重要课题。以上讨论的传热过程大都发生在规则而简单的区域中，但许多实际的热传递现象是在不规则的区域中进行的，例如：套片管中肋片的传热渐扩通道中的流动与换热流体外掠管束以上四种情形中的流动与换热不是直角坐标、圆柱轴对称坐标或极坐标所能方便地予以描述。 虽然有限元法在处理不规则边界方面显示了极大的优越性，但就流动与传热而言，在计算技 巧与方法方面，有限差分法都比有

10、限元法成熟。用有限差分法处理这类问题的方法可以归纳为以下几种。1、采用阶梯形边界（网格） 用阶梯形边界近似代替四分之一圆弧边界。阶梯形边界（网格） 是采用有限差分法计算不规则区域的最普通的方法。缺点：程序缺少通用性；曲面边界上网格必须划得比较细密（否 则会引起较大误差）。2、采用区域扩充法当计算区域的边界不规则程度不很严重时，可以采用区域扩充法，把计算区域扩充为直角坐标，流动：把计算区域扩充到图中虚线所示的整个圆形通道，从而可以应用圆柱轴对称坐标系中 的控制方程加以描述。如何保证原来的情况不变？孔板区视为粘性无限大的“流体”，而其余区域的流体粘性值就等 于真实值，边界上流速赋为零，计算中零边值

11、将迅速传到孔板区域内，有效地模拟了孔板的存在。传热：对三种不同的边界条件，具体的处理方式不同（1）均匀壁温边界条件令扩充区域中的导热系数为无限大，而扩充后的区域边界温度则等于已知值。（2）绝热的边界条件令扩充区域中的材料导热系数为零即可实现此条件（3）均匀热流边界条件可应用附加源项法来实现，真实边界上均匀热流可以 附加源项的形式置于与真实边界相邻的控制容积中去，而 扩充区域则处于绝热状态。周期性二维渐扩、渐缩通道中的换热，倾斜的边界上 作用有均匀的热流。采用左下所示阶梯形网格，并把计算 区域扩充到一个长方形，以便利用直角坐标系求解。看一个网格单元的放大后的情形，P控制容积的附加 源项为S 二

12、qL / AVad of abcdLe实际边界与控制容积P的两条边界相交部分的长度;AV 控制容积P的体积abcd扩充区域九=0,则控容P中的附加源项S不会向扩充区域传递，从而实现了实附边界上 扩ad的均匀热流加热条件。（4）外部对流换热边界条件，l -1f根据附加源项法，此时P控制容积的两个附加源项为Sc,adtLf. of-1/a +6 /九 AVp,ad1Le1/a +6 /九 AVdcb126网格节点P到实际边界的距离九扩二0区域扩充法的优点：可以用按规则区域编制的通用程序来计算非 规则区域的问题；易于实施。缺点：浪费一些计算机内存及计算时间。3、采用三角形网格 外节点法对于不规则区域

13、中的导热问题，采用三角形网格可以得 到比较满意的结果。从不规则区域的三角形网格中划出围绕节点P的多边形来分析。 确定节点P所代表的控制容积：三角形外心作为控容的顶点，要求： 三角形为锐角三角形，以保证外心一定在三角形内；以三角形重心作为控容的顶点，对三角形形状无限制；外心为控制容积顶点，P控制容积如图所示，各部分平均导热系数分别为九，九，九，九，用控制容积能量 a b cf平衡法建立离散方程。P-1之间的热导c =九/L /L +九L /Lp1aa1p1b b1p1La1所以CP111=L cot P L = L cot Pb12 p1211cot P + 九 cot P =(九 cot P

14、+ 九 cot P )1 2 b 2 2 a 1 b 21 二一九(cot P + cot P )2122 p 11二-九2a常物性时，Cp1其它各节点(2-6)与P节点之间的热导亦按上式计算，而只需把仲P-改换成相应顶角即可,P 控制容积的热容：变物性时，图示各部阴影部分的热容量之和Y (PcAV ) p,i ii由控制容积 P 的热平衡可得t -to乙(pc AV )-p p Ai常物性时，上式变为二工 C (t -1 ) + S (AV )pi ippt -to pcAV p p 二乙 C p Atpii以上二式右端温度值的时层未标出，它取决于所采用的格式。边界条件的处理：第一类边界条件

15、，无需专门处理； 第二、三类 B、 C 边界上的节点，需要仿上述方法建立补充方程 和图第二、三类边界上的节点P的控制容积1二一九 cot P2 c 3(t -1 ) + S(AV )i pp1C 二一九 cotP , Cp12 b1p 31C =(九 cot P + 九p 22 b 2 c3cot P )4边界加入热量则计入到P控制容积的源项中，对第三类边界条件,传入热量a (t -1 )L，其中t的时层取决于所采用的格式。f p adpWinslow 指出，当导热系数为常数时，三角形网格所建立的离散方 程的系数矩阵具有正定性，可用 SOR 求解；导热系数与温度有关时，为求迭代收敛，不宜采用S

16、OR,且两次迭代之间的C可用欠松弛。pi 根据实际问题的需要，可以采用三角形网格与矩形网格的组合结构 例如二维纵向肋片管。三角形网格的缺点：节点位置的确定、编号，节点间距的计算比较费时程序比较复杂。 内节点法：以每个三角形的外心为节点，其余同上。4、组合坐标法 二维平直一弯曲段组成的通道，可采用坐标系组合的方法，例如直角坐标与极坐标的组合 把该计算系统分成三个区域，在区域I、III中采用直角坐标系，在区域II中采用极坐标系。y速度U的控制容积三个区域中垂直于主流方向的坐标增量是一致的，dy二dr，但在主流方向上的坐标不一致，为 求得一致，在I、III区域中定义二x/R，则I、iIII区域中的也

17、是弧度，与II区中一致，从而可以 从入口截面开始就统一地连续计量坐标。在各区域中分别按相应的控制方程建立其离散方程，然后在整个区域内统一地求解代数方程。三个区域的两个分界线作为主控制容积的界面线，所以当采用交错网格时，主流方向速度的 控制容积将跨越两个区域，如图所示。该速控制容积的a ,a分别Ew按级坐标和直角坐标系计算，而a ,a可按并联处理：先分别按两NS区域中的半控制容积计算，后相加，源项亦仿此。5、采用特殊的正交坐标系三维正交曲线坐标系中，稳态对流一扩散问题的通用控制方程为h h h 6x1123(pu h h ) +1231ah h h ax1 2 3 2(pu h h ) +1a2

18、13h h h ax1 2 3 3(pu h h )312i a hh a、 i a hh a=（厂一13） +（厂一12） + Sh h h axh ax 2 h h h ax3h ax31 2 3121 2 33控制方程写成这种形式时，一切由曲线坐标，滞流模型引起的附加项都包括在S中去了，因 此， S 的表达式一般比较复杂。曲面边界平面通道中层流充分发展对流换热，设x3主为主流方向，且为直线坐标，再利用充 分发展条件，可得a (厂h aa h a)(厂 2) +(厂 +) = - Sh haxh axaxh ax1 21 1 1 2 2 2式中，r,c例于表中rCU3uh h dp / d

19、x1 23tU/Pra(h h pu t)/ax1 233在该流道截面的二维正交曲线坐标系中取出一个控制容积p来分析，x3方向控制容积厚度为AX3,将上面的控制方程在该控制容积上作积分，等式左端第一项为n e dh deI = Ax JnJe(r屮3 s W dxh11)dx dxdx1 21f h帥 =Ax fn (r -2)3 一 h dx11)dxh dx 211利用度规系数的性质，上式进一步变为e eI 仝 Ax rE p J3 e (6s )1 PEn hdx2r2e eJh dx 2 一W=Ax3re (eE eP)罟1 PEr (e e)(Asi)e W p/、w (6s )2

20、SP类似地，左端第二项积分为：11 Wn (eN -ep)計-r(e2 NP-e炸2 SP等号右端的积分为：c = (s + s e )hhc p p 1 2III = Ax3JJeCdxdx =(S + S e Jf2W1e h h dx dxc p p s W 1 212= ( S+s epp)AVp其中控容体积AV为:p_ (As ) + (As )2 e2 W (As ) + (As ) _ 2 一一_ 2 一将 I、II、III 代入积分式，整理得到Ax3AV =pb+SeSa+NeNa+wewa+EeEa=pepa中其1 7a = a + a + a + a S AV , b =

21、S AVp E w N S p pC p前面导出的三种典型坐标系中的二维导热的离散方程只不过是上式的特殊形式，从上式出 发可以编制正交坐标系中二维导热的通用程序不同坐标系中只是规定hi, h2的不同计算式即计算空间可。（对流动问题，采用原始变量法时，亦需应用交错网格） 对于扩散（导热）问题，当计算区域的边界正好与正交曲线坐标系的坐标相吻合时，采用相 当导热体的概念可以使这类问题的计算得到简化，无内热源时，能量方程为dhhhdt .dhhhdthhhdt ._（九 12 3） +（九 12 3）（九 12 3） = 0dxh 2dxdxh 2dx h 2dx11122333或y Q / hhh

22、dt .门乙（九） = 0dxh2 dxi=1 iiihh h如果把x , x及x3视为一新直角坐标系中的三个坐标，把九T 作为相当导热系数，记着 1 23h21九*，则上式化为iy3 ddt乙 （九*） = 0dx i dxi=1 ii这样，在原直角坐标系中的一个曲面物体就被转换成了新直角坐标系中的规则长方体，导热 系数由九变成了各向异性的相当导热系数九*，原直角坐标系所代表的空间称为物理空间，而新直 i角坐标系所代表的则称为计算空间，如图所示：x3V7x1 物理空间 x于是，该导热问题可先在计算空间中求解，由于它在计算空间中是规则区域，数值计算易于 进行，然后把计算空间中的计算结果传送到物

23、理空间中去，两个空间中点的对应关系由曲线坐标 的定义所规定。在由物理空间向计算空间转换时，边界条件亦应作相应的转换而成为计算空间长方体的边界 条件。第一类边界条件：把给定温度赋给计算空间中相应的点即可。 对第二、第三类边界条件，相应的模式列于表中类别 空间第二类边界条件第三类边界条件物理空间九一qQn n/、Qta(t t ) =X f BQnB计算空间Qthh h九 *=12 3 qi Qxhniihh h ,、1 2 3 a (t t ) hfiaQt=九*Bi QxiBx =const把物理空间中的曲线坐标系当作计算空 间中的直解坐标系，从而把物理空间中不规则 形状的求解区域变换成计算空

24、间中规则的长方 体（三维问题）或矩形（二维问题），这种做不xyx2qn .a ,tf仅适用于导热问题，也可用于对流换热问题，于是，有限差分法的一个主要弱点不易处理 不规则边界问题，可以说原则上已根本解决了。现在的问题是，现成的曲线坐标系的数目x1x2q; consta* = h h a2 1 3h h h=123 q =hh q h n 13n2毕竟是有限的，而不规则形状的物体则是千变万化的，不可能满足上述变换的需要。要使上述变 换思想付诸实现，必须发展一套方法，对于那些没有现成曲线坐标系可以利用的复杂形状，可以利用这套方法来建立一个与该物体相适应的坐标系，这就是适体坐标系的思想。5-5 适体

25、坐标系(Body-Fitted coordinates)一、适体坐标系的基本概念 如上所述，最理想的坐标系是其各坐标轴与所计算物体的边界一一相符合的坐标系，称为适 体坐标系，又称贴体坐标系，附体坐标系，无现成的坐标系可资利用时，希望通过计算来构造这 样的坐标系，如图所示。x-y坐标系中的不规则区域ABCD,把其相交的两个边界作为曲线坐标系的两个轴，记为g ,耳, 且：（1）一条边界上，只能一个坐标单值地发生为化，另一个坐标保持为常数；（2）在两条对应边界上，另一组曲线坐标的最大值与最小值对应相等，以便在计算平面上 得出矩形区域。问题：物理平面上，点（x, y）与计算平面上（g，耳）的对应关系？

26、求解结果的传送?如果把g，耳视为物理平面上的两个未知函数，那么上述确定g，耳的问题就是物理平面上的一 个边值问题，因此，从物理平面上来说，所谓要生成一个适体坐标系，实际上相当于要求解物理 平面上的一个边值问题。反之，首先把物理平面上的区域ABCD按已规定的边界上的g，耳值，画成计算平面上的g，耳坐标系中对应的矩形，然后以均匀网格离散化，于是问题变为：已知计算区域边界上各节点（g，耳） 相应的（g，耳），问在计算区域内部任意一点（g，耳）对应的（g，耳）值是多少？这样，如把X、 y视为计算平面的未知函数，则生成适体坐标系的问题即为计算平面中的一个边值问题。从数值计算的观点，对生成的适体坐标系有以

27、下几点要求：（1）物理平面上的点与计算平面上的点一一对应，同一簇中的曲线不能相变，不同簇中的 两曲线只能相交一次；（2）在适体坐标系中，每一个节点应当是一系列曲线坐标轴的交点，而不是一群三角形元 素的顶点或无序的点群，以便设计有效，经济的算法及程序，采用矩形网格即可；（3）物理平面求解区域内网格疏密程度易于控制；（4）在物理平面的适体坐标的边界上，网格线最好与边界线正交或接近正交，以便于边界 条件的离散化。二、生成适体坐标系的方法分类大致可分为三类1、复变函数法 利用复变函数的映射，可以把相当一批二维不规则区域变换成矩形区域，而且可以得出解极 或部分解析的变换关系式。g = (r - r )/

28、(r -1 )1 2 1耳=0/2兀相应地：x（gn）=r + (r - r )gcos(2兀耳)1 2 1y (g，n)=r + (r - r )gxin(2兀耳)1 2 1复变函数映射W 二 l z 二 l (reiQ)二 l r + i0n nng + in所以g = l rnn =e控制方程的变化： 字+专+丄学 =0or 2 or r 2 d0 2因为n = eo21/00 2 = o21/Oy2ot og o 2 tg = l ,Ot/Or = -+= 0仍然是各向同性！nrog or on 2代入原控制方程，得o 2t o 2t（条件心0）+= 0仍然是各向同性。og 2 oy

29、2可见,正交坐标系法与复变函数法并不相同！2、代数变换法利有一些代数关系式来进行区域变换,具体实施方法颇多,最简单的是边 界规范化的方法。3、解微分方程的方法通过求解边值问题的微分方程建立物理平面与计算平面上各点间的 对应关系。对该边值问题的控制方程的类型,物理问题本身无任何限定,可以比较自由地按照对所生成 网格的要求来选择控制方程。三、计算平面上求解区域形状的选择1、物理平面上的区域由四条两两相交曲线构成的单连通域一计算平面上取为正方形或矩形;2、物理平面的 L 型区域：有两种选择3、物理平面上的型区域：两种选择 特别注意：节点间一一对应关系不要受破坏,计算平面中求解节点方程时,对重迭线上的

30、节 点,每一个变量应有两套数组,以分别存放从左边和右边计算而得之值。注意：计算平面上的尖角点未必对应于 物理平面上的尖角点，反之亦然，例如物理 平面上一条连续封闭的曲线所围成的区域可 以转换成计算平面上的一个正方形。选择计算平面上求解区域形状的一条 件：生成的网格要适用于计算的问题，通常， 在计算平面上都采用正方形均匀网格。四、用适体坐标系求解流动与换热问题 的总体步骤1、网格生成：在计算平面上选择与物xx394理平面上复杂区域相应的求解区域，并找出两个区域内部节点之间的对应关系,X二x（g，耳）,y二y（g，耳）或其反函数。这些关系可以是解析的，也可以是数值的。2、控制方程的生成与离散，把物

31、理平面上的控制方程及边界条件转换成计算平面上的相应 形式，并利用控制容积积人法律立离散方程。3、离散方程的求解及解的传递。在计算平面上获得收敛的解后，根据节点间的对应关系， 传送到物理平面上去。5-6 控制方程的转换及离散化一、方程的转换讨论如何把物理平面上的稳态对流一扩散方程的通用形式转换成一般曲线坐标系g -n中的形式，物理平面上的通用形式为字*字等（碍）+ ?（碍）+ S （x, y）oxoyox ox oyoy帥o ogo on=y+oxog oxon oxoo ogo onM y+oyog oyon oy设已有函 数 x M x（g，n） y M y（g，n） 要求找出 反函 数 g

32、，n 的上述导数（用 已知函 数x M x（g，n）， y M y （ g ， n ）的导数来表示），这涉及到多元函数的微分关系：og1 oy1=yoxJ onJ non _ -1 oy ox J ogygon _ i ox oy J og1_ xJgog1 ox1oy _-J on _- Jxn，一， ” Ix Iy Ix ly式中，j代有Jacobi行列式：j=igm-iyitxg冋-xnyg.=丄住空处空=1( y y ) = 1 沁址 lxjog ininigj g n Jg j ogin1=7(y )g -(yg)J n g g =(-g x + xg) = (-X )g + (xg

33、)IyJ gg Jg g同理，对函数pu, pu可以写出:迪=打衍叫）-1（puyg） oxJogoI(pu) 1 I(pux ) I(pux )=一 g ly J扩散项的转换：，-等 yg)o gI1 I由：=r （丽yIx J Ig有：IX哙)唱r J(gyn-nyg)二 J忌J(gyn-nJ叫珥-叮叩J怙J (g V -nygyn )-InJ(n弋-gX叨同理有：(哼)二-J 捋J ( Xg X -g X2)-In J ( xg - X X )ly oyJ log J ngn g n In J n g gng源项：S (x, y) Sx(g，耳),y (g，耳)S (g，耳)两对流项相加

34、：得:-冲+二冲-2竺 J ogoJ o11P叫-匕)+ lP(匕-塔) = !沁 + 迪 jiginjigi式中，U二uy - ux , V二ux -uy，可分别视为计算平面上g及耳方向上的速度分量。两扩散项相加：得1 a ra r1 a ra ryU(x 2 + y2)-丽(y+ xg x) + +yf)0 (yg y+ X x)eJ ag Jgag J g gJ a J g g a J g g g=I 岂J 他-叫)+J 备J (-fV 叫)式中 a = x2 + y 2,r|r|r = x 2 + y 2ggax ax ay ayag ari+ ag ari计算平面上的守恒方程为：1 aJ ag(pu)+1 aj ar(pv) =1J agL J-J (agg1r一叫)+J J(-(V 叫)+S (g，r)可见，式中各项仍然保持了直角坐标中相应各项的意义，其中源项S完全由直角坐标中的源 项转换而来，其它各项在转换过程中并不给源项增添新的成分。生成的曲线坐标在物理平面上是否正交？0 = 0的点处，曲线坐标正交，且a和r即为相应度规系数的平方；0的点处，曲线坐标斜交，0偏离零的多少反应了局部曲线坐标偏离正交的程度。从物理平面到计算平面区域的简化（规则化）是以控制方程的复杂化为代价的。二、方程的离散化