Русский
Русский
English
Статистика
Реклама

Температура

Как получить температуру в -50C (и ниже!) на дому или вихревая трубка Ранка-Хилша под микроскопом

28.05.2021 12:18:53 | Автор: admin

А что, если я скажу Вам, что можно подуть в Т-образную трубку и с двух других её концов пойдет воздух сильно минусовой и сильно плюсовой температуры? Похоже на какой-то бред, не так ли? Тем не менее, такое замечательное устройство вполне себе существует и известно очень давно. Ученые до сих пор расходятся во мнениях относительно того, как же оно всё-таки работает?!. Предлагаем и Вам ознакомиться с этим любопытным эффектом

В 1931 г. Жозефом Ранком был открыт вихревой эффект энергетического разделения газов, называемый часто эффектом Ранка.

После доклада Ранка Французскому физическому обществу о своем открытии о нем забыли и только с 1946 года вихревой эффект стал объектом исследований ученых разных стран.

В Советском Союзе первые широкие исследования вихревого эффекта были проведены профессором Мартыновским В. С. и доцентом Алексеевым В. П. в Одесском технологическом институте пищевой и холодильной промышленности.

Большую работу по созданию вихревых вакуум-аппаратов провела группа под руководством д-ра техн. наук Дубинского М. Г.; им опубликовано несколько теоретических работ по закрученным потокам.

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

Несмотря на длительное изучение вихревого эффекта до сих пор отсутствует общепризнанная физико-математическая модель феномена.

Наиболее обоснованной считается гипотеза взаимодействия вихрей, но и она не лишена некоторых неточностей и противоречий. Пока не проведено строго обоснованного прямого эксперимента, способного полностью подтвердить ее достоверность.

На основе проведенных исследований разработаны полуэмпирические методики расчета вихревого эффекта и некоторых видов вихревых аппаратов.

На этом период первоначального изучения вихревого эффекта закончился и начался период освоения и внедрения его в производство.

С 1953 года в Куйбышевском авиационном институте велась работа по исследованию вихревого эффекта и его промышленному использованию. В 1958 году в институте создана лаборатория промышленного применения вихревого эффекта, в которой разработано несколько промышленных образцов вихревых холодильно-нагревательных установок и созданы различные типы вихревых аппаратов.

Благодаря своим особенностям вихревой эффект нашел практическое применение в самых различных областях техники и производства.

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

На многих предприятиях СССР работали вихревые холодильные камеры, вихревые термостаты и вихревые вакуум-насосы.

В данный момент, вихревые трубы широко используются в разных отраслях промышленности по всему миру:

image


image


image


image


Сущность вихревого эффекта и конструкции вихревых труб


Вихревой эффект, или эффект Ранка, проявляется в закрученном потоке вязкой сжимаемой жидкости и реализуется в очень простом устройстве, называемом вихревой трубой (трубой РанкаХилша, вихревым энергоразделителем, вихревым холодильником), схематичная конструкция которой изображена на рис. 1. и рис. 2.

То есть, можно сказать, что Вихревой эффект (эффект Ранка Хилша, англ. Ranque-Hilsch Effect) эффект температурного разделения газа при закручивании в цилиндрической или конической камере при условии, что поток газа в трубке проходит не только прямо, но и обратно.

image
Рисунок 1.

image
Рисунок 2.

Вихревая труба представляет собой гладкую цилиндрическую трубу 1, снабженную тангенциальным соплом 2, улиткой 3, диафрагмой 4 с осевым отверстием и дросселем 5.

При втекании газа через сопло образуется интенсивный круговой поток, приосевые слои которого заметно охлаждаются и отводятся через отверстие диафрагмы в виде холодного потока, а периферийные слои подогреваются и вытекают через дроссель в виде горячего потока.

По мере прикрытия дросселя общий уровень давления в вихревой трубе повышается и расход холодного потока через отверстие диафрагмы увеличивается при соответствующем уменьшении расхода горячего потока. При этом температуры холодного и горячего потоков также изменяются.

На периферии образуется закрученный поток с большей температурой, а из центра в противоположную сторону выходит охлажденный поток. Существует распространённое заблуждение, что температурное разделение происходит путём перемещения молекул газа на прямом проходе вихря (в одну сторону). Но объяснимых физикой причин для такого разделения нет, как нет причин и для вращения центрального жгута в противоположную сторону относительно периферии. В противоположную сторону вращаются микровихри между центральным жгутом и периферией, так как жгут вращается с более высокой скоростью относительно периферии. Но катятся они, как ролики в подшипнике, в ту-же сторону, в которую вращаются внешний слой и центральный жгут. Температурное разделение происходит путём теплопередачи от сжатого (и потому горячего) кумулятивным эффектом или имплозией центрального жгута к несжатой периферии, имеющей температуру как на входе. По мере движения к горячему концу периферия нагревается от двигающегося ей навстречу сжатого горячего центрального жгута, который в свою очередь наоборот остывает. Т. о. образующийся в трубке вихрь является тепловым насосом компрессионного типа с противоточным теплообменником, способным передать до 100 % разницы температур. Поэтому для терморазделения необходим не только прямой, но и обратный проход, как на рисунке. Так как после выхода из трубки жгут расширяется до давления окружающей среды (атмосферного), выходящий из холодного конца трубки газ имеет температуру намного ниже температуры окружающей среды (если горячий конец не заглушен), а всё утерянное им тепло уносится газом с горячего конца.

- Достоинства вихревой трубки:

  • Нет движущихся частей;
  • Не требуется обслуживание;
  • Надежность;
  • Без электричества или химикатов;
  • Компактная, легкая;
  • Бюджетная;
  • Мгновенное получение результата;
  • Прочность конструкции;
  • Регулируемая температура;


Парадоксальность эффекта Ранка заключается в том, что центробежные силы во вращающемся потоке направлены наружу. Как известно, более тёплые слои газа или жидкости имеют меньшую плотность и должны подниматься вверх, а в случае центробежных сил стремиться к центру, более холодные имеют большую плотность и, соответственно, должны стремиться к периферии. Между тем, при большой скорости вращающегося потока всё происходит с точностью до наоборот!

Эффект Ранка проявляется как для потока газа, так и для потока жидкости, которая, как известно, является практически несжимаемой и потому фактор адиабатического сжатия / расширения к ней неприменим. Тем не менее, в случае жидкости эффект Ранка обычно выражен значительно слабее возможно, именно по этой причине, да и очень малая длина свободного пробега частиц затрудняет его проявление. Но это верно, если оставаться в рамках традиционной молекулярно-кинетической теории, а у эффекта могут быть и совсем другие причины.

В результате многочисленных экспериментальных исследований создано несколько конструктивных вариантов вихревых труб.

Основным их различием является конструктивное выполнение тангенциального соплового входа сжатого газа и длина цилиндрической части (вихревой зоны) трубы в калибрах.
На рисунке 3 приведены конструктивные особенности исследованных вихревых труб, показавших хорошие результаты.

image
Рисунок 3.

Разработанные Р. Хилшем вихревые трубы имеют один круглый утопленный тангенциальный сопловой вход и входную улитку прямоугольного сечения, ввиду чего на срезе сопла имеется площадка, создающая зону завихрения.

Этот недостаток устранен в конструкции В. С. Мартыновского и В. П. Алексеева, которые создали тангенциально-лотковый сопловой вход, имеющий два круглых наполовину утопленных в теле трубы сопловых входа, переходящих в лоток.

Подобная конструкция позволила устранить зону завихрения и улучшить эффективность вихревой трубы, хотя и усложнила изготовление соплового входа.
Оптимальной длиной вихревой зоны трубы (расстояние от соплового входа до дросселя горячего конца) авторы двух рассмотренных выше конструкций считают 50 калибров.

Конструкция Меркулова А. П. имеет сопловой вход прямоугольного сечения и прямоугольную входную улитку, построенную по спирали Архимеда, что обеспечивает устранение зоны завихрения на сопловом срезе при сохранении простоты конструкции.

Второй особенностью этой конструкции является сокращение длины вихревой зоны до 9 калибров, осуществляемое за счет ограничения вихревой зоны спрямляющей четырехлопастной крестовиной, устанавливаемой перед дросселем горячего потока (рис 4. и рис 5. )

Введение этих двух элементов позволило еще более повысить эффективность вихревой трубы и сделать ее значительно компактнее.

image
Рисунок 4.

image
Рисунок 5.

Влияние диаметра отверстия диафрагмы


С полным основанием можно ожидать, что на характеристики вихревых труб различных диаметров будут влиять не абсолютные значения диаметров отверстия диафрагмы, а их отношение к диаметру вихревой трубы:

image


Экспериментальные данные трех различных исследований хорошо укладываются на прямую линию (График А).
С изменением соотношения, характеристики вихревой трубы заметно меняются (График Б, где -относительный весовой расход холодного потока).

image


Влияние длины вихревой зоны трубы


В ряде работ за оптимальную длину L вихревой зоны трубы принималась длина, соответствующая 50 ее калибрам.

Это делало вихревую трубу очень громоздкой, поэтому многие исследователи пытались уменьшить ее.

В результате было найдено конструктивное решение, обеспечивающее резкое сокращение вихревой зоны.

За счет установки на горячем конце вихревой зоны четырехлопастной крестовины удалось сократить длину этой зоны до 9 калибров и при этом несколько улучшить эффективность вихревой трубы как по эффекту охлаждения, так и по холодопроизводительности.

На рис. 6 приведены опытные кривые зависимости максимальных значений эффекта охлаждения от относительной длины вихревой зоны для двух отличных друг от друга давлений и проходных сечений соплового входа.

При сокращении длины вихревой зоны до оптимальной наблюдается небольшой рост максимального значения эффекта охлаждения, а при (L/D)<9 последний резко снижается.

Выявлено, что и максимальная холодопроизводительность изменяется аналогично.

image
Рисунок 6.

Следующим шагом по пути сокращения вихревой зоны явилось запатентованное В. Хендалом предложение делать вихревую зону в виде усеченного конуса с расширением в сторону дросселя и с углом конусности около 7.

Это позволило сократить длину вихревой зоны до 4-х калибров при сохранении максимального эффекта охлаждения вихревой трубы.

Влияние площади проходного сечения сопла


Экспериментами установлено, что для получения оптимальных характеристик имеется определенное соотношение между диаметром вихревой трубы и площадью Fc проходного сечения соплового входа.

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

Эти два противоположно влияющих фактора дают оптимальное значение величины проходного сечения сопла, которое удобнее выразить относительной величиной, являющейся отношением площади Fc проходного сечения сопла к площади поперечного сечения вихревой трубы:

image


Согласно указанным выше экспериментам оптимальное значение относительной площади проходного сечения сопла лежит в пределах:

image


При повышении давления питающего сжатого газа оптимальные значения приближаются к нижнему пределу.

Для случая получения максимального эффекта охлаждения при малых значение Fc увеличивается.

Для большинства случаев расчета вихревой трубы принимается:

image


Сопловой ввод конструктивно должен быть выполнен так, чтобы при сохранении тангенциальности он обеспечивал плавный вход сжатого газа в вихревую зону и осевую симметрию формирующегося вихря.

Отсутствие зон завихрения и плавность входа струи обеспечивает конструкция 3 (см. рис. 3.), имеющая прямоугольное сечение сопла и входную улитку.

Хорошую осевую симметрию и плавность входа обеспечивают конструкции 2 и 6, но они труднее в изготовлении и требуют усложненной системы подвода сжатого газа, поэтому чаще используется конструкция 3 с прямоугольным сечением сопла.

Влияние давления сжатого газа


При условии сохранения постоянным давления холодного потока за отверстием диафрагмы величина давления, подаваемого к соплу сжатого газа, играет существенную роль в получаемом эффекте охлаждения.

С повышением давления входа растут эффект охлаждения и холодопроизводительность.
На рис. 7 изображены кривые эффектов охлаждения tx и холодопроизводительности tx для различных давлений p.

image
Рисунок 7.

Влияние температуры сжатого газа


Результаты экспериментов показали, что с ростом температуры, растет эффект охлаждения холодного потока и эффект подогрева горячего потока, значение же температурной эффективности при этом практически остается неизменным.

Влияние масштаба вихревой трубы


С уменьшением масштаба вихревой трубы возрастает влияние эффекта подмешивания пограничного слоя диафрагмы к холодному потоку, что приводит к ухудшению получаемого эффекта охлаждения последнего.

Влияние масштабного фактора впервые экспериментально было обнаружено Хилшем, а затем подтверждено Мартыновским и Алексеевым. Автором были исследованы вихревые трубы с диаметрами D= 10, 16, 20, 25, 30 и 33 мм.

В этих исследованиях наблюдался рост эффективности вихревой трубы с ее масштабом.
Несмотря на то, что исследованные различными авторами серии вихревых труб имели различный диапазон диаметров и различное конструктивное исполнение при несколько отличающихся геометрических соотношениях, полученные максимальные эффекты охлаждения имеют довольно четкую зависимость от диаметра D вихревой трубы.

Эта зависимость для двух давлений входа представлена на рис. 8. Как видно из рисунка, экспериментальные точки могут быть аппроксимированы в прямые линии, несколько расходящиеся с ростом D.

image
Рисунок 8.

Возможность местного получения холода является важным достоинством вихревых труб, масштаб которых может быть как угодно малым.

Например, вихревая труба диаметром 1 мм и длиной 10 мм (близкая к размерам полупроводниковых элементов), расходуя около 15 л воздуха в минуту может обеспечить отвод от полупроводникового элемента рассеиваемой мощности порядка 5 вт.

Это делает ее незаменимым устройством для стабилизации температуры отдельных ответственных элементов полупроводниковых электронных устройств летательных аппаратов.

Теперь, возвращаясь к декларированной в названии статьи возможности получения температуры в -50 градусов, можно сказать следующее. Опираясь на результаты эксперимента, проведённого в видео, приложенном ниже, даже без каких-либо расчётов, можно грубо прикинуть минимально возможную температуру, которую мы можем получить. Предположим, что в качестве компрессора, будет использован компрессор от холодильника. Основным плюсом которого является поршень достаточно малого диаметра (около 10 мм). Однако у данного поршня есть свои плюсы: благодаря своему малому диаметру, он может качать воздух под достаточно большими давлениями.


Автором данной статьи экспериментально было выявлено, что компрессор холодильника может качать воздух под давлением в 20 атмосфер. Так как в эксперименте, показанном в ролике выше была озвучена цифра в 4 атмосферы, необходимая для получения температуры в -10 градусов, возьмём эти цифры в качестве ориентира. Основываясь на них, по формуле пропорции, высчитаем необходимую величину давления, которая позволит нам получить нужную отрицательную температуру (конечно, зависимость между давлением и температурой,- скорее всего, не будет настолько линейной, но мы можем взять за некий базис эти числа, чтобы было от чего оттолкнуться):

4 = -10
20 = х, отсюда х = -200/4 = -50 градусов.

Таким образом, используя широкодоступные компрессора от бытовых холодильников (которые можно купить в магазинах запчастей для холодильников), мы можем достичь температуры в -50 градусов.

В инструкции ниже, вкратце рассмотрено, как переделать компрессор от бытового холодильника в источник высокого давления (20 атмосфер) для экспериментов с вихревыми трубками.

Нам для этого понадобятся следующие компоненты:

image


Далее, мы делаем следующее:

image

В итоге должно получиться нечто такое:

image


Некоторые дополнительные пояснения к процессу переделки компрессора от холодильника:

1. Не стоит сильно беспокоиться за компрессор, что он должен работать постоянно в смазке, без нее-он выйдет из строя. У автора данный компрессор от холодильника работал для целей аэрографии в течение 3 лет. Только периодически (раз в месяц) доливал во входной патрубок 3-4 см3 машинного масла. Но, если хотите сделать всё более аккуратно, можно на входной патрубок установить лубрикатор специальное устройство, насыщающее всасываемый воздух распыленным маслом. Используется для смазки компрессора. Тогда вся система в сборе будет выглядеть так: лубрикатор-->компрессор-->ресивер(сойдет даже бутылка 2л от кока-колы это по моему личному, многолетнему опыту). Ресивер нужен для накопления давления, чтобы избежать пульсаций воздух и отстоя паров масла и воды в подаваемом воздухе-->вихревая трубка;

2. На напечатанной с помощью 3d печати вихревой трубке следует нарезать резьбу М10х1(если вы будете использовать рекомендованные мной компоненты). Так как печатная такая мелкая резьба получается оплывшей и лучше ее нарезать плашкой;

3. Ниже вы найдете ссылки на требующиеся в процессе переделки компоненты:



4. Во многом, вихревые трубки изготовляются, основываясь на эмпирических данных, предыдущих исследователей (примерно такого размера и таких параметров вроде работало...). Отсутствие общепризнанной теории вихревых труб сильно осложняет процесс их изготовления согласно заданных параметров (о чём уже было сказано выше).

Однако, не всё так плохо, если:

а) компрессор справляется с подачей нужного количества воздуха в трубку;

б) трубка изготовлена, опираясь на имеющиеся эмпирические данные (например, рисунок 3.) такие трубки работают сразу, почти все. Они отличаются только своими возможностями получения нужных низких/высоких температур.

Самые экстремальные исследователи, могут попробовать запитать вихревую трубку воздухом высокого давления (до 300 атмосфер!) от появившихся в изобилии в последние годы компрессоров высокого давления для PCP пневматических винтовок.
Полагаем, результат будет весьма впечатляющим.

Автор данной статьи, в процессе её подготовки, заметил еще один интересный факт, который должен быть озвучен: все публично доступные результаты экспериментов с этими трубами, основываются на стандартных технологиях.

Под этим подразумевается, что для производства данных вихревых труб используется стандартные металлические детали, изготовленные с применением стандартных металлообрабатывающих станков.

Налицо полное игнорирование возможностей изготовления данных трубок с применением 3D печати!

Однако именно 3D печать может дать новые возможности в изготовлении данных устройств!
Например: изготовление микромассива данных трубок, объединённых в миниатюрный кластер.

Данный кластер может быть установлен в дальнейшем прямо на критичные электронные компоненты, требующие охлаждения, например процессор компьютера.

А с появлением общедоступных фотополимерных принтеров по недорогой цене, например, Anycubic Photon,- возможно изготовление поистине миниатюрных и сложных кластерных систем, базирующихся на вихревых трубках (в буквальном смысле микронного масштаба)!

Таким образом, любой исследователь данного эффекта, имеющий в наличии соответствующей трехмерный принтер и программу инженерного трехмерного проектирования, например, solidworks, может с лёгкостью смоделировать и распечатать вихревую трубку или их кластер любого размера, соответствующего возможностям компрессора, имеющегося в наличии у конкретного исследователя!

Производители на aliexpress тоже не остались в стороне и выпускают весьма интересное устройство Мобильный кондиционер сварщика, основанный на данной трубке Ранка-Хилша. Данное устройство позволяет обеспечить комфортным микроклиматом рабочих в цеху персонально каждого. Это устройство достаточно подключить к цеховой сети сжатого воздуха:

image

При желании провести ряд собственных экспериментов, там же имеются и готовые вихревые трубки.

Используя 3d печать, любой исследователь может легко и быстро создать целую серию вихревых трубок, подогнав их размеры и параметры получаемой температуры под требующиеся именно ему!

С 1960-х годов вихревое движение является темой множества научных исследований. Регулярно проводятся специализированные конференции по вихревому эффекту, например, в Самарском аэрокосмическом университете.

Следует заметить, что возможности некоторых установок на эффекте Ранка внушают уважение например, рекордное охлаждение, которого удалось достигнуть на одной ступени, составляет более 200С!

Желающие более подробно ознакомиться с данным эффектом и вихревыми трубами могут прочитать следующую литературу:

  • А.П. Меркулов Вихревой эффект и его применение в технике;
  • Ш.А. Пиралишвили, В.М.Поляев, М.Н. Сергеев Вихревой эффект. Эксперимент, теория, технические решения.



Подробнее..

О том, как мы температуру в ЦОД мерили

21.04.2021 14:19:32 | Автор: admin

Если у вас большой и серьезный ЦОД, то параметрия температурных режимов не является проблемой. Существуют проверенные решения, например, программируемые контроллеры TAC Xenta, которые работают через LonWorks. Именно так мы собираем данные в московском ЦОД Datahouse. Но непосвящённому смертному весьма непросто собрать правильные показатели из этой связки и выводить их в мониторинг в нужном виде. К тому же решение промышленное и достаточно дорогостоящее. Поэтому при строительстве новой гермозоны
в Екатеринбурге мы решили поэкспериментировать и внедрить альтернативное решение по измерению температуры в холодных и горячих коридорах.

Ничто не предвещало беды

Так как множество систем в этом ЦОДе завязано на открытом коммуникационном протоколе Modbus мы решили заказать температурные датчики, работающие поэтой шине, и собирать данные с дальнейшей интерпретацией в интерфейсе мониторинга. Недорогие датчики быстро нашлись на известном китайском сайте и были заказаны партиями в количестве 20 и 40 штук.

Первая партия из 20 штук пришла достаточно оперативно, но при детальном рассмотрении стало понятно, что датчики незначительно отличаются корпусами. Важно ли это, как выяснилось да.

Из первой партии завелось 15 датчиков. Так как острой потребности в остальных не было, пока работали с ними. К моменту прихода второй партии выявили, что часть уже установленных в шину датчиков имеет поведение новогодней елки: показывают некорректные данные, отдают checksum error или отваливаются по таймауту.

Анализ второй партии показал, что в ней процент брака в разы выше.
В итоге из 60 датчиков мы получили всего 8 стабильных с правильными показаниями.

Вот так это выглядело:

Самое простое решение отказаться от этой затеи или заказать кучу датчиков и выбрать
из них не бракованные, но это не наш метод.

Не ищем легких путей

Вскрыли, посмотрели: микросхемы идентичны. Значит дело в прошивке.

Кстати, обратите внимание на заводское крепление кабеля оно было лишь на половине устройств. На остальных пришлось заливать кабель из термопистолета получилось прочней
и надежней.

Пока разбирались с работоспособностью датчиков поняли, как без особых сложностей обнаружить кривую версию прошивки. Если кроме Modbus работают обычные текстовые команды READ, PARAM, AUTO, STOP значит прошивка хорошая. В мертвых прошивках текст не отдается.

Решили взять прошивку с этих живых 8 датчиков, купили программатор Nu-Link,
но хитрые китайцы заблокировали чтение прошивки. То есть перезалить можно, а считать - нельзя. Запросы правильных прошивок у поставщика потерпели фиаско:
Я продаван, я не разработчик.

И тут я психанул, схватил крепкие напитки и закрылся у себя.
Через пару дней вышел с прошивками и программой.

За основу был взят Keil, пакет С51, позволяющий работать с 8-битными MCU.

В начале я научил читать сенсор SHT 20 (который собственно и снимает температурные данные), потом научил передавать эти данные по Modbus. В виду того, что этот MCU ни что иное как Nuvoton N76E003AT20, то вся база знаний, видимо, сосредоточена в руках наши китайских друзей.

В итоге i2c и Modbus сделал быстро, а вот с таймерами пришлось повозиться. Чтобы не было коллизий в шине, добавил возможность смены SLAVE_ID без выключения датчика в Китайской версии прошивки после смены адреса его нужно было обязательно выключать, что не очень удобно.

Начали шить новой прошивкой, из семи прошитых получили семь успешных и стабильных устройств. Так понемногу реанимировали все датчики, присвоили им номера и подключили витой парой в шину.

Так стало:

Если говорить про результаты измерений, то нормальные показания появляются только
в помещении с явной циркуляцией и движением воздуха. Без продува датчик показывает 30С. Это объясняется тем, что внутри установлен стабилизатор напряжения, который преобразует 24В в 3.3В, переводя разность в тепло.

Относительно влажности, пока можно утверждать, что требуется дополнительная калибровка, но способа ее произвести мы пока не придумали. Теме не менее, прямое воздействие, например, дыханием увеличивает показания влажности. В любом случае, через пару регистров это можно сделать по месту установки.

Несмотря на возникшие осложнения, данное решение имеет два очевидных преимущества: стоимость и гибкость. Датчики можно установить в любом удобном месте, точечно или объединить в гирлянды. Можно измерять как общую температуру, так и частные показания отдельных приборов и устройств. И самое важное, что все это отлично работает по Modbus.

Программа выложена на GitHub кому интересно, можно забрать и поиграться.
Стоимость датчика всего 300 , правда, нужен программатор.

Подробнее..

Умная одежда устройство модуляции температуры на основе графена

24.06.2020 10:13:26 | Автор: admin


У природы нет плохой погоды, как поется в знаменитой песне из кинофильма Служебный роман. Однако далеко не все готовы согласиться с этим утверждением. Кому-то нравится холод, кто-то предпочитает жару, кому-то все равно. Я же отношусь к тем людям, которые будут жаловаться и на жару, и на холод, нам подавай комфортные +20 C. К сожалению, не всегда и не у всех любителей нейтрального климата есть возможность жить в регионах, где он есть. Сейчас лето в самом разгаре, удушающая жара лишь изредка прерывается кратковременными грозами, которые не особо помогают. Если природа не готова идти нам навстречу, значит стоит делать что-то самим. Сегодня мы познакомимся с исследованием, в котором ученые из Манчестерского университета (Великобритания) разработали умную адаптивную ткань, способную снижать температуру тела человека ее носящего в жаркие дни. Что легло в основу умной ткани, как протекал процесс разработки, и какие дополнительные свойства и варианты применения имеются у этого изобретения? Ответы на эти вопросы мы найдем в докладе ученых. Поехали.

Основа исследования


Прежде, чем рассказать нам о своем творении, ученые отмечают, что прогресс в области пользовательской электроники за последние годы идет семимильными шагами. Буквально каждый день появляется что-то новенькое и необычное. Тем не менее, существует ряд ограничений, которые мешают тем или иным разработкам перейти от стадии лабораторных тестов к стадии массового производства. В аспекте производства умных тканей основной проблемой является сложность интеграции электронных / оптических материалов внутрь волокон ткани. Самый простой вариант в носимой электронике это создание отдельных гаджетов (браслеты, часы и т.д.), которые не требуют внедрения в другую систему (в данном случае, ткань), но спектр возможностей этих устройств будет ограничен.

По словам ученых, чтобы достичь вразумительных результатов в сопряжении электроники и ткани, необходимо либо изменить технологию производства ткани, либо использовать нестандартные материалы для электронной части носимого устройства.

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

Тепловое излучение от многослойного графена может модулироваться электрически через интеркалирование* ионов.
Интеркаляция* обратимое внедрение молекулы или группы молекул между другими молекулами или группами молекул.
В данном исследовании ученые представляют нашему вниманию технологию оптического текстиля, основанную на интеграции в текстиль динамических инфракрасных устройств на базе электрически перестраиваемого графена, образованного методом химического осаждения из паровой фазы (ХОПФ).

Результаты исследования


Устройства состоят из объединенных слоев инфракрасно прозрачного полимерного слоя, многослойного графена, выращенного с использованием метода ХОПФ, слоя тканевого разделителя и проводящей ткани (схема устройства на ).


Изображение 1

Изготовление начинается с выращивания многослойных графеновых пленок на никелевой фольге. Тонкая полиэфирная (PE) пленка, которая функционирует как прозрачный для инфракрасного излучения защитный слой, ламинируется на многослойный графен перед травлением Ni-фольги. Графен на полиэфирном листе прикрепляется к ткани с помощью термоплавкого клея.

Одним из важных моментов данной разработки является удобство использования и практичность, потому необходимо было удостовериться в хорошей адгезии между графеном и подложкой (тканью). Это было сделано посредством нескольких циклов стирки и посредством испытаний на механическое сжатие.

Далее на задний электрод (проводящая ткань) был нанесен ионный жидкий электролит (BMIMPF6), который впоследствии диффундировал в текстильную подложку. Текстиль действует как разделитель и ионопроводящий слой, обеспечивая ионное движение, когда разность напряжений приложена к графену и заднему электроду.

На 1b показаны примеры изготовленных устройств на натуральных (хлопок) и синтетических текстильных материалах (полиэфир).

Электрохимическая стабильность заднего электрода играет решающую роль в долговременной стабильности устройства. В качестве основы для заднего электрода тестировались разные материалы: проводящий текстиль на основе серебра, сетка из нержавеющей стали, золотое напыление, графен и восстановленный оксид графена.

Массив задних электродов и проводку на текстиле изготовили с помощью фотолитографии с последующей металлизацией и процессом отрыва*.
Отрыв* в технологии микроструктурирования представляет собой способ создания структур целевого материала на поверхности подложки с использованием жертвенного материала (например, фоторезиста).
Полученные пиксельные электроды позволяют определять динамические инфракрасные структуры на непрерывном графеновом слое с помощью выборочной интеркаляции.

Принцип работы устройств основан на обратимой интеркаляции ионов в графеновые слои и модулировании его электрических и оптических свойств. При 0 В многослойный графен имеет высокое инфракрасное поглощение, что приводит к высокой излучательной способности, раскрывая фактическую температуру устройства ().


Изображение 2

При подаче достаточной разности напряжений (> 2.5 В) ионная жидкость интеркалирует в слои графена, увеличивая оптическую проводимость и подавляя излучательную способность, тем самым скрывая фактическую температуру устройства. Термографы устройства записывались с помощью длинноволновой инфракрасной камеры, которая визуализирует изображения по закону Стефана-Больцмана:
P = T4
где P количество падающего теплового излучения на матрице болометров*; излучательная способность поверхности; постоянная Стефана-Больцмана; T температура поверхности в Кельвинах.
Болометр* тепловой приемник излучения (преобразует энергию поглощенного электромагнитного излучения в тепловую).
Текстильные устройства находясь непосредственно в тепловом контакте с источниками тепла, такими как тело человека, для предотвращения ложного экранирования температуры источника. Кроме того, графен функционирует как слой с высокой теплопроводностью, который удваивает температуропроводность в плоскости текстиля, улучшая теплопроводность от источника к поверхности.

Временной отклик устройств был получен путем записи видео тепловизором, чтобы получить изменение видимой температуры поверхности (2b).


Динамическое изменение инфракрасного излучения на хлопковом устройстве.

Полная интеркаляция (подавление излучательной способности) занимает ~5 с, когда ток устройства не ограничен. Стоит отметить, что эти измерения проводились в лабораторных условиях (21 C), ограничивающих минимальную кажущуюся температуру.

Устройства могут многократно циклически переключаться между состояниями с высокой и низкой излучательной способностью (2c), однако превышение электрохимического окна электролита ухудшает рабочие характеристики устройства.

Модуляция излучательной способности определялась количественно с помощью измерений отражения в инфракрасном и ближнем инфракрасном диапазонах с использованием инфракрасного Фурье-спектрометр (FTIR), оборудованного интегрирующей сферой. При 0 В коэффициент отражения внутреннего устройства почти плоский (2d) и составляет около 30%, за исключением поглощения в верхней полиэфирной пленке на длинах волн ~3.4, ~6.8, ~13.9 мкм и поглощения в атмосфере (например, CO2, H2O).


Демонстрация работы адаптивного инфракрасного текстильного устройства.

В диапазоне спектральной чувствительности тепловой камеры (8-13 мкм) такие поглощения минимизируются благодаря тщательному выбору верхней защитной пленки. Коэффициент излучения (или коэффициент поглощения) рассчитывается как 1 R, где R коэффициент отражения, поскольку свет не проходит через устройство. По мере того, как ионы интеркалируют графеновые слои, энергия Ферми и оптическая проводимость графена увеличиваются, тем самым увеличивая коэффициент отражения инфракрасного излучения.

Средняя излучательная способность устройства в диапазоне длин волн 8-13 мкм достаточно высока ( 0.7) для 0 В и поддерживается в таком значении до порогового напряжения ( 2.5 В) с последующим резким падением до 0.35 при > 4 В (), что отлично согласуется с термограммами на 2а.

Модуляция излучательной способности охватывает как длинноволновый инфракрасный (8-13 мкм), так и средневолновый инфракрасный (MWIR, 3-5 мкм) диапазон. В MWIR, тем не менее, полиэфирная пленка демонстрирует значительное поглощение из-за режима растяжения C-H связей, который не зависит от приложенного напряжения, ограничивая диапазон модуляции излучательной способности до 0.7-0.5 (2e). Из этого следует, что любые устройства, работающие в этом диапазоне длин волн, нуждаются в нестандартном защитном слое.

Другим эффектом полиэфирного слоя является повышенная излучательная способность поверхности благодаря термической экстракции полиэфиром, у которого показатель преломления больше, чем у воздуха.

Также наблюдалась модуляция излучательной способности (0.2-0.4) и в коротковолновом инфракрасном диапазоне (SWIR, 0.9-1.7 мкм). А вот модуляция в видимом спектре была незначительной из-за недостаточного легирования графена.

Улучшить модуляцию в SWIR и видимом диапазоне возможно за счет использования ионной жидкости с большим электрохимическим окном, которая будет совместима с текстилем.

Учитывая, что ткань должна быть растяжимой и гибкой, используемые графеновые элементы должны действовать соответственно. Однако многослойный графен не растягивается и не гнется ввиду механического воздействия. Потому в разработке была использована нестандартная изогнутая конструкция графена, что обеспечило уровень деформации до 60%.


Изображение 3

Решить все проблемы с гибкостью и механическим напряжением можно за счет использования массивов электродов, в не единого элемента. На показан пример такого варианта конфигурации с массивом из 25 индивидуально адресуемых электродов и датчиком термобатареи. В качестве активного слоя использовался большой цельный лист многослойного графена на хлопчатобумажной ткани (3b). Каждый электрод контролирует излучательную способность площадью 2х2 см. Внешняя электронная схема была запрограммирована реагировать на тепловую сигнатуру от датчика. Графики 3c и 3d показывают сигналы датчика и кажущуюся температуру активного пикселя (область контроля 2х2).

Мультипиксельное текстильное устройство отображает буквы C или H (обозначающие cold и hot), настраивая излучательную способность соответствующих пикселей, реагирующих на наличие/отсутствие горячего объекта над датчиком. На 3e показаны тепловые изображения работы устройства при взаимодействии с рукой человека.

Далее ученые провели фактическое практическое испытание устройства, внедренного в обычную футболку. Из-за естественной температуры тела, в условиях окружающей среды, человеческое тело излучает около 100 Вт инфракрасного света в основном в LWIR диапазоне. Этот спектральный диапазон также совпадает с окном атмосферного пропускания, которое позволяет распространять излучаемый LWIR свет на большие расстояния.

Устройство для футболки было изготовлено путем ламинирования пленки графен/полиэфир размером 6х6 см непосредственно на поверхности футболки из 100% хлопка и сеткой из нержавеющей стали на обратной стороне ().


Изображение 4

Для передачи закодированного сигнала был использован микроконтроллер, который был запрограммирован передавать буквы N, G и I азбукой Морзе. Тире и точки создавались путем подавления кажущейся температуры на длительное (9 с) и короткое (3 с) время.

На 4b и показаны инфракрасные снимки футболки в состояниях с высокой и низкой излучательной способностью, а шкала справа показывает зарегистрированную с расстояния в 3 м кажущуюся температуру.

Использование микроконтроллера позволяет строить более сложные схемы на текстиле, что, в свою очередь, обеспечивает более безопасные протоколы связи, например, инициирование связи при получении внешних запускающих стимулов. А человеческое тело в данной конфигурации может служить источником энергии.

В данном опыте скорость связи с использованием одной заплатки ограничена процессом интеркаляции/деинтеркаляции, который масштабируется в зависимости от площади устройства.


Передача букв N, G и I азбукой Морзе.

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


Увеличение отражающей способности полиэфирного устройства в ближнем инфракрасном диапазоне.

Разработанное устройство требует низкого напряжения (~ 3 В) и совсем немного энергии (5.5 х 10-4 мАч/см2 на одно событие интеркаляции, что соответствует плотности заряда ~ 1014 см-2 для каждого слоя графена). Следовательно, обычная дисковая батарейка на 1000 мАч может активировать устройство размером с футболку (1 м2) около 180 раз. Кроме того, энергия потребляется исключительно во время цикла зарядки (интеркаляции). А средняя мощность в режиме ожидания практически равна нулю, что позволяет значительно продлить использования одного устройства без замены внешнего источника питания. Это, конечно, если не рассматривать идею с использованием человека в качестве источника энергии.

Для более подробного ознакомления с нюансами исследования рекомендую заглянуть в доклад ученых и дополнительные материалы к нему.

Эпилог


Мода переменчива, как и погода. А вот наука, хоть иногда и кажется хаотичной, но все же следует одним и тем же естественным законам.

В данном труде ученые использовали оптические свойства графена в своей разработке графенового адаптивного оптического текстиля. Это устройство позволяет не только модулировать его температуру, но и дает возможность лучше понять термические и механические свойства графена. Успешная демонстрация модуляции оптических свойств на различных типах текстиля может дать толчок более широкому использованию волокнистых архитектур. Спектр применения подобных технологий не ограничивается элементами гардероба, она может быть крайне полезна и в технологиях связи, и даже в адаптивных скафандрах.

Сами же ученые намерены шагнуть еще дальше. В дальнейшем они планируют использовать свою разработку в спутниках на околоземной орбите. Спутники, как никто другой, испытывают экстремальные перепады температуры: в тени Земли они замерзают, а обращаясь к Солнцу очень нагреваются. Использование данной технологии в теории позволяет получить контроль над тепловым излучением, следовательно, и над температурой самого спутника. От обычной футболки к спутникам на орбите остается лишь надеяться, что амбиции ученых будут подкреплены успешными результатами их дальнейших исследований.

Благодарю за внимание, оставайтесь любопытствующими и хорошей всем рабочей недели, ребята. :)

Немного рекламы


Спасибо, что остаётесь с нами. Вам нравятся наши статьи? Хотите видеть больше интересных материалов? Поддержите нас, оформив заказ или порекомендовав знакомым, облачные VPS для разработчиков от $4.99, уникальный аналог entry-level серверов, который был придуман нами для Вас: Вся правда о VPS (KVM) E5-2697 v3 (6 Cores) 10GB DDR4 480GB SSD 1Gbps от $19 или как правильно делить сервер? (доступны варианты с RAID1 и RAID10, до 24 ядер и до 40GB DDR4).

Dell R730xd в 2 раза дешевле в дата-центре Equinix Tier IV в Амстердаме? Только у нас 2 х Intel TetraDeca-Core Xeon 2x E5-2697v3 2.6GHz 14C 64GB DDR4 4x960GB SSD 1Gbps 100 ТВ от $199 в Нидерландах! Dell R420 2x E5-2430 2.2Ghz 6C 128GB DDR3 2x960GB SSD 1Gbps 100TB от $99! Читайте о том Как построить инфраструктуру корп. класса c применением серверов Dell R730xd Е5-2650 v4 стоимостью 9000 евро за копейки?
Подробнее..

Нет реальности без боли электронный эквивалент рецепторов кожи человека

11.09.2020 10:13:26 | Автор: admin


Кожа это не только самый большой орган нашего тела, но и самая большая сенсорная система, ежесекундно собирающая информацию о внешних раздражителях и уровнях их воздействия на наш организм. С механической точки зрения, воссоздание кожи человека не является столь сложной задачей, но это будет лишь искусственный защитный слой, лишенный нейронной активности. Ученые из Мельбурнского королевского технологического университета (Мельбурн, Австралия) разработали систему искусственных датчиков, имитирующих различные сенсоры кожи человека. Что потребовалось для создания столь сложной имитации, каков принцип работы устройства, какие раздражители оно способно воспринимать, и где может быть применена данная разработка? Ответы на эти и другие вопросы мы найдем в докладе ученых. Поехали.

Основа исследования


Сенсорная система кожи человека может быть разделена на несколько подсистем, каждая из которых отвечает за определенные раздражители. Самыми распространенными и важными считаются рецепторы давления (тельца Пачини), температуры (терморецепторы) и боли (ноцицепторы).

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

Логично, что подобная биологическая система крайне сложна для воспроизведения, даже при учете современных тактильных сенсоров и КМОП (комплементарная структура металл-оксид-полупроводник).

Существуют разработки, в которых реализован искусственный ноцицептор на основе диффузионного мемристора*, который может демонстрировать нормальное состояние, состоящее из напряженной и релаксационной стадии ноцицептора, а также аномальное состояние с аллодинией* и гипералгезией* ноцицептора, использующего внешние стимулы в качестве напряжения.
Мемристор* пассивный элемент в микроэлектронике, способный изменять свое сопротивление в зависимости от прошедшего через него заряда.
Аллодиния* аномальная боль, вызванная раздражителем, который обычно не вызывает болевых ощущений (например боль при легком прикосновении).

Гипералгезия* аномально высокая чувствительность организма к болевым стимулам.
По словам ученых, данные разработки крайне важны, поскольку механизм переключения мемристора зависит от проводящих нитей, которые имеют примерно субнанометровый диаметр. Используя термоэлектрический модуль и пьезоэлектрический модуль давления, можно успешно достичь напряженного и релаксационного состояния среди четырех основных функций ноцицептора.

На данный момент подобная методика используется для создания искусственного глаза, но реализация ее в формате искусственной кожи пока не была достигнута.

В данном труде ученые демонстрируют рабочий прототип искусственных электронных рецепторов, которые имитируют тельце Пачини, терморецептор и ноцицептор. Достичь этого удалось за счет комбинации нескольких функциональных составляющих:

  • мемристор для принятия решений на основе титаната стронция SrTiO3 (STO) с дефицитом кислорода;
  • датчик давления на основе золота на растяжимом эластомере (полидиметилсилоксан, т.е. PDMS);
  • температурный триггер на основе оксида ванадия (VO2) с фазовым переходом.

Основным отличием данной концепции от предыдущих является отсутствие необходимости в отдельных и сложных термоэлектрических модулях и пьезоэлектрических датчиках давления для практической реализации соматосенсоров. Следовательно, в разработке используются исключительно недорогие и легкодоступные тонкие оксидные пленки, а также носимые датчики давления на биосовместимом PDMS.

Результаты исследования


Прежде чем создавать что-либо, необходимо обдумать что и как будет работать в финальном варианте. С целью создания искусственных рецепторов кожи была разработана гипотетическая основа для реализации функциональных телец Пачини, терморецепторов и ноцицепторов (схема ниже).


Изображение 1

Человеческие соматосенсоры через спинные рога (выступы серого вещества) соединены со спинным мозгом, который передает информацию в мозг (). Существуют специальные пути для определения давления (синий цвет на 1a) и температуры (красный цвет на 1a).

Тельца Пачини это слои мембран, заполненных жидкостью. Отпечатки пальцев являются ярким примером тельца Пачини. Когда на тело оказывается местное давление, часть тельца деформируется, вызывая сдвиг химических ионов (например, натрия или калия) и, как следствие, возникает рецепторный потенциал на нервном окончании кожи. Этот рецепторный потенциал при достижении достаточной энергии (порог) генерирует электрический импульс внутри тельца, который проходит через центральную нервную систему, чтобы активировать двигательный ответ через нервные волокна (синий цвет на 1a).

Когда температура кожи поднимается выше 30 C, терморецептор определяет тепло и запускает потенциалы действия. Частота возбуждения увеличивается с увеличением температуры стимула, пока не достигнет значения насыщения. Кроме того, тепловые ноцицепторы, которые улавливают болевые сигналы, начинают срабатывать при температуре около 45 C. Эти клетки специализируются на обнаружении вредоносного тепла и ожогов.

Когда вредоносный стимул получен тепловым нейроном, расположенным на свободном нервном окончании, электрический ответ отправляется ноцицептору, чтобы сравнить, пересекает ли амплитуда стимула пороговое значение для генерации потенциала действия и отправки в центральную нервную систему через спинной мозг (красный цвет на ).

Чтобы создать аналогичные искусственные рецепторы, были использованы датчики давления на основе золота и PDMS, которые переключаются между состоянием низкого сопротивления (LRS) и состоянием высокого сопротивления (HRS) без и с приложенным давлением для имитации тельца Пачини (1b и 1c).

Чтобы воспроизвести поведение терморецепторов и ноцицепторов, использовался фазовый переход VO2, который может переходить от HRS при комнатной температуре к LRS при температуре выше температуры перехода (68 C).

Кроме того, в качестве элемента принятия решений для оценки пороговых уровней была использована резистивная коммутационная память на основе STO (титанат стронция).

Для искусственного тельца Пачини, когда нет определяемого давления, ток через мемристор принятия решения (I1) недостаточен из-за напряжения смещения, чтобы инициировать двигательную реакцию (1b). При приложении давления датчик переходит в режим HRS, блокирующий I2, что позволяет максимальному току проходить через мемристор. Из-за более высокого I1 мемристор на основе STO переключается на LRS. Следовательно, через тельце протекает более высокий ток, вызывающий двигательную реакцию ().

В случае терморецептора и ноцицептора VO2 может демонстрировать изменение сопротивления на три-четыре порядка при температуре перехода. Если же температура ниже температуры перехода, то VO2 является изолятором.

Таким образом, через рецептор протекает незначительное количество тока, а напряжение, которое появляется на мемристоре, недостаточно для его включения (1d). По достижении температуры перехода VO2 переключается на LRS, в результате чего на мемристоре появляется более высокий потенциал, что вызывает его переключение на LRS. Когда и VO2, и STO находятся в LRS, через рецептор протекает повышенный ток (1e).

Искусственное тельце Пачини


После создания концепции будущего устройства ученые приступили к поэтапной реализации. На первом этапе было создание искусственного тельца Пачини, для чего использовался мемристор на основе кислорододефицитного STO со стековой структурой: Pt (100 нм) / Ti (10 нм) / STO (55 нм) / Pt (25 нм) / Ti (7 нм) и подложка SiO2.

Нижний слой Ti используется как адгезионный слой нижнего слоя Pt, а верхний слой Ti используется как резервуар для кислорода, а также как адгезионный слой верхнего слоя Pt. Нижний слой Pt принимает участие в процессе переключения, а вот верхний служит в качестве инертного материала, предотвращающего TiO2 из-за воздействия кислорода окружающей среды.

Архитектура датчика давления вдохновлена биологическим тельцем Пачини, которое имеет спиральную форму с шириной дорожки и зазором 100 мкм. Диаметр всей спирали целиком составляет 7.8 мм. Для создания сенсора на PDMS толщиной 300 мкм был нанесен Au (200 нм) / Cr (20 нм).


Изображение 2

На изображении показан искусственный эквивалент тельца с интеграцией мемристора и датчика давления.

Сеть датчиков давления работает таким образом, что позволяет рецептору активировать мемристор, который работает как компонент принятия решений.

В биологических системах при достижении достаточного рецепторного потенциала компонент принятия решения может создать электрический импульс для активации мотора центральной нервной системы. Чтобы создать реплику этой функции с определенным пороговым значением, от датчика давления требуется восприятие определенного спектра значений давления. Для простоты демонстрации работоспособности системы ученые решили упростить этот момент до всего двух значений: есть сильное давление и давления нет вообще.

Фиксированное сопротивление в 100 кОм было выбрано для ограничения тока, проходящего через сеть датчика давления, которая имеет сопротивление всего 0.6 кОм. Это гарантирует, что система показывает очень низкий ток при отсутствии давления. На изображении 2b показан отклик и повторяемость автономного датчика давления.

При приложении давления датчик давления переходит в очень HRS с сопротивлением около 1 ГОм из-за деформации и трещин, которые очень часто встречаются для датчиков на основе PDMS. Из-за деформаций и трещин характеристики датчика давления могут ухудшиться после многократных циклов, однако это не мешает демонстрации самой концепции искусственного датчика. Когда давление сбрасывается, зазоры из-за трещин снова закрываются, создавая LRS, в результате чего датчик возвращается в исходное состояние.

Подобная картина наблюдается и в биологических датчиках, которые также деформируются, что приводит к сдвигу химических ионов при приложении давления.

Следует отметить, что компонент принятия решения, состоящий из мемристорного элемента на основе STO, должен быть первоначально подвергнут гальванопластике путем приложения напряжения смещения при очень низком токе 1 мкА к верхнему и нижнему электродам. Этот этап создает локализованный канал для образования проводящих нитей через STO. После этого требуется развертка напряжения для переключения устройства между состояниями HRS и LRS.

Важно и то, что без приложенного давления ток, протекающий через мемристор, недостаточен для его переключения. Однако при приложении давления датчик, содержащий ответвление, переходит в состояние HRS, что приводит к максимальному потенциалу рецептора на мемристоре (). При достижении порога рецепторного потенциала мемристор, принимающий решение, переключается из состояния HRS в состояние LRS (). В этом состоянии примененная последовательность 0 +0.85 В 0 1.12 В 0 переключает устройство в состояние LRS для положительного цикла и в состояние HRS для отрицательного полупериода (2d и ).

Чтобы перевести устройство в LRS, учитывается только положительный полупериод. В соответствии с изображением 2d, когда давление не подается, цепь датчика давления имеет общее сопротивление 100.6 кОм, тогда как параллельный компонент принятия решения (мемристор) имеет сопротивление 70 кОм. Таким образом, эквивалентное сопротивление всего тельца Пачини составляет 41.2 кОм.

Это эквивалентное сопротивление пропускает ток всего 0.02 мА через всю цепь, что можно рассматривать как расслабленное состояние. Приложение давления преобразует цепь датчика давления в состояние с чрезвычайно высоким сопротивлением 1 ГОм, в то время как сопротивление мемристора составляет всего около 2.5 кОм, изменяя эквивалентное сопротивление всего тельца Пачини примерно на 2.5 кОм. Это низкоомное состояние допускает ток 0.35 мА по всей цепи.

Таким образом, стимул давления генерирует ответный сигнал, который почти на 18 раз выше, чем в расслабленном состоянии, что может позволить центральной нервной системе инициировать свой двигательный ответ. После того как двигательная реакция завершена, для инициализации тельца Пачини можно применить обратную полярность к мемристору, используя неиспользованные электроды.

Искусственный терморецептор


Для создания терморецептора в основе мемристора была использована такая же стек-структура, как и для тельца Пачини, т.е. металл-изолятор-металл (МИМ).


Изображение 3

Часть верхнего электрона использовалась совместно с поверхностью VO2 (3а и 3b) для последовательного подключения теплового датчика. Для смещения всего устройства этот электродный слой, состоящий из Pt (100 нм) / Ti (10 нм), был нанесен на поверхность VO2. Между исходным электродом и верхним электродом мемристора поддерживалось существенное расстояние в 100 мкм.

На изображении показана схема подключения терморецептора, в которой смещение приложено через металл к тепловому датчику, а земля (GND на схеме) подключена к нижнему электроду мемристора, принимающего решение.

График 3d показывает кривую зависимости сопротивления от температуры для перехода диэлектрикметалл на тонкой пленки VO2. Очевидно, что при достижении температуры перехода наблюдается падение удельного сопротивления на четыре порядка. Очевидный тепловой гистерезис также наблюдается в циклах нагрева и охлаждения. Также было установлено отсутствие какого-либо заметного влияния температуры на резистивное переключение ().

Мемристор может показывать изменение сопротивления от 100 кОм до 2 кОм в процессе переключения. Однако для более понятного анализа было решено рассматривать сопротивление состояния HRS как 93 кОм, а сопротивление состояния LRS как 9 кОм при 80 мВ напряжения считывания (VREAD), так как при этом напряжении наблюдается максимальное отношение переключения ROFF / RON. Последовательность переключения напряжения готового автономного устройства составила: 0 +0.65 В 0 0.80 В 0. Когда та же последовательность применяется ко всему терморецептору, сопротивление уменьшается, и, следовательно, ток рецептора увеличивается с повышением температуры (3f).

Чтобы обеспечить необходимое коммутируемое напряжение, температура приемника поддерживалась на уровне 70 C. Это необходимо для гарантии того, что VO2 находится в состоянии LRS. Затем приложение напряжения смещения от 0 до 2 В полностью устанавливает и сбрасывает устройство ().


Изображение 4

Для мемристора, принимающего решение, исходное сопротивление 93 кОм намного ниже, чем HRS теплового датчика, которое составляет 11 МОм. Таким образом, частичное напряжение, которое появляется на мемристоре, принимающем решение, не может достичь порогового значения VSET для преобразования его из HRS в LRS. Следовательно, и термодатчик, и мемристор находятся в состоянии HRS, что позволяет минимальному току протекать через терморецептор [4b(i)].

Когда применяется критическая температура 70 C, сопротивление термодатчика уменьшается на четыре порядка, а частичное напряжение мемристора постепенно увеличивается до напряжения SET с увеличением отклика рецептора [4b(ii)]. Как только VSET включает мемристор, он переходит в LRS от HRS с сопротивлением 9 кОм [4b(iii)].

На этом этапе формируется максимальный рецепторный ответ. LRS памяти будет сохраняться в течение длительного времени, даже если тепловой стимул полностью отключен. Чтобы перепрограммировать мемристор, отрицательное напряжение VRESET может перевести его с LRS в HRS [4b(iv)]. Для этого можно подавать отрицательное напряжение от неиспользуемых электродов ().

Искусственный ноцицептор


Можно с уверенностью сказать, что ноцицепторы значительно отличаются от своих собратьев. Ноцицепторы имеются по всему человеческому телу и расположены на конце аксона сенсорного нейрона.

Чтобы избежать воздействия вредоносных раздражителей, ноцицептор реагирует двумя способами: нормальными и анормальными.

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

Анормальный ответ возникает, когда организм сталкивается со стимулами, близкими к порогу повреждения ноцицептора, и в этом состоянии ноцицептор работает как обычный рецептор, чтобы избежать дальнейшего повреждения. Если же травма все же была получена, то уязвимость пораженной ткани возрастает. Ноцицептивная система адаптируется к этой повышенной уязвимости, локально понижая ноцицептивный порог и облегчая ноцицептивный ответ, тем самым обеспечивая адекватную защиту тканей.

Ноцицептор демонстрирует два различных поведения в анормальных условиях: аллодиния и гипералгезия.

Аллодиния дает ответный сигнал при нижнем пороговом значении, тогда как гипералгезия создает более сильный ответный сигнал при превышении порогового значения, указывая на то, что при аномальном состоянии порога для ноцицептора нет.

Для наблюдений за поведением искусственного терморецептора в нормальных условиях, который работает как ноцицептор во время вредоносных стимулов, устройство было переключено на LRS, а показания считывались при VREAD80 мВ.

Поскольку срабатывание биологического ноцицептора в значительной степени зависит от интенсивности стимула, на искусственный ноцицептор воздействовали серией температурных стимулов с различной интенсивностью в диапазоне от 66 до 82 C (5a).


Изображение 5

График 5b демонстрирует ответный сигнал относительно интенсивности приложенного теплового стимула. Следует отметить, что ноцицептор не включается, пока температурный импульс не достигнет 68 C, что является температурой перехода используемого VO2. Таким образом, когда VO2 попадает в LRS из-за температурно-индуцированного перехода, более высокий ток начинает течь по всей цепи. Подобное повеление напоминает биологические системы, в которых ноцицептор генерирует запускающий мозг потенциал действия, когда сила стимула достигает значений выше критического.

Дальнейшее увеличение интенсивности стимула выше порогового значения приводит к большему току, что также согласуется с реакцией биологического аналога: чем выше интенсивность стимула, тем выше интенсивность ответной реакции. На графике 5c показаны стимулы нагрева и охлаждения и соответствующий ответный сигнал при 68 C.

График 5d показывает ослабление ответного сигнала с течением времени после того, как вредоносный стимул был отключен, т.е. процесс релаксации, который определяется VO2, так как нет влияния температуры на мемристор на основе STO.

Из-за тенденции к увеличению сопротивления VO2 по мере снижения температуры с течением времени искусственный ноцицептор ограничивает ток через цепь, и, следовательно, наблюдается уменьшение интенсивности ответных сигналов.

Более сильный ответный сигнал из-за более высоких стимулов требует относительно большего времени для полной релаксации. Например, ответному сигналу при 68 C требуется 100 с, чтобы достичь базового тока в 0.5 мкА, тогда как ответный сигнал при 80 C не может полностью релаксировать за 100 с.

Для наблюдения за поведением устройства в анормальных условиях, к искусственному ноцицептору был применен стимул, интенсивность которого была гораздо выше, чем в нормальных условиях.

Ноцицептор нагревали до 90 C со скоростью 20 градусов в минуту и охлаждали до 60 C, что ниже порогового значения (68 C) в нормальных условиях.

Далее следовал повторный нагрев с 60 до 90 C, необходимый для определения наличия/отсутствия генерации пониженного порога и усиленного ответа, которые являются основными свойствами аллодинии и гипералгезии.

Последовательность 60 90 60 90 применялась к VO2 части рецептора и ко всему рецептору, содержащему VO2 и стек металл-изолятор-металл ().


Изображение 6

На графиках отчетливо видно, что ответный сигнал намного более линейный в VO2 () по сравнению с сигналом всего ноцицептора (6b). Это вполне ожидаемо, так как при высокой интенсивности стимула VO2 находится в почти металлическом состоянии с относительно низким сопротивлением 5 кОм после перехода. Более того, приложенное напряжение смещения VREAD (80 мВ) электрически настраивает VO2, чтобы сделать его еще более металлическим, что приводит к линейному отклику.

В то же время, когда подобное смещение VREAD появляется на всем ноцицепторе, максимальное падение напряжения происходит на мемристоре, который находится в состоянии LRS (9 кОм). На этом этапе напряжения на VO2 недостаточно, чтобы показать линейный отклик. Следовательно, данное поведение является нелинейным.

На схеме показано поведение при аллодинии и гипералгезии. В биологической системе интенсивность ответа выше в аномальном состоянии для подпороговой (аллодиния) и сверхпороговой (гипералгезия) интенсивности стимула.

На 6d показан отклик по отношению к двум циклам нагрева последовательности с последовательностью 60 90 60 90 C. Тут видно, что отклик для второго цикла нагрева усиливается, а порог снижается.

В поведении искусственного рецептора четко видны аллодиния ниже пороговой интенсивности и гипералгезия выше пороговой интенсивности (70 C). Из этого следует, что снижая порог и усиливая интенсивность ответа, ноцицептор активирует и усиливает защитные реакции, такие как отстранение или избегание острых болезненных стимулов.


Демонстрация принципа работы электронного эквивалента тельца Пачини.

Для более детального ознакомления с нюансами исследования рекомендую заглянуть в доклад ученых и дополнительные материалы к нему.

Эпилог


Мозг человека является одной из самых сложных биологических систем. Но нельзя отрицать и того факта, что кожа человека не менее сложна, особенно учитывая немалый список выполняемых ею функций.

Воссоздать некоторые из функций кожи в искусственном эквиваленте не сложно, но вот касательно рецепторов, собирающих информацию об окружающей среде, так сказать не получится.

Тем не менее, ученым все же удалось достичь некоторых результатов в создании искусственных рецепторов, улавливающих давление, температуру и боль.

Как заявляют авторы сего труда, их устройство способно различать легкое прикосновение и, например, укол иглы. На первый взгляд, это весьма банальные вещи, однако ранее такой точности не было в электронных рецепторах.

В будущем ученые, естественно, намерены продолжить свой труд, дабы расширить спектр воспринимаемых внешних стимулов, что позволит сделать их устройство еще более точным. Подобные разработки однозначно найдут свое применение не только в протезировании, но и даже в робототехнике.

Благодарю за внимание, оставайтесь любопытствующими и отличных всем выходных, ребята! :)

Немного рекламы


Спасибо, что остаётесь с нами. Вам нравятся наши статьи? Хотите видеть больше интересных материалов? Поддержите нас, оформив заказ или порекомендовав знакомым, облачные VPS для разработчиков от $4.99, уникальный аналог entry-level серверов, который был придуман нами для Вас: Вся правда о VPS (KVM) E5-2697 v3 (6 Cores) 10GB DDR4 480GB SSD 1Gbps от $19 или как правильно делить сервер? (доступны варианты с RAID1 и RAID10, до 24 ядер и до 40GB DDR4).

Dell R730xd в 2 раза дешевле в дата-центре Equinix Tier IV в Амстердаме? Только у нас 2 х Intel TetraDeca-Core Xeon 2x E5-2697v3 2.6GHz 14C 64GB DDR4 4x960GB SSD 1Gbps 100 ТВ от $199 в Нидерландах! Dell R420 2x E5-2430 2.2Ghz 6C 128GB DDR3 2x960GB SSD 1Gbps 100TB от $99! Читайте о том Как построить инфраструктуру корп. класса c применением серверов Dell R730xd Е5-2650 v4 стоимостью 9000 евро за копейки?
Подробнее..

Квантовый нанотермометр измерение температуры нематоды длиной 1 мм

18.09.2020 10:20:14 | Автор: admin


Одним из основных показателей состояния биологической системы является температура. Если у человека развивается какая-то инфекция, то температура его тела повышается (как правило, но не всегда), что является признаком ответной реакции иммунной системы на угрозу. Другими словами, по температуре можно определить примерное состояние организма. Проблема в том, что человек большой (буквально), а вот, например, нематоды в длину всего лишь около 1 мм. Измерить температуру столь малого организма было крайне сложно, однако ученые из университета Осаки (Япония) разработали методику, позволяющую решить эту проблему. Какие средства были использованы для реализации нанотермометра, что показали практические опыты, и где можно использовать данную разработку? Ответы на эти вопросы мы найдем в докладе ученых. Поехали.

Основа исследования


Температура тела живого организма варьируется в зависимости от степени воздействия внутренних и внешних факторов. Мы привыкли, что температура окружающей среды напрямую влияет на температуру холоднокровных, посему ее значения меняются с завидной регулярностью. Однако даже у теплокровных при нормальных физиологических условиях наблюдаются температурные колебания, которые можно связать с гомеостатической терморегуляцией и энергетическим обменом.

Другими словами, тут отлично подходит шутка: я не бездельничаю, я очень занятой человек на клеточном уровне. Если точно измерить температуру и ее динамику в субмикронном масштабе, то можно получить много информации касательно клеточной и молекулярной активности. Проблема в том, что с уменьшением объекта измерения увеличивается сложность его проведения (сложно засунуть в нематоду обычный термометр из аптеки).

Авторы исследования отмечают, что обычные электрические термометры не имеют субмикронного разрешения, а термография в ближнем инфракрасном диапазоне обычно помогает определять температуру поверхности биологических образцов, но не внутреннюю температуру.

Конечно, сейчас уже есть светоизлучающие нанотермометры (например, термочувствительные молекулярные зонды), которые способны преодолеть это ограничение. Но у такой методики также есть недостатки. Основной это долговременная устойчивость, а точнее ее отсутствие. Подобные устройства не могут точно измерять изменения температуры, которые протекают длительное время (скажем пару часов). Не говоря уже о токсичном воздействии на образец со стороны такого термометра.

В данном труде ученые описывают концепцию наноалмазного (ND от nanodiamond) квантового термометра, который обладает высокой точностью, устойчивостью и низкой токсичностью. Принцип его работы таков: датчик считывает температуру как сдвиг частоты оптически детектируемого магнитного резонанса (ODMR от optically detected magnetic resonance) дефектных центров азотных вакансий (NV от nitrogen-vacancy), который в основном возникает из-за теплового расширения решетки. Сенсорное ядро NV глубоко встроено в решетку алмаза и невосприимчиво к различным биологическим факторам окружающей среды. Внедрение этого квантового датчика в более сложные организмы позволяет считывать их тепловую активность на конкретном участке в режиме реального времени. Но процесс реализации такой техники сопряжен с рядом сложностей.


Нематода (круглый червь) вида Caenorhabditis elegans.

Многоклеточные модельные организмы, такие как черви Caenorhabditis elegans, нуждаются в специальной камере, способной вместить тело миллиметрового размера, а сами образцы необходимо быстро анализировать, чтобы сохранить их физиологическое состояние. Квантовые ND термометры движутся намного быстрее, чем в культивируемых клетках, даже если тело обезвожено, что требует использования алгоритма быстрого отслеживания частиц. Кроме того, позиционное перемещение ND и сложная структура тела вызывают существенные колебания обнаруженной интенсивности флуоресценции, что, вероятно, вызовет артефакты измерения температуры. Решение этих проблем на данном этапе исследования сопряжено с подгонкой устройство под индивидуальные особенности анализируемого образца. Вопрос универсальности и легкости в настройке будущего нанотермометра планируется рассматривать в дальнейших работах, а пока внимание было уделено самой концепции и основным принципам работы.


Небольшой ролик, рассказывающий о нематодах.

Результаты исследования


Основой нанотермометра является конфокальный флуоресцентный микроскоп, оборудованный установкой для микроволнового облучения (1А).


Изображение 1

ODMR азотных вакансий можно измерить как уменьшение интенсивности лазерно-индуцированной флуоресценции при применении спин-резонансного микроволнового возбуждения, поскольку спиновое возбуждение активирует нефлуоресцентный путь релаксации из возбужденного состояния в основное состояние ().

Камера, куда помещаются образцы, представляет собой одноразовую чашу со стеклянным дном, интегрированную в антенну, которая обеспечивает оптический доступ большой площади (диаметр 12 мм) и простоту использования (1C), подходящую для деликатных образцов, таких как стволовые клетки. Время от захвата червя Caenorhabditis elegans до начала фактического измерения составляет всего 15 минут. Это помогает сохранить жизнеспособность червя и способствует получению большего объема данных о его состоянии.

Кроме того, данная система эффективно объединяет быстрое отслеживание частиц и высокоточную оценку температуры в реальном времени по ODMR смещению NV центров.

При отслеживании частиц система измеряет интенсивность флуоресценции ND вдоль осей xyz микроскопа и фокусируется на соответствующем максимуме флуоресценции каждые 4 секунды (возможен более короткий интервал отслеживания), в течение которых температура оценивается со временем выборки от 0.5 до 1.0 секунды. ().


Изображение 2

Методов квантовой термометрии существует несколько, однако в данном труде был использован метод четырехточечных измерений ODMR. Этот метод предполагает, что количество фотонов, зарегистрированных на всех четырех выбранных частотах, линейно масштабируется в соответствии с изменениями обнаруженной интенсивности флуоресценции.

Однако, было обнаружено, что каждый последующий фотон показывает разницу в светочувствительности около 0.5%, что фактически создает существенные артефакты в оценке частотного сдвига (т.е. 300 кГц, что соответствует нескольким градусам Цельсия), особенно при низко-фотонном режиме.

Эти артефакты, скорее всего, возникают из-за зависимой от оптической мощности асимметрии спектра ODMR. Для точного измерения температуры сложных оптических динамических систем (т.е. биологических систем) от подобных артефактов необходимо избавляться. Поэтому в метод четырехточечных измерений был добавлен фильтр коррекции ошибок.

Для оценки работы системы, сопряженной с коррекцией ошибок, в реальном времени были проведены измерения температуры ND во время ступенчатых тепловых событий. Резкие изменения температуры использовать нельзя было, так как внезапные изменения температуры вызывают большую расфокусировку фокальных пятен и связанные с ними флуктуации интенсивности флуоресценции.

На 2B показаны временные профили общего количества фотонов (Itot) и температурная оценка ND (TNV), когда температура образца (TS) изменяется от 44.3 30.4 44.3 с шагом в 2.8. Система точно выдает TNV, соответствующий TS, при этом положение фокуса существенно перемещалось, особенно вдоль оси z на расстояние более 30 мкм (2C).

При шаге в 3 проявляется позиционный сдвиг по оси z на 6 мкм в течение 3-4 минут, но скорость слежения достаточно высока, чтобы следовать динамике 105 нм/с в течение 96 минут ().

Кроме того, TNV четко демонстрирует антикорреляцию с Itot. Статистическое исследование этого типа температурной зависимости определяет средние значения для SD: Itot-1dItot/dT = -3.9 0.7 %/С и dD/dT = 65.4 5.5 кГц/С (2D). При этом точность измерения температуры составляет 0.29 и < 0.6 C, соответственно, что дает чувствительность 1.8 C/Гц.

После достижения надежной и точной термометрии в реальном времени в рамках этапа разработки, был проведен тестовый локальный мониторинг температуры на живых червях.


Изображение 3

На снимке показаны ND внутри анестезированных червей, помещенных рядом с микроволновыми антеннами. Эти ND хорошо диспергируются в воде за счет поверхностной функционализации полиглицерина (PG от polyglycerol) и вводятся путем микроинъекции в гонады (половые железы подопытного червя).

На графике показан ODMR спектр одиночного ND (отмечен стрелкой на ). демонстрирует временные профили Itot и TNV за период в 1 час при изменении температуры TS.

Сначала проводилось измерение Tobj при 33.2 С, через 6 минут было выполнено уменьшение до 25.3 С. В результате Tobj достиг уровня 28.6 С на 35.2 минуте. TNV показал точное изменение температуры между двумя стационарными состояниями: 33.2 и 28.6 С.

Отображение реальной динамики температуры внутри червей между этими двумя стационарными состояниями отображается за счет того, что TNV всегда отстает от TS и демонстрирует немного заниженный отклик из-за конечной теплоемкости объектива микроскопа и окружающей среды. Itot также показывает постепенные изменения интенсивности флуоресценции, вызванные температурой.

Отслеживание частиц также выполнялось на удовлетворительном уровне (). В течение 0-15 минут в подсчитанных фотонах появляются частые всплески, возникающие из-за позиционных флуктуаций ND приблизительно на 400 нм в течение нескольких секунд.

Результаты теста отчетливо свидетельствуют о высокой точности измерения температуры внутри наноразмерной биологической системы в реальном времени. Далее было решено провести дополнительные тесты, перед которыми подопытные черви прошли фармакологическую обработку с помощью C10H5F3N4O (FCCP от карбонил цианид-4- (трифторметокси) фенилгидразон), вызывающую неподвижный термогенез (грубо говоря, повышение температуры ввиду увеличения метаболизма и без дополнительной мышечной активности).


Изображение 4

На снимке показаны ND у червей, стимулированных FCCP. А на графике показан временной профиль TNV ND, отмеченного стрелкой на снимках.

На седьмой минуте после начала измерения был использован раствор FCCP. На 32-ой минуте TNV начинает постепенно увеличиваться, а на 48-ой наблюдается еще большее дополнительное увеличение, когда уровень изменения температуры повышается от 4 до 7 С. Состояние повышенной температуры длилось около 80 минут.

Во время стимуляции ND медленно перемещаются на несколько микрометров в течение часа, что подтверждает результаты отдельных экспериментов, в которых ND непрерывно наблюдались под микроскопом.

Контрольная группа червей ( и 4D), которым не вводили FCCP, показала равномерный отклик TNV во время всего теста без каких-либо явных изменений температуры.

Для дополнительного подтверждения того, что FCCP реально вызывает повышение температуры тела червей, была проведена количественная оценка червей с помеченными ND как в контрольной, так и в подопытной группе (). График явно говорит о повышении температуры у червей из подопытной группы по сравнению с контрольной.

Другой контрольный эксперимент, в котором буферный раствор не добавлялся, а TNV отслеживался статически, показывает, что добавление допанта вызывает колебания TNV на определенном уровне либо из-за изменения температуры, либо из-за артефактов сдвига ODMR. Однако наблюдение подобного сдвига невозможно при добавлении FCCP, что дополнительно подтверждает повышение температуры за счет FCCP у подопытной группы червей (4F).

Для более детального ознакомления с нюансами исследования рекомендую заглянуть в доклад ученых и дополнительные материалы к нему.

Эпилог


В данном исследовании ученым удалось разработать методику, позволяющую точно измерить температуру внутри наноразмерной биологической системы в реальном времени. Утрировано говоря, им удалось измерить температуру тела червя Caenorhabditis elegans, длина которого составляет примерно 1 мм.

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

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

Благодарю за внимание, оставайтесь любопытствующими и отличных всем выходных, ребята! :)

Немного рекламы


Спасибо, что остаётесь с нами. Вам нравятся наши статьи? Хотите видеть больше интересных материалов? Поддержите нас, оформив заказ или порекомендовав знакомым, облачные VPS для разработчиков от $4.99, уникальный аналог entry-level серверов, который был придуман нами для Вас: Вся правда о VPS (KVM) E5-2697 v3 (6 Cores) 10GB DDR4 480GB SSD 1Gbps от $19 или как правильно делить сервер? (доступны варианты с RAID1 и RAID10, до 24 ядер и до 40GB DDR4).

Dell R730xd в 2 раза дешевле в дата-центре Equinix Tier IV в Амстердаме? Только у нас 2 х Intel TetraDeca-Core Xeon 2x E5-2697v3 2.6GHz 14C 64GB DDR4 4x960GB SSD 1Gbps 100 ТВ от $199 в Нидерландах! Dell R420 2x E5-2430 2.2Ghz 6C 128GB DDR3 2x960GB SSD 1Gbps 100TB от $99! Читайте о том Как построить инфраструктуру корп. класса c применением серверов Dell R730xd Е5-2650 v4 стоимостью 9000 евро за копейки?
Подробнее..

Что такое температура и как её учитывать в молекулярном моделировании? Реализация на CUDA

06.01.2021 20:05:30 | Автор: admin

Методы численного моделирования молекулярных систем, такие как молекулярная динамика, рассматривают эти системы как механические (что-то вроде набора шариков на пружинках). Однако, в отличие от механических систем, для молекулярных существует понятие температура. Вещество не может существовать без температуры, а температура без вещества (на счет последней части утверждения есть и другое мнение). Из опыта мы знаем, что очень многие свойства вещества кардинально зависят от температуры, и, естественно, что её надо как-то учитывать. Для поддержания температуры в молекулярной динамике используются специальные алгоритмы-термостаты. Наиболее известные среди них это термостаты Андерсена, Берендсена и Нозе-Гувера. Все они основаны на молекулярно-кинетической теории газов, где температура есть просто величина пропорциональная среднекинетической энергии молекул. Соответственно, работа данных термостатов осуществляется путём умножения скоростей частиц на некоторую величину.

Что же не так?

  1. Масштабирование скоростей изменяет их модуль, но не направление. А ведь тепловое движение, согласно большинству утверждений, должно быть хаотическим. Оно определяет броуновское движение. Термостатирование в молекулярной динамике же не привносит никакого хаоса, частицы движутся в том же направлении, как и без учета температуры, только быстрее/медленнее. Метод остается полностью детерминированным, начальные координаты и скорости частиц полностью определяют конечное состояние.

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

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

Модель

Чтобы сконструировать корректный термостат нам нужно разобраться, а что же такое температура на самом деле? Несмотря на то, что мы этим понятием оперируем чуть ли не с дошкольного возраста, суть температуры, как мне кажется, до конца не ясна.

Соберем известные факты о температуре:

  • Температура не может быть ниже абсолютного нуля, т.е. 0 К.

  • Температура понятие интенсивное, оно не зависит от размеров системы, количества вещества и т.д.

  • Температура имеет свойство выравниваться если мы приведем в контакт два тела разной температуры (даже можно и без контакта, тела могут обмениваться излучением), то рано или поздно температура обоих тел станет одинаковой.

  • Скорости химических реакций и коэффициенты диффузии с температурой растут экспоненциально.

  • При повышении температурах сложные вещества постепенно разлагаются на более простые.

  • Агрегатные состояния с повышением температуры меняются как правило в порядке: твердое, жидкое, газообразное. При этом тепло почти всегда поглощается (единственное известное исключение гелий, его энтальпия плавления отрицательная). При очень высоких температурах начинается ионизация газа (молекулы/атомы газа теряют электроны).

  • С увеличением температуры тела обычно расширяются (газы все, твердые тела почти все, но есть исключения, см. отрицательный КТР).

  • Теплоемкость с температурой как правило растет (для некоторых жидкостей, такие как расплавленные оксиды она постоянная или чуть падает).

  • С температурой связан интересный эффект, адиабатическое расширение в пустоту (эффект Джоуля-Томсона).

  • Другой интересный эффект, сонолюминисценция воздействие ультразвуком приводит к разогреванию жидкости (воды) до температуры в десятки тысяч кельвинов. На хабре, кстати, была статья про неё.

  • С увеличением температуры растет сжимаемость, т.е. при высоких температурах чтобы сжать тело в то же число раз, нужно приложить меньшее давление. По крайней мере, это верно для расплавов.

  • Электропроводность металлов падает с температурой.

  • Спектр излучения абсолютно черного тела определяется исключительно его температурой.

Как мы видим, количественных характеристик температуры не так уж и много. Привязываться к каким-то ключевым точкам типа температуры кипения воды не имеет смысла, поскольку: они очень частные; зависят и от других условий, в первую очередь, давления. Тем не менее, есть характеристика, кажущаяся довольно перспективной. Это спектр излучения абсолютно черного тела. Оказывается, если тело ничего не отражает, а все излучение от него является собственным, то спектр этого излучения не зависит от природы объекта, а определяется исключительно его температурой, согласно формуле Планка. Тот факт, что все тела излучают и характеристики этого излучения зависят от температуры используется, например, в оптической пирометрии. Таким образом, температуру можно определить через спектр излучения. В условии равновесия с термостатом, т.е. средой, имеющей заданную температуру, излучение, испускаемое объектом, должно сравняться с излучением термостата. Тогда можно реализовать термостатирование следующим образом: атомы поглощают и испускают фотоны в случайном направлении, при этом их скорости изменяются согласно закону сохранения импульса. Интуитивно кажется, что если мы будем прибавлять/убавлять поправку к скорости атома в случайном направлении, то атом сначала наберет какую-то скорость, а затем эта скорость выйдет на насыщение, поскольку все выше вероятность того, что новая скорость будет направлена против одной из предыдущих добавок. Однако легко показать, что это не так, начнем с закона сохранения импульса:

m\vec{v_1}=m\vec{v_0}+\vec{p_\gamma},\qquad\qquad\qquad (1)

где m и v масса и скорость атома, индексы 1 и 0 означают состояния до и после поглощения фотона, p импульс фотона. Перейдем к скоростям:

\vec{v_1}=\vec{v_0}+\frac\epsilon{mc}\vec{r},\quad\quad\quad(2)

где энергия фотона, c скорость света в вакууме, r единичный вектор случайного направления. Из этого, квадрат скорости равен:

v_1^2 = (\vec{v_0} + \frac{mc}\vec{r})^2 = v_0^2 + 2 \frac{mc}\vec{r}\vec{v_0} + \frac{^2}{m^2c^2}, \quad\quad\quad(3)

или, используя определение скалярного произведения векторов:

v_1^2 = v_0^2 + 2 \frac{mc}v_0 \cos \phi + \frac{^2}{m^2c^2},\quad\quad(4)

где угол между начальной скоростью атома и направлением движения поглощенного фотона. Косинус принимает значения от -1 до 1, а в среднем ноль, тогда мы видим, что среднее приращение квадрата скорости (что пропорционально кинетической энергии)

\langle\Delta v^2\rangle = \frac{^2}{m^2c^2}\quad\quad(5)

положительное. Т.е. действуя таким образом, мы будем накачивать систему кинетической энергией, независимо от того будем излучать или поглощать фотоны, поскольку знак будет стоять перед косинусом, который в среднем нулевой. Чтобы скомпенсировать неограниченный рост скоростей в системе нужно, чтобы среднее изменение квадрата скорости атома при излучении фотона было равно среднему значению изменения квадрата скорости при поглощении фотона (5) с обратным знаком:

\langle\Delta v^2\rangle = -\frac{^2}{m^2c^2}\quad\quad(6)

для этого косинус между направлением движения атома и излучаемым фотоном, в среднем, должен равняться:

\langle\cos\phi\rangle = -\frac{}{2mcv_0}\quad\quad(7)

Это выражение устанавливает связь между скоростями атомов и энергией фотонов. Можно задать косинус равномерной случайной величиной, так чтобы указанное выше значение было серединой интервала, а -1 его нижней границей.

Вторая проблема возникает с законом сохранения энергии. Если мы допустим, что вся энергия фотона полностью перейдет в кинетическую энергию атома, то должно выполняться равенство:

\frac12mv_1^2 = \frac12mv_0^2 + ,\quad\quad(8)

однако, если мы попытаемся решить уравнение в системе с законом сохранения импульса (2), мы получим нефизичный результат. Например, покоящийся атом должен разогнаться до световой скорости. Причина этого заключается в том, при одном порядке кинетической энергии атома и энергии фотона их импульсы отличаются на несколько порядков. Решение напрашивается следующее: импульс сохраняем, вычисляем изменение кинетической энергии, оставшийся избыток энергии фотона добавляем к какой-нибудь другой энергии атома, назовём её тепловой (обозначим U):

U_1 + K_1 = U_0 + K_0 + ,\quad\quad(9)

где K кинетическая энергия. Остался вопрос, фотоны какой энергии излучать атому? Давайте немного упростим модель, пусть тепловое излучение описывается не формулой Планка, а будет монохроматическим, с длиной волны / частотой / энергией однозначно определяемыми температурой. Для этого подходит, например, закон Вина, определяющий длину волны, дающую наибольший вклад в спектр излучения:

[м] = 0.002898 / T,\quad\quad(10)

где длина волны в метрах, от длины волны несложно перейти к энергии фотона, дробь при этом перевернется и в общем случае можно записать:

 = T,\quad\quad(11)

где некоторый постоянный коэффициент. Т.е. температура T, это такое состояние вещества, когда излучаются фотоны энергией T.

Кроме того, у нас есть ещё одна численная характеристика, касающаяся температуры. Мы совершенно точно знаем, какое количество теплоты (энергии/тепловой энергии) нужно сообщить телу, чтобы нагреть его до заданной температуры от нуля кельвинов. В условиях постоянного давления эта величина называется приращением энтальпии / инкрементом энтальпии / энтальпией (не путать с энтальпией образования), её значение зависит не только от температуры, но и от природы вещества. Если приводить образную аналогию, то вещество можно представить как сосуды с некой жидкостью, где её уровень это температура, а объём энтальпия. Производная энтальпии по температуре называется изобарной теплоемкостью, cp. Для идеального газа (и, с некоторым приближением, для благородных газов) cp постоянна, тогда:

H = c_p T.\quad\quad(12)

Объединяя (11) и (12) можно исключить температуру и получить энергию фотона как функцию от текущей тепловой энергии, запасенной в веществе:

 = \frac{c_p} H.\quad\quad(13)

Мы получили, что энергия излучаемых фотонов определяется энтальпией вещества. Продолжая нашу аналогию с сосудами, можно представить, что они теряют жидкость пропорционально высоте столба. Напоминаю, что выражение (13) справедливо только для систем с постоянной теплоемкостью, для реальных систем нужно находить температуру как функцию от энтальпии и подставлять её в (11). При переходе к числам меня ждал сюрприз, выразив из (10) значение и подставив его в (13), а cp взяв равным 5/2k (изобарная теплоемкость одноатомного идеального газа) получилось, что коэффициент перед энтальпией больше 1, а это означает, что система должна излучать больше энергии, чем имеется, что, конечно, не физично. Немного покопавшись, я нашел, что при переходе от длины волны к частоте (и затем к энергии) нужно использовать несколько отличное выражение, чем (10). В конечном итоге, с использованием данных по теплоемкости гелия, у меня получилось, /cp 0.9, что я и взял для своих численных экспериментов. В качестве энтальпии я взял величину U, описанную выше, формула (9).

Немного кода

Основная часть кода базируется на изложенном ранее. Однако вместо термостата Нозе-Гувера используется вышеописанная схема. Функция термостатирования:

__global__ void tstat_radi(int atPerBlock, int atPerThread, cudaMD* md){    __shared__ int indEng;    __shared__ float engTemp;    if (threadIdx.x == 0)    {        indEng = atomicAdd(&(md->curEng), 1);   // get current index in photon energies array        engTemp = 0.f;    }    __syncthreads();    uint4 randVar = md->ui4rnd;    float pe;        // photon energy    int i, e0, ei;// indexes    // atomic range    int id0 = blockIdx.x * atPerBlock + threadIdx.x * atPerThread;    int N = min(id0 + atPerThread, md->nAt);    e0 = indEng * atPerBlock + threadIdx.x * atPerThread - id0;    for (i = id0; i < N; i++)    {        ei = e0 + i;        if (ei >= md->nAt)            ei -= md->nAt;        pe = md->engPhotons[ei];        adsorb_rand_photon(&(md->vls[i]), &(md->engs[i]), md->masses[i], pe, randVar, md);        radiate_photon(&(md->vls[i]), &(md->engs[i]), md->masses[i], randVar, md);    }    atomicAdd(&engTemp, teng);    if (threadIdx.x == 0)    {        rnd_xor128(md->ui4rnd);        atomicAdd(&(md->engTemp), engTemp);    }}

Функция принимает на вход число атомов в потоке и блоке, а также указатель на структуру, где хранятся все МД-данные. Для каждого атома в обрабатываемом диапазоне находится энергия случайного фотона, вызываются функции его поглощения и функция излучения. Энергия поглощаемого фотона берется из предварительно сгенерированного массива энергий, которое соответствует Планковскому распределению (эта часть вычисляется на хосте и загружается в массив). Я здесь использовал ту фичу CUDA, что порядок выполнения блоков не определён. Первый запустившийся блок берет первый кусок массива энергий фотонов и т.д., таким образом получается что-то вроде рандомизации. Таким же образом я хотел сделать и направления фотонов, но рандомизация выходит слабенькой, поскольку порядок выполнения блоков хоть и не определен, но все-таки даёт очень небольшой разбег, для направлений это критично, даёт заметные артефакты.

Функция поглощения фотона принимает на вход указатель на скорость атома, указатель на тепловую энергию атома / энтальпию, массу атома, энергию фотона, переменную для генератора случайных чисел и указатель на главную структуру:

__device__ void adsorb_rand_photon(float3 *vel, float *int_eng, float mass, float eng, uint4 &rand_var, cudaMD *md){    float u0 = *int_eng;    float v02 = float3_sqr(*vel);   // square of initial velocity    float3 rand_vect = rand_uvect(rand_var, md);    // momentum conservation:    float ermc = eng * revLight / mass;    vel->x += ermc * rand_vect.x;     vel->y += ermc * rand_vect.y;    vel->z += ermc * rand_vect.z;    float v12 = float3_sqr(*vel);    // energy conservation: old kinetic energy + photon energy = new kinetic energy + 'internal' energy    *int_eng += eng + 0.5f * mass * (v02 - v12);}

функция rand_uvect берет единичный вектор случайного направления, revLight константа, обратная скорость света во в единицах программы. Функция излучения фотона принимает на вход указатель на скорость атома, указатель на тепловую энергию атома / энтальпию, массу атома, переменную для генератора случайных чисел и указатель на главную структуру:

__device__ void radiate_photon(float3* vel, float* int_eng, float mass, uint4& rand_var, cudaMD* md){    float u0 = *int_eng;    float v02 = float3_sqr(*vel);   // square of initial velocity    float v0 = sqrt(v02);           // module of initial velocity    float ph_eng = 0.9f * u0;    float ermc = ph_eng * revLight / mass;    // random radiation    float3 rand_vect;    // random cosine between   (1-2a/v0) .. -1 with mean = -a/v0, where a = e/mc    float ermcv0 = ermc / v0;       // ph_eng/(m*c*v0)    float cos_phi;    if (ermcv0 >= 1.f)        cos_phi = -1.f;    else    {        int rnd = rnd_xor128(rand_var) % 2048;        float cos_phi = (float)rnd / 1024.f * (1.f - ermcv0); // 1024, because factor of 2          cos_phi -= 1.f;        rnd = rnd_xor128(rand_var) % 2048;        float theta = (float)rnd / 1024.f * numPi;  // 0 .. 2PI        rand_vect = get_angled_vector(*vel, cos_phi, theta);    }    // momentum conservation:    vel->x += ermc * rand_vect.x;       vel->y += ermc * rand_vect.y;    vel->z += ermc * rand_vect.z;    float v12 = float3_sqr(*vel);    // energy conservation:    *int_eng -= (ph_eng + 0.5f * mass * (v12 - v02));}

Тут учитывается, что выражение (7) может дать и косинус меньше -1, что нужно пофиксить. Поскольку в 3-мерном пространстве векторов, находящихся под заданным углом к данному, целый конус, функция get_angled_vector вычисляет один из них, по выбранному углу вращения:

__device__ float3 get_angled_vector(float3 invec, float cos_phi, float theta){    float3 v1, v2, v3;    float leng1 = float3_length(invec);    // only for non-zero input vectors!    v1 = make_float3(invec.x / leng1, invec.y / leng1, invec.z / leng1);    // find v2 perpendicular to v1    if (v1.x != 0.f)    {        v2.y = 1.f; v2.z = 1.f; // any coordinates        v2.x = -(v1.y * v2.y + v1.z * v2.z) / v1.x;    }    else        if (v1.y != 0.f)        {            v2.x = 1.f; v2.z = 1.f; // any coordinates            v2.y = -(v1.z * v2.z) / v1.y;     // a1 = 0 !        }        else // a1=0, b1=0, c1 <> 0:        {            v2.x = 1.f; v2.y = 0.f; v2.z = 0.f;        }    // v3 is perpendicular to both v1 and v2    v3.x = v1.y * v2.z - v1.z * v2.y;    v3.y = -v1.x * v2.z + v1.z * v2.x;    v3.z = v1.x * v2.y - v1.y * v2.x;    // convert v2 and v3 into unit vectors    float leng2 = float3_length(v2);    float leng3 = float3_length(v3);    v2.x /= leng2;    v2.y /= leng2;    v2.z /= leng2;    v3.x /= leng3;    v3.y /= leng3;    v3.z /= leng3;    float sinPhi, sinTh, cosTh;    sinPhi = sqrt(1 - cos_phi * cos_phi);    sincos(theta, &sinTh, &cosTh);    // parameteric equation of circle    float x = v1.x * cos_phi + sinPhi * (cosTh * v2.x + sinTh * v3.x);    float y = v1.y * cos_phi + sinPhi * (cosTh * v2.y + sinTh * v3.y);    float z = v1.z * cos_phi + sinPhi * (cosTh * v2.z + sinTh * v3.z);    float3 res = make_float3(x, y, z);    return res;} 

Генератор случайных чисел на CUDA не тривиальная задача, спасибо пользователю @denglide за обзор генераторов и пользователю @maratyszcza за xorShift генератор, который и был использован в данной работе.

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

Результаты

Ну а теперь к численным экспериментам. В качестве модельного объекта я взял 40000 атомов аргона, взаимодействующих согласно парному потенциалу из этой работы. Размеры бокса выбраны таким образом, чтобы плотность системы соответствовали плотности аргона при стандартных условиях (T = 298 К и P = 1 атм). Все приведенные результаты даны для 1 млн шагов интегрирования с шагом интегрирования 1 фс. Вычисления проводил на видеокарте GeForce RTX 2080 Ti.

Рис. 1. Зависимость кинетической энергии от времени эволюции системы

Ещё раз напоминаю, что целевую температуру я задаю через набор энергий фотонов, излучаемых термостатом. На данном этапе я следил за изменением кинетической энергии системы (рис. 1) при различных температурах термостата и разных вариантов начальных скоростей атомов. В ряде случаев начальные скорости были нулевыми (Ekin = 0), в других случаях были для всех атомов одинаковыми по модулю, но направленными хаотично (Ekin > 0). Независимо от величины начальных скоростей, кинетическая энергия выходила на примерно постоянное значение, зависящее от температуры термостата, т.е. такой термостат умеет разгонять или замедлять атомы, в зависимости от текущего состояния системы и целевой температуры, что в первую очередь от него и требуется.

Рис. 2. Распределение атомов по скоростям после 1млн шагов эволюции системы.

Второе, на что я обратил внимание, это распределение скоростей атомов, после 1 млн шагов, рис. 2. Несмотря на то, что начальные скорости были заданы одинаковыми по модулю, в результате работы термостата распределение атомов по скоростям приобрело вид, очень похожий на распределение Максвелла. Точно также как и в распределении Максвелла, меньшая температура соответствует более узкому пику. То, что такое распределение получилось само собой, при выполнении простых правил, я считаю определенным успехом.

Заключение

В качестве заключения кратко суммирую плюсы и минусы полученного способа задания температуры.

Плюсы предложенного термостата:

  1. На атомы воздействуют случайные импульсы, так что тепловое движение приобретает хаотический характер, в соответствии с общими представлениями.

  2. Распределение скоростей частиц в системе приобретает вид распределения Максвелла.

  3. Вычислительная сложность порядка O(N), как у других термостатов.

  4. Температура определяет скорости частиц, но не сводится к ним.

  5. Есть возможность учитывать индивидуальные свойства разных веществ, используя их температурную зависимость энтальпии (вернее энтальпийную зависимость температуры).

Минусы и недоработки:

  1. Сумма скоростей системы не равна нулю, возможно проблема в генераторе случайных чисел.

  2. Не удалось воспроизвести правильное давление (и его температурную зависимость), возможно решается подбором коэффициента между энтальпией и энергией фотона.

  3. Формула (7) выведена исключительно из математических соображений и имеет слабое физическое обоснование.

  4. Непонятно, откуда берется начальное распределение фотонов по энергиям, здесь оно просто постулируется на основании формулы Планка, а должно получатся "естественным образом".

Ну и напоследок, вот как выглядят 40000 атомов:

Подробнее..

Категории

Последние комментарии

  • Имя: Макс
    24.08.2022 | 11:28
    Я разраб в IT компании, работаю на арбитражную команду. Мы работаем с приламы и сайтами, при работе замечаются постоянные баны и лаги. Пацаны посоветовали сервис по анализу исходного кода,https://app Подробнее..
  • Имя: 9055410337
    20.08.2022 | 17:41
    поможем пишите в телеграм Подробнее..
  • Имя: sabbat
    17.08.2022 | 20:42
    Охренеть.. это просто шикарная статья, феноменально круто. Большое спасибо за разбор! Надеюсь как-нибудь с тобой связаться для обсуждений чего-либо) Подробнее..
  • Имя: Мария
    09.08.2022 | 14:44
    Добрый день. Если обладаете такой информацией, то подскажите, пожалуйста, где можно найти много-много материала по Yggdrasil и его уязвимостях для написания диплома? Благодарю. Подробнее..
© 2006-2024, personeltest.ru