Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Технология композиционных материалов

..pdf
Скачиваний:
2
Добавлен:
12.11.2023
Размер:
6.98 Mб
Скачать
Рис.2.28. Сечение ячейки кар­ каса горизонтальной плос­ костью ( <? - длина ячейки)

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

d = d (T,t)

При описании процесса будем использовать следующие предполо­ жения:

1) осаждение пироуглерода про­ исходит послойно с повторением по­ верхности тела;

2 ) филаменты и нити двойного сложения имеют цилиндрическую форму.

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

первом л втором этапах происходит заполнение межфиламентного прост

ранства

в нитях, образующих каркас. Третий и четвертый

этапы соот­

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

каркаса и

между ними. Перейдем к

описанию каждого из

этапов.

 

Ча перьсм отапе свободное распределение филаментов внутри ни­

ти позволяет

написать

простые формулы для удельной поверхности S

и плотности

/

 

 

 

 

 

 

 

S(d) =

1йс( r0 +d)n L

 

(2.44)

 

 

 

 

 

 

 

 

 

# ( d ) = Г0 + Т с ! 5 ( z ) d x = r 0 * Гс S ( d )

d =

 

 

О

 

 

 

 

 

= r

+ ,

r.Jr.nL

,

d

. d z

 

(2.45)

Oo

oc

V

{

0 °

 

 

 

Здесь использованы следующие

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

 

 

V -

объем каркаса;

 

 

 

 

 

ft -

начальная (кажущаяся)

плотность изделия;

 

f

- плотность пироуглерода;

 

 

 

 

 

L

- общая дайна одинарной нити, образующей каркас;

 

 

rQ - радиус поперечного

сечения филамента.

 

 

 

 

Формулы

(2.44),

(2.45)

справедливы до начала

смыкания боко­

вых поверхностей филаментов,

образующих нить. Предельное

значе­

ние осадка

dj , при котором еще действуют формулы

(2.44),

(2.45),

определим позднее.

 

 

 

 

 

 

 

 

 

 

Ко второму этапу осаждения отнесем этап дальнейшего запол­

нения пироуглеродом межфиламентного пространства. Ввиду

неравно­

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

внутри нити идет

уменьшение

удельной поверхности

от

значения

S(c/f)

до значения

S(d^) ,

когда

произойдет полное

заполнение межфиламентного

пространства,

^удем считать,

что этот

процесс протекает

по линейному

закону.

Запишем условие полного

заполнения

 

 

 

 

 

 

 

•* * * =

п. л-

г* + Pf

 

 

 

(2.46)

где

R -

радиус поперечного

сечения двойной нити;

 

 

 

 

Dj - суммарная площадь в поперечном сечении нити,

занятая

 

 

пирсуглеродом.

 

 

 

 

 

 

 

Формулу (2.46)

удобно представить

в виде

 

 

 

 

 

 

 

= nsr-rf =

п S' (rQ + d2 )2'

 

(2 *47)

Здесь

под

 

понимается предельное (фиктивное) значение

радиу­

са фшгамонта, покрытого осевшим пироуглеродом. Из формулы (2-47) находим

 

 

 

d2 = rf~ Го

(2.48)

Значение

^

позволяет записать

выражение для плотности

изД'^шя

в конце

второго этапа:

 

 

^ = L n s - r f f a + L n j r ( л / - г д ) Г г

 

 

_

L - n - s r [ r t 2Гг.~П>г ( Г г - Г ф ) 1

(2.49)

 

 

V

>

 

где fcp ” плотность филаментов.

Вернемся к

нахождению величины df

Она определяется из ус­

ловия равенства

значения

в конце

второго этапа, вычислен­

ного по формуле

(2.49):

 

 

f ( d t ) =

Го

*

 

+ < £ / s ( x ) d x

Гг

(2.50)

 

$(х)

 

О

 

di

 

 

 

Здесь

в

первом интервале

определяется по формуле (2.44), а

во втором интервале SCX)

- уравнение

прямой,

проходящей

через

точки

( ^ ,

Sj (df) ) и (

d2 ,

S2 (d2)

). Из

(2.50) находим вы­

ражение для

d]

 

 

 

 

 

и -

2 (Г} -

,

c/9 (S(0)+

S(dt))

______ ________________________________

,

, ,

ЬЗГп-Ь

.

 

5(0)

S(dt )

+ ---р-----d l

В этом случав S(0)

находится по формуле (2.44), a S(dz) лег“

ко определяется по формуле (2.51) для третьего этапа в стыковоч­

ной точке d2 .

Запишем расчетные формулы для второго этапа:

dt>

* fc l ‘sfz >dx ft *&

 

1

 

 

(dt * d * d f )

 

 

Сделаем одно

замечание:

^

* подсчитанное по формуле

(2.49), да­

ет завышенное

значение,

так

как при

нахождении /* но

формуле

(2.48) учитывалось полное заполнение нити пироуглеродом, а

в дей­

ствительности и з-за неравномерного

расположения филаментов

в ни­

ти всегда имеются пустоты,

недоступные для метана.

 

 

На третьем этапе происходит осаждение пироуглерода на

внеш­

нюю поверхность нитей,

образующих каркас. В этом случав

 

 

 

d)-T - }

d2 * d * d 3

 

(2.51)

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

Эта формула справедлива до смыкания горизонтальных плос­

костей, образующих ячейку,

_ .

когда К + и = -у- . Отсюда и

находится предельное значеьие 04 и удельная поверхность S(c/^) для конца четвертого этапа:

Рис.2.30. Вид половины ячейки

к началу четвертого этапа

S'dM)=

5срмула (2.53) отражает квадратичную зависимость удельно;! поверх­

ности

S

от

толщины осадка d

. Однако для

простоты

расчетов

ее

можне

с достаточной

точностью

аппроксимировать прямой линией,

про­

ходящей

через

точки

{ d$ 9 S (d^) ) и (

,

S(c/^)

).

 

 

Таким образом,

окончательные формулы для

четвертого этапа

име­

ют виц

 

 

 

 

 

 

 

 

S(d) =

 

d P Sfdj-

 

 

 

 

d

+ & f sW d z = f s + y c SCd) + S(dj)

(2.55)

(d3 -a d a dn)

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

SCcfq +о) =о

я пиролиз прекращается. Окончательная плотность изделия определя­ ется по формуле (2.55) для d = d^

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

Г-0 = 3,5* I0-4

см;

L

= 6*10® см;

R = 5.65* КГ2

см;

п = 5000;

У = 6Д*104 см3;

 

= tt = 0,3 см;

<2 = 1,6»10“*

см;

 

= 2,2 г/см3;

fo = 0.393 г/cu 3;

у

= 1,74 г/см3.

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

 

 

 

 

Таблица 2.1

 

Значения удельной поверхности S

и плотности У

 

 

в граничных точках

 

С

dL • Ю"4,

см

, см

, г/см3

0

0

 

1081

0,393

I

0,033

 

1091,7

0,401

2

3

 

17,53

0,763

3

185,8

 

23,18

1,581

4

235,8

 

9,71

1,762

2.7.4. Исследование процесса пиролиза метана, осуществляемого термоградиентнш методом на одномерной модели

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

задачи можно ограничиться

системой координат л , £

, где

л -

радиальная координата, a

Z - вертикальная. Анализ

эксперимен­

тальных данных показывает,

что имеет место слабая

зависимость

характеристик процесса от

вертикальной координаты.

Зто дает

воз­

можность рассматривать процесс в рамках одномерной модели вдоль

координаты г*,

так как именно вдоль этой координаты

быстро изме­

няются параметры процесса - температура и плотность,

-

на чем и

основана идея термоградиентного метода.

г

 

LRQ)RJ ] -

Исследуемую область разбиваем на

три части:

 

нагреватель,

Л

 

~ заготовка

изделия,

 

 

J -

теплоизоляция.

 

 

 

 

I Rf? /?2 ]

Схема управления процессом пиролиза в области

следующая: в точке нахождения рабочей термопары

г= |

поддержи­

вается заданная температура

Т *, что обеспечивается

подачей соот­

ветствующего количества электроэнергии на нагреватель,

который

выделяет даоулевое тепло.

 

 

 

 

 

На основе обработки справочных данных построена аппроксими­

рующая формула,

отражающая зависимость

скорости осаждения пиро­

углерода от

температуры:

 

 

 

 

 

К(Т)

= 1151 ехр

^ То

г/см 1-атм

 

Тогда кинетическое уравнение, отражающее изменение плотности ма­ териала изделия, можно записать так (предполагается, что давление метана в реакторе равно атомосферноцу):

dIf*(ft)

S^cr) K(T(r,t>r)),i>0, reZef,Rgl , (2.56)

полагая ln = п

At

где

п - номер шага по времени,

 

tn -

вре­

менной интервал.

 

 

 

 

 

 

tn

 

 

 

 

 

 

 

Для фиксированного момента времени

будем считать

функ­

цию jA(r')tn)

.известной

(например,

при

£ ° = 0

имеем

#*(г,0) =

= fo). Вычисление остальных параметров процесса

распадается

на

следующие этапы:

 

 

 

 

 

 

 

<э(г)#')

 

 

 

 

 

1) по формулам (2.59) находим значения

и

 

 

 

для

г G. [ Rf9 R2] ;

 

 

 

 

 

 

 

 

 

 

 

 

 

2 ) используя функцию

 

 

 

определяем полное

сопротив­

ление

заготовки

R(in)

 

3(in)>

 

 

 

E ( t n) =

 

 

 

3) зная

величину

тока

 

находим

 

 

= j ( t n) R ( i a) / H

;

 

 

 

 

 

 

 

 

 

 

 

 

4) решая краевую

задачу

(2.57),

(2.58),

находим распределение

температуры ТС r ?t n )

;

 

 

 

 

 

 

 

 

 

 

 

 

5) из уравнения (2.56)

вычисляем распределение плотности в об­

ласти

[ R j ,R2

J

для следующего момента времени -

t п*1

 

 

 

 

Опишем подробнее

каждый

из этапов. При численной

реализации

производится дискретизация исследуемой области по

переменной

г

,

для чего строится разностная

сетка

R

= rQ <

<

< гн = ^

£у-

дем использовать

обозначения

J>п =

 

 

} in)

 

 

 

 

 

 

В разностной форме на

этапе I

имеем

 

 

 

 

 

 

 

 

<3. =

<з0 + (6М

 

й П - Гр

 

 

 

 

 

 

 

 

6°} Гт

- Го

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

*1 - ж

Го

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

±

Г 2 , ) 6 *

R

J - 1 s

Для решения краевой задачи (2.57), (2.58) воспользуемся консер­ вативной разностной схемой, которая имеет вид

~ ( г Х Т г )"

J V

T ’n ‘ 0.

* E - e f ’ O, J - 1,4-1

(2.60)

'

 

n ТЛ ___

(2.61)

~ - ^ ( T Z - 8 > ,

где конечно-разностные операторы, аппроксимирующие первую произ­ водную (разности "вперед” и"назад" соответственно),

 

г

Л * ' - . 5

г. =•

Ъ

- Т;

 

Г.

'J-±-

 

V

 

 

у

 

- л._,

а конечно-разностная аппроксимация второй производной

( г з е т Х

= —

г -

(п

Т а -

л

 

Т*)

X j

£ / - £ ■ /

1

'/г у "

 

 

'*

Систему разностных уравнений (2.60), (2.61) относительно неизвест­ ных Т.п , у = 0,/V можно решить методом скалярной прогонки.

Реализацию этапа 5 проведем следующим образом. Б!удем считать, что за единицу времени At характеристики среды существенно не меняются. Тогда приращение толщины осадка пироуглерода ( AOSn

= OSj 1 -

OSj

) в точке

у

за

вРемя tn ,

прошедшее

с момен­

та

tn, может быть вычислено с использованием формулы

 

 

 

 

aosn=

osn* ~

os?

=

. ? ч At

,

 

 

 

 

кСт?уг

 

 

где

-

плотность пироуглерода.

 

 

 

 

 

Соответственно приращение плотности

 

 

 

 

 

 

 

- $ 0 « к 4Г ' ) л о а ;

Гс

 

 

 

В тех

случаях,

когда известна

зависимость S

(у) ,

дщ

р0а.

лизации этапа 5 можно воспользоваться методом Эйлера, применив

ко­

торый к уравнению (2.56), получим

 

 

 

 

 

Г?' = Г/ +Ai syd(rp«(tp