Главная » Просмотр файлов » Диссертация

Диссертация (1145283), страница 18

Файл №1145283 Диссертация (Математическое моделирование нестационарных неизотермических процессов в движущихся многофазных средах) 18 страницаДиссертация (1145283) страница 182019-06-29СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 18)

Он базируется на использовании аналитического решения плоскойзадачи Стефана при соответствующих граничных и начальных условиях.1072.4. Оледенение цилиндрической поверхностив морской воде (модель Л1),численные и аналитические решения2.4.1. Математическая модель Л1, численное решениеЗапишем математическую модель оледенения внешней боковой поверхности цилиндра в морской воде.Модифицированное условие Стефана для морского льдаВ отличии от пресной воды, в соленой морской воде при нарастании льдапроисходит приток солей в слой воды, прилегающий к фронту замерзания,кроме того, теплота плавления, температура фазового перехода становятся сложными функциями солености морского льда и его температуры.

Всеэто приводит к искажению динамики оледенения, которая имела бы место впресной воде, в частности, может измениться температура замерзания, увеличиться плотность прилегающей воды и возникнуть дополнительные конвективные течения. Суммарный эффект, как отмечается в большинстве исследований, выражается в замедлении процесса оледенения по сравнению спроцессом оледенения в пресной воде.В настоящей работе эти эффекты учитываются введением в условие Стефана дополнительного притока тепла к фронту оледенения. В модели принимается линейная зависимость этого дополнительного притока тепла от скорости нарастания льда. Гипотеза о линейности, как и оправданность введениядополнительного притока тепла, должны быть подтверждены (или опровергнуты) сопоставлением расчетов по модели Л1 с экспериментальными данными.Итак, введем в условие Стефана дополнительное слагаемое α dy , проdtпорциональное скорости нарастания льда.

Общий поток тепла q от воды кфронту оледенения положим равным сумме потока тепла от воды q3 , который был бы при отсутствии этого эффекта, и дополнительного притока теплаα dy :dtdyq = q3 + α ,dt108α — эффективный параметр модели, имеющий размерность (кДж/м3 ) удельного объемного источника тепла. О выборе величины потока тепла q3 говорилось выше в замечании 1 параграфа 2.2. В пресной воде эффективныйпараметр α равен нулю.Модифицированное условие Стефана записывается следующим образом:dydy∂T ∂T − q = γρ− q3 = (γρ + α) .⇒ λ (2.16)λ ∂r R+y(t)dt∂r R+y(t)dtАрифметически введение в условие Стефана дополнительного притока теплаα dy эквивалентно замене эффективной теплоты плавления морского льда γdtна некую новую величину γ∗ , равную: γ∗ = γ + αρ .

Однако значение величиныγ∗ может оказаться необъяснимым с физической точки зрения, тогда как впредлагаемой модификации условия Стефана явно отражено, что увеличениемножителя (γρ + α) перед скоростью нарастания льда может быть связаноне только с механизмами, которые приводят к увеличению эффективной теплоты плавления морского льда γ. Параметр α является эффективным параметром модели, учитывающим, в частности, и осреднение по углу, посколькуодномерная модель предполагает аксиальную симметрию процессов, котораяочевидно нарушается, если на динамику оледенения оказывает влияние стекание рассола в образующемся морском льде под действием силы тяжести.Выбор эффективного параметра α рассмотрен в п.

2.4.2.Один момент следует подчеркнуть. Далее показано, что для рассматриваемой одномерной однофазной задачи Стефана решение существует и единственно. Если при выбранном значении α рассчитанное по модели Л1 решениезадачи совпадает с экспериментальными данными, то, в силу единственностирешения однофазной одномерной задачи Стефана, можно считать выбор αудовлетворительным.Модель Л1λ ∂∂T=ρc∂tr ∂r∂Tr∂r,T (r, t0 ) = T 0 (r),r ∈ (R, R + y),r ∈ [R, R + y0 ];t > t0 ;(2.17)(2.18)109T (R, t) = T0 (t),t > t0 ;T (R + y(t), t) = T∗ , t > t0 ;dydy∂T − q = γρ , q = q3 + α ,λ ∂r R+y(t)dtdty t 0 = y 0 .(2.19)(2.20)t > t0 ;(2.21)(2.22)В модели Л1 обозначено: R — радиус цилиндра; ρ, λ, c, — средние по слоюплотность и коэффициенты теплопроводности и удельной теплоемкости морского льда соответсвенно; γ — средняя удельная теплота плавления морскогольда, (для пресного льда эта величина в параграфе 2.3 обозначалась Q); t, r— время и координата вдоль радиуса в цилиндрической системе координат;T = T (r, t) — температура льда; y = y(t) – толщина слоя льда на боковойповерхности цилиндра; q3 — радиальная составляющая осредненного по углувектора потока тепла от воды к фронту оледенения при заданной неизменнойсолености морской воды; α dydt — дополнительный осредненный по углу притоктепла к фронту оледенения, обусловленный перечисленными выше факторами; r = R + y(t) — координата фронта оледенения; y0 — толщина слоя льда вначальный момент времени t0 ; T0 (t) — заданный закон изменения температуры поверхности цилиндра; T 0 (r) — распределение температуры в слое льдав начальный момент времени; T∗ — температура замерзания морской воды.В модели Л1 условие (2.20) выражает неизменность температуры нафронте оледенения; (2.21) является модифицированным условием Стефана,выражающим баланс между потоками тепла на фронте оледенения и количеством тепла, выделяющимся при образовании льда; (2.17) — уравнение теплопроводности в слое льда; (2.18), (2.22) — начальные условия; (2.19) — граничное условие первого рода на поверхности цилиндра.Величины λ, c, ρ, T∗ , γ, α в модели Л1 считаются постоянными.

Еслидиапазон изменения температуры в законе T0 (t) большой, может возникнутьнеобходимость в учете зависимостей этих величин от температуры льда.О существовании и единственности классического решения однофазной одномерной задачи Стефана (модель Л1). Как отмечалосьв параграфе 2.3, в работе С.Л.

Каменомостской [58] доказано существованиеи единственность обобщенного решения задачи Стефана в общей постановке110(многомерность, многофазность, зависимость коэффициентов от температуры) при определенных требованиях на гладкость и ограниченность функций.Из этих доказательств следует, что если классическое решение задачи Стефана существует, то оно совпадает с обобщенным и, следовательно, оно единственно. В работах Ф.П. Васильева [60] и Дугласа [170] доказано существование классического решения одномерной однофазной задачи Стефана. Такимобразом, для модели Л1 решение существует и единственно при ограниченности величины потока q3 и при непрерывности и ограниченности функцийT0 (t) и T 0 (r).Оценка максимальной толщины y∗ слоя льда в модели Л1При неизменной температуре T0 (t) = T0 < T∗ поверхности цилиндратолщина слоя льда монотонно увеличивается со временем до максимальнойвеличины y∗ в равновесном состоянии при∂∂t= 0.

Модифицированное условиеСтефана (2.21) в состоянии равновесия переходит в равенствоdT − q3 = 0,λ dr R+y∗(2.23)которое совместно со стационарным вариантом уравнения (2.17) позволяетрассчитать максимальную толщину y∗ слоя льда. Из уравнения (2.17) при∂∂t= 0 следует:rdTB dr(T∗ − T0 )= B ⇒ dT =⇒ B=.drrln(1 + yR∗ )Модифицированное условие Стефана (2.23) с учетом найденного выражениядля dT преобразуется в следующее трансцендентное уравнение относительноdry∗ :(T∗ − T0 )λ− q3 = 0.(2.24)ln(1 + yR∗ )(R + y∗ )Решение уравнения (2.24) не представляет трудности, оно может быть получено методом половинного деления, методом касательных, методом хорди другими хорошо известными методами. Из уравнения (2.24) определяетсямаксимальная толщина y∗ слоя льда в модели Л1.111Если температура поверхности цилиндра убывает со временем, стремяськ постоянной величине, для оценки y∗ в качестве T0 следует взять минимальное значение температуры поверхности цилиндра.Численное решение системы уравнений модели Л1Система уравнений модели Л1 относится к классу нелинейных задач Стефана с подвижным фронтом.

Как отмечалось в параграфе 2.3, для численного решения одномерных задач стефановского типа наиболее эффективныметоды с явным выделением границы оледенения.Принцип построения схемы численного решения задачи (2.17)–(2.22), являющейся однофазной одномерной задачей Стефана, взят из работы Ф.П.Васильева [60]. Численное решение с явным выделением фронта оледенениястроится по неявной разностной схеме решения уравнения теплопроводности (2.17) с использованием итерационной процедуры решения возникающейнелинейной системы уравнений. В этом методе переменный шаг по времениτ [n+1] на (n + 1)-м временном слое определяется таким образом, чтобы слойльда за время τ [n + 1] увеличился на заданную постоянную величину h, равную шагу сетки по радиусу.

Сходимость такого метода численного решения,как и существование решения этого варианта задачи Стефана, доказаны вработе [60].Введем равномерную сетку по радиусу с шагом h, ri = R + ih, i =0, 1, . . . , N . Число N определяется как целая часть от величины (y∗ /h), гдеy∗ — максимальная толщина слоя льда при выбранных параметрах модели,которая рассчитывается из уравнения (2.24).Обозначим n — номер шага по времени; толщина слоя льда за n шагов повремени равна (y0 + nh); вводим массив шагов по времени τ [n], n = 1, . .

. , N ;τ s — величина шага по времени на (n + 1)-м временном слое в s-й итерации;после окончания итерационного процесса искомому шагу τ [n + 1] присваивается значение τ s+1 в последней итерации; tn+1 = t0 + τ [1] + τ [2] + . . . + τ [n + 1]— суммарное время нарастания льда толщиной y, (y = y0 + (n + 1)h); Tis —температура льда в i-м узле на (n+1)-м временном слое в s-й итерации, Tin+1— температура льда в i-м узле на (n + 1)-м временном слое, Tin+1 присваива-112ется значение Tis+1 в последней итерации при выполнении условия окончанияитерационного процесса.

Величины Tin , τ [n] считаются известными, величины Tin+1 , τ [n + 1] рассчитываются по следующему алгоритму.(0)Первый шаг по времени в нулевом приближении τ1 (при отсутствиислоя льда в начальный момент) может быть найден, например, следующимобразом. На первом шаге слой льда содержит два расчетных узла: i = 0 иi + 1 = 1. Положим в них температуру равной: T01 = T0 , T11 = T∗ . Модифи-цированное условие Стефана (2.21) в разностном виде записывается в видеравенства:λhT∗ − T0− q3 = (γρ + α) (0) ,hτ1из которого следует:(0)τ1 =(0)(γρ + α)h2.(λ(T∗ − T0 ) − q3 h)(2.25)Далее величина τ1 уточняется по аналогу итерационного алгоритма (пункты9)–12)), приведенного ниже.Начиная со второго шага по времени динамика оледенения рассчитывается в цикле по времени.Цикл по времени по n от 1 до N1) n := n + 1.Находится время процесса:t := 0, в цикле по i от 1 до (n − 1) выполняется:t := t + τ [i].

После этого цикла текущее время t∗ определяется равен-ством: t∗ = t + t0 .В цикле по итерациям используются следующие обозначения: τ s+1 обо-значается τ 1 , а τ s обозначается τ 0 .Перед входом в цикл по итерациям шагу по времени τ 1 присваиваетсязначение шага по времени на предыдущем временном слое:τ 1 := τ [n − 1].113Цикл по итерациям по s от 0 до s∗ , (s∗ — заданное максимальноечисло итераций.)2) s := s + 1, τ 0 := τ 1 , t := t∗ + τ 0 .Уравнение теплопроводности (2.17) аппроксимируется по неявной монотонной разностной схеме Самарского и решается методом прогонки в областиr ∈ (R, R + nh) с использованием граничных условий (2.19), (2.20).λTis − Tin−1=τ0ρ c h2hhss)(Ti+1− Tis ) − (1 −)(Tis − Ti−1) ,(1 +2ri2rii = 0, .

Характеристики

Список файлов диссертации

Математическое моделирование нестационарных неизотермических процессов в движущихся многофазных средах
Свежие статьи
Популярно сейчас
Как Вы думаете, сколько людей до Вас делали точно такое же задание? 99% студентов выполняют точно такие же задания, как и их предшественники год назад. Найдите нужный учебный материал на СтудИзбе!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6417
Авторов
на СтудИзбе
307
Средний доход
с одного платного файла
Обучение Подробнее