Программно-инструментальные средства численного моделирования процессов адсорбции, хранения и десорбции водорода статическими и динамическими наноструктурами тема диссертации и автореферата по ВАК РФ 05.13.18, кандидат технических наук Суетин, Михаил Валерьевич

• Суетин, Михаил Валерьевич
• кандидат технических наук
• 2007

Оглавление диссертации кандидат технических наук Суетин, Михаил Валерьевич

ОБОЗНАЧЕНИЯ.

ВВЕДЕНИЕ.

ГЛАВА 1. МАТЕМАТИЧЕСКАЯ ПОСТАНОВКА ЗАДАЧИ МОДЕЛИРОВАНИЯ АДСОРБЦИИ, ХРАНЕНИЯ И ДЕСОРБЦИИ ВОДОРОДА НАНОСТРУКТУРАМИ.

1.1 Моделирование адсорбции, хранения и десорбции водорода наноструктурами методом молекулярной динамики.

1.1.1 Метод молекулярной динамики.

1.1.2 Периодические граничные условия.

1.2 Алгоритмы численного расчёта.

1.2.1 Алгоритм предиктор-корректор.

1.2.2 Предиктор - корректор алгоритм Нордсика.

1.2.3 Алгоритм Верле.

1.2.4 Скоростной алгоритм Верле.

1.2.5 Алгоритм Верле «прыжков лягушки» (Leapfrog).

1.3 Вычисление термодинамических параметров моделируемой системы.

1.4 Программный комплекс.

1.4.1 Блок подготовки.

1.4.2 Вычислительные модули.

1.4.3 Блок согласования.

1.4.4 Блок обработки результатов.

1.5 Тестовые расчёты.

1.5.1 Моделирования адсорбции водорода в нанотрубку с параметрами (9,9).

1.5.2 Сравнение результатов моделирования структуры фуллерена и фуллерита с опытными данными.

1.5.3 Устойчивость, сходимость и точность.

Выводы по главе 1.

ГЛАВА 2. МОДЕЛИРОВАНИЕ АДСОРБЦИИ, ХРАНЕНИЯ И ДЕСОРБЦИИ ВОДОРОДА СТАТИЧЕСКИМИ НАНОСТРУКТУРАМИ: ФУЛЛЕРЕНАМИ И ДЕФЕКТНЫМИ ФУЛЛЕРЕНАМИ.

2.1 Адсорбирование молекул водорода фуллеренами.

2.1.1 Постановка задачи.

2.1.2 Результаты моделирования.

2.2 Адсорбирование молекул водорода дефектными фуллеренами.

2.2.1 Постановка задачи.

2.2.2 Результаты расчётов.

2.3 Адсорбирование молекул водорода в кластеры, образованные фуллеренами.

2.3.1 Постановка задачи.

2.3.2 Результаты моделирования

Выводы по главе 2.

ГЛАВА 3. МОДЕЛИРОВАНИЕ АДСОРБЦИИ, ХРАНЕНИЯ И ДЕСОРБЦИИ ВОДОРОДА ДИНАМИЧЕСКИМИ НАНОСТРУКТУРАМИ: НАНОКАПСУЛАМИ.

3.1 Постановка задачи моделирования адсорбции, хранения и десорбции водорода нанокапсулами при изменении термодинамических условий.

3.1.1 Результаты моделирования адсорбции, хранения и десорбции водорода в нанокапсуле, состоящей из нанотрубки (10.10) и фуллерена Сбо.

3.1.2 Результаты моделирования адсорбции, хранения и десорбции водорода в нанокапсуле, состоящей из нанотрубки (10.10) и двух фуллеренов Сбо.

3.1.3 Результаты моделирования адсорбции, хранения и десорбции водорода в нанокапсуле, состоящей из нанотрубки (15.15) и фуллерена С240.

3.2 Постановка задачи моделирования адсорбции, хранения и десорбции водорода нанокапсулами при изменении электрического поля

3.2.1 Результаты моделирования адсорбции, хранения и десорбции водорода в нанокапсуле, состоящей из нанотрубки (10.10) и фуллерена Сбо5+.

3.2.2 Результаты моделирования адсорбции, хранения и десорбции водорода в нанокапсуле, состоящей из нанотрубки (10.10) и фуллерена Cso5+.

Выводы по главе 3.

Введение диссертации (часть автореферата) на тему «Программно-инструментальные средства численного моделирования процессов адсорбции, хранения и десорбции водорода статическими и динамическими наноструктурами»

Актуальность темы.

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

Для хранения водорода используются четыре основные технологии, позволяющие хранить водород в виде сжатого газа, жидкости, металлогидридов и при помощи адсорбции водорода углеродными наноматериалами, [73,58,101].

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

• объёмная плотность водорода - р = , W где тн - масса атомов водорода, W - объём, занимаемый водородом;

• массовая плотность водорода - Wt = ——--100%, тс+тн где тс - масса атомов углерода; Е

• энергетическая объёмная плотность - Ew = —, где Е - энергия, содержащаяся в атомах водорода Е

• энергетическая массовая плотность - Ет =н тн

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

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

Плотность водорода при давлении 0,1 МПа и температуре 273 К равна о

0,089 кг/м , что обуславливает создание систем хранения, в которых водород имеет более высокую плотность. На рис. 2. представлена диаграмма, характеризующая энергетическую объёмную плотность водорода в различных формах хранения, и другие виды топлива. Здесь, в отличие от рис. 1., показатели водорода резко отстают от традиционных видов топлива. Наилучшую энергетическую объёмную плотность среди форм хранения о водорода показали металлогидриды (10,78 ГДж/м ), затем жидкий водород

•5 о

9 ГДж/м ), газообразный водород при 70 МПа (4,7 ГДж/м'), и водород в наноструктурах (3,5 ГДж/м ). дизель бензин этан метан жидкий газообразный водород в водород в водород водород (70 металлах наноструктурах МПа) Л

Рис. 2. Энергетическая объёмная плотность (ГДж/м ) традиционных видов топлива и водорода в различных способах хранения.

Сравнивая четыре приведённых способа хранения по объёмной плотности водорода внутри систем хранения, (рис. 3.) видим, что металлогидриды занимают лидирующие позиции (150,8 кг/м ), даже превосходя жидкий водород (70,8 кг/м3), далее следует газообразный водород под давлением 70 МПа (36

3 3 кг/м ) и водород в наноструктурах (30 кг/м ).

На рис. 4. приведены показатели массовой плотности водорода (%) в зависимости от существующих способов хранения. Хранение водорода в сжатом состоянии при давлении 70 МПа является наилучшим с точки зрения гравиметрических показателей (10 масс.%). Далее следует хранение водорода в металлах (4 5 масс.%) при температуре 298 К и давлении 0,1 МПа, затем в сжиженном состоянии (1,6 масс.%), система хранения находится при температуре 298 К и давлении 0,1 МПа и в наноструктурах (0,01 масс.%) при тех же термодинамических параметрах. жидкий водород газообразный водород водород в металлах водород в

70 МПа) наноструктурах

Рис. 3. Зависимость объёмной плотности водорода от типа хранения жидкий водород газообразный водород водород в металлах водород в

70 МПа) наноструктурах

Рис. 4. Зависимость гравиметрической плотности водорода от типа хранения

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

Рассмотрим более подробно различные способы хранения водорода.

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

Проблема прочности баллонов очень актуальна, т.к. их разгерметизация > освобождает огромную энергию сжатого газа, которая может привести к значительным разрушениям. Предполагается, что баллоны будущего будут иметь три слоя: внутренний, состоящий из полимеров, средний, состоящий из углеродного волокна, обеспечивающего хорошие деформационные характеристики всей конструкции, и внешний слой, обеспечивающий защиту от механических повреждений и коррозии. В отличие от сжижения, процесс сжатия водорода требует гораздо меньшее количество энергии. Например, сжатие водорода от атмосферного давления до 45 МПа требует 25кДж/моль, что составляет 9% от наибольшей тепловой величины водорода (285 кДж/моль).

Жидкий водород хранится в криогенных системах при температуре 21,2К и атмосферном давлении. Вследствие низкой критической температуры водорода (ЗЗК) жидкий водород может храниться только в открытых системах, т.к. водород не имеет жидкой фазы при температурах выше критической. Давление в закрытых системах при нормальной температуре достигает величины в 100 МПа. С жидким водородом связаны две значительные проблемы: для его сжижения требуются большие энергетические затраты, составляющие до 40% от общего энергетического содержания водорода, и при его хранении происходит вскипание, вследствие чего он должен стравливаться, это приводит к потерям водорода как топлива на 1.5 - 2% в день, что тем более не безопасно в замкнутых помещениях. Величина стравливаемого водорода из-за притока тепла пропорциональна отношению площади системы хранения к объёму. Очевидно, что наилучшее отношение будет у сферических систем. Однако, такие системы в изготовлении более трудны, чем цилиндрические. Массовой содержание водорода в таких системах достигает 1,6 масс.%. Существует также другая разновидность водорода - кашеобразный водород. Это смесь жидкого и твёрдого водорода, она рассматривается как форма хранения водорода на воздушных судах, но также может быть использована в автомобильном транспорте. Два преимущества очевидны: более высокая плотность (примерно на 15% выше чем у жидкого водорода), т.к. присутствуют твёрдые составляющие и меньшая склонность к испарению, т.к. поступающая энергия в первую очередь используется для перевода водорода из одного агрегатного состояния в другое, т.е. энергия прежде всего тратится на плавку твёрдого * водорода и уже затем на процесс испарения. В настоящее время кашеобразный водород используется в космических программах в качестве топлива.

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

Гидриды формируются из многих металлов: палладия (Рс1Но.б), редкоземельных металлов (REH2 и REH3) и магния (MgH2). Однако больше всего для хранения водорода подходит сплав LaNi5, который обладает такими свойствами, как быстрая и полностью обратимая сорбция с малым гистерезисом. Объёмная плотность водорода в LaNi5H6 при давлении 0,2 МПа такая же, как у газообразного водорода при давлении 180 МПа. La и Ni имеют высокую молекулярную массу, вследствие чего содержание водорода в LaNi5H6 менее 2 мачх.%. Более высокая массовая плотность водорода может быть достигнута при использовании лёгких элементов, таких, как кальций и магний. Гидрид MgH2 содержит 7.6 масс.% водорода, но его формирование из газообразного водорода и магния происходит чрезвычайно медленно.

Mg2FeH6 Mg2NiH4 FeTiH2 ZrV2H5.5 LaNi5H6 PdH0.6 Sn(AIH4)4 Ce(AIH4)3 Zr(AIH4)3 ln(AIH4)3 Ti(AIH4)4 ScAIH4 Ga(AIH4)3 Ti(AIH4)3 AgAIH4 Fe(AIH4)2 Mn(AIH4)2 Ca(AIH4)2 Cu(AIH4)2 Mg(AIH4)2 Na2LiAIH6 Be(AIH4)2 KAIH4 NaAIH4 LiAIH4

0 ' 1 2 3 4 5 6 7 8 9 10 11 12 13

Массовое содержание водорода (%)

Рис. 5. Массовое содержание водорода в различных металлогидридах

И | | ^ШЛДД^Д! 3,6

М^ДР^^ИЧН—»;?. ' : • "У У :: У У; . У/ . : : .;. ; : ;/ ; ;ч:уу , у :.ууУ| 9 g

-■'■-■ ' .^-"""ifey'j 2,2 д— y^vg;" -А, - а - ЯШ■--'■'■- 33ШВШШШЯЖ 6,7

-.,. - ■■■■■■. ■ ■ ■: 42

Г.:. У::.,.,, -Vft : ' ■ ■■ -j,- : ~У„У, У5,3

Для того, чтобы улучшить кинетику формирования гидрида MgH^, были проведены различные исследования, которые показали необходимость существования катализатора для ускорения процесса диссоциации. Подобным катализатором может быть никель. Такой подход заключается в использовании двух или более компонентов для нивелирования слабых сторон друг друга. Комплексный гидрид Mg2NiH4 может хранить в себе до 3.6 масс.% водорода. Другой гидрид Mg2FeH6 содержит 5.5 масс.% водорода. На рис. 5. показано массовое содержание водорода в зависимости от типа металлогидрида.

Недостатками металлогидридных систем являются резкое падение массового содержания водорода после прохождения нескольких десятков циклов зарядки - разрядки, медленное формирование гидрида с большим выделением тепла, которое необходимо отводить. В настоящее время всесторонне протестированы только системы Mg2FeH6, Mg2NiH4, LaNi5H6 и PdHo.6, не отличающиеся большими показателями содержания водорода.

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

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

Lee S. et al. [69], моделируя взаимодействия системы различных нанотрубок при температуре 77 К и давлении 15 МПа, получили массовую плотность водорода - Wt равную 6.88%. Результат YeY. et al. [100] для системы нанотрубок при температуре 80 К и давлении 10 МПа составляет

Wt=8.25 %. В [96] при температуре 77 К и давлении 5 МПа Wang Q. и Johnson К. в результате моделирования получили Wt=5 % для системы нанотрубок. Zuttel A. et al [101] опытным путём получили Wt=5,5 %. В работе [97] Wang Q. и Johnson К. выполнили моделирование пространственного расположения нанотрубок с целью получения наилучшей структуры для повышения количества хранимого водорода. Показано, что расстояние между нанотрубками в большой степени влияет на количество адсорбированного водорода, которое достигает максимального значения, когда влияние соседних нанотрубок мало, и адсорбцию водорода можно рассматривать как на одну изолированную нанотрубку. После оптимизации пространственной структуры системы нанотрубок при давлении 5МПа и температуре 77К Wt составило 10%. Эти данные соответствуют [96], где те же авторы, при таких же термодинамических параметрах для одиночной нанотрубки получили Wt=10.5%. •

Кроме того, необходимо отметить, что в [97] значение плотности водорода вблизи стенок нанотрубки превосходит плотность жидкого водорода в 6 раз, что является ошибочным вследствие пренебрежения размерами молекул водорода. Шаг распределения плотности был выбран менее 1 А, тогда как диаметр атома водорода составляет 2,4 А. В результате был неправильно построен профиль плотности водорода внутри нанотрубки.

В [88] Simonyan et al. осуществили моделирование адсорбции водорода на заряженные одностенные нанотрубки. Это дало увеличения показателя адсорбции Wt всего на несколько десятых долей процента даже при использовании улучшенной конфигурации расположения нанотрубок.

Следует учитывать, что на внутреннюю поверхность нанотрубок адсорбировалось только Wt=l,6% при давлении 5МПа и температуре 77К, а при давлении 5МПа и температуре 293К результат составил Wt=0,2%. Этот результат согласуется с опытными данными Anson A. et al. [38], где Wt=0,l% и результатами моделирования Guay P. et al. [53], где авторы получили результат адсорбции менее одного массового процента. Ими был получен интересный результат, что металлические примеси между нанотрубками значительно увеличивают адсорбцию водорода, при этом возрастает показатель объёмной адсорбции, но массовый показатель практически не увеличивается вследствие большой массы примесей.

Нечаев Ю.С. поднимает вопрос о существовании суперадсорбентов в связи с недостаточно полным пониманием в мире природы адсорбции водорода углеродными наноструктурами [19,20]. Следы суперадсорбции проявляются в статьях многих авторов из различных лабораторий. В 2006г. появилась работа Yang R. et al., вошедшая в отчёт американского департамента энергетики по водородной тематике, в которой описывается механизм адсорбции, способный достичь показателя Wt=4% при нормальной температуре и давлении 10 МПа [98]. Проведённые исследования Chen J. et al. по адсорбции водорода в неорганические нанотрубки прорывных результатов не принесли [43]. Массвое и объёмное содержание водорода почти равно приведённым раннее показателям адсорбции в углеродные нанотрубки. Численный эксперимент Jhi S. et al. [63] с участием неорганических нанотрубок также показал результаты, сравнимые с адсорбцией водорода в углеродные нанотрубки.

Покрытые одноатомным слоем никеля углеродные нанотрубки, как это показал Lee J. et al. [70] способны хранить при нормальных условиях (Т=300 К, Р=0,1 МПа) до Wt=10% водорода. При этом десорбция водорода происходит при температуре 328 К, которая является идеальной для использования в топливных элементах. Предприняты попытки аналитически вычислить предельную физическую адсорбцию водорода. Например, Богданов А.А. в [3] г выполнил расчёты, согласно которым двухсторонняя адсорбция водорода на графитовый лист при 77 К составляет Wt=5%, а при 293 К всего Wt=l% при рассмотрении диапазона давлений от 1 МПа до 10 МПа. На рис. 6. - диаграмма результатов работ различных авторов, показывающая Массовое содержание водорода в системах нанотрубок, одиночных нанотрубках и во внутреннем пространств;; нанотрубок.

Zuttel А. [100] (77К 5МПа) Wang Q. [96] (77К 5МПа) Ye V. [100] (80К ЮМПа) LeeS. [69] (77К 15МПа) Yang R. [98] (ЗООКЮМПа) Guay Р. [53] (298 ЗМПа) Anson А. [38] (3Г0К 2МПа) Wang Q. [96] (77К 5МПа) Wang Q. [97] (77К 5МПа) Lee J. [70] (298К ЮМПа) Wang Q. [96] (77К 5МПа) Wang Q. [97] (77K 5МПа)

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

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

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

Для одиночной нанотрубки /

10,5 • • - , I 10 :: ' :' : ' ■ ■ . . | 10

Для внутреннего пространства нанотрубки

1,6

4 6 8

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

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

Объектом исследования являются процессы адсорбции, хранения и десорбции молекул водорода статическими наноструктурами - фуллеренами и дефектными фуллеренами, а также динамическими наноструктурами -нанокапсулами; математическая модель адсорбции, хранения и десорбции водорода наноструктурами.

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

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

Для реализации поставленной цели необходимо решить следующие задачи:

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

- провести вычислительные эксперименты по моделированию адсорбции, хранению и десорбции водорода статическими наноструктурами: фуллеренами и дефектными фуллеренов при различных термодинамических параметрах; определить температуру и давление, при которых адсорбция водорода проявляется наиболее сильно;

- численно выявить динамику адсорбции, хранения и десорбции водорода фуллеренами и дефектными фуллеренами различного диаметра;

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

- изучить степень влияния термодинамических параметров на структуру нанокапсул; определить параметры, при которых происходит десорбция водорода;

- численно исследовать характеристики нанокапсул, состоящие из нанотрубки (10.10) и фуллерена с6о, нанотрубки (10.10) и двух фуллеренов с6о, нанотрубки (15.15) и фуллерена С24о> по адсорбции, хранению и десорбции водорода при изменении термодинамических параметров; определить адсорбционные характеристики нанокапсул;

- выявить величины напряжённости электрического поля, необходимой для перемещения заряженных фуллеренов Сбо5+ и Cgo5+ в нанотрубке (10.10);

- численно исследовать характеристики нанокапсул, состоящие из нанотрубки (10.10) и заряженного фуллерена с6о5+, нанотрубки (10.10) и заряженного фуллерена Cso5+ по адсорбции, хранению и десорбции водорода под действием электрического поля.

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

Программно-инструментальные средства реализованы с помощью сред визуального программирования Borland С++ Builder 6 и Borland Delphi 7.

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

Задача моделирования адсорбции методом молекулярной динамики исследована на устойчивость, сходимость и точность.

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

- математическая модель для решения задач адсорбции, хранения и десорбции водорода статическими и динамическими наноструктурами.

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

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

- временные пределы по адсорбции и десорбции водорода фуллеренами и дефектными фуллеренами различного диаметра;

- величины выявленных энергетических параметров нанокапсул, изменяющихся при смене фаз адсорбции, хранения и десорбции водорода;

- зависимость влияния термодинамических параметров на структуру нанокапсул; значения температуры и давления наиболее выгодные для адсорбции, хранения и десорбции водорода;

- адсорбционные характеристики нанокапсул, состоящие из нанотрубки

10.10) и фуллерена С6о, нанотрубки (10.10) и двух фуллеренов С6о> нанотрубки »

15.15) и фуллерена С24о;

- величины напряжённости электрического поля, необходимые для перемещения заряженных фуллеренов Сбо5+ и Cg05+ в нанотрубке (10.10);

- адсорбционные характеристики нанокапсул, состоящие из нанотрубки (10.10) и заряженного фуллерена С6о5+, нанотрубки (10.10) и заряженного фуллерена Cgo5+ по адсорбции, хранению и десорбции водорода под действием электрического поля;

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

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

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

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

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

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

- выявлена динамика закрытия нанотрубок с адсорбированным водородом фуллеренами, вычислено изменение плотности водорода в процессе закрытия; определена температура выхода фуллеренов из нанотрубок;

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

Практическая полезность исследования состоит в том, что оно связано с моделированием процессов адсорбции, хранения и десорбции водорода статическими и динамическими наноструктурами при помощи программного комплекса также осуществляющего интерпретацию полученных результатов. При помощи компьютерного моделирования производятся точные расчёты свойств новых материалов для хранения водорода. Ценность полученных решений заключается в получении результатов при моделировании адсорбции, хранения и десорбции молекул водорода отдельными наноструктурами и их образованиями, что является невозможным при опытном исследовании. Полученные результаты являются новыми и дают представление об адсорбции, хранении и десорбции водорода статическими и динамическими наноструктурами. Разработанные методики и программный комплекс позволяют моделировать процессы адсорбции, хранения и десорбции молекул водорода на поверхности наноструктур различных видов и детально исследовать, особенности процессов в зависимости от термодинамических параметров и влияния электрического поля. Работа выполнена в рамках госбюджетной научно-исследовательской темы государственный регистрационный номер 012006097787, осуществлённой Институтом прикладной механики УрО РАН.

Диссертационная работа выполнялась при поддержке комплексной программы фундаментальных исследований Президиума РАН № 26 «Водородная энергетика», гос. контракт № 10002-251/П-26/117-383/290404-138, грантов молодых учёных и аспирантов УрО РАН 2005-2006гг. и стипендии президента Удмуртской Республики 2006-2007гг.

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

Апробация работы Основные результаты диссертационной работы и отдельные её части были представлены и доложены на следующих российских и международных конференциях: Первый Всемирный Конгресс "Альтернативная энергетика и экология" - (2006), Вторая Российская конференция «Физические проблемы водородной энергетики» (Санкт-Петербург, - 2005), Третья Российская конференция «Физические проблемы водородной энергетики (Санкт-Петербург, - 2006), 14-я зимняя школа по механике сплошных сред (Пермь, - 2005), The International Conference on Carbon (Aberdeen, - 2006), Конференция "Математическое моделирование и краевые задачи" (Самара, - 2005), Третья Всероссийская научная молодёжная конференция «Под знаком Сигма» (Омск, -2005), 14-я всероссийская школа-конференция молодых ученых "Математическое моделирование в естественных науках" (Пермь, - 2005), Школа-семи.чар "Нанотехнологии и наноматериалы" - "КоМУ 2005" (Ижевск, -2005), Конференция "ICOC" (Москва, - 2005), Первый Российский научный форум - "Демидовские чтения" (Екатеринбург, - 2006), Конференция "Теория управления и математическое моделирование" (Ижевск, - 2006), Научная конференция молодых учёных - «Поздеевские чтения» (Пермь, - 2006), Конференция "Материалы и технологии XXI века" (Пенза, - 2006), IV Международная научно-техническая конференция "Материалы и технологии

XXI века" (Пермь, - 2006), VI конференция молодых учёных «КоМУ - 2006» (Ижевск, - 2006), Зимняя школа по механике сплошных сред (Пермь, - 2007), Nanoizh 2007 (Ижевск, - 2007).

Публикации. Основные научные результаты по теме диссертации опубликованы в 16 научных работах, из них 8 - статей, 8 материалов конференций и 3 научно-технических отчёта, получен патент на изобретение. Автор имеет 4 научных труда в изданиях, выпускаемых в РФ и рекомендуемых ВАКом для публикации основных результатов диссертации.

Структура и объём диссертации.

Диссертация состоит из введения, трёх глав, заключения и списка цитируемой литературы. Работа изложена на 131 странице, включая 79 рисунков. Список литературы содержит 102 наименования.

Заключение диссертации по теме «Математическое моделирование, численные методы и комплексы программ», Суетин, Михаил Валерьевич

Выводы по главе 3

1. Численно показано, что нанокапсулы сохраняют при нормальных условиях Wt=l,6% водорода, которое было адсорбировано при температуре кипения азота (77 К) и давлении 5-10 МПа.

2. Численно показано, что нанокапсул а, состоящая из нанотрубки (10,10) и фуллерена Сбо> способна хранить при нормальных условиях Wt=l,6%, что на два порядка больше, чем у нанотрубок с теми же параметрами.

3. Численно выявлена температура десорбции водорода из нанокапсулы, состоящейиз нанотрубки (10,10) и фуллерена Сбо> равная 550К при внешнем давлении 0,1 МПа.

4. Численно рассчитанные процессы закрытия-открытия дефекта нанокапсулы заряженным фуллереном Сбо5+ изменяющимся электрическим полем Е напряженностью 1.044*109 В/м показали слабое влияние изменения термодинамических параметров на нанокапсулу: при повышении температуры с 77 до 300 К и понижении давления с 5 до 0,1 МПа плотность водорода снизилась ha 9 масс.%.

ЗАКЛЮЧЕНИЕ

В ходе выполнения диссертационной работы были достигнуты следующие цели и результаты:

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

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

3. Выяснено влияние температуры и давления на адсорбцию водорода фуллеренами. При уменьшении температуры и повышении давления количество адсорбированного водорода возрастает, достигая при температуре 60 К и давлении 10 МПа достигая 11,98 масс.%.

4. Исследовано, что дефектные фуллерены С46, С]67 и С220 с микропорами насыщают своё внутреннее пространство водородом при малых давлениях (0,1-2 МПа) до 0,71 масс.%, 1,28 масс.%, 1,27 масс.% соответственно, а дефектные фуллерены с мезопорами С505 при высоких (более 10 МПа) до 2,04 масс.%

5. Исследовано, что нанокапсулы способны хранить при нормальных условиях 1,6 масс.% водорода, адсорбированных при температуре 77 К и давлении 5 МПа.

6. Проведено исследование процессов закрытия-открытия дефекта нанокапсулы заряженным фуллереном Сбо5+ изменяющимся электрическим полем Е напряженностью 1.044*109 В/м, которое показало существование слабого влияния изменения термодинамических параметров на нанокапсулу: при повышении температуры с 77 К до 300 К и понижении давления с 5 МПа до 0,1 МПа плотность водорода снизилась на 9 %.

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

Список литературы диссертационного исследования кандидат технических наук Суетин, Михаил Валерьевич, 2007 год

1. Армстронг Т.М., Хейр М.Д. Исследования по хранению водорода в Ок-Риджской национальной лаборатории // Альтернативная энергетика и экология. 2004. № 2. С. 15-20.

2. Берлин А.А., Балабаев Н.К. Имитация свойств твёрдых тел и жидкостей методами компьютерного моделирования // Соросовский образовательный журнал № 11. 1997. С. 85-92.

3. Богданов А.А. О предельной физической адсорбции водорода в углеродных материалах // ЖТФ. 2005. Т. 75, Вып. 9. С. 139-142.

4. Липанов A.M., Суетин М.В. Моделирование механизма для хранения водорода, состоящего из нанотрубки и фуллерена // Материалы Третьей Российской Конференции "Физические проблемы водородной энергетики". Санкт-Петербург, 2006. С. 16-17.

5. Вахрушев А.В., Липанов A.M., Суетин М.В. Моделирование процессов адсорбирования водорода наноструктурами // Альтернативная энергетика и экология. 2007. № 1. с. 13-20.

6. Суетин М.В. Численное исследование процессов поглощения, хранения и выделения водорода фуллеренами. Проблемы термогазодинамики и прочности механических систем. Ижевск.: ИПМ УрО РАН, 2005. С.157-167.

7. Вахрушев А.В., Липанов A.M., Суетин М.В. Исследование процессов поглощения, хранения и выделения водорода в нанообъектах. Отчёт о НИР к государственному контракту №10002-251/П-26/117-383/290404-138. Ижевск, 2004.-20 с.

8. Суетин М.В. Численное исследование процессов поглощения и выделения водорода фуллеренами // Материалы Зимней школы по механике сплошных сред (четырнадцатая). Пермь, 2005. С .60.

9. Вахрушев А.В., Суетин М.В. Моделирование процессов водородной адсорбции углеродными наноструктурами // Альтернативная энергетика и экология. 2006. № 6. С. 64-66.

10. Грег С., Синг К. Адсорбция, удельная поверхность, пористость: Пер. с англ. 2-е изд.-М.: Мир, 1984.

11. Гулд X., Тобочник Я. Компьютерное моделирование в физике: Часть 1.С. 14-16.

12. Зайцев Общая химия. М. Высшая школа, 1983. - 248 с.

13. Зефиров Ю.В., Зоркий П.М. Журн. структур, химии. 1974. Т. 15, № 1. С. 118-122.

14. Кобелев Н.П., Николаев Р.К., Сидоров Н.С., Сойфер Я.М. Температурная зависимость упругих модулей твёрдого Сбо // Физика твёрдого тела. 2001. Т. 43, Вып. 12. С. 2244-2250.

15. Коротеев А.С., Миронов В.В., Смоляров В. А. Перспективы использования водорода в транспортных средствах // Альтернативная энергетика и экология. №1. 2004. С. 5-13.

16. Липанов A.M., Вахрушев А.В., Суетин М.В., Вахрушев А.А. Моделирование процессов поглощения и выделения водорода фуллеренами и фуллеритами // Химическая физика и мезоскопия. 2005. Т. 7,№3. С. 309-318.

17. Суетин М.В. Моделирование процессов адсорбции и десорбции молекул водорода фуллеренами и их кластерами // Материалы Второй Российской Конференции "Физические проблемы водородной энергетики". Санкт-Петербург, 2005. С. 15.

18. Морозов А.И. Физика твердого тела кристаллическая структура Учебное пособие. М.: 2005 34 - с.

19. Нечаев Ю.С., Алексеева O.K., Гусев А.Л., Филиппов Г.А. «Открытые» вопросы о природе и характеристиках сорбции водорода углеродными наноматериалами и пути их решения // Альтернативная энергетика и экология. 2006. № 4. С. 15-18.

20. Нечаев Ю.С. О природе, кинетике и предельных значениях сорбции водорода углеродными наноструктурами // УФН. Т. 176, № 6. С. 581-610.

21. Никитина Е.А. Рентгенографическое исследование и компьютерное моделирование углеродных материалов Петрозаводск, 2003. 22 с.

22. Романова Т.А., Краснов П.О., Качин С.В., Аврамов П.В. Теория и практика компьютерного моделирования нанообъектов, Справочное пособие.: Красноярск ИПЦ КГТУ, 2002.

23. Сидоров Л.Н., Юровская М.А. и др. Фуллерены: Учебное пособие. -М.: Экзамен, 2005. 688 с.

24. Суетин М.В. Моделирование процессов адсорбирования водорода фуллеренами // Сборник статей IV Международной научно-технической конференции "Материалы и технологии XXI века". Пенза, 2006. С. 241245.

25. Суетин М.В. Численное исследование процессов поглощения, хранения и выделения водорода фуллеренами // Труды Второй Всероссийской научной конференции "Математическое моделирование и краевые задачи". Самара, 2005. С. 234-237.

26. Суетин М.В. Исследование процессов поглощения, хранения и выделения водорода фуллеренами // Материалы Третьей Всероссийской научной молодёжной конференции «Под знаком Сигма». Омск, 2005. С. 245.

27. Суетин М.В., Вахрушев А.В. Исследование процессов поглощения, хранения и выделения водорода фуллеренами // Труды 14-ой Всероссийской школы-конференции молодых ученых "Математическое моделирование в естественных науках". Пермь, 2005. С. 68

28. Суетин М.В. Численное исследование процессов поглощения, хранения и выделения водорода фуллеренами // Труды Второй

29. Всероссийской научной конференции "Математическое моделирование и краевые задачи". Самара, 2005. С. 234-237.

30. Суетин М.В. Хранение водорода в фуллеренах, их модификациях и полученных из них образований // Материалы Школы-семинара "Нанотехнологии и наноматериалы КоМУ 2005". Ижевск, 2005. С. 57.

31. Суетин М.В. Хранение водорода в фуллеренах, их модификациях и составленных из них кластеров // Материалы Первого Российского Научного Форума "Демидовские Чтения". Екатеринбург, 2006. С. 205206.

32. Суетин М.В. Моделирование адсорбирования водорода фуллеренами // Сборник научных трудов Научной конференции молодых учёных -"Поздеевские чтения". Пермь, 2006. С. 125-127.

33. Суетин М.В. Адсорбция водорода углеродными кластерами на основе фуллеренов // Материалы VI конференции молодых учёных "КоМУ -2006". Ижевск, 2006. С. 56.

34. Харрис П. Углеродные нанотрубы и родственные структуры. Новые материалы XX века. М.: Техносфера, 2003.-336 с.

35. Хеерман Д.В. Методы компьютерного эксперимента в статистической физике. М.: Наука, 1990. - 176 с.

36. Эткинс П. Физическая химия том 1: Пер. с англ. М.: Мир, 1980. -584 с.

37. Эткинс П. Физическая химия том 2: Пер. с англ. М.: Мир, 1980. -584 с. '

38. Allen М. P. and Tildesley А. К. Computer Simulation of Liquids. -Oxford: Clarendon Press. 1987.

39. Bischof C., Bucker M. Computing derivatives of computer programs // Modern Methods and Algorithms of Quantum Chemistry, J. Grotendorst (Ed.), John von Neumann Institute for Computing, Julich, NIC Series, 2000. Vol. l.P;287-299.

40. Brown D., Clarke J., Okuda M., Yamazaki Т., A domain decomposition parallel processing algorithm for molecular dynamics simulations of polymers. // Сотр. Phys. Comm. 1994. Vol. 83 P. 1-13.

41. Buczek В., Czepirski L., Zietkiewicz J. Improvement of Hydrogen Storage Capacity for Active Carbon // Adsorption. 2005. Vol. 11. P. 877-880.

42. Cabria I., Lopez M., Alonso J. Adsorption of hydrogen on normal and pentaheptite single wall carbon nanotubes // Eur. Phys. J. 2005. D 34. P. 279282.

43. Cagin Т., Che J., Qi Y., Zhou Y., Demiralp E., Gao G., Goddard III W. Computational materials chemistry at the nanoscale // Journal of Nanoparticle Research. 1999. N 1. P. 51-69.

44. Chen J. Wu F. Review of hydrogen storage in inorganic fullerene-like nanotubes // Applied Physics A: Materials Science & Processing. 2004 P. 989-994.

45. Derylo-Marczewska A., Goworek J., Swiaztkowski A., Buczek B. Influence of differences in porous structure within granules of activated carbon on adsorption of aromatics from aqueous solutions // Carbon. 2004. N 42. P. 301-306.

46. Diaz-Tendero S. Structure and fragmentation of neutral and positively charged carbon clusters and fullerenes. Madrid. 2005. 240 p.

47. Fischer D. Dissertation: Theoretical investigation of nanoscale solid state and cluster structures on surfaces. 144 p.

48. Fonseca A., Pierard N., Tollis S., Bister G., Konya Z., Nagaraju N., Nagy J.B. Hydrogen storage in carbon nanotubes produced by CVD // J. Phys. IV France. 2002. Vol. 12. P. 4-129.

49. Furuya Y., Hashishin Т., Iwanaga H., Motojima S., Hishikawa Y. Interaction of hydrogen with carbon coils at low temperature // Carbon. 2004. N42. P. 331-335.

50. Garberoglio G., Skoulidas A., Johnson K. Adsorption of Gases in Metal Organic Materials: Comparison of Simulations and Experiments // J. Phys. Chem. B. 2005. Vol. 109. P. 13094-13103.

51. Gauss J. Molecular Properties // Modern Methods and Algorithms of Quantum Chemistry, Proceedings, Second Edition, J. Grotendorst (Ed.), John von Neumann Institute for Computing, Julich, NIC Series. 2000. Vol. 3. P. 541-592.

52. Gerndt M. Parallel programming models, tools and performance analysis // Modern Methods and Algorithms of Quantum Chemistry, J. Grotendorst (Ed.), John von Neumann Institute for Computing, Julich, NIC Series. 2000. Vol. 1. P. 9-27.

53. Guay P., Stansfield В., Rochefort A., On the control of carbon nanostructures for hydrogen storage applications // Carbon. 2004. N 42. P. 2187-2193.

54. Gueorguiev G. K., Pacheco J.M. Tomanek D. Quantum size effects in the polarizability of carbon fullerenes. // Physical Review Letters. 2004. Vol. 92, N. 21, P. 215501-1 -215501-4.

55. Gulseren O., Yildirim Т., Ciraci S. Effects of hydrogen adsorption on single-wall carbon nanotubes: Metallic hydrogen decoration // Physical Review B. 2002. Vol. 66. P. 121401-1 121401-4.

56. Han S., Kim H., Han K., Lee J., Lee H., KangaSeong J. Nanopores of carbon nanotubes as practical hydrogen storage media // Applied physics letters. 2005. Vol. 87. P. 213113-1 213113-3.

57. Hassel M. Thesis for the degree of doctor of philosophy: Vibrations and Rotations of Diatomic Molecules at Surfaces. 62 p.

58. Hoover W. G. Isomorphism linking smooth particles and embedded atoms // Physica A. 1998. V. 260. .3-4. P. 44-254.60. Humphrey, Dalke, Schulten, VMD Visual Molecular Dynamics // J.

59. Molec. Graphics 1996, Vol. 14, P. 33-38.

60. Iijima S. Helical microtubules of graphitic carbon // Nature. 1991. Vol. 354. N6348. P. 56-58.

61. Iijima S., Ichihashi T. Single shell carbon nanotubes of 1-nm diameter // Nature. 1993. Vol. 363. N 6340. P. 603-605.

62. Jhi S., Kwon Y. Hydrogen adsorption on boron nitride nanotubes: A path to room-temperature hydrogen storage // Physical Review B. 2004. Vol. 69. P. 245407-1* 245407-4.

63. Kim, Choi, Chang, Tomanek Defective fullerenes and nanotubes as molecular magnets: An ab initio study // Physical Review B. 2003. Vol. 68. P. 125420-1 125420-4.

64. Krivtsov A. I. Second Order Equation of State for Lennard-Jones Chain // Proceedings of the XXVIII Summer School Actual Problems in Mechanics. St.-Petersburg. Russia. 2001. V. 1. P. 79-90.

65. Kroto H.W., Heath J.R., Curl R.F., Smalley R.E. C60: Buckminsterfullerene //Nature. 1985. Vol. 318. P. 162.

66. Lau K., Hiu D. Effectiveness of using carbon nanotubes as nanoreinforcement for advanced composite structures // Carbon. 2002. Vol. 40. N. 9. P. 1605-1606

67. Lee S., Kay A., Lee Y., Seifert G., Frauenheim T. Novel Mechanism of Hydrogen Storage in Carbon Nanotubes // Journal of the Korean Physical Society. ?001. Vol. 38. N 6. P. 686-691.

68. Lee J., Kim H., Kang J. Hydrogen storage and desorption properties of Ni-dispersed carbon nanotubes // Appl. Phys. Lett. 2006. Vol. 88. P. 143126-1 143126 -3.

69. Levesque D., Gicquel A., Darkrim F., Kayiran B. Monte Carlo simulations of hydrogen storage in carbon nanotubes // J. Phys.: Condens. Matter. 2002 Vol. 14. P. 9285-9293.

70. Li J., Furuta Т., Goto H., Ohashi Т., Fujiwara Y., Yip S. Theoretical evaluation of hydrogen storage capacity in pure carbon nanostructures // Journal о!'chemical physics. Vol. 119, N 4 P. 2376- 2385.

71. Maeland A.J., Hydrogen storage for automobile transport survey and approbation. - Institute of Energy Technology, 2000. - 45 p.

72. Maruyama S., Kimura T. Molecular dynamics simulation of hydrogen storage in single-walled carbon nanotubes // 2000 ASME International Mechanical Engineering Congress and Exhibit. Orland. 2000.

73. Marx D., Hutter J. Ab initio molecular dynamics: Theory and Implementation // Modern methods and algorithms of quantum chemistry, J. Grotendorst (Ed.), John von Neumann Institute for Computing, Julich, NIC Series. 2000. Vol. 1. P. 301-449.

74. Metropolis N., Rosenbluth A., Rosenbluth M., Teller A., Teller E., Equation of state calculations by fast computing machines. // J. Chem. Phys. 1953. Vol. 21. P. 1087-1092.

75. Nechaev Y.S. New Challenges in Superconductivity: Experimental Advances and Emerging Theories // The NATO Science Series II — Mathematics, Physics and Chemistry. 2005. Vol. 183 P. 91

76. Neskovicl O., Djustebek J., Djordjevic V., Cveticanin J., Velickovic S., Veljkovic M., Bibic N. Hydrogen storage on activated carbon nanotubes // Digest Journal of Nanomaterials and Biostructures. 2006. Vol. 1. No. 4. P. 121 127.

77. Nose S. Constant-Temperature Molecular-Dynamics // Journal of Physics -Condensed Matter. 1990. V. 2. P. SA115-SA119.

78. Pol V., Motiei M., Gedanken A., Calderon-Moreno J., Yoshimura M. Carbon spherules: synthesis, properties and mechanistic elucidation. // Carbon. 2004. N42. P. 111-116.

79. Qian t>., Liu W., Ruoff R. Mechanics of C60 in Nanotubes // J. Phys. Chem. B. 2001. Vol. 105. P. 10753-10758.

80. Ramachandran G.N., Sasisekharan V., Adv. Prot. Chem. 1968. Vol. 23. P 283.

81. Patent 6,473,351 US. Nanocapsules containing charged particles, their uses and methods of forming same / Tomanek et al. (2002)

82. Schimmel H.G., Nijkamp G., Kearley G.J., Rivera A., de Jong K.P., Mulder F.M. Hydrogen adsorption in carbon nanostructures compared // Materials Science and Engineering B. 2004. Vol. 108. P. 124-129.

83. Scott R. A., Scheraga H. A. // J. Chem. Phys. 1965. Vol. 45. P. 2091.

84. Shelimov K.B., Jarrold M.F. Carbon Clusters Containing Two Metal Atoms: Structures, Growth Mechanism, and Fullerene Formation // J. Am. Chem. Soc. 1996. Vol. 118. P. 1139-1147.

85. Simonyan, Diep, Johnson Molecular simulation of hydrogen adsorption in charged single-walled carbon nanotubes // Journal of Phys. Chem. B. 1999. N 111. P. 9778-9783.

86. Srivastava D., Menon M., Cho K. Computational nanotechnology with carbon nanotubes and fullerenes // Computing in science & engineering American Institute of Physics. 2001. P. 42-55.

87. Svensson K., Bengtsson L., Bellman J., Hassel M., Persson M., Andersson S. Two-Dimensional Quantum Rotation of Adsorbed H2 // Physical Review Letters. 1999. Vol. 83, N 1. P. 124-127.

88. Sumanasekera G.U., Adu C.K., Pradhan В., Chen G., Romero H., Eklund P.C. Thermoelectric study of hydrogen storage in carbon nanotubes // Physical Review B. 2001. Vol. 65. P. 035408-1 035408-5.

89. Vakhrouchev A.V., Suyetin M.V. Molecular simulation of hydrogen adsorption on fullerenes and in fullerites // Abstract "Carbon 2006" Aberdeen, 2006. P. 108.

90. Vakhrouchev A.V. Simulation of nano-elements interactions and self-assembling // Modeling and simulation in materials science and engineering 2006. N 14. P. 975-991

91. Verlet L. Computer "experiments" on classical fluids. I. Thermo dynamical properties of Lennard-Jones molecules // Phys. Rev. 1967. Vol. 159, N 1. P. 98-103.

92. Verlet L. Computer 'experiments' on classical fluids. I. Thermodynamical properties of Lennard-Jones molecules // Phys. Rev. 1967. V. 159. .98. P. 103.

93. Wang Q., Johnson K. Molecular simulation of hydrogen adsorption in single-walled carbon nanotubes and idealized carbon slit pores // Journal of chemical physics. 1999. Vol. 110, N 11. P. 577-586.

94. Wang Q., Johnson K., Optimization of Carbon Nanotube Arrays for Hydrogen Adsorption//Journal Phys. Chem. B. 1999. N 103. P. 4809-4813.

95. Yang R., Li Y., Lachawiec A. Hydrogen Storage in Graphite Nanofibers and the Spillover Mechanism // 2006 Annual progress report U.S. D.O.E. Hydrogen program. 2006. P. 484-485

96. Yoon, Berber, Tomanek. Energetics and packing of fullerenes in nanotube peapods//Physical Review B. 2005. Vol. 71. P. 155406-1 155406-4.

97. Ye Y. et al. Hydrogen adsorption and cohesive energy of single-walled carbon nanotubes // Applied physics letters. 1999. Vol. 74, N 16. P. 23072309.

98. Zhu H.W., Chen A., Mao Z.Q., Xu C.L., Xiao X., Wei B.Q., Liang J. Wu D.H. The effect of surface treatments on hydrogen storage of carbon nanotubes // Journal of Materials Science Letters 2000 Vol. 19. N. 14 pp. 1237-1239.t