Научная литература
booksshare.net -> Добавить материал -> Физика -> Адзерихо К.С. -> "Лекции по теории переноса лучистой энергии" -> 29

Лекции по теории переноса лучистой энергии - Адзерихо К.С.

Адзерихо К.С. Лекции по теории переноса лучистой энергии — БГУ, 1975. — 192 c.
Скачать (прямая ссылка): lekciipoteoriiperenosaluchistoyenergii1972.djvu
Предыдущая << 1 .. 23 24 25 26 27 28 < 29 > 30 31 32 33 34 35 .. 44 >> Следующая


П -f і + 1

2

п + і ’

Таким образом, неизвестные постоянные An можно определить из следующей системы алгебраических уравнений:

Функция q (т), полученная методом Мензела—Сена для т = 5, имеет вид

q(т) = 0,71044936 — 0,27186451 E2 (т) -f 0,41860309 Е3(т) —

— 0,31239946 Е4(т) + 0,13433965 E6 (т). (5.27)

Расчет по (5.27) дает погрешность лишь IO-5—10~4% для всего интервала изменения оптических толщин.

В заключение отметим, что если в качестве пробных функций выбирать используемые остаточные разности, то мы приходим к вариационным методам, которые подобны методам, рассмотренным в предыдущем параграфе.

т

= _А,_________L

J I о м I

п + 2 п + 3

(5.26)

124
§ 4*. Метод Монте-Карло

Изложенные выше методы иг позволяют решать задачи теории переноса .излучения, в которых наряду с пере* распределением по частотам и неизотропностью актов рассеяния необходимо учитывать такие важные факторы, как неоднородность среды и ее произвольную трехмерную конфигурацию. В настоящее время одним из основных методов решения подобных задач является метод Монте-Карло, или метод статистических испытаний. Этот метод, являясь чисто численным и индуктивным,' принадлежит к вероятностным методам (об этом говорит и само название метода). Основная идея метода Монте-Карло заключается в моделировании реального физического процесса (в частности, траектории фотона), которое производится с использованием таблиц случайных чисел. При моделировании траектории фотона последовательно учитываются акты взаимодействия фотона с элементарными объемами вещества вплоть до исчезновения рассматриваемого фотона (до выхода его из среды или поглощения его средой).

При решении задач теории переноса излучения методом Монте-Карло с помощью прямого моделирования необходимо сначала выбрать некоторое число N пробных фотонов. Прослеживая историю распространения этих фотонов в среде и статистически усредняя результаты, можно отыскать искомые характеристики выходящего излучения. Так как при этом относительная ошибка расчетов пропорциональна IlVN, то число N должно выбираться достаточно большим. Из фотонов выделяется t-й фотон. Согласно таблицам случайных чисел (см., например, [44]), выбирается некоторый элементарный телесный угол (і = I, 2, ... , N\ / = 1, 2, ... ... , ti), в котором распространяется і-й фотон. Набор случайных чисел можно получить с помощью генератора случайных чисел и таблиц псевдослучайных чисел. Последние достаточно случайны, хотя и получаются по вполне определенному закону. Например, в [45] они определяются как остаток при делении 23 ап на (IO8 + -г 1), где п = 0, 1, 2, ... — некоторое произвольное

восьмизначное число. Набор чисел, получаемых таким образом, имеет повторяемость через 5882352 числа, и поэтому его можно считать достаточно случайным. Если эти

12*
числа разделить на IO81 получается набор псевдослучайных чис^л между 0 и 1.

После выбора dQ\l) необходимо определить величину пробега /-го фотона R\l) (I = 1,2......т) с учетом оптиче-

ских свойств среды. Обычно при определении R\l) используют закон Бугера. При этом следует учитывать возможность выхода кванта из среды или его поглощения. Если R0 является граничной поверхностью среды, то при R}1) ^R0 фотон покидает среду и его историю можно считать оконченной, а при R}1) < R0 фотон или испытывает поглощение (конец истории), или распространение его в среде продолжается. Таким образом, перед каждым актом рассеяния (или поглощения) состояние фотона можно охарактеризовать некоторой совокупностью чисел

Sh — Rh I ^.fc). k = 1, 2......kg,

где под Xh подразумевается вероятность выживания кванта, причем при k<k0 Xh = 1, а при k = kg Xh = 0 (факт исчезновения фотона). Набор этих чисел образует цепь Маркова, которую можно охарактеризовать некоторой плотностью 1Ir {sh, sh+1). Функция 1Ir определяет вероятность перехода фотона из k-ro в (k H- 1)-е состояние. Принимая во внимание закон Бугера, функцию 1Ir можно представить в следующем виде [24]:

при Xh = Xm = 1

ч (Vh+i) = Kp !*)(* + о)

при Xh = I, Х*+1 = 0

V{sh,sh+1}=l- a

х + а и + а при Xh = 0

1JrK. sfc+1} = о.

Функция" плотности 1F (sfc, Sh+1} позволяет находить распределение фотонов при (k 1)-м столкновении, если известно распределение при k-іл столкновении. Такое последовательное определение распределения фотонов аналогично обычному процессу итераций и при достаточно большом количестве шагов (k0, п, т= 10 — 100) эквивалентно интегральному уравнению переноса излучения. Кроме того, функция 4{sh, sft+1} позволяет определить начальное распре-

126
деление фотонов в среде, от которого существенно зависит сходимость рассматриваемого метода.

Так, начальное распределение фотонов можно задать в виде

OO

F(S1)= J^s11 S0) ео (г — х1, I)dx =

О

OO

= J (х + О) е-(х+а,\ (г — л: I1 I) dx. (5.28) о

Для точечного коллимированного источника излучения е0(г, 1) = б (г — г0) б (1 —10)

и поэтому

F (S1) = (х + <т) e-(*+0,(r-V!o6 (I —10).

Следует отметить, что применение метода Монте-Карло связано с большим объемом вычислительной работы, который доступен лишь для очень быстродействующих машин. Поэтому в настоящее время уделяется большое внимание вопросам ускорения сходимости метода Монте-Карло. В сильной степени ускорение метода достигается использованием па промежуточных этапах некоторых простейших аналитических решений. Сходимость существенно зависит и от числа выбираемых фотонов. В пределах задаваемой точности можно выбрать оптимальное число систем и, кроме того, выделить наиболее существенные области. Для выделения последних обычно пользуются теоремой перестановки источников [38]: если р — решение интегрального уравнения
Предыдущая << 1 .. 23 24 25 26 27 28 < 29 > 30 31 32 33 34 35 .. 44 >> Следующая
Реклама

c1c0fc952cf0704ad12d6af2ad3bf47e03017fed

Есть, чем поделиться? Отправьте
материал
нам
Авторские права © 2009 BooksShare.
Все права защищены.
Rambler's Top100

c1c0fc952cf0704ad12d6af2ad3bf47e03017fed