Тепло и массообмен в пограничных слоях - Патанкар С.
Скачать (прямая ссылка):
Подпрограмма COEFF. Является по существу сердцевиной программы. Она предназначена для введения коэффициентов Аи, Ви, Си уравнений (2.4-35), а также коэффициентов А, В, С уравнений (2.4-31) для каждого Ф-уравнения и всех промежуточных точек сетки. Она начинается с вычисления g-коэффициентов уравнения скорости с последующим определением С-величин, входящих в конвективный член. Затем различные члены группируются согласно равенствам (2.4-36) — (2.4-38) и (2.4-32) — (2.4-34). Для этой цели используются различные накопительные устройства машинной памяти. Полученные такие образом совокупности величин дополняются коэффициентами с подстрочными индексами 2 и JV+2, которые вырабатываются с помощью специальной подпрограммы SLIP. К ее рассмотрению мы сейчас и переходим.
Подпрограмма SLIP. Представленные в § 2.5 соотношения связывают величины скольжения с соседними истинными величинами и содержатся в виде соответствующих коэффициентов в подпрограмме SLIP. Этн коэффициенты вычисляются для обеих границ I и Е в зависимости от характера границы. В случаях, когда границей является линия симметрии, расчеты плоского и осесимметричных течений различаются между собой. Если границей служит стенка, то граничное условие вводится подпрограммой FBC. Тогда для величин скольжения используются различные соотношения в зависимости от того, задается ли на стенке поток либо Ф-велпчина.
Подпрограмма SOLVE. Содержит математические операции (описанные в §2.7) непосредственно решения системы уравнений типа (2.7-1). Подпрограмма вводится командой
CALL SOLVE (А, В, С, F, NP3),
где А. В, С — таблицы заданных коэффициентов, F — таблица Ф-величин, a NP3 определяет размер таблицы F, включающей известные значения на границах.
Подпрограмма READY. После каждого интегрирования мы получаем новые значения и и других Ф-величин для значений со. Подпрограмма READY производит вычисление соответствующих расстояний по нормали ординат у н радиусов поперечной-кривизны г для каждой точки сетки. Это позволяет осуществить переход к следующей операции интегрирования.
Подпрограмма READY требует введения величин плотности р во всех точках сетки и задания характеристик г, и а. Они заимствуются 'из подпрограмм DENSTY и RAD, Формулы, используемые в подпрограмме READY, приведены в разд. 2.8-2.
П1-3-2. Подпрограммы ввода физических гипотез
Подпрограмма VEFF. Предназначена для ввода величины [хЭф в соответствии, с используемой гипотезой.
Ввод осуществляется командой
CALL VEFF (I + I + l, EMU),
которая означает, что требуемое значение [х0ф, EMU должно быть вычислено посередине между линиями сетки I, 1 + 1.
В программе, приведенной в § Ш-5, используется гипотеза о пути смешения в форме уравнения (1.3-5). Характерная длина yi согласно определению разд. 2.8-3 вводится подпрограммой LENGTH.
( CALL LENGTH
, r~
( callentrnУ
Bbiffop продольного шага: 1
(CALL MASS ? t ;
CALL WALL V —----------
(call output)^
" I
( CALL SOEFF
' i "
( CALL SOLVE V;
CONST
Вдавятся величины констант
BEGIN
Ввод шналрных профилей ЩнмЬление величины
Д&пали ctt.окончание дцсёграш^/
LENGTH Вычисление у1
ENTRN
вычисление т" и ib'E'дця свободной границы и линии симметрии
PRE
Вводятся величины
MASS
Вводятся величины т’‘ на стенке
Детали см. окончание диаграммы
OUTPUT Печатаются результаты
Детали см. продолжение диаграммы
SOLVE j
Решаются алгебраические уравнения последовательной подста-\ ___________новной________________j
Рис. П1-3.1. Начало.
Подпрограмма LENGTH. Предназначена для вычисления характерной толщины слоя в соответствии с определением разд. 2.8-3. Вначале разыскиваются максимальное и минимальное значения скорости, затем обследуются области вблизи обеих границ. В задачах ламинарного движения можно обойтись без этой подпрограммы.
