LINUX.ORG.RU

Linux как инструмент для решения задач моделирования


0

0

В продолжение разговора о разгерметизации салона самолета на высоте 10.000 м

Задачка решалась на маленьком 2-х нодном кластере примерно в течение 2-х часов.

>>> Просмотр (1024x768, 194 Kb)

★★★★★

Проверено: maxcom

А что это за Linux и что такое FreeCFD?

Rolex ★★
()

Моделирование ...

Это сделано на основе Slackware 8.0/8.1 

Кластер сделан на MOSIX версии 1.8.0 (ядро 2.4.19) 

FreeCFD это одна из _моих_ CFD софтинок - используется для моделирования
сжимаемых потоков - параллельная версия использует PVM.  

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

Моделирование ...

Вдогонку: со временем решения задачки в комментарии я чуток прогнал ...
время решения задачи составило чуть больше одного часа
Узлы кластера: Athlon 1000Mhz TB /0.5Gb и Athlon 1700+ XP /O.5Gb
Среда передачи: обыкновенная 100-мегабитка

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

хто на рабочих машинах ставит КДЕ? ствят Виндовмакир

кде - дома и в офиссе

anonymous
()

Можно только позавидовать :)

Интереснае задачи на интересном железе, что еще нужно для счастливой жизни :)

А что, софт и правда сам писал?

Dead ★★★★
()

Как сенсоры прикрутил? А то у меня, блин, ничего с ними не выходит 8( sensors-detect находит сенсоры и предлагает вставить модуля, которых у меня нет и в ядре нет (я все i2c скомпилил - ну нет там такого модуля)

Zulu ★★☆☆
()

А чего у молодого человека шрифты такие отстойные?

Да, и что за widget set был использован для FreeCFD? То, что Мотиф - это я вижу, а ещё что? Просто там tabs какие-то непривычные.

ДВ

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

Моделирование ...

>хто на рабочих машинах ставит КДЕ?

_Я_ ставлю - по юзабилити (ImHO) пока ничего лучше не придумано.

>ствят 
>Виндовмакир 

лично _МНЕ_ сей WM не нравится категорически ...

PS: ImHO у нас такие "советы" называются "пионерская распальцовка"
потому как выбор персонального десктопа дело абсолютно личное
потому как на нем _МНЕ_ работать а не _ВАМ_



PS2: ничего личного ;)

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

Моделирование ...

2Dead: Да сам. Правда софту уже много лет. Отдельные куски кода писаны еще годике эдак в 94-95. 2Zulu: lm_sensors 2.6.3 - поставленные согласно инструкции тонкость существует в порядке наложения патчей (MOSIX в этом отношении довольно капризен) рецепт такой - сначала на чистом 2.4.19 генерируется lm_sensors.diff - затем ставится MOSIX а потом прикладывается этот самый lm_sensors - остальные патчи прикладываются после этих 2-х в данном случае это lv_cool и на рабочей машинке alsa (шоб нескучно работалось ;) ) 2ДВ: Ну вопервых я не так уж и молод :) Ну а во вторых (и третьих) - мне некогда заниматься декорацией - мне нравятся именно _такие_ шрифты как максимально юзабельные на моем десктопе Насчет Widget-ов - это не Motif (просто я по старинке стараюсь делать внешний вид Motif-like - так привычнее.) а FLTK - в данном случае древней версии 1.0.6 к сожалению в более поздних поломали кучу полезных фишек типа fl_positioner и тех же самых Fl_Tabs :( Так что сейчас я нахожусь в стадии перехода на FOX Toolkit :(

sS ★★★★★
() автор топика
Ответ на: Vopros от SolarWind

Моделирование ...

2SolarWind:

С самолетом скорее всего ничего...

По предварительной прикидке на данных взятых с потолка (кроме разве
что параметров атмосферы на высоте 10 км и скорости в 900 км/ч) будет
следующее:

Наименее приятное место это возле разрушенного илюминатора - за пару
миллисекунд там давление падает до 0.7 атмосферного и скорость газа поперек салона возле ближайших кресел будет порядка 120 м/сек (ноут с коленок унесет запросто :)) вдоль салона в проходе между креслами
"ветерок" будет около 50-60 м/сек. Ближайшие к илюминатору кресла
действительно могут быть отломаны потоком воздуха - перепад на
них составит порядка 0.1 атм что при размере спинки кресла 1x0.6 метра создаст усилие в 600 кг

Правда точность сего решение +/- полкилометра потому как прикидка
делалась в 2D постановке (облом было рисовать салон в 3D да и времени
на расчет ушло бы на порядок больше :)) и не учитывается падение
давления в салоне - хотя и время процесса не слишком велико ...

Если интересны более точные данные - цену можно обсудить (начиная с 10 килобаксов) :)



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

А какие методы расчета использует эта самая FreeCFD? Конечных элементов? Или что-то свое?

S@nek...

anonymous
()

Насчёт шрифтов прошу пардону за замечание, это я сослепу в первый раз не разглядел, что там, на самом деле, всё нормально. Так что вопрос снят :)

Однако появились ещё два вопроса:

1) Что это за виджеты такие, которые векторное поле и графики внизу рисуют? Сам я пользуюсь библиотечками от Sitraka (бывшая KL Group) и очень ими доволен; но дело не в этом, а в том, что охота знать, чем другие пользуются.

2) Опять-таки прошу пардону за прямоту - но SlickEdit какой версии и где брал, если не секрет?

ДВ

DmVo
()

2 sS (*) (2002-11-04 12:16:49.001), sS (*) (2002-11-05 17:57:02.821)

Блин. Ну вы что, сомневались в моей правоте? Я же своими глазами кадры документальных съемок видел. :) Какой-то (немецкий?) науч-поп фильм. :)
Жалко ноут. И коленки. И кресло. :)
Завидую ученым - интересным делом заняты...

Успехов!

Виндузятник

anonymous
()

А можно вопрос к автору софта -почему был выбран MOSIX? Как был организован обмен нитками между считающими процами.
И последнее - насколько тормознее было-бы считать на двух-головом атлоне? И на сколько дешевле?
Вопрос задаю не из праздного интереса... Возможно это поможет чуток продвинуть нашу советскую науку. Надеюсь на ответ...

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

Моделирование ...


2S@nek...: Конечные элементы это скорее из solid dynamic хотя
есть варианты применения их и в CFD...данная софтина использует обыкновенный сеточный метод (в данном случае на регулярной сетке) установления (time dependence) схема a-la Мак-Кормак но со своими фишками (адаптивная численная вязкость).

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

Моделирование ...

2DmVos: Векторное поле сделано на основе fl_positioner (от него только
красное перекрестье) а остальное програмилось вручную через функции
FLTK - это довольно просто сам код занял около 1000 строк, это не только вектора но и поля распределений параметров (5 видов отображения) и изолинии. VSE 6 версии (есть уже 7 но я пока использую старый - мне его хватает) брать триал-версию можно у производителя на http://www.slickedit.com - а потом пропатчить ее до нормальной и найти в Inet-е keygen ну или выложить $$ за 7 версию :)
Если уж _ОЧЕНЬ_ нужно давай мыло - пришлю ссылку на комплект (для 6 версии)

2anonymous: Выбран MOSIX по следующим причинам:
- нужен был полноценный batch pool (не все мои задачи расчитаны на параллельное исполнение) хотя есть такие штуки как Condor но ... так сложилось что пошел MOSIX
- MOSIX+PVM гораздо проще и эффектиыней реализует балансировку в гетерогенном кластере чем чистый PVM (потому как СРАЗУ поиметь $$ на N-ное количество одинаковых узлов "простому советскому ученому" проблематично, то есть мощности наращиваются постепенно) а с MPI
получается больше возни.
- MOSIX наиболее (ImHO) динамично развивающаяся кластерная софтина для
Linux

Когда только плантровалась кластеризация в природе не существовало
не то что 2-х головых Athlon-ов а Athlon-ов вообще - в смысле они
были еще "нераспробованы" - первый кластер был построен на 2-х целиках
(433 и 500) - сейчас рассматиривается замена части кластера (всего он 8 узловой - просто в данном случае используется только 2 узла) в том числе на 2-х головые Tyan-ы. Наверное на 2-х головом считалось бы
быстрее. Кстати MOSIX не умеет (пока) никакие виды shared memory
так что задачи использующие pthread в нем параллелится не будут.
Данная софтина сделана хитро. Сама рисовалка малтитредовая,
а вот solver это пачка PVM-овских slave-ов которые могут свободно мигрировать на другие узлы (данные синхронизируются через pipe-ы)


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

А можно поподробнее о скрещивании PVM с MOSIX?

Deimos
()

Моделирование ...

>А можно поподробнее о скрещивании PVM с MOSIX?
>Deimos

Это проще чем кажется. MOSIX для PVM выглядит так же как SMP.
То есть особых телодвижений не нужно, все коммуникации идут
по схеме process2process минуя TCP транспорт pvmd.

sS ★★★★★
() автор топика
Ответ на: Моделирование ... от sS

sS (*) (2002-11-05 17:57:02.821):
>... и скорости в 900 км/ч
А на что скорость - то влияет? Типа, высоту эффективно увеличивает?

Die-Hard ★★★★★
()

2sS:
Кстати, а не мог бы чуть подробнее про метод? Просто в 2-х словах, для общего развития.
Я с CFD, в общем, не знаком, но с гиперболическими системами сталкивался достаточно
плотно, что такое предиктор-корректор знаю.

Die-Hard ★★★★★
()

2 sS

Не могли бы Вы ответить так на вскидку, ориентировочно, на основе Вашего опыта, сколько в процентном отношении по времени относительно указанного Вами длился бы расчет:

- отдельно на одном Athlon 1000Mhz TB - отдельно на одном Athlon 1700+ XP - при замене в Вашей схеме 100 Мбит на 1 Гбит - на 8-ми нодном 100 Мбит кластере, состоящим из одинаковых Athlon1000 - при добавлении к Вашей схеме одного старого P166MMX

Эти ответы были бы интересны для оценки потенциала кластерных вычислений. Допустим есть класс из разнородных 12 машин, на лучшей из которых Celeron1000-256MB расчетная задача считается 11 часов. Мне интересно, при организации кластера чего можно добиться. Если вместо указанных 11 часов получится минимум 10 - то не стоит и заморачисаться. Насколько узким местом является 100 Мбит-й интерфейс? Очень старые машины - будут тормозом, никак не проявятся, или добавят минимум к производительности?

anonymous
()

Моделирование ...

2Die-Hard:
> А на что скорость - то влияет? Типа, высоту эффективно увеличивает?
Если рассуждать в терминах полное давление-температура то скорее
уменьшает ;) На самом деле внешняя скорость в данном случае нужна
для расчета параметров пограничного слоя вблизи точки вдува - это сильно
влияет на параметры возле разрушенного илюминатора 

Про схему - К. Флетчер - Вычичслительные метода в динамике жидкостей Т2
            стр. 480 (Cхема явная)
            или C.A.J. Fletcher - "Computational Techniques for Fluid Dynamics 2" - страницу не скажу ;) 

Единственное отступление от классики это схема с весами с адаптивными 
весовыми коэффициентами.

2anonymous:

Навскидку тяжело но попробую:

>- отдельно на одном Athlon 1000Mhz TB - отдельно на одном Athlon >1700+ XP 

на 1000 TB примерно 2 часа 
на на XP1700+ пропорционально частоте - 1.43 часа   


>- при замене в Вашей схеме 100 Мбит на 1 Гбит - на 8-ми нодном 100 >Мбит кластере,

Для данной задачи изменение вряд ли было бы очень значительным... 

Количество данных - передаваемых по сети прямо пропорционально
размеру задачи и слабо зависит от числа узлов. Повторю что речь 
идет о данной конкретной задаче - фокус в том что для конкретной
программы существует некое оптимальное число slave-процессов
которое в данном конкретном случае равно 2-м - со стороны железа оно определяется производительностью узлов и скоростию среды передачи
в 8 узловом кластере эти 2 процесса просто влезут на 2 самых быстрых
узла и будут периодически менятся местами для обеспечения равномерности загрузки - это кстати хорошо видно на Load-диаграммах
gkrellm-ов 



> состоящим из одинаковых Athlon1000 - при добавлении к Вашей схеме >одного старого P166MMX

Скорее всего никак - хотя повторюсь, все зависит от задачи -
на memory expensive задаче все будет зависеть в большей степени от скорости среды передачи а на cpu expensive от скорости узлов.

>Эти ответы были бы интересны для оценки потенциала кластерных >вычислений. Допустим есть класс из разнородных 12 машин, на лучшей из >которых Celeron1000-256MB расчетная задача считается 11 часов. Мне >интересно, при организации кластера чего можно добиться. Если вместо >указанных 11 часов получится минимум 10 - то не стоит и >заморачисаться. Насколько узким местом является 100 Мбит-й интерфейс? 

Это опять же зависит от задачи - но в любом случае кластер с более
чем 10 узлами я бы не стал делать на 100 мегабитке - она годится разве 
что на "попробовать" - ImHO Ethernet вообще слабо годится как среда
передачи для сильносвязанных задач НО существует огромное число других
факторов (например финансовый) которые заставляют рассматривать такое 
решение ...

>Очень старые машины - будут тормозом, никак не проявятся, или добавят >минимум к производительности?

По опыту они ничего не добавляют (кроме опыта построения кластеров :))

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