Перейти к содержимому
UzScite
  • НСИ
    • Новости События
    • Методическая информация
    • Нормативные документы
  • Каталог журналов
  • Указатель авторов
  • Список организаций

Численный метод для решения нелинейных уравнений трубопроводного транспорта реального газа

Хужаев И.К.

Мамадалиев Х.А.

Вестник ТУИТ

  • № 3 (51) 2019

Страницы: 

68

 – 

79

Язык: русский

Открыть файл статьи
Открыть страницу статьи в Интернет

Аннотация

Уравнения сохранения импульса и массы, используемые при моделировании процесса трубопроводного транспорта реальных газов, содержат члены с третьим порядком неизвестных. В инженерной практике ограничиваются использованием известных формул, которые получены для стационарного режима течения. Нестационарные задачи решаются обычно линеаризацией уравнений. Для адекватного описания объекта требуется более полный учет факторов, определяющих потери давления в трубопроводе. Степень неизвестных в уравнениях сохранения импульса и массы понижена введением вспомогательной функции в виде натурального логарифма от приведенной плотности. Применением подхода бегущих волн составлены автономные уравнения для прямых и обратных волн, в которые входят члены нелинейного закона сопротивления и силы гравитации. Данный подход учитывает наличие двух масштабов скорости: обычная гидродинамическая скорость газа и скорость распространения малых возмущений, значения которых отличаются на один порядок. Разработан численный метод решения уравнений трубопроводного транспорта, сверхсжимаемого газа с учетом силовых факторов: силы сопротивления в рамках квадратичного закона сопротивления, силы гравитации при постоянном уклоне трассы от горизонта, локальной и конвективной составляющих силы инерции газа. Для случая задания входного давления и выходного массового расхода использован итерационный процесс, позволяющий определить неизвестные граничные и внутренние значения газодинамических параметров. Для удовлетворения граничных условий постоянства массового расхода реализован метод касательных Ньютона. Исходя из особенностей составленного трансцендентного уравнения, составлена рекуррентная формула. Согласно предложенного численного метода разработаны алгоритм и программное средство в среде Pascal ABC. Проведены серии расчетов при различных длинах и уклонах участка, входного давления и выходного массового расхода газа. В настоящей работе ограничились анализом результатов вычислительного эксперимента, относящихся к путевому изменению скорости потока при различных уклонах трассы газопровода,которые демонстрировали особенности постперевального режима течения газа.Метод и алгоритм расчета допускают адаптацию для различных вариантов начальных и граничных условий, а также путевого изменения уклона трассы трубопровода.

Табиий газларни қувур орқали узатиш жараёнини моделлаштиришда фойдаланиладиган импульс ва масса сақланиш тенгламалари учинчи тартибли номаълумлардан ташкил топган.Муҳандислик амалиёт ишларида оқимнинг стационар режими учун олинган формулалардан фойдаланиш билан чегараланмоқда. Ностациоар масалалар одатда тенгламаларни чизиқлилаштириб ечилади. Объектни адекват ифодалаш учун қувурда босим йўқотишларини ифодаловчи омилларни янада тўлиқ ҳисобга олиш талаб этилади.Импульс ва масса сақланиш тенгламаларида номаълумлар даражаси зичлик билан берилган натурал логарифм кўринишидаги ёрдамчи функция киритиш билан пасайтирилган. Қаршилик қонуниятлари ва гравитация кучларининг ночизиқли ҳадлари билан ифодаланувчи тўғри ва тескари тўлқин автоном тенгламалари югурувчи тўлқинлар ёндашувини тадбиқ этиб тузилган. Бундай ёндашув тезликнинг икки масштаби мавжуд бўлганида ҳисобга олинади: қийматлари битта тартибда фарқланувчи оддий гравитацион тезлик билан кичик кўчишларнинг тарқалиш тезлиги.Куч омиллари: қаршиликнинг квадратик қонуни доирасида қаршилик кучлари, горизонтга нисбатан ўзгармас бурчак остида бўлганда гравитация кучлари, газ инерция кучларининг локал ва конвектив ташкил этувчиларини ҳисобга олган ҳолда ўта сиқилувчан газларни қувур орқали ташиш тенгламаларини ечишнинг сонли усули ишлаб чиқилган.Киришдаги босим ва чиқишдаги масса сарфи берилган ҳоллар учун газодинамик кўрсаткичларнинг номаълум чегаравий ва ички қийматларини аниқлашга имкон берувчи итерацион жараёндан фойдаланилган. Масса сарфи ўзгармаслик чегаравий шартини қаноатлантириш учун Ньютоннинг уринмалар усули тадбиқ этилган. Тузиб олинган транцендент тенглама хусусиятларидан келиб чиқиб рекуррент формула ишлаб чиқилган. Таклиф этилган сонли усулга асосан ҳисоблаш алгоритми ва Pascal ABC муҳитида дастурий восита ишлаб чиқилган. Газнинг кириш босими ва чиқишдаги масса сарфи берилганда участканинг турли узунликлари ва оғиш бурчаклари қийматлари учун сонли ҳисоблашлар ўтказилди. Мазкур ишда газ оқимининг довонорти режими хусусиятларини ифодаловчи қувур турли оғмаликлари учун оқим тезлигининг йўл-йўлакай ўзгаришларига бағишланган ҳисоблаш тажрибалари натижаларини таҳлил қилиш билан чекланилган.Ҳар хил бошланғич ва чегаравий шартлар, қувур трассаси оғиш бурчагининг ўзгаришларини хисоблаш учун келтирилган усул ва алгоритм мослаштириш имконини беради.

The equations of momentum and mass conservation, used in modeling the process of real gases pipeline transport, contain terms with third-order unknown variables.In practice of engineering, there are widely used known formulas that were obtained for the stationary flow. Unsteady problems are usually solved by linearizing the equations. In order to describe the object more adequate, more complete consideration of the factors are important determining the pressure loss in the pipeline.The degree of the unknowns in the equations of momentum and mass conservation is reduced by introducing an auxiliary function in the form of the natural logarithm of the reduced density. By using the traveling-wave approach, autonomous equations for forward and backward waves are composed, which include terms of the non-linear law of resistance and gravity. This approach takes into account the presence of two speed scales: the usual hydrodynamic gas velocity and the propagation velocity of small perturbations whose values differ by one order of magnitude.A numerical method was developed for solving the equations of supercompressible gas pipeline transport taking into account such power factors as quadratic law of resistance, gravity forces with a constant slope of the route from the horizon, and local and convective components of the gas inertia force.For the case of specifying the input pressure and output mass flow, an iterative process was used, which allows us to determine the unknown boundary and internal values of the gas-dynamic parameters. To satisfy the boundary condition of the mass flow constancy, the Newton tangent method was implemented. Based on the characteristics of the compiled transcendental equation,a recurrent formula was obtained.According to the proposed numerical method, an algorithm and software were developed using Pascal ABC environment. There was conducted a series of calculations with different lengths and slopes of the pipeline section, inlet pressure and output mass flow rate of gas. In the present work, we restricted ourselves to an analysis of the results of a computational experiment relating to the path change in the flow velocity at various slopes of the gas pipeline, which demonstrated the features of the post-cross mode of gas flow.The method and the algorithm for calculation allow adaptations for different variants of initial and boundary conditions, as well as the route variation of the slope of the pipeline.

Список использованных источников

  1. Селезнев В.Е., Алешин В.В., Прялов С.Н. Математическое моделирование магистральных трубопроводных систем: дополнительные главы. – М.: МАКС Пресс, 2009. – 356 с.
  2. Чарный И.А.Неустановившееся движение реальной жидкости в трубах.–М.:Недра,1975.–296 с.
  3. Figueiredo A.B. et al. Numerical simulation of stratified-pattern twophase flow in gas pipelines using a two-fluid model // International Journal of Multiphase Flow. – 2017, Vol. 88, PP.30-49. doi https://doi.org/10.1016/j.ijmultiphaseflow.2016.09.016
  4. He G.etal. Numerical simulation of vapor condensation in gas-water stratified wavy pipe flow with varying interface location //International Journal of Heat & Mass Transfer. Dec2017 Part B, Vol. 115, PP.635-651. DOI: 10.1016/ j.ijheatmasstransfer.2017.08.069
  5. Krivovichev G.V. A computational approach to the modeling of the glaciation of sea offshore gas pipeline // International Journal of Heat and Mass Transfer Volume 115, Part B, December 2017, PP. 1132-1148. https://doi.org/10.1016/j.ijheatmasstransfer. 2017.08.117
  6. Shin Ch. A numerical study on the characteristics of transient flow in a pressure regulator resulting from closure of the pressure control valve // Journal of Mechanical Science & Technology. Feb2013, Vol.27 Issue 2, PP.443-449. DOI: 10.1007/s12206-012-1257-y
  7. Ермолаева Н.Н. Математическое моделирование нестационарных неизотермических процессов в движущихся неизотермических многофазных средах: Дисс… д-ра физ.-мат. наук (05.13.18), СПб,2017. – 323 с.
  8. Xu H., Kong W., Yang F. Decomposition characteristics of natural gas hydrates in hydraulic lifting pipeline // Natural Gas Industry B Volume 6, Issue 2, April 2019, PP.159-167. https://doi.org/10.1016/j.ngib.2018.07.005
  9. Dong H. Numerical study on the thermal characteristics and its influence factors of crude oil pipeline after restart // Case Studies in Thermal Engineering Available online 26 April 2019, 100455. https://doi.org/10.1016/j.csite.2019.100455
  10. Yuan Q. Study on the restart algorithm for a buried hot oil pipeline based on wavelet collocation method // International Journal of Heat and Mass Transfer Volume 125, October 2018, PP. 891-907. https://doi.org/10.1016/j.ijheatmasstransfer.2018.04.127
  11. Hafsia Z., Elaouda S., Mishra M. A computational modelling of natural gas flow in looped network: Effect of upstream hydrogen injection on the structural integrity of gas pipelines // Journal of Natural Gas Science and Engineering Volume 64, April 2019, PP.107-117. https://doi.org/10.1016/j.jngse.2019.01.021
  12. Hafsia Z.,Mishra M.,Elaouda S.Hydrogen embrittlement of steel pipelines during transients Procedia Structural Integrity Volume 13,2018, PP. 210-217. https://doi.org/10.1016/j.prostr.2018.12.035
  13. Yang J.etal. Numerical Simulation of Sudden Gas Pipeline Leakage in Urban Block // Energy Procedia Volume 105, 2017, PP. 4921-4926.https://doi.org/10.1016/j.egypro.2017.03.1049
  14. Demirci H.E. et al. Experimental and numerical modelling of buried pipelines crossing reverse faults // Soil Dynamics and Earthquake Engineering Volume 114, November 2018, PP. 198-214. https://doi.org/10.1016/j.soildyn.2018.06.013
  15. Mamadaliev X. A., Khujaev I. Q., Boltibaev Sh.K. Modelling the Propagation of mass consumption waves in the Pipeline with Damper of pressure Disturbances //Florence (Italy), International Journal of Sciences and Research.Vol.74 No.8.1 Aug 2018. DOI:10/21506/j.ponte 2018.8.12.–PP. 163-170.
  16. Mamadaliev X.A., Khujaev I.Q. Mathematical model of the pipeline connected to the ends of an area with dampers of pressure American Journal of Mathematical and Computational Sciences, 2016; 1(1) РP. 43-49. http://www.aascit.org/journal/ajmcs.
  17. Бозоров О.Ш., Маматкулов М.М. Аналитические исследования нелинейных гидродинамических явлений в средах с медленно меняющимися параметрами. – Ташкент, ТИТЛП, 2015. – 96 с.
  18. Калиткин Н.Н. Численные методы. – М.: Наука, 1978. – 512 с.
  19. Самарский А.А., Вабищевич П.Н. Вычислительная теплопередача. – М.: Едиториал УРСС, 2003. – 784 с.
  20. Садуллаев Р. и др. Расчёт магистрального газопровода с учётом рельефа местности // Газовая промышленность, Москва, 2003,№8. – С. 58-59.

Список всех публикаций, цитирующих данную статью

Copyright © 2025 UzScite | E-LINE PRESS