LINUX.ORG.RU

Считаем плазму

 , , ,


0

3

На фоне в ssh крутится самописная (c++/openmp) программа для расчёта горения разряда в холловском двигателе в аксиально-азимутальном приближении. На графиках — нарисованная с помощью matplotlib визуализация градиентно-дрейфовой неустойчивости. Эта неустойчивость возникает при сравнительно небольших концентрациях плазмы и в отсутствие столкновений. В результате, азимутальная симметрия разряда нарушается, а электроны формируют вихри и ненулевой ток в сторону анода. В целом, это похоже на формирование сгустков в магнетронах. При большой концентрации плазмы (и наличии столкновений) вихри разрушаются, но возбуждаются более короткие ионно-звуковые волны.

В остальном, обычный Debian8.1 с Xfce4; на числодробилке — ScientificLinux6.1.

>>> Просмотр (1920x1080, 576 Kb)

★★★★★

Проверено: JB ()
Последнее исправление: thunar (всего исправлений: 3)

Ответ на: комментарий от Deleted

чего с icewm на крысу переехал?

xfce4-clipman & xfce4-notes

thunar ★★★★★
() автор топика

В результате, азимутальная симметрия разряда нарушается, а электроны формируют вихри и ненулевой ток в сторону анода.

Я так и знал!

enso ★★★★★
()
Ответ на: комментарий от BceM_IIpuBeT

Это Х окна по ssh?

Нет. Монтирую через sshfs и работаю как с локальными файлами.

ШГ такое, шо ну нафиг.

Вот как ни странно, что на домашнем, что на рабочем компьютере они выглядят отлично. А поглядел как-то с компьютера брата — глаза заболели. Видимо, от монитора всё зависит.

thunar ★★★★★
() автор топика
Последнее исправление: thunar (всего исправлений: 2)
Ответ на: комментарий от BceM_IIpuBeT

К сожалению, крыса не умеет отдельный шрифт на панели. А что с общим видом не так-то?

thunar ★★★★★
() автор топика
Последнее исправление: thunar (всего исправлений: 1)

Считаем плазму

ту, которая KDE?

w1nner ★★★★★
()

О, а как реализовано несколько окон с матплотлибом? Просто разные скрипты запускаются, или это внутри одного несколько окон рисуются?

И кстати, как сделать цветную шкалу такую? А то стандартная rainbow рисует с фиолетовым.

knovich
()
Последнее исправление: knovich (всего исправлений: 1)

В остальном, обычный Debian8.1

Почему не Slackware?

densss ★★★★
()
Ответ на: комментарий от knovich

как реализовано несколько окон с матплотлибом

fig1 = plt.figure(); fig2 = plt.figure() же

цветную шкалу

https://matplotlib.org/examples/color/colormaps_reference.html (jet)

thunar ★★★★★
() автор топика
Последнее исправление: thunar (всего исправлений: 2)
Ответ на: комментарий от thunar

А, jet я и проглядел, спасибо!

Что касается нескольких окон, то, как я понимаю, такое решение блокирует весь скрипт. Просто я пишу некоторую программу, в которой пользователь (я) выбирает произвольные данные, строит график, потом выбирает другие данные, тоже строит, не закрывая предыдущую картинку. А если как у вас, то графики не отрисуются без plt.show, который, в свою очередь, не даст вернуться в mainloop интерфейса, пока не закроешь окошки с графиками. Или я неправ? Мне приходится делать через multiprocessing.

knovich
()
Ответ на: комментарий от thunar

Так я заранее не знаю, сколько будет окон с графиками.

knovich
()

За Redmond-тему (Win95) в заголовках окон — плюс.

EXL ★★★★★
()

Не сказать что ШГ, но панель предельно вырвиглазна из-за отступов, которых нет.

bread
()
Ответ на: комментарий от knovich

Не, мне интерактивность не нужна же. Там один полный кадр считается несколько минут (а иной раз и по паре часов).

thunar ★★★★★
() автор топика
Ответ на: комментарий от thunar

jet

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

i-rinat ★★★★★
()

Теперь понятно, почему плазма падает.

KOT040188 ★★
()
Ответ на: комментарий от thunar

Для меня необходимость перерисовывать картинки для печати в ч/б была неприятным сюрпризом. Если только смотреть на экране, то не так страшно. Но другие палитры всё равно попробовать стоит. Особенно те равномерные. Они не зря в списке не той страничке первыми указаны.

i-rinat ★★★★★
()
Ответ на: комментарий от i-rinat

Вот сколько встречал статей с ЧБ рисунками, всегда градиенты там были совершенно нечитаемыми.

thunar ★★★★★
() автор топика
Последнее исправление: thunar (всего исправлений: 1)
Ответ на: комментарий от thunar

Ну так мало кто парится. Берут и конвертируют готовую картинку в ч/б. По дефолту много где стоит похожая на jet палитра, вот и получается жуть.

А ещё стоит делать изолинии вместо градиентов, чтобы были видны мелкие детали.

i-rinat ★★★★★
()
Последнее исправление: i-rinat (всего исправлений: 1)

Круто!

Интересно было бы посчитать плазму в лазерном двигателе, перспективная имхо штука (есть ссылочка на недосчитанную работу, которая хорошо отработала в эксперименте)

michwill ★★★★★
()

У тебя оно 2Д принципиально?

Если планируешь переходить в 3Д, то там и «гляделки» уже есть широко используемые, типа VisIt или ParaView или Vapor (кажется).

Это оно как двухжидкостная МГД считается?

sshestov
()
Ответ на: У тебя оно 2Д принципиально? от sshestov

Это оно как двухжидкостная МГД считается?

PiC/MCC по всем компонентам. МГД там из серии «на гранцие применимости / не работает».

Если планируешь переходить в 3Д, то там и «гляделки» уже есть широко используемые, типа VisIt или ParaView или Vapor (кажется).

В ближайших проектах как раз перейти в 3d3v-геометрию (и перенести код на CUDA).

Ну да, действительно. Надо будет освоить формат смотрелок что бы велосипед не писать.

thunar ★★★★★
() автор топика
Последнее исправление: thunar (всего исправлений: 4)
Ответ на: комментарий от t184256

говорю ж

это VisIt и ParaView. Специфические программы для просмотра а) больших; б) 3D данных разнообразных форматов. Спроектированные и сделанные с учетом этого, типа клиент-сервер и многопроцессорность изкоробки (искаропки? как там правильно?), поддержка разных хитрых форматов типа vtk или HDF5.

А чего такое MCC? Сходу не соображу аббревиатуру?

sshestov
()
Ответ на: говорю ж от sshestov

Ну вот есть у меня многогиговые hdf5-файлы с compound datatypes, как их отображать с помощью этого Paraview? И что я с этого выиграю? Как-то непонятно.

MCC

Google говорит, что Monte-Carlo Collisions. У меня такой фигни нет.

t184256 ★★★★★
()

Годная работа. А при каких температурах (эВ) это у тебя там крутится? Какие взаимодействия частиц просчитываешь? Диаграмму направленности переизлучения сам рассчитываешь?

pacify ★★★★★
()
Ответ на: комментарий от t184256

так ты чего именно спрашиваешь?

ни там, ни там - интерфейс не очевидный. Глянь help какой-нибудь. В ParaView нужно сделать (кажется) pipeline, куда сверху пихаешь данные, потом на него фильтры для вывода конкретной величины. Он чуть более заморочен чем визит.

VisIt в этом смысле чуть попроще: данные загружаешь через меню file, потом в основном окошке делаешь какой-либо тип вывода через «add», например pseudocolor. Далее можно ему сечение прорисовать - это Operators->Slicing->Slice.

Многогиговые - это сколько? У меня при 200Gb снапшоте лаги терпимые; при 400Gb уже хуже. Причем visit в бинарниках скомпилирован без mpi. Прересборка с ним ситуацию улучшила, особенно если файл умеет читаться в mpi-режиме.

sshestov
()
Ответ на: комментарий от pacify

А при каких температурах (эВ) это у тебя там крутится?

Там нет термодинамического равновесия, так что, если под «температурой» понимать среднюю энергию, то ~10эВ в спокойных областях, и до ~300эВ в ускорительном слое (разброс большой).

Какие взаимодействия частиц просчитываешь?

Электронов с фоном: кулоновское рассеяние на ионах, упругие и неупругие рассеяния на нейтралах.

Диаграмму направленности переизлучения сам рассчитываешь?

Слишком медленные процессы — безизлучательное приближение: ∂ₜE⟶0 и ∂ₜB⟶0.

thunar ★★★★★
() автор топика
Последнее исправление: thunar (всего исправлений: 2)
Ответ на: комментарий от Odalist

На десктопе надёжность важнее мифической сэкурности.

thunar ★★★★★
() автор топика
Ответ на: так ты чего именно спрашиваешь? от sshestov

так ты чего именно спрашиваешь?

Что я теряю в сравнении с просмотром через самописные велосипеды c matplotlib.

Видимо ответ — с этими штуками не надо писать код.

Многогиговые - это сколько?

Не больше 10.

t184256 ★★★★★
()
Последнее исправление: t184256 (всего исправлений: 1)
Ответ на: комментарий от thunar

Внутри debian'овского релиза это не вовсе не принято. Лень обновить минорщинку внутри oldstable? Не надо отмазок, это точно лень.

t184256 ★★★★★
()
Ответ на: комментарий от t184256

не теряешь вроде особенно ничего

но когда я пытался ковырять в IDL относительно большие данные c AMR сеткой, то мне показалось это неудобно. Да и не быстро. А там сразу какая угодно визуализация и преобразования, которые почти что из коробки могут быть в mpi.

Я в последнем своем посте в попугаях запутался. У меня снапшоты в мегабайтах, конечно же. При Gb (такие, впрочем, тоже есть) оно даже считывается долго. Поэтому, конечно, удобнее пускать в клиент-сервер режиме, чтобы рабочий комп не лагал.

sshestov
()
Ответ на: комментарий от t184256

Внутри debian'овского релиза это не вовсе не принято. Лень обновить минорщинку внутри oldstable? Не надо отмазок, это точно лень.

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

Odalist ★★★★★
()
Ответ на: комментарий от Odalist

Что делает? Если в ведре или фуррифоксе обнаружиат уязвимость, то она мне прикатит в уведомлениях в секурити-апдейтах. Но это совершенно не повод всю остальную систему потенциально ломать и обновлять софт всем скопом.

thunar ★★★★★
() автор топика
Ответ на: не теряешь вроде особенно ничего от sshestov

Ну маленькие файлы визуализировать я могу и коротким сочинением на numpy/matplotlib в Jupyter. А про большие надо покрутить/подумать, вся ли обработка у меня упирается в IO.

t184256 ★★★★★
()
Ответ на: комментарий от thunar

Но это совершенно не повод всю остальную систему потенциально ломать и обновлять софт всем скопом.

Вы плохо ознакомлены с внутри debian'овскими релизами....это вам не Arch, или какая то Федорка, где одним чихом все ломается. Ну да ладно. Это ваше личное дело.

Odalist ★★★★★
()
Вы не можете добавлять комментарии в эту тему. Тема перемещена в архив.