Тепло и массообмен в пограничных слоях - Патанкар С.
Скачать (прямая ссылка):
В примере, приведенном ,в § П1-5, все необходимые величины считываются с перфокарт. Это удобно, если имеются экспериментальные профили. С другой стороны для построения начальных профилей можно использовать алгебраические выражения. В этом случае для показателей (3 и у принимаются одинаковые значения, равные 0,143 (профиль одной седьмой), т. е. хорошо известные степенные профили, характеризующие турбулентное пристеночное течение. Для ламинарного течения выбираются другие профили. Соотношения для величин скольжения, используемые в примере, взяты из § 2.5, а значения (i|)e—грj) и и подсчитываются по формулам (1.1-2) и (2.1-1).
Подпрограмма OUTPUT. Содержит команды на печатание результатов. Остальная программа не нуждается в этой подпрограмме, и поэтому вычислитель обладает свободой выбора печатать по своему желанию любую информацию, относящуюся к расчету.
Целочисленную переменную INTG, учитывающую количество интегрирований, можно успешно использовать для команды иа воспроизведение результатов после желаемого числа интегрирований. В примере, приведенном в П1-5, предусмотрено печатание величин абсциссы х и соответствующих профилей скорости и температуры через каждые пять интегрирований.
П1-4. КОНТРОЛЬ ПРАВИЛЬНОСТИ СОСТАВЛЕНИЯ ПРОГРАММЫ
Приведенное нами выше описание деталей вычислительной программы и сводка отдельных подпрограмм, представленная в § П1-5, составляет информацию, достаточную для самостоятельного выполнения машинного счета. Однако после приобретения общего понимания программы расчетчик может не пожелать использовать всю эту информацию целиком всякий раз при составлении заново вычислительной программы для решения конкретной задачи.
В этом случае нужна уверенность в том, что соблюдены следующие десять операций (они выражают ту же самую информацию в сжатой и удобной для использования форме).
1. Подпрограммы COEFF, SLIP, SOLVE и READY сохраняются без изменений:
2. В подпрограмме .MAIN определение длины шага (xD—Хс) и условия для окончания вычислений выбираются в соответствии с поставленными целями.
3. Выбирается формула для эффективной вязкости. В случае приемлемости гипотезы о пути смешения подпрограммы
VEFF, LENGTH, ENTRN, WALL, WF1, WF2
сохраняются неизменными. Для любой другой гипотезы течения эти подпрограммы МО' дифицируются соответствующим образом.
4. Выбираются величины KRAD, NEQ, KIN, КЕХ, N и вводятся в подпрограмму BEGIN совместно с начальными значениями х н начальными профилями. Выясняется возможность согласования типа границы (т. е. KIN и КЕХ) с характером движения в рассматриваемой задаче. Если нужны изменения, то величины KL\ пли КЕХ вводятся в подходящем месте программы. Следует изменить величины [3 и у BEGIN, есля течение не турбулентно. Нужно, насколько это возможно, ограничить все команды READ, относящиеся к этой подпрограмме.
5. Подпрограмма SOURCE приводится в соответствие с используемыми ©-уравнениями. Если граница совпадает со стенкой, то подпрограммы Р'ВС и MASS рассматриваются таким образом, чтобы обеспечить задание желаемых Ф-граничных условий1 и интенсивности массоперепоса через стенку.
6. Проводится подготовка для вычислений соответствующих свойств жидкости в подпрограммах DENSTY и VISCO.
7. Подготавливаются подпрограмма PRE для задания величины продольного градиента давления в соответствии с условиями задачи и подпрограмма RAD, вводящая информацию о геометрии 1-границы.
8. Нужная для расчета информация о постоянных вводится в CONST.
9. Предусматривается вывод и воспроизведение (печатание) нужных результатов в OUTPUT.
10. Подготавливается перфокарта с данными с указанием последовательных ко манд READ по всей программе.
П1-4. ВЫЧИСЛИТЕЛЬНАЯ ПРОГРАММА
В этом параграфе мы приводим полный текст вычислительной программы. Подпро-граммы, за исключением MAIN, которая приводится первой, расположены в алфавитном порядке.
MAIN
COMMON/GEN/PEI,AMI,AME,DPDX,PREF (2) ,PR(2) ,Р (2),DEN,AMU, IU,XD,XP,XL,DX,INTG,CSALFA 1/I/N,NP1,NP2,NP3,NEQ,NPH,КЕХ,KIN,KASE,KRAD I /В/В ETA,GAMA (2) ,TAUI,TAUE,AJI (2) ,AJE (2) ,IND [ (2) ,INDE (2)
1IV/V (43) ,F (2,43) ,R (43) ,RHO (43) ,OM (43) ,Y (43) iI/C/SC (43) ,AU(43) ,BU (43) ,CU (43) ,A (2,43) ,B (2,43) ,C (2.43) COMMON/L/AK,ALMG !G Продолжение INTG = 0 XL = 20.
CALL CONST CALL BEGIN AMI = 0.
AME = 0.
Переход к 25 15 CALL READY
25 Продолжение INTG= INTG+ 1 CALL LENGTH CALL ENTRN С Выбор предыдущего шага FRA --=.05
DX = FRA*PEI/(R(1)*AMI—R(NP3)*AME)
IF(DX.GT..5*Y(NP3))DX = .5*Y(NP3)
XD=XU+DX
77 Продолжение
CALL PRE(XU,XD,DPDX)
IF(KASE.EQ,2) Переход к 26 IF (KIN.EQ. 1) CALL MASS (XU,XD,AMI)
1 Обычно трактовки граничных условий в подпрограмме SLIP оказывается достаточно, но при рассмотрении сложных граничных условий необходима некотопая модификация. ' F
