【正文】
i j Fj i j j i j j i j j i j V i Fj i j j i j j i j j i j L i Fj i jL x V y V yL x L x V y V y j NL x L x V y V y f j NL x L x V y V y f j NLx? ? ? ?? ? ? ?? ? ? ???? ? ? ?? ? ? ? ? ? ?? ? ? ? ? ? ?? ? ? ? ? ??, , 1 , 11 , 1 , ,( ) 0 , . . . . . . . . . 1 1( ) 0 , . . . . . . . . . . . . . . . . . . . . . . . . . .j i j j i j j i j F Tj i j j i j j i j TL x V y V y N j NL x L x V y j N????? ? ? ? ? ? ?? ? ? ?, , ,i j i j i jy k x?變形后的 M方程組 上式方程組是一個有三對角矩陣的線性方程組 第一個方程僅有 1, 2下標(biāo)的變量,最后一個方程 僅有, NT- 1, NT下標(biāo)的變量,其他都是三個變量。 同時如把平衡方程 代入后消去 y后 可得下述方程組。 下面是計算框圖: 22( ) ( ( ) ( 1 ) ) 0 . 0 1( ) ( ( ( ) ( 1 ) ) / ( ) )j j Tj j jF T T k T k NF T T k T k T k e s p? ? ? ?? ? ? ???計算框圖 輸入已知條件 設(shè) Tj, Vj Ki, j= f( Tj, Pj) 物料衡算計算 Lj 三對角矩陣計算 xi, j Xi, j正規(guī)化 泡點方程修正 Tj, ki, j F(T)esp 計算個 Hj和 hj 二對角方程修正 Vj 是否收斂 輸出 計算 M次 第一次用理想體系 選擇迭代方法 選擇迭代方法