Android-приложение для поиска дешевых авиабилетов: play.google.com
Главная -> Применение эвм

0 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


SUBROUTINE MKE(X,Y,IHD,M,N,AAA,BBB.li!S,G.F)

ПРИВЕПРШЕ КООРДИНАТ УЗЛОВ К ЦЕИТРУ

ПОЦПРОГРАМЫА ФОПЙ1Р0ВАНИЯ МАТРИЦЫ МЕТОДА КОНЕЧНЫХ

Х1=Х(1)-ХС

3JEHEHT05 ДЛЯ ДВУМЕРНОЙ СТАЦИОНАРНОЙ ЗАДАЧИ ТЕПЛОПРОВоДносп,

5 56

X2=X(J)-XC ХЗ=Х(К)-ХС

ОПИСАНИЕ ПАРАМЕТРСВ:

fi=ftI)-YC

М - ЧИСЛО УЗЛОВ

.Y2=Y(J)-YC

а - число КОНЕЧНЫХ элйентов

ВХОДНЫЕ

Y3=WK)~YC

)((M),Y(M) - HACCHHJ КООРДИНАТ УЗЛОВ

ПАРАМЕГГРЫ

IND{4«N) - ИНДЕКСНАЯ МАТРИЦА

ВЫЧИСЛЕНИЕ КОЭСФИЦИНТГОВ ФУНКЦИИ ФОРШ

ААА.ВВБ - ИМЕНА ВНЕШНИХ ПОДПРОГРАММ

&4

MS - ШРИНА ЛЕНТЫ ЛЕНТОЧНОЙ МАТРИЦЫ

ВЫХОШЬЕ

B1=Y2-Y3

G - ВЫХОДНАЯ ЛЕНТОЧНАЯ МАТРИЦА

С1=ХЗ-Х2

F - ВЫХОДНСЙ ВЕНДОР-СТОЛБЕЦ ПРАВЫХ ЧАСТЕЙ \ "AFMEWbl

B2=Y3-Y1

С2=Х1-ХЗ

DIMENSION X(l),Y{l),IND(l),G(l),F(i)

B3=Y1-Y2

СЗ=Х2-Х1

ВЫЧИСЛЕНИЕ ШРИНЫ ЛЕНТЫ MS

ВЫЧИСЛЕНИЕ ПЛОЩАДИ ТРЕУГОЛЬНИКА

MS=e

LL=I

А=(Х2*УЗ-ХЗ«У2+Х1*У2-

№ ! L=1,N

«XbY3+X3«Yl-X2.Yl )/2

J-lNt-(LL+l)

ШЧИСЛЕНИЕ ЗЛЕМЕтоВ ЛСЖАЛЬНОЙ МАТРИЦЫ

K=IKD(a+2)

И ВЕКТОРА ПРАВЫХ ЧАСТЕЙ

MS=MAXe(MS,IABS(I-J),

»IABS(I-K),IABS(J-K))

CALL AAA(ALX.ALy,ALFV,QV,XC,YC)

1 LULL+4

Н=4«А

D=ALFV*A/6

ОЧИС™ МАССИВОВ G,F

Gll=(ALX*BI*BItALY«Cl*!CI)/H-*-])

G12=(ALX»Bl»B2+ALY«Cl»C2)/H+D/2

J=M+MS»(2*M-MS-l)/£!

G13-- {ALX»B I .B3t ALY »C 1 »C3) /H+Il/2

DO 2 bl.J

G2£=(ALX»B2»B2>ALY.C2 »C2)/Н+D

2GtO=0

G23=(ALX»B2«B3tALY»C2»C3)/H+D/2

DO 3 1=1,M

G3?=(AU«B3*B3+ALY«C3«C3)1+D

3 F(l)-e

f;-5v.a/3

F2=F!

с НАЧАЛО ЦИКЛА ПС ТРЕУГОЛЬНИКАМ

F3=F1

LL=1

ВЫЧИСЛЕНИЕ ДОПОЛНИТЕЛЬНЫХ СЛАГАЕМЫХ

DO 12 L=1,N

ДЛЯ ГРАНИЧНЫХ ЭЛЕМЕНТОВ

НОМЕРА УЗЛОВ

LL=LL+3

IF(IND(LL).EQ,»)GO ТО 5

I = IND(a)

IFCIND(LL).EQ.1)G0 ТО 4

J=IND{LL+1)

K=IND(LU2)

ВКЛАД СТОРОНЫ JK

ВЫЧИСЛЕНИЕ КООРДИНАТ ЦЕНТРА ТРЕУГОЛЬНИКА

CALL BBB(ALFL.QL.(X(J)+X(K))/2,(Y(J)+Y(K)

H=SQRT((X2-X3)»«2+(Y2-Y3)«»2)

XC=(X{l)+X(J)+X(K))/3

IfcALFL»H/3

ir-(y(I)ty(J)ty(K))/3.

G22=G224D

Рис. 4.15

G33=G33+D

Рнс. 4.15. Продолжение





G23-G23-t-D/2

КМЩ ЦИКЛА го ТРЕЭТОЛЫШАЫ"

D=QL*iH/2

Fl=Fl+I>

12 Ш-LL+l

F3=F3+D

ВКЛАД СТОРОНЫ и

ФОРМИРОВАНИЕ МАТРИЦЫ В ВИДЕ, ПРИГОДКШ

ДЛЯ 0БРАЩЕМ1Я К ПОДПРОГРАММЕ ИОВ

4 CALL BBB(ALFL.3L,CX(I)+X(J))/2.(y(I)+y(J))/2)

H=SQRT((Xi-X2)**i2+(YI-Y2)«2)

13 M1=M1N0(N1.H2)

D-ALR.H/a

M2=UAX0(N1,N2)

GJJ=GJ]+D

Ш=]4М2-М1

G22=G22+D

HS1=MS+1

G1212+D/2

MMSKS

12»

D=QL»H/2

IF(M1.GE.HMS)G0 TO 14

Fl-FUD

NK=MS1<(M1-1 )+Ив1

F2=F2+D

GO TO 15

14 HM№=M-MU1

Л0БАВЛ£К;1Е ЭШЕШЪ ЛОКАЛШИ

№(=MSl#)lftlS+ie»KSi/2-««*(MMM+l)/2t«M

МАТРИЦЫ К ЭЛЕМЕЖАН ГЛОБАЛЫК)

15 G(NN)=G{KN)+D

GO TO LBL,(6,7,6,9,1®,11)

5 N1-1

123 130

D=G11

ASSIGN е ТО LBL

Рис. 4.15. Продолжение

GO TO 12

6 N2=J

133 134 135 136 137 138 139 14в 141 142 143 144 145 146 147 146 149 150 151 152 153 154 155 356 157 156 159

I>=Gi2

ASSIGH 7 TO LBL GO TO 13 7 N2=K IM]13

ASSIGK 8 TO LBL GO TO 13 6 N1=J K2=J IbG22

ASSIGN 9 TO LM-GO TO 13 9 N2=K D=G23

ASSIGK 10 TO LBL GO TO 13 le N1=H D=G33

ASSIGN U TO LBL GO TO J3

TO IE ДЛЯ ВЕКТОРА ПРАВЫХ ЧАШЯ

и F(l)=F(I)tFI F(J)=F(J)+F2 F(K)=F(K)+F3

Рис. 4.15. Продолжение

no индексной матрице определяются глобальные номера узлов i, /, k в данном элементе и принадлежность сторон границе;

вычисляются коэффициенты функций формы элемента по (4.14) и рассчитываются коэффициенты локальных матрицы н столбца по (4.29), (4.30),

значения всех коэффициентов локальных матрицы н вектор-столбца прибавляются к текущим значениям коэффициентов глобальных матрицы н вектор-столбца, индексы которых определяются глобальяыми номерами узлов i, /, k.

S результате обхода всех треугольников оказываются полостью сформированными глобальная матрица G н столбец правых астеи ф системы уравнений.

цднм, что прн формировании матрицы G необходимо учнты-*ан записи матрицы в машинной памяти для используемой В Дартной подпрограммы решения системы линейных уравнений. CHr случае предполагается использование подпрограммы щ математического обеспечения ЕС ЭВМ 1151, реализующей эт, "ДРйтного корня для симметричных ленточных матриц. Прн MaccHsP" матрицы должны быть записаны в одномерный г.11авцд<."Ум последовательного обхода верхней части ленты над Matpfjjj-"ональю по строкам. Такой пересчет индексов элемента S8-l индекс одномерного массива реализован операторами



ГЛАВА

КОНЕЧНО-РАЗНОСТНЫЕ МЕТОДЫ РЕШЕНИЯ ЗАДАЧ КОНВЕКТИВНОГО ТЕПЛООБМЕНА

§ 5.1. КОНЕЧНО-РАЗНОСТНЫЕ СХЕМЫ ДЛЯ УРАВНЕНИЯ ЭНЕРГИИ

Класс задач конвективного теплообмена весьма широк [19, 23] В связи с этим в данной книге нет возможности остановиться на сколько-нибудь подробном разборе особенностей даже основных постановок задач, н поэтому мы рассмотрим методы численного расчета только дтя некоторых наиболее простых, но довольно часто встречающихся иа практике задач решения уравнения энергии при заданном поле скоростей.

Для однофазного потока уравнение энергии можно записать в виде

(5.1)

где в левой части записана полная производная от энтальпии, включающая локальную н конвективную производные, а в правой части- члены, соответствующие тепловым потокам в единице объема, вызванным различными причинами. Первый член учитывает перенос теплоты за счет теплопроводности, второй - изменение энтальпии за счет работы сил давления, третий - выделение энергии за счет вязкого трения.

При моделировании процессов конвективного теплообмена уравнение энергии должно рассматриваться совместно с уравнениями неразрывности, движения и состояния. Прн анализе многих процессов, например в случае свободной конвекции или прн необходимости учета зависимости вязкости от температуры, необходимо все эти уравнения решать совместно. Численные схемы для уравнений гидродинамики гораздо сложнее, чем рассмотренные в главе 3 схемы для уравнения теплопроводности. С ними можно познакомиться по книгам [19-21, 23. Мы будем считать, что поле скоростей известно, т. е. либо была решена гидродинамическая задача при условии отсутствия влияния на нее температуры, либо было проведено соответствующее экспериментальное исследование.

Поясним основные особенности численных схем для уравнения энергии, введя еще ряд упрощений в уравнение (5.1). В случае малых скоростей течения (числа Маха М 4; 1) ш":ледние два члеи

1-56

Р правой части (5.1) можно не учитывать н рассматривать уравнение энергии вида

/ дТ

-vVT]V(kVT).

(5.2)

На первый взгляд схемы! рассмотренные в главе 3, легко пере-jjecTH на уравнение (5.2). Действительно, оно отличается от уравнения теплопроводности только членом IVT, содержащим первые производные от температуры по координатам, которые можно аппроксимировать конечными разностями. Однако некоторые варианты такого «естественного» подхода приводят к неудачным численным схемам. Поэтому новый конвективный член вносит ряд существенных особенностей в процедуру выбора вида разностной схемы. Рассмотрим их на примере простейшего одномерного стационарного уравнения энергии

(5.3)

Введем равномерную пространственную сетку Хп {п - 1) Я, ft = 1, Л, Конечно-разностное уравнение для внутренней точки будем строить методом баланса, выбрав элементарные объемы вида [д; Я/2, Хп + hl2\. Сеточную функцию численного решения обозначим, как обычно, через «„, п 1, Л. Уравнения баланса п-го элементарного объема (рис. 5.1) для единичного промежутка времени записывается так:

Энтальпия 0 {Xn + h/2) потока жидкости, выходящего из объема через сечение д:,1 + Я/2

/ h \

Энтальпия Ха-

потока

жидкости, входящего в объем через сечение х„, - h/2 Количество теплоты (Хп-гfil2), выносимое путем теплопроводности из элементарного объема через сечение Xn-j-h/2

Количество теплоты (х - h/2), вносимое путем теплопроводности в элементарный объем через сечение Хп - hf2

Q" (л-„ + л/2)-Q« ix„-h/2) = {Хп-Л/2)- (хп + (5.4)



0 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



0.0078
Яндекс.Метрика