Главная  Журналы 

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 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

ЗНАЧЕНИЯ ПЕРЕМЕННЫХ, ЗАДАВАЕМЫХ ПО УМОЛЧАНИЮ

Имя переменной

ALAM(I, J) BIG

CRIT(NF) DT

F(I,J,N) FLXCIl(J) FLXCJl(I) FLXCLl(J) FLXCMl(I) FLXPll(J) FLXPJl(I) FLXPLl(J) FLXPMl(I) GAM(I, J)

IBLOCK(I, J)

ITER

KBCIl(J) KBCJl (I) KBCLl(J) KBCMl (I) KBLOC(NF) KORD

KOUT

KPGR

KPLOT(NF) KPRINT(NF)

Значение

LE20

l.E-5

LE20

Комментарий

Значение X = 1

Очень большое число (аналог бесконечности)

Критерий выхода в алгоритме прогонки Шаг по времени для стационарных задач Все Г1еремеш1ые изначально равны нулю /(. = О на границе с I = 1 Го же J = I То же I = L1 То же J = М1

= О па rpatmue с I = 1 То же J = 1 Го же I = LI То же J = М1

Везде одинаковое значение Г, равное единице

В расчетной области отсутствуют

заблокированные КО

Начальное значение счетчика итераций

Задано значение ф при 1 = 1

То же J = 1

То же I = L1

То же J = М1

Выполняется блочная коррекция На границах реализуется аппроксимация второго порядка Результаты выводятся на экран и в файл с именем PRINTF

Распечатываются координаты расчетных точек

Данные в файл для графики не записываются

Поле переменной не выводится на печать



KSOLVE(NF)

Уравнение для переменной не решается

KSTOP

Выполнение расчетов не прекращается

LAST

Выполняются пять итераций

MODE

Декартова система координат

NTIMES(NF)

В подпрограмме SOLVE выполняются 10 повторов алгоритма

PLOTF

PLOTl

Имя файла для графической обработки

POWERX

Равномерная сетка по оси х

POWERY

То же по оси у

POWRX(NZ)

Во всех зонах равномерная сетка по оси х

POWRY(NZ)

То же по оси у

PRINTF

PRINTl

Имя файла для вывода результатов

R(l)

Радиус при J = 1

RELAX(NF)

Релаксации не используются

SC(I,J)

Коэффициент S, = 0

SMALL

l.E-20

Малое (практически нулевое) число

SP(I,J)

Коэффициент Sp = 0

TIME

Начальный момент времени



ПАМЯТКА ПО ПОСТРОЕНИЮ ПОДПРОГРАММЫ ADAPT

Общая часть. Используйте выражения с операторами DIMENSION и EQUIVALENCE с целью введения удобных имен для частей массива F.

GRID. Определите переменную HEADER и имена файлов для вывода результатов и графической обработки. Затем задайте MODE, Ll, Ml, XU(I) (I = 2,L1)hYV{J) (J = 2,M1). Если MODE ?4 1, задайте также R (1). Помните, что 4 < Ll < NI и 4 < Ml < NJ. Для построения равномерной сетки необходимо задать MODE, NCVLX, NCVLY, XL и YL, a также R{1) для MODE Ф 1, затем вызвать подпрограмму EZGRID. С использованием переменных POWERX и POWERY можно ввести простую неравномерность сетки. Для более сложных сеток следует использовать подпрограмму ZGRID с соответствующими входными данными.

BEGIN. Задайте значения для массивов переменных KSOLVE (NF), KPRINT (NF) , TITLE (NF), RELAX (NF), NTIMES (NF), CRIT (NF) если для задачи неприемлемы их значения по умолчанию. Заполните массивы F (1, J,NF) их начальными значениями, используя корректные значения на границах, если они известны.

OUTPUT. Реализуйте вывод любых желаемых результатов. Пересчитайте величины, которые изменяются после каждой итерации или щага по времени. Для распечатки полей переменных F (I, J, NF) можете вызвать подпрограмму PRINT.

PHI. Для каждого используемого NF задайте массив ALAM (I, J), GAM(I, J) , SC (I, J) hSP{I,J) при I = 2, L2 и J = 2,M2. Определите подобающие значения КВС и соответствующие FLXC и FLXP для граничных точек, в которых не заданы значения зависимой переменной.





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 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