: Мда ...
: Если эффект действительно имеет место в том числе и для большого числа частиц, то это означает фундаментальные проблемы у статистической физики.
Да, ну! Я на такое не замахиваюсь. Это просто маленькое хакерство: поиск проблем.
: Вы не пробовали варьировать число частиц и смотреть за величиной эффекта? Может быть эффект будет при этом уменьшаться. А, следовательно, есть надежда, что при достаточно большом числе частиц - полностью исчезнет (предсказания статистической физики справедливы лишь для большого числа частиц).
Надежда конечно есть. Но в поле 1-D надежд маловато IMHO. А вот чем большая размерность пространства, тем плотнее центральные области.
: Вы хорошо протестировали программу? А то может где-то ошибка?
Так вот, обращаюсь:
: Предлагаю протестировать на распределение Максвелла. Оно должно устанавливаться после нескольких времен столкновения (имеется в виду в отсутствии гравитационного поля) и независимо от начального распределения. Если же оно не устанавливается, то либо имеет место ошибка в программе, либо взято слишком малое число частиц, чтобы можно было сравнивать с предсказанием статистической физики.
Кто спорит? Конечно, малое число частиц. Но кроме этого есть еще логика. Мое дело показать проблему, заострить внимание.
: Для скольких частиц Вы пробовали?
Ну, ясное дело, не для мильена. :)
Для большого числа частиц требуется «Эксперимент» (BvaguMup)
: Только не понимаю как же вы до центра то добрались? Ведь там - сингулярность. Или же Вы этими словами лишь озвучили тенденцию увеличения плотности на поверхности внутреннего ящика при уменьшении его радиуса?
Не понял? Сингулярность я использую как термояд в отсутствии 'внутреннего ящика'.
: Скажите, а эффект зависит от радиуса внутреннего ящика. Именно - не исчезает ли он при достаточно большой его величине? Если бы не исчезал, то модель может быть бы была адекватна просто атмосфере планеты.
Нет сингулярности, нет нагрева - нет и эффекта. Модель именно адекватна атмосфере или, как в сверхразреженном газе, адекватна планетарной модели.
: Тяжко мне это (разрываюсь).
: А ситуация здесь "проще простой" - основное из того, что нужно сделать - Вами уже сделано: ящик в ящике и частицы между ними. Все что теперь нужно - это дополнительно к численному решению уравнений движения всех частиц - параллельно решать еще одно линейное уравнение, которое будет зависеть от единственной функции - некоторой функции координат всех частиц (А именно - от суммы потенциалов взаимодействия всех частиц). В свою очередь это решение никак не входит в уравнения движения частиц. Так что вроде довольно просто. Жаль, что я не программирую на Дельфи (да уже и СИ забыл порядком). Впрочем нельзя ли как-нибудь Вашу программу просто перевести на мэйпл (в нем я еще - более мнение).
Что такое 'мейпл'? Что-то я не слышал о таком. А на Дельфи научиться - проще простого. Лучше, конечно, если модель будете делать Вы сами. Короче, устанавливайте Дельфи; заодно проверите и мою программу. |