ddddddddddddddddddddddddddsssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssss
КИПЕНИЕ Кипение - превращение жидкости в пар Кипение можно создать: нагревом поверхности тепловыделением в объеме Условие начала кипения T > Ts Кипение - процесс парообразования в толще жидкости Теплообмен при кипении интенсивнее, чем в однофазной жидкости Три причины повышенной интенсивности теплообмена: сильная турбулизация пограничного слоя за счет роста и отрыва паровых пузырей перенос тепла от поверхности при испарении микрослоя жидкости у основания пузыря в паровой пузырь перенос скрытой теплоты парообразования внутри пузыря Классификация процессов кипения - по роду кипения : пузырьковое, пленочное; - по типу конвекции кипение в большом объеме (кипение при свободной конвекции) кипение в каналах (при вынужденной конвекции) - по соотношению средней температуры жидкости (Tf) и температуры насыщения (Ts) жидкость недогрета до температуры насыщения (поверхностное кипение) Т f < Т s Т ≈ Тs жидкость догрета до температуры насыщения f - по характеру кипения: развитое, неразвитое, неустойчивое Кривая кипения Зависимость плотности теплового потока q от температурного напора (Tw - Ts) при кипении в большом объеме Нукияма Nukiyama (1934) 5 – свободная конвекция жидкости; В – начало кипения – неразвитое пузырьковое кипение; - развитое пузырьковое кипение; С' – DNB (departure from nucleate boiling - отклонение от пузырькового кипения); С – точка кризиса теплообмена; – переход от пузырьк к пленочн.; – пленочное кипение; D – точка Лейденфроста Явление Лейденфроста q ~ ∆T 3÷ q ~ ∆T \\\\\ - неустойчивое кипение α ~ ∆T ÷1 q = α ⋅ ∆T ~ ∆T ÷4 Зарождение парового пузыря Зарождение парового пузырька на твердой поверхности происходит в центрах парообразования Работа образования парового пузыря объемом V с полной площадью поверхности F, часть которой f опирается на твердую стенку L = −∆p V + σ П − Ж ( F − f ) + σТв − П f − σТв − Ж f F σ П−Ж σТв − П = σТв − Ж + σ П − Ж сosθ V σТв− П Условия равновесия в точке σТв− Ж L = −∆p V + σ f f F L П−Ж f F 1 − (1 − cos θ) F 10 Расчетные соотношения Развитое пузырьковое кипение воды на чистой поверхности в большом объеме Формула ЦКТИ α = 4,34q 0, (p 0,14 + 0,0135 p ) Вт/м2К q - Вт/м2; p - МПа α = 10,45q 0, [ 3,3 − 0,0113 ⋅ (Ts − 373)] −1 [Ts] = K Формула Лабунцова 3,4 ⋅ Р α= ⋅q − 0,0045Р ,18 18 Кризис теплообмена АВ – свободная конвекция; ВС – пузырьковое кипение; СD – переходное пленочное кипение; DF – устойчивое пленочное кипение; ВМ – затянутая свободная конвекция Точка С – переход пузырькового кипения в пленочное Причина ухудшения теплообмена: с наступлением пленочного кипения большая часть поверхности и в течение большего времени омывается паром Т.к интенсивность теплообмена к пару значительно меньше, чем при испарении, температура поверхности возрастает 19 КУТАТЕЛАДЗЕ САМСОН СЕМЕНОВИЧ (18.07.1914 - 20.03.1986) Гидродинамическая теория кризиса теплообмена С.С.Кутателадзе Переход пузырькового кипения в пленочное носит черты кризиса и в большинстве случаев определяется устойчивостью паровой пленки (гидродинамический процесс) Первая критическая плотность теплового потока qкр1 переход от пузырькового кипения к пленочному Вторая критическая плотность теплового потока обратный переход qкр2 21 Гидродинамическая теория кризиса теплообмена С.С.Кутателадзе qкр1, соответствующая переходу пузырькового кипения в пленочное, определяется в основном физическими свойствами жидкости и пара Возникновение пленочного режима кипения, является результатом такого нарушения устойчивости структуры пристенного двухфазного слоя, при котором жидкость оказывается оттесненной образующимся паром от нагреваемой поверхности Такой же эффект наблюдается при барботаже газа через жидкость 22 Гидродинамическая теория кризиса теплообмена С.С.Кутателадзе Стабилизирующий фактор - силы тяжести и поверхностного натяжения, возмущающий фактор - динамический напор образующегося пара - ρ"w"2 Динамический напор равен работе оттеснения жидкости от поверхности ρ "W "2 ≈ g ⋅ ( ρ '− ρ ") ⋅ δ δ - средняя толщина паровой пленки, которая определяется капиллярной постоянной l = σ / g ( ρ'−ρ") ρ "W "k р = g ( ρ '− ρ ")δ ρ "W "k р gσ ( ρ '− ρ ") = K ≅ const 23 Гидродинамическая теория кризиса теплообмена С.С.Кутателадзе Скорость паровой фазы для небольших давлений имеет порядок Тогда W "к р = qк р r ⋅ ρ" qк р ≅ K ⋅ r ⋅ ρ" ⋅ g σ (ρ'−ρ" ) Критерий устойчивости K ≈ 0,13 ÷ 0,16 r - удельная теплота парообразования, КДж/кг ρ' , ρ" - плотности жидкости и пара на линии насыщения, кг/м σ - коэффициент поверхностного натяжения, Н/м 24 Кризис теплообмена 25 Термодинамическая теория кризиса теплообмена При 0,1 МПа предельный перегрев составляет ~200 K, при 14 МПа ~15 K Т2-Т1 максимальный перегрев Наступление кризиса связывается с условиями, когда температура поверхности достигает предельной температуры состояния жидкости (температуры Лейденфроста) 26 Теплообмен при пленочном кипении – теплопроводность; – конвекция; – излучение; δ - средняя толщина паровой пленки; lкр – критическая длина волны (колебания границы раздела фаз) имеет порядок капиллярной постоянной 27 Теплообмен при пленочном кипении Допущение: пар движется под действием только подъемной силы (справедливо при δ [...]... большом объеме Давление Р = 5, 5 МПа, тепловая нагрузка q = 0,2 МВт/м2 Найти: Коэффициент теплоотдачи α и температуру поверхности нагрева tw При Р =5, 5 МПа ts = 270 оС Коэффициент теплоотдачи Температура стенки Формула ЦКТИ α = 4,34 ⋅ q 0, 7 ( P 0,14 + 0,0135P 2 ) = = 3,74.104 Вт/м2К tw = ts + q = α 276,8 оС Формула Лабунцова 2 3,4 ⋅ Р 0,18 α= ⋅ q 3 = 2,93.104 Вт/м2К 1 − 0,00 45 tw = 2 75, 3 оС 32 ... + 0,5c p ,п (t w − t s ) эффективная теплота парообразования, учитывающая перегрев пара в пленке Физические свойства перегретого пара по t = 0 ,5( t w + t s ) Пленочное кипение на вертикальной поверхности (турбулентное течение пленки): α = 0, 25 3 (λп ) 2 ⋅ с р ,п ⋅ g ⋅ ( ρ '− ρ п ) νп 30 Теплообмен при пленочном кипении Суммарный коэффициент теплообмена: За счет излучения: α рад εw (T ) − Ts4 = 5, 67... Расчетные соотношения Развитое пузырьковое кипение воды на чистой поверхности в большом объеме Формула ЦКТИ α = 4,34q 0, 7 (p 0,14 + 0,01 35 p 2 ) Вт/м2К q - Вт/м2; p - МПа α = 10,45q 0, 7 [ 3,3 − 0, 0113 ⋅ (Ts − 373)] −1 [Ts] = K Формула Лабунцова 2 3,4 ⋅ Р α= ⋅q 3 1 − 0,00 45 0 ,18 18 Кризис теплообмена АВ – свободная конвекция; ВС – пузырьковое кипение; СD – переходное пленочное кипение; DF – устойчивое... парообразования, КДж/кг ρ' , ρ" - плотности жидкости и пара на линии насыщения, кг/м σ - коэффициент поверхностного натяжения, Н/м 3 24 Кризис теплообмена 25 Термодинамическая теория кризиса теплообмена При 0,1 МПа предельный перегрев составляет ~200 K, при 14 МПа ~ 15 K Т2-Т1 максимальный перегрев Наступление кризиса связывается с условиями, когда температура поверхности достигает предельной температуры состояния... молекулярным сцеплением между жидкостью и твердым телом) При хорошей смачиваемости (краевой угол мал (θ→0)) углубления заполняются жидкостью, что затрудняет вскипание и приводит к большим перегревам поверхности 11 Рост парового пузыря Чтобы сферический пузырек существовал в жидкости, а не схлопывался под действием сил поверхностного натяжения, Рп > Рокр, R > Rкр Критический (минимальный) радиус парового зародыша... rρ" (Tw − Ts ) скорость парообразования q W "= rρ" 14 Влияние давления При π = p > 0,2 влияние давления резко увеличивается pкр Причина - характер изменения физических свойств жидкости ( r, σ, ρ', ρ" 15 ) Влияние шероховатости Увеличение шероховатости все каверны одного размера нормальное распределение каверн по размерам специально подобранное распределение каверн по размерам 16 Влияние отложений Отложения ... ЦКТИ α = 4,34q 0, (p 0,14 + 0,01 35 p ) Вт/м2К q - Вт/м2; p - МПа α = 10,45q 0, [ 3,3 − 0, 0113 ⋅ (Ts − 373)] −1 [Ts] = K Формула Лабунцова 3,4 ⋅ Р α= ⋅q − 0,00 45 ,18 18 Кризис теплообмена АВ –... поверхностного натяжения, Н/м 24 Кризис теплообмена 25 Термодинамическая теория кризиса теплообмена При 0,1 МПа предельный перегрев составляет ~200 K, при 14 МПа ~ 15 K Т2-Т1 максимальный перегрев Наступление... ⋅ ν п ⋅ ∆t ⋅ d r* = r + 0,5c p ,п (t w − t s ) эффективная теплота парообразования, учитывающая перегрев пара в пленке Физические свойства перегретого пара по t = 0 ,5( t w + t s ) Пленочное кипение