<?xml version="1.0" encoding="UTF-8"?>
<rss version="2.0"
	xmlns:content="http://purl.org/rss/1.0/modules/content/"
	xmlns:wfw="http://wellformedweb.org/CommentAPI/"
	xmlns:dc="http://purl.org/dc/elements/1.1/"
	xmlns:atom="http://www.w3.org/2005/Atom"
	xmlns:sy="http://purl.org/rss/1.0/modules/syndication/"
	xmlns:slash="http://purl.org/rss/1.0/modules/slash/"
	>

<channel>
	<title>Электронный научно-практический журнал «Современные научные исследования и инновации» &#187; фазовые переходы</title>
	<atom:link href="http://web.snauka.ru/issues/tag/fazovyie-perehodyi/feed" rel="self" type="application/rss+xml" />
	<link>https://web.snauka.ru</link>
	<description></description>
	<lastBuildDate>Sat, 18 Apr 2026 09:41:14 +0000</lastBuildDate>
	<language>ru</language>
	<sy:updatePeriod>hourly</sy:updatePeriod>
	<sy:updateFrequency>1</sy:updateFrequency>
	<generator>http://wordpress.org/?v=3.2.1</generator>
		<item>
		<title>Пирометрические исследования структурообразования в процессе стационарного нагрева вольфрама</title>
		<link>https://web.snauka.ru/issues/2016/12/76197</link>
		<comments>https://web.snauka.ru/issues/2016/12/76197#comments</comments>
		<pubDate>Thu, 29 Dec 2016 18:20:51 +0000</pubDate>
		<dc:creator>Долматов Алексей Викторович</dc:creator>
				<category><![CDATA[05.00.00 ТЕХНИЧЕСКИЕ НАУКИ]]></category>
		<category><![CDATA[argon]]></category>
		<category><![CDATA[brightness]]></category>
		<category><![CDATA[diagnostic]]></category>
		<category><![CDATA[heat]]></category>
		<category><![CDATA[phase transition]]></category>
		<category><![CDATA[pyrometry]]></category>
		<category><![CDATA[spectral]]></category>
		<category><![CDATA[structure]]></category>
		<category><![CDATA[temperature]]></category>
		<category><![CDATA[tungsten]]></category>
		<category><![CDATA[вольфрам]]></category>
		<category><![CDATA[диагностика]]></category>
		<category><![CDATA[нагрев]]></category>
		<category><![CDATA[пирометрия]]></category>
		<category><![CDATA[спектральная]]></category>
		<category><![CDATA[среда аргона]]></category>
		<category><![CDATA[структурообразование]]></category>
		<category><![CDATA[температура]]></category>
		<category><![CDATA[фазовые переходы]]></category>
		<category><![CDATA[яркостная]]></category>

		<guid isPermaLink="false">https://web.snauka.ru/?p=76197</guid>
		<description><![CDATA[1. Введение Формирование нового производственного уклада во многом связано с развитием аддитивных технологий на базе металлов и металлокерамик. В них синтез материалов и изделий осуществляется из мелкодисперсных компонентов в ходе быстропротекающих  высокотемпературных процессов. Здесь актуальность контроля структурных и фазовых изменений обусловлена широким спектром научных и технических задач [1-6]. Высокая температура аддитивного синтеза металлов и металлокерамик [...]]]></description>
			<content:encoded><![CDATA[<p><strong>1. Введение</strong></p>
<p style="text-align: justify;">Формирование нового производственного уклада во многом связано с развитием аддитивных технологий на базе металлов и металлокерамик. В них синтез материалов и изделий осуществляется из мелкодисперсных компонентов в ходе быстропротекающих  высокотемпературных процессов. Здесь актуальность контроля структурных и фазовых изменений обусловлена широким спектром научных и технических задач [1-6].</p>
<p style="text-align: justify;">Высокая температура аддитивного синтеза металлов и металлокерамик (1500 – 4000 К) открывает возможность косвенного контроля структурообразования  методами оптической пирометрии с разрешением до 1 мкм и 10 мкс [7-9]. Яркостная температура линейно зависит от спектрального коэффициента излучения, который в свою очередь определяется микроструктурой и фазовым составом поверхностного слоя наблюдаемого тела. Погрешность современных средств яркостной пирометрии в диапазоне температур 1500 – 5000 К достигает 0.05 % [10, 11]. Метод спектральной пирометрии, напротив, позволяет выбрать такой оптический диапазон, где измеряемая температура практически не зависит от свойств материала, а отклонение спектральной температуры от действительной не превышает 0.5 %. Случайная погрешность цифровых спектральных пирометров на 1-2 порядка ниже систематической составляющей [12]. Таким образом, сочетание в измерительном комплексе средств яркостной и спектральной пирометрии теоретически дает возможность контроля излучательной способности неподвижных объектов с высокой точностью.</p>
<p style="text-align: justify;">Принципиально система оптической диагностики структурообразования может быть построена на основе результатов пирометрических и рентгенографических исследований. Первые позволяют определить зависимость от температуры спектрального коэффициента излучения материала, а вторые изучить поведение его структурно-фазового состава. Цель настоящей работы заключалась в обнаружении физических явлений, влияющих на измеряемую величину спектральной излучательной  способности в процессе синтеза материалов, и оценке точности подхода пирометрической диагностики структурно-фазовых превращений.</p>
<p><strong>2. Эксперимент</strong></p>
<p>Экспериментальные исследования выполнялись с помощью оригинального микропирометрического комплекса (рис. 1).</p>
<p style="text-align: center;"><a href="https://web.snauka.ru/issues/2016/12/76197/fig-1-2" rel="attachment wp-att-76200"><img class="aligncenter size-full wp-image-76200" src="https://web.snauka.ru/wp-content/uploads/2016/12/Fig-1.bmp" alt="Рис. 1. Экспериментальная установка (а) и крепление образца в реакторе высокотемпературного нагрева (б) (1 – микроскоп МБС-10, 2 – полосовой светофильтр в тепловизионном канале, 3 –  камера ВидеоСпринт, 4 – спектрометр LR1-T, 5 – источник тока PSH-2035, 6 – реактор высокотемпературного нагрева (РВТН), 7 – вакуумный насос, 8 – охлаждаемые проточной водой  медные электроды, 9 – закрепленный образец фольги, 10 – область визирования спектрометра)" width="756" height="466" /></a></p>
<p style="text-align: center;">Рис. 1. Экспериментальная установка (а) и крепление образца в реакторе высокотемпературного нагрева (б) (1 – микроскоп МБС-10, 2 – полосовой светофильтр в тепловизионном канале, 3 –  камера ВидеоСпринт, 4 – спектрометр LR1-T, 5 – источник тока PSH-2035, 6 – реактор высокотемпературного нагрева (РВТН), 7 – вакуумный насос, 8 – охлаждаемые проточной водой  медные электроды, 9 – закрепленный образец фольги, 10 – область визирования спектрометра)</p>
<p style="text-align: justify;">В нем микроскоп МБС-10 (ЛЗОС, Россия) (1) служит оптической системой, которая совмещает в пространстве два измерительных канала: тепловизионный и спектральный. В тепловизионном канале установлен полосовой светофильтр SL-725-40 (PhotoОptic-Filters, Россия) (2) с центральной длиной волны 725 нм и полушириной полосы пропускания 40 нм. Регистратором здесь служит цифровая камера ВидеоСпринт (НПК «Видеоскан», Россия) (3). Микроскоп обеспечивает пространственное разрешение тепловизионной съемки 2.9 мкм, а быстродействие камеры позволяет фиксировать кадры с частотой от 2 до 250000 Гц [13]. С помощью образцовой температурной лампы ТРУ-1200-2350 выполнена калибровка тепловизионного канала на доступном множестве значений времени экспозиции, и построены таблицы соответствия дискретного уровня сигнала яркостной температуре [14]. Погрешность измерения яркостной температуры оценена на уровне 0.1 %. Спектральный канал комплекса организован с помощью цифрового спектрометра LR1-T (ASEQ Instruments, Канада) (4), соединенного оптоволокном с микроскопом. Обработка данных теплового спектра и определение температуры выполнялись на компьютере в режиме реального времени [15]. Случайная погрешность измерения спектральной температуры оценена с помощью образцовой лампы на уровне 0.07%.</p>
<p style="text-align: justify;">Образцы для пирометрических исследований изготавливались из вольфрамовой фольги толщиной 100 мкм. Длина образца составляла 50 мм, ширина около  2 мм. Нагрев образцов выполнялся электрическим током от программируемого источника PSH-2035 (Good Will Instrument, Тайвань) (5) в оригинальном реакторе высокотемпературного нагрева (РВТН) (6). РВТН имеет цилиндрическую кварцевую стенку, которая позволяет удалять воздух с помощью вакуумного насоса (7), замещать его инертной средой (аргон) и контролировать излучение с поверхности образца.  Медные электроды (8) РВТН во время опыта охлаждались проточной водой. Область визирования спектрометра (10) имела форму круга с диаметром 230 мкм и размещалась в центре тепловизионного кадра размерами 3700 х 2970 мкм.</p>
<p style="text-align: justify;">Экспериментальная зависимость коэффициента излучения от температуры на длине волны 725 нм строилась для общей области визирования измерительных каналов комплекса. При этом допускалось, что во время наблюдения в пределах данной области излучательная способность однородна. Дополнительно объективными и субъективными средствами детектировалось наличие процесса структурообразования в образце, и производилась его тепловизионная съемка.</p>
<p style="text-align: justify;">Методика эксперимента включала следующие шаги. С помощью источника PSH-2035 с шагом 0.2 А изменялся ток. Мощность тепловыделения в образце определялась по фактическим данным тока и напряжения, абсолютная погрешность измерения которых составляла 0.01 ампера и вольта соответственно. После установки нового уровня тока экспериментатор отслеживал состояние поверхности образца, а напряжение на выходе источника контролировалось автоматически. Если через 5 секунд после смены уровня тока изменение мощности источника продолжалось или наблюдались структурные изменения на поверхности фольги, то запускалась длительная тепловизионная съемка с частотой от 2 до 50 Гц. В противном случае считалось, что переходные процессы в образце завершены. Тогда  в общей области визирования измерялась спектральная температура, а по данным тепловизионной съемки оценивалась яркостная температура. Спектральный коэффициент излучения рассчитывался по формуле:</p>
<p style="text-align: center;"><a href="https://web.snauka.ru/issues/2016/12/76197/formula-1-18" rel="attachment wp-att-76203"><img class="aligncenter size-full wp-image-76203" src="https://web.snauka.ru/wp-content/uploads/2016/12/formula-1.bmp" alt="Формула определения излучательной способности в эксперименте" width="372" height="56" /></a></p>
<p style="text-align: left;"><a href="https://web.snauka.ru/issues/2016/12/76197/formula-1-note" rel="attachment wp-att-76204"><img class="size-full wp-image-76204" src="https://web.snauka.ru/wp-content/uploads/2016/12/formula-1note.bmp" alt="Примечание к формуле (1)" width="683" height="44" /></a></p>
<p>Для изготовления экспериментальных образцов использовался прокат вольфрама, химический состав которого представлен в таблице 1.</p>
<p style="text-align: center;" align="right">Таблица 1. Химический состав вольфрамовой фольги</p>
<div align="center">
<table border="1" cellspacing="0" cellpadding="0">
<tbody>
<tr>
<td rowspan="2" width="73">
<p align="center">W, %</p>
</td>
<td colspan="7" valign="top" width="427">Содержание примесей, не более мкг/г</td>
</tr>
<tr>
<td width="47">
<p align="center">Fe</p>
</td>
<td width="59">
<p align="center">Al</p>
</td>
<td width="59">
<p align="center">Ni</p>
</td>
<td width="59">
<p align="center">Si</p>
</td>
<td width="59">
<p align="center">C</p>
</td>
<td width="70">
<p align="center">O</p>
</td>
<td width="74">
<p align="center">N</p>
</td>
</tr>
<tr>
<td valign="top" width="73">
<p align="center">99,95</p>
</td>
<td valign="top" width="47">
<p align="center">80</p>
</td>
<td valign="top" width="59">
<p align="center">50</p>
</td>
<td valign="top" width="59">
<p align="center">50</p>
</td>
<td valign="top" width="59">
<p align="center">50</p>
</td>
<td valign="top" width="59">
<p align="center">100</p>
</td>
<td valign="top" width="70">
<p align="center">100</p>
</td>
<td valign="top" width="74">
<p align="center">100</p>
</td>
</tr>
</tbody>
</table>
</div>
<p>Создание аргоновой среды выполнялось в ходе четырех итераций процесса, включающего откачку газа из РВТН до давления 2 Па и заполнения ее аргоном при давлении 400 кПа. На последней итерации давление аргона устанавливалось на заданном уровне 100 кПа. В эксперименте РВТН не отключался от магистрали подачи аргона, что способствовало стабилизации давления во время нагрева образца вольфрама. Таким образом, содержание примесей кислорода и азота в среде аргона перед началом эксперимента оценивалось величиной &lt;10<sup>-8 </sup>мкг/г. Так  как концентрация данных элементов в самом образце на 8-10 порядков выше, то процессы химического взаимодействия образца с газовой средой не рассматривались. Более того, при анализе структурно-фазовых изменений считалось, что атомы примесей или их химические соединения, покинувшие поверхность образца, обратно не возвращается.</p>
<p style="text-align: justify;">В ходе опыта ступенчатый нагрев образца продолжался либо до достижения верхнего предела генерации источника тока или сигнала измерительных приборов, либо прерывался специально. Во всех случаях выполнялось мгновенное отключение тока, что вело к закалке образца с сохранением микроструктуры и фазового состава последнего этапа нагрева. Все образцы подвергались только однократному нагреву. По окончанию пирометрических исследований над образцами произведен микроскопический, рентгенофазовый и микроэлементный анализ.</p>
<p><strong>3. Обсуждение результатов</strong></p>
<p>В результате опытов с давлением аргоновой среды 100 кПа по данным тепловизионной и спектральной съемки и формуле (1), получены зависимости, представленные на рисунке 2.</p>
<p style="text-align: center;"><a href="https://web.snauka.ru/issues/2016/12/76197/fig-2" rel="attachment wp-att-76205"><img class="aligncenter size-full wp-image-76205" src="https://web.snauka.ru/wp-content/uploads/2016/12/Fig-2.bmp" alt="Результаты пирометрических исследований вольфрамовых образцов в области с максимумом температурного поля при давлении аргоновой среды 100 кПа" width="700" height="419" /></a>Рис. 2. Результаты пирометрических исследований вольфрамовых образцов в области с максимумом температурного поля при давлении аргоновой среды 100 кПа</p>
<p style="text-align: justify;">На зависимости спектральной излучательной способности (спектральный коэффициент излучения) всех образцов отчетливо видны перепады, по которым были определены температуры фазовых превращений.  Кроме того, микротепловизионная съемка зафиксировала изменения структуры материала на поверхности образцов в области температур близких к фазовому переходу (рис. 3).</p>
<p style="text-align: center;"><a href="https://web.snauka.ru/issues/2016/12/76197/fig-3-2" rel="attachment wp-att-76206"><img class="aligncenter size-full wp-image-76206" src="https://web.snauka.ru/wp-content/uploads/2016/12/Fig-3.bmp" alt="Микротепловизионные изображения вольфрамового образца с разрешением 5.9 мкм" width="749" height="541" /></a></p>
<p style="text-align: center;">Рис. 3. Микротепловизионные изображения вольфрамового образца с разрешением 5.9 мкм</p>
<p>Для интерпретации экспериментальных зависимостей были привлечены справочные данные (в том числе таблица 2) и диаграмма фазовых состояний системы &#8220;вольфрам – кислород&#8221; (рис. 4)[16-21].</p>
<p style="text-align: center;"><span style="text-align: justify;">Таблица 2. Физические свойства оксидов вольфрама</span></p>
<p style="text-align: justify;"><a href="https://web.snauka.ru/issues/2016/12/76197/tabl-2-5" rel="attachment wp-att-76207"><img class="aligncenter size-full wp-image-76207" src="https://web.snauka.ru/wp-content/uploads/2016/12/Tabl-2.bmp" alt="Справочные данные по системе &quot;Вольфрам-Кислород&quot;" width="426" height="196" /></a></p>
<p style="text-align: justify;">В итоге возникло следующее объяснение результатов опытов. Во время эксперимента температура на концах образца за счет водяного отвода тепла от электродов не превышала 400-500 К. Максимум температурного поля образцов наблюдался примерно на их середине. В отсутствие структурно-фазовых изменений данное поле являлось стационарным. Причем с весьма высоким градиентом до 10<sup>5</sup> К/м. Таким образом, были созданы хорошие условия для диффузии углерода, азота и кислорода (и других примесей) к поверхности кристаллитов вольфрама и в высокотемпературную область образца, где их растворимость больше. Следствием подобной диффузии является повышенная концентрация  примесей на межзеренных границах с образованием там оксидов, нитридов и карбидов вольфрама (и других соединений).</p>
<p style="text-align: center;"><a href="https://web.snauka.ru/issues/2016/12/76197/fig-4-2" rel="attachment wp-att-76208"><img class="aligncenter size-full wp-image-76208" src="https://web.snauka.ru/wp-content/uploads/2016/12/Fig-4.bmp" alt="Диаграмма фазовых состояний системы &quot;вольфрам-кислород&quot;" width="382" height="364" /></a></p>
<p style="text-align: center;">Рис. 4. Диаграмма фазовых состояний системы &#8220;вольфрам-кислород&#8221;</p>
<p style="text-align: justify;">Взаимодействие вольфрама с кислородом ведет к образованию двух основных фазы: WO<sub>2</sub> и WO<sub>3</sub>. При температуре выше 1000 К обе фазы начинают заметно сублимировать. Однако газовая фаза WO<sub>2 </sub>распадается на W и WO<sub>3</sub>. В диапазоне температур 1400 &#8211; 1550 К в кристаллической структуре WO<sub>2</sub> постепенно нарушается дальний порядок, что может выражаться в пластификации оксида и высвобождению кислорода как в чистом виде, так и в виде соединения с углеродом. Микротепловизионные наблюдения и &#8220;заморозка&#8221; образцов показали, что при температурах 1350 – 1550 К на поверхности вольфрамовых образцов возникают сферические и полусферические образования, подобные пузырям (рис. 5). Только парциальное давление насыщенных паров WO<sub>3</sub> на температуре около 1500 К по справочным данным составляет почти 50 кПа (таблица 2). Поэтому при температурах чуть меньше 1550 К в реакторе, наполненном аргоном с давлением 100 кПа, создаются неплохие условия для возникновения пузырей, наполненных паром WO<sub>3</sub>+O<sub>2</sub>+CO<sub>x</sub> с пластичной оболочкой из WO<sub>2</sub>.</p>
<p style="text-align: center;"><a href="https://web.snauka.ru/issues/2016/12/76197/fig-5" rel="attachment wp-att-76209"><img class="aligncenter size-full wp-image-76209" src="https://web.snauka.ru/wp-content/uploads/2016/12/Fig-5.bmp" alt="Полусферические образования на поверхности вольфрамовых образцов" width="397" height="354" /></a></p>
<p style="text-align: center;">Рис. 5. Полусферические образования на поверхности вольфрамовых образцов</p>
<p style="text-align: justify;">Образование подобных пузырей согласуется с поведением излучательной способности экспериментальных образцов в диапазоне температур 1350-1500 К, где рост коэффициента излучения происходит  по логарифмическому закону от значения 0,4 до 1 (рис. 2). На образцах фольги, закаленных при максимуме температуры 1500-1750 К, обнаружено, что в областях, где локальная температура была выше, радиус кривизны пузырей и их количество на единице поверхности больше (рис. 5). Каждый пузырь может работать как полость, в которой происходит многократное отражение теплового излучения. Тоже самое можно сказать про полости между пузырями. Излучательная способность любой полости стремиться к единице с ростом количества отражений в ней. Таким образом, увеличение количества пузырей на единице поверхности должно вести к росту ее излучательной способности, что и показывают экспериментальные данные, в которых спектральный коэффициент излучения определялся как среднее значение по площадке диаметром 230 мкм.</p>
<p style="text-align: justify;">При достижении температуры плавления WO<sub>2</sub> (~1500 К) на экспериментальной кривой излучательной способности наблюдается спад со значения 1 до 0,48 (рис. 2). Его можно объяснить разрушением структуры указанного оксида и исчезновением пузырей с поверхности образца. Причем на границах кристаллитов возможно резкое ускорение процесса восстановления вольфрама углеродом. Также следует отметить, что снижение коэффициента излучение не повлекло в эксперименте изменение электрической мощности, подводимой к образцу (т.е. электрическое сопротивление образцов не изменилось), но вызвало резкий рост (на 450-470 градусов) спектральной (термодинамической) температуры (рис. 2), которая практически достигала уровня закипания оксидов вольфрама. Таким образом, фактически плавление WO<sub>2</sub> инициировало структурно-фазовые превращения, в результате которых с поверхности образца исчезли оксиды, а ее излучательная способность с погрешность 4% стала равна справочному значению для чистого вольфрама. Высокую вероятность реализации подобного сценария подтверждает рентгенофазовый анализ образцов, которые в эксперименте находились при температуре выше 2000 К. На их дифрактограммах обнаруживаются только линии нитрида вольфрама с гексагональной решеткой, но отсутствуют  явные признаки оксидов и карбидов (рис. 6).</p>
<p style="text-align: center;"><a href="https://web.snauka.ru/issues/2016/12/76197/fig-6" rel="attachment wp-att-76210"><img class="aligncenter size-full wp-image-76210" src="https://web.snauka.ru/wp-content/uploads/2016/12/Fig-6.bmp" alt="Дифрактограмма вольфрамового образца после эксперимента с PAr=100 кПа и Tmin=2050К" width="571" height="373" /></a></p>
<p style="text-align: center;">Рис. 6. Дифрактограмма вольфрамового образца после эксперимента с P<sub>Ar</sub>=100 кПа и T<sub>min</sub>=2050К</p>
<p style="text-align: justify;">На графиках спектральной и яркостной температуры (рис. 2) при наблюдается небольшой излом. Однако излучательная способности образцов здесь остается неизменной. Данная температура хорошо совпадает с температурой диссоциации нитрида вольфрама. Причем азот может оставаться в вольфраме в виде твердого раствора внедрения вплоть до температуры около 2400 К, когда его тепловая энергия сравнивается с работой выхода из металла. Таким образом,  после разрушения нитридов вольфрама концентрация азота в вольфраме останется на прежнем уровне, что объясняет постоянство излучательной способности. Но происходит увеличение теплоемкости, которое вызывает уменьшение температуры.</p>
<p style="text-align: justify;">В диапазоне температур 2150 &#8211; 2400 К наблюдается снижение излучательной способности нагретых образцов до значения 0,27. Затем вплоть до температуры 2450 К излучательная способность возвращается к уровню 0,48, а температура падает на 150 градусов (рис. 2). В основе такого поведения может лежать процесс диффузии азота в вольфраме после разрушения нитридных связей, приводящий к концентрации его атомов в межзеренных границах (в том числе на поверхности образца).  По достижении температурного предела существования твердого раствора азот улетучивается с поверхности вольфрама и происходит восстановление излучательной способности металла, которое влечет снижение температуры при условии постоянства мощности тепловыделения в образце.</p>
<p><strong>4. Выводы</strong></p>
<p style="text-align: justify;">Оптический контроль излучательной способности вольфрамовых образцов позволил: идентифицировать температуру плавления WO<sub>2</sub>, температуру кипения оксидов вольфрама, температуру диссоциации нитридов вольфрама и температуру, при которой атомы азота покидают вольфрам; выявить на поверхности материала возникновение и разрушение макроструктур (полусферические образования). Оригинальное объяснение механизмов структурообразования в экспериментальных образцах стоит воспринимать на уровне рабочей гипотезы, для повышения достоверности которой необходимо провести более тщательные структурные исследования. Однако выполненные эксперименты показали, что чувствительности средств яркостной и спектральной пирометрии достаточно, чтобы с погрешностью 1-3% оценивать температуры фазовых переходов и выявлять структурные изменения в веществе как на микро-, так и на макромасштабе.  С точки зрения диагностики фазовых превращений влияние макроструктур на коэффициент излучения является отрицательным моментом, так как оно вносит искажения или маскирует зависимость этого оптического свойства от микроструктуры материала. В таком случае помогает тепловизионная съемка, по которой можно обнаружить появление макроструктур на поверхности материала. Кроме того, большое увеличение в микротепловизионной съемке позволит контролировать локальные значения коэффициента излучения и использовать селективный подход для устранения влияния поверхностных макроструктур на пирометрическую диагностику фазовых переходов. В области температур фазовых переходов от 1500 до 6000 К пространственное разрешение оптических тепловизионных систем может составлять от 10 до 1 мкм соответственно.  Особенно важным здесь является возможность регистрировать динамику процесса с временным разрешением от 10 мс до 10 мкс. Это позволяет средствами пирометрической диагностики вести наблюдение структурообразования в быстропротекающих процессах.</p>
<p style="text-align: justify;">Работа выполнена при финансовой поддержке РФФИ (грант № 15-48-00100).</p>
]]></content:encoded>
			<wfw:commentRss>https://web.snauka.ru/issues/2016/12/76197/feed</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>Модели динамики гетероструктур снеголедовых масс</title>
		<link>https://web.snauka.ru/issues/2017/04/80911</link>
		<comments>https://web.snauka.ru/issues/2017/04/80911#comments</comments>
		<pubDate>Tue, 18 Apr 2017 14:07:33 +0000</pubDate>
		<dc:creator>Шорин Владимир Алексеевич</dc:creator>
				<category><![CDATA[05.00.00 ТЕХНИЧЕСКИЕ НАУКИ]]></category>
		<category><![CDATA[dynamics]]></category>
		<category><![CDATA[heterostructures]]></category>
		<category><![CDATA[models]]></category>
		<category><![CDATA[phase transitions]]></category>
		<category><![CDATA[snow-ice masses]]></category>
		<category><![CDATA[гетероструктуры]]></category>
		<category><![CDATA[динамика]]></category>
		<category><![CDATA[модели]]></category>
		<category><![CDATA[снеголедовые массы]]></category>
		<category><![CDATA[фазовые переходы]]></category>

		<guid isPermaLink="false">https://web.snauka.ru/issues/2017/04/80911</guid>
		<description><![CDATA[Введение Динамика гетерогенных структур с фазовыми переходами на границах играет важную роль в природных и искусственных системах [1]. Например, каждую зиму в г.Москве от сорвавшихся с крыш снеголедовых масс и сосулек страдает более 50 человек и до 300 автомобилей. Самолеты, вертолеты, морские суда, портовые сооружения, нефтегазовые, космические объекты и др. подвержены обледенению. Эти процессы чреваты [...]]]></description>
			<content:encoded><![CDATA[<p><strong><span>Введение</span></strong><br />
<span>Динамика гетерогенных структур с фазовыми переходами на границах играет важную роль в природных и искусственных системах [1].</span><br />
<span>Например, каждую зиму в г.Москве от сорвавшихся с крыш снеголедовых масс и сосулек страдает более 50 человек и до 300 автомобилей.</span><br />
<span>Самолеты, вертолеты, морские суда, портовые сооружения, нефтегазовые, космические объекты и др. подвержены обледенению. Эти процессы чреваты серьезными финансовыми потерями, гибелью людей.</span><br />
<span>Математическое моделирование динамики гетерогенных структур снеголедовых масс является актуальной задачей.</span><br />
<strong><span>1. Фундаментальная модель нестационарного теплообмена</span></strong><br />
<span>Модель описывается в рамках двухмерного нестационарного уравнения теплопроводности, граничные условия задают конвективный теплообмен с окружающей средой.</span><br />
<span>Рассматривается бесконечная структура прямоугольного сечения, составленная из разнородных материалов. Материалы различаются теплофизическими свойствами: теплопроводностью , удельной теплоемкостью </span><em><span>С</span></em><span>, коэффициентом теплообмена с окружающей средой </span><img src="http://content.snauka.ru/web/80911_files/0EXHGYS4E.gif" alt="" width="22" height="24" /><span>. На верхних и нижней гранях происходит конвективный теплообмен с окружающей средой, температура которой </span><img src="http://content.snauka.ru/web/80911_files/0ILECHRB1.gif" alt="" width="25" height="25" /><span>. Две другие грани теплоизолированы.</span><br />
<span>В начальный момент времени структура нагрета до температуры </span><img src="http://content.snauka.ru/web/80911_files/08HDV1OVU.gif" alt="" width="70" height="25" /><span>. Задача состоит в нахождении поля температур </span><img src="http://content.snauka.ru/web/80911_files/001G15Y3D.gif" alt="" width="76" height="24" /><span>.</span><br />
<span>В математическом отношении задача сводится к решению нестационарного уравнения теплопроводности в двухмерной области </span><em><span>G</span></em><span> c соответствующими краевыми и начальными условиями:</span><br />
<img src="http://content.snauka.ru/web/80911_files/0UAFA4Y79.gif" alt="" width="304" height="44" /><span>;</span><br />
<img src="http://content.snauka.ru/web/80911_files/0JV3RQBA1.gif" alt="" width="181" height="24" /><span>;</span><br />
<img src="http://content.snauka.ru/web/80911_files/0PGO9VRNK.gif" alt="" width="121" height="25" /><span>, </span><img src="http://content.snauka.ru/web/80911_files/0H1IJ78TE.gif" alt="" width="52" height="25" /><span>, </span><img src="http://content.snauka.ru/web/80911_files/0N23UDRRC.gif" alt="" width="90" height="25" /><span>, </span><img src="http://content.snauka.ru/web/80911_files/0R3TLE00R.gif" alt="" width="38" height="20" /><span>.</span><br />
<span>Здесь </span><img src="http://content.snauka.ru/web/80911_files/05JHH0DT2.gif" alt="" width="165" height="29" /><span> тепловой поток; </span><img src="http://content.snauka.ru/web/80911_files/0A721DM12.gif" alt="" width="72" height="25" /><span> – коэффициент теплопроводности, </span><img src="http://content.snauka.ru/web/80911_files/0KMVQ2BX8.gif" alt="" width="84" height="24" /><span> – удельная теплоемкость; </span><img src="http://content.snauka.ru/web/80911_files/0EG8ZUMXH.gif" alt="" width="69" height="25" /><span> – температура в момент </span><img src="http://content.snauka.ru/web/80911_files/04ELCJV8E.gif" alt="" width="38" height="20" /><span>. Решение </span><img src="http://content.snauka.ru/web/80911_files/0Q3SDIH4X.gif" alt="" width="68" height="25" /><span> ищется в цилиндре </span><img src="http://content.snauka.ru/web/80911_files/0T097VHE7.gif" alt="" width="149" height="29" /><span>, основанием которого является прямоугольник </span><em><span>G</span></em><span> с границей </span><em><span>дG</span></em><span>.</span><br />
<strong><span>2. Конечно-разностная модель гетероструктуры</span></strong><br />
<span>Задача решается методом конечных разностей, в </span><em><span>G</span></em><span> вводится пространственная сетка</span><br />
<img src="http://content.snauka.ru/web/80911_files/0TUQG7GDO.gif" alt="" width="340" height="24" /><br />
<img src="http://content.snauka.ru/web/80911_files/0X7H9FSQ7.gif" alt="" width="330" height="28" /><br />
<img src="http://content.snauka.ru/web/80911_files/0XODC9NWN.gif" alt="" width="316" height="25" /><br />
<span>и сетка по времени</span><br />
<img src="http://content.snauka.ru/web/80911_files/1.gif" alt="" width="460" height="29" /><span>. </span><br />
<span>Здесь </span><img src="http://content.snauka.ru/web/80911_files/1(1).gif" alt="" width="113" height="25" /><span> – переменные шаги сетки по пространству в направлениях </span><em><span>Х</span></em><span> и </span><em><span>Y</span></em><span> и по времени, соответственно.</span><br />
<span>Задача решается на сетке </span><img src="http://content.snauka.ru/web/80911_files/1(2).gif" alt="" width="98" height="25" /><span>, вводится сеточная функция температуры </span><img src="http://content.snauka.ru/web/80911_files/1(3).gif" alt="" width="125" height="32" /><span>, определенная на </span><img src="http://content.snauka.ru/web/80911_files/1(4).gif" alt="" width="30" height="25" /><span>. Разностная схема во внутренней области </span><em><span>G</span></em><span> записывается на крестообразном шаблоне с центром в узле </span><img src="http://content.snauka.ru/web/80911_files/1(5).gif" alt="" width="340" height="28" /><br />
<span>На сетке </span><img src="http://content.snauka.ru/web/80911_files/1(6).gif" alt="" width="24" height="25" /><span> рассматривается ячейка с центром в узле </span><img src="http://content.snauka.ru/web/80911_files/1(7).gif" alt="" width="58" height="28" /><span> и вершинами в полуцелых узлах, то есть образованными пересечением прямых, проходящих через середины отрезков, соединяющих узлы шаблона, параллельно направлениям </span><em><span>X</span></em><span> и </span><em><span>Y</span></em><span>. Размеры ячейки по этим направлениям:</span><br />
<img src="http://content.snauka.ru/web/80911_files/2.gif" alt="" width="284" height="25" /><br />
<img src="http://content.snauka.ru/web/80911_files/2(1).gif" alt="" width="308" height="28" /><br />
<span>Площадь ячейки </span><img src="http://content.snauka.ru/web/80911_files/2(2).gif" alt="" width="68" height="28" /><span>, разностные производные определяются</span><br />
<img src="http://content.snauka.ru/web/80911_files/2(3).gif" alt="" width="412" height="58" /><br />
<img src="http://content.snauka.ru/web/80911_files/2(4).gif" alt="" width="268" height="60" /><br />
<img src="http://content.snauka.ru/web/80911_files/3.gif" alt="" width="406" height="24" /><br />
<span>Здесь введены безиндексные обозначения для размеров ячейки и производных </span><img src="http://content.snauka.ru/web/80911_files/3(1).gif" alt="" width="208" height="30" /><span>Вводятся также безиндексные обозначения для потоков через грани ячейки:</span><br />
<img src="http://content.snauka.ru/web/80911_files/3(2).gif" alt="" width="290" height="28" /><br />
<img src="http://content.snauka.ru/web/80911_files/3(3).gif" alt="" width="297" height="30" /><br />
<img src="http://content.snauka.ru/web/80911_files/3(4).gif" alt="" width="264" height="24" /><br />
<span>В соответствии с интегро-интерполяционным методом построения разностных схем уравнение теплового баланса для ячейки имеет вид:</span><br />
<img src="http://content.snauka.ru/web/80911_files/4.gif" alt="" width="326" height="34" /><span> (2)</span><br />
<span>Здесь использовано обозначение </span><img src="http://content.snauka.ru/web/80911_files/4(1).gif" alt="" width="174" height="29" /><span>, </span><img src="http://content.snauka.ru/web/80911_files/4(2).gif" alt="" width="66" height="20" /><span>. </span><br />
<span>Узлам, лежащим на границе </span><em><span>дG</span></em><span>, будут соответствовать шаблоны и ячейки несколько иного вида. Если дополнить шаблон в этих точках фиктивными узлами, то уравнение баланса в них запишется так же, требуется только положить нулевым соответствующие фиктивные шаги. В общем случае размеры ячейки будут:</span><br />
<img src="http://content.snauka.ru/web/80911_files/4(3).gif" alt="" width="340" height="57" /><span> (3)</span><br />
<span>где </span><img src="http://content.snauka.ru/web/80911_files/5.gif" alt="" width="244" height="54" /><span>, </span><img src="http://content.snauka.ru/web/80911_files/5(1).gif" alt="" width="262" height="54" /><span>.</span><br />
<span>Потоки через пограничные грани ячейки определяются в случае граничного узла из краевых условий исходной задачи:</span><br />
<img src="http://content.snauka.ru/web/80911_files/5(2).gif" alt="" width="246" height="28" /><br />
<img src="http://content.snauka.ru/web/80911_files/5(3).gif" alt="" width="285" height="28" /><br />
<img src="http://content.snauka.ru/web/80911_files/5(4).gif" alt="" width="356" height="29" /><br />
<img src="http://content.snauka.ru/web/80911_files/6.gif" alt="" width="357" height="26" /><br />
<span>Разностная схема на сетке </span><img src="http://content.snauka.ru/web/80911_files/6(1).gif" alt="" width="30" height="25" /><span>:</span><br />
<img src="http://content.snauka.ru/web/80911_files/6(2).gif" alt="" width="260" height="68" /><span>, (4)</span><br />
<span>где</span><br />
<img src="http://content.snauka.ru/web/80911_files/6(3).gif" alt="" width="409" height="60" /><br />
<img src="http://content.snauka.ru/web/80911_files/6(4).gif" alt="" width="410" height="60" /><br />
<img src="http://content.snauka.ru/web/80911_files/7.gif" alt="" width="426" height="62" /><br />
<img src="http://content.snauka.ru/web/80911_files/7(1).gif" alt="" width="449" height="60" /><br />
<span>Таким образом, на каждом временном слое </span><img src="http://content.snauka.ru/web/80911_files/8.gif" alt="" width="140" height="25" /><span> получена нелинейная система алгебраических уравнений (4), в которую решение подобной системы на предыдущем слое входит как неизвестная функция. На нулевом слое задано начальное распределение:</span><br />
<img src="http://content.snauka.ru/web/80911_files/8(1).gif" alt="" width="112" height="32" /><span>.</span><br />
<strong><span>3. Модель организации внешнего и внутреннего итерационных процессов.</span></strong><br />
<span>Для организации внешнего интеграционного процесса вводится вектор-функция </span><img src="http://content.snauka.ru/web/80911_files/8(2).gif" alt="" width="305" height="26" /><span> и операторы </span><img src="http://content.snauka.ru/web/80911_files/8(3).gif" alt="" width="62" height="25" /><span>, определенные равенствами:</span><br />
<img src="http://content.snauka.ru/web/80911_files/8(4).gif" alt="" width="252" height="58" /><span>, тогда система (4) запишется так:</span><br />
<img src="http://content.snauka.ru/web/80911_files/8(5).gif" alt="" width="177" height="29" /><span>. (5)</span><br />
<span>Для решения ее используется линейно-квадратический процесс. Система (5) переписывается в виде: </span><img src="http://content.snauka.ru/web/80911_files/8(6).gif" alt="" width="57" height="24" /><img src="http://content.snauka.ru/web/80911_files/8(7).gif" alt="" width="222" height="38" /><br />
<span>(здесь </span><img src="http://content.snauka.ru/web/80911_files/9.gif" alt="" width="16" height="33" /><span> – значение функции на верхнем временном слое), применяется линейно-квадратичный итерационный процесс:</span><br />
<img src="http://content.snauka.ru/web/80911_files/9(1).gif" alt="" width="314" height="29" /><span>.</span><br />
<span>Верхним индексом помечается номер итерации. Здесь </span><img src="http://content.snauka.ru/web/80911_files/9(2).gif" alt="" width="25" height="29" /><span> и </span><img src="http://content.snauka.ru/web/80911_files/9(3).gif" alt="" width="25" height="29" /><span> линейные операторы</span><br />
<img src="http://content.snauka.ru/web/80911_files/9(4).gif" alt="" width="269" height="66" /><span>,</span><br />
<span>где например,</span><br />
<img src="http://content.snauka.ru/web/80911_files/9(5).gif" alt="" width="404" height="70" /><br />
<img src="http://content.snauka.ru/web/80911_files/10.gif" alt="" width="372" height="88" /><br />
<img src="http://content.snauka.ru/web/80911_files/10(1).gif" alt="" width="226" height="37" /><span>.</span><br />
<span>Внутренний итерационный процесс организуется следующим образом. На каждой итерации внешнего процесса линейная система разностных уравнений записывается без индексов внешней итерации, вводится диагональный оператор </span><em><span>D</span></em><span> так, что </span><img src="http://content.snauka.ru/web/80911_files/10(2).gif" alt="" width="94" height="24" /><span>, и в операторном виде система имеет вид:</span><br />
<img src="http://content.snauka.ru/web/80911_files/10(3).gif" alt="" width="182" height="24" /><span>.</span><br />
<span>Для приближения решения этого оператора уравнения применяется двухслойная итерационная схема с чебышевским упорядоченным набором параметров</span><br />
<img src="http://content.snauka.ru/web/80911_files/10(4).gif" alt="" width="300" height="62" /><span>.</span><br />
<span>Здесь </span><img src="http://content.snauka.ru/web/80911_files/11.gif" alt="" width="37" height="24" /><span> – упорядоченный чебышевский набор параметров. В пространственной сеточной функции </span><img src="http://content.snauka.ru/web/80911_files/11(1).gif" alt="" width="28" height="24" /><span> в смысле некоторого скалярного произведения удовлетворяются условия самосопряженности, положительной определенности и ограниченности оператора </span><br />
<img src="http://content.snauka.ru/web/80911_files/11(2).gif" alt="" width="121" height="25" /><span>,</span><br />
<img src="http://content.snauka.ru/web/80911_files/11(3).gif" alt="" width="268" height="25" /><span>,</span><br />
<img src="http://content.snauka.ru/web/80911_files/11(4).gif" alt="" width="226" height="25" /><span>,</span><br />
<img src="http://content.snauka.ru/web/80911_files/11(5).gif" alt="" width="84" height="25" /><span>,</span><br />
<span>гарантирующие сходимость внутреннего итерационного процесса.</span><br />
<span>Границы спектра </span><img src="http://content.snauka.ru/web/80911_files/11(6).gif" alt="" width="18" height="25" /><span> и </span><img src="http://content.snauka.ru/web/80911_files/11(7).gif" alt="" width="21" height="25" /><span> оператора </span><img src="http://content.snauka.ru/web/80911_files/11(8).gif" alt="" width="118" height="24" /><span> эффективно оцениваются по теории Гершгорина.</span><br />
<strong><span>4. Импульс градиента температур</span></strong><br />
<span>Конкретные исследования по описанию алгоритму проведены для следующих значений параметров линейной задачи:</span><br />
<img src="http://content.snauka.ru/web/80911_files/12.gif" alt="" width="200" height="60" /><br />
<img src="http://content.snauka.ru/web/80911_files/12(1).gif" alt="" width="194" height="56" /><br />
<img src="http://content.snauka.ru/web/80911_files/12(2).gif" alt="" width="312" height="57" /><br />
<span>где </span><img src="http://content.snauka.ru/web/80911_files/12(3).gif" alt="" width="329" height="25" /><br />
<span>Температурное поле и характер его изменения во времени имеют общий для всех вариантов расчета характерный вид. Температура металла практически постоянна по объему, и при переходе из стали в наледь поверхности теплообмена со средой круто падает, то есть в точках наледи, лежащих на его поверхности по границе с металлом, возникает значительный градиент температуры, направленный вдоль поверхности наледи. Кривая зависимости градиента температуры в угловой точке от времени имеют характерную форму, близкую к форме импульса. Близость кривой к импульсной форме определяется величинами </span><img src="http://content.snauka.ru/web/80911_files/13.gif" alt="" width="28" height="26" /><span> и </span><img src="http://content.snauka.ru/web/80911_files/13(1).gif" alt="" width="30" height="26" /><span>.</span><br />
<span>При уменьшении </span><img src="http://content.snauka.ru/web/80911_files/13(2).gif" alt="" width="42" height="26" /><span>и </span><img src="http://content.snauka.ru/web/80911_files/13(3).gif" alt="" width="30" height="26" /><span> импульс сглаживается. Такая зависимость от </span><img src="http://content.snauka.ru/web/80911_files/13(4).gif" alt="" width="28" height="26" /><span> и </span><img src="http://content.snauka.ru/web/80911_files/13(5).gif" alt="" width="30" height="26" /><span> сохраняется при всех исследованных отношениях </span><img src="http://content.snauka.ru/web/80911_files/13(6).gif" alt="" width="48" height="26" /><span>, причем от </span><img src="http://content.snauka.ru/web/80911_files/13(7).gif" alt="" width="30" height="26" /><span> острота и амплитуда импульса зависят существенно сильнее, чем от </span><img src="http://content.snauka.ru/web/80911_files/13(8).gif" alt="" width="42" height="26" /><span>, отношение </span><img src="http://content.snauka.ru/web/80911_files/13(9).gif" alt="" width="49" height="30" /><span> слабо влияют на форму зависимости и определяет преимущественно амплитуду импульса.</span><br />
<span>Градиент температуры в угловой точке достигает своего максимального значения примерно в одно и то же время (около одной секунды с начала остывания) при различных значениях </span><img src="http://content.snauka.ru/web/80911_files/13(10).gif" alt="" width="162" height="26" /><span>из исследованных интервалов. Крутизна фронта пропорциональна его амплитуде и растет с увеличением как </span><img src="http://content.snauka.ru/web/80911_files/13(11).gif" alt="" width="42" height="26" /><span> и </span><img src="http://content.snauka.ru/web/80911_files/13(12).gif" alt="" width="30" height="26" /><span>, так и с увеличением отношения </span><img src="http://content.snauka.ru/web/80911_files/13(13).gif" alt="" width="49" height="30" /><span>. В то же время крутизна спада слабо зависит от отношения </span><img src="http://content.snauka.ru/web/80911_files/13(14).gif" alt="" width="49" height="30" /><span> и определяется главным образом значениями </span><img src="http://content.snauka.ru/web/80911_files/14.gif" alt="" width="28" height="26" /><span> и </span><img src="http://content.snauka.ru/web/80911_files/14(1).gif" alt="" width="30" height="26" /><span>.</span><br />
<span>Полученные результаты позволяют объяснить предпочтения в выборе материала кровли.</span><br />
<strong><span>5. Модель процесса замерзания жидкого слоя</span></strong><br />
<span>Математическое моделирование нестационарных тепловых полей состоит в решении нестационарного уравнения теплопроводности в двумерной области </span><em><span>G</span></em><span> с соответствующими граничными и начальными условиями:</span><br />
<img src="http://content.snauka.ru/web/80911_files/14(2).gif" alt="" width="129" height="106" /><span> (6)</span><br />
<span>Особенностью рассматриваемой задачи является необходимость учета фазового перехода из жидкого в твердое состояние.</span><br />
<span>Для сквозного счета таких задач без явного выделения фронта затвердевания нужно учесть, что при температуре фазового перехода </span><img src="http://content.snauka.ru/web/80911_files/14(3).gif" alt="" width="50" height="28" /><span>энергия </span><em><span>Е</span></em><span>, как функция температуры, испытывает переход величины </span><img src="http://content.snauka.ru/web/80911_files/14(4).gif" alt="" width="24" height="33" /><span>, который называется теплотой фазового перехода, поэтому для энергии справедливо:</span><br />
<img src="http://content.snauka.ru/web/80911_files/14(5).gif" alt="" width="213" height="60" /><span>, где</span><br />
<img src="http://content.snauka.ru/web/80911_files/14(6).gif" alt="" width="180" height="52" /><br />
<span>Это выражение подставляется в уравнение энергии:</span><br />
<img src="http://content.snauka.ru/web/80911_files/14(7).gif" alt="" width="120" height="46" /><span>, и учитывая что </span><img src="http://content.snauka.ru/web/80911_files/15.gif" alt="" width="109" height="46" /><br />
<span>есть дельта-функция Дирака, получается уравнение:</span><br />
<img src="http://content.snauka.ru/web/80911_files/15(1).gif" alt="" width="282" height="46" /><span>, справедливое и в области фазового перехода. Выражение </span><img src="http://content.snauka.ru/web/80911_files/15(2).gif" alt="" width="42" height="24" /><span> и </span><img src="http://content.snauka.ru/web/80911_files/15(3).gif" alt="" width="93" height="33" /><span> входят в уравнение одинаковым образом, причем </span><img src="http://content.snauka.ru/web/80911_files/15(4).gif" alt="" width="93" height="33" /><span> представляет собой сосредоточенную теплоемкость на поверхности </span><img src="http://content.snauka.ru/web/80911_files/15(5).gif" alt="" width="50" height="28" /><span>.</span><br />
<span>Для перехода к разностной схеме заменяется дельта-функция приближенно – образной или размазанной функцией </span><img src="http://content.snauka.ru/web/80911_files/15(6).gif" alt="" width="121" height="33" /><span>, где – величина полуинтервала, на котором функция </span><img src="http://content.snauka.ru/web/80911_files/15(7).gif" alt="" width="94" height="33" /><span>отлична от нуля.</span><br />
<span>Таким образом, вводится сглаженная или эффективная теплоемкость </span><img src="http://content.snauka.ru/web/80911_files/16.gif" alt="" width="224" height="34" /><span>, которая удовлетворяет условию </span><img src="http://content.snauka.ru/web/80911_files/16(1).gif" alt="" width="100" height="26" /><span> вне интервала </span><img src="http://content.snauka.ru/web/80911_files/16(2).gif" alt="" width="138" height="33" /><span>.</span><br />
<span>Изменение энтальпии на интервале </span><img src="http://content.snauka.ru/web/80911_files/16(3).gif" alt="" width="137" height="33" /><span> сокращается, т.е.</span><br />
<img src="http://content.snauka.ru/web/80911_files/16(4).gif" alt="" width="244" height="81" /><span>.</span><br />
<span>На интервале </span><img src="http://content.snauka.ru/web/80911_files/16(5).gif" alt="" width="138" height="33" /><span>можно, например, взять </span><img src="http://content.snauka.ru/web/80911_files/16(6).gif" alt="" width="138" height="50" /><span>, что будет соответствовать интерполяции – функции с помощью прямоугольного импульса. На том же интервале производится сглаживание коэффициента теплопроводности . Вводится сглаженный, или, эффективный коэффициент </span><img src="http://content.snauka.ru/web/80911_files/17.gif" alt="" width="48" height="26" /><span>, совпадающий с </span><img src="http://content.snauka.ru/web/80911_files/17(1).gif" alt="" width="45" height="25" /><span> при </span><img src="http://content.snauka.ru/web/80911_files/17(2).gif" alt="" width="80" height="29" /><span> и с </span><img src="http://content.snauka.ru/web/80911_files/17(3).gif" alt="" width="48" height="25" /><span> при </span><em></em><img src="http://content.snauka.ru/web/80911_files/17(4).gif" alt="" width="82" height="24" /><span>.</span><br />
<span>Например, если задавалось</span><br />
<img src="http://content.snauka.ru/web/80911_files/17(5).gif" alt="" width="201" height="62" /><br />
<span>То можно взять</span><br />
<img src="http://content.snauka.ru/web/80911_files/17(6).gif" alt="" width="489" height="118" /><br />
<span>В результате получается задача для уравнения теплопроводности со сглаженными коэффициентами:</span><br />
<img src="http://content.snauka.ru/web/80911_files/18.gif" alt="" width="164" height="42" /><span>,</span><br />
<img src="http://content.snauka.ru/web/80911_files/18(1).gif" alt="" width="154" height="32" /><span>.</span><br />
<span>Моделирование процессов затвердевания границ льда и основы позволяют выбрать композитные тонкие слои по границам гетероструктур [1].</span></p>
<p><strong><span>Выводы</span></strong><br />
<span>1. Рассмотрена фундаментальная модель процесса теплообмена в гетероструктуре в форме нестационарного уравнения теплопроводности.</span><br />
<span>2. Описана конечно-разностная модель гетероструктур с введением пространственной сетки и применением интегро-интерполяционного метода построения разностных схем.</span><br />
<span>3. Предложена эффективная модель организации внешнего и внутреннего итерационных процессов с применением двухслойной итерационной схемы с чебышевским упорядоченным набором параметров.</span><br />
<span>4. Модель процесса замерзания жидкого слоя наледи построена без явного выделения фронта затвердевания с учетом теплоты фазового перехода.</span><br />
<span>5. Проведены системные исследования тепловых процессов в гетероструктурах, установлены закономерности протекания нестационарных процессов.</span><br />
<strong><span>Заключение</span></strong><br />
<span>1. Модели динамики гетероструктур снеголедовых масс с фазовыми переходами на границах позволили установить новые закономерности нестационарных процессов обледенения наиболее распространенных в технике гетероструктур.</span><br />
<span>2. Использование установленных закономерностей в практике проектирования гетероструктур представляет новые возможности в получении более безопасных конструкций.</span><br />
<span>3. Наиболее распространенные гетероструктуры крыш зданий и сооружений целесообразно устраивать с разными коэффициентами трения точечно по всей крыше и узкой (около 1 %), полосе гидрофобного композита по краям крыш для снижения размера сосулек до безопасного. В качестве тонкослойных композитов для разных материалов крыш разработаны и испытаны эффективные, весьма долговечные и недорогие составы, а также технологии их нанесения с учетом конкретных условий применения.</span></p>
]]></content:encoded>
			<wfw:commentRss>https://web.snauka.ru/issues/2017/04/80911/feed</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
	</channel>
</rss>
