Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

ОСОБЕННОСТИ ЗАДАНИЯ ГРАНИЧНЫХ УСЛОВИЙ ТРЕТЬЕГО РОДА НА ПОВЕРХНОСТИ МНОГОЛЕТНЕМЕРЗЛОГО ГРУНТА ПРИ МОДЕЛИРОВАНИИ ЕГО ТЕПЛОВОГО ВЗАИМОДЕЙСТВИЯ С «ГОРЯЧИМ» ТРУБОПРОВОДОМ В ПРОГРАММЕ COMSOL MULTIPHYSICS 4.3b

Марков Е.В. 1 Пульников С.А. 1 Казакова Н.В. 1
1 ФГБОУ ВО «Тюменский государственный нефтегазовый университет»
Создана математическая модель взаимодействия многолетнемерзлого грунта и трубопровода в программном комплексе Comsol Multiphysics 4.3b. Проведен анализ исходных данных для моделирования. Сделано заключение о том, что граничные условия третьего рода на поверхности грунта зависят от времени. Обозначены основные требования, предъявляемые к регрессионной модели, описывающей граничные условия на поверхности грунта. Разработана методика аппроксимации исходных данных рядами Фурье. Определено оптимальное число слагаемых в ряду Фурье, обеспечивающее плавное изменение аппроксимируемой величины при переходе от одного месяца к другому. При помощи инструментов пакета Matlab с интегрированным языком символьного исчисления MuPad сформировано выражение, задающее граничные условия в Comsol Multiphysics 4.3b. Выполнена проверка соответствия использованных символьных обозначений функций и переменных смысловому содержанию тех же символьных обозначений в интерпретируемом языке Comsol Multiphysics 4.3b.
горячий трубопровод
многолетнемерзлый грунт
ореол оттаивания
переменные граничные условия
климатические параметры
солнечная радиация
граничные условия третьего рода
1. Алескерова З.Ш., Пульников С.А., Сысоев Ю.С., Казакова Н.В. Оценка эффективности геотехнического мониторинга магистральных газопроводов по качеству продуцируемой информации // Известия высших учебных заведений. Нефть и газ. – 2015. – № 3. – С. 81–86.
2. Владимиров В.С. Уравнения математической физики. – 4-е изд. – М.: Наука. Главная редакция физико-математической литературы, 1981. – 512 с.
3. Горковенко А.И., Вагнер В.В., Пульников С.А. Результаты исследований взаимодействий вибронагруженных магистральных подземных газопроводов с окружающими грунтами // Известия высших учебных заведений. Нефть и газ. – 2007. – № 4. – С. 73–77.
4. Карнаухов М.Ю., Лазарев С.А., Пульников С.А., Сысоев Ю.С. Исследование процесса формирования эксплуатационного положения протяженных участков МГ Уренгой – Челябинск в сложных гидрогеологических условиях // Газовая промышленность. – 2015. – № S724 (724). – С. 53–58.
5. Кушнир С.Я., Карнаухов М.Ю., Пульников С.А., Сысоев Ю.С. Анализ пространственных перемещений магистральных газопроводов с определением граничных зон // Известия высших учебных заведений. Нефть и газ. – 2011. – № 5. – С. 72–75.
6. Кушнир С.Я., Пульников С.А., Серебренников А.А., Сенив Д.М. Устойчивость призм обвалования подземных магистральных нефтегазопроводов // Известия высших учебных заведений. Нефть и газ. – 2010. – № 4. – С. 65–70.
7. Кушнир С.Я., Пульников С.А., Сысоев Ю.С. Пространственная устойчивость подземного магистрального газопровода на обводненных участках трассы // Известия высших учебных заведений. Нефть и газ. – 2012. – № 1. – С. 72–76.
8. Луканин В.Н., Шатров М.Г. Теплотехника. – М.: Высшая школа, 2000. – 671 с.
9. Марков Е.В., Пульников С.А., Гербер А.Д. Проблемы задания граничных условий при моделировании теплового взаимодействия «горячих» трубопроводов с многолетнемерзлыми грунтами // Фундаментальные исследования. – 2015. – № 2–10. – С. 2106–2110.
10. Марков Е.В., Пульников С.А., Сысоев Ю.С. Анализ типовых методик расчета глубины протаивания вечномерзлых грунтов под трубопроводами // Современные проблемы науки и образования. – 2015. – № 1. – С. 142.
11. Мерданов Ш.М. Совершенствование методов строительства и ремонта трубопроводов в условиях Крайнего Севера: дис. ... канд. техн. наук. – Тюмень, 1996.
12. Пульников С.А. Взаимодействие вибронагруженных магистральных газопроводов с окружающими грунтами: дис. ... канд. техн. наук. – Тюмень, 2007. – 173 с.
13. Пульников С.А., Сысоев Ю.С., Карнаухов М.Ю., Лазарев С.А. Проблемы внедрения систем геотехнического мониторинга на объектах газотранспортной системы ОАО «Газпром» // Геотехника. – 2014. – № 5/6. – С. 77–83.
14. Пульников С.А., Сысоев Ю.С., Гербер А.Д., Карнаухов М.Ю. Методика обработки данных мониторинга пространственного положения эксплуатируемого участка подземного магистрального газопровода // Проблемы сбора, подготовки и транспорта нефти и нефтепродуктов. – Уфа: Изд-во ИПТЭР, 2013. – Вып. 3 (93). – С. 58–66.
15. Тарасенко А.А., Сильницкий П.Ф., Тарасенко Д.А. Проектирование теплоизолированных фундаментов резервуаров на вечной мерзлоте // Известия высших учебных заведений. Нефть и газ. – 2012. – № 2. – С. 84–88.

При проектировании трубопроводных систем, перекачивающих горячие продукты в условиях сплошной или островной многолетней мерзлоты, особенно важным является правильный выбор способа прокладки и обоснование толщины теплоизоляционного покрытия. При подземном способе прокладки определяющий фактор – это допустимая величина просадки оттаявшего многолетнемерзлого грунта (ММГ), зависящая от глубины его оттаивания и жесткости трубопровода [1; 3–7; 9–15].

Цель работы – разработать наиболее рациональную методику задания периодических граничных условий третьего рода на поверхности моделируемого грунтового массива, воспроизводящую среднестатистические значения основных климатических параметров в районе проектируемого трубопровода.

Задачи:

– выполнить анализ исходных данных задачи о моделировании теплового взаимодействия трубопроводов с ММГ;

– определить основные возможности программного комплекса Comsol Multiphysics 4.3b в части задания граничных условий третьего рода;

– разработать математическую формулировку граничных условий третьего рода, позволяющую использовать ее в Comsol Multiphysics 4.3b.

Для создания модели трубопровода и грунта в программе Comsol Multiphysics 4.3b необходимо выбрать физический режим Heat Transfer with Radiation in Participating Media, позволяющий моделировать теплопроводность в средах, испытывающих фазовые переходы типа «плавление – затвердевание» на основании нелинейного дифференциального уравнения теплопроводности [2]:

markov01.wmf (1)

где ρ(T) – плотность грунта; T – температура; cр(T) – теплоемкость грунта; t – время; K(T) – теплопроводность грунта; G – внутренние объемные источники тепла.

Среди решателей необходимо выбрать Time Dependent, позволяющий решать задачи с нестационарным теплообменом. Во вкладке Geometry необходимо создать геометрическую модель исходной задачи (рис. 1).

Во вкладке Materials, секции Material Contents задаются теплофизические характеристики грунтов, теплоизоляции и материала стенки трубопровода. В секции Geometric Entity Selection производится выбор расчетных областей, к которым принадлежит тот или иной материал.

pic_34.tif

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

pic_35.tif

Рис. 2. Добавление граничного условия «Heat Flux». Осуществляется вызовом контекстного меню во вкладке «Heat Transfer with Radiation in Participating Media»

Во вкладке Heat Transfer with Radiation in Participating Media необходимо добавить граничное условие Heat Flux (рис. 2), задающее тепловой поток в общем виде. Затем в секции Heat Flux в строке General inward heat flux вводится имя переменной, которая будет задавать величину теплового потока на поверхности грунта, например «Full_Heat_Flux». В секции Boundary Selection необходимо выбрать верхнюю границу расчетной области (рис. 3).

pic_36.tif

Рис. 3. Ввод переменной, задающей величину теплового потока и выбор границы № 3 (выделена синим цветом в окне Graphics)

Далее во вкладке Defenitions необходимо добавить новую таблицу с переменными и записать в столбец Name ранее определенное имя переменной, задающей тепловой поток на поверхности грунта: «Full_Heat_Flux». В столбце Expression вводится выражение в символьном виде, задающее непосредственно сам тепловой поток. Удобные возможности для выполнения операций с символьными выражениями предоставляет программный комплекс Matlab с интегрированным языком символьного исчисления MuPad.

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

markov02.wmf (2)

где A – альбедо деятельной поверхности; a – коэффициент теплоотдачи от поверхности земли; Tв – температура воздуха; ε – степень черноты поверхности земли в длинноволновом диапазоне; σ0 – постоянная Стефана – Больцмана; b – коэффициент излучения атмосферы [8]; T – температура моделируемого массива грунта.

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

Для определения мощности теплового потока Q необходимо разделить величину суммарной солнечной радиации на количество секунд в месяце. В этом случае размерность Q будет Вт/м2.

Перейдем к выводу выражения, аппроксимирующего величину Q в зависимости от времени года. Для этих целей можно использовать различные регрессионные модели. Однако выбор должен быть сделан в сторону той модели, которая позволит учесть две важные особенности:

1. Регрессионная модель должна представляться периодической функцией времени. Это обеспечит повторяемость значений основных климатических характеристик в одно и то же время для каждого расчетного года.

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

Разработанная в начале XIX века французским математиком Фурье теория сходящихся тригонометрических рядов предоставляет очень удобные возможности для вычисления выражения, аппроксимирующего Q в соответствии со всеми требованиями, приведенными выше:

markov03.wmf (3)

markov04.wmf (4)

markov05.wmf (5)

markov06.wmf (6)

где ty – продолжительность года, с; Qi – мощность теплового потока солнечной радиации в i-м месяце, Вт/м2; t0 = 0 с, соответствует началу года; ti – момент окончания i-го месяца в году, с; k – число слагаемых в ряду Фурье; t – расчетный момент времени.

Аналогичные операции выполняются для всех переменных выражения (2), зависящих от времени года. Для плавного изменения аппроксимируемой величины при переходе от одного месяца к другому рекомендуется брать не более 7–9 слагаемых ряда Фурье. После расчета коэффициентов ряда Фурье для всех переменных, окончательное выражение для Full_Heat_Flux записывается в колонку Expression.

Следует отметить, что использованная переменная «t», «T» функции «cos» и «sin» интерпретируются в Comsol Multiphysics 4.3b в соответствии со смыслом, положенным в формулах (3)–(6).

Выводы

Создана математическая модель взаимодействия многолетнемерзлого грунта и трубопровода в программном комплексе Comsol Multiphysics 4.3b. Проведен анализ исходных данных для моделирования. Сделано заключение о том, что граничные условия третьего рода на поверхности грунта зависят от времени. Обозначены основные требования, предъявляемые к регрессионной модели, описывающей граничные условия на поверхности грунта. Разработана методика аппроксимации исходных данных рядами Фурье. Определено оптимальное число слагаемых в ряду Фурье, обеспечивающее плавное изменение аппроксимируемой величины при переходе от одного месяца к другому. При помощи инструментов пакета Matlab с интегрированным языком символьного исчисления MuPad сформировано выражение, задающее граничные условия в Comsol Multiphysics 4.3b. Выполнена проверка соответствия использованных символьных обозначений функций и переменных смысловому содержанию тех же символьных обозначений в интерпретируемом языке Comsol Multiphysics 4.3b.

Рецензенты:

Якубовский Ю.Е., д.т.н., заведующий кафедрой «Прикладная механика», ФГБОУ ВО «Тюменский государственный нефтегазовый университет», г. Тюмень;

Тарасенко А.А., д.т.н., профессор кафедры «ТУР», ФГБОУ ВО «Тюменский государственный нефтегазовый университет», г. Тюмень.


Библиографическая ссылка

Марков Е.В., Пульников С.А., Казакова Н.В. ОСОБЕННОСТИ ЗАДАНИЯ ГРАНИЧНЫХ УСЛОВИЙ ТРЕТЬЕГО РОДА НА ПОВЕРХНОСТИ МНОГОЛЕТНЕМЕРЗЛОГО ГРУНТА ПРИ МОДЕЛИРОВАНИИ ЕГО ТЕПЛОВОГО ВЗАИМОДЕЙСТВИЯ С «ГОРЯЧИМ» ТРУБОПРОВОДОМ В ПРОГРАММЕ COMSOL MULTIPHYSICS 4.3b // Фундаментальные исследования. – 2015. – № 11-2. – С. 287-291;
URL: https://fundamental-research.ru/ru/article/view?id=39326 (дата обращения: 28.03.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674