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

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 пользователей

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




  • Сообщения

    • nicomed
      Поднял старые записи. Из того что было максимально похоже на то, что можно было назвать "получалось" : ModelDoc2 modelDoc = swApp.ActiveDoc; DrawingDoc drawing = modelDoc as DrawingDoc; MathUtility mathUtils = swApp.IGetMathUtility(); SelectionMgr selectionMgr = modelDoc.SelectionManager; Edge edge = selectionMgr.GetSelectedObject6(1,-1); View view = selectionMgr.GetSelectedObjectsDrawingView2(1, 0); //Не уверен - надо ли ... drawing.ActivateView(view.Name); double[] viewpositionArray = view.Position; // Здесь по идее можно и просто получить координаты вершины выбранного EDGE var curveParams = edge.GetCurveParams2(); var startPoint = new double[3]; startPoint[0] = curveParams[0]; startPoint[1] = curveParams[1]; startPoint[2] = curveParams[2]; MathPoint modelPoint =mathUtils.CreatePoint(startPoint); MathTransform modelTransform = view.ModelToViewTransform; MathPoint drawingPoint = modelPoint.MultiplyTransform(modelTransform); double[] drawingPointArray = drawingPoint.ArrayData; SketchManager sketchManager = modelDoc.SketchManager; // Первые три координаты даны для простоты тестирования sketchManager.CreateLine( 0, 0, 0, (drawingPointArray[0] - viewpositionArray[0]) / modelTransform.ArrayData[12], (drawingPointArray[1] - viewpositionArray[1]) / modelTransform.ArrayData[12], drawingPointArray[2] / modelTransform.ArrayData[12]);   Самому не совсем понятно поведение СкетчМенеджера при отрисовке линии от координат 0,0,0 - частенько слетает в произвольное место на виде. 
    • maxx2000
      на втором  скрине PLC сообщает что-то типа "старт программы запрещён потому что оси ХУ не связаны". Возможно станок в режиме точения должен заблокировать перемещение У, включить тормоз или что-то ещё для произвольного смещения под нагрузкой. Возможно косяк в тексте программы.
    • maxx2000
      @sneg0vik как это? Если ХУ не перпендикулярны то это уже не У, а скажем ось В. Нет, конечно в теории можно построить станок с неперпендикулярными ХУ. Вопрос зачем? Gildemeister CTX 420 это же токарно-фрезерный с осью У
    • sneg0vik
      На станке ось "X" является наклонной по отношению к оси "Y" (т. е. они не перпендикулярны). Поэтому надо полагать (учитывая текст ошибки) у вас отключена связь оси "Y" с осью "X".   Проверьте меняются ли координаты оси "X", когда вы двигаете ось "Y". И наоборот. Если - нет, то ищите как включить (наверное через параметры) взаимосвязь оси "Y" с осью "X".
    • ДОБРЯК
      Если ваш Заказчик и такую работу примет, то необходимо в шпоночное отверстие добавить и массу воды. Примерно 1/4 от той массы которую вы будете добавлять.  Да и еще и не забыть добавить момент инерции от воды.
    • Soprin
      В функции MOVE по F7.3 должно же из R270 переносить в R278?
    • niagara39
      Причина все еще не найдена... Но заметил, что сразу после включения в меню диагностики появляется следующее предупреждающее сообщение: 10155 M: Y-axis: Y-axis and X-axis not coupled , но причина его появления и его смысл неизвестны
    • maxx2000
      скажите, Вам уже делали эти детали? Вы себе представляете как и главное  чем, можно обработать внутренние углы без скругления на вашей детали? Разве что проволокой. НО оно же стоить будет как крыло от боинга. 
    • Sivolap36
      Спасибо за информацию. А как определили, что нижний  инкрементный? Спасибо. Как то я пролистал наверное. Благодарю
    • Jesse
      deleted этого я и хотел узнать у товарищей форумчан :-) вот для сферы оценка есть: присоед. масса воды равна половине от объёма сферы может кто в справочнике каком видел сколько для рабочих колёс присоединять надо..) мне тоже кажется там серьёзная отстройка должна получиться..   не спец в насосах и лопаточных машинах..)) у меня вот тяк: да это геморр... ещё и приложить его надо как-то к валу. Кстати, я вообще говоря собирался для модальника взять модель только вала, с колёса учесть в виде сосредоточ. масс с моментами инерции. Законно ли это, не сильно исказит расчёт СЧ? И ещё непонятка: куда присоединять эту массу? Прям в шпоночное отверстие вала?: Или вдоль всей окружности куда садится колесО?   не не , там не биссектриса. \Гляньте мою картинку выше из Вики. И все картинки такие для ДК что в сети нашлись. Там явно BW и  FW явно не зеркальные отражения друг друга относ-но гориз-й прямой, проведённой через ординату СЧ для нулевой частоты вращения ротора...   откуда вы знаете что сбалансировано? идеальной балансировки не существует   вот прям так в отчёте и напишу!:D точка Меня другой вопрос волнует.. Что центробеж. сила и гироскоп. момент пропорц-ны квадрату частоты вращенеия, т.е. на больших частотах может появиться вилка... Мне нужно где-то найти оценку, возможно опять таки относительно СЧ когда нет вращения..
×
×
  • Создать...