ГОСТ 25645.301—83 Стр. 81
н готовят соответствующие массивы
2.4. Подпрограмма TELA обрабатывает кнформацзю о типовых поверхностях
информации обо всех поверхностях второго порядка ИСЗ
2.5. Уравнение поверхности второго порядка относительно осей Ххi имеет вид
(i.* —1,2,3).
(2)
Р.
где Ид Ai, Л,* — коэффициенты уравнения второго порядка, вычисляемые по значениям параметров, приведенных в
табл. 2.
2.6. Границы области D поверхности второго порядка рассчитывают в системе координат Уу у. ориентированной от
носительно системы координат Xxt направляющими косинусами aJt углов между осями Х( п у/ н координатами х 0(
точки Y.
26.1. Для задания области
D
используют:
в прямоугольной системе координат (у/) — переменные(j>i, рг);
в полярной системе (ф, г, у») — переменные (<р, г);
в цилиндрической системе (<р.
г
, уа)—переменные («р, р,);
в
сферической системе (ф, в,
г)
—
переменные (Ф,
»).
В указанных системах координат почти для всех поверхностей второго порядка ИСЗ, встречающихся в практике,
форма области D простейшая: прямоугольная с пределами V, ипо внешней псрсмсзвой (у, и <р) и с пределами Yi
и
Ya по внутренней переменной (х/я. г, у* и в ).
2.G.2. Для обозначении систем координат вводят признакравный:
1 — для прямоугольной;
2 — для цилиндрической;
3 — для полярной;
4 — для сферической.
При £=. t к £= »3 в расчетах учитывают только те участки поверхности,для которых£fr>0, при Р<*2 зада
ют у, = 0.
2.7. На поверхностях второго порядка могут быть области исключения, имеющие сложную границу. Области исклю
чения задают в соответствии с п. 2.2.
2.7.1. Для задания области исключения используют два массива:
NPD — массив, элементы которого равны покерам областей исключении;
MD2 — массив, порядок расположения элементов которого соответствует расположению элементов NPD, а значения
раины номерам поверхностей второго порядка, к которым принадлежит область исключения.
Границы областей исключения данной поверхности могут быть односвязными, мкогосвязными. пересекаться друг с
другом и с границей поверхности.
2.8. Все поверхности второго порядка могут быть пористыми. Коэффициент пористостиMDI — часть молекул газа,
проходящая через поверхиость, не изменяя векгора скорости. Оставшаяся часть молекул газа (1—MDI) отражается от по
верхности, как от сплошной стенки.
2.9. Поверхность ИСЗ окружаюг условной контрольной поверхностью (КП), имеющей форму параллелепипеда.
2.10. Выбор случайных значений координат точки начала движения к компонент вектора скорости молекул набегаю
щего потока газа производят методом Мойте—Карло на поверхности КП с помощью функций распределения плотности
вероятности.
2.11. Для поверхности второго порядка проводят анализ решения системы двух уравнений:
уравнения поверхности второго порядка по формуле (2) и уравнения первого порядка траектория движения молеку
лы газа по формуле
flT ju + t- k i* •~ (3)
где x ,i,\ i — хоординаты точки начала движения и проекции единичного вектора скорости молекулы газа;
t
— расстояние, пройденное молекулой газа;
?. — индекс, значение которого равно номеру точки начала движения молекулы газа.
В результате решения определяют, столкнулась молекула с поверхностью или нет.
2.12. При столкновении молекулы газа с поверхностью ИСЗ проводят выбор случайных значений параметров закона
отражения молекулы от поверхности. Закон отражения молекулы от поверхности задают значением исходного признака
2.12.1. При зеркальном отражении (LOT а —-1<) компоненты скорости | ./отраженной молекулы газа равны
2(е+ |-п,)ль 4—1,2.3. ,(4)
где п/ — компоненты единичного вектора внутренней нормали к поверхности;
«+» и «—» — индексы, обозначающие соответственно сталкивающиеся к отраженные от поверхности молекулы.
2.12.2. При зеркально-диффузном отражении молекулы газа от поверхности (LO T=l) задают исходные данные:
коэффициент г IN — часть диффузно отражающихся молекул;
функцию <р,;
зависимость коэффициента аккомодации энергии а, от угла столкновения молекулы газа вс поверхностью при
характерном значении энергии столкновения молекул £ , . , вычисляемом по формуле
£ ,.- 0 ,8 6 2 9 10~ ‘ S i T „ ,(5)
где S . - отношение скорости набегающего потока газа к наиболее вероятной скорости молекул набегающего потока;
Гд,-- температура набегающего потока газа.
Параметры <р, и о, имеют п программе идентификаторы соответственно FIE и SIGF,
Функцию ф /, описывающую изменение о,, и зависимости от энергии молекулы £ + при характерном значении угла
столкновения молекулы 6 +.( 0 ° < б.,* < 9 0 “). вычисляют по формуле
где о ,» = 1 — значение коэффициента о . для случая равновесного отражениямолекул от поверхности
6
3
jk
. и «
°«Я-—
?,<*, . е : .)
(<•)