Разделы сайта
Читаемое
Обновления Mar-2024
|
Промышленность Ижоры --> Динамика жидкости: уравнения 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
лежащие внутри расчетной области (рис. 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
|
© 2003 - 2024 Prom Izhora
При копировании текстов приветствуется обратная ссылка |