Сравнение гидравлических расчетов Опции

[Масюня]

Добрый вечер господа хорошие! Прошу помочь в таком вопросе, он не нов, но все еще актуален! Гидравлика газопроводов. Ручной расчет (экселька) Сравнивала я расчет по СП , по ДБН по Ионину по СНиП 66 года и по Борисову и Даточному. А так же на глаза попалась статейка одного автора.
Кто какой литературой пользуется? сравнивать начала потому что для низкого давления цифры по ДБН и СП сошлись несмотря на немного отличные формулы, а вот со средним давлением которым занимаюсь сейчас значительная разница в результатах. Оговорюсь что для меня руководящий документ это ДБН. (но в нем нашла в формулах гидравлики опечатки) поэтому доверие к нему слегка упало... В Ионине и в Борисове нет единиц измерения которые используют в конкретных формулах, поэтому слегка теряюсь... Прошу совета и помощи тех кто может помочь установить истину.

Давайте все кого этот вопрос не оставил равнодушным разберем конкретный пример, и выясним где собака зарыта! Данные беру из примера статьи Г.Г. Крицкого, так как он там приводит свое сравнение которое тоже весьма интересно и хочется обсудить https://www.politerm.com/articles/gaz/gaz/

Давайте для примера возьмем такие исходные данные:

Рн (входное избыточное) 0,3 МПа
плотность 0,717 кг/м3
Расход при н.у. 50004 нм3/ч (13,89 м3/с)
Температура газа 373.15 К (100 С)
длинна участка 500 м
диаметр 0,35м (35 см)
как пишут в статье примем лямбда 0,015

1. Какую брать величину кинематической вязкости? 14,3*10-6 степени в м2/с или высчитывать по формуле 14 из Борисов и даточный :v= дин. вязкость/плотность

От этого коэффициента пляшет число Рейнольдса,

в ДБН Rе = 0,0354 Qo/dυ = 3536746,45 (Q=50004, d=35см , u=14.3*10-6)

где Qo - расход газа, м3/ч, при температуре 0 °С и давлении 0,10132 МПа;
d - внутренний диаметр газопровода, см;
υ - коэффициент кинематической вязкости газа, м2/с (при температуре 0 °С и давлении 0,10132 МПа).

В СП 42-101-2003 более шире представлена эта формула

Q Q
0 0
Re = ────────── = 0,0354 ───────, что как пишет Крицкий коэфициент высчитан с небольшой погрешностью.
9 пи d v d v
При использовании полного значения Пи = 3533525,93 (Q=50004, d=35см , u=14.3*10-6)

в Борисове и даточном Re = скорость*диаметр/кин. вязкость = 3533525,93 (Q=50004, d=35см , u=14.3*10-6)

что за собой тянет небольшую погрешность в вычислении лябда при использовании полного значения пи лямбда = 0,014536267 при 0,0354 лямбда =0,01453606


Среднее давление
по ДБН формула Е5
2 2 -5 0,25 2 5
Р - Р / L = 1,4*10 *((n/d + 1922* ((кин. вязкость*d)/Q) ) * Q /d * плотность (0,717)= 7.1682*10-6 умножить на длинну чтоб получить разность квадратов (L=500м) = 0,003584
1 2 о о

а искомая потеря (корень из пред. числа)= 0,0598 МПа

СП

2 2 -4 2 5
Р - Р = 1,2687*10 *лямбда * Q /d * плотность *длину = 0,0032479
1 2 о

а искомая потеря (корень из пред. числа)= 0,18022 МПа

Как видите разница абалдеть....


по СНиП 66 и 62 года та же формула

Р1^2 - Р2^2 / L = 1,45*10^-3 *((n/d + 1922* ((кин. вязкость*d)/Qo) )^0.25 * Qo^2 /d^5 * уд. вес(который при переводе получается равен плотности (0,717)= 0,00654 умножить на длину чтоб получить разность квадратов =4,0576 , корень = 2,0144 кгс/см2 или 0,20144 МПа

----- добавлено через ~2 мин. -----
В Ионине не обозначена размерность используемая в формуле поэтому что в итоге там получу не знаю, кто может еще как считал?

[Масюня]

Добрый вечер господа хорошие! Прошу помочь в таком вопросе, он не нов, но все еще актуален! Гидравлика газопроводов. Ручной расчет (экселька) Сравнивала я расчет по СП , по ДБН по Ионину по СНиП 66 года и по Борисову и Даточному. А так же на глаза попалась статейка одного автора.
Кто какой литературой пользуется? сравнивать начала потому что для низкого давления цифры по ДБН и СП сошлись несмотря на немного отличные формулы, а вот со средним давлением которым занимаюсь сейчас значительная разница в результатах. Оговорюсь что для меня руководящий документ это ДБН. (но в нем нашла в формулах гидравлики опечатки) поэтому доверие к нему слегка упало... В Ионине и в Борисове нет единиц измерения которые используют в конкретных формулах, поэтому слегка теряюсь... Прошу совета и помощи тех кто может помочь установить истину.

Давайте все кого этот вопрос не оставил равнодушным разберем конкретный пример, и выясним где собака зарыта! Данные беру из примера статьи Г.Г. Крицкого, так как он там приводит свое сравнение которое тоже весьма интересно и хочется обсудить https://www.politerm.com/articles/gaz/gaz/

Давайте для примера возьмем такие исходные данные:

Рн (входное избыточное) 0,3 МПа
плотность 0,717 кг/м3
Расход при н.у. 50004 нм3/ч (13,89 м3/с)
Температура газа 373.15 К (100 С)
длинна участка 500 м
диаметр 0,35м (35 см)
как пишут в статье примем лямбда 0,015

1. Какую брать величину кинематической вязкости? 14,3*10-6 степени в м2/с или высчитывать по формуле 14 из Борисов и даточный :v= дин. вязкость/плотность

От этого коэффициента пляшет число Рейнольдса,

в ДБН Rе = 0,0354 Qo/dυ = 3536746,45 (Q=50004, d=35см , u=14.3*10-6)

где Qo - расход газа, м3/ч, при температуре 0 °С и давлении 0,10132 МПа;
d - внутренний диаметр газопровода, см;
υ - коэффициент кинематической вязкости газа, м2/с (при температуре 0 °С и давлении 0,10132 МПа).

В СП 42-101-2003 более шире представлена эта формула

Q Q
0 0
Re = ────────── = 0,0354 ───────, что как пишет Крицкий коэфициент высчитан с небольшой погрешностью.
9 пи d v d v
При использовании полного значения Пи = 3533525,93 (Q=50004, d=35см , u=14.3*10-6)

в Борисове и даточном Re = скорость*диаметр/кин. вязкость = 3533525,93 (Q=50004, d=35см , u=14.3*10-6)

что за собой тянет небольшую погрешность в вычислении лябда при использовании полного значения пи лямбда = 0,014536267 при 0,0354 лямбда =0,01453606


Среднее давление
по ДБН формула Е5
2 2 -5 0,25 2 5
Р - Р / L = 1,4*10 *((n/d + 1922* ((кин. вязкость*d)/Q) ) * Q /d * плотность (0,717)= 7.1682*10-6 умножить на длинну чтоб получить разность квадратов (L=500м) = 0,003584
1 2 о о

а искомая потеря (корень из пред. числа)= 0,0598 МПа

СП

2 2 -4 2 5
Р - Р = 1,2687*10 *лямбда * Q /d * плотность *длину = 0,0032479
1 2 о

а искомая потеря (корень из пред. числа)= 0,18022 МПа

Как видите разница абалдеть....


по СНиП 66 и 62 года та же формула

Р1^2 - Р2^2 / L = 1,45*10^-3 *((n/d + 1922* ((кин. вязкость*d)/Qo) )^0.25 * Qo^2 /d^5 * уд. вес(который при переводе получается равен плотности (0,717)= 0,00654 умножить на длину чтоб получить разность квадратов =4,0576 , корень = 2,0144 кгс/см2 или 0,20144 МПа

----- добавлено через ~2 мин. -----
В Ионине не обозначена размерность используемая в формуле поэтому что в итоге там получу не знаю, кто может еще как считал?