Re: Выставлю на обсуждение работу akayomov_andrew.pioneer14May 5 2013, 17:20:47 UTC
Объясните, пожалуйста, какое отношение этот пример имеет к гравитации? В коде считаются только столкновения частиц, у которых есть какая-то начальная скорость. Ни притяжения, ни отталкивания там не моделируется.
Re: Выставлю на обсуждение работу akayomov_andrew.akayomov_andrewMay 6 2013, 22:48:05 UTC
Вот это то самое главное! Только сталкивание частиц и ничего больше. Но Большие шары сближаются! и с ускорением. То есть никаких гравитационных сил для объяснения сближения шаров не нужно. Вполне достаточно 2 закона Ньютона действие равно противодействию. 3 не надо 3 выводится.
Гравитация - это поток гравитационных частиц, который исходит из ядра, а потом возвращается. Гравитация действует механическим образом, избегая (других видов) взаимодействий со всеми, кроме некоторых субатомных частиц, так что будучи независимыми, гравитационные частицы движутся нисходящим потоком и восходящей струёй, сталкивая встречающиеся на своём пути составные части атомов. Взаимодействие с другими частицами является, по существу, механическим, потому что поток гравитационных частиц имеет очень большую скорость, и поэтому множества взаимодействий с потоком, который выходит из гравитационного гиганта и возвращается обратно в этот гигант, не случается. Причина того, почему восходящий поток более быстр, - не только в сжатии, когда давление потока находит выход, но и в образующейся траектории. Таким образом, на выходе скорость более высокая, в то время как на входе механическое давление распространяется вниз на более обширную площадь, и это требует большего времени
( ... )
Спорить на эту тему не буду мне это не интересно. вот код писать интересно а спорить нет. Да и смысла не вижу. Давайте распишите все подробнее и попробуем вместе сделать модель. Если заработает вы правы если нет...ну нет так нет... У меня модель на идею Виктора http://viictor.livejournal.com/228424.html#t4057416 Там очень мало допущений поэтому она меня и заинтересовала
Re: линк на скачивание Visual и Pythonse_frolovMay 3 2013, 06:51:21 UTC
Завтра посмотрю на ваш код подробнее, но пока масса проблем с кодом и методологией.
(а) Traceback (most recent call last): File "vii.py", line 110, in x,y,z,vv,r =xyzclist[i] IndexError: list index out of range
Цикл по Natoms, но размер массива xyzclist меньше => buggy code.
(b) Вывести колонку чисел, изменяющихся с 0.24 до 0.23 как-то маловато. Нужно полноценное статистическое исследование результатов. Кроме того, "теория" Виктора утверждает, что результирующее движение тел будет совпадать со следствием привычного гравитационного взаимодействия. Так что нужно проводить более внимательный анализ.
(c) Наконец, выбраны отражающие граничные условия. Почему? Как это скажется на результате? А если стенки будут в 2 раза дальше?
Re: линк на скачивание Visual и Pythonakayomov_andrewMay 6 2013, 22:57:46 UTC
Вы правы конечно. Но здесь речь идет только о тенденции. Теорию Виктора проверить на привычное гравитационное взаимодействие на компе нереально. Нужно напустить в модель слишком много мелких шариков. ну 10 000 000 000 не меньше. А мой уже при 4000 дает ошибку памяти. А багу уже исправил и здесь об этом сообщил и исправленный код выложил.
либо текст программы не актуальный, либо одно из двух ... но кроме статичных картинок с шариками и кривой ничего программа не даёт. интерпретатор выдаёт ошибки индексации и создания массивов в циклах, строки с 100 по 150. строчки с проблемами отмечены ####### дальше ничего не изучал, т.к. не знаю сути алгоритма.
Re: Ошибка в 87 строкеvakomiMay 3 2013, 19:58:52 UTC
благодарю! вроде как заработало. правда есть некоторое различие в картинках по сравнению с вашими скриншотами.
облако гораздо мельче виден куб целиком. второй экран - гистограмма а не такой график, как на скриншоте. Он у вас в конце работы программы формируется? гистограмма за верхний край вылезает. вы не подскажите где какие параметры подвигать?
большое спасибо. я благодаря вам до питона добрался. а то всё не досуг было. :)
Re: Ошибка в 87 строкеakayomov_andrewMay 6 2013, 23:01:23 UTC
а он мышкой крутится. Кнопку нажмите и двигайте. Ну потом все начальные значения случайны и все время разные. Только большие шары задаются. (они кстати в начале неподвижны и импульса не имеют)
2 комментария: 1. Нужно исследовать статистику для многих запусков. 2. Нужно продемонстрировать, что dt выбран достаточно малым. Я попробовал уменьшить dt вдвое и флуктуации результата стали заметно меньше.
Да так и есть. А еще лучше частиц добавить и их размеры уменьшить ну и отдалить стенки чтоб шары не притягивали. Я собственно для того и разместил код чтоб с ним можно было "поиграться" и найти наилучшие начальные условия.
Comments 20
Reply
Reply
В коде считаются только столкновения частиц, у которых есть какая-то начальная скорость.
Ни притяжения, ни отталкивания там не моделируется.
Reply
То есть никаких гравитационных сил для объяснения сближения шаров не нужно. Вполне достаточно 2 закона Ньютона действие равно противодействию. 3 не надо 3 выводится.
Reply
Reply
Да и смысла не вижу.
Давайте распишите все подробнее и попробуем вместе сделать модель.
Если заработает вы правы если нет...ну нет так нет...
У меня модель на идею Виктора
http://viictor.livejournal.com/228424.html#t4057416
Там очень мало допущений поэтому она меня и заинтересовала
Reply
Какая именно библиотека - там есть варианты
Reply
это Python
http://sourceforge.net/projects/vpythonwx/reviews/?source=navbar
это библиотека визуализации ее надо скачать после
(моя) программа
Передел вот эту.
Python27\Lib\site-packages\visual\examples\gas.py
Вообще-то базар пошел отсюда
http://viictor.livejournal.com/228424.html#t4057160
предлагаю здесь и продолжить. Впрочем не настаиваю это как вам удобно.
Reply
(а) Traceback (most recent call last):
File "vii.py", line 110, in
x,y,z,vv,r =xyzclist[i]
IndexError: list index out of range
Цикл по Natoms, но размер массива xyzclist меньше => buggy code.
(b) Вывести колонку чисел, изменяющихся с 0.24 до 0.23 как-то маловато. Нужно полноценное статистическое исследование результатов. Кроме того, "теория" Виктора утверждает, что результирующее движение тел будет совпадать со следствием привычного гравитационного взаимодействия. Так что нужно проводить более внимательный анализ.
(c) Наконец, выбраны отражающие граничные условия. Почему? Как это скажется на результате? А если стенки будут в 2 раза дальше?
В общем, вопросов и проблем пока много.
Reply
Reply
но кроме статичных картинок с шариками и кривой ничего программа не даёт.
интерпретатор выдаёт ошибки индексации и создания массивов в циклах, строки с 100 по 150.
строчки с проблемами отмечены #######
дальше ничего не изучал, т.к. не знаю сути алгоритма.
for i in range(Natoms):
x,y,z,vv,r =xyzclist[i] #############################################
Atoms = Atoms+[sphere(pos=(x,y,z), radius=r, color=vv)]
mass = Matom*r**3/Ratom**3
pavg = sqrt(2.*mass*1.5*k*T) # average kinetic energy p**2/(2mass) = (3/2)kT
if i
Reply
while na
Reply
вроде как заработало.
правда есть некоторое различие в картинках по сравнению с вашими скриншотами.
облако гораздо мельче виден куб целиком.
второй экран - гистограмма а не такой график, как на скриншоте. Он у вас в конце работы программы формируется?
гистограмма за верхний край вылезает. вы не подскажите где какие параметры подвигать?
большое спасибо. я благодаря вам до питона добрался. а то всё не досуг было. :)
Reply
Только большие шары задаются. (они кстати в начале неподвижны и импульса не имеют)
Reply
1. Нужно исследовать статистику для многих запусков.
2. Нужно продемонстрировать, что dt выбран достаточно малым. Я попробовал уменьшить dt вдвое и флуктуации результата стали заметно меньше.
Reply
Я собственно для того и разместил код чтоб с ним можно было "поиграться" и найти наилучшие начальные условия.
Reply
Leave a comment