www.chms.ru - вывоз мусора в Жуковском
Читаемые статьи

Читаемые книги

Ссылки


Главная >  Вычислительная гидроаэродинамика 

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

5.2.3, FIVOL: метод конечных объемов в применении к реилению уравнения Лапласа

Описанный выше метод реализован в компьютерной программе FIVOL, распечатка которой приводится на рис. 5.6. Основные параметры, используемые в программе FIVOL, опи-<:ываются в табл. 5.4.

Таблица 5.4. Параметры, используемые в программе FIVOL

Параметр

Описание

JMAX

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

КМАХ

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

NMAX

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

Tw (рис. 5.4)

гх (рис. 5.4)

гу (рис. 5.4)

гг (рис. 5.4)

THEB

Qwx (рис. 5.4)

THEN

OzK (рис. 5.4)

Критерий сходимости ПНР

Показатель релаксации h в формуле (5.42)

PHIX

XA, YA

Ха, у а

Ах а в

Axk-i,k

110*-011гтз, 110+*-01гтз

Программа FIVOL считывает и выдает исходные управ-.ляющие параметры (строки 1-25). Задается сетка, рассчитывается начальное решение (для проведения ПВР), а также точное решение (строки 38-53). Поскольку интерес представляет только сходящийся процесс решения, то точное решение также использовалось в качестве одного из начальных решений для ПВР. Этот вариант приводит к быстрой сходимости процесса по схеме ПВР. Задаются граничные условия (строки 57-64); строго говоря, эта часть программы не является необходимой, так как граничные условия уже были в неявной форме учтены в строке 5L

На каждой стороне области задаются параметры, связан-яые с сеткой {Qab, Раз и т. д.); это отражено в строках 68-123.



3 С FIYOb APPLIES THE FINITE VOLUHE METHOD TO THE SOLUTIOK OF

4 С LAPUCES EQUATION IN CARTESIAN COORDINATES ON A POLAR GRID*

5 С THE DISCRETISED EQUATION IS SOLVED BY SOR

7 DIMENSION X(21,21) Д(21,21) ,QAB(21,21) ,PAB(21,2.1) QBC(21,21)

8 1PBC(21,21) ,QCD(21,21) ,PCD(21,21) ,QDA(21,21.) ,PDA(21r2X),

9 2PHI(21,21),PHIX(21,21)

10 С

11 OPEN(1,FILE* FIVOL.DAT)

12 OPEN(б,FILE- FIVOL.OUT)

13 READ(1,1)JMAX,KMAX,NMAX

14 READ(1,2)RV,RX,RY,RZ,THEB,THEN,EPS,0H

15 1 FORMAT(815)

16 2 FORMAT(8E10.3) IT С

18 VRITE(6,3)

19 VRITE(6,4)JMAX,KMAX,NMAX,EPS,0M JtQ VRITE (6,5) RV, RX, RY,1tZ, THEB, THEN

21 3 FORMAT(! LAPLACE EQUATION BY FINITE VOLUME METHOD, )

22 4 FORMATC JMAX-M2; mX M2, NMAX *;I5,

23 15Х/ EPS-4E10.3, OM \F5.3)

24 5 FORMATC RW ,F5.3/ RX-\F5.3, RY ,F5.3/ RZ *,F5.3r

25 15Х/ THEB *,F5.1, THEN ,F5a, )

26 С

27 JMAP JMAX 1

28 KMAP Ш KMAX - 1

29 AJM JMAP

30 ARM Ш KMAP

31 DRVX Ш (RX - RV)/AJM

32 DRZY (RY - RZ)/AKM

33 Ш (THEN-THEB)/AKM

34 PI m 3.141592T

35 С

36 С SET X, Y, EXACT AND INITIAL PHI

37 С

38 DO 7 К .1,KMAX

39 AK К - 1

40 ТНК (THEB + AK*DTH) PI/180.

41 CK COS(THK)

42 SK SIN(THK)

43 DR * DRVX -I- (DRZY - DRVX)*AK/AKM

44 RWZ RW + (RZ - RW)*AK/AKM

45 DO б J W l.JMAX

46 AJ J 1

47 R RVZ AJ*DR

48 X(J,K) R*CK

49 Y(J,K) R*SK

50 PHIX(J,K) SK/R

51 PHI(J,K) PHIX(J,K)

52 6 CONTINUE

53 7 CONTINUE

54 С

55 С SET BOUNDARY VALUES OF PHI

56 С

57 DO 8 J 1,JMAX

58 PHI(J,i) 0.

59 PHKJ.KMAX) = PHIX(J,KMAX)

60 8 CONTINUE

Рис. 5.6. Распечатка программы FIVOL (начало).



DO 9 К 1ДМАХ

РН1(1Д) РН1Х(1Д)

PHI (JMAX Д) РН1Х(аМАХД)

9 CONTINUE

65 С

66 С

SET GRID REbATED PARAMETERS

67 С

DO 11 К 2,KMAP

KM К - 1

KP К + 1

DO 10 J 2,JMAP

JM J - 1

JP J + 1

XA 0.25*(X(J,K)

+ X(JM,K) +

X(JM,KM)

+ X(J,KM))

YA 0.25*(Y(J,K)

+ YCJM,K) ♦

Y(JM,KM)

+ Y(J,KM))

XB Ш о.25*(х(ад)

+ X(J,KM) ♦

ХЫРДМ)

+ X(JP,K))

YB 0.25*(Y(J,K)

+ Y(J,KM) +

Y(JPM)

+ Y(JP,K))

XC 0.25*(X(J,K)

+ X(JP,K) +

X(JP,KP)

+ X(J,KP))

YC 0.25*(Y(J,K)

+ Y(JP,K) +

Y(JP,KP)

+ Y(J,KP))

XD 0.25*(X(J,K)

+ X(J,KP) +

X(JM,KP)

+ X(JM,K)) + Y(JM,K))

YD . 0.25*(Y(J,K)

+ Y(J,KP) +

Y(JMP)

Z2 С

83 С

SIDE AB

S4 С

DXA XB - XA

DYA YB - YA

DXK X(J,K) - X(J,Krl)

DYK Y(J,K) - Y(J,K-1)

SAB * ABS(DXA*DYK

- DXK*DYA)

QAB(J,K) (DXA*DXA -l- DYA*DYA)/SAB

PAB(J,K) (DXA*DXK DYA*DYK)/SAB

92 С

93 С

SIDE ВС

94 С

DXB XC - XB

DYB YC - YB

DXJ X(J,K) - X(J+1,K)

DYJ Y(J,K) - Y(J+1,K)

SBC Ш ABS(DYJ*DXB

- DXJ*DYB)

QBC(J,K) (DXB*DXB -1- DYB*DYB)/SBC

PBC(J,K) s (DXB*DXJ -1- DYB*DYJ)/SBC

102 С

103 С

SIDE CD

104 С

DXC XD - XC

DYC YD - YC

DXK X(J,K) - X(J,K+1)

DYK Y(J,K) - Y(J,K-H)

SCD ABS(DXC*DYK

- DYC*DXK)

QCD(J,K) (DXC*DXC ♦ DYC*DYC)/SCD

PCD(J;k) (DXC*DXK DYC*DYK)/SCD

112 С

113 С

SIDE DA

114 С

DXD = XA - XD

DYD YA - YD

DXJ X(J,K) - X(J-1,K)

DYJ Y(J,K) - Y{J-1,K)

Рис. 5.6 (продолжение).



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


Чем хороши многотопливные котлы?



Нетрадиционное отопление



Детище отечественной Оборонки



Что такое автономное индивидуальное отопление?



Использование тепловых насосов



Эффективное теплоснабжение для больших помещений



Когда удобно применять теплые полы
© 1998 - 2024 www.300mm.ru.
При копировании материала обязательно наличие обратных ссылок.
Яндекс.Метрика