Лекция 7 (1133447)

Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла

Глава 3. Математическое моделированиенелинейных объектов и процессов1.Математическиемоделипроцессовнелинейной теплопроводности и горения1. Краевые задачитеплопроводностидляквазилинейногоуравненияРассмотрим квазилинейное уравнение теплопроводности:∂u ∂∂ucρ= (k (u ) ),∂t ∂x∂xгде С – удельная теплоемкость,теплопроводности.ρ(1)- плотность,k- коэффициентПараболическое уравнение (1) с коэффициентомК(u),зависящимот температуры и удовлетворяющим условию K(0)=0, называетсявырождающимся.Автомодельными решениями уравнения (1) мы будем называть такие егочастные решения специального вида, которые могут быть полученыпутем интегрирования некоторых обыкновенных дифференциальныхуравнений, аргументы искомых функций которых представляют собойкомбинацию независимых переменных x и t.Найдем автомодельноеусловиям:решениеуравнения(1)=u (0, t ) u=u21 , u ( x, 0)ПустьТогда:=ξx2удовлетворяющее(2)x ) θ (ξ ).

(3),(,)(u=xtθ=2 tt(k (θ ) θ ′)′ = −2c ρξθ ′, ξ >0,θ (0)=u1 , θ (∞)=u2 .(4)(5)Задача (4),(5) имеет единственное решение, которое в общем случаенаходится численно.Рассмотрим автомодельное решение уравнения (1) типа бегущейволны:ξ =−x Dt , u ( x, t ) =θ ( x − Dt ) =θ (ξ ).Из (1) и (6) получим: ( k (θ )θ ′)′ = − Dc ρθ ′, − ∞Ищем непрерывное решение, обладающее«тепловым потоком»:(6)<ξ <∞(7)непрерывным∂u ( x, t )−k (θ ) θ ′ ⇒ W ( x, t ) =−k (u ( x, t ))∂x(8)θ (ξ ) → 0, ξ → ∞; k (θ )θ ′ → 0, ξ → ∞.(9)Пусть бегущая волна начала движение по нулевому фонутемпературы:Тогда из формул (7) и (9) получим уравнение:k (θ ) dθ= − Dc ρθ dξ(10)) k0u , k0 > 0, σ > 0 (при σ = 2 получаем коэффициентПусть k (u=электронной теплопроводности в полностью ионизованнойплазме).

Тогда решение уравнения (10) имеет вид:σ51σ  σ Dcρ (−ξ )  , ξ ≤ 0,θ (ξ ) =   k0ξ > 0.0,(11)11u0t σ (1 − x ) σ ,Dtu ( x, t ) = 0,(12)Из формул (6) и (11) получим:0 ≤ x ≤ Dt ,x > Dt ,1гдеσ D cρ u0 =  .k02σФункция(12)–обобщенноерешениеквазилинейноговырождающегося параболическогоуравнения, обладающаяследующими свойствами:а) функция (12) финитна по х в любой конечный моментвремени: существует постоянная А>0 такая, что u(x,t)=0 при x ≥ А;б) функция (12) не имеет всюду непрерывных производных,входящих в уравнение (1): при σ = 1 они имеют разрыв первогорода, при σ = 2 - разрыв второго рода на прямой x = Dt.Замечание. Тепловой поток: k0u0σ +1 t σ1 (1 − x ) σ1 , σDDtW ( x, t ) = 0,0 ≤ x ≤ Dt ,x > Dt(13)является непрерывной функцией: W ( Dt − 0, t )= W ( Dt + 0, t )= 0;в) отличительной особенностью вырождающихся уравнений (1)является то, что они могут описывать процессы с конечнойскоростью распространения возмущения, то есть их решениямогут быть финитными функциями.Условие финитности по х решения типа «бегущей волны» имеетвид:1∫0k (η )ηdη < +∞.Данное неравенство является необходимым и достаточнымусловием конечности скорости распространения возмущения впроцессах, описываемых уравнением (1).Из данного условия следует существование конечного «фронта»тепловой волны, где может теряться необходимая гладкостьрешения, что требует определения решения в обобщенном смысле.Обобщенноерешениевырождающегосяуравнения(1)удовлетворяет не самому дифференциальному уравнению, анекоторымпостроеннымнаегоосновеинтегральнымсоотношениям, которые не содержат первых производных по х и по tи вторых производных по х, которые могут быть разрывными.Существенным моментом является то, что эти соотношенияспециальным образом учитывают непрерывность теплового потока(8).Решения вида (3) и (6) – единственные типы нетривиальныхавтомодельных решений, которые допускает уравнение (1) припроизвольных коэффициентах k(u).

Новые типы автомодельныхрешений появляются только при специальном виде функции k(u).Рассмотрим задачу: ∂u ∂σ ∂u=ku( ∂t ∂x 0 ∂x ), x > 0, t > 0,x, 0) 0, x ≥ 0,u (=u (0, t )= u t m , t ≥ 0, m > 0.0Автомодельноерешение: ξ=(14)(15)(16)mx=uxtut,(,)σ 1+ mσ10 θ (ξ )222k0 u0 t(17)Для функции θ (ξ )получаем задачу:(θ σ θ ′)′ + 12 (1 + mσ )ξθ ′ − mθ =0, ξ >0,(18)θ=(0) 1, θ (∞)=0,(19)непрерывность теплового потока−θ θ ′σ.Приm = σ −1решение (17) совпадает с бегущей волной (11).m ≠σ−1Призадача (18),(19) решается численно.Можно показать, что при любом m>0 существует такое ξ 0 (m , σ ) > 0 ,что θ (ξ ) = 0 при ξ > ξ 0 и θ (ξ ) > 0 при 0 ≤ ξ ≤ ξ 0. Решение(17) – тепловая волна, движущаяся по невозмущенному фонутемператур.

Фронт волны xΦ (t ) :σxΦ (t ) = ξ 0 k0 u0 t1221+ mσ2(20)ускоряет своё движение и приt → ∞ нагревает до бесконечнобольших температур всю полупрямую x ≥ 0.Решение асимптотически устойчиво при t → ∞относительноk (u )малых изменений начальных значений и малых отклоненийот степенной зависимости. При незначительных изменений условий(15),(16) основные закономерности процесса нагрева, которые даётпространственно-временная структура автомодельного решения(17), сохраняются.2. Режимы с обострением. Эффект локализациитеплаРежимом с обострением называется такой закон изменениянекоторой величины, который обеспечивает её неограниченноевозрастание в течение конечного времени.Рассмотримзадачутеплопроводности:Кошидляодномерногоut k0 (u 2u x ) x + q0u β , -∞ < x < ∞, 0 < t ≤ T ,=, 0) u0 ( x), -∞ <x <∞, k0 >0, q0 >0.u ( x=уравнения(1)(2)Конкуренциянелинейныхпроцессовтеплопередачиитепловыделения: эффект локализации процесса горения(проявлениепроцессасамоорганизациинелинейнойдиссипативной среды).

Возникновение в среде целого набораразличных структур, не взаимодействующих друг с другом.Сверхинтенсивное горение: β > 1.Введём преобразование переменных:t → qt ,0k0x→x q0и приведём уравнение (1) к безразмерному виду:2βut (u u x ) x + u .=(3)(4)Свойства решений уравнения (4) существенно различаются вслучаях β = 3, β < 3, β > 3.u = (u u ) .2Частный случай уравнения (4):tx xПолучим точное автомодельное решение уравнения (5) вида:2x1u ( x, t ) = 1 θ (ξ ), ξ =t4t(5)(6)Из формул (5),(6) получаем:2θ (ξ )θ ′(ξ =) + 14 θ (ξ ) ψ (ξ ), 2ξψ ′(ξ=) +ψ (ξ ) 0.2Решение Зельдовича – Компанейца – Баренблатта:1 1  2 x2   21=u A ( x, t ) 1  η 0 −   ,(7)24t  t +E0η0Где1zzutη0 2===,max0,,(0,),{ } A1π [ ]+t4 2u A ( x, t ) ≠ 0, x ∈ (−η0t ,η0t ), t > 0.1414Решение (7) описывает тепловые волны.

u A ( x, t ) - это решениетипамгновенноготочечногоисточникамощностиЕ0,действовавшего в точке х=0 в момент t=0.1) Пусть в среде имеетсясоответствующий β = 3 :источник=ut (u 2u x ) x + u 3 .тепловойэнергии,(8)Ищем автомодельное решение (8) в виде:u=1 θ ( x)T0 − tИз формул (8) и (9) следует уравнение для функции(9)θ ( x) :1θ ( x)θ ′′( x) + 2(θ ′( x)) 2 + θ 2 ( x) =2(10)Формулы (9) и (10) определяют локализованный режим: 3Lsπ x  2 c o s  L  , x < 2 , su A ( x, t ) = 1 T0 −t Ls≥x0,,2(11)где Т0 – время существования решения,Lsносителя решения в любой момент времени.= π 3 - длинаLs LsωL = (− , ).

ПриРешение локализовано в области2 23 → ∞, t → T .t → T0 , x ∈ ωL , u A ( x, t ) → ∞=, u A (0, t ) 102 (T0 −t )Тепловая структура (11) называется локализованным S-режимом собострением и представляет собой стоячую температурную волну.2) Рассмотрим уравнение (4) сβ = 2:=ut (u 2u x ) x + u 2(12)Ищем решение в виде:=u A ( x, t )1 θ (ξ ), ξ =x T − t ,0T0 −tгде Т0- время существования тепловой структуры.(13)Из формул (12) и (13) получаем уравнение:ξ ′′′(θ θ ) + θ − θ + θ 2 =0.22(14)Функция θ (ξ ) строго положительна на интервале ( −ξ 0 , ξ 0 )иравна нулю вне его.

Функция θ (ξ ) и ξ 0 определяются численно.Свойства решения:а) режим с обострением:=u A (0, t )1 θ (0) → ∞,T0 −tt → T0 ;б) в любой момент времени тепловая структура имеет конечныйправый x (t ) и левыйx (t ) фронты:−+x+ (t ) =ξ0T0 − t,x− (t ) = −ξ0T0 − t(15)в) фронты движутся с увеличивающейся скоростью и в моментобострения t=T0 тепловая структура охватывает всю прямую,нагревая её всюду до бесконечной температуры.Такой процесс горения, описываемый уравнением (12),называется HS –режимом.3) Рассмотрим уравнение (4) при β = 4 :ut (u 2u x ) x + u 4=(16)Мощность источника энерговыделения Q (u ) = u при большихтемпературах выше, чем в S-режиме ( β = 3) и HS-режиме ( β = 2).Такой тепловой процесс, описывающий сильно локализованныеструктуры, называется LS-режимом.Ищем решение уравнения (16) в виде:4u A ( x, t ) =1 θ (ξ ), ξ = x ,3 T −t6 T −t00где Т0 > 0 –время обострения решения.(17)Из формул (16) и (17) следует уравнение:− 1 θ + θ 4 0,(θ 2θ ′)′ − 1 ξθ ′=63θ (ξ )>0, -∞ <ξ <∞.

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

Тип файла PDF

PDF-формат наиболее широко используется для просмотра любого типа файлов на любом устройстве. В него можно сохранить документ, таблицы, презентацию, текст, чертежи, вычисления, графики и всё остальное, что можно показать на экране любого устройства. Именно его лучше всего использовать для печати.

Например, если Вам нужно распечатать чертёж из автокада, Вы сохраните чертёж на флешку, но будет ли автокад в пункте печати? А если будет, то нужная версия с нужными библиотеками? Именно для этого и нужен формат PDF - в нём точно будет показано верно вне зависимости от того, в какой программе создали PDF-файл и есть ли нужная программа для его просмотра.

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