Перейти к публикации

3D расчеты двухфазной неравновесной теплогидродинамики


Рекомендованные сообщения

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

Одно время на сайте <noindex>http://www.inel.gov/relap5/relap5-3.htm</noindex> можно было скачать исходники программы. Язык программирования - Фортран.

Сайт продолжает активно функционировать и развиваться. В свое время я скачал исходники. Причем тогда не было никаких ограничений по распространению программы или ограничения лицензии.

Если уважаемым участникам форума будет интересно освоить этот мощный комплекс, по которому проводится обоснование безопасности АЭС с водоохлаждаемыми реакторами в различных аварийных и переходных режимах, то я готов предоставить все имеющиеся у меня материалы.

Считаю целесообразным открыть отдельный раздел по обсуждению этого программного комплекса.

Жду Ваших откликов и предложений.

Ссылка на сообщение
Поделиться на других сайтах


Очень интересно покопаться... судя по сайту, там расчеты далеко не ограничиваются двухфазными потоками. много уж вкусного по атомной области

division494 собака gmail точка ru

Изменено пользователем Hren
Ссылка на сообщение
Поделиться на других сайтах

Странно, что так мало откликов.

Вроде бы задач с двухфазной гидродинамикой очень много, а коды CFD такие задачи не могут решать, особено с неравновесными потоками.

Откликнувшиеся форумяне имеют ли опыт в решении подобных задач? Хотя бы в самой простой постановке 1D?

Текст кода на Фортране тоже немалый. Для работы с ним надо изучить мануалы, которые есть на оффсайте.

Вопрос: ознакомился ли кто с мануалами?

Есть ли у него своя задача для расчета?

Ссылка на сообщение
Поделиться на других сайтах

Вроде бы задач с двухфазной гидродинамикой очень много, а коды CFD такие задачи не могут решать, особенно с неравновесными потоками.

А как позиционируются в этой области Phoenics от CHAM и Гиневский из МЭИ?
Ссылка на сообщение
Поделиться на других сайтах

А как позиционируются в этой области Phoenics от CHAM и Гиневский из МЭИ?

Никак :no:

Здесь надо, хотя бы, спинодаль от бинодали отличать!

Ссылка на сообщение
Поделиться на других сайтах

Странно, что так мало откликов.

Я не специалист по многофазным течениям.

Хотя именно <noindex>Phoenics используют в МЭИ</noindex>. Там у них именно такие течения. Это их специфика.

Мы здесь задавали вопрос более простой по горению и также получили мало откликов.

Полагаю большинство решает просто задачи гидродинамики.

Ссылка на сообщение
Поделиться на других сайтах

Я не специалист по многофазным течениям.

Хотя именно <noindex>Phoenics используют в МЭИ</noindex>. Там у них именно такие течения. Это их специфика.

Мы здесь задавали вопрос более простой по горению и также получили мало откликов.

Полагаю большинство решает просто задачи гидродинамики.

К сожалению, МЭИ сейчас уже нет. Есть конгломерат техуниверситетов. Какие "течения" в МЭИ :flush:

я бы не взялся обсуждать. Если иметь в виду кафедру АЭС, которую я имел честь заканчивать (еще в МЭИ), то там Phoenics не пользуют. :bleh:

Вообще Релап считает пароводяные потоки, т.е. с водой в разных агрегатных состояниях (кроме льда).

Здесь и испарение и конденсация в разных видах, с теплообменом и без, с равновесием и без. Плюс специфичная задача нейтронной кинетики в активной зоне реактора.

Горение тоже есть. Циркониевый сплав горит в паровой среде с выделением водорода, который тоже может гореть и взрываться.

Ссылка на сообщение
Поделиться на других сайтах

Ну у вас же есть документация по Фениксу и техническая поддержка. Задайте им вопрос.

Ссылка на сообщение
Поделиться на других сайтах

Серьезность программы я определяю без самого кода.

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

Когда сравниваются сходные программы, напр., CATHARE, ТРАС или Relap, то обращаешь внимание на уровень эмпирических и замыкающих соотношений.

Еще понятно, что студентам не под силу пользоваться вышеупомянутыми программами, это дело непростое.

Ссылка на сообщение
Поделиться на других сайтах

Странно, что так мало откликов.

Вроде бы задач с двухфазной гидродинамикой очень много, а коды CFD такие задачи не могут решать, особено с неравновесными потоками.

Откликнувшиеся форумяне имеют ли опыт в решении подобных задач? Хотя бы в самой простой постановке 1D?

Текст кода на Фортране тоже немалый. Для работы с ним надо изучить мануалы, которые есть на оффсайте.

Вопрос: ознакомился ли кто с мануалами?

Есть ли у него своя задача для расчета?

Еще раз здравствуйте.

Мой опыт решения ограничивается только следующей, собственно говоря, и интересующей меня на данный момент задачей: впрыск воды в горячий скоростной поток газа (воздуха). Задача решалась в одномерной постановке (на базе ОДУ) и в многомерных постановках при использовании пакета Fluent (Discrete Phase Model).

Повторюсь, решение именно этой задачи меня и интересует.

С мануалами именно что ознакомился. Думаю, Вы согласитесь, что аппетит приходит во время еды: хорошо бы начать с простейшего примера. Хотя бы течение в трубе.

Ссылка на сообщение
Поделиться на других сайтах

А с помощью решения уравнения Больцмана методом решеток получают что-то серьезное?

Ускорение расчетов на GPU достигает 60-80 раз.

Ссылка на сообщение
Поделиться на других сайтах

А с помощью решения уравнения Больцмана методом решеток получают что-то серьезное?

Для программ расчета пароводяных потоков давно сделаны проверочные тесты.

Самый примитивный - тест Эдвардса и О'Брайена. Это моделирование открытия отверстия на конце трубки с горячей водой под давлением. На этом тесте "засыпался" не один знаменитый пакет CFD!

Теперь апологеты CFD(с этого форума) старательно избегают обсуждение этого вопроса. :biggrin:

Вот и в Вашем случае хорошо бы потестировать :doctor:

Ссылка на сообщение
Поделиться на других сайтах

На CFX уже повторили тест Эдвардса с уравнением состояния реального газа (не запомнил наименования).

Пики смазали, конечно же, но в целом намного лучше чем раньше.

А вообще, если в коммерческих кодах чего-то нет, то это просто означает, что это клиентам этих кодов эта тема не интересна, а не потому что это нельзя реализовать.

Ссылка на сообщение
Поделиться на других сайтах

Да. Это вопрос приделывания фазового перехода к основному коду.

Где-то получается глубже, где-то мельче.

Обычно в коллективы разработчиков добирают специалистов и они приделывают.

Надо прокачать свою важность и вас могут призвать. :rolleyes:

Возраст в частных фирмах не имеет особого значения, был бы опыт.

Ссылка на сообщение
Поделиться на других сайтах

В частных фирмах все стараются делать не ради красоты, а ради рубля.

Если есть потребность на рынке - сделают. А если нет, то зачем лишние проблемы себе на плечи? Сделаешь ветку кода, а потом ее при каждом релизе тестируй, хотя никто ей не пользуется.

Ссылка на сообщение
Поделиться на других сайтах

Верно. Поэтому если конденсация, то равновесная. И т.п.

Тогда можно самому дописать.

Fluent: через UDS+UDF.

И продать им? :rolleyes:

Ссылка на сообщение
Поделиться на других сайтах

На CFX уже повторили тест Эдвардса с уравнением состояния реального газа (не запомнил наименования).

Пики смазали, конечно же, но в целом намного лучше чем раньше.

А вообще, если в коммерческих кодах чего-то нет, то это просто означает, что это клиентам этих кодов эта тема не интересна, а не потому что это нельзя реализовать.

А что такое - " уравнение состояния реального газа"?

В тесте Эдвардса и О'Брайена есть только вода. Уравнения состояния воды давно известны и есть куча программ расчета свойств воды и пара. Неужели Вы считаете. что "прилепив" нужные свойства вещества, можно решить проблему?

Может быть речь идет о разных тестах?

Да. Это вопрос приделывания фазового перехода к основному коду.

Где-то получается глубже, где-то мельче.

Обычно в коллективы разработчиков добирают специалистов и они приделывают.

Надо прокачать свою важность и вас могут призвать. :rolleyes:

Возраст в частных фирмах не имеет особого значения, был бы опыт.

Если Вы сравните физмодели и исходные уравнения из Relap-5 (в мануалах это подробно расписано) и традиционных здесь пакетов CFD, Вы должны увидеть существенные отличия. И тогда Вы поймете что это НЕ "вопрос приделывания фазового перехода к основному коду"!

Верно. Поэтому если конденсация, то равновесная. И т.п.

Тогда можно самому дописать.

Fluent: через UDS+UDF.

И продать им? :rolleyes:

Вспоминается басня Крылова "Слон и Моська", особенно последняя строка :biggrin:

В частных фирмах все стараются делать не ради красоты, а ради рубля.

Если есть потребность на рынке - сделают. А если нет, то зачем лишние проблемы себе на плечи? Сделаешь ветку кода, а потом ее при каждом релизе тестируй, хотя никто ей не пользуется.

Программы типа Relap делали и продолжают совершенствовать НЕ частные фирмы.
Ссылка на сообщение
Поделиться на других сайтах

А что такое - " уравнение состояния реального газа"?

В тесте Эдвардса и О'Брайена есть только вода. Уравнения состояния воды давно известны и есть куча программ расчета свойств воды и пара. Неужели Вы считаете. что "прилепив" нужные свойства вещества, можно решить проблему?

Может быть речь идет о разных тестах?

Тест Эдвардса он один такой.

Что такое уравнение состояния можно узнать в каком-нибудь учебнике.

;)

Касаемо CFX постановка двухфазная как минимум. Уравнение состояния реального газа помогает при моделировании свойств вещества в зоне вблизи фазового перехода.

Если Вы сравните физмодели и исходные уравнения из Relap-5 (в мануалах это подробно расписано) и традиционных здесь пакетов CFD, Вы должны увидеть существенные отличия. И тогда Вы поймете что это НЕ "вопрос приделывания фазового перехода к основному коду"!

А в Релапе в основе Навье-Стокс считается или что?

Ссылка на сообщение
Поделиться на других сайтах

Тест Эдвардса он один такой.

Что такое уравнение состояния можно узнать в каком-нибудь учебнике.

;)

Касаемо CFX постановка двухфазная как минимум. Уравнение состояния реального газа помогает при моделировании свойств вещества в зоне вблизи фазового перехода.

А в Релапе в основе Навье-Стокс считается или что?

А как в СFX моделируется метастабильное состояние жидкости или пара?

А каково уравнение состояния реального газа в метастабильном состоянии?

И в каком учебнике про эти нетривиальные вещи написано?

А каков механизм перехода из метастабильного состояния в стабильное в CFX?

Да много еще подобных вопросов можно задать!

И еще одно соображение. В атомной отрасли большой прагматизм. Если бы была возможность использовать CFX, то его бы купили в самом полном варианте и и давно бы использовали вместо того, чтобы "городить" Релапы. Но, почему-то картина обратная.

На вопрос: "А в Релапе в основе Навье-Стокс считается или что?"

Ответ - не Навье-Стокс.

Подробно модель напр. в этом документе:RELAP5/MOD3 CODE MANUAL

NUREG/CR-5535

INEL-95/0174

(Formerly EGG-2596)

VOLUME I: CODE STRUCTURE, SYSTEM MODELS, AND SOLUTION METHODS

Ссылка на сообщение
Поделиться на других сайтах

Присоединяйтесь к обсуждению

Вы можете опубликовать сообщение сейчас, а зарегистрироваться позже. Если у вас есть аккаунт, войдите в него для написания от своего имени.
Примечание: вашему сообщению потребуется утверждение модератора, прежде чем оно станет доступным.

Гость
Ответить в тему...

×   Вставлено в виде отформатированного текста.   Вставить в виде обычного текста

  Разрешено не более 75 эмодзи.

×   Ваша ссылка была автоматически встроена.   Отобразить как ссылку

×   Ваш предыдущий контент был восстановлен.   Очистить редактор

×   Вы не можете вставить изображения напрямую. Загрузите или вставьте изображения по ссылке.

  • Сейчас на странице   0 пользователей

    Нет пользователей, просматривающих эту страницу.




  • Сообщения

    • Александр1979
      Информация по номерам для заказа есть в руководстве "B-65515", также информация есть в руководствах "B-65302" и "B-65262"      
    • gudstartup
      а чего тут искать то первый это абсолютный на 1000000 имп\об bi типа второй ai инкрементальный xз но нав 4000000 имп\об разъем подключения описан в схеме. инфа передается по цифровому протоколу обмена. кое что есть тут   в старых описаниях таких датчиков у фанука нет  
    • anykeyto
      У нас AC. Так то я понимаю, что G108 это переключить режим шпинделя из оборотов в координату С. Но я почему про нее спросил, потому что в мануале к G07.1, в примере программки указывалось таким образом, что сначала включали G108 (хотя это может для токарки, а для фрезера и не надо совсем), а потом уже остальное.   G07.1 это как раз цилиндрическая интерполяция, то есть преобразование оси С из угла в миллиметры, то есть как бы разворачивает цилиндр в плоскость. Ну как я это понял. Проблема у меня встала в том, как именно программируются такие действия. Пробовал на основе примеров из мануала, стойка выдает ошибки. в чем проблема непонятно из-за китайской стойки и некорректного перевода. Ну и вообще хочется пощупать, как это вообще работает именно на фрезерном ОЦ.   Тут еще засада в том, что на 5-оснике из-за RTCP станок не дает одновременно перемещаться одновременно по осям XY и С например...  Я понимаю, что такое надо в САМ системах делать, но есть потребность в упрощенном варианте и со стойки писать, у нас с технологами проблема... G68.1, 68.2, 68.3  это преобразование характеристической СК, развороты короче координат G53.1, G53.2 - фактически активирование поворота стола с RTCP и без... 
    • ДОБРЯК
      А собственные вектора вы как будете находить для любой размерности? :=)
    • Sivolap36
      Добрый день. Прошу подсказать в каком мануале есть описание энкодеров моторов, где-то видел, сейчас нем могу вспомнить.  Вот например два энкодера по ним ищу инфу  
    • Рома калужский
      это код подключения шпинделя как ось "С" но при этом он не перемещает револьверку, т.е. не работает по программе
    • Горыныч
      Еще совет: строить нужно только половину, если модель симметричная, тогда Солид меньше "косячит". И соглашусь с коллегами, идеальный вариант для гладкого построения - два сечения и две направляющие кривые, т.е. 4 управляющие кривые (кромки). При необходимости можно добавлять не только промежуточное сечение, но и направление. Не нужно пятаться построить всю поверхность одним фичером.  Пусть поверхности будут лоскутами, но соединяются по касанию или кривизне. 
    • ДОБРЯК
      Откуда появится смещение центра масс колеса, если колесо сбалансировано.  Жесткий сбалансированный вал, это условие данной задачи. На 4000 оборотах эти силы будут минимальны. Это разговор ни о чем в данной задаче. :=) Но даже если вы абсолютно точно определите эти частоты, с учетом всех сил. Что дальше с этими частотами делать? На основании чего делать вывод о работоспособности изделия? Это в учебной работе получили число, работу приняли и сдали в архив. Какая практическая ценность всех этих уточнений в данном случае?:=)  
    • Борман
      Если еще не поздно... я и сам погрузился бы. Физика такая, для вала с одним колесом... F - это понятно, центробежное смещение центра масс колеса M - момент от центробежных сил, который не появится, если колесо не повернется при изгибе вала. Обрати внимание, что BW-момент выворачивает ось в сторону действия силы F, т.е. при обратной прецессии они действуют сообща, и в этом опасность.   Отсюда будут разные собственные частоты при прямой и обратной прецессиях.  Теперь осталось понять, что такое прямая и обратная прецессии :), и почему момент при них направлен в разные стороны. По первой части вопроса см. мою картинку в этой теме от 27 августа 2013. По второй части вопроса - нужно включать воображение и усердно рисовать на бумажке. Мне пока сложно.    Физика именно такая. Нужно как то оценить величину этого момента и приложить его в софте, который не умеет ДК. А с тестовыми значениями можно и так порешать, и на выходе должна быть вилка. Т.е. F - задает типа "биссектрису" вилки (престресс-модальник), а момент - угол вилки.  
    • boomeeeer
      @anykeyto Я не оператор. Настраивал работу этой функции на токарно-фрезерном, где нужно было позиционировать шпиндель, для нарезания пазов. Что такое G07 не знаю. Вот такой пример для фрезерного станка нашёл   G108 B0; переключение на ось B для управления положением; G68.2 X10 Y10 Z20 I0 J90 K0 Положение заготовки (10,10,20) Угол Эйлера для определения характеристической системы координат; G53.2 Управление направлением инструмента   G109 B0 - отключит режим
×
×
  • Создать...