Жидкокристаллический осмос

и другие научные идеи

Содержание:

Жидкокристаллический двигатель
или
акустомеханический преобразователь

А.Ю. Дроздов

4.01.2010 - 17.07.2010

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

word   htm

Ниже приводится полемика по данной работе:

Принципиальные вопросы, поставленные
редактором журнала «Химия и Жизнь»
и ответы на них

17.11.2010

word   htm

Жидкокристаллический осмос
или
о возможности нарушения принципа детального равновесия в жидкокристаллической дисклинации

А.Ю. Дроздов

17.11.2010

В данной работе ставится под сомнение всеохватность принципа детального равновесия. И предлагается идея экспериментальной проверки возможности нарушения этого принципа в системе жидкий кристалл - мембрана.

word   htm

Гравитационно-осмотический кольцар

А.Ю.Дроздов

16.11.2010

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

word   htm

  Чудо-дерево породы «БИОСОЛЯР»

Юный Техник № 8 за 1986 г.

Магнитный двигатель на эффектах скалярного магнитного поля

А.Ю.Дроздов

1.11.2017

В данной работе предлагается конфигурация двигателя на постоянных магнитах, основанная на конфигурации опыта N 8 из книги Николаева.

word   htm

Электродинамический расчёт рельсотрона

с использованием понятий векторного потенциала, магнитного поля и силы Лоренца

А.Ю.Дроздов

02.04.2018

Электродинамический расчёт рельсотрона

с использованием понятий векторного потенциала, магнитного поля и силы Лоренца

а также формулы для силы действующей на движущийся заряд взятой из "Новой электродинамики" Ф.Ф.Менде

уточнённый привлечением понятий скалярного магнитного поля и силы Николаева

с последующим выводом

формулы Взаимодействия элементов тока

 

А.Ю.Дроздов

02-19.04.2018

Опыт Николаева и Дейны с цилиндрами

 

А.Ю.Дроздов

14.08.2019-14.06.2020

'

К вопросу об интерретации результатов эксперимента Майкельсона Морли

Michelson-Morley-0

 

Michelson-Morley-90

 

Michelson-Morley-180

 

Michelson-Morley-270

 

А.Ю.Дроздов

24.11.2020

 

 

Электродинамический расчёт рельсотрона  

с использованием понятий векторного потенциала, магнитного поля и силы Лоренца 

а также формулы для силы действующей на движущийся заряд взятой из "Новой электродинамики" Ф.Ф.Менде 

уточнённый привлечением понятий скалярного магнитного поля и силы Николаева 

с последующим выводом формулы 

Взаимодействия элементов тока 

 

А.Ю.Дроздов 

 

В предыдущей работе [1] было показано, что в конфигурации рельсотрона использование формул Лоренца при вычислении взаимодействия отдельных токовых элементов друг с другом не позволяет корректно вычислить место приложения сил реакции (отдача рельсотрона). Данная проблема была поднята автором на научных форумах [2] и [3].  

 

В результате ее обсуждения был сформулирован план эксперимента, который был осуществлён автором [4]. Впоследствии автором был осуществлён ещё один эксперимент с жидким металлом [5]. Результаты проведенных экспериментов по мнению автора показывают наличие продольной силы реакции приложенной к рельсам.  

 

В связи с этим в данной работе опробована усовершенствованная методика электродинамического расчёта рельсотрона с привлечением понятия скалярного магнитного поля Николаева [7], с помощью которой место приложения сил реакции может быть вычислено. В результате обнаружено, что полученная сила взаимодействия отдельных токовых элементов удовлетворяет третьему закону Ньютона. Обобщая полученный результат была выведена общая формула взаимодействия элементов тока, которая отличается от уравнений Лоренца-Грассмана [8]. Применимость полученной формулы обоснована ссылками на рассуждения  Е.И. Тамма. Кроме того данная формула была протестирована при рассчёте силы взаимодействия бесконечных параллельных проводников.  

 

Векторный потенциал в точке x, y, z, создаваемый участком тока `#mover(mi( , расположенном в точке x__i1, y__i1, z__i1 

dA_ := proc (x, y, z, t, I__1, x__i1, y__i1, z__i1, dl_) options operator, arrow; `+`(`/`(`*`(`/`(1, 4), `*`(I__1, `*`(dl_))), `*`(Pi, `*`(sqrt(`+`(`*`(`^`(`+`(x, `-`(x__i1)), 2)), `*`(`^`(`+`(y, `-`(... 

`#mover(mi( (1)
 


Теперь попробуем решить задачу вычисления отдачи рельсотрона
 

Drawing-Canvas 

 

 

X1_peremychka := 1.1 

X1_peremychka := 1.1 (3)
 

X2_peremychka := 1.9 

X2_peremychka := 1.9 (4)
 

Y_peremychka := 0 

Y_peremychka := 0 (5)
 

Z_relsotron := 0 

Z_relsotron := 0 (6)
 

 

Y_istochnik := -10 

Y_istochnik := -10 (7)
 

 

X_left_rail := 1 

X_left_rail := 1 (8)
 

X_right_rail := 2 

X_right_rail := 2 (9)
 

 

Y1_rails := -10 

Y1_rails := -10 (10)
 

Y2_rails := 0 

Y2_rails := 0 (11)
 

Рассчитаем векторный потенциал создаваемый перемычкой, то есть снарядом, которую мы расположим вдоль оси икс, пусть ток по ней течёт справа налево, длина перевычки равна 1 

 

A_peremychka_ := subs(y__i1 = Y_peremychka, z__i1 = Z_relsotron, int(dA_(x, y, z, t, I__r, x__i1, y__i1, z__i1, `+`(`-`(`#mover(mi( 

`#mover(mi(
`#mover(mi(
(12)
 

 


 

 


Векторный потенциал, создаваемый правой рельсой
 

 

A_right_rail_ := subs(x__i1 = X_right_rail, z__i1 = Z_relsotron, int(dA_(x, y, z, t, I__r, x__i1, y__i1, z__i1, `#mover(mi( 

`#mover(mi( (13)
 

 

Векторный потенциал, создаваемый левой рельсой 

 

A_left_rail_ := subs(x__i1 = X_left_rail, z__i1 = Z_relsotron, int(dA_(x, y, z, t, I__r, x__i1, y__i1, z__i1, `+`(`-`(`#mover(mi( 

`#mover(mi( (14)
 

 

 

 

Суммарный векторный потенциал, создаваемый правой и левой рельсой 

A_rails_ := simplify(`+`(A_right_rail_, A_left_rail_)) 

`#mover(mi(
`#mover(mi(
(15)
 

 


Магнитное поле создаваемое правой рельсой в области перемычки
 

H_of_right_rail_on_peremychka_ := subs(y = Y_peremychka, z = Z_relsotron, Curl(A_right_rail_)); 1; H_of_left_rail_on_peremychka_ := subs(y = Y_peremychka, z = Z_relsotron, Curl(A_left_rail_)); 1
H_of_right_rail_on_peremychka_ := subs(y = Y_peremychka, z = Z_relsotron, Curl(A_right_rail_)); 1; H_of_left_rail_on_peremychka_ := subs(y = Y_peremychka, z = Z_relsotron, Curl(A_left_rail_)); 1
 

 

`#mover(mi(

`#mover(mi( (16)
 

 

 

 

 

Рассчитаем векторный потенциал создаваемый источником тока, который мы расположим вдоль оси икс, ток по ней течёт с лева на право, длина равна 1 

 

A_istochnik_ := subs(y__i1 = Y_istochnik, z__i1 = Z_relsotron, int(dA_(x, y, z, t, I__r, x__i1, y__i1, z__i1, `#mover(mi( 

`#mover(mi( (17)
 

 

Находим также магнитное поле создаваемое источником тока в области перемычки 

 

H_of_istochnik_on_peremychka_ := subs(y = Y_peremychka, z = Z_relsotron, Curl(A_istochnik_)); 1 

`#mover(mi( (18)
 

 

Применяем формулу силы Лоренца - находим таким образом выражение для силы, с которой магнитное поле источника действует на перемычку 

f_lorenz_from_istochnik_to_peremychka_ := simplify(`*`(mu, `*`(Typesetting:-delayCrossProduct(`+`(`-`(`*`(I__r, `*`(`#mover(mi( 

`#mover(mi( (19)
 

 

Интегрируем по длине перемычки 

 

 

F_lorenz_from_istochnik_to_peremychka_ := evalf(simplify(int(f_lorenz_from_istochnik_to_peremychka_, x = X_left_rail .. X_right_rail))); 1 

`#mover(mi( (20)
 

 

 

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

 

f_lorenz_from_right_rail_to_peremychka_ := `*`(mu, `*`(Typesetting:-delayCrossProduct(`+`(`-`(`*`(I__r, `*`(`#mover(mi( 

`#mover(mi( (21)
 

Интегрируем по длине перемычки 

 

F_lorenz_from_right_rail_to_peremychka_ := int(f_lorenz_from_right_rail_to_peremychka_, x = X1_peremychka .. X2_peremychka); 1 

`#mover(mi( (22)
 


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

 

 

f_lorenz_from_left_rail_to_peremychka_ := `*`(mu, `*`(Typesetting:-delayCrossProduct(`+`(`-`(`*`(I__r, `*`(`#mover(mi( 

`#mover(mi( (23)
 

Интегрируем по длине перемычки 

 

F_lorenz_from_left_rail_to_peremychka_ := int(f_lorenz_from_left_rail_to_peremychka_, x = X1_peremychka .. X2_peremychka); 1 

`#mover(mi( (24)
 

 

 

 

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


Магнитное поле, создаваемое перемычкой в области правой рельсы и в области левой рельсы
 

 

H_of_peremychka_on_right_rail_ := subs(x = X_right_rail, z = Z_relsotron, Curl(A_peremychka_)); 1; H_of_peremychka_on_left_rail_ := subs(x = X_left_rail, z = Z_relsotron, Curl(A_peremychka_))
H_of_peremychka_on_right_rail_ := subs(x = X_right_rail, z = Z_relsotron, Curl(A_peremychka_)); 1; H_of_peremychka_on_left_rail_ := subs(x = X_left_rail, z = Z_relsotron, Curl(A_peremychka_))
 

 

`#mover(mi(
`#mover(mi(
`#mover(mi(

`#mover(mi(
`#mover(mi(
`#mover(mi(
(25)
 

 

 

 

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

 

f_lorenz_from_peremychka_to_right_rail_ := `*`(mu, `*`(Typesetting:-delayCrossProduct(`*`(I__r, `*`(`#mover(mi(
f_lorenz_from_peremychka_to_right_rail_ := `*`(mu, `*`(Typesetting:-delayCrossProduct(`*`(I__r, `*`(`#mover(mi(
f_lorenz_from_peremychka_to_right_rail_ := `*`(mu, `*`(Typesetting:-delayCrossProduct(`*`(I__r, `*`(`#mover(mi(
f_lorenz_from_peremychka_to_right_rail_ := `*`(mu, `*`(Typesetting:-delayCrossProduct(`*`(I__r, `*`(`#mover(mi(
 

 

`#mover(mi(
`#mover(mi(
`#mover(mi(
`#mover(mi(

`#mover(mi(
`#mover(mi(
`#mover(mi(
(26)
 

 

 

F_lorenz_from_peremychka_to_right_rail_ := int(f_lorenz_from_peremychka_to_right_rail_, y = Y1_rails .. Y2_rails); 1 

 

 

F_lorenz_from_peremychka_to_left_rail_ := int(f_lorenz_from_peremychka_to_left_rail_, y = Y1_rails .. Y2_rails); 1 


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

 

Использование подхода Ф.Ф.Менде 

 

В работе Ф.Ф.Менде [6] приводится формула (2.6), где определены все силы, действующие на заряд, движущийся в поле векторного потенциала. Применим эту формулу для решения задачи.  

 

Рассмотрим третье слагаемое формулы 2.6 : градиент потенциальной энергии заряда, движущегося в поле векторного потенциала, как градиент скалярного произведения векторного потенциала на вектор тока второго проводника. Вычислим градиент потенциальной энергии тока рельс и источника в поле векторного потенциала перемычки 

grad_P_of_right_rail_on_peremychka_A_ := `*`(mu, `*`(I__r, `*`(Gradient(Typesetting:-delayDotProduct(`#mover(mi(
grad_P_of_right_rail_on_peremychka_A_ := `*`(mu, `*`(I__r, `*`(Gradient(Typesetting:-delayDotProduct(`#mover(mi(
grad_P_of_right_rail_on_peremychka_A_ := `*`(mu, `*`(I__r, `*`(Gradient(Typesetting:-delayDotProduct(`#mover(mi(
grad_P_of_right_rail_on_peremychka_A_ := `*`(mu, `*`(I__r, `*`(Gradient(Typesetting:-delayDotProduct(`#mover(mi(
 

 

 

`#mover(mi(

`#mover(mi(

`#mover(mi(
`#mover(mi(
`#mover(mi(
`#mover(mi(
`#mover(mi(
`#mover(mi(
`#mover(mi(
`#mover(mi(
`#mover(mi(
(27)
 

 

Как видим рельсы в поле векторного потенциала перемычки не испытывают влияние градиента потенциальной энергии тока в поле векторного потенциала, но противоположная перемычка (источник тока) испытывает. Поэтому интегрируя по длине источника тока этот градиент, получаем в интеграле две компоненты силы продольную и поперечную  

 

Int_dx_grad_P_of_istochnik_on_peremychka_A_ := int(grad_P_of_istochnik_on_peremychka_A_, x) 

`#mover(mi(
`#mover(mi(
`#mover(mi(
(28)
 

 

Подставляем пределы интегрирования 

 

Int_on_istochhik_len_of_grad_P_of_istochnik_on_peremychka_A_ := `+`(subs(x = X_right_rail, int_dx_grad_P_of_istochnik_on_peremychka_A_), `-`(subs(x = X_left_rail, int_dx_grad_P_of_istochnik_on_peremyc...
Int_on_istochhik_len_of_grad_P_of_istochnik_on_peremychka_A_ := `+`(subs(x = X_right_rail, int_dx_grad_P_of_istochnik_on_peremychka_A_), `-`(subs(x = X_left_rail, int_dx_grad_P_of_istochnik_on_peremyc...
 

`#mover(mi( (29)
 

 

После упрощения видим, что осталась только поперечная компонента силы.  

 

evalf(Int_on_istochhik_len_of_grad_P_of_istochnik_on_peremychka_A_) 

0. (30)
 

 

 

 

Исследуем теперь второе слагаемое формулы 2.6. Записываем выражение для конвективной производной 

 

ConvectiveDerivative := proc (A_, B_) options operator, arrow; `+`(`*`(Typesetting:-delayDotProduct(A_, Gradient(Component(`#mover(mi(
ConvectiveDerivative := proc (A_, B_) options operator, arrow; `+`(`*`(Typesetting:-delayDotProduct(A_, Gradient(Component(`#mover(mi(
 


(31)
 

 

 

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

 

Conv_deriv_to_istochnik_on_peremychka_A_ := `+`(`-`(`*`(mu, `*`(I__r, `*`(ConvectiveDerivative(`#mover(mi( 

 

Conv_deriv_to_right_rail_on_peremychka_A_ := `+`(`-`(`*`(mu, `*`(I__r, `*`(ConvectiveDerivative(`#mover(mi( 

 

Conv_deriv_to_left_rail_on_peremychka_A_ := `+`(`-`(`*`(mu, `*`(I__r, `*`(ConvectiveDerivative(`+`(`-`(`#mover(mi( 

 

 

Conv_f_to_istochnik_on_peremychka_A_ := subs(z = Z_relsotron, y = Y_istochnik, Conv_deriv_to_istochnik_on_peremychka_A_); -1 

 

Conv_f_to_right_rail_on_peremychka_A_ := subs(z = Z_relsotron, x = X_right_rail, Conv_deriv_to_right_rail_on_peremychka_A_) 

`#mover(mi(
`#mover(mi(
`#mover(mi(
`#mover(mi(
(32)
 

 

Conv_f_to_left_rail_on_peremychka_A_ := subs(z = Z_relsotron, x = X_left_rail, Conv_deriv_to_left_rail_on_peremychka_A_) 

`#mover(mi(
`#mover(mi(
`#mover(mi(
(33)
 

 

Интегрируем по длине источника тока 

Int_dx_Conv_f_to_istochnik_on_peremychka_A_ := int(simplify(Conv_f_to_istochnik_on_peremychka_A_), x); -1 

 

Подставляем пределы интегрирования 

 

Int_on_istochhik_len_of_Conv_f_to_istochnik_on_peremychka_A_ := evalf(`+`(subs(x = X_right_rail, Int_dx_Conv_f_to_istochnik_on_peremychka_A_), `-`(subs(x = X_left_rail, Int_dx_Conv_f_to_istochnik_on_p...
Int_on_istochhik_len_of_Conv_f_to_istochnik_on_peremychka_A_ := evalf(`+`(subs(x = X_right_rail, Int_dx_Conv_f_to_istochnik_on_peremychka_A_), `-`(subs(x = X_left_rail, Int_dx_Conv_f_to_istochnik_on_p...
 

`#mover(mi( (34)
 

Конвективную силу приложенную к источику можно считать равной нулю.

Интегрируем по длине рельс
 

 

Int_dy_Conv_f_to_right_rail_on_peremychka_A_ := int(Conv_f_to_right_rail_on_peremychka_A_, y); -1 

 

Int_dy_Conv_f_to_left_rail_on_peremychka_A_ := int(Conv_f_to_left_rail_on_peremychka_A_, y); -1 

 

Int_on_rail_len_of_Conv_f_to_right_rail_on_peremychka_A_ := evalf(`+`(subs(y = Y2_rails, Int_dy_Conv_f_to_right_rail_on_peremychka_A_), `-`(subs(y = Y1_rails, Int_dy_Conv_f_to_right_rail_on_peremychka...
Int_on_rail_len_of_Conv_f_to_right_rail_on_peremychka_A_ := evalf(`+`(subs(y = Y2_rails, Int_dy_Conv_f_to_right_rail_on_peremychka_A_), `-`(subs(y = Y1_rails, Int_dy_Conv_f_to_right_rail_on_peremychka...
 

Error, numeric exception: division by zero
 

 

 

Int_on_rail_len_of_Conv_f_to_left_rail_on_peremychka_A_ := evalf(`+`(subs(y = Y2_rails, Int_dy_Conv_f_to_left_rail_on_peremychka_A_), `-`(subs(y = Y1_rails, Int_dy_Conv_f_to_left_rail_on_peremychka_A_...
Int_on_rail_len_of_Conv_f_to_left_rail_on_peremychka_A_ := evalf(`+`(subs(y = Y2_rails, Int_dy_Conv_f_to_left_rail_on_peremychka_A_), `-`(subs(y = Y1_rails, Int_dy_Conv_f_to_left_rail_on_peremychka_A_...
 

Error, numeric exception: division by zero
 

 

К сожалению, рассчёт не удался, то тем не менее видно, что конвективная сила, приложенная к рельсам направлена вдоль оси x, а в рамках данного исследования эта компоненты силы нас не интересует. 

 

 

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

Вычислим градиент потенциальной энергии тока верхней в поле векторного потенциала рельс и источника тока 

 

grad_P_of_peremychka_on_rails_A_ := `*`(mu, `*`(I__r, `*`(Gradient(`+`(`-`(Typesetting:-delayDotProduct(`#mover(mi( 

`#mover(mi( (35)
 

grad_P_of_peremychka_on_istochnik_A_ := subs(z = Z_relsotron, y = Y_peremychka, `*`(mu, `*`(I__r, `*`(Gradient(`+`(`-`(Typesetting:-delayDotProduct(`#mover(mi( 

`#mover(mi(
`#mover(mi(
(36)
 

 

Интегрируем 

 

Int_dx_of_grad_P_of_peremychka_on_istochnik_A_ := int(grad_P_of_peremychka_on_istochnik_A_, x) 

`#mover(mi(
`#mover(mi(
(37)
 

 

Применяя пределы интегрирования - длину перемычки, находим потенциальную компоненту силы, действующую на перемычку со стороны источника тока 

 

Int_on_peremychka_len_of_grad_P_of_peremychka_on_istochnik_A_ := evalf(`+`(subs(x = X2_peremychka, Int_dx_of_grad_P_of_peremychka_on_istochnik_A_), `-`(subs(x = X1_peremychka, Int_dx_of_grad_P_of_pere...
Int_on_peremychka_len_of_grad_P_of_peremychka_on_istochnik_A_ := evalf(`+`(subs(x = X2_peremychka, Int_dx_of_grad_P_of_peremychka_on_istochnik_A_), `-`(subs(x = X1_peremychka, Int_dx_of_grad_P_of_pere...
 

`#mover(mi( (38)
 

 

Вычисляем конвективную компоненту силы, действующую на перемычку со стороны источника и рельс 

 

Conv_deriv_of_peremychka_on_istochnik_A_ := subs(z = Z_relsotron, y = Y_peremychka, `+`(`-`(`*`(mu, `*`(I__r, `*`(ConvectiveDerivative(`+`(`-`(`#mover(mi( 

`#mover(mi( (39)
 

 

 

Conv_deriv_of_peremychka_on_rails_A_ := simplify(subs(z = Z_relsotron, y = Y_peremychka, `+`(`-`(`*`(mu, `*`(I__r, `*`(ConvectiveDerivative(`+`(`-`(`#mover(mi( 

`#mover(mi(
`#mover(mi(
(40)
 

 

Интегрируем по длине перемычки 

 

Int_dx_Conv_deriv_of_peremychka_on_istochnik_A_ := int(Conv_deriv_of_peremychka_on_istochnik_A_, x) 

`#mover(mi( (41)
 

 

Int_dx_Conv_deriv_of_peremychka_on_rails_A_ := int(Conv_deriv_of_peremychka_on_rails_A_, x) 

`#mover(mi( (42)
 

Применяем пределы интегрирования 

 

Int_peremychka_len_Conv_deriv_of_peremychka_on_istochnik_A_ := evalf(`+`(subs(x = X2_peremychka, Int_dx_Conv_deriv_of_peremychka_on_istochnik_A_), `-`(subs(x = X1_peremychka, Int_dx_Conv_deriv_of_pere...
Int_peremychka_len_Conv_deriv_of_peremychka_on_istochnik_A_ := evalf(`+`(subs(x = X2_peremychka, Int_dx_Conv_deriv_of_peremychka_on_istochnik_A_), `-`(subs(x = X1_peremychka, Int_dx_Conv_deriv_of_pere...
 

`#mover(mi( (43)
 

 

 

 

Int_peremychka_len_Conv_deriv_of_peremychka_on_rails_A_ := evalf(`+`(subs(x = X2_peremychka, Int_dx_Conv_deriv_of_peremychka_on_rails_A_), `-`(subs(x = X1_peremychka, Int_dx_Conv_deriv_of_peremychka_o...
Int_peremychka_len_Conv_deriv_of_peremychka_on_rails_A_ := evalf(`+`(subs(x = X2_peremychka, Int_dx_Conv_deriv_of_peremychka_on_rails_A_), `-`(subs(x = X1_peremychka, Int_dx_Conv_deriv_of_peremychka_o...
 

`#mover(mi( (44)
 

 

Таким образом мы видим, что применение формулы 2.6 из электродинамики Менде показало что конвективная сила приложенная к перемычке в поле векторного потенциала рельс численно равна силе Лоренца, приложенной к перемычке в магнитном поле, создаваемом рельсами. Потенциальная сила приложенная к перемычке в поле векторного потенциала источника тока численно равна и противоположна по направлению потенциальной силе, приложенной к источнику в поле векторного потенциала перемычки. Потенциальная сила, приложенная к рельсам в поле векторного потенциала перемычки равна нулю, а конвективная сила, приложенная к рельсам направлена перпендикулярно рельсам. То есть формула 2.6 не позволяет рассчитать отдачу рельсотрона приложенную к рельсам, что указывает на неполноту формулы 2.6 

 

 

Учёт скалярного магнитного поля 

 

Попытаемся дополнить полученное решение, применив подход Геннадия Николаева.  

 

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

 

 

`+`(`-`(Divergence(A_peremychka_))) 

`/`(`*`(`*`(0.7957747157e-1, `+`(`*`(10., `*`(sqrt(`+`(`*`(100.0000000, `*`(`^`(x, 2))), `*`(100.0000000, `*`(`^`(y, 2))), 361.0000000, `*`(100.0000000, `*`(`^`(z, 2))), `-`(`*`(380.0000000, `*`(x))))...
`/`(`*`(`*`(0.7957747157e-1, `+`(`*`(10., `*`(sqrt(`+`(`*`(100.0000000, `*`(`^`(x, 2))), `*`(100.0000000, `*`(`^`(y, 2))), 361.0000000, `*`(100.0000000, `*`(`^`(z, 2))), `-`(`*`(380.0000000, `*`(x))))...
`/`(`*`(`*`(0.7957747157e-1, `+`(`*`(10., `*`(sqrt(`+`(`*`(100.0000000, `*`(`^`(x, 2))), `*`(100.0000000, `*`(`^`(y, 2))), 361.0000000, `*`(100.0000000, `*`(`^`(z, 2))), `-`(`*`(380.0000000, `*`(x))))...
`/`(`*`(`*`(0.7957747157e-1, `+`(`*`(10., `*`(sqrt(`+`(`*`(100.0000000, `*`(`^`(x, 2))), `*`(100.0000000, `*`(`^`(y, 2))), 361.0000000, `*`(100.0000000, `*`(`^`(z, 2))), `-`(`*`(380.0000000, `*`(x))))...
(45)
 

 

 

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

 

with(plots); -1; contourplot(subs(z = Z_relsotron, I__r = 1, `+`(`-`(Divergence(A_peremychka_)))), x = 0 .. 3, y = -10 .. 1); 1 

Plot_2d
 

 

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

 

DivA_peremychka_right := subs(x = X_right_rail, z = Z_relsotron, Divergence(A_peremychka_)); 1; DivA_peremychka_left := subs(x = X_left_rail, z = Z_relsotron, Divergence(A_peremychka_))
DivA_peremychka_right := subs(x = X_right_rail, z = Z_relsotron, Divergence(A_peremychka_)); 1; DivA_peremychka_left := subs(x = X_left_rail, z = Z_relsotron, Divergence(A_peremychka_))
 

 

DivA_peremychka_right := `+`(`-`(`/`(`*`(`*`(0.7957747157e-1, `+`(`*`(10., `*`(sqrt(`+`(1, `*`(100.0000000, `*`(`^`(y, 2))))))), `-`(`*`(10., `*`(sqrt(`+`(81.0000000, `*`(100.0000000, `*`(`^`(y, 2))))...

DivA_peremychka_left := `+`(`-`(`/`(`*`(`*`(0.7957747157e-1, `+`(`*`(10., `*`(sqrt(`+`(81.0000000, `*`(100.0000000, `*`(`^`(y, 2))))))), `-`(`*`(10., `*`(sqrt(`+`(1, `*`(100.0000000, `*`(`^`(y, 2)))))... (46)
 

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

plot(subs(I__r = 1, DivA_peremychka_right), y = Y1_rails .. Y2_rails) 

Plot_2d
 

 

plot(subs(I__r = 1, DivA_peremychka_left), y = Y1_rails .. Y2_rails) 

Plot_2d
 

 

 

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

 

f_nikolaev_from_peremychka_to_right_rail_ := `+`(`-`(`*`(mu, `*`(I__r, `*`(`#mover(mi(
f_nikolaev_from_peremychka_to_right_rail_ := `+`(`-`(`*`(mu, `*`(I__r, `*`(`#mover(mi(
f_nikolaev_from_peremychka_to_right_rail_ := `+`(`-`(`*`(mu, `*`(I__r, `*`(`#mover(mi(
f_nikolaev_from_peremychka_to_right_rail_ := `+`(`-`(`*`(mu, `*`(I__r, `*`(`#mover(mi(
 

 

`#mover(mi(

`#mover(mi( (47)
 

 



Интегрируем силу Николаева по длине рельсы
 

 

F_nikolaev_from_peremychka_to_right_rail_ := int(f_nikolaev_from_peremychka_to_right_rail_, y = Y1_rails .. Y2_rails); 1 

`#mover(mi( (48)
 

F_nikolaev_from_peremychka_to_left_rail_ := int(f_nikolaev_from_peremychka_to_left_rail_, y = Y1_rails .. Y2_rails); 1 

`#mover(mi( (49)
 

 

 

 

Сравниваем с силой Лоренца, с которой магнитное поле рельс действует на перемычку 

 

F_right_peremychka_ := int(f_right_peremychka_, x = X1_peremychka .. X2_peremychka); 1 

`#mover(mi( (50)
 

 

F_left_peremychka_ := int(f_left_peremychka_, x = X1_peremychka .. X2_peremychka); 1 

`#mover(mi( (51)
 

 

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

 

Расчёт для малых элементов тока 


 

Мы можем уточнить этот результат расчётом не в интегральной, а в дифференциальной форме 

 

 

X_center_peremychka := `/`(3, 2) 

X_center_peremychka := `/`(3, 2) (52)
 


Y_center_relsy := -5
 

Y_center_relsy := -5 (53)
 

 


Векторный потенциал, создаваемой бесконечно малым отрезком тока центральной части перемычки
 

 

dA_from_peremychka_center_ := subs(y__i1 = Y_peremychka, z__i1 = Z_relsotron, x__i1 = X_center_peremychka, dA_(x, y, z, t, I__r, x__i1, y__i1, z__i1, `+`(`-`(`#mover(mi( 

`#mover(mi( (54)
 

 

 


Магнитное поле создаваемое этим бесконечно малым отрезком тока центральной части перемычки
 

dH_from_peremychka_center_ := Curl(dA_peremychka_center_) 

`#mover(mi( (55)
 


 


Бесконечно малый отрезок тока в центральной части перемычки создаёт магнитное поле в центральной точке  правой и левой рельсы
 

 

dH_from_peremychka_center_on_right_rail_center_ := subs(z = Z_relsotron, x = X_right, y = Y_center_relsy, dH_from_peremychka_center_); 1; dH_from_peremychka_center_on_left_rail_center_ := subs(z = Z_r...
dH_from_peremychka_center_on_right_rail_center_ := subs(z = Z_relsotron, x = X_right, y = Y_center_relsy, dH_from_peremychka_center_); 1; dH_from_peremychka_center_on_left_rail_center_ := subs(z = Z_r...
 

 

`#mover(mi(

`#mover(mi( (56)
 

 


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

 

df_lorenz_from_peremychka_center_to_right_rail_center_ := `*`(mu, `*`(Typesetting:-delayCrossProduct(`*`(I__r, `*`(`#mover(mi(
df_lorenz_from_peremychka_center_to_right_rail_center_ := `*`(mu, `*`(Typesetting:-delayCrossProduct(`*`(I__r, `*`(`#mover(mi(
 

 

`#mover(mi(

`#mover(mi( (57)
 


 

 

Скалярное магнитное поле Николаева, создаваемое этим бесконечно малым отрезком тока центральной части перемычки 

dS_from_peremychka_center_ := `+`(`-`(Divergence(dA_from_peremychka_center_))) 

`#mover(mi( (58)
 


 

Бесконечно малый отрезок тока в центральной части перемычки создаёт скалярное магнитное поле Николаева в центральной точке  правой и левой рельсы 

 

dS_from_peremychka_center_on_right_rail_center_ := subs(z = Z_relsotron, x = X_right_rail, y = Y_center_relsy, dS_from_peremychka_center_); 1; dS_from_peremychka_center_on_left_rail_center_ := subs(z ...
dS_from_peremychka_center_on_right_rail_center_ := subs(z = Z_relsotron, x = X_right_rail, y = Y_center_relsy, dS_from_peremychka_center_); 1; dS_from_peremychka_center_on_left_rail_center_ := subs(z ...
 

 

`#mover(mi(

`#mover(mi( (59)
 

 


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

 

df_nikolaev_from_peremychka_center_to_right_rail_center_ := `*`(mu, `*`(I__r, `*`(`#mover(mi(
df_nikolaev_from_peremychka_center_to_right_rail_center_ := `*`(mu, `*`(I__r, `*`(`#mover(mi(
 

 

`#mover(mi(

`#mover(mi( (60)
 


 

 

 

 

 


Векторный потенциал, создаваемый бесконечно малыми участками тока центральной части правой и левой рельсы
 

 

dA_from_center_of_right_rail_ := subs(x__i1 = X_right_rail, z__i1 = Z_relsotron, y__i1 = Y_center_relsy, dA_(x, y, z, t, I__r, x__i1, y__i1, z__i1, `#mover(mi( 

`#mover(mi( (61)
 

dA_from_center_of_left_rail_ := subs(x__i1 = X_left_rail, z__i1 = Z_relsotron, y__i1 = Y_center_relsy, dA_(x, y, z, t, I__r, x__i1, y__i1, z__i1, `+`(`-`(`#mover(mi( 

`#mover(mi( (62)
 


 


Магнитные поля создаваемые этими участками тока в центральной точке перемычки
 

 

dH_from_center_of_right_rail_on_peremychka_center_ := subs(x = X_center_peremychka, z = Z_relsotron, y = Y_peremychka, Curl(dA_from_center_of_right_rail_)); 1; dH_from_center_of_left_rail_on_peremychk...
dH_from_center_of_right_rail_on_peremychka_center_ := subs(x = X_center_peremychka, z = Z_relsotron, y = Y_peremychka, Curl(dA_from_center_of_right_rail_)); 1; dH_from_center_of_left_rail_on_peremychk...
dH_from_center_of_right_rail_on_peremychka_center_ := subs(x = X_center_peremychka, z = Z_relsotron, y = Y_peremychka, Curl(dA_from_center_of_right_rail_)); 1; dH_from_center_of_left_rail_on_peremychk...
dH_from_center_of_right_rail_on_peremychka_center_ := subs(x = X_center_peremychka, z = Z_relsotron, y = Y_peremychka, Curl(dA_from_center_of_right_rail_)); 1; dH_from_center_of_left_rail_on_peremychk...
 

 

`#mover(mi(

`#mover(mi( (63)
 


Скалярные магнитные поля Николаева создаваемые этими участками тока в центральной точке перемычки
 

 

dS_from_center_of_right_rail_on_peremychka_center_ := subs(x = X_center_peremychka, z = Z_relsotron, y = Y_peremychka, `+`(`-`(Divergence(dA_from_center_of_right_rail_)))); 1; dS_from_center_of_left_r...
dS_from_center_of_right_rail_on_peremychka_center_ := subs(x = X_center_peremychka, z = Z_relsotron, y = Y_peremychka, `+`(`-`(Divergence(dA_from_center_of_right_rail_)))); 1; dS_from_center_of_left_r...
dS_from_center_of_right_rail_on_peremychka_center_ := subs(x = X_center_peremychka, z = Z_relsotron, y = Y_peremychka, `+`(`-`(Divergence(dA_from_center_of_right_rail_)))); 1; dS_from_center_of_left_r...
dS_from_center_of_right_rail_on_peremychka_center_ := subs(x = X_center_peremychka, z = Z_relsotron, y = Y_peremychka, `+`(`-`(Divergence(dA_from_center_of_right_rail_)))); 1; dS_from_center_of_left_r...
 

 

`#mover(mi(

`#mover(mi( (64)
 


Применяя формулу силы Лоренца - находим таким образом выражение для силы, с которой Магнитные поля создаваемый бесконечно малыми участками тока центральной части правой и левой рельсы действует на бесконечно малый участок тока в точке центральной части перемычки
 

 

df_lorenz_from_center_of_right_rail_to_peremychka_center_ := `*`(mu, `*`(Typesetting:-delayCrossProduct(`+`(`-`(`*`(I__r, `*`(`#mover(mi(
df_lorenz_from_center_of_right_rail_to_peremychka_center_ := `*`(mu, `*`(Typesetting:-delayCrossProduct(`+`(`-`(`*`(I__r, `*`(`#mover(mi(
 

 

`#mover(mi(

`#mover(mi( (65)
 


 

 

 

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

 

df_nikolaev_from_center_of_right_rail_to_peremychka_center_ := `+`(`-`(`*`(mu, `*`(I__r, `*`(`#mover(mi(
df_nikolaev_from_center_of_right_rail_to_peremychka_center_ := `+`(`-`(`*`(mu, `*`(I__r, `*`(`#mover(mi(
 

 

`#mover(mi(

`#mover(mi( (66)
 


 

 

 


Сравним теперь сумму сил Лоренца и Николаева
 

Взаимодействие перемычки и правой рельсы: 

 

df_lorenz_from_peremychka_center_to_right_rail_center_ 

0 (67)
 

df_nikolaev_from_peremychka_center_to_right_rail_center_ 

`+`(`-`(`*`(mu, `*`(`^`(I__r, 2), `*`(`#mover(mi( (68)
 

df_lorenz_from_center_of_right_rail_to_peremychka_center_ 

`*`(mu, `*`(`^`(I__r, 2), `*`(sqrt(101), `*`(sqrt(4), `*`(`#mover(mi( (69)
 

 

df_nikolaev_from_center_of_right_rail_to_peremychka_center_ 

`+`(`-`(`*`(5, `*`(mu, `*`(`^`(I__r, 2), `*`(`#mover(mi( (70)
 

 

 

Взаимодействие перемычки и левой рельсы: 

 

 

df_lorenz_from_peremychka_center_to_left_rail_center_ 

0 (71)
 

df_nikolaev_from_peremychka_center_to_left_rail_center_ 

`+`(`-`(`*`(mu, `*`(`^`(I__r, 2), `*`(`#mover(mi( (72)
 

df_lorenz_from_center_of_left_rail_to_peremychka_center_ 

`*`(mu, `*`(`^`(I__r, 2), `*`(sqrt(101), `*`(sqrt(4), `*`(`#mover(mi( (73)
 

df_nikolaev_from_center_of_left_rail_to_peremychka_center_ 

`+`(`*`(5, `*`(mu, `*`(`^`(I__r, 2), `*`(`#mover(mi( (74)
 

 

 

 


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

 

Взаимодействие элементов тока 

 

 

Радиус вектор, направленный от элемента тока `#mover(mi( , расположенном в точке x__i1, y__i1, z__i1 к точке x, y, z 

 

R_ := proc (x, y, z, x__i1, y__i1, z__i1) options operator, arrow; `+`(`*`(`+`(x, `-`(x__i1)), `*`(_i)), `*`(`+`(y, `-`(y__i1)), `*`(_j)), `*`(`+`(z, `-`(z__i1)), `*`(_k))) end proc 

`#mover(mi( (75)
 

 

 

 

Векторный потенциал в точке x, y, z, создаваемый участком тока `#mover(mi( , расположенном в точке x__i1, y__i1, z__i1 

dA_ := proc (x, y, z, t, I__1, x__i1, y__i1, z__i1, `#mover(mi( 

`#mover(mi( (76)
 

 

 

 

Магнитное поле в точке x, y, z, создаваемое участком тока `#mover(mi( , расположенном в точке x__i1, y__i1, z__i1 

 

dH := simplify(Curl(dA_(x, y, z, t, I__1, x__i1, y__i1, z__i1, `#mover(mi( 

dH := `+`(`-`(`*`(I__1, `*`(Typesetting:-delayCrossProduct(`+`(`*`(`+`(x, `-`(x__i1)), `*`(`#mover(mi( (77)
 

 

 

 

 

Скалярное магнитное поле в точке x, y, z, создаваемое участком тока `#mover(mi(, расположенном в точке x__i1, y__i1, z__i1 

 

dS := simplify(`+`(`-`(Divergence(dA_(x, y, z, t, I__1, x__i1, y__i1, z__i1, `#mover(mi( 

dS := `*`(I__1, `*`(Typesetting:-delayDotProduct(`+`(`*`(`+`(x, `-`(x__i1)), `*`(`#mover(mi( (78)
 

 

Сила Лоренца, действующая на элемент тока `#mover(mi( , расположенном в точке x, y, z со стороны магнитного поля dH 

df_lorenz_ := simplify(`*`(mu, `*`(Typesetting:-delayCrossProduct(`*`(I__2, `*`(`#mover(mi( 

`#mover(mi( (79)
 

 

 

Сила Николаева, действующая на элемент тока `#mover(mi( , расположенном в точке x, y, z со стороны скалярного магнитного поля Николаева dS 

`*`(I__1, `*`(Typesetting:-delayDotProduct(`+`(`*`(`+`(x, `-`(x__i1)), `*`(`#mover(mi( (80)
 

df_nikolaev_ := `*`(mu, `*`(I__2, `*`(`#mover(mi( 

`#mover(mi( (81)
 

Суммарная сила, действующая со стороны элемента тока `#mover(mi( на элемент тока `#mover(mi( 

df_ := simplify(`+`(df_lorenz_, df_nikolaev_)) 

`#mover(mi( (82)
 

 

Приводим выражение суммарной силы к удобочитаемому векторному варианту 

 

 

`#mover(mi( 

`#mover(mi( (83)
 

Мы видим, что полученное выражение для выражение силы взаимодействия двух токовых элементов отличается от приведенной в книге Тамма [8], формулы (43.1)  

 

`*`(mu, `*`(I__1, `*`(I__2, `*`(Typesetting:-delayCrossProduct(`#mover(mi( 

 

введением дополнительного слагаемого `*`(`#mover(mi( 

 

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

 

 

Image 

 

 

Основываясь вышеприведенном рассуждении Тамма о неоднозначности определения силы взаимодействия токовых элементов, легко показать соответствие полученной в данной работе формулы формуле  (43.1) у Тамма. Действительно дивергенция векторного потенциала, согласно формуле 46.7 [Тамм, 8] вне поверхности всех обтекаемых током проводников равна нулю. Следовательно, интеграл по замкнутому токовому контуру для силы Николаева вне данного токового контура равен нулю. 

 

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

dF := int(`*`(mu, `*`(I__1, `*`(I__2, `*`(`+`(`*`(`#mover(mi( 

 

выражения для радиус вектора и векторов тока, установив в ноль координату z для обоих проводов, получаем 

 

assume(x__i2, 'real', y__i2, 'real', z__i2, 'real') 

 

 

dF := int(subs(`#mover(mi( 

dF := `/`(`*`(`#mover(mi( (84)
 

 

правильное значения силы взаимодействия параллельных проводников. 

 

Резюмируя данную работу, хотелось бы отметить, что взятая за основу формула векторного потенциала подвергается ревизии, например, в работе [9]. Следовательно, полученная в данной работе формула взаимодействия двух токовых элементов может быть пересмотрена вместе с пересмотром формулы для векторного потенциала. 

 

Литература. 

 

1. А.Ю.Дроздов. Электродинамический расчёт рельсотрона с использованием понятий векторного потенциала, магнитного поля и силы Лоренца. 02.04.2018. http://liquidcrystalosmos.narod.ru/railgunelectrodynamic.htm 

2.Научно-технический форум SciTecLibrary - Точные науки и дисциплины - Физика,  астрономия,  математические решения - Отдача в рельсотроне с позиций электродинамики http://www.sciteclibrary.ru/cgi-bin/yabb2/YaBB.pl?num=1522662090/new  

3.Форум движения за возрождение отечественной науки - Главная категория - Физика, техника, технологии - Отдача в рельсотроне  http://www.forum.za-nauku.ru/index.php/topic,4756.0.html  

4.Отдача в рельсотроне приложена к рельсам 2. https://youtu.be/bGBGk_iJtZQ 

5.Опыт Николаева номер 3 с продольным движением жидкого металла https://youtu.be/2oVKJBspu0E https://youtu.be/s3MfFzIBRHE 

6. Ф.Ф.Менде. Новая электродинамика, 2011 

7. Геннадий Николаев. Непротиворечивая Электродинамика. Теория, эксперименты, парадоксы. Книга 1. 1997 

8. И.Е.Тамм. Основы теории электричества. М.1957 

9. А.С. Чуев. О ВЕКТОРНОМ ПОТЕНЦИАЛЕ ОДИНОЧНО ДВИЖУЩЕГОСЯ ЭЛЕКТРИЧЕСКОГО ЗАРЯДА  http://www.sciteclibrary.ru/texsts/rus/stat/st3067.pdf 

 



02.04.2018 -19.04.2018
 

 

 

 

 

 

 

 

 

 

 

Обратная связь

Стена

 

© liquidcrystalosmos