Вычисление температуры нагрева головы человека при пользовании сотовым телефоном
Моделирование эффекта поглощения мощности в голове человека при пользовании сотовым телефоном. Численный анализ электромагнитного поля и поглощаемой мощности внутри головы пользователя. Расчет соответствующего этому поглощению повышение температуры.
Рубрика | Физика и энергетика |
Вид | статья |
Язык | русский |
Дата добавления | 04.11.2018 |
Размер файла | 1,3 M |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
ВЫЧИСЛЕНИЕ ТЕМПЕРАТУРЫ НАГРЕВА ГОЛОВЫ ЧЕЛОВЕКА ПРИ ПОЛЬЗОВАНИИ СОТОВЫМ ТЕЛЕФОНОМ
А. А. Курушин
Широкое распространение сотовой телефонной связи усилило беспокойство, связанное с облучением головы человека радиоволнами, излучаемыми антенной телефона при его использовании. Комплексные исследования физиологических изменений, происходящих в биологической ткани во время разговора по сотовому телефону, установили, что через десять минут разговора температура кожи человека в непосредственной близости от телефона возрастает на 0.5…2° [2]. В соответствии с требованиями американского стандарта, мощность, излучаемая антенной сотового телефона, не должна превосходить некоторого предельного уровня, при котором поглощаемая мощность внутри головы человека не превышает 1.6 мВт на 1г веса биологической ткани.
Удельная поглощаемая мощность 1.6 мВт/г - это характерная величина, которая имеет место при мощностях, излучаемых антеннами современных сотовых телефонов. Стремление увеличить чувствительность, дальность и устойчивость связи входят в противоречие с необходимостью уменьшать поглощаемую в голове мощность. Чтобы этого достичь, можно идти двумя путями: уменьшить мощность передатчика телефона, что может привести к увеличению вероятности сбоя связи (“звонок не доходит”), или разработать антенную структуру телефона так, чтобы уменьшить облучение головы, не уменьшая эффективности излучения антенны в дальней зоне.
За последние 10 лет антенны сотовых телефонов [6-8, 12] значительно усовершенствовались по характеристикам и усложнились по конструкции (рис. 1). При этом в некоторых моделях в качестве излучающей системы задействован и корпус телефона [11].
a) cпиральная антенна вне корпуса сотового телефона |
b) PIFA- антенна в корпусе сотового телефона |
c) планарная патчевая антенна сотового телефона |
Рис.1. Возможные модели сотового телефона и его антенны.
В данной работе рассматривается сотовый телефон с наружной спиральной антенной (рис. 1а), находящийся вблизи головы пользователя (talk position). Хотя эта модель телефона является несколько устаревшей, тем не менее, она позволяет определить основные качественные и количественные характеристики взаимодействия поля излучения телефона и головы пользователя. В модели используется фантом головы человека, который хорошо описывает биологическую структуру и электродинамические свойства реального тела. На основе численного анализа электромагнитного поля и поглощаемой мощности внутри фантома головы рассчитано соответствующее этому поглощению повышение температуры.
Характеристики воздействия радиоизлучения на биологическое тело
Для оценки воздействия радиоизлучения на организм человека можно применить, по крайней мере, два критерия: 1) критерий увеличения температуры тела при длительном воздействии облучения и 2) величину поглощенной мощности на единицу веса. Повышение температуры тела зависит от мощности, поглощаемой и рассеиваемой в виде тепла в теле, поэтому второй параметр - т. н. величина SAR (Specific Absorption Rate) - получил более широкое распространение [3-5]. Величина SAR в ячейке анализируемого пространства определяется как отношение мощности, поглощаемой в данной ячейке, к весу биологической ткани в ней [2]:
(1)
где у - проводимость биологической ткани, заполняющей ячейку, См/м,
E - напряженность электрического поля внутри ячейки, В/м
r - плотность ткани внутри ячейки, кГ/м3.
В 1948 г. H.H.Pennes [1] предложил математическую модель для передачи тепла в биологическом теле, обусловленной диффузией крови. Хотя были разработаны и более сложные модели процесса теплопередачи, этот подход был уточнен и всё еще используется в настоящее время. Модель описывает влияние тока крови на температуру ткани в стационарном состоянии. Поэтому источник и канал передачи тепла введены в уравнение теплопроводности:
, (2)
где r - плотность массы биологического материала, cp - теплоемкость биологического материала, l--- теплопроводность материала, MR - скорость генерации тепла в соответствии с метаболическими процессами, x--- скорость диффузии тепла, rb -плотность массы крови, cb - теплоемкость биологического материала, Ta - артериальная температура, SAR - удельная мощность поглощения. Каждое слагаемое в (2) имеет смысл и размерность мощности.
Если все механизмы теплопроводности в рамках простейшей диэлектрической модели головы свести к потерям, то из (2) можно получить упрощенное уравнение, в котором увеличение температуры связано с величиной SAR [2]:
, (3)
Где С - коэффициент теплопроводности [Дж/кГ/град],
DT - увеличение температуры в заданной точке тела, град,
Dt - время действия облучения, мин.
Таким образом, есть два способа определения SAR, основанные 1) на измерении или расчете увеличения температуры тела на заданном отрезке времени, и 2) на измерении или расчете напряженности электрического поля в заданной точке внутри фантома головы.
Для того, чтобы определить поле E в системе, состоящей из излучающего телефона и расположенного рядом биологического объекта, нужно решить систему уравнений для электромагнитного поля вместе с граничными условиями на границах раздела. Так как законы распространения тепла и электромагнитного поля описываются похожими уравнениями, то методы решения уравнения теплопроводности и уравнений Максвелла на этапе их дискретизации могут использовать одну и ту же сетку разбиения.
Погрешность в расчете SAR зависит от погрешности в расчете поля и от точности задания параметров модели головы. Важным фактором является также степень упрощения корпуса и конструкции телефона, от чего зависит точность расчета его ближнего поля. Очевидно, что подробное моделирование корпуса телефона является довольно сложной задачей, требующей больших вычисительных ресурсов. Однако, основные детали корпуса, особенно близкие к антенне и к голове человека, желательно описать как можно точнее.
Необходимо также учитывать, что напряжение источника, действующего на входе антенны, изменяется во времени. Поэтому для каждого временного шага необходимо рассчитывать значения электрического и магнитного полей во всех точках анализируемого пространства, в соответствии с методом расчета во временной области (метод FDTD) [1].
Голова человека как биологическое тело имеет достаточно сложную структуру и состоит из несколько видов биологической ткани с потерями, в которых мощность поглощается неравномерно из-за разницы в проводимости каждого слоя. Поэтому важным условием точности расчетов внутреннего поля является достаточно точное построение модели головы. Следует отметить, что электромагнитное моделирование тела человека в настоящее время, - это одна из выделившихся научных областей на стыке радиотехники и медицины. Усредненные значения основных параметров часто используемой трехслойной модели головы человека представлены в таблице 1.
Табл.1. Параметры трехслойной модели головы человека (в скобках указаны параметры для частоты 1.9 ГГц)
Вещество |
Толщина, мм |
Отн.диэлектрическая проницаемость |
Проводимость слоя См/м |
tg d [6] |
Плотность слоя кГ/м3 |
Теплопроводность Вт/К/м |
|
Мозг |
заполнение |
53 (0.9 ГГц) 46 (1.9 ГГц) |
1.1 (1.7) |
0.415 (0.369) |
1030 |
0.3 |
|
Кость |
3 |
9 ( 8 ) |
0.06 ( 0.1 ) |
0.133 (0.125) |
1800 |
0.01 |
|
Кожа |
1 |
59 (46) |
1.3 (1.9) |
0.44 (0.41) |
1100 |
0.5 |
В данной работе численное исследование взаимодействия электромагнитного излучения телефона и головы пользователя было проведено с помощью программы CST MicrowaveStudio (CST MS). Эта программа основана на методе FDTD и позволяет найти все требуемые характеристики поля в ближней и дальней зонах.
Влияние головы пользователя на работу сотового телефона
Сначала рассчитали частотную характеристику сотового телефона, т.е. зависимость согласования его антенны от частоты, в свободном пространстве в отсутствие головы пользователя. Модель телефона состоит из спиральной антенны и металлического корпуса (см. рис. 1а). Полученная зависимость представлена на рис.2 синим цветом. Можно предположить, что находящееся в ближнем поле антенны большое диэлектрическое тело (каковым является модель головы пользователя) должно влиять на характеристики антенны, в частности на ее согласование. Численные расчеты, проведенные с помощью CST MS в присутствии фантома головы человека, подтверждают это (см. рис. 2).
Рис. 2. Рассчитанные частотные характеристики антенны в отсутствие и в присутствии фантома головы пользователя. Параметр Udal - расстояние между телефоном и головой.
В программе CST MS расчет частотной характеристики выполняется как преобразование Фурье от временной зависимости волны, отраженной от антенны. Основное время в вычислениях занимает расчет во временной области. Программа может выполнять этот расчет сравнительно быстро на двухядерных компьютерах с частотой процессора 2 ГГц и RAM> 2 ГГб.
Рис. 3. Диаграмма направленности дальнего поля в системе телефон-голова для ряда значений параметра удаления телефона от головы.
На рис. 3 показаны рассчитанные диаграммы направленности в дальней зоне для нескольких значений параметра удаления телефона от головы. Изменение интенсивности излучения при разных значениях Udal (изменение коэффициента усиления антенны) объяснятся рассогласованием, см. рис. 2. Голова человека экранирует поле антенны в направлении ц=180 град. Поэтому поле излучения имеет значительно меньшую величину в этом направлении по сравнению с противоположным направлением «от головы».
Рис. 4. Распределение поля Е вдоль линии внутри головы.
Удельное поглощение рассчитывается на основе вычисленных значений электрического поля и заданных потерь в принятой 3-слойной модели головы, см. ф. (1). Интересно отметить, что рассчитанное распределение модуля электрического поля вдоль линии, проходящей внутри головы, имеет скачок на границе раздела сред, и, следовательно, в этой области можно ожидать наибольший нагрев.
Рис. 5. Различные опции при расчета величины SAR.
В диалоговых окнах, относящихся к вычислению величины SAR (рис. 5), можно выбрать метод вычисления, а именно: SAR в ячейке с массой усреднения 10 г или 1 г, SAR в точке и другие характеристики, а также стандарт (например, IEEE C95.3), в соответствии с которым рассчитывается эта величина. Средняя мощность излучения антенной системы, от которой зависит величины SAR, принимаем равной 200 мВт.
Рис. 6. Зависимость максимального значения SAR от расстояния между телефоном и головой.
Зависимость SAR от расстояния между телефоном и головой (рис. 6) имеет несколько неожиданный характер (уменьшение SAR при приближении телефона к голове от 2 мм до 0 мм) и, это может быть объяснено частичной экранировкой определенных областей головы близко расположенным корпусом телефона. Значение SAR=2.7 Вт/кГ достигается при расстоянии 2 мм от головы, а при значительном удалении снижается до 1.5 Вт/кГ, что является близким к допустимой норме значением [5].
Расчет распределения температуры внутри головы пользователя
По рассчитанным значениям внутреннего поля на следующем шаге можем решить задачу вычисления температурного распределения внутри головы. Для этого в программе CSTMWS необходимо задать соответствующие граничные условия.
a) распределение температуры в сечении со стороны телефона |
б) распределение температуры в сечении с обратной стороны головы |
Рис. 7. Рассчитанные распределения температуры в сечениях фантома головы пользователя сотового телефона.
Рис. 8. Рассчитанная температура вдоль линии максимального SAR для нескольких значений расстояния телефона от головы.
Из графиков на рис. 8 следует, что с обратной (по отношению к телефону) стороны головы температура на 1.3є меньше, чем со стороны телефона, но тем не менее она остается на 1.3є выше температуры окружающей среды. Если сравнить рис. 8 и 4, то видно, что величина SAR внутри головы уменьшается значительно быстрее, чем температура, которая даже на стороне, противоположной телефону, превышает температуру окружающей среды. Из этого результата заключаем, что значения и распределение SAR не полностью отражают влияние поля и поглощения на распределение температуры внутри головы. Это объясняется тем, что установившая температура зависит, в частности, и от теплопроводности биологической среды. Отсюда можно сделать вывод, что расчет температуры внутри головы является важным при оценке степени воздействия СВЧ излучения на организм человека.
Ранее [6, 7] с помощью программы CST MWS была смоделирована сложная конфигурация корпуса сотового телефона, с крышкой и традиционной спиральной антенной, в свободном пространстве. Прототипом является телефон, измеренная величина SAR у которого равна 1.5 Вт/кг [6]. Расчет SAR с помощью программы CST MS показал хорошее совпадение с экспериментом [6].
Отметим одно из перспективных приложений данной технологии моделирования электромагнитных и тепловых процессов в биологическом теле. Известно, что в ряде современных медицинских приборов в области гипертермии нагрев биологического тела используется для достижения положительного эффекта. Программа CST MS позволяет моделировать эти сложные тепловые эффекты нагрева и рассчитывать распределения напряженности электрического поля и плотности потерь мощности внутри тела с соответствующим этому поглощению распределением температуры.
Выводы
Показано, что теплопередача биологического материала играет важную роль при моделировании эффекта поглощения мощности в голове человека при пользовании сотовым телефоном. При задании установившегося значения излучаемой мощности 200 мВт рассчитанное значение нагрева составляет примерно 2є по сравнению с температурой окружающей среды.
Данный подход является перспективным для проектирования современных медицинских приборов, например, в гипертермии (медицинская терапия), которая использует сфокусированное СВЧ излучение для локального нагрева ткани внутри тела до определенной температуры. Компьютерное моделирование и медицинское исследование в комплексе позволяют перевести лечение СВЧ излучением на качественно новый более высокий уровень.
сотовый телефон температура электромагнитный
Литература
1. H.H. Pennes, Analysis of tissue and arterial blood temperatures in the resting forearm, J. Appl. Physiol., 1, 93-122 (1948).
2. Paolo Bernardi, etc. Specific Absorption Rate and Temperature Increases in the Head of a Cellular-Phone User. IEEE Trans. on MTT, Vol. 48, No.7, July 2000.
3. Hanus, M. Luong, F. Lethimonnier, "Electromagnetics Fields and SAR Computations in a Human Head with a Multi-port Driven RF Coil at 11.7 Tesla", Proc. Intl. Soc. Mag. Reson. Med. 13 (2005).
4. Biological Effects in the cm/mm Wave Range, Part II/III.
5. IEEE C95.1-1991, “IEEE Standard for Safety Levels with Respect to Human Exposure to Radio Frequency Electromagnetic Fields, 3 kHz to 300 GHz', New York: The Institute of Electrical and Electronics Engineers, Inc., April 1992.
6. Курушин А., Титов А. Расчет мощности, поглощаемой в голове пользователя сотового телефона//CHIP NEWS. 2001. 8. С.40 45.
7.Курушин А.А., Титов А.П. Расчет антенной системы сотового телефона в присутствии биологического тела. LVII научная сессия посвященная дню радио. 15-16 мая 2002 г., Москва. Том 2, стр.97-99.
8. Банков С.Е., Курушин А.А., Разевиг В.Д. Анализ и оптимизация трёхмерных СВЧ-структур с помощью HFSS М.: СОЛОН-Пресс, 2005. - 216 с.
9. Antonio Faraone, Q. Baszano, D. Simunic. Experimental Dosimetry in a Sphere of Simulated Brain Tissue Near a Half-Wave Dipole Antenna. IEEE Trans. On Veh. Tech., Vol. 44, pp. 390-403.
10. J. Wang and O. Fujiwara, "Reduction of electromagnetic absorption in the human head for portable telephones by a ferrite sheet attachment", IEICE Trans. Commun., vol. E80-B, no.12, pp.1810-1815, Dec. 1997.
11. Калиничев В.И., Курушин А.А. Микрополосковые антенны в сотовых телефонах. New Chips (Инженерная микроэлектроника), 2001, № 7, стр. 5-10.
12. M.Pelosi, O.Franek, M.B.Knudsen, G.F.Pedersen, and J.B.Andersen. Antenna Proximity Effects for Talk and Data Modes in Mobile Phones. IEEE Antennas and Propagation Magazine, Vol.52, No 3, June 2010, pp. 15-27.
Размещено на Allbest.ru
...Подобные документы
Оценка влияния атмосферной термической неоднородности на атомное поглощение электромагнитного излучения. Основные сведения о спектроскопии. Эффекты Зеемана и Штарка. Профиль атомного поглощения в условиях градиента температуры. Канал передачи данных.
дипломная работа [610,6 K], добавлен 21.04.2016Основные параметры электромагнитного поля и механизмы его воздействия на человека. Методы измерения параметров электромагнитного поля. Индукция магнитного поля. Разработка технических требований к прибору. Датчик напряженности электромагнитного поля.
курсовая работа [780,2 K], добавлен 15.12.2011Измерение поглощаемой мощности как наиболее распространенный вид измерения СВЧ мощности. Приемные преобразователи ваттметров проходящей мощности. Обзор основных методов для измерения импульсной мощности, характеристика их преимуществ и недостатков.
реферат [814,2 K], добавлен 10.12.2013Методика и характеристика основных этапов расчёта ребристого радиатора при естественном воздушном охлаждении для транзистора 2Т808А заданной мощности 15 Вт. Определение необходимого напора внутри радиатора, температуры среды и коэффициента теплоотдачи.
реферат [12,0 K], добавлен 10.11.2009Принцип действия трансформатора, элементы его конструкции. Вычисление мощности фазы, номинальных токов и короткого замыкания. Расчет основных размеров трансформатора и обмотки. Определение размеров магнитной системы, массы стали и перепадов температуры.
курсовая работа [649,9 K], добавлен 25.06.2011Расчет и определение режимов работы двигателя. Выбор мощности двигателя для продолжительного режима работы с повторно-кратковременной нагрузкой, проверка на перегрузочную способность, пусковые условия. Вычисление потребляемой мощности, расшифровка марки.
контрольная работа [248,7 K], добавлен 07.02.2016Установление методами численного моделирования зависимости температуры в точке контакта от угла метания пластины при сварке взрывом. Получение мелкозернистой структуры и расчет параметров пластины с применением программного расчетного комплекса AUTODYN.
дипломная работа [6,2 M], добавлен 17.03.2014Определение мощности лазерного излучения, подаваемого на образец. Вычисление размеров лазерного пучка на образце. Разработка системы измерения мощности излучения и длительности лазерного импульса, системы измерения температуры в зависимости от времени.
лабораторная работа [503,2 K], добавлен 11.07.2015Характеристики потребителей электроэнергии. Расчет электрических нагрузок. Определение мощности компенсирующего устройства реактивной мощности. Выбор числа и мощности трансформаторов подстанции. Вычисление параметров и избрание распределительной сети.
курсовая работа [884,2 K], добавлен 19.04.2021Определение линейного теплового потока методом последовательных приближений. Определение температуры стенки со стороны воды и температуры между слоями. График изменения температуры при теплопередаче. Число Рейнольдса и Нусельта для газов и воды.
контрольная работа [397,9 K], добавлен 18.03.2013Определение мощности теплового потока при конвективной теплопередаче через трубу заданного диаметра. Расход пара на обогрев воды в пароводяном теплообменнике, превращение пара в конденсат. Изменение температуры теплоносителей вдоль поверхности нагрева.
контрольная работа [308,7 K], добавлен 13.05.2015Выбор метода измерения СВЧ в зависимости от необходимой точности, диапазона находимой мощности, сложности конструкции измерителя и стоимости прибора. Диагностирование поглощаемой и проходящей силы источника электромагнитных колебаний (генератора).
доклад [233,0 K], добавлен 01.02.2011Определение физических величин, явлений. Изменение температуры углекислого газа при протекании через малопроницаемую перегородку при начальных значениях давления и температуры. Сущность эффекта Джоуля-Томсона. Нахождение коэффициентов Ван-дер-Ваальса.
контрольная работа [231,7 K], добавлен 14.10.2014Анализ модели температуры в радиально бесконечном пласте. Моделирование давления и температуры сигнала, связанного с переменной скоростью. Определение сигнала температуры отдельного слоя связанного с постоянной скоростью добычи слабо сжимаемой жидкости.
курсовая работа [770,7 K], добавлен 20.02.2021Основные шкалы измерения температуры. Максимальное и минимальное значение в условиях Земли. Температура среды обитания человека. Температурный фактор на территории Земли. Распределение температуры в различных областях тела в условиях холода и тепла.
доклад [1,0 M], добавлен 18.03.2014Анализ квантовой теории полей. Способ получения уравнения Клейна-Гордона-Фока для электромагнитного поля и его классическое решение, учитывающее соответствующие особенности. Процедура квантования (переход к частичной интерпретации электромагнитного поля).
доклад [318,7 K], добавлен 06.12.2012Расчет основных величин трансформатора станции. Определение потерь короткого замыкания, механических сил в обмотках и их нагрева. Вычисление размеров магнитной системы и потерь холостого хода трансформатора. Расчет превышения температуры устройствами.
курсовая работа [2,5 M], добавлен 22.06.2015Способы регулирования температуры воды в электрических водонагревателях. Методы интенсификации тепломассообмена. Расчет проточной части котла, максимальной мощности теплоотдачи конвектора. Разработка экономичного режима работы электродного котла в Matlab.
магистерская работа [2,5 M], добавлен 20.03.2017Расчет мощности и выбор соответствующего оборудования для трансформаторной электрической подстанции двух предприятий - потребителей энергии первой и третьей категории. Определение мощности и числа трансформаторов, расчет токов короткого замыкания.
курсовая работа [413,2 K], добавлен 18.05.2012Определение объема магнитопровода, оптимальной магнитной индукции, потерей мощности, плотности тока в проводах обмоток, добавочных потерь. Выбор сечений проводов. Расчет тепловых режимов, схемы замещения трансформатора. Его моделирование в среде OrCAD.
курсовая работа [696,4 K], добавлен 05.12.2012