Разработка метода расчета пластинчатого скребкового теплообменного аппарата Текст научной статьи по специальности «Механика и машиностроение»

ВестникВГУИТ, № 3, 2012

УДК66.015.24(045)

Аспирант К.А. Рашкин, профессор С.А. Бредихин,

(Московский гос. ун-т пищ. произ-в) кафедра процессов и аппаратов, тел. (926) 152-70-88

доцент В.М. Чесноков

(Московский гос. ун-т пищ. произ-в) кафедра высшей математики, тел. (495) 676-48-13

Разработка метода расчета пластинчатого скребкового теплообменного аппарата

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

The current work presents a method of calculating the required heat exchange area for the scraper plate- type heat exchanger for high viscosity products. Calculations are per-formedusing theappliedcomputingpackageMathCAD. Based on results of calculations, tem-peraturegraphs were constructedforstraight-flowandcountercurrentflowpatterns in liquid environment.

Ключевые слова: теплообмен, скребковые теплообменники, методика расчета.

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

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

Расчет скребкового пластинчатого тепло-обменного аппарата проводим по следующей методике:

- задаем производительность аппарата;

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

© Рашкин К.А., Бредихин С.А., Чесноков В.М., 2012

Рис. 1. Схема движения продукта в аппарате

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

Рис. 2. Расчетная схема: 1 - теплообменные пластины; 2 - продукт; Б0 - размер, ограничивающий продуктовую зону, Б0=2Я0; Б1 - диаметр входного отверстия, Б1=2Я1; Б2 - диаметр, на котором расположены выходные отверстия, Б2=2Я2; Р1 - давление продукта на входе в теплообменный элемент; Р 2 - давление продукта на выходе из теплообменного элемента, То - температура продукта; Т1 - температура пластины 1; Т 2 - температура пластины 2

Для расчета температуры продукта на выходе из каждого теплообменного элемента используем уравнения(1), (2) и граничные условия (3):

дТ дг

+ V,

дТ дг

^д 2Т

= а

1 дТ +--+

д 2Т^

дг г дг дг2

(1)

где Т - температура в точках продукта, °С ; г , г - цилиндрические координаты точки продукта; vr, vz - проекции скорости точек продукта на оси г ; г, а - коэффициент температуропроводности;

г

= V

1 др

р дг

дг2

д2 V,,

(2)

1 ду +--- + - 2

г дг дг

где р - давление в точках продукта; V - коэффициент кинематической вязкости продукта;

г = Ях, г = Я2 ,Р = р ,

р = Р2,Т = Т, Т = Т, , (3)

Т = Т3, г = 0= к, vr(г,0) = 0, /(0) = 0 где Р\ - модифицированное давление на входе продукта в междисковое пространство; Р2 -модифицированное давление на выходе

продукта; к - расстояние между дисками; Т3 и Т4 - температура продукта на стенках дисков.

Решение уравнений (1) и (2) при граничных условиях (3) имеет вид

2А к2

Т (г,г )=_ ап _ Т3 _(Т' _ Т4) - * 2 _ * _ г 4 _ *> _ 2 - п

_г ) . кя I

81П-г)

к I

+Т3 _ -(Тз _ Т4) к

(4)

Однако каждый теплообменный элемент должен рассчитываться, исходя из своих граничных условий (температура на входе и температуры охлаждающих стенок). По этой причине непосредственное применение выражения (4) для практических расчетов малопригодно, так как нужно знать температуру продукта на входе в каждый элемент и температуру стенок дисков в этом элементе. В силу этого расчет теплообменного элемента будет состоять из двух основных этапов: 1 этап -определение температуры стенок каждого элемента охладителя и 2 этап - определение температуры продукта на выходе из охладителя.

Температура стенок охладительного элемента, или температура хладагента, не может быть рассчитана аналитически ввиду сложности канала, по которому течет хладагент. Поэтому найдем ее приближенно, считая, что она одинакова на всех стенках данного элемента. Для этого используем формулу, аналогичную формуле В.Г. Шухова, приме ни-

V

г

ВесжнмкВГУИТ, № 3, 202

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

Т = Ти + Г21

2

где 2£К

+

ТЖ, 1_1

Т + Т Л 1 11 ^ 1 21

2

е , (5)

(6)

Ср РЧх

I - средняя длина канала для хладагента в каждом элементе; £ - средняя площадь сечения канала для хладагента в каждом элементе; К - коэффициент теплопередачи через стенку между хладагентом и продуктом; р- плотность хладагента; с- теплоемкость хладагента при постоянном давлении; т и Т21 - соответственно температура продукта на входе в элемент с номером 1 и температура продукта на выходе из этого элемента; Тх г-1 - температура хладагента на входе в 1 -й элемент; Т + Т

_21 - средняя температура продукта в

2

каждом элементе; дх - расход хладагента.

Коэффициент р не может быть рассчитан по формуле (6), поскольку для каналов сложной формы коэффициент теплопередачи К может быть определен только экспериментально. Поэтому определим весь комплекс

_в, входящий множителем в показатель

к

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

Тх1 = Тх0 ±А, где Тх1 - средняя температура хладагента в первом элементе; Тх0 - температура хладагента на входе в первый элемент (при прямотоке) или на выходе из первого элемента (при противотоке), А - разность температур хладагента в первом элементе и на входе (выходе) из него, которую можно найти экспериментально. Знак плюс соответствует случаю прямотока, а знак минус - случаю противотока. Тогда из (5) получим уравнение для определе-

ния комплекса

Ё! к

Тхо ±А = Т11 + Т21. + (Тх0 _ Т11 + Т21 1 ев

(7)

2 ^ хи 2 в котором температура продукта Т21 на выходе из первого элемента рассчитывается по

формуле

(4)

при

г = Я,

Т = Т4 = Тх1 = Тх0 ±А , Т1 = Тп и осреднении функций, зависящих от z по толщине зазора между дисками.

Из (7) находим

_В- = 1п

к

1 ±-

2 А

2Т _Т _Т

х0 21 У

(8)

На основании (8) формула (5) примет

вид

Г = Т1, + Т2г +ГГ _ Т11 + Т21

Тх1 ~ „ "Ч Тх,1 _1 „

1 ±-

2 А

2Т _Т _Т

х0 -41 21 У

(9)

Таким образом, по формуле (9) можно приближенно найти температуру стенок для элементов 1 = 1,2,... N.

Преобразуем теперь формулу (4) так, чтобы с ее помощью можно было бы найти температуру продукта на выходе из любого охлаждающего элемента. Для этого положим, что температура всех стенок данного элемента с номером 1 примерно одинакова и определяется по формуле (9). Тогда в (4) следует положить Т3 = Т4 = ТХ1, Т1 = Ть-, а истинную температуру в каждой точке охлаждающего элемента заменить средней температурой по толщине зазора между дисками при радиусе г = Я2, т. е. Т(г, z) « Т21, где Т и Т21 - температура продукта на входе в элемент с номером 1 и температура продукта на выходе из этого элемента. После этого преобразованная формула примет следующий вид:

Т2, =_П[Т _Тх, )(1 _со*кп)

Т<Я1_Я2 >1 _ С08 кп

Я22 _ Я,2

^Я24 _ Я-4) _

(10)

По формуле (10) последовательно под-считывается температура на выходе из любого охлаждающего элемента, начиная с первого по ходу течения продукта. Подсчет прекращается для такого количества теплообменных элементов 1 = N, при котором температура продукта на выходе из этого элемента станет равной заданной конечной температуре Т2N. Следует иметь в виду, что температура про-

ВестникВГУИТ, № 3, 2012

дукта на выходе из любого элемента Т2г присутствует как в равенстве (9), так и в равенстве (10). Поэтому ее определение будет связано с решением системы уравнений (9) и (10).

В равенство (10) вместо Т^ подставляется правая часть (9) и решается уравнение относительно Т2г, при этом температура Т на входе

в элемент с номером г равна температуре продукта на выходе из предыдущего г _ 1 элемента.

Пример расчета. Определить количество теплообменных элементов для охладителя сгущенного молока для следующих исходных данных: температура продукта на входе -Т = 80° С, температура продукта на выходе -Т2 N = 16° С, температура хладагента на входе в первый элемент (при прямотоке) или на выходе из первого элемента (при противотоке) -Тх0 = 2° С, радиус входного отверстия -= 0,041 , радиус, на котором расположены выходные отверстия - Я2 = 0,21 , расстояние между пластинами - k = 0,0151 , коэффициент температуропроводности продукта -a = 10_71 2/с, расход продукта -

q = 4 -10_413/с, разность температур хладагента на входе и выходе в первом элементе -

д = 0,8°й.

Пример расчета приведен для движения продукта и хладоносителя по схеме прямотока и противотока.Все численные расчеты количества теплообменных элементов проведены в прикладном вычислительном пакете МЛТНСЛБ.

Прямоток. Из расчета следует, что заданная температура продукта Т2Н = 16° С примерно достигается при количестве теплооб-менных элементов i = N = 20 шт и составляет 15, 9 °С. Зависимость температуры продукта на выходе из последовательно установленных теплообменных элементов при прямоточной схеме течения сред показана на рис. 3. При этом температура хладоносителя на выходе из аппарата равна Тх20 = 8,4 °С.

Для противоточной схемы течения сред (рис. 4) задаем температуру хладоносителя Тх0 = 12 °С и разность температур Д = 1 °С

Зависимость температуры продукта на выходе из последовательно установленных теплообменных элементов при противоточной схеме течения сред показана на рис. 4. При этом температура хладоносителя на выходе из аппарата равнаТх20 = 2,3 °С.

80 ¡50 40

О

о

п

м е

Количество пластин, шт

Рис. 3. Зависимость температуры продукта от количества теплообменных элементов при прямоточной схеме течения сред: Т& - температура хладагента в г-м элементе, - выходная температура продукта из i - го элемента; г - количество охладительных элементов

30

60

о ей

а

Ь 20 S а

<D

е

S n

(D 0

[— 0 5 10 15 20

Рис. 4. Температурный график скребкового тепло-обменного аппарата при противоточной схеме течения сред

Из расчета следует, что заданная температура продукта T2N = 16 °C примерно достигается при количестве теплообменных элементов i = N = 20 шт и составляет 15, 9 °С.

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

ЛИТЕРАТУРА

1. Бредихин, С.А. Технологическое оборудование предприятий молочной промышленности [Текст] / С.А. Бредихин. - М.: Колос, 2010. - 408 с.

2. Коршак, А.А. Трубопроводный транспорт нефти, нефтепродуктов и газа [Текст]: учебное пособие / А.А. Коршак, А.М. Нечваль. - Уфа: «ДизайнПолиграфСервис», 2005. - 516 с.

Количество пластин, шт