книги / Технология композиционных материалов
..pdfДля перехода от этих зависи мостей к зависимостям от темпера туры необходимо учитывать, что толщина осадка 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, |
см |
S£ , см |
, г/см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