Введение в технологию программирования. Часть 3. Проект OLYMPUS-D. Коржов Е.Н - 66 стр.

UptoLike

ВУЗ: 

ВГУ | Воронеж

Составители: 

Рубрика: 

66
W[i,Nx] := 3/sqr(Hx)*(Psi[i,Nx]-Psi[i,Nx-1])-0.5*W[i,Nx-1];
end;
for
j := 0 to Nx do
begin
W[0,j] := 2/sqr(Hy)*(Psi[0,j]-Psi[1,j])-0.5*W[1,j];
W[Ny,j] := 2/sqr(Hy)*(Psi[Ny,j]-Psi[Ny-1,j])-0.5*W[Ny-1,j];
end;
end
; { WoodFormulas }
procedure JensenFormulas;
{
Назначение: Вычисление функции вихря на границах области реше-
ния задачи по формулам Йенсена
}
var
i, j : byte;
begin
for
i := 0 to Ny do
begin
W[i,0] := 0.5/sqr(Hx)*(7*Psi[i,0]-8*Psi[i,1]+Psi[i,2]);
W[i,Nx] := 0.5/sqr(Hx)*(7*Psi[i,Nx]-8*Psi[i,Nx-1]+Psi[i,Nx-2]);
end;
for
j := 0 to Nx do
begin
W[0,j] := 0.5/sqr(Hy)*(7*Psi[0,j]-8*Psi[1,j]+Psi[2,j]);
W[Ny,j] := 0.5/sqr(Hy)*(7*Psi[Ny,j]-8*Psi[Ny-1,j]+Psi[Ny-2,j]);
end;
end
; { JensenFormulas }
Аналогично создаются процедуры для других способов вычисления
функции вихря на границах области. Отметим, что можно без особого тру-
да наращивать функциональные возможности комплекса без перепрограм-
мирования существующего кода, расширяя лишь список возможных аль-
тернативных ветвей оператора case в процедуре вычисления вихря в гра-
ничных точках области решения
Calculate_Vorticity_in_Boundary_Points. 
          W[i,Nx] := 3/sqr(Hx)*(Psi[i,Nx]-Psi[i,Nx-1])-0.5*W[i,Nx-1];
      end;
   for j := 0 to Nx do
      begin
          W[0,j] := 2/sqr(Hy)*(Psi[0,j]-Psi[1,j])-0.5*W[1,j];
          W[Ny,j] := 2/sqr(Hy)*(Psi[Ny,j]-Psi[Ny-1,j])-0.5*W[Ny-1,j];
      end;

end; { WoodFormulas }


procedure JensenFormulas;
{
   Назначение: Вычисление функции вихря на границах области реше-
                ния задачи по формулам Йенсена

}
var
   i, j : byte;
begin
   for i := 0 to Ny do
        begin
            W[i,0] := 0.5/sqr(Hx)*(7*Psi[i,0]-8*Psi[i,1]+Psi[i,2]);
            W[i,Nx] := 0.5/sqr(Hx)*(7*Psi[i,Nx]-8*Psi[i,Nx-1]+Psi[i,Nx-2]);
        end;
   for j := 0 to Nx do
        begin
            W[0,j] := 0.5/sqr(Hy)*(7*Psi[0,j]-8*Psi[1,j]+Psi[2,j]);
            W[Ny,j] := 0.5/sqr(Hy)*(7*Psi[Ny,j]-8*Psi[Ny-1,j]+Psi[Ny-2,j]);
        end;

end; { JensenFormulas }

    Аналогично создаются процедуры для других способов вычисления
функции вихря на границах области. Отметим, что можно без особого тру-
да наращивать функциональные возможности комплекса без перепрограм-
мирования существующего кода, расширяя лишь список возможных аль-
тернативных ветвей оператора case в процедуре вычисления вихря в гра-
ничных точках области решения Calculate_Vorticity_in_Boundary_Points.


                                       66

Страницы