СЕЛФ

Новогодний вопрос Лео

С.Б. Каравашкин

Письмо Лео ко мне

От: xu xu <xuszxu@yahoo.com>

Кому: Сергей Каравашкин selflab@go.com, selflab@mail.ru

Дата: Thu, 18 Dec 2003 09:56:31 -0800 (PST)

Тема: проблема с формулой дивергенции

 

Уважаемый Сергей,

У меня вопрос по твоей формуле дивергенции. Посмотри пожалуйста прикрепленный файл.

С Рождеством!

Лео

Приложение: problem.pdf

 

figa.gif (3372 bytes)

 

В сферических координатах электрическое поле в точке Р, создаваемое диполем Il , имеет вид [1]:

Image1940.gif (1973 bytes)
где Image1941.gif (1139 bytes)

Естественно, что направление распространения электрического поля  vectorn.gif (845 bytes) = vectorr.gif (839 bytes).. В соответствии с твоей формулой [2]:

Image1943.gif (1107 bytes) (1)
мы имеем
Image1944.gif (1421 bytes)

Однако легко убедиться, что

Image1945.gif (1551 bytes)

Я думаю, что проблема в том, что (1) неверно. Ссылка [2] не вносит вклада стороны поверхности в дивергенцию вектора, как описано в (1).

Литература

[1] James R. Wait. "Introduction to antennas & propagation", Peter Peregrious Ltd, London, Unided Kindom, 1986, pp 89. Equation (5.30) - (5.31)

[2] S.B. Karavashkin "Transformation of divergence theorem in dynamic fields"

 

Мой ответ Лео

Дата: 22.12.03

Уважаемый Лео,

Ваше представление о сферическом диполе, которое Вы почерпнули из Джеймса Вейта, безусловно, похвально. Но я был бы Вам ещё более признателен, если бы Вы прежде, чем сравнивать приведенные Вами выражения для дивергенции потока от сферического диполя с выражением для дивергенции вектора в динамических полях из моей работы, проанализировали бы вывод самих уравнений из Джеймса Вейта. При этом Вы нашли бы много интересного - того, что не видно из конечных формул. Проверьте внимательно, причём не чисто механически, а феноменологически, анализируя, что откуда берётся и почему. Тогда поймёте причину. У меня, к сожалению, нет книги именно Джеймса Вейта, но я могу Вам помочь в анализе, основываясь на решении аналогичной задачи у Кугушева [1]. В этой задаче также исследуется электрическое и магнитное поле провода с током при длине провода  l << lumbdacut.gif (841 bytes) , только в непроводящей среде, т.е. при sigmacut.gif (843 bytes) = 0 . Конечное выражение для Е-компоненты поля у него имеет вид [1, с. 98]

Image1955.gif (1777 bytes)

(1)

где betacut.gif (852 bytes) = omegacut.gif (838 bytes)/c  - волновое число, Im   - амплитуда тока в проводе, epsiloncut.gif (833 bytes)a  - диэлектрическая постоянная,  mycut.gif (843 bytes)a   - магнитная проницаемость. Сравнив (1) с Вашими выражениями для Er  и ETetabigbottom.gif (826 bytes)  при sigmacut.gif (843 bytes) = 0, Вы можете убедиться, что по структуре эти выражения полностью соответствуют друг другу. Поэтому дивергенция, взятая от (1), так же обращается в ноль, как и в Вашем случае.

Чтобы понять причину, обратим сначала наше внимание на тот факт, что решение (1), как и Ваши решения, получено путём подстановки значения векторного потенциала Abigdot.gif (850 bytes)  в стандартную систему уравнений [1, с. 98]

Image1964.gif (1294 bytes)

(2)

т.е. окольным путём. Сначала находим вектор Hbigdot.gif (847 bytes) , а потом уже вектор Ebigdot.gif (842 bytes) . При этом мы дважды берём операцию ротора над выражением для векторного потенциала, поскольку (2) можно представить в виде

Image1967.gif (1191 bytes)

(3)

Таким образом, мы дважды пропускаем выражение (1) через фильтр, отсекающий продольную компоненту. Факт, что в результирующем выражении эта компонента сохраняется, обусловлен тем, что, вдобавок ко всему, мы теряем очень важное слагаемое в выражении для электрического поля. Действительно, если согласно [2, с. 41], до модельных расчётов мы подставим первое выражение из (2) в закон Фарадея, то при этом получим:

Image1968.gif (1154 bytes)

(4)
или

Image1969.gif (1239 bytes)

(5)

"Последнее уравнение показывает, что вектор

Image1970.gif (1136 bytes)

остаётся потенциальным вектором, т.е. может быть представлен в виде

Image1971.gif (1279 bytes)

(6)

где ficut.gif (844 bytes) - функция координат и времени, которую мы будем именовать скалярным потенциалом" [2, с. 42]. Обратим внимание, что правая часть (6) равна не нулю, но некоторой функции КООРДИНАТ и ВРЕМЕНИ. Поэтому мы не имеем права обращать произвольно для себя в ноль ВЕКТОР  grad ficut.gif (844 bytes) , как это принято в настоящее время (см. например [3, с. 72]). Полное обоснование этого тезиса приведено в моей работе [4]. Но даже и без строгого математического доказательства понятно, что пренебрегать указанным членом нельзя. Действительно, если мы представим себе ЭМ волну с частотой в несколько сотых герца (такие волны реально существуют в пространстве и учитываются астрофизиками при исследовании межзвёздных туманностей), и поднимем вопрос о присутствии скалярного потенциала в данном поле, то учитывая огромную протяжённость длины волны, мы обязаны признать его наличие, как должны признать и его переменный во времени характер. Можем ли мы пренебрегать этим членом, если он является частью электрического поля, причём той частью, которая исключительно определяет потенциальный характер этого поля? Безусловно, нет. Те компоненты поля, которые действительно малы или не влияют на общий характер поля, сами "отсеиваются" и не требуют некоторых насильственных калибровок. Вместе с тем, возвращаясь к стандартному выводу на основе системы (2), мы теперь видим, что двойная операция взятия ротора от выражения для скалярного аргумента изъяла из функции электрического поля именно grad ficut.gif (844 bytes), что существенно исказило результат. А второй путь через закон Фарадея не уничтожает указанного вектора, поскольку самой операции взятия ротора от функции при этом мы не производим! И это я сейчас покажу.

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

Image1974.gif (1152 bytes)

(7)

При этом второе слагаемое, связывающее электрическое поле с магнитным, выражает закон электромагнитной индукции Фарадея" [2, с. 42].

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

Чтобы на этой основе определить правильную зависимость электрического поля от времени и координат, мы должны сделать ещё одно отступление от изложения и обратить внимание на саму запись векторного потенциала. Согласно Кугушеву, выражение для векторного потенциала, на основе которого он получил выражение (1) для напряжённости электрического поля, имеет вид [1, с. 98]:

Image1975.gif (1385 bytes)

(8)

При этом следует отметить, что данное выражение не является точным решением волнового уравнения для векторного потенциала, что очень важно для дальнейшего изложения. Действительно, как известно, выражение для векторного потенциала находится на основе волнового уравнения типа [1, с. 88]

Image1976.gif (1054 bytes)

(9)

где mywave.gif (850 bytes)a    - комплексная магнитная проницаемость.

Но решением этого уравнения является не функция (8), которая нами исследуется, а интегральная функция типа [1. с. 89], [3. с. 95]:

Image1978.gif (1246 bytes)

(10)

Однако "для области r >> l  можно считать все точки диполя равноотстоящими от исследуемой точки поля, вследствие чего (10) можно переписать в виде

Image1980.gif (1292 bytes)

(11)

где  e3 орт, параллельный направлению тока. Так как jm = Im/S и  Sl = V , S   - сечение проводника, получаем

Image1985.gif (1187 bytes)

(12)

Переходя к сферической системе координат, в которой

Image1986.gif (1086 bytes)

(13)

получаем формулу (8), определяющую векторный потенциал на расстоянии r >> l " [1. с. 97]. Таким образом, (8) определяет векторный потенциал только для дальней зоны, когда согласно [5], все параметры волнового процесса стабилизировались. Поэтому, когда исследователи (в том числе и сам Кугушев) применяют для ближней зоны результаты, полученные на основе данного векторного потенциала, то это они делают неправомерно, поскольку этим они нарушают условия задачи. Отсюда возникает сомнение в правильности выражения (1). Если выражение (1) описывает процесс в дальней зоне, то никаких суперпозиций напряжённости поля с разными характеристиками затухания и разными сдвигами фаз в выражении присутствовать не должно. Это дополнительно подтверждает высказанное выше утверждение о том, что вывод выражения (1) был произведен некорректно.

Чтобы получить правильное выражение для напряжённости электрического поля на основе выражения (7), нам нужно определить скалярный потенциал ficut.gif (844 bytes). Это мы можем сделать на основе зависимости, полученной в работе [2, с. 106]:

Image1987.gif (936 bytes)

(14)

Учитывая, что в нашей задаче

Image1988.gif (916 bytes)

(15)

мы можем записать

Image1989.gif (922 bytes)

(16)

или с учётом (8)

Image1990.gif (1217 bytes)

(17)

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

Image1991.gif (1532 bytes)

(18)

Image1992.gif (1521 bytes)

(19)

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

Image1993.gif (1544 bytes)

(20)

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

Теперь, на основе скорректированного, а главное, полного выражения для напряжённости электрического поля, сравним дивергенцию напряжённости этого поля с результатами доказанной мною теоремы о дивергенции вектора в динамических полях [6] (странно только, что, хорошо оформив свою заметку, Вы не заметили, что эта статья опубликована в международном журнале Archivum mathematicum, 37 (2001), 3, с.233-243. Ссылку без издания неприлично представлять, как Вы думаете? ;-)).

Для проверки воспользуемся стандартной формулой, приведенной в Вашей заметке:

Image1994.gif (1512 bytes)

(21)

С учётом (20), первый член правой части будет иметь вид

Image1995.gif (1468 bytes)

(22)

Соответственно второй член примет вид

Image1996.gif (1563 bytes)

(23)

В результате имеем

Image1997.gif (1310 bytes)

(24)

Как мы видим из (24), дивергенция электрического поля не равна нулю и фаза её изменения во времени смещена на 90о по отношению к фазе продольной составляющей, что вполне соответствует феноменологии распространения волнового процесса в пространстве.

Несложно проверить, что и амплитудное значение дивергенции вектора напряжённости электрического поля в Вашей задаче тоже соответствует доказанной мною теореме о дивергенции вектора в динамических полях. Для этого воспользуемся выражением (1) Вашей заметки, которое Вы вполне правильно интерпретировали, и продолжим его до результата:

Image1998.gif (1704 bytes)

(25)

Так что можете с превеликим удовольствием сравнивать (24) и (25). Как видите, в действительности проблем с теоремой нет. Есть проблемы только с существующими решениями для конкретных моделей.

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

Если же вернуться к разговору о справедливости теоремы, доказанной мною, то её можно опровергнуть только, найдя ошибку в самом доказательстве, поскольку теорема первична по отношению ко всем решениям моделирующих уравнений. И это не моё желание. Это объективность построения любой теории. Сначала формулируются базовые теоремы (в данном случае - теоремы сохранения для динамических полей), потом строится математический аппарат, а потом решаются задачи по моделированию конкретных процессов. Так строил свою систему Максвелл, так нужно и дальше развивать теорию ЭМ полей.

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

С новым годом Вас и с новым знанием.

Сергей.

 

Литература:

1. Кугушев А.М., Голубева Н.С. Основы радиоэлектроники. Линейные электромагнитные процессы. Москва, Энергия, 1969, 880 с.

2. Левич В.Г. Курс теоретической физики, т.1. М., Государственное издательство физико-математической литературы, 1962, 695 с.

3. Ландау Л.Д., Лифшиц Е.М. Теория поля. - // Теоретическая физика, т. 2. М., Наука, 1973, 504 с.

4. Каравашкин С.Б. К вопросу о продольных электромагнитных волнах. Глава 1. Снятие запретов. SELF Transactions, 1 (1994), 15-47,

http://angelfire.lycos.com/la3/selftrans/archive/archive.html#long

5. Каравашкин С.Б., Каравашкина О.Н. Сравнение характеристик ЭМ и акустического полей. SELF Transactions, 2 (2002),

http://angelfire.lycos.com/la3/selftrans/v3_1/contents3.html#taew

6. Каравашкин С.Б. Трансформация теоремы о дивергенции в динамических полях. Archivum mathematicum, 37 (2001), 3, 233-243. SELF Transactions, archive,

  http://angelfire.lycos.com/la3/selftrans/archive/archive.html#div

 

Hosted by uCoz