_Vova: Огромное спасибо за ссылку с соленоидом!!! Я как-то повидимому мимо прошел и не обратил внимания, так как та ссылка с витком в моей первом сообщении была из той же конторы и той же оперы....
С соленоидом формулы проще, возможно осилю аналитически последний набор по радиусу, тогда получится то, что надо, если нет, то буду как есть считать, формула в общем-то очень быстро вычисляется и с адаптивным интегрированием должна давать хорошие результаты.
Цитата(_pv @ Sep 1 2017, 15:50)
ну и не нужен тут fem, всмысле разбиения всего простанства на элементы и решение граничных условий.
та же радия граничными интегралами считает
http://accelconf.web.cern.ch/accelconf/pac...s/pdf/9P027.PDFдля таких задач тоже вроде как быстрее должно быть.
да, бем - вещь хорошая, я за свое время несколько бем солверов, в том числе один максвеллосвкий написал. Есть большое подозрение, что конкретно для той геометрии, что у меня, он слостно продувает в производительности на точность, так как я аппроксимирую rot-rot в максвелле тензорными матрицами, для которых комплексити получается линейная, и даже использование доморощенно вылизанного пакета арш-матриц (H-matrix) не позволит бему выйти по скорости вперед.
По сути у меня куча катшек, мне надо их правильно расставить. Катушки могут иметь разные токи, даже (тут еще не обсуждалось) градиент тока. Для каждой такой задачи я совместно с задачей теплопереноса решаю, что у меня получается, и ищу оптимум. В этой задаче получается около 50 оптимизационных параметров и оптимизирующая функция обладает жудким числом локальных минимумов. Пробный старт позволил с приемлимой точностью считать за одну минуту примерно 20 различных наборов катушек на моем 24 терафлопнике (примерно 5 простых катушек в секунду), но сходимость будет примерно на миллионе наборов, что в общем-то не сильно радует, ибо около полутора месяцев. А вдруг мне еще что-то пересчитать захочется?
_pv: могу войсом рассказать зачем и почему все так, вы же мои контакты знаете