Привет, ты когда-нибудь задумывался, как формирование кристалликов сахара в карамели почти то же самое, что и математика, объясняющая, как пыль в космосе собирается в звезду? Меня это безумно завораживает, особенно то, как температура влияет на этот процесс. Уверен, тебе это тоже очень интересно, но уже в космических масштабах. Давай обменяемся мыслями?

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

Я тут с этим возилась, как-то комбинированно: решаю идеальные уравнения МГД явно на сетке, но потом добавляю порог электропроводности, который увеличивается с локальной плотностью. Получается, как будто магнитное поле надеваю защитные очки для плотных областей. Ещё я использую полунеявный метод для жёстких членов диффузии, чтобы шаг по времени не превратился в какой-нибудь бесформенный, никогда не поднимающийся пирог. Это много бумажной работы, но если я вынесу связь между магнитным полем в отдельный модуль, который потом можно будет заменить на более сложную модель плазмы с заряженными пылинками, то хотя бы код не развалится от перегрузки.