WWW.KONFERENCIYA.SELUK.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА - Конференции, лекции

 

Pages:   || 2 |

«Посвящается 100-летию со дня рождения профессора Лебедева Ивана Кирилловича ТЕПЛОФИЗИЧЕСКИЕ ОСНОВЫ ЭНЕРГЕТИЧЕСКИХ ТЕХНОЛОГИЙ Сборник научных трудов II Всероссийской научно-практической ...»

-- [ Страница 1 ] --

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ

Федеральное государственное бюджетное образовательное учреждение

высшего профессионального образования

«НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ

ТОМСКИЙ ПОЛИТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ»

ЭНЕРГЕТИЧЕСКИЙ ИНСТИТУТ

Посвящается 100-летию со дня рождения

профессора Лебедева Ивана Кирилловича

ТЕПЛОФИЗИЧЕСКИЕ ОСНОВЫ

ЭНЕРГЕТИЧЕСКИХ ТЕХНОЛОГИЙ

Сборник научных трудов II Всероссийской научно-практической конференции с международным участием 06 – 08 октября 2011 г.

Томск УДК 621.1(063) ББК 31.3л Т Теплофизические основы энергетических технологий:

сборник научных трудов II Всероссийской научно-практической Т конференции / под ред. Г.В. Кузнецова, А.С. Заворина, К.В. Бувакова; Томский политехнический университет. – Томск:

Изд-во Томского политехнического университета, 2011. – 358 с.

ISBN 978-5-98298-959- Сборник содержит теоретические и практические работы, отражает аспекты консолидации возможностей и усилий научно-инженерного сообщества на инновационном направлении в решении теплофизических проблем технологий энергетического производства. Особое внимание уделено вопросам энергосбережения и эффективности технологий сжигания органического топлива, надежности и безопасности технологических систем теплоэнергетики.

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

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

УДК 621.1(063) ББК 31.3л Редакционная коллегия Г.В. Кузнецов, доктор физико-математических наук

, профессор ТПУ;

А.С. Заворин, доктор технических наук, профессор ТПУ;

В.В. Литвак, доктор технических наук, профессор ТПУ;

В.С. Логинов, доктор физико-математических наук, профессор ТПУ;

П.Г. Яковенко, доктор технических наук, профессор ТПУ;

К.В. Буваков, кандидат технических наук, доцент ТПУ.

Конференция организована и проведена при финансовой поддержке Российского фонда фундаментальных исследований № 11-08-06045-г © ФГБОУ ВПО НИ ТПУ, ISBN 978-5-98298-959- УДК 621.182.002 (571.16) (092)

ОСНОВАТЕЛЬ ТОМСКОЙ НАУЧНО-ПЕДАГОГИЧЕСКОЙ

ШКОЛЫ КОТЛОСТРОЕНИЯ

Заворин А.С., д.т.н.

Томский политехнический университет, г. Томск E-mail: zavorin@tpu.ru Развитие теплоэнергетического образования в Томском политехническом институте во второй половине двадцатого века и становление в Томске и в Сибири в целом научно-педагогической школы по котлостроению связаны с именем профессора Ивана Кирилловича Лебедева, 100 лет со дня рождения которого исполнится 8 октября 2011 года. Результаты его деятельности в этих направлениях не утратили актуальности до сих пор и востребованы в современных условиях реформирования высшей школы.

Жизнь И.К. Лебедева охватила почти все эпохальные для России события двадцатого века. Родился 8 октября 1911 г. в Пензенской губернии в крестьянской семье. После гибели отца на фронте, в самом начале первой мировой войны, был взят на воспитание в семью дяди. В 1928 г. с мечтой стать инженером-энергетиком поступил на работу учеником электрика на шахте в г. Новошахтинске Ростовской области и получил возможность обучаться на рабфаке – так именовались курсы по подготовке представителей рабочего класса к учебе в вузах и техникумах. Успешное окончание И.К. Лебедев рабфака позволило не только поступать в институт без вступительных экзаменов, но и выбирать вуз.

Предпочтение было отдано старейшему техническому учебному заведению Сибири. Так в 1931 г. И.К. Лебедев стал студентом первого курса Сибирского механического института, который был выделен наряду с другими из Томского технологического института в ходе проходившего тогда разукрупнения многопрофильных учебных заведений. Вскоре начался обратный процесс преобразования системы высшего образования, и путем слияния мелких вузов образовался Томский индустриальный институт, который И.К. Лебедев и закончил в 1936 г. по кафедре профессора И.Н. Бутакова с дипломом инженера-теплотехника по специальности «Производство, распределение и промышленное использование тепловой энергии».

Самостоятельная деятельность для И.К. Лебедева началась в Новокузнецке в должности инженера ТЭЦ Кузнецкого металлургического комбината. Ранее здесь проходила студенческая производственная практика, и с тех пор комбинат привлекал своими масштабами индустриального гиганта и размахом созидания нового металлургического центра страны. Молодой специалист с энтузиазмом включился в работу, с большим желанием перенимая опыт старших коллег и участвуя в решении разнообразных и сложных производственных задач. В последующем, уже на преподавательском поприще, Иван Кириллович часто обращался к этому периоду, приводя многочисленный примеры практических ситуаций в своих лекциях, что делало их очень убедительными и запоминающимися.

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



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

После возвращения на КМК в 1940 г. И.К. Лебедев был назначен на должность заместителя начальника котельного цеха ТЭЦ комбината и оставался на этом посту до 1944 г. По воспоминаниям Ивана Кирилловича, это были годы неимоверно напряженного, поистине фронтового труда, когда для обеспечения военных заданий по выпуску металла энергетики порой сутками и неделями оставались на рабочих местах.В 1944 г. И.К. Лебедева по личному распоряжению министра черной металлургии И.Ф. Тевосяна переводят начальником котельного цеха на ТЭЦ Новотагильского металлургического комбината, а в 1947 г. по решению Свердловского обкома КПСС назначают начальником котельного цеха Богословского алюминиевого завода.

В 1948 г. началась научно-педагогическая деятельность Ивана Кирилловича: поступив в аспирантуру к профессору И.Н. Бутакову по его приглашению, он, как оказалось, окончательно связал свою жизнь с Томском и родным институтом. Здесь уместно отметить, что Иннокентий Николаевич Бутаков, который сам прошел богатую производственную школу до того как стать преподавателем, первым оценил высокую профессиональную квалификацию и организаторские способности своего ученика. Являясь в то время заместителем ректора (проректором) по учебному и научному управлению, он вскоре рекомендовал И.К. Лебедева ректору А.А. Воробьеву как потенциального лидера в деле подготовки и осуществления актуальной тогда для развития политехнического института задачи создания самостоятельного направления научной и образовательной деятельности на стыке теплофизики, энергетического машиностроения и теплотехники.

В 1951 г. И.К. Лебедев защитил кандидатскую диссертацию на тему «Борьба с золовым износом водяных экономайзеров и дымососов котельных установок при сжигании высокозольных топлив». Проведенное исследование, отличавшееся большой практической значимостью, было в русле систематических исследований теплотехнических свойств энергетического топлива, проводившихся школой профессора И.Н. Бутакова, и вместе с тем дало начало новому направлению – исследованию влияния минеральной части топлив на условия работы котельных агрегатов.

Заняв по конкурсу должность доцента кафедры теплосиловых установок, И.К. Лебедев стал совместно с доцентом В.Н. Смиренским курировать деятельность кафедры по котельным установкам. Одновременно с этим он с 1952 г. становится деканом энергетического факультета, сменив на этом посту профессора В.Т. Юринского.

В этот период И.К. Лебедев проявил себя как талантливый и целеустремленный организатор высшего профессионального образования.

В 1952 г. по его инициативе и при безоговорочной поддержке ректора А.А. Воробьева открыт прием по специальности «Котлостроение».

Однако уже в 1955 г. пришлось отстаивать его перед вышестоящими инстанциями, так как там преобладала точка зрения, что подготовку специалистов по энергетическому машиностроению, к которому относилась специальность, надлежит вести в индустриальных центрах, где действуют крупные производственные мощности соответствующего профиля. Иван Кириллович не отступил перед авторитетными решениями, он аргументировал тем, что главенствующую роль играет наличие высококвалифицированных кадров для фундаментальной, общеинженерной и специальной подготовки, а отсутствие привязки к конкретной производственной базе дает выгодные предпосылки готовить универсальных специалистов, ориентированных на проблематику всей отрасли, а не замкнутых на специфику конкретного предприятия. В стремлении закрепить в Сибири подготовку инженеров для энергомашиностроения он обратился за поддержкой на котлостроительные заводы, которые инициировали пересмотр планов потребности в специалистах Министерством тяжелого машиностроения СССР и направление их в Госплан и Министерство высшего и среднего специального образования. Полагая, что на нужное решение уйдет слишком много времени, И.К. Лебедев обратился с личным письмом к Первому секретарю Центрального комитета КПСС Н.С. Хрущеву [1]. В итоге плановый прием на специальность был не только сохранен, но даже увеличен в два раза.

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

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





В 1956 г. И.К. Лебедев подготовил и осуществил разделение энергетического факультета и стал первым деканом теплоэнергетического факультета, оставаясь на этом посту до 1959 г. О том, каким он был деканом, одновременно масштабно целеустремленным на перспективу и дотошно внимательным в ситуациях со студентами, известно немало из воспоминаний выпускников и коллег. О значении его личности как декана немало говорит, например, тот факт, что два, пожалуй, наиболее выдающихся из выпускников-энергетиков тех лет, являющиеся ныне действительными членами Российской Академии наук – лауреаты премии «Глобальная энергия» академик Геннадий Андреевич Месяц и академик Владимир ЕлиферьевичНакоряков – оба в своих официальных интервью назвали Ивана Кирилловича Лебедева среди тех, кто оказал наибольшее влияние на их судьбу и на путь в науку.

За время своего пребывания в должности декана теплоэнергетического факультета, а до того энергетического факультета, И.К. Лебедев организовал и возглавил работу по открытию специальностей «Котлостроение» (1952 г.), «Промышленная теплоэнергетика» (1955 г.), «Автоматизация теплоэнергетических процессов» (1959 г.), «Атомные электрические станции» (1959 г.) и кафедр «Парогенераторостроение и парогенераторные установки» (1958 г.), «Промышленная теплоэнергетика» (1959 г.), «Автоматизация теплоэнергетических процессов»

(1959 г.) [2].Он вникал буквально во все вопросы методического обеспечения учебного процесса, подготовки кадров преподавателей и создания материальной базы.

Особенно показателен принцип егопонимания ответственности и роли руководителя на примере организации подготовки по специальности «Атомные электрические станции». Принимая во внимание, что в подготовке этих инженеров определяющее значение имеет теплотехнический базис образования, ректорат не смог опираться на возможности физико-технического факультета. После консультаций с авторитетными специалистами (в частности, для изучения возможностей на месте, в ТПИ была командирована профессор МЭИ Т.Х. Маргулова) декану ТЭФ И.К. Лебедеву было предложено разработать решение вопроса на базе теплоэнергетического факультета. В итоге подготовка по специальности «Атомные электрические станции» была им организована пригод назад созданной и на деле близкой по профилю кафедре «Парогенераторостроение и парогенераторные установки», которой он сам и заведовал. Понятно, что это решение существенно осложняло работу кафедры и прежде всего заведующего кафедрой, однако не смогло повлиять на системный и последовательный характер проводимых им мероприятий. Практически одновременно с приемом на первый курс студентов, начались предметные стажировки преподавателей в Московском энергетическом институте, на первой АЭС в г. Обнинске, ежегодно из лучших выпускников проходил отбор для направления в целевую аспирантуру на профильные кафедры московских вузов. Многое в учебном процессе переносилось из проверенного опыта подготовки по специальностям «Тепловые электрические станции» и «Котлостроение».

Результат не заставил себя ждать: уже ко времени первого выпуска инженеров по атомным станциям (1964 г.) во всех деталях работала система их воспроизводства, которая в последующие годы только проходила дальнейшую огранку. Кафедра парогенераторостроения и парогенераторных установок профилировала подготовку по атомным электростанциям вплоть до 1971 г., затем она была передана на образованную тогда кафедру «Теплофизика и атомная энергетика».

Отдаленные во времени результаты образовательной деятельности дают наиболее объективную её оценку. В этом отношении интересен такой факт: в середине восьмидесятых годов, т. е. через 20 лет после первого выпуска кафедрой И.К. Лебедева инженеров для атомной энергетики, из 18-ти АЭС, задействованных в энергетике, на 12-ти в числе первых руководителей (директор или главный инженер, либо сразу оба) были выпускники Томского политехнического института, а из них на 9ти АЭС выпускники теплоэнергетического факультета.

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

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

В начале шестидесятых годов, когда завершилось формирование контингента студентов и преподавателей на теплоэнергетическом факультете, обозначились масштабы учебной и научной работы, которые превзошли возможности имевшихся площадей. Политехническому институту вновь пригодились организаторский талант и инженерный опыт И.К. Лебедева. К тому времени институтская электростанция, расположенная в четвертом корпусе, утратила уже свое значение для обеспечения электроэнергией и оставалась только источником теплоснабжения для группы корпусов, включая главный. И.К. Лебедев обосновал необходимость реконструкции 4-го корпуса. После принятия положительного решения во всех инстанциях института и города он лично на общественных началах возглавил разработку проекта реконструкции, а затем и строительство. А до того требовалось перевести пять корпусов института на централизованное теплоснабжение от ГРЭС-2. И.К. Лебедев стал общественным прорабом и техническим руководителем форсированного строительства теплотрассы вдоль улицы Усова. Были привлечены со всех факультетов студенты-добровольцы, которые в зимнее время круглосуточно, смена за сменой, отогревая кострами замерзшую землю, вырыли траншею для прокладки трубопроводов теплотрассы. Силами студентов велись и строительные работы на корпусе, которые, как и поставку материалов, И.К. Лебедев каждодневно контролировал. В 1966 г.

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

Хотя И.К. Лебедев приступил к научно-педагогической деятельности уже в зреломвозрасте, он, несмотря на выпавшую ему большую занятость организационными делами, сумел внести заметный вклад в науку. Практически с самого начала освоения для энергетики крупнейшего в мире Канско-Ачинского буроугольного бассейна он вместе с коллективом кафедры занялся научно-техническимипроблемами этого направления, актуального для развития Сибири и всей страны. Вскоре были получены результаты, которые привлекли внимание ведущих отраслевых институтов, инициировали финансирование исследований кафедры и включение их в государственные координационные планы. Всё это позволило привлекать средства для расширения экспериментальной базы, что, в свою очередь, создавало условия для расширения спектра и объема научных исследований, развития на этой основе кадрового состава. Уже к концу шестидесятых годов коллектив кафедры был признан одной из ведущих научных школ страны в области технологии и техники энергетического использования низкосортных топлив. Свидетельством тому стало приглашение И.К. Лебедева с научным докладом на VIII Конгресс Мировой энергетической конференции (1968 г.) [3].

В 1971 г. И.К. Лебедевым была защищена докторская диссертация на тему «Особенности сжигания углей Канско-Ачинского бассейна в топках энергетических котельных агрегатов большой паропроизводительности». Его рекомендации использованы отечественными котлостроительными заводами при разработке ряда конструкций парогенераторов. По тематике научного направления им подготовлено 20 кандидатовнаук, многие из которых впоследствии защитили докторские диссертации.

В соавторстве с учениками И.К. Лебедевым опубликовано более 150 трудов, он подготовил несколько учебных пособий, среди которых выделяется изданное в «Энергоатомиздате» с грифом Министерства высшего и среднего образования СССР пособие «Гидродинамика паровых котлов», которое до сих пор является основным учебным изданием по дисциплине «Гидродинамика энергоустановок» для студентов направления «Энергетическое машиностроение».

И.К. Лебедев оставался бессменным заведующим кафедрой парогенераторостроения и парогенераторных установок до 1988 г.Он курировал подготовку научных кадров, консультировал ученых и специалистов предприятий Барнаула, Бийска, Иркутска, Красноярска, Новокузнецка, Новосибирска, Омска, Прокопьевска, Читы, содействовал становлению учебно-методической базы для подготовки котлостроителей в Алтайском политехническом институте.

За высокие производственные достижения и результативную научно-педагогическую деятельность И.К. Лебедев был награжден тремя орденами «Знак Почета» (1945, 1961, 1981 гг.), четырьмя медалями, Почетной грамотой Верховного Совета РСФСР, почетными знаками «Отличник энергетики» и «Отличник черной металлургии»,дипломом почета ВДНХ СССР. В 1982 г. ему было присвоено почетное звание «Заслуженный деятель науки и техники РСФСР», в 1993 он был избран членом-корреспондентом Сибирского отделения Международной Академии наук высшей школы.

До последних дней жизни Иван Кириллович оставался действующим преподавателем. В знак признания заслуг профессора перед Томским политехническим университетом по решению Ученого совета на учебном корпусе № 4 установлена мемориальная доска.

1. Заворин А.С., Беляев С.А. Томская школа котло- и реакторостроения. – Томск:

Изд-во «S-press», 2008. – 192 с.

2. Беляев Л.А., Беляев С.А., Заворин А.С. Итоги и традиции теплоэнергетического образования в ТПУ: 120 лет со дня рождения И.Н. Бутакова и 90 лет со дня рождения И.К. Лебедева // Известия Томского политехнического университета. – Т. 305. – Вып. 2. – 2002. – С. 7–14.

3. Заворин А.С., Загромов Ю.А., Беляев Л.А. Развитие энергоэффективных технологий в теплоэнергетике // Известия Томского политехнического университета. – Т. 303. – Вып. 2. – 2003. – С. 78–84.

ТЕПЛОФИЗИЧЕСКИЕ И ЭНЕРГО-ЭКОЛОГИЧЕСКИЕ

ХАРАКТЕРИСТИКИ ГОРЕНИЯ ВОДОРОДА В ДВИГАТЕЛЯХ

ВНУТРЕННЕГО СГОРАНИЯ

Ассад М.С., д.т.н., Пенязьков О.Г., д.ф-м.н.

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

Использование водорода в тепловых двигателях, в частности – двигателях внутреннего сгорания, приводит к существенным изменениям параметров рабочего процесса и влияния энергоустановки на окружающую среду. Для выявления характера и величин этих изменений в лаборатории физико-химической гидродинамики ИТМО НАН Беларуси выполнен цикл экспериментов на автоматизированном стенде, описанном в [1], ядром которого является двигатель HONDA – D15B6. В ходе работы снимались теплофизические, топливно-энергетические и экологические характеристики двигателя в зависимости от содержания водорода в смеси. Контролировались следующие параметры: давление и интенсивность излучения пламени в камере сгорании, объем поступающего в двигатель воздуха, объемное и массовое содержание водорода и синтезгаза в смеси, расход бензина, мощность и крутящий момент двигателя, количественный состав (CO, CO2, O2, CHx и NO) и температура продуктов сгорания, частота вращения коленчатого вала, коэффициент избытка окислителя и другие (в общей сложности 27 параметров). Методика проведения эксперимента описана в [2, 3]. Суть эксперимента заключалась в выходе на заданный режим работы на основном топливе (бензине) и последующее увеличение подачи водорода от 0 до 20 % от объема поступающего в двигатель воздуха.

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

Обогащение бензина водородом в количестве 20 % от объема поступающего в двигатель воздуха приводит к некоторому увеличению осредненного максимального давления (по 70 циклам) примерно на 8 % при частоте вращения коленчатого вала n = 1000 мин –1 и на 11 % при n = 3000 мин –1 по сравнению с работой на чистом бензине, что объясняется высокой скоростью горения водорода, которая, по нашим данным [4, 5], почти на порядок выше, чем у бензина.

Зависимости, характеризующие содержание СО, CO2, O2, CHx и NO в продуктах сгорания исследуемого двигателя, который работал на бензине с добавками водорода 0–8 % от объема поступающего в двигатель воздуха, представлены на рисунках 1 и 2. Там же даны зависимости изменения удельного расхода топлива Gb = f(H2), мощности двигателя Ne = f(H2), частоты вращения коленчатого вала n = f(H2) и коэффициента избытка окислителя = f(H2).

Рис. 1. Зависимость эмиссии CO и NO, а также динамики изменения расхода бензина Gb и мощности двигателя Ne от размера добавки H2 к воздуху Рис. 2. Зависимости эмиссии CHx, O2 и CO2, а также изменения коэффициента избытка окислителя и частоты вращения коленчатого вала n от размера добавки Н2 к воздуху С увеличением доли H2 в топливно-воздушной смеси, т. е. по мере уменьшения отношения С/Н в элементном составе топлива, содержание CO и CHx существенно падает и при 8 %-й добавке водорода концентрации этих компонентов ниже на 17 и 27 % соответственно, чем при работе двигателя на бензине (кривая 1 рис. 1 и кривая 2 рис. 2). Однако наличие водорода в смеси ведет к возрастанию содержания NO в продуктах сгорания (кривая 2 рис. 1). Например, при 8 %-й добавке H2 значение концентрации NO достигает 585 ppm, что в 2,15 раза выше, чем при работе на бензине.

Интерес представляет анализ влияния водорода на топливноэнергетические показатели двигателя. При питании двигателя бензином с добавкой водорода до 12 % по объему поступающего в двигатель воздуха и коэффициенте избытка воздуха 1 мощность практически сохраняется на уровне мощности базового двигателя (отклонение в пределах 4 %). Дальнейшее обогащение H2 до 20 % приводит к снижению мощности до 9 %. При этом существенно возрастает топливная экономичность. Так, при добавлении 6 %, 12 % и 20 % H2 расход бензина снижается примерно на 25 %, 43 % и 70 % соответственно. Объясняется это кинетикой и высокой скоростью горения водорода, которые положительно влияют на процесс горения всей смеси и способствуют уменьшению тепловых потерь через стенки цилиндра.

Используя методы математического планирования, описанные в [6], получены адекватные уравнения регрессии второго порядка, позволяющие достаточно точно оценить состав продуктов сгорания двигателя в зависимости от частоты вращения коленчатого вала, коэффициента избытка окислителя и размера добавки водорода в смеси.

Для установления динамики процентного прироста выброса (повышения или понижения), соответствующего изменению доли водорода на 1 % вычислена эластичность концентрации продуктов сгорания как предел отношения относительного приращения концентрации i-го компонента ri к относительному приращению доли водорода в смеси, когда приращение аргумента стремится к нулю (если производная исследуемой функции существует):

Таким образом, определены термодинамические характеристики горения внутри цилиндра, установлены особенности протекания рабочего процесса при использовании водородсодержащих топлив и изучено их влияние на мощность, топливную экономичность и состав продуктов сгорания (CO, CO2, CHx, O2 и NO) двигателя с распределенным впрыском при работе на разных скоростных и нагрузочных режимах.

1. Ассад М.С., Пенязьков О.Г. Особенности работы двигателя внутреннего сгорания на водородсодержащих топливах // ИФЖ. – 2010. – Т. 83. – № 4. – С. 814–820.

2. Ассад М.С., Пенязьков О.Г. Продукты сгорания жидких и газообразных топлив.

Образование, расчет, эксперимент: Монография. – Минск: Беларус. навука, 2010. – 305 с.

3. Ассад М.С., Пенязьков О.Г. Методика измерения состава продуктов сгорания водородсодержащих топлив в двигателе внутреннего сгорания // Тепло- и массоперенос – 2008. – Минск: ИТМО им. А.В. Лыкова НАН Беларуси, 2009. – С. 76–81.

4. Ассад М.С., Лещевич В.В., Миронов В.Н., Пенязьков О.Г., Севрук К.Л. Горение модифицированных топлив в модели камеры сгорания ДВС // Тепло- и массоперенос – 2005. – Минск: ИТМО им. А.В. Лыкова НАН Беларуси, 2005. – С. 100–105.

5. Ассад М.С., Лещевич В.В., Миронов В.Н., Пенязьков О.Г., Севрук К.Л., Скилондь А.В. Горение водородсодержащих смесей в модели камеры сгорания ДВС // ИФЖ. – 2009. – Т. 82. – № 6. – С. 1031–1045.

6. Ассад М.С. Оценка влияния добавок водорода на экологические показатели ДВС методами математического планирования // Весцi НАН Беларусі. Сер. фіз.-тэхн.

навук. – 2010. – № 3. – С. 84–89.

МОДЕЛИРОВАНИЕ ФИЗИКО-ХИМИЧЕСКИХ ПРОЦЕССОВ

ПРИ ЗАЖИГАНИИ ТОНКИХ ПЛЕНОК ЖИДКИХ

КОНДЕНСИРОВАННЫХ ВЕЩЕСТВ ИСТОЧНИКАМИ НАГРЕВА

С ОГРАНИЧЕННЫМ ТЕПЛОСОДЕРЖАНИЕМ

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

Установлена многофакторность процессов тепломассопереноса при воздействии на горючие жидкости типичных источников ограниченной энергоемкости – нагретых до высоких температур частиц металлов и неметаллов. Разработаны плоская [3] и пространственная [4] модели процессов зажигания жидкостей одиночными частицами. Но и плоская [3], и особенно пространственная [4] модели являются слишком сложными. Для их реализации необходимы большие вычислительные ресурсы. По этим причинам представляется целесообразным разработка упрощенных моделей и методов численного расчета условий зажигания горючих и легковоспламеняющихся жидкостей.

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

Предполагалось, что в начальный момент времени разогретая частица находится на поверхности жидкого вещества. Считалось, что источник нагрева не погружается в вещество, а на границе «частица – жидкость» обеспечивается идеальный контакт. За счет энергии частицы жидкость разогревается. Интенсифицируется испарение вещества.

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

Рис. 1. Схема области решения задачи на первом (a) и втором 2 – частица, 3 – жидкость уравнение диффузии паров горючего в воздухе (y2 < y < l):

уравнение баланса парогазовой смеси (y2 < y < l):

уравнение теплопроводности для частицы (y1 < y < y2):

уравнение теплопроводности для пленки жидкости (0 < y < y1):

Здесь C – удельная теплоёмкость, Дж/(кг·К); Cг – концентрация паров Принимая во внимание, что для исследуемой конструкции имеет задачу (15)–(21) удается разделить на две самостоятельные более простые задачи. Первая, из которых запишется:

Данная система уравнений равноценна задаче нестационарной теплопроводности неограниченной пластины при нулевых граничных условиях и неравномерном распределении температуры в начальный момент времени. Аналитическое решение системы (22)–(25) известно и приводится, в частности, в монографии А.В. Лыкова [3].

Используя функцию температуры 1 (l1, ) из этого решения, можно составить систему уравнений для нахождения 2 путем её подстановки в оставшуюся часть задачи (15)–(21):

Задача (26)–(29) также может рассматриваться как процесс нестационарной теплопроводности плоского тела при соответствующих граничных и начальном условиях. Математическое решение такой системы тоже приводится в [3].

1. Керн Д., Краус А. Развитые поверхности теплообмена. – М.: Энергия, 1977. – 461 с.

2. Видин Ю.В., Бойков Г.П., Колосов В.В., Ромащенко А.С. Краткий справочник по тепломассообмену. – Красноярск: Сибирский федеральный ун-т, 2007. – 169 с.

3. Лыков А.В. Теория теплопроводности. – М.: Высшая школа, 1967. – 600 с.

НЕСТАЦИОНАРНЫЙ РАДИАЦИОННО-КОНВЕКТИВНЫЙ

ТЕПЛООБМЕН В ГАЗОПРОНИЦАЕМОМ МАТЕРИАЛЕ

C ИСТОЧНИКАМИ ТЕПЛОВЫДЕЛЕНИЯ

Бурка А.Л., к.т.н., Емельянов А.А., к.т.н., Синицын В.А., к.т.н.

Институт теплофизики им. С.С. Кутателадзе СО РАН, г. Новосибирск В настоящей работе предлагается расчетно-экспериментальное исследование механизма накопления тепловой энергии с высокой удельной плотностью в объеме газопроницаемого материала. Особый интерес представляют системы частиц двуокиси ванадия. Известно, что нестехиометрические оксиды VO2-x состоят из большого ряда различных структур, называемых фазы с переменной валентностью. В этих структурах решеточный кислород слабо связан и может быть хемодесорбирован при довольно низких температурах, что приводит к изменению кристаллической решетки оксида. Структурные изменения с образованием новых капилляров, внутренних границ раздела фаз и новых фаз переменной валентности в частицах могут быть достигнуты с помощью механической активации материала.

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

С целью изменения количества фаз переменной валентности и адсорбированного слабосвязанного кислорода порошок двуокиси ванадия производства фирмы «Aldrich», выдержанный на воздухе при комнатной температуре в течение десяти месяцев, подвергался обработке на планетарно-центробежной мельнице в воздушной среде. Удельная поверхность частиц двуокиси ванадия была измерена методом тепловой десорбции аргона. Рентгенофазовый анализ образцов произведен на приборе «Дрон-4» [1].

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

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

Рис. 1. Эффективная теплоемкость двуокиси ванадия при различных значениях плотности энергии активации. Образец 1 (Aldrich) выдерживался на воздухе, остальные образцы находились без В работе проводится исследование процесса переноса тепла и массы в пористом материале в условиях, когда теплообмен излучением является сопоставимым с теплопроводностью в твердом веществе и конвекцией газа в порах. Плоский слой газопроницаемого полупрозрачного материала из двуокиси ванадия нагревается излучением от внешних серых источников. Происходит теплообмен с окружающей средой излучением и конвекцией. В условиях нагрева или охлаждения твердый материал испытывает фазовый переход второго рода «металлдиэлектрик», происходит хемосорбция кислорода и перестройка кристаллической структуры. Эти процессы сопровождаются отводом или подводом тепла в материале, что приводит к появлению объемных источников или стоков массы и тепла в слое. При моделировании радиационного переноса в газопроницаемом материале учитываются процессы поглощения, испускания и рассеяния теплового излучения. Оптические свойства материала не зависят длины волны. Теплоемкость материала зависит от температуры. Решение задачи проводится в приближении сплошной среды с эффективными теплофизическими и оптическими свойствами.

Уравнение энергии с граничными условиями записываются в виде:

Удельный расход газа g находится из уравнения с граничным условием Дивергенция интегральной плотности радиационного потока E / в уравнении (1) определяется из решения системы интегродифференцальных уравнений переноса энергии излучения [2].

С помощью функции Грина G (, z ), которая является решением однородной краевой задачи:

начально-краевая задача (1)(4) сводится к нелинейному интегральному уравнению относительно безразмерной температуры, которое записывается в виде:

Это уравнение на каждом временном шаге решается итерационным методом Ньютона-Канторовича [3].

Расчеты проводились при следующих значениях параметров, соотР = 0,8, ветствующих условиям эксперимента:

E 4.4 Вт/мК, g 0.026 Вт/мК, c g 1005Дж/кгК, 4339 кг/м3.

Результаты расчета, представленные на рисунке 2, позволяют выяснить влияние температуры внешнего источника на формирование температурного поля в процессе нагрева материала. Температура источника, облучающего границу 0, принимала значения Т1 = (127, 227) °С, а температура источника, облучающего границу 1, Т2 = 127 °С. Результаты расчета относятся к фиксированному моменту времени t =720 с от начала нагрева. Видно возрастание температурного уровня всех образцов с увеличением температуры источника, что объясняется ростом потока падающего излучения. При этом наблюдается большое расслоение температурных кривых по величине активации материала.

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

В температурном интервале 50–140 °С экспериментально исследовано калориметрическим методом влияние активации двуокиси ванадия в воздушной среде на важнейшую термодинамическую характеристику – эффективную теплоемкость материала. Показано, что с увеличением энергии активации происходит смещение максимума эффективной теплоемкости в область более высоких температур с одновременным существенным ростом его абсолютной величины. Расчеты показали, что аккумулирующий эффект в материале проявляется наиболее сильно при температурах, соответствующих высоким значениям эффективной теплоемкости.

Работа проводилась при поддержке Российского фонда фундаментальных исследований по проекту 06-08-00361-а.

Обозначения:

cр эффективная теплоемкость материала, Дж/(кгK); cg теплоемкость газа, Дж/(кгK); Е плотность потока излучения, Вт/м2; L толщина слоя, м; mV массовая скорость газовыделения в порах, кг/(м3с); Р пористость; Т температура, К; T1*, T2* температуры внешних источников излучения, K; t время, с; х координата в перпендикулярном направлении, м; 1, 2 коэффициенты теплоотдачи на поверхностях слоя, Вт/(м2 K); объемный коэффициент поглощения, 1/м; 1, степени черноты поверхностей слоя; коэффициент эффективной теплопроводности материала, Вт/(мK); g коэффициент теплопроводности газа, Вт/(мK); E коэффициент теплопроводности твердого каркаса, Вт/(мK); плотность материала, кг/м3.

1. Емельянов А.А., Полубояров В.А., Бурка А.Л., Коротаева З.А., Великанов Е.В., Лапин А.Е., Бан Бонг-Чан. Теплообмен при фазовом переходе и адсорбциидесорбции кислорода в дисперсном материале // Промышленная теплотехника. – 2006. – Т. 28. – № 1. – С. 22–30.

2. Оцисик М.Н. Сложный теплообмен. – M.: Мир, 1976. – 616 с.

3. Канторович Л.В. О методе Ньютона // Труды Математического института АН СССР. – 1949. – Т. 28.

МОДЕЛИРОВАНИЕ ТЕПЛООБМЕНА В ПОЛУПРОЗРАЧНОМ

КОМПОЗИТЕ С МИКРОСФЕРИЧЕСКИМИ ЧАСТИЦАМИ

Бурка А.Л., к.т.н., Емельянов А.А., к.т.н.

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

Композитный материал является эффективной теплозащитой только при высокой объемной концентрации микросфер, вследствие чего не выполняется условие независимости рассеивающих и поглощающих излучение частиц. Поэтому использование классической теории Ми [1] в этом случае может привести к значительным ошибкам при расчете радиационных характеристик композитного материала и температурных полей в нем. Для определения ослабляющих и рассеивающих характеристик композитного материала при наличии взаимодействий между частицами могут быть использованы приближенные решения [2]. В предлагаемой работе учет влияния взаимодействия между микросферами на коэффициенты ослабления и рассеяния в композитном материале осуществляется применением эффективных оптических свойств полученных на основании приближения Максвелл-Гарнетта. Радиационнокондуктивный теплообмен в рассматриваемом материале описывается уравнением энергии с граничными условиями в безразмерном виде:

Дивергенция плотности радиационного потока dE / d, которая находится в уравнении (1), определяется из решения системы интегродифференциальных уравнений для переноса энергии излучения относительно интенсивностей I (,), I (,) [3]. С помощью функции Грина [3] начально-краевая задача (1)(4) о совместном переносе тепла теплопроводностью и радиацией в полупрозрачном материале сводится к нелинейному интегральному уравнению относительно безразмерной температуры (,), которое имеет вид Эффективные значения плотности и коэффициента теплопроводности композита определялись [4]. Интегральное уравнение (5) на каждом временном шаге решается итерационным методом НьютонаКанторовича [5]. Расчеты нестационарных температурных полей в слое композита проводились с учетом переноса тепла в материале за счет теплопроводности и излучения. Были выбраны следующие значения параметров: * = 1, 1 2 5 Вт/м2К, T* = 1000K, L = 0,1 м,, получены с учетом фактора заполнения, = +, g = 0,196 Вт/мК, E = 0,343 Вт/мК, 0 () 0.333, c1 = 783Дж/кгК, c2 = 1465Дж/кгК, 1 = 118 кг/м, 2 = 2200 кг/м, температуры источников внешнего излучения: T1* 500 0 C, T2* 300 0 C.

На рисунке 1 представлены температурные поля в композитах в отсутствии взаимодействия между сферами (N = 2,5· 10 6 м 3 ). Из рисунка видно, что величина внешнего диаметра микросферы заметно влияет на температурное распределение. В композите с малыми микросферами (кривые 2) практически отсутствует температурный перепад между горячей ( 0 ) и холодной ( 1 ) границами слоя, а в композите с большими микросферами (кривые 1) температурный перепад заметен и растет в процессе нагрева. В этом случае рассматриваются образцы с одинаковыми факторами заполнения, поэтому их эффективные теплопроводности практически совпадают, а различие в поведении температурных профилей будет определяться в основном радиационным переносом тепла. Рассеяние излучения в композите с малыми частицами приводит к меньшему перепаду температур на границах слоя. На рисунке представлены температурные кривые для композитов с факторами заполнения, при которых расстояние между микросферами становится соизмеримым с расстоянием начала их взаимодействия (расстояние между микросферами в этом случае составляет 3d). Из рисунка 2 видно, что при нагреве композитов с микросферами d = 3,5 мкм наблюдается начало расслоения температурных кривых для 0,0056, 0,00056 и тем заметнее, чем выше температура.

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

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

Рис. 1. Температурное распределение в композитах, содержащих полые Рис. 2. Температурное распределение в композитах, содержащих полые Обозначения:

c, c1, c2 – эффективная теплоемкость композита, кварца, полиметилметакрилата, Дж/(кг·К); Е – плотность потока излучения, Вт/м2; Ib – функция Планка абсолютно черного тела при температуре T, Вт/(м3стер); L – толщина слоя, м; T* – характерная температура, К; R – радиус сферической частицы, м; qR – радиус ядра сферической частицы; T1*, T2* – температуры внешних источников излучения, К; 1, 2 – коэффициенты теплоотдачи на граничных поверхностях слоя, Вт/(м2K);

– объемная концентрация частиц (фактор заполнения); – объемный коэффициент поглощения композита, 1/м; – объемный коэффициент рассеяния композита, 1/м; – степень черноты слоя композита; 1*, 2 – степени черноты внешних источников излучения; 1, 2 – приведенные степени черноты границ слоя композита.

1. Дейерменджан Д. Рассеяние электромагнитного излучения сферическими полидисперсными частицами. – М.: Мир, 1971. – 300 с.

2. Петров Ю.И. Физика малых чтиц. – М.: Наука, 1982. – 360 с.

3. Бурка А.Л., Емельянов А.А. Влияние взаимодействия микросфер на формирование температурных полей в полупрозрачном материале // Промышленная теплотехника. – 2011. – Т. 33. – № 1. – С. 24–30.

4. Полежаев Ю.В., Юревич Ф.Б. Тепловая защита. – М.: Энергия, 1976. – 392 с.

5. Канторович Л.В. О методе Ньютона // Труды Математического института АН СССР. – 1949. – Т. 28. – С. 135–139.

ИЗМЕРИТЕЛЬНАЯ СИСТЕМА, РЕАЛИЗУЮЩАЯ МЕТОД

ОПРЕДЕЛЕНИЯ КОМПЛЕКСА ТЕПЛОФИЗИЧЕСКИХ

СВОЙСТВ МАТЕРИАЛОВ

Рогов И.В., к.т.н., Кацуба Д.С., Полунина Н.Ю., Полунин Е.П.

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

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

На практике для определения ТФС применяют различные приборы и средства измерения. Наиболее перспективными являются измерительные системы, которые обеспечивают автоматический сбор и обработку измерительной информации об исследуемом материале [1].

Разрабатываемая измерительная система предназначена для определения комплекса ТФС: теплопроводности (), температуропроводности (а) и теплоемкости (с).

Схема измерительной системы представлена на рисунке 1.

Рис. 1. Схема измерительной системы теплофизического эксперимента:

1 – теплоизмерительная ячейка, 2 – термостат, 3 – циркуляционный насос, 4 – электродвигатель, 5 – блок реле, 6 – блок управления, 7 – персональный компьютер, 8 – блок холодных спаев термопар Управление и регистрация измерительной информации осуществляется с помощью персонального компьютера (ПК) 7, оснащенного многофункциональной платой сбора данных.

Образцы испытуемого материала располагаются с обеих сторон от сердечника теплоизмерительной ячейки (ТИЯ) 1 и прижимаются с наружной стороны металлическими блоками, в которых с помощью насоса 3 циркулирует из термостата 2 нагреваемая жидкость. Управление электродвигателем 4 насоса осуществляется через контакты реле а блока реле 5.

На аналого-цифровой преобразователь (АЦП) платы поступают сигналы с измеренными значениями термо-ЭДС от термопар, расположенных на внутренних стенках металлического блока и сердечника ТИЯ.

Поддержание заданного температурного режима в термостате осуществляется при помощи контактов реле б блока 5. Значение температуры жидкости в термостате измеряется при помощи термопары, подключенной через блок холодных спаев 8 к АЦП платы. Напряжение на электрический нагреватель термостата подается от блока управления 6.

Реле а и б, входящие в состав блока 5, программно управляются ПК через дискретные выходы платы (ДВ).

Теоретическую основу определения ТФС материалов составляет метод, основанный на решении краевой задачи теплопроводности [2].

Исходные данные для такой задачи формируются на основе результатов измерений и включают в себя граничные условия (первого и четвертого рода) и зависимость разности значений температуры в двух точках исследуемого образца от времени вида Т f ().

Исследуемые образцы имеют форму диска. За счёт конструктивного исполнения ТИЯ и равномерного нагрева поверхности образцов реализуется одномерный процесс теплопереноса. Температура обратной поверхности образца (со стороны металлического блока ТИЯ) используется для задания граничного условия.

На основании решений краевой задачи теплопроводности получено расчетное выражение, описывающее термограмму на рабочем участке:

где Тс – значение температуры на поверхности исследуемого образца, К;

Т0 – начальное распределение температуры внутри исследуемого образца, К; с11h1 ch – относительная теплоемкость материала сердечника ТИЯ;, с,, а – теплопроводность, Вт/(мК), удельная теплоемкость, Дж/(кгК), плотность, кг/м3, температуропроводность, м 2/с, образца; с1, 1 – удельная теплоемкость, Дж/(кгК) и плотность, кг/м3, материала сердечника; h, h1 – значения толщины образца и сердечника, м;

– время, с.

Для расчета ТФС по экспериментальным данным представим выражение (1) в виде:

Здесь параметры модели, описывающей термограмму на рабочем участке;

сп – теплоёмкость, отнесённая к единице площади поверхности сердечника, Дж/(м2К), c п с11 h1.

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

Представленная измерительная система реализует современный метод определения ТФС материалов.

Создание подобных систем является перспективным направлением в практике теплофизических исследований и тепловых испытаний теплоизоляционных материалов.

1. Виноградова Н.А., Гайдученко В.В., Карякин А.И. и др. Системы автоматизации теплофизического эксперимента: Учеб. пособие / Под ред. С.Г. Свиридова. – М.:

Издательский дом МЭИ, 2007. – 252 с.

2. Алифанов О.М. Обратные задачи теплообмена. – М.: Машиностроение, 1988. – 280 с.

ДВА РЕЖИМА ТЕПЛОМАССОПЕРЕНОСА ПРИ ПОДЗЕМНОЙ

ГАЗИФИКАЦИИ БУРОГО УГЛЯ

Мазаник А.С., Субботин А.Н., к.ф.-м.н.

Томский политехнический университет, г. Томск E-mail: nucleariya@pochta.ru; subbot@vtomske.ru Переход к рыночной экономике с особой остротой обнаружил отставание угольной отрасли как по технико-технологическому уровню, так и в области экономического и организационного управления. Интенсивная эксплуатация недр ставит проблему развития направления комплексного их освоения и разработку ресурсосберегающих технологий, позволяющих экономить минеральные ресурсы при одновременном снижении трудовых, энергетических и экономических затрат. Этим требованиям полностью удовлетворяет технология подземной газификации угля.

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

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

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

С учетом испарения влаги, разложения угля, гетерогенных и гомогенных реакций [1, 2]:

уравнения, моделирующие рассматриваемые процессы будут иметь вид [3]:

Система дифференциальных уравнений (1)–(5) решалась при краевых условиях Здесь уравнения (1) описывают законы сохранения массы исходного конденсированного вещества (угля), влаги, кокса и золы соответственно; (2) – уравнение сохранения массы газовой фазы; (3) – уравнения сохранения массы отдельных компонентов газовой фазы ( 1 – кислорода, 2 – оксида углерода, 3 – диоксида углерода, 4 – водорода, 5 – метана, 6 – водяного пара); последнее уравнение в (2) – алгебраическое соотношение необходимое для замыкания системы уравнений; (4) – уравнение сохранения энергии; (5) – уравнение движения в форме Дарси в проекциях на цилиндрические оси координат r, z и уравнение состояния; (6) – начальные условия; (7), (8) граничные условия вдоль координатных осей z и r соответственно.

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

Установлено, что меняя концентрацию кислорода в подаваемой в угольный пласт паровоздушной смеси можно менять режим протекания процесса подземной газификации. При сравнительно небольшой массовой концентрации кислорода равной 0,1–0,15 газификация протекает при температуре не превышающей 800 К (рис. 2), а массовая концентрация горючего газа (СО) в образующихся газообразных продуктах не превышает 27 %. Концентрации водорода и метана при этом малы: концентрация водорода не превышает 0,27 %, а метана менее 1,2 %.

Рис. 2. Типичное распределение температуры в угольном пласте при массовой концентрации кислорода равной 0,1–0, Увеличение массовой концентрации кислорода до 0,25 и выше, в закачиваемой в угольный пласт паровоздушной смеси, приводит к смене режима протекания процесса (рис. 3).

Температура газификации для данного режима возрастает до 1300 К даже при условии подачи холодной паровоздушной смеси. Массовая концентрация горючего газа (СО) в образующихся газообразных продуктах возрастает до 45 %. Концентрации водорода и метана также возрастают, но незначительно: концентрация водорода повышается до 1,8 %, а метана до 1,3 %.

Рис. 3. Типичное распределение температуры в угольном пласте при массовой концентрации кислорода равной 0, Таким образом, рассмотренный подход анализа процессов подземной газификации угля на базе математического моделирования тепломассопереноса и физико-химических превращений [3], протекающих при локальном нагреве угольных пластов и подводе окислителя, позволяет делать оценку состава образующейся газообразной смеси, прогнозировать режимы протекания данных процессов в условиях ПГУ и показывает существование двух режимов: низкотемпературный, протекающий при подаче подогретой паровоздушной смеси с низким содержанием кислорода; высокотемпературный, реализующийся при большом содержании кислорода в подаваемой паровоздушной смеси.

1. Померанцев В.В., Арефьев К.М., Ахмедов Д.Б. и др. Основы практической теории горения: Учебное пособие для вузов; 2-е издание переработанное и дополненное. – Л. : Энергоатомиздат, 1986. – 312 с.

2. Щетинков Е.С. Физика горения газов. – М.: Наука, 1965. – 739 с.

3. Кузнецов Г.В., Субботин А.Н. Тепломассоперенос в условиях подземной газификации угля // Тепловые процессы в технике. – 2010. – Т. 2. – № 9. – С. 422–426.

ИЗМЕРИТЕЛЬНАЯ СИСТЕМА НЕРАЗРУШАЮЩЕГО

КОНТРОЛЯ ТЕПЛОФИЗИЧЕСКИХ СВОЙСТВ

Жуков Н.П., д.т.н., Рогов И.В., к.т.н., Попов О.Н., Бардадымова О.Н.

Тамбовский государственный технический университет, г. Тамбов Сложность и большой объем экспериментальных исследований по определению качества, долговечности и надежности изделий, как из традиционных, так и вновь синтезируемых материалов, требуют создания новых эффективных методов и средств контроля. Среди них особое место занимают тепловые методы неразрушающего контроля (НК) и диагностики, позволяющие определять теплофизические свойства (ТФС) таких материалов. Известно, что современные методы теплового НК наиболее эффективно реализуются измерительными системами (ИС), позволяющими автоматизировать проведение измерений, адаптивно изменять режимные параметры, обеспечивать адекватность физической и математической моделей теплового процесса, обеспечивать оперативность и точность измерений при сохранении целостности и эксплуатационных характеристик объектов исследования [1].

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

Структурная схема ИС, реализующей НК ТФС материалов, представлена на рисунке 1. В состав ИС входят: персональный компьютер (ПК), измерительно-управляющая плата PCI-1202H, сменные измерительные зонды (ИЗ), регулируемый блок питания (БП). Зонд обеспечивает создание теплового воздействия на исследуемое изделие, фиксирование температуры в заданных точках контроля термоэлектрическими преобразователями (ТП). При измерениях ИЗ устанавливают контактной стороной на поверхность исследуемого изделия. Тепловое воздействие осуществляется с помощью нагревателя (Н). Мощность и длительность теплового воздействия БП задаются программно через интерфейс (И), контроллер К1, цифро-аналоговый преобразователь (ЦАП).

Сигналы с ТП поступают через мультиплексор (П), усилитель (У), аналого-цифровой преобразователь (АЦП), буфер обмена (Б) и интерфейс в ПК. Контроллер К2 обеспечивает необходимый порядок опроса каналов и различные диапазоны измерения на каждом из них. Сбор информации производится при нагреве и остывании исследуемого тела. В качестве ПК в ИС используется IВМ-совместимый компьютер. Программное обеспечение ИС включает системное (СПО), прикладное (ППО) и вспомогательное (ВПО) обеспечения.

Зонд состоит из двух основных узлов: измерительной ячейки 1 и корпуса 2. Измерительная ячейка, в свою очередь, состоит из основания 3, разъема 9. С контактной стороны измерительной ячейки на поверхности теплоизолятора 4 размещены микротермопары 5, сваренные встык. Линейный нихромовый нагреватель 6 изготовлен в виде узкой полосы. Выводы термопар и нагревателя проходят через специальные отверстия. В качестве микротермопар 5 используются дифференциальные термопары, горячие спаи которых устанавливаются в плоскости контакта ИЗ с исследуемым объектом. Холодные спаи микротермопар устанавливаем на подложку ИЗ. На крышке 11 измерительной ячейки расположен разъем 9, предназначенный для коммутации сенсора со встраиваемой в него платой PCI-1202H. Между измерительной ячейкой и крышкой 10 корпуса расположена пружина 8, которая обеспечивает одинаковое усилие прижима измерительной ячейки к поверхности исследуемого материала, что обуславливает равенство контактных тепловых сопротивлений при каждом измерении. Разностная ЭДС, полученная на зажимах микротермопар 5, поступает в плату PCI-1202H.

На рис. 2, в качестве примера, представлена конструкция ИЗ с линейным нагревателем.

Рис. 1. Структурная схема измерительной системы ИС с линейным нагревателем использовалась при экспериментах на материалах: полиметилметакрилат (ПММА), политетрафторэтилен (ПТФЭ), Капролон-В. Градуировку ИС проводили на ПММА. Результаты определения теплопроводности () и темратуропроводности (a) указанных материалов представлены в таблице 1. Для сравнения в таблеце представлены значения теплопроводности ( ИТ-3), полученные стационарным методом на приборе ИТ-3, и справочные данные ( спр., aспр).

Рис. 2. Схема измерительного зонда с линейным нагревателем Таблица 1. – Результаты определения теплопроводности ИС реализует метод неразрушающего контроля ТФС изделий из твердых неметаллических материалов [1, 2] и обеспечивает повышение точности, достоверности результатов и оперативность измерений за счет использования рабочих участков экспериментальных термограмм и самоконтроля ИС при каждой реализации эксперимента.

Работа выполнена при финансовой поддержке Фонда содействия развитию малых форм предприятий в научно-технической сфере по Программе «У.М.Н.И.К.», проект № 13962.

1. Патент 2161301 РФ, G 01 N 25/18. Способ неразрушающего определения теплофизических свойств материалов / Жуков Н.П., Муромцев Ю.Л., Майникова Н.Ф., Рогов И.В., Балашов А.А. № 99104568, заявл. 03.03.1999; опубл. 27.12.2000, 2. Жуков Н.П. Многомодельные методы и средства неразрушающего контроля теплофизических свойств твердых материалов и изделий / Н.П. Жуков, Н.Ф. Майникова // Монография. – М.: Машиностроение-1, 2004. – 288 с.

МОДЕЛИРОВАНИЕ ТЕПЛОПЕРЕНОСА ПРИ

ТЕПЛОФИЗИЧЕСКИХ ИЗМЕРЕНИЯХ

Майникова Н.Ф., д.т.н., Рогов И.В., к.т.н., Пудовкина Е.В., Антонов А.О.

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

В результате точность и надежность определения ТФС во многом обуславливается тем, насколько адекватно математическая модель описывает тепловые процессы, проходящие при измерении [1].

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

Исследуемое тело представляет собой конструкцию, состоящую из двух слоев: первый – низкотеплопроводный с теплофизическими свойствами: теплопроводностью (1), температуропроводностью (а1), плотностью (1); второй – высокотеплопроводный с теплофизическими свойствами: теплопроводностью (2), температуропроводностью (а2), плотностью (2). Толщина первого слоя (относительно тонкой пластины, например, полимерного покрытия) – h1. Второй слой – полуограниченное металлическое тело.

В начальный момент времени на свободную поверхность пластины (при х = 0) воздействуют постоянным тепловым потоком q = const. Боковые поверхности тел имеют тепловую изоляцию (q = 0).

Рис. 1. Схема системы, состоящей из ограниченного и В математическом виде задача записывается следующим образом:

Здесь Т 1 ( х, ), Т 2 ( х, ) – соответственно температурные поля исследуемых покрытия и основания в точках с координатами х в моменты времени ; q – тепловой поток, подводимый к свободной поверхности исследуемой пластины в точке х = 0.

После ряда преобразований и упрощений при больших решение задачи (1)(7) имеет следующий вид:

Здесь 1, 2 – тепловые активности первого и второго тел, соответственно.

Данное решение может быть использовано для определения толщины полимерного покрытия (h1), если известны ТФС материалов двухслойных объектов.

Для определения оптимальных условий метода неразрушающего контроля (НК) толщины покрытий изделий воспользуемся численным моделированием температурных полей методом конечных элементов с помощью пакета программ COMSOL/FEMLAB.

Построены термограммы для системы ограниченного и полуограниченного тел с различными толщинами второго слоя (рис. 2) для объекта, состоящего из металлического основания (сталь3) с полимерным покрытием (политетрафторэтилен (ПТФЭ)).

На термограммах (2–6) определены рабочие участки, на которых второе тело можно считать полуограниченным. Рабочий участок характеризуется регуляризацией теплового режима в локальной области нагревателя и тепловой поток, проходящий через точку контроля, становится практически постоянным. Это позволяет использовать математическую модель (8) описания процесса теплопереноса для рабочего участка термограммы.

Зная границы рабочих участков термограмм для различных h2, найдена зависимость f (h2 ) (рис. 3). По оси ординат на рисунке 3 представлены значения, соответствующие времени окончания рабочего участка, на термограммах (2–6) рисунка 2.

Зависимость, представленная на рисунке 4 для конкретного двухслойного объекта, позволяет определить время эксперимента. Например, при толщине второго слоя 20 мм модель (8) описывает теплоперенос в объекте исследования до момента времени = 15 с с начала эксперимента.

Рис. 2. Термограммы, построенные для системы ограниченного и 1 кривая получена по выражению (8); 26 кривые получены численным моделированием по методу конечных элементов при толщине второго слоя h Рис. 3. Зависимость = f (h2) 1. Жуков Н.П. Многомодельные методы и средства неразрушающего контроля теплофизических свойств твердых материалов и изделий / Н.П. Жуков, Н.Ф. Майникова // Монография. – М.: Машиностроение-1, 2004. – 288 с.

МОДЕЛИРОВАНИЕ ТЕПЛОПЕРЕНОСА ПРИМЕНИТЕЛЬНО К

МЕТОДУ ТЕПЛОФИЗИЧЕСКОГО ЭКСПЕРИМЕНТА

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

На рисунке 1 представлена измерительная схема метода. Тела 1 и характеризуются различными ТФС (теплопроводность и температуропроводность: 1, а1 и 2, а2). Размеры первого тела (исследуемый ПМ):

длина – l1; толщина – d1; ширина – 2H1. Размеры второго тела: длина – l2; толщина – d2; ширина – 2H2. В плоскости контакта тел действует линейный источник тепла. Второе тело – подложка измерительного зонда (ИЗ). Начальная температура тел одинакова и равна Т0. Мощность нагревателя – q 0. В ходе эксперимента фиксируется термограмма – зависимость разности значений температуры (Т) в точках контроля от времени (). Если на поверхности полуограниченного в тепловом отношении тела расположен нагреватель в виде полосы шириной 2h, через которую в материал подается тепловой поток q0, то температура T в точке с координатами (x, 0) поверхности (рис. 1) в момент времени [1]:

где – время, с; х – координата, м.

Рис. 1. Измерительная схема метода с линейным источником тепла:

1 – изделие; 2 – измерительный зонд; 3 – нагреватель; 4 – термоприемник После ряда упрощений в безразмерной форме:

где Fo a h 2 – время в безразмерном представлении; x x h – безразмерная координата; Т ( x, 0, ) q0 h – температура, в безразмерном представлении.

Для точки с координатами (0; 0), выражение (2) примет вид:

После ряда преобразований получено выражение для температуры поверхности z = 0 в точке с координатой x для больших значений [2]:

где – число Эйлера.

Выражение (4) в безразмерной форме применительно к точке с координатами (0; 0) для больших значений примет вид Таким образом, получено безразмерное представление температурной функции, которое имеет преимущества по сравнению с первоначальной (1). Действительно, вместо семи размерных величин в данном случае оказываются только две – безразмерная температура () и число Фурье (Fo). Число Фурье заняло место времени и может рассматриваться как безразмерное время.

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

На рисунке 2 представлены термограммы = f(Fo). Кривая 1 получена численным решением при соотношении d 2 h 5 методом конечных элементов. Кривые 2 и 3 получены по выражениям (3) и (5) с помощью системы аналитических вычислений Maple.

На термограмме выделены три участка – I, II и III.

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



Pages:   || 2 |
Похожие работы:

«Министерство образования и наук и Российской Федерации Санкт-Петербургский государственный политехнический университет Неделя Науки СПбГПу Материалы научно-практической конференции с международным участием 2–7 декабря 2013 года НаучНо-образовательНый цеНтр возобНовляемые виды эНергии и устаНовки На их осНове Санкт-Петербург•2014 УДК 621.31:627:502.63 ББК 31.6:31.15; 38.77 Н 42 Неделя науки СПбГПУ : материалы научно-практической конференции c международным участием. Научно-образовательный центр...»

«V ВСЕРОССИЙСКАЯ НАУЧНО – ПРАКТИЧЕСКАЯ КОНФЕРЕНЦИЯ Инновационные технологии в обучении и производстве Камышин 4-6 декабря 2008 г. МАТЕРИАЛЫ КОНФЕРЕНЦИИ Том 3 Вузы и организации, участвующие в конференции 1. Волгоградский государственный технический университет 2. Волжский политехнический институт (филиал) Волгоградского государственного технического университета 3. Камышинский технологический институт (филиал) Волгоградского государственного технического университета 4. Волгоградский...»

«РЕЦЕНЗИИ обсуждениях: Глобальное управление и безопасность: коллективная безопасность в Европе и Энергетическая безопасность: диалог Востока и Запада, за которыми последовали заседания рабочих групп, рассматривавших соответствующие вопросы в интерактивном режиме. Второй день был отмечен пленарными обсуждениями по темам Инвестиции и развивающиеся рынки: модели развития рынков и экономик в период финансовой нестабильности и Корпоративное управление: эффективные стратегии во времена глобальных...»

«Жизнин Станислав Захарович д.экон.н. Кафедра международных проблем ТЭК, профессор Доктор экономических наук, профессор кафедры международных проблем ТЭК МИЭП МГИМО (У) МИД России/ Работает на кафедре международных проблем ТЭК с сентября 2002 г. В 1969 г. окончил Харьковский авиационный институт по специальности инженерэлектрик. В 1977 г. - Дипломатическую академию МИД СССР по специальности международные экономические отношения. В 1998 г. защитил кандидатскую диссертацию Энергетическая...»

«Министерство образования и наук и Российской Федерации Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования ИВАНОВСКИЙ ГОСУДАРСТВЕННЫЙ ЭНЕРГЕТИЧЕСКИЙ УНИВЕРСИТЕТ ИМЕНИ В.И.ЛЕНИНА _ СОВРЕМЕННАЯ ЭКОНОМИЧЕСКАЯ НАУКА VIII МЕЖДУНАРОДНАЯ НАУЧНО-ТЕХНИЧЕСКАЯ КОНФЕРЕНЦИЯ СТУДЕНТОВ, АСПИРАНТОВ И МОЛОДЫХ УЧЕНЫХ ЭНЕРГИЯ – 2013 ИВАНОВО, 23-25 апреля 2013 г. МАТЕРИАЛЫ КОНФЕРЕНЦИИ ТОМ 6 _ ИВАНОВО ИГЭУ УДК 330. ББК 65. СОВРЕМЕННАЯ ЭКОНОМИЧЕСКАЯ НАУКА //...»

«Министерство образования Российской Федерации Ухтинский государственный технический университет Институт социально-экономических и энергетических проблем Севера Коми научного центра УрО РАН Сыктывкарский государственный университет Институт управления, информации и бизнеса Научно-исследовательский и проектный институт ПечорНИПИнефть ООО ВНИИгаз – филиал СеверНИПИгаз Межрегиональная научно-практическая конференция ПРОБЛЕМЫ ЭФФЕКТИВНОГО ОСВОЕНИЯ ПРИРОДНЫХ РЕСУРСОВ В УСЛОВИЯХ РЫНКА (29–30 октября...»

«Министерство образования и наук и Российской Федерации Тульский государственный университет Администрация Тульской области Академия горных наук Российская академия архитектуры и строительных наук Международная академия наук экологии и безопасности жизнедеятельности Совет молодых ученых Тульского государственного университета Международная научно-практическая конференция молодых ученых и студентов ОПЫТ ПРОШЛОГО – ВЗГЛЯД В БУДУЩЕЕ Конференция посвящена 150-летию со дня рождения учёного,...»

«ни-' ‘ in ± ь -Q > X НX S шу - mо нх оs Q. d >s ТЕХНОЛОГИЯ И АВТОМАТИЗАЦИЯ оы оо ш АТОМНОЙ ЭНЕРГЕТИКИ т S >5: 1_ sо п; ОО Q. ШX ШX Шш Он Материалы отраслевой научно-технической конференции 12-14 мая 2004г. ьо МИНИСТЕРСТВО РОССИЙСКОЙ ФЕДЕРАЦИИ ПО АТОМНОЙ ЭНЕРГИИ АДМИНИСТРАЦИЯ ЗАТО СЕВЕРСК СИБИРСКИЙ ХИМИЧЕСКИЙ КОМБИН АТ ТОМСКИЙ...»

«Министерство сельского хозяйства РФ Департамент научно-технологической политики и образования Министерство сельского хозяйства Иркутской области Иркутская государственная сельскохозяйственная академия НАУЧНЫЕ ИССЛЕДОВАНИЯ СТУДЕНТОВ В РЕШЕНИИ АКТУАЛЬНЫХ ПРОБЛЕМ АПК Материалы студенческой научно-практической конференции с международным участием, посвященной 80-летию ФГБОУ ВПО ИрГСХА (19-20 марта 2014 г., г. Иркутск) Часть II Иркутск, 2014 1 УДК 001:63 ББК 40 Н 347 Научные исследования студентов в...»

«МЕЖДУНАРОДНАЯ НАУЧНО-ПРАКТИЧЕСКАЯ КОНФЕРЕНЦИЯ УРАЛЬСКАЯ ГОРНАЯ ШКОЛА – РЕГИОНАМ 11-12 апреля 2011 г. БЕЗОПАСНОСТЬ ЖИЗНЕДЕЯТЕЛЬНОСТИ УДК 621.039 ВЫШЕДШИЕ ИЗ УПОТРЕБЛЕНИЯ ПЕСТИЦИДЫ КАК УГРОЗА ЭКОЛОГИЧЕСКОЙ БЕЗОПАСНОСТИ ТЕРРИТОРИЙ И ЗДОРОВЬЮ ЧЕЛОВЕКА БОЛТЫРОВ В. Б. ГОУ ВПО Уральский государственный горный университет Проблема пестицидов в общем и СОЗ в частности особенно актуальна для России и связана с развитым сельскохозяйственным производством, высоким удельным весом энергетического и...»

«Министерство образования и наук и Российской Федерации Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования Тюменский государственный нефтегазовый университет ЭНЕРГОСБЕРЕЖЕНИЕ И ИННОВАЦИОННЫЕ ТЕХНОЛОГИИ В ТОПЛИВНО-ЭНЕРГЕТИЧЕСКОМ КОМПЛЕКСЕ Материалы региональной научно-практической конференции студентов, аспирантов, молодых учёных и специалистов Тюмень ТюмГНГУ 2012 УДК 338.45 (06)+656.5(06) ББК 65.301 Э653 Редакционная коллегия: А. Л. Портнягин...»

«РОССИЙСКАЯ АКАДЕМИЯ НАУК КАРЕЛЬСКИЙ НАУЧНЫЙ ЦЕНТР ИНСТИТУТ ГЕОЛОГИИ НАУЧНО-ПРОИЗВОДСТВЕННЫЙ КОМПЛЕКС ООО НПК КАРБОН-ШУНГИТ ПРОМЫШЛЕННО-СТРОИТЕЛЬНАЯ КОМПАНИЯ ООО АЛЬФА-ПОЛ ШУНГИТЫ И БЕЗОПАСНОСТЬ ЖИЗНЕДЕЯТЕЛЬНОСТИ ЧЕЛОВЕКА МАТЕРИАЛЫ ПЕРВОЙ ВСЕРОССИЙСКОЙ НАУЧНО-ПРАКТИЧЕСКОЙ КОНФЕРЕНЦИИ (3-5 октября 2006) Под редакцией д.т.н. Ю.К.Калинина Петрозаводск 2007 УДК Шунгиты и безопасность жизнедеятельности человека. Материалы Первой Всероссийской научнопрактической конференции. Петрозаводск:, 2007....»









 
2014 www.konferenciya.seluk.ru - «Бесплатная электронная библиотека - Конференции, лекции»

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