Существующие методы численного моделирования электромагнитного поля в среде обладают, к сожалению, каждый своими недостатками. Авторы поставили себе задачу проанализировать наиболее популярные методы. В качестве модельной задачи авторами рассматривается линза Люнеберга. В данной работе авторы рассматривают метод конечных разностей во временной области, программное средство openEMS и его применимость к задачам численного моделирования распространения электромагнитных волн в среде на примере сферической линзы Люнеберга. Благодаря своей простоте и широким возможностям метод конечных разностей во временной области (Finite-Difference Time-Domain method, FDTD) применяется для решения широкого спектра задач. Существует достаточно большое количество программных инструментов, как с открытым исходным кодом, так и проприетарных, позволяющих производить расчёт этим методом. Программный комплекс openEMS является набором функций для MATLAB или Octave, с помощью которых можно произвести расчёт характеристик электромагнитного поля методом EC-FDTD в декартовых или цилиндрических координатах. Программный комплекс openEMS является бесплатным и имеет открытый исходный код. Поддерживаются параллельные технологии вычисления (MPI). В данной работе на примере моделирования прохождения электромагнитных волн сквозь сферическую линзу Люнеберга показан процесс работы с openEMS, его установка и настройка, а также даны общие сведения о работе метода FDTD и алгоритма Йи. Приведён пример работы алгоритма Йи. Показан способ анализа и визуализации результатов моделирования с помощью программы ParaView. Приведён исходный код скрипта для моделирования. Исследованы возможности openEMS и метода FDTD при моделировании распространения электромагнитных волн в среде.