On a class of vector fields

Full Text

Настоящая статья представляет собой расширенный вариант доклада [1], сделанного автором на Четвёртой международной конференции «Математическая физика и её приложения» (Россия, Самара, 25 августа - 1 сентября 2014). Введение. Аналогия между уравнениями гидродинамики Гельмгольца и уравнениями электродинамики Максвелла отмечалась в разных работах (см., например, [2-4]). Этот важный факт находит явное подтверждение в кинематической таблице физических величин - LT -таблице Брауна [5] и Бартини [6]. Эта таблица использовалась рядом авторов (П. Г. Кузнецов, О. Л. Кузнецов, Б. Е. Большаков, А. Е. Петров) при анализе и поиске общих законов природы. Аналогия разных физических величин, попадающих в одну клетку Lm T n кинематической таблицы, где L и T есть соответственно размерности длины и времени, заключается в том, что они имеют одинаковую качественную определённость. При заполнении этой таблицы можно исходить из разных физических формул. Однако во всех LT -таблицах для размерности электрического заряда и массы принята формула Максвелла [q] = [m] = L3 T -2 . Кроме того, общей оказывается и формула для размерности силы [F ] = L4 T -4 . Из выражения для лоренцевой силы F = q(E +v ×B), где E и B есть соответственно напряжённость и индукция электромагнитного поля, а v - скорость движения заряда ([v] = LT -1 ), находим выражения для размерностей [E] = = LT -2 , [B] = T -1 . Длина и время - это две физические величины, для измерения которых созданы наиболее совершенные приборы: лазерный радар и атомные часы. Поэтому можно предположить, что значение кинематической таблицы будет возрастать в процессе дальнейшего научного исследования природы и общества. 1. Поле смещений вакуума. Основу математического исследования составляет трёхпараметрическое поле (смещений вакуума): u(x; t) = P (x) cos k(x)t + Q(x) sin k(x)t. Здесь t - время; k(x) - частота колебаний в точке x трёхмерного евклидова пространства; P (x), Q(x) - ортонормированная пара стационарных векторных полей; (k, P, Q) - список параметров смещения. Оно получено из простого постулата: поле смещений вакуума есть нормированное электрическое поле. При этом нормировочный коэффициент k 2 (x) имеет размерность T -2 . Он обеспечивает единичную норму смещения u(x; t) при любом t. Вычислим скорость поля смещений v(x; t) = ∂u(x; t) = k(x)(Q(x) cos k(x)t - P (x) sin k(x)t). ∂t Напряжённость электрического поля, отвечающего указанному распределению поля смещения вакуума, даётся формулой E(x; t) = - ∂v(x; t) ∂ 2 u(x; t) =- = k 2 (x)u(x; t). ∂t ∂t2 При этом магнитная индукция B(x; t) = rot v(x; t) = rot(k(x)(Q(x) cos k(x)t - P (x) sin k(x)t)). Динамика электромагнитного поля, отвечающего паре (E(x; t), B(x; t)), описывается известными уравнениями Максвелла (1861). Математической 681 И с л а м о в Г. Г. теории этих уравнений посвящено огромное число работ. Мы упомянем лишь одну работу [7]. Исторический аспект возникновения уравнений Максвелла рассматривался разными авторами (см., например, [8-10]). Известные виды симметрии уравнений Максвелла подробно рассмотрены в [11]. Независимость и полнота уравнений Максвелла, а также теоремы единственности в теории электромагнетизма, изучались в работе [12]. Непосредственным применением аппарата векторного исчисления [13] проверяется справедливость следующего утверждения. Теорема 1. Пусть H есть векторная проекция индукции B(x; t) на векторное произведение P (x)×Q(x). В каждой точке x евклидова пространства тройка векторов (E, H, v) ортогональная. Вектор E вращается в плоскости векторов P (x) и Q(x) от первого вектора ко второму с частотой k(x). Вектор v вращается в той же плоскости от Q(x) к -P (x) с той же частотой. Вектор H перпендикулярен плоскости вращения векторов E, v. Следующее утверждение представляет собой переформулировку теоремы разложения Гельмгольца-Ходжа для векторных полей [14]. Теорема 2. Поле смещений u(x; t) в любой ограниченной области D с гладкой границей ∂D может быть единственным образом разложено в сумму скалярной и вихревой компоненты u(x; t) = ψ(x; t) + w(x; t), причём div w = 0, ∆ψ = div u, на границе области D вихревое поле w перпендикулярно к нормали границы. Для решения трёхмерной задачи Неймана в предыдущей теореме существует множество методов (см., например, [15]). В приводимой ниже теореме раскрывается физический смысл скалярного и векторного потенциалов электромагнитного поля [16] , а также скалярного магнитного поля, рассмотренного в работе [17]. Теорема 3. Скалярный потенциал ϕ электромагнитного поля (E, B) есть 2 ускорение ∂ ψ . ∂t2 Векторный потенциал A поля (E, B) есть скорость изменения ∂w вих∂t ревой компоненты поля смещений u(x; t). 3 Скалярное магнитное поле B ∗ есть рывок - c1 ∂ ψ , c - скорость света в 2 ∂t3 вакууме. Применим трёхпараметрическое поле смещений вакуума к отысканию локальных решений уравнений Максвелла. Предположим сначала, что круговая частота не зависит от выбора точки x : k(x) ≡ ω. Ниже будет показано, что это предположение выполнено для уравнений Максвелла без источников: div E = 0, div B = 0, rot E = - ∂B , ∂t rot B = 1 ∂E , c2 ∂t где c - скорость света. Подстановка выражений для E(x; t) и B(x; t) даёт div P = 0, 682 div Q = 0, rot rot P = ω c 2 P, rot rot Q = ω c 2 Q. Об одном классе векторных полей Типичная схема решения спектральной задачи для квадрата ротора основана на векторном тождестве rot rot F = grad div F - ∆F, содержащем векторный оператор Лапласа ∆. В евклидовой системе координат он распадается на тройку скалярных операторов Лапласа, что не имеет места в криволинейных системах координат [18]. Внимательный анализ выписанных выше уравнений показывает, что достаточно решить спектральную задачу для ротора rot F = λF , где в нашем случае λ = ω/c. Это известная задача, которой посвящён ряд работ (см. цитирования ниже и приведённую там литературу). Как правило, эту задачу рассматривают в классическом гильбертовом пространстве квадратично-суммируемых в некоторой области Ω гладких вектор-функций с определёнными граничными условиями (см., например, [19, 20]). В приложениях обычно область Ω ограничена (случай локальных решений уравнений Максвелла). При этом требуют, чтобы векторное поле F (x) было «приглажено» на границе области Ω (см. [21, 22]). В частности, для шара радиуса R с центром в начале координат ненулевое собственное значение спектральной задачи для ротора с условием перпендикулярности поля к нормали границы шара даётся формулой λ= ρm,n , R где ρm,n - любой нуль функции Бесселя дробного порядка Jn+1/2 (z), m, n 1 (случай n = 0 функции Бесселя, как будет показано ниже, приводит к решениям, сингулярным на полюсах сферы). При этом кратность собственного значения ρm,n равна 2n + 1. Доказательство этого факта приведено в работе Сакса [23]. Там же найдена соответствующая система собственных функций. При этом отмечено, что ρm,n совпадают с нулями функции ψn (z) = (-z)n d zdz n sin z z и, таким образом, лежат на действительной оси и симметричны относительно точки z = 0. Отсюда получаем важную для нас формулу, связывающую круговую частоту колебаний локального электромагнитного поля с радиусом охватывающего его шара: c ω = ρm,n . R При n = 1 имеем sin z - z cos z ψ1 (z) = . z2 Нас интересует минимальный положительный нуль этой функции ρ1,1 . Вычисления в системе Mathematica показывают, что ρ1,1 ≈ 4.49341. Применим полученную формулу для частоты электромагнитного поля к шару радиуса протона. Из физического справочника берём радиус протона 683 И с л а м о в Г. Г. R = 8.79 · 10-16 м. Скорость света в вакууме c = 299792458 м/с. Отсюда получаем круговую частоту ω = 1.53253 · 1024 . Произведём сравнение энергий, вычисленных по двум разным формулам: де Бройля ( ω) и Эйнштейна (mc2 ). Произведение приведённой постоянной Планка = h = 1.054571628 · 10-34 2π на найденную частоту ω даёт величину ω = 1.61616 · 10-10 . Масса покоя протона m = 1.672621637·10-27 . Вычисляя энергию по второй формуле, получим mc2 = 1.50328·10-10 . Относительная погрешность расхождения результатов вычисления по двум известным формулам для вычисления энергии (де Бройля и Эйнштейна) равна 0.0750899, т.е. меньше восьми процентов. Можно явно выписать для нуля ρ1,1 ≈ 4.49341 все три собственные функции, образующие базис подпространства собственных функций спектральной задачи для ротора, и тем самым получить систему трёх электрических полей, образующих «остов» предположительно протона. Заметим, что эти поля колеблются с одной и той же круговой частотой ω = 1.53253 · 1024 . Но мы рассмотрим наиболее интересный случай, когда n = 0. Он не рассмотрен в работе [23] по понятной причине: не выполнены наложенные автором граничные условия. Хотя кратность нулей функции ψ0 (z) = sin z/z равна единице, тем не менее соответствующая спектральная задача для ротора имеет два линейно независимых решения, правда, сингулярных. Если ограничиться только положительными нулями, то наименьшее значение равно π. Поэтому нас будет интересовать следующая спектральная задача: π rot F = F R в шаре радиуса R. Вычислим частоту ω = πc/R электрического поля в шаре R радиуса электрона. Относительно значения последней величины есть множество мнений. Мы же рассчитаем значение R исходя из равенства ω = mc2 , которое мы уже использовали выше. Но в этой формуле теперь возьмём массу покоя электрона m = 9.10938215 · 10-31 кг. Имеем ω = mc2 / . Отсюда находим радиус шара R = π /(mc) = 1.21316× × 10-12 м. Видим, что найденное значение на три порядка больше классического радиуса электрона, равного 2.8179402894 · 10-15 м. При этом частота колебаний ω = πc/R = 7.76344 · 1020 , т.е. на четыре порядка ниже частоты колебаний электрического поля для шара радиуса протона. 2. Задача для ротора в криволинейных координатах. Так как трёхмерная область Ω, отличная от шара, может иметь свою симметрию, рассмотрим спектральную задачу rot F = λF в общей криволинейной системе координат. Если в области Ω введены криволинейные координаты (ξ1 , ξ2 , ξ3 ), то радиус-вектор A текущей точки (x1 , x2 , x3 ) евклидовой системы координат является функцией от криволинейных координат (ξ1 , ξ2 , ξ3 ). Через произвольно 684 Об одном классе векторных полей взятую точку криволинейной системы координат проходят три координатные поверхности и три координатные линии. В общем случае координатные поверхности и линии будут криволинейными [13]. ∂A ∂A ∂A Предполагается, что система векторов ∂ξ1 , ∂ξ2 , ∂ξ3 линейно независима, за исключением конечного числа особых точек области Ω. В этом случае система называется подвижным репером криволинейной системы координат. Смешанное произведение векторов подвижного репера отлично от нуля. Отсюда следует линейная независимость попарных векторных произведений векторов этой системы. Поэтому подвижному реперу отвечает сопутствующий базис, составленный из попарных векторных произведений векторов этого репера. По каждому такому базису криволинейной системы координат можно разложить произвольно взятый вектор. Иногда бывает удобно разложить скалярное поле или векторное поле не по векторам подвижного репера, а по векторам сопутствующего базиса. В [13] для градиента скалярного поля Φ(ξ1 , ξ2 , ξ3 ) в терминах скалярного a · b, векторного a × b и смешанного (a, b, c) произведений векторов даётся вывод следующих разложений: grad Φ= 1 ∂A ∂A ∂A ∂ξ1 , ∂ξ2 , ∂ξ3 ∂Φ ∂A ∂A ∂Φ ∂A ∂A ∂Φ ∂A ∂A · × + · × + · × . ∂ξ1 ∂ξ2 ∂ξ3 ∂ξ2 ∂ξ3 ∂ξ1 ∂ξ3 ∂ξ1 ∂ξ2 Выражение для дивергенции векторного поля F (ξ1 , ξ2 , ξ3 ) в криволинейных координатах получено в следующей удобной форме:  div F = 1 ∂A ∂A ∂A ∂ξ1 , ∂ξ2 , ∂ξ3 ∂A ∂A ∂ F, ∂ξ2 , ∂ξ3  ∂ξ1 + ∂A ∂A ∂ F, ∂ξ3 , ∂ξ1 + ∂ξ2 ∂A ∂A ∂ F, ∂ξ1 , ∂ξ2 ∂ξ3  . Для ротации этого поля получена формула в виде следующего символического определителя: rot F = 1 ∂A ∂A ∂A ∂ξ1 , ∂ξ2 , ∂ξ3 · ∂A ∂ξ1 ∂ ∂ξ1 ∂A F · ∂ξ1 ∂A ∂ξ2 ∂ ∂ξ2 ∂A F · ∂ξ2 ∂A ∂ξ3 ∂ ∂ξ3 ∂A F · ∂ξ3 . ∂A ∂A ∂A Если векторы подвижного репера ∂ξ1 , ∂ξ2 , ∂ξ3 образуют ортогональную систему, то обычно её нормируют. Получаемый таким способом новый репер e0 , e0 , e0 будет ортонормированным. Предположим, что векторы этого ре1 2 3 пера образуют правую систему. Тогда будут иметь место следующие соотношения: e0 × e0 = e0 , 2 3 1 e0 × e0 = e0 , 3 1 2 Обозначив h1 = ∂A , ∂ξ1 h2 = e0 × e0 = e0 , 1 2 3 (e0 , e0 , e0 ) = 1. 1 2 3 ∂A , ∂ξ2 ∂A ∂ξ3 h3 = 685 И с л а м о в Г. Г. (эти скалярные величины называются коэффициентами Ламе ортогональной криволинейной системы координат), получим следующее представление подвижного репера через ортонормированный: ∂A = h1 e0 , 1 ∂ξ1 ∂A = h2 e0 , 2 ∂ξ2 ∂A = h3 e 0 . 3 ∂ξ3 Дифференциальные операции в ортогональной криволинейной системе координат в терминах ортонормированного репера e0 , e0 , e0 и коэффициентов 1 2 3 Ламе {h1 , h2 , h3 } принимают компактный вид: 1 ∂Φ 0 1 ∂Φ 0 1 ∂Φ 0 e1 + e2 + e , h1 ∂ξ1 h2 ∂ξ2 h3 ∂ξ3 3 ∂(h2 h3 F · e0 ) ∂(h3 h1 F · e0 ) ∂(h1 h2 F · e0 ) 1 1 2 3 · F = div F = + + , h1 h2 h3 ∂ξ1 ∂ξ2 ∂ξ3 h2 e 0 h3 e0 h1 e0 1 2 3 1 ∂ ∂ ∂ × F = rot F = · . ∂ξ1 ∂ξ2 ∂ξ3 h1 h2 h3 0 h F · e0 h F · e0 h1 F · e1 2 3 2 3 Φ = grad Φ = Решение F спектральной задачи λF = rot F разложим по ортонормированному базису: F = F1 e0 + F2 e0 + F3 e0 . 1 2 3 Разлагая символический определитель в выражении для rot F по элементам первой строки, получим ∂G3 ∂G2 0 1 ∂G3 ∂G1 0 1 ∂G2 ∂G1 0 1 e1 - e2 + e , - - - h1 g1 ∂ξ2 ∂ξ3 h2 g2 ∂ξ1 ∂ξ3 h3 g3 ∂ξ1 ∂ξ2 3 где приняты следующие замены переменных: G1 = F1 h1 , G2 = F2 h2 , G3 = F3 h3 и следующие преобразования коэффициентов Ламе: g1 = h2 h3 , h1 g2 = h1 h3 , h2 g3 = h1 h2 . h3 Нетрудно видеть, что h1 h2 h3 = g1 g2 g3 , а обратное преобразование имеет вид √ √ √ h1 = g2 g3 , h2 = g1 g3 , h3 = g1 g2 . Теперь видно, что скалярная форма спектральной задачи для ротора преобретает простой вид относительно новой группы переменных G1 , G2 , G3 : λG1 g1 = ∂G3 ∂G2 - , ∂ξ2 ∂ξ3 λG2 g2 = ∂G1 ∂G3 - , ∂ξ3 ∂ξ1 λG3 g3 = ∂G2 ∂G1 - . ∂ξ1 ∂ξ2 В силу симметрии шара радиуса R, где решается спектральная задача для ротора, целесообразно перейти в сферическую систему координат относительно переменных ξ1 = r, 0 686 r < ∞; ξ2 = θ, 0 θ π; ξ3 = ϕ, 0 ϕ < 2π, Об одном классе векторных полей которые называются соответственно «радиусом», «широтой» и «долготой». Формулы перехода от сферической системы к евклидовой задаются соотношениями x1 = r sin θ cos ϕ, x2 = r sin θ sin ϕ, x3 = r cos θ. Поэтому для радиус вектора A текущей точки (x1 , x2 , x3 ) евклидовой системы координат имеем A = {r sin θ cos ϕ, r sin θ sin ϕ, r cos θ} . Отсюда находим, что векторы подвижного репера выражаются формулами ∂A = {sin θ cos ϕ, sin θ sin ϕ, cos θ} , ∂r ∂A = {r cos θ cos ϕ, r cos θ sin ϕ, -r sin θ} , ∂θ ∂A = {-r sin θ sin ϕ, r sin θ cos ϕ,0} . ∂ϕ Параметры Ламе следующие: h1 = 1, h2 = r, h3 = r sin θ. Их преобразования имеют вид g1 = r2 sin θ, g2 = sin θ, g3 = 1 . sin θ Поэтому в сферической системе координат спектральная задача для ротора λF = rot F в новых переменных G1 = F1 , G2 = F2 r, G3 = F3 r sin θ примет следующий вид: λG1 r2 sin θ = ∂G3 ∂G2 - , ∂θ ∂ϕ λG2 sin θ = ∂G1 ∂G3 - , ∂ϕ ∂r 1 ∂G2 ∂G1 = - . sin θ ∂r ∂θ Непосредственной подстановкой проверяется, что эта система имеет два решения. Первое решение: λG3 G1 = 1 sin(λr), r2 G2 = λ cos(λr) cot θ, G3 = -λ sin(λr) cos θ. Второе решение: G1 = 0, G2 = - cos(λr) 1 , sin θ G3 = sin(λr). Заметим, что найденные решения зависят только от радиуса и широты θ, но не зависят от долготы ϕ. 687 И с л а м о в Г. Г. Переходя к основным переменным и вспоминая, что λ = π решения спектральной задачи rotF = R F : F1 = π R, получим два πr 1 π πr π πr sin , cos cot θ, - sin cot θ , 2 r R rR R rR R F2 = 1 πr 1 πr 0, - cos csc θ, sin csc θ . r R r R Эти сингулярные решения указаны и в статье [24], при этом автор отмечает, что он находил их в работах 1966 года. Возможно, они приводятся в более ранних публикациях по электромагнетизму. Для устранения сингулярности введём скалярное произведение с весом (r sin θ)2 : R π (F 1 , F 2 ) = F 1 · F 2 (r sin θ)2 d r d θ. 0 0 Найденные решения ортогональны относительно этого скалярного произведения и после нормировки дают нам стационарные поля P = F 1/ (F 1 , F 1 ), Q = F 2 / (F 2 , F 2 ) изменяющегося с частотой ω = πc/R = 7.76344 · 1020 локального электрического поля (предположительно) электрона E(x; t) = ω 2 (P (x) cos ωt + Q(x) sin ωt). 3. О независимости частоты от положения точки в пространстве. Изучим характер зависимости частоты k(x) электрического поля E(x; t) от положения x в пространстве R3 . Известная формула rot f (x)F (x) = f (x) rot F (x) + f (x) × F (x) (здесь f (x) - скалярное, а F (x) - векторное поле в трёхмерном пространстве), применённая к равенству B(x; t) = rot[k(x)(Q(x) cos k(x)t - P (x) sin k(x)t)], даёт линейный рост зависимости индукции магнитного поля от времени t. При вычислении дивергенции электрического поля E(x; t) = k 2 (x)(P (x) cos k(x)t + Q(x) sin k(x)t), воспользовавшись формулой div f (x)F (x) = f (x) div F (x) + f (x) · F (x) (здесь точка “·” обозначает скалярное произведение в R3 ), опять получим линейный рост зависимости от времени t. 688 Об одном классе векторных полей Отсюда следует, что первое уравнение системы Максвелла ∂B 1 ∂E , rot B = 2 ∂t c ∂t выполняется тогда и только тогда, когда имеют место следующие тождества: div E = 0, div P (x) = 0, div B = 0, rot E = - k(x) · P (x) = 0, div Q(x) = 0, k(x) · Q(x) = 0. Второе и третье уравнение автоматически выполнены в силу представления трёхпараметрического поля смещения вакуума. Осталось получить эквивалент последнего уравнения. Имеем rot B(x; t) = rot rot v(x; t) = k(x) c 2 v(x; t). В силу самосопряжённости оператора ротора у него нет присоединённых векторов. Поэтому переходим к спектральной задаче для ротора rot v(x; t) = k(x) v(x; t). c Так как v(x; t) = k(x)(Q(x) cos k(x)t - P (x) sin k(x)t), получим эквивалентную последнему уравнению Максвелла систему тождеств rot P (x) = k(x) P (x), c k(x) × P (x) = 0, k(x) Q(x), k(x) × Q(x) = 0. c Все условия на функцию частоты k(x) приводят к простому равенству rot Q(x) = k(x) = 0, что в случае отсутствия источников электромагнитного поля означает независимость частоты колебаний электрического поля от рассматриваемой точки x. Однако, если применить трёхпараметрическую форму с переменной частотой k(x) к системе уравнений Максвелла с источниками, то обнаружим линейный по времени рост плотности электрического заряда и плотности тока. 4. Малые течения физического вакуума. Течение физического вакуума описывается уравнением с нестационарным полем скоростей: dy = v(y; t), dt y(0) = x. Характер течения зависит от величины времени наблюдения за процессом. Предполагая, что интервал наблюдения мал, можно надеяться, что приращение h = y(t) - x будет мало. В окрестности точки (x; 0) имеет место разложение ∂v(x; 0) v(y; t) = v(x; 0) + v(x; 0)h + t + O(h2 + t2 ). ∂t 689 И с л а м о в Г. Г. Введём обозначения для тензоров: A= v(x; 0), 1 D = (A + AT ), 2 1 S = (A - AT ) 2 и векторов: b = B(x; 0) = rot v(x; 0), f (x) = v(x; 0) + ∂v(x; 0) t = v(x; 0) - E(x; 0)t. ∂t Заметим, что D - симметрический тензор скоростей деформаций в точке x, S - кососимметрический тензор с аксиальным вектором индукции магнитно1 го поля в точке x : Sh = 2 b × h. Так как dy = dh, линеаризованное уравнение относительно приращения h запишется в виде dh 1 = Ah + f = Dh + b × h + f, dt 2 h(0) = 0. Классическая формула Коши решения этого неоднородного уравнения t eA(t-s) f (s) ds, h(t) = 0 где используется матричная экспонента eA(t-s) , не позволяет увидеть физический смысл результата. Поступим следующим образом. Так как тензоры D и S есть функции тензора A, они коммутируют. Значит, матричная экспонента интегрального представления решения расщепляется в произведение матричных экспонент eA(t-s) = eD(t-s) eS(t-s) = eS(t-s) eD(t-s) . Имеем t eS(t-s) g(t, s) ds, h(t) = g(t, s) = eD(t-s) f (s) = eD(t-s) (v(x; 0) - E(x; 0)s), 0 где E(x; 0) - вектор электрической напряжённости в точке x в начальный момент времени t = 0. Пусть D = T (x)diag {ν1 (x), ν2 (x), ν3 (x)}T ∗ (x) - спектральное разложение симметрического тензора. Тогда g(t, s) = T (x)diag{eν1 (x)(t-s) , eν2 (x)(t-s) , eν3 (x)(t-s) }T ∗ (x)(v(x; 0) - E(x; 0)s). Отсюда видно влияние тензора скоростей деформаций на процесс интегрирования линеаризованного уравнения, определяющего малые приращения при формировании течения эфира. Остаётся выяснить влияние кососимметрического тензора. Для этого разложим в степенной ряд матричную экспоненту eS(t-s) : eS(t-s) g = g + 690 1 1 (t - s)Sg + · · · + (t - s)n S n g + · · · . 1! n! Об одном классе векторных полей Заметим, что Sg = 1 b×g, где b = B(x; 0) - вектор индукции магнитного поля 2 в точке x в начальный момент времени t = 0. 1 Далее, S 2 g = 4 b × (b × g). Так как S 2 g = 1 b(b · g) - 1 gb2 , имеем 4 4 1 b2 b2 S 3 g = b × S 2 g = - b × g = - Sg. 2 8 4 Сумму сходящегося степенного ряда найдём, суммируя сначала все нечётные степени отображения S, а потом все чётные степени этого отображения. С учётом найденных закономерностей образования степеней получим выражение для eS(t-s) g в виде eS(t-s) g = α(t, s)b × g + β(t, s)b × (b × g), где скалярные функции α(t, s) и β(t, s) могут быть явно выписаны. Теперь остаётся воспользоваться теоремой о среднем для интеграла t eS(t-s) g(t, s) ds, h(t) = 0 согласно которой для некоторой точки s∗ = s(t) можем написать выражение для зависимости от времени приращения в выражении для функции y(t) = = x + h(t), описывающей течение физического вакуума: h(t) = g(t, s∗ ) + α(t, s∗ )b × g(t, s∗ ) + β(t, s∗ )b × (b × g(t, s∗ )). 5. Нестандартное решение системы уравнений Максвелла. Примеры решений уравнений Максвелла без источников, отличные от найденных Максвеллом плоских волн, широко известны [24-26]. Здесь мы дополним этот список. Для этого рассмотрим «вырожденную форму» трёхпараметрического представления, когда Q(x) = iP (x). Умножение на комплексную единицу i приводит к повороту на угол π/2 вектора P (x) вокруг своей оси. По формуле Эйлера имеем экспоненциальный вид поля смещения эфира: u(x; t) = P (x)eik(x)t . Далее ограничимся случаем k(x) ≡ iµ, где µ - константа любого знака. Стационарный вектор P (x) выберем в виде произведения скалярной функции µ f (x) = e c s , s = αx1 + βx2 + γx3 , α2 + β 3 + γ 2 = 1 и постоянного вектора G = (G1 , G2 , G3 ). Очевидно, что f (x) = Теперь имеем µ f (x) {α, β, γ} c и rot G = 0. u(x; t) = f (x)Ge-µt . Предыдущие выкладки показывают, что для нахождения решения системы уравнений Максвелла без источников достаточно проверить выполнение равенства iµ P (x) = rot P (x). c 691 И с л а м о в Г. Г. Так как rot P (x) = rot(f (x)G) = f (x) rot G + f (x) × G, получим спектральную задачу iµ f (x)G = c После сокращения на f (x) × G. µ f (x) эта задача запишется в простом виде: c iG - {α, β, γ} × G = 0. Это обыкновенная однородная система линейных уравнений относительно координат G1 , G2 , G3 постоянного вектора G. Определитель этой системы i γ -β -γ i α β -α i = i(-1 + α2 + β 2 + γ 2 ) = 0. Значит, эта система имеет ненулевое решение. При γ 2 = 1 подпространство решений одномерно и вектор G= iβ + αγ iα + βγ , ,1 1 - γ2 γ2 - 1 образует базис этого подпространства. В противном случае базис будет двумерным. Его также можно выписать явно. Очевидно, что спектральный параметр µ может быть любым вещественным числом (спектр непрерывный, сплошной, занимает всю числовую прямую). Заметим, что скалярная функция f (x) электрического поля E(x; t) = -µ2 f (x)Ge-µt вещественна, чего не скажешь относительно компонент постоянного вектора G. Таким образом, существует электрическое поле, которое удовлетворяет системе уравнений Максвелла без источников и «заряжается» и «разряжается» по экспоненциальному закону с любой скоростью µ. Заключение. Система уравнений Максвелла без источников нуждается в дальнейшем изучении. Новые локальные и глобальные решения этой системы открывают путь к изучению электромагнитного строения физического вакуума и элементарных частиц.

About the authors

Galimzian G Islamov

Udmurt State University

Email: ggislamov@gmail.com
(Dr. Phis. & Math. Sci.; ggislamov@gmail.com), Head of Department, Dept. of HPC & Parallel programming 1, Universitetskaya str., Izhevsk, 426034, Russian Federation

References

Supplementary files