Числовое моделирование электромагнитного поля в теле человека давно перешло из академической дисциплины в рабочий инструмент разработчика сотовых модулей, радиочастотной медицинской аппаратуры и систем беспроводной передачи данных. Метод конечных разностей во временной области, известный как FDTD, решает уравнения Максвелла напрямую, продвигая поле во времени малыми шагами, и выдаёт то, чего не способна дать ни одна аналитическая формула, точную картину распределения поля в слоистой неоднородной среде с реальной геометрией. Именно эта картина определяет, на какой глубине поглощается энергия, где возникают максимумы удельной поглощённой мощности и как конструкция антенны меняет нагрев тканей. В статье разобраны практические стороны такого расчёта, от выбора шага сетки до получения карты SAR.
Сетка Яи и шаг по времени из условия Куранта
В основе метода лежит разбиение пространства на прямоугольные ячейки, в которых напряжённость электрического поля и магнитная индукция хранятся со сдвигом на половину шага друг относительно друга и по координатам, и по времени. Такое расположение, предложенное Яи, превращает уравнения Максвелла в разностные соотношения второго порядка точности, где каждая компонента поля вычисляется из значений соседних компонент на предыдущем полушаге. Производные по времени заменяются центральными разностями, что даёт вторую точность одновременно по всем переменным. На каждую ячейку приходится шесть компонент поля, три электрических и три магнитных. Практическое следствие простое: сетка из трёх миллионов ячеек с одинарной точностью и шестью компонентами требует порядка 72 Мбайт памяти только под поля, а с материалами, источниками и служебными массивами счёт легко уходит за 300 Мбайт. Отсюда первый инженерный вывод, каждая лишняя тысяча ячеек в поперечнике оплачивается кубическим ростом общего объёма задачи, поэтому экономию начинают с правильного выбора шага, а не с оптимизации кода.
Явная схема устойчива не при любых параметрах. Связь между шагом по пространству и шагом по времени задаёт условие Куранта, для трёхмерной равномерной сетки оно записывается как Δt ≤ Δx/(c√3), где c скорость света в вакууме. Для шага сетки Δx = 1 мм предельный временной шаг составляет 1e-3/(2,998e8 × 1,732) ≈ 1,92 пс. Если взять шаг крупнее, решение не сходится, поле начинает расти экспоненциально и через несколько тысяч итераций задача разваливается. Для рабочей частоты 2,45 ГГц период колебаний равен 408 пс, значит при Δt = 1,92 пс на один период приходится примерно 213 шагов, что даёт хороший запас по точности. В расчётах с импульсными источниками общее число шагов достигает нескольких тысяч, для установившихся режимов счёт ведут до затухания переходного процесса, поэтому реальное время вычислений определяется именно произведением числа ячеек на число шагов.
Электрические параметры кожи, жира и мышц на рабочих частотах
Биологические ткани ведут себя на СВЧ как диэлектрики с высокими потерями, и их параметры сильно зависят от частоты. На 2,45 ГГц мышца имеет относительную диэлектрическую проницаемость около 52,7 и проводимость около 1,74 S/m, жировая ткань проницаемость около 5,3 и проводимость около 0,10 S/m, сухая кожа проницаемость около 38 и проводимость около 1,44 S/m, кровь проницаемость около 58 и проводимость около 2,54 S/m. На 900 МГц проницаемость мышцы возрастает примерно до 55, а проводимость падает до 0,97 S/m, на 5,8 ГГц она снижается до 48 при проводимости около 4,9 S/m. Такую частотную зависимость задают дебаевскими и коул-коуловскими моделями по базе параметров, собранной в ходе многолетних измерений образцов тканей. Длина волны в мышце на 2,45 ГГц составляет λ = c/(f√εr) = 3e8/(2,45e9 × 7,26) ≈ 16,9 мм, поэтому для разрешения поля правило десятой доли длины волны требует шага сетки не крупнее 1,7 мм в среде с максимальной проницаемостью. Именно мышца, а не воздух, диктует размер ячейки, и это удорожает счёт в разы по сравнению с расчётами в свободном пространстве.
Глубина проникновения поля в слоистой модели тканей
Для плоской волны в однородной среде глубину проникновения считают через коэффициент затухания α = ω√(εμ/2)[√(1+(σ/ωε)²) − 1]^0,5, где ω угловая частота, ε абсолютная диэлектрическая проницаемость, μ магнитная проницаемость, σ проводимость. Подстановка параметров мышцы на 2,45 ГГц даёт следующее, отношение σ/ωε = 1,74/(1,539e10 × 8,854e-12 × 52,7) ≈ 0,24, корень из разности даёт 0,17, множитель ω√(εμ/2) равен примерно 263 1/м, итоговое затухание α ≈ 45 1/м. Поле ослабевает в e раз на глубине 1/α ≈ 2,2 см, а плотность мощности, пропорциональная квадрату поля, падает в e раз примерно на 1,1 см. В жировой ткани картина иная, там сочетание малой проводимости и малой проницаемости даёт заметно меньшее поглощение, поэтому жир остаётся почти прозрачным слоем, а основная энергия выделяется в коже и мышце. На границе слоёв возникает отражение из-за разницы волновых сопротивлений, которые для приближённой оценки равны η ≈ 377/√εr, то есть около 164 Ом для жира и около 52 Ом для мышцы. Коэффициент отражения по амплитуде Γ = (52 − 164)/(52 + 164) ≈ −0,52, по мощности это порядка 27 процентов, и в слое жира формируется стоячая волна с максимумом напряжённости у границы. Отсюда характерное для расчётов явление, локальные максимумы нагрева на границах слоёв, которых однородная модель в принципе не покажет.
Удельная поглощённая мощность, пределы SAR и типичные результаты
Поглощение энергии характеризует удельная поглощённая мощность, определяемая как SAR = σ|E|²/(2ρ), где ρ плотность ткани, для мышцы её берут около 1090 кг/м³. Подстановка среднеквадратичного значения напряжённости 100 В/м и проводимости 1,74 S/m даёт SAR = 1,74 × 10⁴/(2 × 1090) ≈ 8,0 Вт/кг. Это в четыре раза выше общественного предела 2 Вт/кг, нормируемого с усреднением по 10 г ткани, и в пять раз выше американского предела 1,6 Вт/кг с усреднением по 1 г. Среднее берут по шестиминутному интервалу, потому что тепловая инерция тканей сглаживает быстрые колебания мощности. В расчёте FDTD локальный SAR вычисляется в каждой ячейке по локальному значению поля, затем усредняется по кубу или сфере нужной массы, и максимальное значение этого среднего сравнивается с пределом. Такой порядок исключает ситуацию, когда одна горячая точка размазана усреднением и ускользает от контроля, а также показывает, насколько глубоко под кожу уходит зона заметного поглощения. Разница между усреднением по 1 г и по 10 г измеримо велика, одна и та же горячая точка при усреднении по меньшему объёму даёт значение на 30-60 процентов выше, и устройство, проходящее по одной норме, может не пройти по другой. Именно поэтому в отчётах по расчёту всегда указывают, по какой массе велось усреднение, иначе число само по себе теряет смысл и не позволяет сравнивать результаты между собой.
Поглощающие границы и объём вычислительной области
Расчётная область должна имитировать открытое пространство, иначе волны отражаются от краёв сетки и искажают решение. Стандартное решение, поглощающие слои PML, которые гасят падающие волны с отражением на уровне минус 60-80 дБ при толщине 8-12 ячеек. Оценка объёма задачи ведётся для фантома головы диаметром 18 см. Плюс воздушный зазор 2 см до антенны, плюс два слоя PML, итоговый габарит около 24 см, при шаге 1,7 мм это примерно 141 ячейка по каждой оси и около 2,8 миллиона ячеек на всю область. Несколько тысяч временных шагов дают доли миллисекунды физического времени, а на счёт уходит от десятков минут на видеокарте до часов на одном процессорном ядре. Существуют гибридные схемы, где основная часть области считается грубой сеткой, а зона интересов уточняется вложенной сеткой с шагом вдвое или вчетверо меньшим, это экономит память без потери точности в критичных местах. Экономия на ячейках вдали от источников почти не меняет результат, поэтому грубую сетку делают крупной, а уточнение ведут только вокруг антенны и тканей.
Отдельная тонкость, согласование шага сетки с толщиной слоёв. Кожа толщиной 2 мм при шаге 1,7 мм попадает в одну-две ячейки, и её вклад в поглощение оценивается грубо, поэтому в критичных задачах шаг в области кожи уменьшают до 0,5-0,8 мм, иначе ошибка в локальном SAR у поверхности достигает десятков процентов. Воксельные модели стандартного разрешения 1-2 мм согласованы с этим ограничением лишь частично, и ответственный расчёт всегда проверяет, сколько ячеек приходится на самый тонкий слой в зоне интересов.
Для сотового модуля, работающего в диапазонах 900 МГц и 1,8-2,6 ГГц, расчёт показывает одну закономерность, максимум SAR локализуется в ухе и прилегающих тканях на глубине первых миллиметров, а до мозга доходит лишь малая доля энергии. При увеличении зазора между антенной и головой с 5 до 15 мм локальный SAR снижается в 3-5 раз, поскольку в ближней зоне поле спадает почти квадратично с расстоянием. Форма корпуса и положение согласующих элементов меняют картину не меньше, чем мощность передатчика, сдвиг точки питания на несколько миллиметров перестраивает распределение тока по корпусу и снижает пиковые значения без потери связи. Разница между диапазонами заметна на глаз, на 900 МГц глубина проникновения примерно вдвое больше, чем на 2,45 ГГц, поэтому там энергия распределяется в большем объёме тканей, а пиковый SAR при той же излучённой мощности выходит ниже. Эти закономерности проектировщик использует прямо, выбирая компромисс между положением антенны, её полосой и запасом по пределам. Для импульсных воздействий картина меняется, SAR усредняется по времени, и при скважности 0,01 пиковые значения поля в тканях в десять раз выше средних при том же среднем поглощении. Расчёт учитывает это автоматически, поскольку сетка видит мгновенные значения поля, а пределы сравниваются со средним за окно, поэтому для импульсных режимов запас по пиковой мощности задают отдельно от среднего.
Сравнение FDTD с другими численными методами
Метод конечных элементов решает ту же задачу в частотной области и выигрывает на гладких геометриях с точным согласованием сетки, но для широкополосного расчёта ему нужен отдельный прогон на каждую частоту. FDTD за один прогон с импульсным источником даёт отклик во всём диапазоне сразу, достаточно преобразования Фурье временных реализаций, поэтому для оценки SAR в нескольких диапазонах одновременно он выигрывает в разы. Метод моментов удобен для антенн над однородным телом, но слоистые неоднородные модели переводят его в разряд громоздких. Лучевые методы работают выше 10-20 ГГц, где длина волны много меньше деталей тела, однако на 900 МГц и 2,45 ГГц они теряют дифракцию и интерференцию, которые как раз формируют горячие точки. На практике связку строят так, лучевые методы дают быструю оценку на высоких частотах, FDTD уточняет критичные случаи, конечные элементы проверяют отдельные узлы с гладкой геометрией. Выбор метода определяет не только скорость, но и достоверность результата, поскольку каждый подход видит свою часть физики и слеп к остальной.
Порядок проведения расчёта от модели до карты SAR
Полный цикл моделирования укладывается в шесть последовательных шагов:
- построение воксельной модели тела или фантома с разрешением, согласованным с шагом сетки;
- присвоение каждой ячейке электрических параметров ткани на рабочей частоте;
- задание источника в виде гауссова импульса для широкополосного расчёта или гармонического сигнала для одной частоты;
- итеративное решение разностных уравнений с шагом Куранта до затухания переходного процесса;
- преобразование Фурье временных реализаций поля и расчёт локального SAR с усреднением по 1 г и по 10 г;
- сравнение максимума SAR с пределами и корректировка конструкции антенны или радиомодуля.
Добросовестный расчёт всегда сверяют с измерением на фантоме, заполненном жидкостью с параметрами тканей, и расхождение в пределах 10-15 процентов считают хорошим результатом. При аккуратной постановке задачи метод FDTD предсказывает и глубину проникновения, и расположение горячих точек с точностью, достаточной для проектирования безопасных радиочастотных устройств задолго до первого прототипа.