|
|
Закрытое акционерное общество Научно-производственная фирма Центральное конструкторское бюро арматуростроения ЗАО "НПФ ЦКБА" РУКОВОДЯЩИЙ ДОКУМЕНТ МЕТОДИЧЕСКИЕ УКАЗАНИЯ ПО РАСЧЕТУ КОЛИЧЕСТВЕННЫХ ХАРАКТЕРИСТИК ТЕПЛООБМЕНА ДЛЯ ЗАДВИЖЕК РД 26-07-38-2001
Дата введения «01 » октября 2001 г. Содержание1. ОБЛАСТЬ ПРИМЕНЕНИЯНастоящий руководящий документ распространяется на задвижки, предназначенные для установки в системах атомных электростанций, опытных и исследовательских ядерных реакторах, в системах и установках химических, нефтехимических и других производств с максимальным рабочим давлением до 20 МПа (200 кгс/см2) и с температурой рабочей среды до 873К (600°С) и устанавливает порядок определения количественных характеристик теплообмена для задвижек DN от 100 до 800 мм включительно, для сред: пар, вода, воздух, азот, гелий, нефтепродукты. УСЛОВНЫЕ ОБОЗНАЧЕНИЯb - характерный геометрический размер (высота или диаметр), м; τ - время, с; t - температура, °С; Т - температура, К; tc - температура поверхности, ºС; tж - температура жидкости, газа, °С; jJ - избыточная температура, ºС; V - объём, м3, или обьёмный расход, м3/с; w - скорость, м/с; g - ускорение силы тяжести, м/с2; β - температурный коэффициент объёмного расширения, 1/°C, 1/К; ν - коэффициент кинематической вязкости, м2/с; λ - коэффициент теплопроводности, Вт/м·ºС; а - коэффициент температуропроводности, м2/с; С- коэффициент излучения, Вт/м2·К4; ε - степень черноты; bП - скорость изменения температуры рабочей среды в проточной части задвижки, К/с; bГ - скорость изменения температуры рабочей среды в зоне горловины задвижки, К/с; αП - коэффициент теплоотдачи в проточной части задвижки, Вт/м2·К; α1- коэффициент теплоотдачи в зоне горловины задвижки, Вт/м2·К; α2 - коэффициент теплоотдачи под крышкой в зоне горловины задвижки, Вт/м2·К; αi - коэффициент теплоотдачи на наружной поверхности, Вт/м2·К; αЛ- коэффициент теплоотдачи излучением, Вт/м2·К; αк- коэффициент теплоотдачи конвекцией, Вт/м2·К; αТ - коэффициент теплоотдачи с торцевой поверхности арматуры,Вт/м2·К; КРИТЕРИИ ПОДОБИЯRe = w·β/ν - критерий Рейнольдса; Рr = μ·ср /λ= ν/a - критерии Прандтля; Nu = α·β/λ - критерий Нуссельта; Gr=g·β·Δt·b3/J - критерий Грасгофа 2. НОРМАТИВНЫЕ ССЫЛКИВ настоящем руководящем документе использованы ссылки на следующие руководящие документы и стандарты: РД 302-07-122-89 «Методика расчёта температурных полей арматуры для сред с температурой 873 К (600°С)»; СТП 07.81-634-96 «Методические указания по проведению теплового расчета задвижек для АЭС»; РД 26-07-25-97 «Методика расчёта температурных полей трубопроводной арматуры»; РД 26-07-31-99 «Методические указания по расчету количественных характеристик теплообмена в арматуре клапанного типа»; Общие технические требования «Арматура для оборудования и трубопроводов АЭС» (ОТТ-87). 3.ОБЩИЕ ПОЛОЖЕНИЯ3.1. Расчёт количественных характеристик теплообмена (коэффициент теплоотдачи на внутренних и наружных поверхностях арматуры и скорости изменения температуры рабочей среды) производится с целью определения исходных данных для теплового расчета, выполняемого по программе на ПЭВМ в соответствии с РД 26-07-25-97. 3.2. Рассчитываются следующие количественные характеристики: а) коэффициент теплоотдачи в проточной части арматуры (на входе) для заданной скорости движения рабочей среды по формулам теплообмена для определённых режимов течения и типа рабочей среды; б) коэффициенты теплоотдачи (локальные) по внутренней поверхности зоны горловины; в) коэффициенты теплоотдачи с наружных поверхностей при различных условиях эксплуатации, в том числе нормальных условиях эксплуатации (НУЭ), нарушении нормальных условий эксплуатации (ННУЭ). г) расчёт скорости изменения температуры рабочей среды в зоне горловины по заданной скорости изменения температуры рабочей среды в проточной части арматуры. 3.3. Расчёт количественных характеристик теплообмена выполняется при наличии всех исходных данных, приведённых в разделе 4. Расчету предшествует составление расчётной тепловой модели. 4. ИСХОДНЫЕ ДАННЫЕ4.1 Для проведения расчета количественных характеристик теплообмена необходимы следующие данные: а) температурное поле арматуры на начальный момент времени τ = 0, К; б) закон изменения температуры рабочей и окружающей сред; в) Тн-температура рабочей среды на начальный момент времени τ = 0, К; г) Тк -температура рабочей среды на конец нестационарного режима, К; д) b -скорость изменения температуры рабочей среды в проточной части К/с; е) ta -температура окружающей среды, К; ж) геометрические размеры расчётного изделия, м; и) марки материалов деталей; к) расход рабочей среды м3/час; л) l- коэффициент теплопроводности материалов деталей, вт/м К; м) ν- кинематическая вязкость сред (рабочей и окружающей), м2/с. 5. РЕКОМЕНДАЦИИ ПО ПРОВЕДЕНИЮ РАСЧЕТА КОЛИЧЕСТВЕННЫХ ХАРАКТЕРИСТИК ТЕПЛООБМЕНА ДЛЯ ЗАДВИЖЕК5.1. Выбор начальных условий расчёта.5.1.1. При проведении теплового расчёта за начальные условия следует принять температурное поле на начальный момент времени τ = 0. Если на момент времени τ = 0 температура рабочей среды равна температуре окружающей среды, то температурное поле задвижки на начальный момент времени принимается одинаковым во всех точках и равным температуре окружающей среды. Если на момент времени τ = 0 температура рабочей среды не равна температуре окружающей среды, то начальное температурное поле задвижки принимается по расчёту соответствующего стационарного распределения температур. 5.2. Расчёт коэффициента теплоотдачи от рабочей среды к внутренней поверхности проточной частиРасчёт коэффициента теплоотдачи от рабочей среды к внутренней поверхности проточной части задвижки (на входе) производится по формуле:
где l - коэффициент теплопроводности рабочей среды, Вт/м·К; d - диаметр входного патрубка, м; Nu - критерий Нуссельта, определяемый в зависимости от характера рабочей среды и режима её течения по соответствующим формулам. 5.2.1. Формулы для жидких сред при различных режимах течения в трубах [1, 2]: а) ламинарный режим течения, когда 1/d >10; Reж>10
где Nuж = aПd/lж ; Reж = w·d/ν; Prж=νжаж ; Prс = νс /aс; ν - коэффициент кинематической вязкости рабочей среды, м2/с; Рrж, Рrc - критерий Прандтля для жидкости и стенки соответственно; 1 - характерный геометрический размер - длина участка, м; w- скорость движения рабочей среды (жидкости), м/с; а - коэффициент температуропроводности рабочей среды (жидкости) , м2/с; a, ν, Pr, - выбираются по справочной литературе. Индексы "ж" и "с" означают, что физические свойства выбираются по средней температуре жидкости и стенки соответственно. За определяющую температуру принимать среднюю температуру жидкости tж (рабочей среды). б) турбулентный режим течения, когда 1/d >50; Reж=1·104...5·106
5.2.2. Формулы для газообразных сред при различных характерах течения в трубах: а) ламинарный режим течения
б) турбулентный режим течения
5.3. Расчёт коэффициента теплоотдачи от рабочей среды к внутренней поверхности горловиныРасчёт коэффициента теплоотдачи в горловине задвижки производится по формуле: aГ = K1´aП где К1 - поправочный коэффициент для проточной части и различных участков горловины (рисунок 1). Выбор поправочного коэффициента производится по таблице 1 в зависимости от параметров режима, используемой рабочей среды. Значения К1 были определены по результатам расчётно-экспериментальных исследований количественных характеристик теплообмена для задвижек. Приведены в приложении А. Значения поправочного коэффициент К1 для проточной части и зон горловины задвижки Таблица 1
I - проточная часть II - зона горловины III - зона горловины под крышкой Рис. 1. Схема определения коэффициентов теплопередачи. 5.4. Расчёт коэффициентов теплоотдачи с наружных поверхностей задвижки.5.4.1. Расчёт при нормальных условиях эксплуатации 5.4.1.1. Коэффициент теплоотдачи на наружной поверхности трубопроводной арматуры ai определяется как сумма коэффициентов теплоотдачи излучением aл и конвекцией aк. 5.4.1.2. Коэффициент теплоотдачи излучением aл рассчитывается по формуле:
где Ti = tw + 273 - абсолютная температура изделия или его ступени (части), К; Та = tа + 273 - абсолютная температура окружающей среды, К; t a - температура окружающей среды, °С; tw- средняя температура изделия или ступени, °С. Определяется как 1/2 суммы температур рабочей и окружающей сред (по экспериментальным или расчётным аналогам); с - 5,77 Вт/м2К - коэффициент лучеиспускания абсолютно чёрного тела; e - степень черноты излучающего тела. 5.4.1.3. Коэффициент теплоотдачи конвекцией aĸ рассчитывается по формулам:
где bi -характерный или определяющий размер (высота или диаметр), м; g = 9,8 м/с2 - ускорение силы тяжести; b = 1/ta+273 - коэффициент объёмного расширения воздуха, 1/град. 5.4.1.4. Коэффициенты "с" и "n" в формуле (10) являются функцией Gr·Pr и выбираются в зависимости от численного значения этого аргумента и положения арматуры. Их значения приведены в таблице 2. Таблица 2
5.4.1.5. Значения физических параметров, входящих в формулы (9), (10), (11) выбираются по температуре окружающей среды tа, из справочной литературы [1, 2, 5... 8]. 5.4.1.6. Для ступеней с характерным размером теплообменной поверхности "b" от 10 до 200 мм коэффициент aк, можно брать из графиков на рисунках 2, 3, 4, 5, 6. При вертикальном расположении арматуры «b» - это высота ступени или изделия; при горизонтальном - диаметр. 5.4.1.7. Коэффициент теплоотдачи с торцевой поверхности арматуры определяется так же, как и ai (7). Полученное по формуле (9) значение коэффициента теплоотдачи конвекцией необходимо увеличить на 30%, если теплоотдающая поверхность обращена кверху, и уменьшить на 30%, если вниз. В качестве определяющего размера для коэффициента теплоотдачи с торцевой поверхности aT берётся наружный диаметр торцевой части изделия. 5.4.2. Расчет коэффициентов теплоотдачи с наружных поверхностей задвижки в условиях аварии (АР). Возможны 2 варианта теплоотдачи: 5.4.2.1. Теплоотдача при турбулентном режиме, когда скорость движения паровоздушной смеси (окружающая среда) равна ω £ 20 м/с, рассчитывается по зависимостям: Рис. 2. Зависимость αк от ν при температуре воздуха ta для комплекса Gr·Pr от 5·102 до 2·107; с = 0,76; n = 0,25 (при вертикальном расположении конструкции). Рис. 3. Зависимость aк от n при температуре воздуха ta =50ºС для комплекса Gr·Pr от 5·102 до 2·107; с = 0,76; n = 0,25 (при вертикальном расположении конструкции). Рис. 4. Зависимость aк от n при температуре воздуха ta = 20ºС для комплекса Gr·Pr от 5·102 до 2·107; с = 0,5; n = 0,25 (при горизонтальном расположении конструкции). Рис.5. Зависимость aк от n при температуре воздуха ta =50ºС для комплекса Gr·Pr от 5·102 до 2·107; с = 0,5; n = 0,25 (при горизонтальном расположении конструкции) . Рис. 6. Зависимость aк от n при температуре воздуха ta =20ºС для комплекса Gr·Pr от 2·107 до 1·1013; с = 0,15; n = 0,33 (при вертикальном расположении конструкции).
где (Prж/Prс)0,25- множитель, представляющий собой поправку, учитывающую зависимость физических свойств рабочей среды от температуры; Re, Pr -рассчитываются аналогично указаний раздела 5.2. За определяющую температуру tж принимается заданная температура окружающей среды, а определяющий размер Н - высота изделия. Индексы "ж" и "с" означают, что физические свойства рабочей среды (ν, Pr, λ,) выбираются по средней температуре рабочей среды Тж и температуры внутренней поверхности стенки Тс. 5.4.2.2. Теплоотдача при естественной конвекции при температуре окружающей среды (паровоздушной смеси) 100< ta < 180°С производится по зависимостям раздела 4 и 5; коэффициент теплоотдачи конвекцией при tа = 100°C и ta = 180ºC может быть определён из графиков на рисунках 7, 8. 5.5. Расчёт скорости изменения температуры рабочей среды в расчётной зоне теплообменаРасчет производится по формуле:
где bП - скорость изменения температуры рабочей среды в проточной части (на входе) арматуры, К/с . Задается в техническом задании на проведение теплового расчёта; К2 - поправочный коэффициент для различных зон проточной части и горловины. Выбор поправочного коэффициента К2 производится по таблице 3 в зависимости от скорости изменения температуры рабочей среды в проточной части задвижки и вида среды. 6. ТЕПЛОВАЯ МОДЕЛЬТепловая модель задвижки составляется согласно выбранному методу расчёта. При использовании метода конечных разностей соблюдать рекомендации, содержащиеся в РД 26-07-25-97 "Методика расчёта температурных полей трубопроводной арматуры". Тепловые модели двух типов задвижек (для газа и жидкости) представлены на рисунках 9 и 10. Таблица 3Значения поправочного коэффициента К2 для проточной части и зон горловины задвижки
Рис.7. Зависимость коэффициента теплоотдачи aк с наружных поверхностей арматуры при свободно-конвективном теплообмене с окружающей средой (паром) от избыточной температуры ν. Температура окружающей среды ta = 100ºC. Рис. 8. Зависимость коэффициента теплоотдачи aк с наружных поверхностей арматуры при свободно-конвективном теплообмене с окружающей средой (паром) от избыточной температуры n. Температура окружающей среды ta = 180ºC. Рис. 9. Тепловая модель задвижки Л13189-200 (Рабочая среда - жидкость) Рис. 10. Тепловая модель задвижки Л13149-200 (Рабочая среда - газ)
ПРИЛОЖЕНИЕ А (рекомендуемое)РЕЗУЛЬТАТЫ РАСЧЁТНО-ЭКСПЕРИМЕНТАЛЬНЫХ ИССЛЕДОВАНИЙ КОЛИЧЕСТВЕННЫХ ХАРАКТЕРИСТИК ТЕПЛООБМЕНА ДЛЯ ЗАДВИЖЕКВ результате экспериментально-расчётных работ по исследованию теплового состояния задвижек, в том числе и количественных характеристик теплообмена, были получены значения коэффициента теплоотдачи по периметру задвижки [СТП 07.81-634-96]. Значения коэффициентов теплоотдачи представлены на рисунке А.1 и в таблице А.1. В работе [3,4] исследовалось изменение скорости разогрева рабочей среды в проточной части задвижки bпр и в горловине задвижки bГ (график на рисунке А.2). Были предложены соотношения К2 = 0,2...0,5. Очевидно, что выполнение тепловых расчетов задвижек с заданием тринадцати значений коэффициентов теплоотдачи в горловине aг по рисунку A.1 крайне трудоёмко и часто нецелесообразно. Практика выполнения тепловых расчётов трубопроводной арматуры, в том числе и задвижек, показала возможность унифицировать ранее полученные значения К1 и К2. В таблице А.2 и A.3 представлены результаты тепловых расчётов задвижек Л13189-200 (рисунок 9) и Л13149-250 (рисунок 10). Расчёты выполнены методом конечных разностей с заданием различных коэффициентов местной теплоотдачи aг, и скоростей разогрева bГ в пределах максимального и минимального значений, рекомендованных в работах [СТП 07.81-634-96; 3,4]. Условия и параметры тепловых режимов представлены в таблицах А.2 и А.З. Для задвижки по чертежу Л13189-200, предназначенной для жидкой рабочей среды, и для задвижки по чертежу Л13149-250, предназначенной для газа, были выполнены расчёты в девяти тепловых режимах. За критерий в оценке влияния изменения принимаемых количественных характеристик теплообмена К1 и К2 на температурное поле задвижек были выбраны значения максимальных перепадов температур во фланцевом соединении: D t4-2 -радиальный перепад, D t1-.3 - осевой перепад, D tст - перепад по стенке горловины. На основании выполненных тепловых расчётов и данных таблиц А.2 и А.З можно сделать следующие выводы: Рис. А.1 Схема определения коэффициентов теплоотдачи Таблица А.1
1. Для задвижек, предназначенных для газа, в пределах изменения скорости разогрева 150 £ bпр £ 300ºС/ч и значения коэффициента теплоотдачи в проточной части a < 3500 вт/м2 ºС принимать одинаковые значения К1=0,5 в горловине (зона цилиндрической образующей) и под крышкой (таблица 1 и рисунок 1). 2. Для задвижек, предназначенных для газа, в пределах изменения скорости разогрева. 10 £ bпр£ 100° С/с и значения коэффициента теплоотдачи aпр < 3500 вт/м2 °С принимать К1=1,4 в горловине (зона цилиндрической образующей) и К1= 0,5 -под крышкой (таблица 1). 3. Для задвижек, предназначенных для жидкой рабочей среды, независимо от скорости разогрева рабочей среды (75 £ bпр£ 300° С/ч) при aпр > 3500 вт/м2 °С, принимать К1 = 0,5 по всей горловине (таблица 1). 4. Для задвижек, предназначенных для жидкой рабочей среды, в пределах изменения скорости разогрева 150 £ bпр£ 300°С/ч и скорости 10 £ bпр £ 100ºC/c принимать К2=1 (таблица 3) по всей зоне горловины. 5. Для задвижек, предназначенных для газа, при малых скоростях разогрева рабочей среды (в пределах 150 £ bпр£ З00ºС/ч) принимать в горловине К2 = 0,5 (график на рисунке А.З). 6. Для задвижек, предназначенных для газа, при скоростях разогрева рабочей среды в пределах 10 £ bпр£ 100°С/с принимать в горловине К2 = 1,0 (график на рисунке А.4). 7. Для задвижек, имеющих изоляцию корпуса (до бугельной стойки), значения количественных показателей теплообмена К1 и К2 в горловине принимать равным K1 = 1,0 и K2 = 1,0 независимо от скорости разогрева и значения коэффициента теплоотдачи в проточной части Рис. А.2. Изменение температуры рабочей среды в проточной части tж.п. и горловине tж.г. макета клиновой задвижки DN300 в режиме горячего термоудара 240ºС®10ºС. Значения максимальных перепадов температур во фланцевом соединении задвижки Л 13189-200 и по стенке горловины в зависимости от значений количественных характеристик теплообмена К1 и К2 Таблица А.2
Значения максимальных перепадов температур во фланцевом соединении задвижки Л13149-250 и по стенке горловины в зависимости от количественных характеристик теплообмена Таблица А.3
Рис. А.3. Изменение максимальных перепадов температур D t1-3 , D t4-2 во фланцевом соединении и по стенке горловины D tст. задвижки Л131189-200 в режимах по табл. 5.
Рис. А.4. Изменение максимальных перепадов температур D t1-3 , D t4-2 во фланцевом соединении и по стенке горловины D tст. задвижки Л1349-250 в режимах по табл. А.3. ПРИЛОЖЕНИЕ Б (справочное)1. М.А. Михеев. Основы теплопередачи. М., Энергия. 1977. 2. В.П. Исаченко и др. Теплопередача. М., Энергоиздат. 1981. 3. Исследование теплового состояния и разработка методики расчёта запорных задвижек Ду от 100 до 800мм для АЭС. Отчёт о НИР, В.И. Лебедевич, № Р01850020304, УДК621.646.5.001.24:536.3, Ленинград, 1988. 4. Исследование и расчет теплонапряжённого состояния трубопроводной арматуры в нестационарных режимах. Сборник научных трудов, ЦКБА, Ленинград,1987. 5. И.Б. Варгафтик. Справочник по теплофизическим свойствам газов и жидкостей. М., Наука. 1972. 6. B.C. Чиркин. Теплофизические свойства материалов ядерной техники. М. Атомиздат. 1968. 7. ВТИ им. Дзержинского. Министерство энергетики и электрификации СССР. Физические свойства сталей и сплавов, применяемых в энергетике. М.-Л., Энергия. 1967. 8. Д.Ф. Гуревич, В.В. Ширяев, И.Х. Пайкин. Арматура атомных электростанций. М., Энергоиздат. 1982. ЛИСТ РЕГИСТРАЦИИ ИЗМЕНЕНИЙ
|
|