Это не официальный сайт wikipedia.org 01.01.2023

Тензоры в физической кинетике — Википедия

Тензоры в физической кинетике

Масса, заряд, импульс и энергия в уравнениях механики сплошной средыПравить

Основные уравнения механики сплошной среды – непрерывности, движения и энергии – демонстрируют причины изменения во времени плотностей трех основных механических величин: массы m  , импульса p   и энергии ε  .

При этом:

  • первое слагаемое левой части каждого из названных уравнений представляет собой изменение плотности (количества в единице объема) соответствующей величины в единицу времени;
  • второе – результат обмена этой величиной выделенного единичного объема с соседними объемами;
  • третье – изменение плотности соответствующей величины в единицу времени под действием внешних сил;
  • правая часть – изменение плотности соответствующей величины в единицу времени в результате столкновений частиц в объеме.

При описании газа заряженных частиц одной из форм уравнения непрерывности является закон сохранения заряда в дифференциальной форме, в котором количество вещества представлено не массой m  , а зарядом q  .

Общими характеристиками массы, заряда, импульса и энергии являются:

Однако, по одному признаку импульс в этом списке стоит особняком. А именно: масса, заряд и энергия – скаляры. Соответственно плотности массы, заряда и энергии – скаляры, плотности потока массы, заряда и энергии – векторы.

Импульс же сам является вектором. Соответственно, плотность импульса есть вектор – полностью описывается тремя величинами. Плотность же потока импульса полностью описывается уже девятью величинами: любая из трех проекций импульса вместе с частицей может переноситься в любом из трех пространственных направлений. Таким образом, плотность потока импульса представляет собой тензор второго ранга (кинетический тензор):

Π = Π 11 Π 12 Π 13 Π 21 Π 22 Π 23 Π 31 Π 32 Π 33  , Π m n = m f ( v ) v m v n d 3 v  ,

где

  • Π m n   – количество m  -й проекции импульса, которое в единицу времени переносится через единицу поверхности в n  -м направлении.
  • f ( v )   – функция распределения частиц по скоростям.

Можно заметить, что половина следа (суммы диагональных компонент) кинетического тензора равна плотности кинетической энергии:

Π 11 + Π 22 + Π 33 2 = m 2 f ( v ) ( v 1 2 + v 2 2 + v 3 2 ) d 3 v = ε K ( V )  

В результате форма записи уравнения движения в традиционном представлении отличается от формы записи уравнений непрерывности:

  ρ   t + p ( V ) = δ   ρ δ   t  

и энергии:

  ε ( V )   t + q V F ( V ) = δ   ε ( V ) δ   t  .

А именно:

  p m ( V )   t + m n   Π m n   r n F m ( V ) = δ   p m ( V ) δ   t  

где:

  • ρ   – плотность массы;
  • p ( V ) = ρ V  – плотность потока массы, математически тождественная плотности импульса;
  • V   – среднемассовая скорость;
  • ε ( V ) = ε K ( V ) + ρ m ε χ   – плотность энергии;
  • q = m 2 f ( v ) v 2 v   d 3 v + f ( v ) ε χ v   d 3 v   – плотность потока энергии;
  • ε χ   – внутренняя энергия частицы;
  • F ( V )   – внешняя сила, действующая на единицу объема газа;
  • δ δ   t   – изменение в единицу времени в результате столкновений (здесь приведены уравнения для компоненты многокомпонентной среды).

Можно заметить три основных неудобства последней записи по сравнению с двумя предыдущими:

  • громоздкость;
  • необходимость записи в трех проекциях;
  • привязанность конкретно к декартовым координатам.

Последнее, например, означает, что в зависимости от системы координат, в которых решается задача, запись уравнения движения в проекциях будет иметь разные формы.

Такими же недостатками обладает и запись кинетического тензора в развернутой форме:

Π m n = ρ V m V n + δ m n P σ m n I  

и запись последнего слагаемого в приведенном выражении:

σ m n I = η (   V m   r n +   V n   r m 2 3 δ m n V )  

где  

Наличие названных недостатков скорее всего можно объяснить тем, что кинетика в физике и тензорный анализ в математике – это сравнительно молодые направления в науке, возникшие уже после того, как натурфилософия, фактически, уже разделилась на отдельные отрасли: физику, химию, математику и т.п. В отличие от Эйлера, Гаусса, Стокса физики уже были только физиками, а математики – только математиками.

В результате тензорный анализ в математике, с одной стороны, оказался достаточно отстраненным от проблем современной физики и, с другой стороны, не сформировал еще общепринятой и достаточно компактной символики.

Необходимость выбораПравить

Развитие математического аппарата любой естественной науки часто ставит исследователя перед выбором:

1. Оставаться в кругу уже определенных категорий, правил и символики, за счет громоздкости, необщности и большого количества выражений.

2. Обобщить понятия, упростить и сократить количество выражений за счет введения новых категорий, правил и символики.

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

Примером преимущества второго выбора в математике и физике является появление векторного анализа, возникшего ввиду трехмерности геометрического пространства.

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

Второй выбор – введение понятия вектора – означает необходимость усвоения немногих новых определений: вектор, скалярное произведение, векторное произведение и т.п., но легко окупается следующими выгодами:

- вектор сохраняет свою целостность в любой системе координат, в то время как значения проекций меняются;

- правила преобразования положения в скорость, скорости в ускорение как векторов, связь между скоростью и импульсом, характеристиками поля и силой как между векторами сохраняются в различных системах координат.

Наиболее существенную новизну сравнительно со скалярным анализом представляет здесь само понятие вектора – нужно просто привыкнуть к тому, что в геометрии и математике одна величина может характеризоваться не одним, а тремя числами – по числу пространственных измерений в нашей Вселенной. В операциях с кинетическим тензором мы сталкиваемся с названной выше необходимостью выбора – по-прежнему оперировать объектами двух типов (векторами и скалярами), описывая перенос импульса девятью скалярами или тремя векторами (со всеми издержками, названными выше) или ввести понятие и правила операций с новыми объектами, характеризующимися девятью числами. Массовость обращений к переносу импульса в механике сплошной среды скорее располагает ко второму выбору. Кроме того, есть соображения, по которым оптимальным здесь становится не просто определение нового класса объектов, но введение некоего "над-класса", к которому равно относятся и скаляры, и векторы и вновь вводимые объекты. Таким "над-классом" в математике и физике являются тензоры соответствующих рангов.

Понятие тензора определенного рангаПравить

В нашем случае тензоры являются математическим представлением конкретных физических величин, но не операторами в матричном анализе. Предлагаемая символика и правила относятся именно к такому случаю и не обязательно полностью соответствуют символике матричного анализа.

Тензором определенного ранга M   в I  -мерном пространстве называют величину, которая полностью описывается I M   числами – элементами тензора. Предметом механики сплошной среды является обычное трехмерное пространство ( I = 3  ), поэтому в дальнейшем мы и будем говорить только о нем. Таким образом, в нашем случае тензором ранга M   является величина, которая полностью описывается 3 M   элементами.

В таком случае:

  • скаляр, имеющий один элемент, есть тензор нулевого ранга;
  • вектор, имеющий три элемента, есть тензор первого ранга.

Появление нового класса объектов требует новой символики. А именно:

  • единственный элемент, из которого состоит скаляр Φ  , не требует индекса в записи значения;
  • каждый из трех элементов A m   вектора A   обозначается индексом m  , изменяющимся от 1 до 3 – соответственно числу измерений геометрического пространства;
  • каждый из 3 M   элементов A ( m 1 m 2 m 3 . . . m M 1 m M )    тензора M  -го ранга A   обозначается M   индексами m 1 m 2 m 3 . . . m M 1 m M  , изменяющимися от 1 до 3 – в дальнейшем для краткости вместо m 1 m 2 m 3 . . . m M 1 m M   будем писать m 1 . . . m M  .

В тензорном анализе, так же, как и в векторном, важным является понятие базиса, основанное на определении единичного тензора.

Единичным тензором M  -го ранга  есть тензор e ( m 1 . . . m M )  , в котором равны нулю все элементы, кроме равного единице m 1 . . . m M  -го элемента.

В таком случае:

  • e (     )  : единичный тензор нулевого ранга есть единица (единичный скаляр);
  • e ( m ) = i m  : единичный тензор первого ранга есть орта (единичный вектор).

Операции с тензорамиПравить

Для облегчения восприятия правила операций с тензорами покажем сравнительно с правилами аналогичных операций с векторами.

 Правило 1. Сложение тензоров и умножение тензора на скаляр

Вектор D   равен сумме векторов A   и B  , если элемент вектора D   равен сумме соответствующих элементов векторов   A   и B  :

1.1. D = A + B                         D m = A m + B m  .

Прямым следствием правила сложения векторов является правило умножения вектора на скаляр: вектор D   равен произведению вектора A   и скаляра C  , если элемент вектора D    равен произведению соответствующего элемента вектора A   и скаляра C  :

1.2. D = C   A                          D m = C   A m  .

Тензор M  -го ранга D   равен сумме тензоров такого же ранга A   и B  , если элемент тензора D   равен сумме соответствующих элементов тензоров A   и B  :

1.3. D = A + B                       D ( m 1 . . . m M ) = A ( m 1 . . . m M ) + B ( m 1 . . . m M )  .

Прямым следствием правила сложения тензоров является правило умножения тензора на скаляр: тензор D   равен произведению тензора A   и скаляра C  , если элемент тензора D   равен произведению соответствующего элемента тензора A   и скаляра C  :

1.4. D = C   A                       D ( m 1 . . . m M ) = C   A ( m 1 . . . m M )  .

Правило 2. Запись тензоров как суммы элементов

Вектор A   может быть представлен как векторная сумма элементов с использованием орт:

1.5.  A = m i m A m  .

При этом нет смысла говорить о результате произведения i m A m   – единственный смысл записи i m A m   состоит в указании, что величине A m   равен именно  m  -й элемент вектора .

Тензор M  -го ранга A   может быть представлен как тензорная сумма элементов с использованием единичных тензоров:

1.6.  A = m 1 . . . m M e ( m 1 . . . m M ) A ( m 1 . . . m M )  .

При этом нет смысла говорить о результате произведения e ( m 1 . . . m M ) A ( m 1 . . . m M )   – единственный смысл записи e ( m 1 . . . m M ) A ( m 1 . . . m M )   состоит в указании, что величине A ( m 1 . . . m M )    равен именно  m 1 . . . m M  -й элемент тензора A  .

Правило 3. Инвариантность произведения скаляра и единичного тензора

"Результат" произведения скаляра и орты не зависит от последовательности сомножителей:

1.7. i m A m = A m i m  .

"Результат" произведения скаляра и единичного тензора не зависит от последовательности сомножителей:

1.8.   e ( m 1 . . . m M ) A ( m 1 . . . m M ) = A ( m 1 . . . m M ) e ( m 1 . . . m M )  .

Правило 4. Тензорное произведение и представление единичных тензоров через орты

Внимание !!! Правило 4 является, фактически, единственным новым правилом тензорного анализа, не представленным в векторном анализе.

Тензорным произведением тензора M  -го ранга A   и тензора N  -го ранга B   является тензор M + N  -го ранга D  , если m 1 . . . m M n 1 . . . n N  -й элемент тензора D   равен произведению m 1 . . . m M  -го элемента тензора A   и n 1 . . . n N  -го элемента тензора B  :

1.9.  D = A   B                       D ( m 1 . . . m M n 1 . . . n N ) = A ( m 1 . . . m M ) B ( n 1 . . . n N )  .

Таким образом, тензорное скаляра и тензора произвольного ранга есть, фактически, "простое" произведение скаляра на тензор (1.4).

С учетом (9) единичный тензор M  -го ранга может быть представлен как кратное тензорное произведение орт:

1.10.  e ( m 1 . . . m M ) = i m 1 . . . i m M  .

Выражения (1.6) и (1.8) с использованием (1.10) можно записать так:

1.11,  A = m 1 . . . m M i m 1 . . . i m M A ( m 1 . . . m M )  ,

1.12.  i m 1 . . . i m M A ( m 1 . . . m M ) = i m 1 . . . i m i A ( m 1 . . . m M ) i m i + 1 . . . i m M  .

В дальнейшем вместо (1.6) можно использовать запись (1.11) как более удобную в случаях, которые будут названы ниже.

Правило 5. Произведение тензоров

Существуют три вида произведений тензоров: тензорное, векторное и скалярное. Каждому из произведений соответствует знак: пробел в тензорном, крестик в векторном и точка в скалярном. Кроме того, удобно использование обобщенного знака произведения "  ", соответствующего трем различным случаям:

1.13.    = (   "    ", " × ", " · " )  .

Правило 5 является прямым следствием (1.11) и (1.12) – скаляры в произведениях можно произвольно перемещать относительно знаков произведений – знак произведения тензоров фактически относится к ближайшим ортам:

1.14. A B = m n ( i m A m ) ( i n B n ) = m n ( i m i n ) A m B n = m n A m B n ( i m i n )  ,

1.15. A B = m 1 . . . m M n 1 . . . n N ( i m 1 . . . i m M A ( m 1 . . . m M ) ) ( i n 1 . . . i n N B ( n 1 . . . n N ) ) =  

= m 1 . . . m M n 1 . . . n N i m 1 . . . i m M 1 A ( m 1 . . . m M ) ( i m M i n 1 ) i n 2 . . . i n N B ( n 1 . . . n N ) =  

= m 1 . . . m M n 1 . . . n N A ( m 1 . . . m M ) B ( n 1 . . . n N ) i m 1 . . . i m M 1 ( i m M i n 1 ) i n 2 . . . i n N  .

Правило 6. Произведение орт

Тензорное, векторное и скалярное произведение орт имеют следующие значения:

1.16. i m i n = e ( m n )  ,

1.17. i m × i m = 0  ,

i 1 × i 2 = i 2 × i 1 = i 3  ,

i 1 × i 3 = i 3 × i 1 = i 2  ,

i 2 × i 3 = i 3 × i 2 = i 1  ,

1.18. i m i n = δ m n  ,

где   δ m n   – символ Кронекера:

1.19. δ m n = { 1 , m = n 0 , m n  .

Правило 7. Извлечение элемента из вектора и тензора

Произвольный вектор A    в различных выражениях может встречаться не в прямой записи (1.5), а как результат операций с другими векторами или скалярами. При необходимости "извлечения" конкретной проекции вектора A   из такой записи можно на основании (1.5) и (1.18) использовать операцию:

1.20. A m = i m A = A i m  

Аналогично, на основании (11) и (18) можно записать для элемента тензора ранга выше нулевого:

1.21. A ( m 1 . . . m M ) = i m M ( . . . ( i m 1 A ) ) = ( ( A i m M ) . . . ) i m 1  .

Дифференциальные операторы в применении к тензорамПравить

Правило 8. Операторы Гамильтона и Лапласа

Любой из трех знаков (1.13) может использоваться не только в произведениях, но и в обозначениях действия оператора Гамильтона  , имеющего, как известно, запись:

1.22. = n i n   r n  .

Результатом тензорного действия оператора   является градиент, векторного – ротор, скалярного – дивергенция. К сожалению,  в большинстве источников отсутствуют общие правила развернутой записи результатов действия оператора   в произвольной (не декартовой) системе координат. Приводятся общие записи в декартовой системе – единственной системе координат с неизменными ортами, а также частные случаи записи для наиболее часто используемых ортогональных криволинейный координат – сферических, цилиндрических.

Правильные результаты в использовании оператора можно, однако, получить, распространив Правило 5 на скалярные дифференциальные операторы – скалярные дифференциальные операторы в произведениях можно произвольно перемещать относительно знаков произведений:

1.23. = ( n i n   r n ) = n i n   r n  .

Таким образом, имеем:

  • градиент скаляра:

1.24. Φ = n i n   Φ   r n  ;

  • ротор вектора:

1.25. × A = n i n ×   A   r n = n m i n ×   r n ( i m A m ) = n m ( i n × i m   A m   r n + i n ×   i m   r n A m )  ;

  • дивергенция вектора:

1.26. A = n i n   A   r n = n m i n   r n ( i m A m ) = n m ( i n i m   A m   r n + i n   i m   r n A m )  .

Аналогично (1.24) – (1.26) правило (1.23) в применении к тензору дает:

  • градиент тензора произвольного ранга:

1.27. A = n i n   A   r n = n m 1 . . . m M i n   r n ( i m 1 . . . i m M A ( m 1 . . . m M ) ) =  

= n m 1 . . . m M ( i n i m 1 . . . i m M   A ( m 1 . . . m M )   r n + i n   r n ( i m 1 . . . i m M ) A ( m 1 . . . m M ) )  ;

  • ротор тензора ранга выше нулевого:

1.28. × A = n i n ×   A   r n = n m 1 . . . m M i n ×   r n ( i m 1 . . . i m M A ( m 1 . . . m M ) ) =  

= n m 1 . . . m M ( i n × i m 1 . . . i m M   A ( m 1 . . . m M )   r n + i n ×   r n ( i m 1 . . . i m M ) A ( m 1 . . . m M ) )  ;

  • дивергенция тензора ранга выше нулевого:

1.29. A = n i n   A   r n = n m 1 . . . m M i n   r n ( i m 1 . . . i m M A ( m 1 . . . m M ) ) =  

= n m 1 . . . m M ( i n i m 1 . . . i m M   A ( m 1 . . . m M )   r n + i n   r n ( i m 1 . . . i m M ) A ( m 1 . . . m M ) )  .

Обобщенно выражения (1.27) – (1.29) можно записать так:

1.30. A = n i n   A   r n = n m 1 . . . m M i n (   r n i m 1 . . . i m M A ( m 1 . . . m M ) ) =  

= n m 1 . . . m M ( i n i m 1 . . . i m M   A ( m 1 . . . m M )   r n + i n   r n ( i m 1 . . . i m M ) A ( m 1 . . . m M ) )  .

Правила дифференцирования производных орт в (1.27) – (1.30) аналогичны правилам дифференцирования произведений скаляров:

1.31.   r n ( i m 1 . . . i m M ) =   i m 1   r n ( i m 2 . . . i m M ) + i m 1   i m 2   r n ( i m 3 . . . i m M ) + i m 1 i m 2   i m 3   r n ( i m 4 . . . i m M ) + . . .  

При этом конкретная система координат представлена просто набором значений производных орт по проекциям координаты   i m   r n  .

В тензорном анализе, как и в векторном, используется также оператор Лапласа Δ  :

1.32. Δ = 2 =  .

Результат действия оператора Δ   на произвольный тензор в произвольной ортогональной системе координат можно получить с использованием тех же, что и выше, правил:

1.33. Δ A = m n i m   r m ( i n   A   r n ) = m ( 2 A   r m 2 + ( n i m   i n   r m )   A   r n )  .

Некоторые характеристики тензоров второго рангаПравить

Тензор A   второго ранга может быть представлен как матрица:

1.34. A = A ( 11 ) A ( 12 ) A ( 13 ) A ( 21 ) A ( 22 ) A ( 23 ) A ( 31 ) A ( 32 ) A ( 33 )  .

Следом T r A   тензора второго ранга A   называют сумму его диагональных элементов:

1.35. T r A = n A ( n n ) = A ( 11 ) + A ( 22 ) + A ( 33 )  .

  • Сопряженный тензор

Тензор A    называют сопряженным тензору A  , если элементы тензора A   получаются перестановкой индексов элементов тензора A  :

1.36. B = A                       B ( m n ) = A ( n m )  .

Можно заметить, что:

1.37. ( A ) A  .

Тензор A   называют симметричным, если его элементы, симметричные относительно главной диагонали, равны друг другу:

1.38. A = A                       A ( m n ) = A ( n m )  .

  • Унитарный тензор

Унитарный тензор δ   есть тензор второго ранга, недиагональные элементы которого равны нулю, а диагональные – единице:

1.39. δ ( m n ) = δ m n                        δ = 1 0 0 0 1 0 0 0 1 = m n i m i n δ m n = n i n i n  .

Можно показать следующие свойства унитарного тензора:

1.40. δ A = A  ,

1.41. ( δ   A ) = A  ,

1.42. ( δ A ) = A  ,

1.43. ( A   δ ) = δ   A  ,

1.44. ( δ A δ ) = ( A )  ,

где тензор, сопряженный градиенту вектора A  :

1.45. ( A ) = n   A   r n i n  

и внутреннее произведение  унитарного тензора и вектора A  :

1.46. δ A δ = n i n A   i n  .

Симметричные тензоры. Операция симметрииПравить

Тензор A   произвольного ранга является симметричным, если перестановка любой пары индексов не изменяет значение компоненты тензора, например:

- для симметричных тензоров 2-го ранга:

1.47. A ( m n ) = A ( n m )  ;

- для симметричных тензоров 3-го ранга:

1.48. A ( k m n ) = A ( k n m ) = A ( m k n ) = A ( m n k ) = A ( n k m ) = A ( n m k )  ;

- для симметричных тензоров 4-го ранга:

1.49. A ( i k m n ) = A ( i k n m ) = A ( i m k n ) = A ( i m n k ) = A ( i n k m ) = A ( i n m k ) =  

A ( i k m n ) = A ( k i n m ) = A ( m i k n ) = A ( m i n k ) = A ( n i k m ) = A ( n i m k ) =  

A ( i k m n ) = A ( k n i m ) = A ( m k i n ) = A ( m n i k ) = A ( n k i m ) = A ( n m i k ) =  

A ( k m n i ) = A ( k n m i ) = A ( m k n i ) = A ( m n k i ) = A ( n k m i ) = A ( n m k i )  

и так далее.

Можно заметить, что число независимых комбинаций индексов для тензора n  -го ранга равно n !  .

Любой тензор A [ n ]    произвольного ранга n   может быть преобразован в симметричный тензор  с помощью операции симметрии:

1.50. A [ n ] = 1 n ! A [ n ]  ,

где A [ n ]    – сумма исходного тензора A [ n ]    и всех тензоров, получаемых путем перестановки индексов его компонент аналогично (1.47) – (1.49).

Можно заметить, что если исходный тензор A [ n ]    уже является симметричным, имеет место A [ n ] = A [ n ]  .

Тензорная степень вектора. Бином и дифференциал тензорной степени вектораПравить

Тензор A [ n ]   произвольного ранга n   может быть результатом кратного тензорного произведения одного и того же вектора a  :

1.51. A [ n ] = a   a . . . a n  .

Для краткости можно использовать символ тензорной степени вектора (итерации вектора):

1.52. a   a . . . a n = a n  .

Показатель тензорной степени нужно записывать в скобках n  , чтобы не путать тензорный квадрат вектора:

1.53. a 2 = a   a  

с принятым в векторном анализе обозначением квадрата вектора (квадрата модуля):

1.54. a 2 = a a = | a | 2 = a 2  .

Можно убедиться, что:

1.55. a n a n  .

В алгебре скаляров существует запись для степени суммы скаляров (бином Ньютона):

1.56. ( a + b ) n = k = 0 n n ! k ! ( n k ) ! a n k b k  .

С использованием принятых здесь обозначений можно показать для тензорной степени суммы векторов:

1.57. ( a + b ) n = k = 0 n n ! k ! ( n k ) ! a n k b k  .

Для дифференциала степени скаляра известно, что:

1.58. d   a n = d ( a   a . . . a n ) = ( d   a ( a   a . . . a n 1 ) + a   d   a ( a   a . . . a n 2 ) + . . . n ) = n   a n 1 d   a  .

Можно показать, что для дифференциала тензорной степени вектора:

1.59. d   a n = d ( a   a . . . a n ) = ( d   a ( a   a . . . a n 1 ) + a   d   a ( a   a . . . a n 2 ) + . . . n ) = n a n 1 d   a  .

Кратное скалярное произведениеПравить

В операциях с тензорами часто встречается необходимость кратного скалярного произведения тензоров, для которого можно использовать следующую символику:

1.60. i m 1 . . . i m m ( k ) i n 1 . . . i n n = i m 1 . . . i m m k ( i m m k + 1 i n k ) . . . ( i m m i n 1 ) i n k + 1 . . . i n n  ,

то есть:

1.61. C [ m + n 2 k ] = A [ m ] ( k ) B [ n ] m 1 . . . m m n 1 . . . n n A ( m 1 . . . m m ) B ( n 1 . . . n n ) i m 1 . . . i m m ( k ) i n 1 . . . i n n  

и

1.62. C ( m 1 . . . m m + n 2 k ) = k 1 . . . k k A ( m 1 . . . m m k k 1 . . . k k ) B ( k 1 . . . k k m m k + 1 . . . m m + n 2 k )  .

Например, с использованием символа кратного скалярного произведения выражение (1.21) приобретает более компактную форму:

1.63. A ( m 1 . . . m m ) = i m m . . . i m 1 ( m ) A [ m ] = A [ m ] ( m ) i m m . . . i m 1  

Уравнение движения. кинетический тензор и тензор вязких напряженийПравить

 Использование предложенной символики и правил позволяет записать уравнение движения в универсальной и компактной форме без привязки к конкретной системе координат:

1.64.    p ( V )   t + Π = δ   p ( V ) δ   t  

Аналогично, для развернутой формы записи кинетического тензора имеем:

1.65. Π = ρ V   V + δ   P σ I  

и для тензора вязких напряжений:

1.66. σ I = 2   η ( V δ 3 V )  .

Следует отметить, что последнее выражение, так же как и выражение для кондуктивной составляющей плотности потока энергии (теплопроводность):

1.67. q κ = κ T  

представляют собой приближенные формы, применимые только при описании относительно плотных средств, когда производными величин σ I   и q κ   по времени и координатам можно пренебречь по сравнению с их изменением в единицу времени в единице объема в результате столкновений.

Одним из часто встречающихся недоразумений, связанных с "урезанностью" (1.66) и (1.67) является представление о том, что причинам вязкого переноса импульса и теплопроводности являются переменность среднемассовой скорости и температуры в пространстве. Аналогично, в диффузионном приближении, когда производными среднемассовой скорости по времени и координатам можно пренебречь по сравнению с ее изменением в единицу времени в единице объема в результате столкновений, причиной течения называют переменность давления в пространстве.

На самом деле, с учетом отброшенных в обоих названных случаях слагаемых течение, теплопроводность, вязкий перенос импульса, переменность скорости, температуры и давления являются следствиями общей в каждом из названных случаев причины. Например: изменение среднемассовой скорости и давления – как следствия переменности сечения канала в реактивных системах.

Уравнения моментов функции распределенияПравить

Механические моменты и моменты функции распределенияПравить

Основные уравнения газодинамики представляют собой уравнения моментов функции распределения частиц по скоростям. Функция f α ( v )  

2.1. d 2 N α = f α ( v ) d 3 r   d 3 v  ,

где

  • d 3 r   – элемент объема в пространстве координат;
  • d 3 v   – элемент объема в пространстве скоростей;
  • d 2 N α   – количество частиц в элементе объема d 3 r   в пространстве координат и элементе объема d 3 v   в пространстве скоростей.

Инструментом для отыскания функции f α ( v )    является кинетическое уравнение:

2.2.   f α ( v )   t + ( f α ( v ) v ) + v ( f α ( v ) F α ( v ) m α ) = δ   f α ( v ) δ   t  ,

где

  • F α ( v )  – сила, действующая на частицу сорта α  , имеющую скорость v  ;
  • v = i x   v x + i y   v y + i z   v z  – оператор Гамильтона в пространстве скоростей;
  • δ   f α ( v ) δ   t  – интеграл столкновений – изменение f α ( v )   в единицу времени в результате столкновений.

Основные механические характеристики частицы представляют собой моменты массы [1], где момент массы m α [ n ] ( v )   порядка n   определяется выражением:

2.3. m α [ n ] ( v ) = m α v n  .

Например:

  • момент массы 0-го порядка m α [ n ] ( v ) = m α v n   представляет собой просто массу частицы;
  • момент массы 1-го порядка m α [ 1 ] ( v ) = m α v 1 = m α v   представляет собой импульс частицы;
  • момент массы 2-го порядка m α [ 2 ] ( v ) = m α v 2 = m α v   v   представляет собой тензор 2-го ранга, не имеющий специального названия, но половина следа которого 1 2 T r ( m α [ 2 ] ( v ) ) = 1 2 m α v 2   есть кинетическая энергия частицы,

Основные газодинамические параметры представляют собой моменты функции распределения [1]:

2.4. M α [ n ] = n α m α [ n ] ( v ) = m α [ n ] ( v ) f α ( v ) d 3 v = m α n α v n  ,

где   

  • n α   – концентрация (количество частиц в единице объема);
  •       – символ осреднения по скоростям.

Например:

  • момент 0-го порядка M α [ 0 ] = m α n α v 0 = m α n α = ρ α   представляет собой плотность массы (массу единицы объема);
  • момент 1-го порядка M α [ 1 ] = m α n α v 1 = m α n α v = p α ( V )   представляет собой плотность импульса (количество импульса в единице объема), количественно равную плотности потока массы (количеству массы, в единицу времени переносимое через единицу поверхности);
  • момент 2-го порядка M α [ 2 ] = m α n α v 2 = m α n α v   v = Π α   представляет собой плотность потока импульса (кинетический тензор, количество импульса, в единицу времени переносимое через единицу поверхности), половина следа которого  1 2 δ ( 2 ) M α [ 2 ] = ε α V   есть плотность энергии (количество энергии в единице объема);
  • момент 3-го порядка M α [ 3 ] = m α n α v 3 = m α n α v   v   v = Q α   представляет собой тензор 3-го ранга, не имеющий специального названия, но половина вектор-следа которого 1 2 δ ( 2 ) M α [ 3 ] = q α   равна плотности потока энергии (количеству энергии, в единицу времени переносимое через единицу поверхности).

Похожее описание приведено в работе R. Fitzpatrick Plasma Physics : An Introduction, но для моментов, отнесенных к единице массы, и с записью только для следа-вектора момента третьего порядка Q α  .

Уравнения моментов функции распределенияПравить

Уравнение момента n  -го порядка функции распределения частиц по скоростям может быть получено умножением всех слагаемых кинетического уравнения (2.2) на момент массы n  -го порядка с последующим интегрированием всех слагаемых по всем значениям скорости.

В результате в применении к заряженной компоненте газа возникает уравнение следующего общего вида:

2.5.   M α [ n ]   t + M α [ n + 1 ] n q α m α M α [ n ] × B + M α [ n 1 ] E = δ   M α [ n ] δ   t  ,

где    δ   M α [ n ] δ   t     – изменение момента   в единицу времени в результате столкновений:

2.6. δ   M α [ n ] δ   t = m α [ n ] ( v ) δ   f α ( v ) δ   t d 3 v  .

В зависимости от порядка n   можно записать следующие случаи для уравнения (2.5):

  • при n  =0 – уравнение непрерывности:

2.7.   ρ α   t + p α ( V ) = δ   ρ α δ   t  ;

  • при n  =1 – уравнение движения:

2.8.   p α ( V )   t + Π α q α m α ( p α ( V ) × B + ρ α E ) = δ   p α ( V ) δ   t  ;

  • при n  =2 – уравнение потока импульса:

2.9.   Π α   t + Q α 2 q α m α Π α × B + p α ( V ) E = δ   Π α δ   t  ;

  • при n  =3 – уравнение моменте третьего порядка:

2.10.   Q α   t + M α [ 4 ] 3 q α m α Q α × B + Π α E = δ   Q α δ   t  .

Незамкнутость системы уравнений моментов функции распределения. Уравнения статических моментовПравить

На основе анализа уравнений (2.5) – (2.10) можно заметить, что система уравнений моментов функции распределения является принципиально незамкнутой – при записи уравнения для очередного неизвестного момента n  -го порядка, во втором слагаемом левой части возникает дивергенция момента порядка n + 1  .

В любом описании система уравнений газодинамики замыкается приближенно с использованием предположений того или иного уровня точности.

Оставляя пока открытым вопрос о незамкнутости системы уравнений моментов функции распределения, можно показать возможность ее иного представления с использованием записей статических моментов.

Сопутствующей системой координат называют инерциальную систему отсчета, в которой в данный момент в данной точке среднемассовая скорость компоненты равна нулю. Скорость частицы в сопутствующей системе  (хаотическая скорость) может быть представлена так:

2.11. v = v V α  .

При этом в соответствии с определением среднемассовой скорости V α = v   имеем:

2.12. v 0  .

Таким образом, первый момент  (плотность потока частиц, плотность импульса) в сопутствующей системе равен нулю по определению.

Моменты функции распределения в сопутствующей системе называются статическими моментами и  могут находиться подстановкой v   вместо v   в выражения для моментов функции распределения:

2.13. P α = m α n α v   v  ,

2.14. G α = m α n α v   v   v  ,

2.15. W α = m α n α v   v   v   v  ,

где 

  • P α   – тензор давления компоненты, равный кинетическому тензору Π α   в сопутствующей системе координат;
  • G α   – третий статический момент, равный третьему моменту Q α   в сопутствующей системе координат;
  • W α   – четвертый статический момент, равный четвертому моменту M α [ 4 ]   в сопутствующей системе координат.

Тензор G α   можно условно называть потоком давления.

Можно показать следующие связи между величинами полных и статических моментов:

2.16. Π α = ρ α V α V α + P α  ,

2.17. Q α = ρ α V α V α V α + 3 V α P α + G α  ,

2.18. M α [ 4 ] = ρ α V α V α V α V α + 6 V α V α P α + 4 V α G α + W α  .

При этом вместо уравнений потока импульса (2.9) и третьего момента (2.10) можно использовать уравнение давления:

2.19.   P α   t + ( V α P α + G α ) 2 q α m α P α × B + 2 P α V α = δ   P α δ   t  

и уравнение потока давления:

2.20.   G α   t + ( V α G α + W α ) 3 q α m α G α × B + 3 G α V α P α P α ρ α = δ   G α δ   t  ,

где  δ   P α δ   t    и δ   G α δ   t   – изменения P α   и G α   в единицу времени в результате столкновений, равные:

2.21. δ P α δ t = δ Π α δ t 2 V α δ p α ( V ) δ t + V α δ ρ α δ t  

2.22. δ   G α δ   t = δ   Q α δ   t 3 V α δ   Π α δ   t + P α V α p α ( V ) ρ α ( δ   p α ( V ) δ   t V α δ   ρ α δ   t )  .

Можно заметить, что половина следа тензора G α   представляет собой кондуктивную составляющую плотности потока энергии:

2.23. 1 2 δ ( 2 ) G α = q α κ  .

Следует также отметить, что понятие тензора вязких напряжений является реликтом, связанным с попыткой представления единичного объема газа как материального тела, изменение импульса которого происходит в результате действия неких сил.

На самом деле все три слагаемые в (1.65) соответствуют переносу импульса вместе с частицами, а не результата обмена импульсами (действия сил). Поэтому предпочтительным является представление кинетического тензора в виде (2.16) или в виде:

2.24. Π α = ρ α V α V α + δ   P α + π α  ,

где π α   – тензор вязкости, равный:

2.25. π α = P α δ   P α = σ I  .

Поток вектора в математике и физике принято считать положительным, если он направлен наружу из выделенного объема, а силу в физике - положительной, если она направлена внутрь. Этим и объясняется разница в знаках π α = σ I   между тензором вязкости (как составляющей плотности потока импульса вместе с молекулами и тензором вязких напряжений как "силой", действующей на объем.

Уравнение вида (2.19) приведено, например в книге Б. Росси и С. Ольберта "Introduction to the physics of space" [2], но в форме уравнения для компоненты тензора давления, а не в нашем компактном виде и без каких-либо рекомендаций о способах отыскания тензора G α  .

Приближение третьего рангаПравить

Как уже сказано, система уравнений моментов функции распределения принципиально незамкнута. Замыкание достигается приближенно в зависимости от степени анизотропии функции распределения в сопутствующей системе координат, которую считают необходимым учесть в конкретной решаемой задаче.

Уравнения (2.7), (2.8), (2.19) - (2.22) содержат тензоры от 0-го до 3-го ранга, необходимые для расчета характеристик устройства с учетом, в том числе и диссипативного переноса импульса и энергии. Проблему представляет тензор 4-го ранга W α  . При этом, моменты четных рангов не равны нулю даже при изотропном по v   распределении.

Например, при максвелловском распределении имеет место равенство:

2.26. W α = 3 δ δ P α P α ρ α  .

В работе [1] предложено следующее приближенное обобщение зависимости (2.26):

2.27. W α = 3 P α P α ρ α  .

В таком представлении уравнение (2.20) приобретает вид

2.28.   G α   t + ( V α G α ) 3 q α m α G α × B + 3 G α V α + 3 P α ( P α ρ α ) = δ   G α δ   t  ,

то есть не содержит уже новых неизвестных, что позволяет приближенно замкнуть систему на уровне моментов от 0-го до 3-го ранга.

Найдя след каждого слагаемого в (2.19) можно показать для скаляра давления P α = 1 3 δ ( 2 ) P α  :

2.29.   P α   t + ( V α P α + 2 3 G α ) + 2 3 P α V α + ( π α ) V α = δ   P α δ   t  .

С учетом (2.19), (2.29) для тензора вязкости можно записать:

2.30.   π α   t + ( V α π α + G α δ 3 q α κ ) 2 q α m α π α × B + 2 P α ( V α δ 3 V α ) = δ   π α δ   t  .

Найдя след каждого слагаемого в (2.28) с учетом (2.23) и подставляя P α = δ   P α   (как при максвелловском распределении) можно показать для теплопроводности:

2.31.   q α κ   t + ( V α q α κ ) + 3 2 q α κ V α 3 2 q α m α q α κ × B + 5 2 P α k ρ α T α = δ   q α κ δ   t  .

Для однородного газа правые части (2.30) и (2.31) могут быть представлены так[3]:

2.32. δ   π α δ   t = 3 2 π α τ α α d  ,

2.33. δ   q α κ δ   t = q α κ τ α α d  ,

где τ α α d  – эффективное время передачи давления[3].

Подстановка (2.32), (2.33) в (2.30) и (2.31) дает для однородного газа:

2.34.   π α   t + ( V α π α + G α δ 3 q α κ ) 2 q α m α π α × B + 2 P α ( V α δ 3 V α ) = 3 2 π α τ α α d  ,

2.35.   q α κ   t + ( V α q α κ ) + 3 2 q α κ V α 3 2 q α m α q α κ × B + 5 2 P α k ρ α T α = q α κ τ α α d  .

Граничные условия и характеристики столкновенийПравить

Можно заметить, что "классические" выражения для вязкости и теплопроводности можно получить из (2.34) и (2.35), оставляя в левых частях только последние слагаемые. Однако в разреженных газах, где роль столкновений в объеме мала, нельзя пренебрегать слагаемыми, связанными с изменением диссипативных характеристик в пространстве и во времени. В отличие от "классических" полные записи являются дифференциальными по искомым параметрам, а значит для их решения необходимы граничные условия.

Названные условия должны отражать факторы изотропии или анизотропии в процессах на границе газа (или плазмы).

Например, на границе плазмы с поверхностью или окружающей нейтральной средой существует пространственный слой заряда (ленгмюовский слой), по своей природе неоднородный и нестационарный. Отражение электронов от потенциального барьера в этом слое не является зеркальным - происходит релаксация импульса и трансформация моментов функции распределения высоких рангов.

Вопрос формулировки граничных условий с учетом названного процесса решался в работах[4][5][6].

Актуальным остается вопрос о записи правых частей для уравнений моментов порядка выше 1 с учетом упругих и неупругих столкновений в объеме.

См. такжеПравить

ЛитератураПравить

ПримечанияПравить

  1. 1 2 3 С. Ю. Нестеренко, Ш. Рошанпур. Система уравнений моментов функции распределения частиц по скоростям в разреженной среде индукционных источников плазмы, электронов и ионов // Авиационно-космическая техника и технология. — 2013. — № 7 (104). — С. 117-120. — ISSN 1727-7337. Архивировано 28 августа 2016 года.
  2. Б. Росси, С. Ольберт. Введение в физику космического пространства // Москва, Атомиздат. — 1974. — ISSN 5-9648-0006-8.
  3. 1 2 Е.М. Лифшиц, Л.П. Питаевский. Физическая кинетика. Теоретическая физика. Т. 10. // Москва: Наука. Архивировано 28 августа 2016 года.
  4. S h . R o s h a n p u r. Electron gas parameters change inside Langmuir layer in electric propulsion devices // Восточно-Европейский журнал передовых технологий. — 2013. — № 4/5 ( 64 ). — С. 36-39. — ISSN 1729-3774. Архивировано 16 сентября 2016 года.
  5. Le Quang Quyen, Ngo Dai Phong, S. Yu. Nesterenko, S. Roshanpour. Effect of electrons non-mirror reflection from potential shield on plasma borders inside helicon and Hall effect thrusters // IEPC-2013-411. Архивировано 28 августа 2016 года.
  6. А.В. Лоян, С.Ю. Нестеренко, Ш. Рошанпур, А.И.Цаглов. Математическое моделирование процессов в индукционных высокочастотных источниках плазмы и электронов // Авиационно-космическая техника и технология. — 2011. — № 10(87). — С. 203-206. — ISSN 1727-7337. Архивировано 17 сентября 2016 года.