[i, j] + apt * f [i, j, nf]) * vol; //коеф-т b [i, j]:=(apt-sp [i, j]) * vol;// коеф-т ар
end;
//коефіцієнти для поширення по Х
for j:=2 to M2 doi:=2 to L3 do:=arx[j]*2.*gam[i,j]*gam[i+1,j]/((Xcv[i]*gam[i+1,j]+Xcv[i+1]*gam[i,j]+small)*sx[j]);[i,j]:=diff+small;// коеф-т aE [i + 1, j]:=aip [i, j];// коеф-т aW; j:=2 to M2 do
//Розглядається межа i=1:=gam [2, j]/(0.5 * xcv [2] * sx [j]) + small;// дифузійна [2, j]:=beta * diff;// коеф-т aW [1, j]:=aim [2, j];// коеф-т aE [2, j]:=aim [2, j] * arx [j]; [1, j]:=(beta - 1) * aip [2, j]/arx [j]; [2, j]:=aip [2, j] + aim [1, j] * arx [j]; (kbci1 [j]=1) then [2, j]:=sc [2, j] + aim [2, j] * f [1, j, nf] [1, j]:=aip [1, j] -flxpi1 [j];// коеф-т AР [1, j]:=flxci1 [j];// коеф-т b:=aim [2, j]/sp [1, j]; [2, j]:=sp [2, j] -aip [1, j] * temp; [2, j] :=aip [2, j] -aim [1, j] * temp;// коеф-т aE [2, j]:=sc [2, j] + sc [1, j] * temp ;; [2, j]:=sp [2, j] + aim [2, j] ;
aim [2, j]:=0.;
//Розглядається межа i=L1
diff:=gam[L2,j]/(0.5*xcv[L2]*sx[j])+small;[L2,j]:=beta*diff;[L1,j]:=aip[L2,j];[L2,j]:=aip[L2,j]*arx[j];[L1,j]:=(beta- 1) * aim [L2, j]/arx [j]; [L2, j]:=aim [L2, j] + aip [L1, j] * arx [j]; (kbcL1 [j]=1) then[L2,j]:=sc[L2,j]+aip[L2,j]*f[L1,j,1][L1,j]:=aim[L1,j]-flxpL1[j];[L1,j]:=flxcL1[j];:=aip[L2,j]/sp[L1,j];[L2,j]:=sp[L2,j]-aim[L1,j]*temp;[L2,j]:=aim[L2,j]-aip[L1,j]*temp;[L2,j]:=sc[L2,j]+sc[L1,j]*temp;;[L2,j]:=sp[L2,j]+aip[L2,j];
aip [L2, j]:=0. ;;
//Коефіцієнти для поширення Y
for j:=2 to M3 doi:=2 to L2 do
begin:=rv [j + 1] * xcv [i];// площа грані КО
diff:=area*2.*gam[i,j]*gam[i,j+1]/(ycv[j]*gam[i,j+1]+ycv[j+1]*gam[i,j]+small);[i,j]:=diff+small;// коеф-т aN [i, j + 1]:=ajp [i, j];// коеф-т aS; i:=2 to L2 do
begin
//Розглядається межа j=1:=rv [2] * xcv [i];
diff:=gam[i,2]/(0.5*ycv[2])+small;[i,2]:=beta*diff;[i,1]:=ajm[i,2];[i,2]:=ajm[i,2]*area;[i,1]:=(beta- 1) * ajp [i, 2]/(rv [3] * xcv [i]); [i, 2]:=ajp [i, 2] + ajm [i, 1] * area; (kbcj1 [i ]=1) then[i,2]:=sc[i,2]+ajm[i,2]*f[i,1,nf][i,1]:=ajp[i,1]-flxpj1[i];[i,1]:=flxcj1[i];:=ajm[i,2]/sp[i,1];[i,2]:=sp[i,2]-ajp[i,1]*temp;[i,2]:=ajp[i,2]-ajm[i,1]*temp;[i,2]:=sc[i,2]+sc[i,1]*temp;;[i,2]:=sp[i,2]+ajm[i,2];
ajm [i, 2]:=0.;
//Розглядається межа j=M1
area:=rv[M1]*xcv[i];:=gam[i,M2]/(0.5*ycv[M2])+small;[i,M2]:=beta*diff;[i,M1]:=ajp[i,M2];[i,M2]:=ajp[i,M2]*area;[i,M1]:=(beta- 1) * ajm [i, M2]/(rv [M2] * xcv [i]); [i, M2]:=ajm [i, M2] + ajp [i, M1] * area; (kbcM1 [i ]=1) then[i,M2]:=sc[i,M2]+ajp[i,M2]*f[i,M1,2][i,M1]:=ajm[i,M1]-flxpM1[i];[i,M1]:=flxcM1[i];:=ajp[i,M2]/sp[i,M1];[i,M2]:=sp[i,M2]-ajm[i,M1]*temp;[i,M2]:=ajm[i,M2]-ajp[i,M1]*temp;[i,M2]:=sc[i,M2]+sc[i,M1]*temp;;[i,M2]:=sp[i,M2]+ajp[i,M2];[i,M2]:=0.;;
//Побудова масивів sp [i, j] і sc [i, j] j:=2 to M2 doi:=2 to L2 do:=aip [i, j] + aim [i , j] + ajp [i, j] + ajm [i, j];// сума сусідніх коеф-тів:=anb * rlx;// інерція I [i, j]:=sp [i, j] + anb + ainr;// коеф-т aP [i, j]:=sc [i, j] + ainr * f [i, j, 1];// коеф-т b
end;
//Виклик процедури Solve для отримання рішення дискретних аналогів
end ;;; ::=time + dt;:=iter + 1; (iter gt;=last) then kstop:=1 ;;
TForm1.Button1Click (Sender: TObject) ;; not TryStrToFloat (Form1.Edit1.Text, t1) or not TryStrT (Form1.Edit2.Text, t2) (StrToFloat (Form1.Edit2.Text) lt; 0) or (StrTFloat (Form1.Edit1.Text) lt; 0) then (handle, PChar ( Введіть коректні температури! ), PChar ( Помилка ), MB_OK + MB_ICONWARNING) t1 gt;= t2 then
begin (handle, PChar ( Температура на кордонах внутрішньої труби повинна + # 10 + бути не рівний і менше температури дільниці нагріву! ), PChar ( Помилка ), MB_OK + MB_ICONSTOP)
endStrToFloat (Form1.Edit3.Text) gt;=StrToFloat (Form1.Edit1.Text) then (handle, PChar ( Температура рідини повинна бути менше температури стінки каналу! ), PChar ( Помилка ), MB_OK + MB_ICONSTOP) (StrToFloat (Form1.Edit2.Text) - StrToFloat (Form1.Edit1.Text) lt; 10) then.Label4.Caption:= Обрано матеріал кожуха: + Material + # 10 + На внутрішньому кордоні + Form1.Edit1.Text + C + # 176 + # 10 + По каналу тече рідина + Form1.Edit3.Text + C + # 176 + # 10 + Нагрівання йде з нижньої частини труби + Form1.Edit2.Text + C + # 176; (Form1.RadioGroup1.ItemIndex=2) then.Label1.Caption:= Матеріал кожуха: + Material + . Raquo ;;. Caption:= Нагрівання кожуха raquo ;;; (Form1.RadioGroup1.ItemIndex=1) then.Label1.Caption:= Обрана рідина: + Gidkost + . Raquo ;;. Caption:= Нагрівання рід...
