Главная Журналы ЗНАЧЕНИЯ ПЕРЕМЕННЫХ, ЗАДАВАЕМЫХ ПО УМОЛЧАНИЮ Имя переменной 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 Распечатываются координаты расчетных точек Данные в файл для графики не записываются Поле переменной не выводится на печать
ПАМЯТКА ПО ПОСТРОЕНИЮ ПОДПРОГРАММЫ 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 |