Распределённая МДС: источник коэффициентов 1/2 и 1/3

Методика: О. Б. Буль, 2005, печатные с.18–21. Геометрия и размеры
эталона созданы для курса. Оригинальные модели и страницы не включены.

1. Сохраните reproduce.py в отдельную папку и выполните Run Python Script
   в FEMM 5. Результаты каждого повтора появятся в новом calculated-NN.
2. Откройте FEM базового режима. N=300, I=1 А, ширина g=50 мм,
   высота h=40 мм, Depth=20 мм. Вся область немагнитная, mu_r=1.
   Группа 1 — распределённая обмотка. Имена цепи и свойств русские.
3. Верхняя граница A=0 подразумевает обратный поверхностный ток.
   На трёх остальных сторонах dA/dn=0. Это локальный предельный слой,
   не самостоятельный полный электромагнит с внешним воздухом.
4. При равномерном токе Hx(y)=-NI*y/(g*h), Bx=mu0*Hx.
   Геометрическая проводимость Lambda_g=mu0*b*h/g.
   Поток через слой: Phi=Lambda_g*NI/2.
   Потокосцепление: Psi=Lambda_g*N^2*I/3.
   Энергия: W=I*Psi/2. Все длины в этих формулах в метрах.
5. Сопоставьте три сетки 1;0.5;0.25 мм. Поток определяется через
   b*Delta A; точки взяты на 0.001 мм внутрь от горизонтальных границ,
   поэтому даже при идеальной сетке коэффициент немного меньше 1/2.
   Контролируйте эту известную погрешность пути отдельно от сеточной.
6. При I=2 А Psi удваивается, W увеличивается в четыре раза. Результат
   относится к тому же линейному слою. Разные распределения витков
   требуют новых интегралов, коэффициенты 1/2 и 1/3 не универсальны.

Проверенный режим на сетке 0.25 мм:
72559 элементов; k_Phi=0.4999754844; k_Psi=0.3333299566;
Psi=0.000603179679104 Вб; W=0.000301589839552 Дж.
Четыре заново сохранённых и решённых поля, 31 численная проверка.
verified/result.json хранит числа, publication.json — SHA-256 точных файлов.