Разделы сайта

Читаемое

Обновления Mar-2024

Промышленность Ижоры -->  Динамика жидкости: уравнения 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 [ 27 ] 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182

S6 Гл. 12. Обобщенные криволинейные координаты 1

2 SUBROUTINE TRAPA(JMAX,KMAX)

4 С FROM GRID COORDINATES CALCULATES TRANSFORM PARAMETERS

€ DIMENSION XG<23,23),YG<23,23),PHIX(21,21),PHI<21,21)

7 DIMENSION GW(21,21) ,GWT(21,21) ,GTT(21,21) ,DELZI(21,21),

8 1DELET(21,21)

9 COMMON /GRIDP/XG,YG,PHIX,PHI

10 COMMON /TRAPP/GWW,GWT,GTT,DELZI,DELET

11 DO 2 К > 1ДМАХ

12 KP = K+1

13 KPP К + 2

14 DO 1 J > 1,JMAX

15 JP J + 1

16 JPP J + 2

17 С

IS С BASIC TRANSFORM PARAMETERS

19 С

20 XZI 0.5*(XG(JPP,KP) - XG(J,KP))

21 YZI 0.5*(YG(JPP,KP) - YG(J,KP))

22 XET 0.5*(XG(JP,KPP) - XG(JP,K)) i3 YET 0.5*(YG(JP,KPP) - YG(JP,K))

24 XZZ XG(J,KP) > 2.*XG(JP,RP) * XG(JPP,KP)

25 YZZ - YG(J,KP) - 2.*YG(JP,KP) + YG(JPP,KP)

26 XEE XG(JP,K) - 2.*XG(JP,KP) + XG(JP,KPP)

27 YEE YG(JP,K) - 2.*YG(JP,KP) + YG(JP,KPP)

28 XZE 0.25*(XG(JPP,KPP)-XG(J,KPP)+XG(J,K)-XG(JPP,K))

29 YZE 0.25*(YG(JPP,KPP)-YG(J,KPP)+YG(J,K)-YG(JPP,K))

30 AJ XZI*YET - XET*YZI

31 С

32 С MODIFIED METRIC TENSOR COEFFICIENTS, G11,G12,G22

33 С

34 GWW(J,K) = (XZI*XZI + YZI*YZI)/AJ

35 GWT(J,K) -2.*(XZI*XET + YZI*YET)/AJ

36 GTT(J,K) (XET*XET + YET*YET)/AJ

37 С

38 С MODIFIED DEL**2ZI AND DEL**2ETA

39 С

40 DUM GTT(J,K)*(XET*YZZ-YET*XZZ)/AJ

41 DUM = DUM + GWT(J,K)*(XET*YZE-YET*XZE)/AJ

42 DELZI(JД) = DUM + GWW{J,K)*{XET*YEE-YET*XEE)/AJ

43 DUM = GTT(J,K)*(YZI*XZZ - XZI*YZZ)/AJ

44 DUM = DUM + GWT(J,K)*(YZI*XZE-XZI*YZE)/AJ

45 DELET(J,K) = DUM + GWW(J,K)*(YZI*XEE-XZI*YEE)/AJ

46 1 CONTINUE

47 2 CONTINUE

48 RETURN

49 END

Рис. 12.11. Распечатка подпрограммы TRAPA.

- (GWT.);+ + [(QWW..)/,- 2 (GWW ,. * +

+ (GWW.sb)/.ft+,] = 0. (12.84)

При записи (12.84) предполагалось, что в расчетной области Л = Ат1= 1. Уравнение (12.84) можно решить методом верхней



2 SUBROUTINE ITER(JMAX,КИАХ,NMAX.N,ОМ,EPS)

4 с ITERATE USING SOR APPLIED TO DISCRETISED EQUATIONS

6 DIMENSION GWW(21,21),GWT(21,21),CTT(21,21),DELZI(21,21).

7 1DELET(21,21),PHI(21,21),PHIX(21,21),XG(23,23),YG(23,23)

8 COMMON /GRIDP/XG,YG,PHIX,PHI

9 COMMON /TRAPP/GWW,GWT,GTT.DELZI,DELET

10 KMAP KMAX-1

11 JMAP = JMAX-1

12 АХИ = KMAP-1

13 AJM = JMAP-1

14 с

15 DO 3 N = 1,NMAX

16 SUM = 0.

17 DO 2 К = 2,KMAP

18 KM К - 1

19 KP * К + 1

20 DO 1 J = 2,JMAP

21 JM = J - 1

22 JP = J + 1

23 PHD = -0.5*(DELZI(JP,K)*PHI(JP,K)-DELZI(JM,K)*PHI(JM,K))

24 PHD - PHD-0.5*(DELET(J,KP)*PHI(J,KP)-DELET(J,KM)*PHI(J,KM)>

25 PHD - PHD + GTT(JM,K)*PHI(JM,K) + GTT(JP,K)*PHI(JP,K)

26 PHD = PHD + GWW(J,KM)*PHI(J,KM) + GWV(J,KP)*PHI{J,KP)

27 PHD = PHD + 0.25*(GWT(JP,KP)*PHI(JP,KP)-GWT(JM,KP)*PHI(JM,KPr

28 1 + GWT(JM,KM)*PHI(JM,KM) - GWT(JP,km)*PHI(JP,KM))

29 PHD = 0.5*PHD/(GTT(J,K)+GWW(J,K))

30 DIF = PHD - PHI(J,K)

31 SUM SUM + DIF*DIF

32 PHI(J,K) PHKJ.K) + OM*Dir

33 1 CONTINUE

34 2 CONTINUE

35 RMS SQRT(SUM/AJM/AKM)

36 IF(RMS .LT. EPS)RETURN

37 3 CONTINUE

38 VRITE(6,4)NMAX,RMS

39 4 FORMATC CONVERGENCE NOT ACHIEVED INM5, STEPS, RMS=,

40 1Б12.5)

41 RETURN

42 END

Рис. 12.12. Распечатка подпрограммы ITER.

релаксации (§ 6.3). Значение на следующей итерации получается из (12.84) в виде

k = {-0.5 [(DELZI./+, , - (DELZI.Л;-,..] -- 0.5 [(DELET.0 ),.- (DELET. )/, + + KGTT. );. , + (QTT..);+,., + (GWW.r) + + (GWW.n/.*+.] +

+ 0.25 [(GWT.- (GWT.r)/-!, + + (GWT. ); ,.- (GWT.n+i. ft-.]} X

X {2 (GTT/., + GWW/.,)}-. (12.85)> <i>tk =Ф1н + ifl k - ф . k). (I2.86>



где X - параметр релаксации. Данный алгоритм реализован в программе LAGEN. Распечатка программы LAGEN и подпрограмм GRID, TRAPA и ITER приведена на рис. 12.9-12.12. Основные параметры, используемые в программе LAGEN, описаны в табл. 12.1.

После считывания различных параметров программа LAGEN вызывает подпрограмму GRID для определения координат точек сетки, которые хранятся в массивах XG, YG. Точки сетки,

Таблица 12.1. Идентификаторы, используемые в программе LAGEN

Идентификатор

Описание

JMAX

Число точек в радиальном направлении

KMAX

Число точек в окружном направлении

NMAX

Максимальное число итераций

rw (рис. 12.8)

гх (рис. 12.8)

гу (рис. 12.8)

rz (рис. 12.8)

THEB

Qwx (рис. 12.8)

THEN

QzY (рис. 12.8)

Точность сходимости SOR

Параметр релаксации X в (12.86)

PHIX

Точное решение для ф (12.78)

0* в (12.85) (подпрограмма ITER)

XG, YG

Координаты расширенной сетки (подпрограмма GRID)

ф* 0 \\фП + 1 ф

лежащие внутри расчетной области (рис. 12.8), дополняются рядами точек (/ = 0, JMAX+1, й = 0, KMAX-fl), лежащих за границей расчетной области. Это позволяет использовать для определения параметров и т. д. на границе области те же формулы (12.30), (12.31), что и внутри.

Различные параметры сетки GTT и т. д., входящие в (12.85), определяются в подпрограмме TRAPA, а итерационный процесс решения (12.85), (12.86) реализуется подпрограммой ITER. Впоследствии программа LAGEN формирует выдачу решения.

Типичная выдача для сетки 6X6 приведена на рис. 12.13. В данном случае сетка ортогональная и члены р и V/J



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 [ 27 ] 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182

© 2003 - 2024 Prom Izhora
При копировании текстов приветствуется обратная ссылка