Научная литература
booksshare.net -> Добавить материал -> Физика -> Рабинович М.И. -> "Введение в теорию колебаний и волн." -> 656

Введение в теорию колебаний и волн. - Рабинович М.И.

Рабинович М.И. Введение в теорию колебаний и волн. — НИЦ, 2000. — 564 c.
Скачать (прямая ссылка): vvedenievteoriuvoln2000.djvu
Предыдущая << 1 .. 650 651 652 653 654 655 < 656 > 657 658 659 660 661 662 .. 942 >> Следующая

Х3 = А - (В + 1) Аз + А3У3 + Dx (А2 - 2А3 + А4), (Р8-11)
Уз = ВАз - АзУз + D2 (У2 - 2У3 + У4), (Р8-12)
Х4 = А - (В + 1) А4 + X4Y4 + Dx (Аз - А4), (Р8-13)
У4 = ВХ4 - X4У4 + Д2 (Уз - У4). (Р8-14)
Эти дифференциальные уравнения можно представить в форме
(4.2.1), положив
Х = (АЬ У" А2, У2, Аз, Уз, А4, У4), р = (А, В, Dx, D2). (Р8-15)
108
Глава 4
Совершенно аналогично можно построить модели для систем реакторов,
соединенных различными способами, например, для линейной цепочки из пяти
и более реакторов, для замкнутого кольца из реакторов, для гексагональной
структуры из шести реакторов, расположенных по окружности с одним
центральным реактором, и т. д. Читатель может легко вывести
соответствующие уравнения, используя формулу (Р8-1).
4.2.9, Задача 9. Система типа "брюсселятор" с изменяющимся во времени
впуском
Схема модельной реакции "брюсселятор" может быть использована при
исследовании свойств нелинейных неавтономных систем, т. е. систем, для
которых время входит в правые части соответствующих дифференциальных
уравнений явным образом [4.37].
Будем предполагать, что концентрация исходного вещества А не постоянна, а
меняется во времени по следующему закону:
сА = А + a sin <оС (Р9-1)
Тогда изменение концентраций промежуточных продуктов X и Г в одном
реакторе описывается уравнениями, аналогичными (Р7-1) и (Р7-2):
^.==X2Y-(B+ 1)А + A + asinmf, (Р9-2)
~ = ВХ - X2Y, (Р9-3)
или
х, р), x = (X,Y), р = (А, В, а, <о). (Р9-4)
4.2.10. Задача 10. Модель Лоренца
Рассмотрим слой жидкости, изображенный на рис. 4.8. В направлении
вертикальной оси слой имеет толщину Я, в горизонтальных направлениях
размеры слоя считаются бесконечными. Будем предполагать, что жидкость
имеет свободную поверхность и что тепло, возникающее за счет внутреннего
трения в жидкости, пренебрежимо мало. Будем предполагать, далее, что
распределения скоростей и температуры, возникающие из-за разности
температур АТ, могут быть описаны системой дифференциальных уравнений с
частными производными в приближении Буссинеска (в частности, являются
постоянными плотность, коэффициенты вязкости и теплопроводности, см.
4.2. Задачи с сосредоточенными параметрами
109
[4.38, 4.40, 4.52]). Эту систему можно аппроксимировать системой
обыкновенных дифференциальных уравнений, получаемых в результате
разложения полей температур и скоростей течения в ряды Фурье [4.52].
Лоренц [4.40] применил очень простую аппроксимацию, в которой
используется всего три
Рис. 4.8. Слой жидкости, нагреваемый снизу (к модели Лоренца).
члена этих разложений. При этом он получил следующую систему уравнений:
х = оу - ах,
y = - xz + rx - y, (Р10-1)
z = xy - bz.
Здесь переменная х отвечает одной из компонент скорости, а переменные у,
г соответствуют членам разложения температуры в ряд Фурье. Параметр а
представляет собой число Прандтля, а параметр г - число Рэлея. В
векторной форме
(4.2.1) мы имеем
х = (х, у, z), р = (а, г, Ь).
Тривиальное стационарное решение системы (Р10-1), х - = у - z = 0,
соответствует случаю передачи тепла с помощью теплопроводности, когда
жидкость покоится и температурный профиль оказывается линейным. Два
других стационарных решения системы (Р10-1) имеют вид
x==y==±[b(r- 1)]1/2, z = г - 1. (Р10-2)
Эти решения соответствуют простому конвективному течению, изображенному
на рис. 4.8.
110 Глава 4
За последние десять лет модель Лоренца превратилась в наиболее изученную
модель возникновения хаотического движения жидкости [4.40]. При этом в
большинстве работ исследуется обычно изменение характера решения
(возникновение периодического решения определенного типа или же появление
хаотического решения) в зависимости от изменения параметра г (числа
Рэлея).
4.3. ЗАДАЧИ С РАСПРЕДЕЛЕННЫМИ ПАРАМЕТРАМИ
В этом параграфе будут сформулированы нелинейные задачи, приводящие к
решению систем дифференциальных уравнений с частными производными,
которые в дальнейшем, в гл. 6, используются для иллюстрации различного
рода численных подходов. Принимая во внимание сложности численного
анализа, мы будем рассматривать только системы с одной пространственной
координатой (параметрические исследования для систем с большим числом
пространственных переменных много труднее и в настоящее время только
начинают широко применяться на практике).
4.3.1. Системы типа "реакция - диффузия"
Рассмотрим s компонент, реагирующих между собой в ходе R независимых
реакций при постоянной температуре. Изменение концентрации во времени и
пространстве может быть описано системой локальных уравнений баланса
массы для s выбранных компонент
дс.
-jj- + div jj- = h. (4.3.1)
Здесь Ci - молярная концентрация компоненты i, ft описывает
возникновение г-й компоненты в результате R независимых
/ я
реакций (при этом fi-Yu^uri^ ГДе v<7 - стехиометрический V /=1
коэффициент t-й компоненты в /-й реакции и г,- - скорость
^ .
/-и реакции I, j; -молярная плотность потока компоненты и
Если учитывать только диффузионные и конвективные составляющие потока, то
Предыдущая << 1 .. 650 651 652 653 654 655 < 656 > 657 658 659 660 661 662 .. 942 >> Следующая

Реклама

c1c0fc952cf0704ad12d6af2ad3bf47e03017fed

Есть, чем поделиться? Отправьте
материал
нам
Авторские права © 2009 BooksShare.
Все права защищены.
Rambler's Top100

c1c0fc952cf0704ad12d6af2ad3bf47e03017fed