ТЕОРИЯ ИНТЕНСИВНЫХ ПРОЦЕССОВ
1. Получение и свойства неравновесных эффективных характеристик
Скачать:
Вложение | Размер |
---|---|
![]() | 563 КБ |
Предварительный просмотр:
Получение и свойства неравновесных эффективных характеристик
© ГЕ. Скворцов
С.-Петербургский государственный университет
Поступило в Редакцию 9 июля 1999 г
дается универсальный проекционный метод сокращения описания и демонстрируется его применение для всех уровней. Рассматриваются получающиеся в результате сокращения эффективные характеристики. Обсуждаются их свойства и поведение для широкого диапазона мер действия, а также обусловленные ими эффекты.
1. Изучение процессов возрастающей интенсивности неизбежно связано с учетом все более глубоких уровней структуры систем и увеличением числа вовлеченных в процесс величин. Аналогичное этому положение возникает для режимов структурных переходов, когда система “открывается” и реагирует на множество внешних факторов. В том и другом случае существенно увеличивается размерность “пространства действия”.
Вместе с тем любое наблюдение ограничено конечным числом величин. В результате возникает несоответствие между возрастанием размерности процесса и ограниченными возможностями его наблюдения и описания.
Естественным и, по существу, единственным способом решения этой проблемы является целесообразное сокращение описания (СО). При этом лишь несколько увеличивается число определяющих величин с общим уменьшением размерности описания и одновременно осуществляется переход к эффективным характеристикам (ЭХ).
Основным признаком последних служит зависимость от меры действия, т. е. это — существенно неравновесные характеристики. Таковыми являются эффективная масса, квазиэнергия, неравновесные восприимчивости. В общем случае, как известно, ЭХ зависят от мер действия операторным образом [1,2]. В конкретных задачах из разных областей (теория ядра, динамика молекулы, кинетика, гидродинамика) переход к ЭХ осуществляется обычно посредством использования феноменологических моделей. Таковые, как правило, не позволяют произвести интерполяцию для режимов существенной неравновесности, т. е. не дают возможности решить основную проблему.
Следует отметить, что в ряде случаев приходится иметь дело изначально с неравновесными характеристиками, поскольку их пределы при стремлении меры действия к нулю не существуют. Примеры такого рода указаны в [3,4].
Заметим, что в общем плане, согласно постулату всеобщей связи, все используемые характеристики неравновесных систем являются эффективными.
2. Представим универсальный метод СО и перехода к эффективным характеристикам.
Метод проекции исходных величин и уравнений на. подходящие подтгространства и получения замкнутых уравнений для интересующих компонент является универсальным и весьма удобным. Он с успехом используется для СО разного рода и уровня [5—7].
Приведем общую схему проекционного метода и выделим получающиеся при этом ЭХ.
Исходное описание динамики системы дается эволюционным уравнением для набора величины
(1)
набор общих координат, Q= {qs} — член источника. Это уравнение описывает систему на уровне более подробном, нежели нас интересует в итоге. Оно может быть неконструктивным, как например уравнение Лиувилля для 1020 частиц. Важно, что на его основе переходим посредством СО к следующему, менее подробному и более конструктивному уровню описания.
Сокращение описания любого уровня представляет собой проектирование на подпространство РF меньшей размерности
(2)
Это представление аналогично разложению в ряд Фурье: Fp - конечное число первых членов, - остаток.
Осуществляя разбиение (2) в уравнении (1) и производя проектирование на подпространства Fp и , получаем два уравнения
О —произвольный оператор (принято для сокращения выражений, что P u коммутируют). Первое уравнение (система) дает основу СО нужного уровня; второе — используется для замыкания. Подставляем операторное решение уравнения (4)
(5)
в уравнение (3) и получаем систему СО замкнутого вида
Величину принято называть “производством”, а выражения (7) — замыкающими соотношениями. Последние в гйдродинамяке определяют общие законы переноса и операторы переноса.
Соотношения (7) определяют ЭХ разного рода. С учетом вида резольвенты (5) они указывают качественную зависимость ЭХ от мер действия. В частности, в низшем порядке по малым мерам действия получаются равновесные характеристики (коэффициенты пропорциональности между “потоками” и “силами”). В противоположном предельном случае очень больших мер действия резольвента задает ЭХ в виде “убывающих операторов”. Например, при большой скоростной мере мер действия имеем ЭХ 1?
Актуализация правых частей уравнений (6), главным образом соотношений (7), для различных систем и уровней описания — ключевая проблема динамической теории.
Простейшим примером общей схемы СО служит сведение двух линейных уравнений
к одному (для x1y)
(8)
В данном случае эффективный оператор релаксации принимает вид Очевидно, при малой скорости - О он pавен , а при большой —.
Подобным образом Больцман и Максвелл, а затем Леонтович и Мандельштам получали из моментной системы релаксационную теорию.
3. Примёняя приведенную общую схему к динамическому описанию системы частиц с рождением и уничтожением, получаем общее кинетическое уравнение.
Такое уравнение описывает эволюцию набора функций распределения взаимодействующих комплексов структурно-кинетических элементов (s=1,2,… n), —. импульсы, координаты и внутренние характеристики s комплексов. Оно имеет вид
(9)
где V — скорость s-комплекса, — градиент, .[] — оператор сил, действующих на комплекс. Оператор сильных взаимодействий i[Е] учитывает всевозможные процессы рождения и уничтожения s-комплексов, а также запаздывание и нелокальность в области столкновений.
Простейшим вариантом уравнения (9) для разреженного газа бесструктурных частиц является уравнение Больцмана.
Уравнение (9) чрезвычайно сложное и конкретизация его для определенных систем непростая задача, а решение — открытая проблема.
Ввиду этого подходящим его предназначением является использование его как основы макроскопической теории, описывающей режимы широкого диапазона Неравновесности.
Такая теория для жидких систем (супергидродинамика) представлена в [5], а линейный ее вариант с достаточной полнотой рассмотрен в [4,8,9]. Основное место в этой теории занимают эффективные восприимчивости: операторы диффузии, термодиффузии, вязкости, вязкотеплопроводности и т. д., которые определяются при подстановке в соотношения (7) оператора в виде правой части (9).
4. Использование ЭХ (восприимчивостей) в полной макроскопической теории (супергидродинамике) выявляет многообразие эффектов, которые “не видны” в кинетической теории и отсутствуют в классических теориях.
Ввиду практической ценности укажем их основные классы.
Прежде всего это — перекрестные эффекты, которые в сувергидродинамике представлены полностью: тепловязкие, гидровязкотеплопроводные и т. д. [4]. Как правило, они малы при малой неравновесности, а при больших градиентах, Например, сравниваются с диагональными эффектами (вязкостью, теплопроводностью).
Хорошо известны релаксационные эффекты последействия, которые одним из первых наблюдал Больцман.
За последние десятилетия получили признание эффекты нелокальности, указываемые супергидродинамикой, для низкотемпературных условий, электромагнитоактивных сред, разреженных газов.
Эффекты сверхбольших восприимчивостей предписываются видом резольвенты (5) при условии малости обращаемого оператора (что соответствует неустойчивости структуры). В частности, эффект бесконечной проводимости для порогового значения электрического поля указан в [9].
Целый ряд эффектов, связанных с особенностями в зависимостях ЭХ (восприимчивостей) от мер действия, представлен в [10].
Таким образом, неравновесность ЭХ и ее адекватное отражение весьма существенны.
5. Наряду со свойствами указания эффектов СО играет важную вычислительную роль.
Известны вычислительные преимущества при сведении задач динамики систем с бесконечным числом степеней к конечной системе обыкновенных дифференциальных уравнений.
Как правило, это делается посредством общих рядов Фурье, и обычно при этом они обрываются, т. е. принимается= 0.
Это способ замыкания в общем случае недостаточен, и использование его при приближении даже очень большой размерности может оказаться малоэффективным. Так показывает пример задачи о звуке при частотах, превышающих характерную внутреннюю частоту [11]. В ней использование более пятисот моментов [12] не привело к согласию с опытом, в то время как выбор в качестве Е, первых пяти моментов (для одномерной задачи) в рамках проекционного метода дал полное совпадение.
Таким образом, следует рекомендовать всегда применять при сведении к конечномерной задаче полную проекционную схему, а затем использовать подходящие приближения при актуализации соотношений (7).
6. При рассмотрении динамики конечномерных систем осуществление СО в форме уменьшения числа уравнений или выделения регулярной компоненты процесса полезно в ряде отношений.
Продемонстрируем СО на примере системы трех уравнений с квадратичной нелинейностью
(10)
— постоянные коэффициенты скорости. Система (10) включает в себя большинство известных моделей (Вольтерра, Лоренца, Белоусова—Жаботянского и др.).
для СО — уменьшения числа уравнений, примем в соответствии с общей рекомендацией, что величина Х2 имеет наибольшее время жизни и наименьшую частоту по линейному приближению. Ввиду этого целесообразно сокращать величины Х1 и Х3.
Будем рассматривать, как это принято, нормальную реакцию, т. е. считать x1,3(t=0)=0. для уменьшения громоздкости итоговых выражений оставим в уравнениях по одному нелинейному недиагональному члену с v112, v2ЗI, VЗ2З и примем q1= 0.
С учетом этих упрощений итоговое уравнение для величины принимает вид
(11)
(12)
(13)
Полученное уравнение (11 )—( 13) существенно сложнее своего прототипа (10)2, однако как единый объект оно имеет ряд преимуществ перед системой (10). На его примере удобно продемонстрировать все известные особенности: автоколебания, временные бифуркации, странный аттрактор и переход к динамическому хаосу. Удобно оно и для применения приближенных методов, особенно с использованием усреднений и выделений вибрационных эффектов.
Этим вопросам будет посвящена специальная работа.
7. Завершая обсуждение последовательности сокращений, выделим основное звено проекционного метода и укажем его роль в различных задачах.
Основным звеном является учет характерных особенностей резольвенты СО (5) при актуализации замыкающих соотношений и ЭХ.
Характерные особенности резольвенты: асимптотическое поведение при больших степенях неравновесности, скачки и обращение в бесконечность на спектре, играют важную роль в фундаментальных и прикладных задачах.
В проблеме перехода от обратимого описания к необратимому (Боголюбов, Фок и Крылов, Пригожин) ключевая роль принадлежит особенности резольвенты, обусловленной непрерывным спектром. Учет непрерывного спектра в сочетании с СО усреднением обеспечивают переход к необратимости.
В вычислительных задачах, в которых используется метод общего ряда Фурье, как показывает пример задачи о звуке, целесообразно использовать резольвентную форму замыкания. Это позволяет учесть спектральные особенности более адекватно, нежели обычное “зануление” остатка.
Успешное применение ЭХ в широком диапазоне условий предполагает подходящее интерполирование на основе точных операторных выражений. даже весьма приближенное представление резольвенты позволяет получать более удачные выражения ЭХ, как например выражение эффективной проводимости турбулентной плазмы в магнитном поле [13].
Список литературы
[1] Скворцов ЛЕ. /1 Письма в ЖТФ. 1997. Т. 23. В. 22. С. 7.
[2] Скворцов ЛЕ. 1/ Письма в ЖТФ. 1998. Т. 24. В. 19. С. 7.
[3] Скворцов ЛЕ. /1 ЖЭТФ. 1969. Т. 57. В. 12. С. 2054. [4] Скворцов ТЕ. /1 ЖЭТФ. 1972. Т 63. В. 8. С. 502.
[5] Скворцов ЛЕ. 1/ Вестник ЛГУ 1979. В. 13. С. 94.
[6] Блехман ИИ Что может вибрация? М. 1988. 208 с.
[7] Скворцов ЛЕ. /1 ЖТФ. 1989. Т 59. ВЭ. С. 62.
[8] Скворцов ЛЕ. / ЖЭТФ. 1975. Т. 68. В. 3. С. 956.
[9] Перевозников Е.’-!., Скворцов ЛЕ. /1 ЖТФ. 1991. Т. 61. В. 9. С. 1
[10] Скворцов ТЕ. /1 Письма в ЖТФ. 1997. Т 23. В. 6. С. 85; В. 7. С. 23.
[11] Скворцов ЛЕ. 1/ ЖЭТФ. 1965. Т 49. В. 10. С. 1248.
[12] Pekeris C, Alterman Z, Finkeistein L., Frankowski K. II Phys. Fluids. 1962. V. 5. P. 1608.
[13] дмитриев АС., Синкевич О.А. /1 ТВТ 1997. Т 15. В. 3. С. 489.