МЕТОД МОНТЕ-КАРЛО ДЛЯ РАСЧЕТА ХАРАКТЕРИСТИК СВЕТОВЫХ ПОЛЕЙ В МОРСКОЙ ВОДЕ

  • Д. И. Глуховец Институт океанологии им. П. П. Ширшова РАН
DOI: 10.29006/1564-2291.JOR-2024.52(3).1
Ключевые слова: метод Монте-Карло, уравнение переноса излучения, коэффициент диффузного ослабления света, коэффициент диффузного отражения, морская вода

Аннотация

Развитие численных методов решения интегро-дифференциального уравнения переноса излучения остается актуальной задачей на сегодняшний день. Среди них можно выделить востребованный в различных нишах современной оптики океана метод Монте-Карло. Цель данной работы – доступное и лаконичное изложение основ метода прямого Монте-Карло моделирования световых полей в морской воде, сопровождающееся подробным описанием способа его программной реализации. Изложены основы метода, описаны процедуры выбора типа взаимодействия, величины длины свободного пробега и направления движения фотонов. Рассмотрен простой случай, соответствующий бесконечно удаленному точечному источнику неполяризованного света, отсутствию влияния атмосферы, гладкой поверхности раздела и отсутствию стратификации первичных гидрооптических характеристик. При этом использованы реалистичные значения показателей поглощения и рассеяния, рассчитанные в соответствии с моделью вод первого оптического типа для концентрации хлорофилла 1 мкг/л, и сильно вытянутая индикатриса Хеньи-Гринстейна с параметром g = 0.95. Выполнен учет френелевского отражения света от поверхности раздела «воздух – морская вода». Относительные ошибки значений коэффициента диффузного ослабления нисходящего потока излучения Kd и коэффициента диффузного отражения R, рассчитанные в спектральном интервале 400–700 нм с использованием 106 фотонов, в сравнении с результатами HydroLight составили 1.5 % и 0.4 % соответственно. Расчет спектров на одном ядре мобильного процессора Intel Core i5-8250U 2017 г. в среде MATLAB занимает 6 минут. Выполнена оценка выбора оптимального количества фотонов, необходимого для получения искомых величин с заданной точностью. Реализованный метод полезен для знакомства с основными принципами, применяемыми для численного решения уравнения переноса излучения в морской воде статистическими методами, и используется в курсе «Оптика океана», читаемом автором студентам 4 курса кафедры термогидромеханики океана МФТИ.

Опубликован
2024-11-06
Раздел
Физика океана и климат