расчет интенсивности ультрафиолетовой радиации методом Монте-Карло

  • Автор темы Komandor
  • Дата начала
Статус
Закрыто для дальнейших ответов.
K

Komandor

#1
Нужна помощь в написании программы, производящей расчет интенсивности ультрафиолетовой радиации методом Монте-Карло. Программа написана, но результаты не совпадают с эталоном. Если кто занимался данной тематикой или может в ней разобраться, буду рад Вашей помощи, естественно, не за спасибо. Программу и документацию по проблеме и методу могу предоставить в любой момент.
 

Kmet

Java Team
25.05.2006
1 036
8
#2
метод Монте-Карло относится к стохастическим. говорить в таком случае о совпадении значений некоректно. при большом количестве итераций результат должне стремится к эталону, не более того
 
K

Komandor

#3
Про стохастичность метода Монте-Карло я знаю, но несовпадение значений очень грубое. К тому же увелечение количества экспериментов с одними и теми же начальными параметрами не ведут к увеличению точности.
Вы вроде знаете метод... Сможете ли помочь?
 

Kmet

Java Team
25.05.2006
1 036
8
#4
метод то я знаю... но вот каким образом вы применяете его к измерению радиации понятия не имею
 

morpheus

скриптописец
07.08.2006
3 915
1
#5
<!--QuoteBegin-Kmet+8:10:2007, 17:23 -->
<span class="vbquote">(Kmet @ 8:10:2007, 17:23 )</span><!--QuoteEBegin-->метод Монте-Карло относится к стохастическим. говорить в таком случае о совпадении значений некоректно. при большом количестве итераций результат должне стремится к эталону, не более того
[snapback]81012" rel="nofollow" target="_blank[/snapback]​
[/quote]
*[ B) пашёл нервно курить кофе... ]
 
K

Komandor

#6
По применению метода Монте-Карло к расчету радиации.
Представим определенный объем воздуха на верхнюю границу которого падает излучение. Расчет ведется для одного фотона -- прослеживаем его путь в атмосфере. Для фотона известно -- его энергия, углы направления и параметры атмосферы.
Считается, что фотон пролетает определенное растояние -- длину свободного пробега -- и рассеивается на молекуле. Для него рассчитываются новые углы и очередная длина свободного пробега и энергия. На каждом рассеянии фотон теряет энергию. Если она становиться 0 или меньше, фотон убирается. И так далее.
Метод Монте-Карло используется для расчета длины свободного пробега. Формулы есть, литература есть. Углы расчитываются по вероятностной кривой, там все просто. Вес вообще можно убрать. Основная проблема именно в расчете длины свободного пробега.
 
K

Komandor

#7
Мое описание ужаснуло? Ну хоть кусок по методу Монте-Карло посмотрите, остальное я уж как нибудь сам...
 

Kmet

Java Team
25.05.2006
1 036
8
#8
Вес?! Каким весом можно принебречь?!
Могу предположить что следующие ошибки в реализации:
1. не учитывается возрастаение плотности воздуха.
2. как я понял рассчитывается радиация для некоторого разреза расположеного на растоянии Н от верхней границы. В таком случае рассчитать длину свободного пробега не достаточно да и не нужно, определющим значением является длина проекции свободного пробега на вертикальную ось
 
K

Komandor

#9
Про вес прошу прощения. Вес -- это энергия фотона, которая уменьшается при каждом столкновении. Весом можно принебречь, заставив фотон летать вечно, пока не вылетит за верхнюю или нижнюю границу объема.

Кроме длины свободного пробега определяются углы к горизонту -- по верояностной кривой (индикатрисе рассеяния). После можно найти точку следующего столкновения и опять рассчитать длину свободного пробега и углы. и тд
 
Статус
Закрыто для дальнейших ответов.