Главная » Просмотр файлов » Н.А. Боголюбов - ОММ Методическое пособие

Н.А. Боголюбов - ОММ Методическое пособие (1133432), страница 5

Файл №1133432 Н.А. Боголюбов - ОММ Методическое пособие (Н.А. Боголюбов - ОММ Методическое пособие) 5 страницаН.А. Боголюбов - ОММ Методическое пособие (1133432) страница 52019-05-12СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

По колонкеслева направо движется поток жидкости q . Газ переносит растворенное внем изучаемое вещество - компонент , осаждающийся на поверхности гранулсорбента.qΔx0xu(x, t) - концентрация компонента в межзерновом пространстве( количество компонента в межзерновом пространстве на единицу объема колонки) .a(x, t) - концентрация компонента , сорбированного на поверхности гранул( количество сорбированного компонента на единицу объема колонки).70Функция a = φ(u) , связывающая концентрацию компонента в фазесорбента с его концентрацией в фазе раствора при условия их равновесия(равновесные концентрации ) называется изотермой сорбции.С ростом концентрации компонента количество «свободных мест» наповерхности сорбента ( куда могут осаждаться его молекулы ) убывает икривая φ(u) выходит на насыщение :В случае малых концентрацийрассматривают изотерму Генри:a = Γu ,Γ - константа, характеризующая наклоннашей кривой φ(u) на начальном участке;φuв случае значительных концентраций:Kua = φ(u) = 1 - изотерма Ленгмюра.u+KТак как изотрема является взаимно однозначной функцией , мы можемввести изотермы, обратные изотермам Генри и Ленгмюра :Обратная изотерма:Обратная изотерма Генри:u = ψ (a) ,au = ψ (a) = ,ΓKaОбратная изотерма Ленгмюра: u = ψ (a) =.K1 − a71Рассмотрим некоторый малый участок Δx и запишем уравнениебаланса рассматриваемого компонента за время Δt :x+Δxa η, t + Δt) + u(η, t + Δt) − a(η, t) − u(η, t)]dη =∫ [ (xt+Δt=∫q[u(x, τ) − u(x + Δx, t)]dτ .tПусть функции u(x, t) и a(x, t) имеют непрерывные частныепроизводные первого порядка.Тогда поделим полученное равенство на ΔxΔt и устремляя Δx и Δt кнулю, мы получим уравнение баланса вещества в дифференциальном виде:at + ut + qux = 0 .Сорбция вещества происходит не мгновенно, а за определенное время.Скорость изменения концентрации сорбированного веществапропорциональна разности текущей и равновесной концентраций.Дополнив полученное уравнение, уравнением кинетики и задавначальные и граничные условия, мы получаем постановку задачи сорбции:at + ut + qux = 0,at = β(u − ψ (a)),at=0=uux=0= u0 .t=0β - коэффициент кинетики.= 0,Для того, чтобы упростить решение поставленной задачи, нам будетудобно ввести понятие локального времени:xτ =t−.q72Значение τ « запаздывает» относительно t на время переноса веществаот границы до рассматриваемой точки x .Если мы рассмотрим u , как функцию переменных τ , x ,то :∂u(τ, x)∂u(τ, x)∂u dτ∂u dτ∂u+q=+q+=( ∂τ d x∂t∂x∂τ dt∂x )=∂u∂u1∂u∂u+q−+=q.()∂τq∂x )∂x( ∂τЗначит: при переходе к локальному времени , наше уравнение упрощается:ut + qux = qux .Введение понятия локального времени, стандартный способ упрощения решенияподобных задач по переносу вещества.Перейдем к решению поставленной задачи сорбции:введем новые переменные через старые переменные:xxτ =β t−, ξ=β .(q)qТогда область x > 0, t > 0 переходит в область ξ > 0, τ > − ξ .Рассматриваемая задача преобразуется в систему следующего вида:aτ + uξ = 0,aτ = u − ψ (a),aτ=−ξuξ=0=uτ=−ξ= 0,= u0 .Продифференцируем второе уравнение полученной системы по ξ иполучим:73aτξ + aτ + ψ′(a)aξ = 0,aaττ+ξ=0= 0,τ+ξ=0= (u − ψ (a))τ+ξ=0= 0.Будем искать решение нашей задачи в заштрихованной области−ξ ≤ τ ≤ 0 : τξM0ABnτnξnτ+ξ =0Для любой фиксированной функции a , мы можем рассматриватьпроизводную обратной функции ψ′(a) , как заданный переменныйкоэффициент k.В тоже время, на прямой τ + ξ = 0 выполняется следующее условие:ππaτ = ancos+ alcos,(4)(4)где l - вектор вдоль прямой , а n - нормаль к ней.74С учетом граничных условий рассматриваемой задачи, мы получаем ,= 0.что anξ+τТаким образом, в заштрихованной области, мы имеем следующуюзадачу с данными на характеристике:aξτ + kaξ + aτ = 0,{a = an = 0, при τ + ξ = 0 .Применяя обобщенную формулу Даламбера мы приходим к выводу:значение a в любой точке M0 , принадлежащей заштрихованной области,определяется заданными ( нулевыми в данном случае ) условиями наотрезке AB - основании характеристического треугольника M0 AB .Следовательно, a равно нулю во всей заштрихованной области.По непрерывности a = 0 при τ = 0 , а на границе ξ = 0 выполняетсяследующее условие:aτ = u0 − ψ (a), aτ=0 = 0 .Обозначим решение полученной задачи через Q(τ) и получаем вобласти ξ > 0, τ > 0 следующую задачу:aτξ + aτ + ψ′(a)aξ = 0,{aτ=0= 0, aξ=0= Q(τ) .Решая полученную задачу сорбции , мы рассмотрим и сравним накачественном уровне линейный и нелинейный случаи данной задачи.В линейной задаче сорбции мы будем исследовать задачу сорбции сизотермой Генри , в нелинейном случае - задача сорбции с изотермойЛенгмюра.75a) Линейный случай задачи сорбции ( изотерма Генри) .Итак, изотерма Генри φ(u) = Γu , и обратная ей изотерма ψ (a) =aΓсоответственно.Подставляем выписанную изотерму в полученную нами ранее задачу сданными на характеристиках ( в решение Q(τ) ) .Получаем задачу Гурса для линейной задачи сорбции:aτξ + aτ +aaξΓ= 0,= 0, aτ=0−Γ=Γu1−e0().ξ=0τЧтобы упростить решение, как и в случае с локальным временем , введемτновую переменную θ =и подберем не ограничивая общности значениеΓu0 так , чтобы выполнялось условие Γu0 = 1 .Тогда уже в области θ > 0, τ > 0 мы имеем более упрощенныйвариант предыдущей задачи Гурса:aθξ + aθ + aξ = 0,{a= 0, aθ=0−θ=1−e().ξ=0Сделаем замену функции a = ze −(θ + ξ) , вычисляем производные z по ξи θ и подставляем их в нашу задачу с данными на характеристиках:zξθ − z = 0,{z= 0, zθ=0θ=e− 1) .(ξ=0Как мы уже хорошо знаем, решением полученной задачи являетсяфункция Римана и имеет следующий вид:76v(ξ, θ, ξ1, θ1) = I0 2((ξ − ξ1)(θ − θ1) ) .Применяем обобщенную формулу Даламбера и получаем:(zv)A + (zv)Bz = (θ1, ξ1) =eθ1−12−1zvξ − vzξ)dξ + (vzθ − zvθ)dθ =([]∫2A01[vzθ − zvθ]∫2−θ1θ1=2Be θ I0 2(∫0ξ=0dθ =eθ1−12+[zv]θ=0, ξ=0θ=θ1, ξ=02θ1+ [vzθ]∫0ξ1(θ1 − θ) dξ .)Возвращаемся к старым переменным и получаем решениерассматриваемой задачи сорбции:a(θ1, ξ1) = ze −(θ1 + ξ1) =θ1e −(θ1 − θ)e −ξ1I0 2(∫0θ1= e −ξ1 e −θ I0(2 ξ1θ )dθ .∫0a1t4 > t3t2 > t1t1t3ξξ1(θ1 − θ) dθ =)графический вид решениязадачи сорбции в различныймомент времени.Значение концентрациицелевого компонента в фазесорбента на границе ( входе всорбционную колонку) современем возрастает до 1,фронт распространениякомпонента продвигаетсявдоль колонки, расплываясь.77ξ=0dθ =б) нелинейный случай задачи сорбции ( изотерма Ленгмюра ) .k1uИтак, изотерма Ленгмюра a = φ(u) =, коэффициент k1 - этоu+kмаксимальное количество компонента в сорбированной фазе, соответствующееu → ∞ (0 ≤ a < k1) .Обратная изотерма ( обратная зависимость равновесного значения u отa имеет следующий вид:kau = ψ (a) =.k1 − aψu0 = ψ (a0)a0 k1aВычислим первую и вторую производные обратной изотермы сорбцииЛенгмюра:ψ′(a) =k1k(k1 − a)2> 0, ψ″(a) =782k1k(k1 − a)3> 0.Подставляем вид изотермы в рассматриваемую задачу и получаемквазилинейное уравнение переноса:aξτ + aτ + ψ′(a)aξ = 0 , где ψ′(a) - скорость переноса со временем значения aвдоль оси ξ .Аналитически полученное уравнение не решается!Будем рассматривать его решение на качественном уровне.Пусть профиль a , как функции ξ , имеет вид плавно убывающейфункции.Тогда u ≈ ψ (a) .Итак, ψ′(a) растущая функция a , то большие значения a переносятсябыстрее.

Фронт становится более крутым.Но из уравнения кинетики следует, что aτ = − uξ = u − ψ (a) < umax .Значит крутизна фронта ограничена. Решение выходит на режимбегущей волны a = f (ξ − vτ) , когда образуется некоторыйустановившийся профиль, переносимый вдоль оси ξ . Численноемоделирование процесса,подтверждает приведенныеaкачественные соображения.ψ′(a)ξ79Теперь определим вид профиля волны . Для этого рассмотрим виднашего решения в форме бегущей волны a = f (ξ − vτ) следующей задачина бесконечном участке −∞ < ξ < ∞:aτξ + aτ + ψ′(a)aξ = 0,{aξ=−∞= a0, a= 0.ξ=∞Подставляем заявленный вид решения в рассматриваемую задачу:−vf″ − − vf′ + ψ′( f )f′ = 0,{f (−∞) = a0, f′(−∞) = f (∞) = f′(∞) = 0 .Проинтегрируем полученное выражение от −∞ до ∞ и получим:va0 = ψ (a0) .Отсюда, легко можно выразить скорость переноса:v=ψ (a0)a0.Теперь проинтегрируем вид нашего решения от α до ∞ и получим:vf′(α) + vf (α) − ψ ( f (α)) = 0 .Подставим вычисленное значение скорости переноса и получим:a0df=ψ(f ) − f .dα ψ (a0)Полученное уравнение определяет профиль бегущей волны f (α) .80Значение правой частиуравнения равно длиневертикального отрезка междупрямой и кривой,соединяющих точку 0 c точкойav(a0, ψ (a0)) .ξdfdαотрицательна, мала при a ≈ 0 иa ≈ a0 и возрастает по модулюПроизводнаяв средней части интервала значений a .Cравнивая расмотренные : линейный и нелинейный случаи задачи сорбции легкозаметить:в линейном случае: фронт волны концентрации представляет собой продвигающийся посорбционной колонке расплывающийся профиль.В нелинейном случае образуется бегущая волна постоянного профиля ( образованиеподобных структур имеет место в нелинейных задачах).81Задачи:1) Пусть через слой 0 ≤ x ≤ l пористого сорбента в направлении оси x проходитпоток воздуха со скоростью q , переносящий газообразное химическое вещество.Пусть C(x, t) - концентрация вещества в порах, a(x, t) - концентрация наповерхности сорбента.

Эти концентрации связаны соотношением a = f (C) , где f изотерма сорбции. Ее вид зависит от свойств сорбента.Пусть при t = 0 сорбент чистый, а при t > 0 на вход подается концентрацияC(0, t) = αt .И пусть f (C) = kln(1 + C).Необходимо определить C(l, t) - концентрацию вещества на выходе из слояqсорбента. Будем предполагать, что l ≤.kαРешение:Рассмотрим некоторый участок (x, x + Δx) за период времени от t доt + Δt и запишем соответствующее ему уравнение баланса вещества:[C(x, t + Δt) + a(x, t + Δt) − C(x, t) − a(x, t)]Δx = q[C(x, t) − C(x + Δx, t)]Δt .Следовательно, процесс сорбции описывается следующем уравнением:∂∂C(C + a) + q= 0.∂t∂xУчитывая явный вид зависимости a от C , начальные и граничныеусловия, получаем следующую задачу:∂C∂t∂a∂t∂C∂x{C(x,0) = 0, C(0, t) = αt .++q= 0, a = kln(1 + C), 82Чтобы упростить поставленную задачу, перейдем к локальномувремени:xτ =t− .qПерейдем к новым переменным:xτ =t− , ξ =x .qПодставим и получим:∂∂1 ∂∂∂.=−, =∂x∂ξ q ∂τ ∂t∂τТогда для функции C(ξ, τ) мы имеем следующую задачу:∂a∂τ+q∂C∂ξ= 0, a = kln(1 + C), 0 ≤ ξ ≤ l, C(ξ, − q ) = 0, C(0, τ) = ατ, τ ≥ − q .ξξПодставляем a = kln(1 + C) в наше уравнение и получаемквазилинейное уравнение переноса:∂Ck∂C⋅+q= 0.∂τ 1 + C∂ξДанное уравнение решается уже хорошо известным нам методомхарактеристик:dτk(1 + C)dξ dC==.q0Запишем первые интегралы полученного уравнения характеристик:ξkqτ −= const, C = const .1+C83Множество S , на котором задаются дополнительные условия, состоитиз двух лучей (τ > 0, ξ = 0) и (τ < 0, ξ = − qτ) .Обозначим τ* , ξ* координаты пересечения проекции характеристики,проходящей через точку (τ, ξ, C(τ, ξ)) с множеством S .С учетом граничных условий рассматриваемой задачи, мы получаем:qτ −ξk1 + C(τ, ξ)= qτ* −ξ*k1 + Φ(τ*, ξ*)C(τ, ξ) = Φ(τ*, ξ*) .,Здесь возможны две ситуации :либо τ* ≥ 0, ξ* = 0, либо τ* ≤ 0, ξ* = − qτ* .ατ*, τ* ≥ 0, ξ* = 0, Φ(τ*, ξ*) ={0, τ* ≤ 0, ξ* = − qτ* .

а) τ* ≥ 0 :qτ −ξk1 + ατ*= qτ*,C(τ, ξ) = αt* .( c учетом граничных условий)Следовательно, проекция характеристики на плоскость ξ , τопределяется соотношением:τdξ q= (1 + ατ*) . - прямая линия.dτkqξ= τkξξ = − qτ84( линии при различныхзначениях τ ) .Значение τ* определяем из полученного нами квадратного уравнения:qτ −ξk= qτ* .1 + ατ*Выражаем τ* :11τ* =τ− ±α214kξτ−+τ−(α)α(q)2.kТаким образом, при τ ≥ ξсуществует единственныйqположительный корень τ* ≥ 0 .kk1Если τ < ξ < l ≤, то оба корня полученного решенияqq αотрицательны.б) τ* < 0 :qτ − ξk = qτ ∗ − ξ ∗k,{C(ξ, τ) = 0 .Запишем наклон проекции характеристик:dξ q=.dτkКак видно из приведенного выше рисунка, характеристикипредставляют собой параллельные линии, расположенные ниже линииqξ= τ.k85Итак, мы нашли решение при всех значениях τ ∗ и возвращаясь кстарым переменным, можем записать ответ:C(x, t) =12α(t − qx ) − 1 +0, t ≤1+kqα(t − qx ) + 1 −()24kxαq, t ≥1+kx,qx.2) Пусть по реке вместе с водой переносятся сброшенные в реку отходы химическогопроизводства.

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

Тип файла
PDF-файл
Размер
7,78 Mb
Тип материала
Высшее учебное заведение

Список файлов книги

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