Комментарии 27
CFD по определению динамическое моделирование, а вот для установившихся процессов любой студент посчитает, это типовая задача матфизики.
Сравнение результатов, полученных по авторской программе, с результатами CFD моделирования по установившимся (после нанесенного возмущения) значениям…
Ну да, это все равно что сказать, что автор решил на пальцах задачу движения трех тел - после того, как они все упали на планету и движение прекратилось.
Вы совершенно правы. В цитате, приведенной Вами, говорится о совпадении стационарных температурных полей в конце переходного процесса, полученных при CFD моделировании и при использовании изобретенного подхода. Получение стационарных полей, с использованием уравнения теплопроводности, действительно не представляет сложности и совпадение “стационаров” лишь подтверждает адекватность работы изобретенного подхода. Я бы хотел обратить Ваше внимание на рисунки 4-5, на которых представлено сравнение изменений температурных полей при протекании переходных процессов. На данных рисунках видно, что “динамическое моделирование” также было проведено и было получено хорошее совпадение результатов.
Я не понял, сравнения проводились для всей активной зоны с подробным ее описанием или на том маленьком упрощённом фрагменте?
Сравнение проводилось на упрощенном фрагменте, поскольку CFD моделирование всей зоны я не могу себе позволить. Однако мой подход возможно масштабировать для анализа зоны целиком. Осталось только найти результаты CFD моделирования всей зоны для верификации
CFD моделирование всей зоны я не могу себе попозволить.
Какое количество элементов сетки получается? Если не учитывать радиационный теплообмен ( а вы его как-нибудь его учитываете?), то 3 млн за 4000 интеграцией реально посчитать за неделю.
Но есть несколько но.
Во-первых во fluent есть возможность ускорить вычисления за счёт использования видеокарты. Мне неизвестно насколько это эффективно, сами разработчики обещают чуть ли не в сотни раз на А100.
Во-вторых в CFD постановке это расчёт сопряжённого теплообмена жидкости и твёрдого тела. Для сходимости таких рачётов требуется на порядок больше итераций, чем для обычных CFD расчётов. Во ansys CFX есть есть возможность задать разные шаги по времени для жидкости и твёрдого тела, что позволяет в ряде случаев уменьшить необходимое количество итераций.
Однако мой подход возможно масштабировать для анализа зоны целиком.
Абсолютно неочевидно, что при масштабировании и увеличении сложности не упадет драматически точность. Чем проще модель, тем точнее ее обсчитывают простые решения.
Кстати, курьез из жизни. В одной фирме мы продавали нашу крутейшую прогу расчета трехмерного поля силы радиосигнала от сотовых вышек, с учетом зданий их типов, размеров, лесов и прочего. Картинка получалась замечательная, очень красивая. Прогу продавали за очень дорого для сотовых компаний, для оптимизации и предсказания покрытия. На формулах расчета поля была сделана докторская диссертация. А потом как-то раз я с коллего ради прикола посчитал поле по простейшией формуле для вакуума и данные сравнили и той крутейшей моделью. Точнее, сравнили с данными по реальнвм замерам поля в том районе. И средняя ошибка крутой модели оказалась примерно такой же, как и для простой формулы для вакуума.
Мы тихо шепнули начальнику про результаты, он посмотрел, нахмурился и сказал нам никогда и никому это больше не показывать :) А фирму потом продали за хорошие бабки и все получили плюшек ;)
Можно уточнить, натурные испытания будут, и если где, то когда? Планирую отпуск в другом полушарии.
Число Рейнольдса чему равно?
Число Рейнольдса примерно равно 2,3*10^5
Число Рейнольдса примерно равно 2,3*10^5
Да, не взялся бы я предсказывать, какая там эффективная теплопровоность будет. Ну и, судя по началу статьи, речь про воду или вообще про газ? Т.е. то, что называется semi-transparent media? В этом случае, если оптическая толщина маленькая, о теплопроводности вообще лучше не говорить.
Норм работа, от/для молодого учёного, аспиратнта. Зачот!
Что то мне подсказывает, что Вы в качестве "расчетной модели активной зоны" используете один ряд ТВЭЛ из "квадратной" ТВС Вестингауза.
Так ли это в действительности?
Прошу не минусовать, ОФФТОП...
Я просто мечтаю сесть в компании пьяных баб и выдать с глубокомысленным видом:
"На протяжении последних шести лет я участвовал в исследованиях, посвященных моделированию стационарных и нестационарных температурных полей активных зон ядерных реакторов..."
Спасибо, очень интересный подход!
Вопрос по поводу функции P - функция, отражающая влияние граничных условий.
Какие у нее аргументы?
Благодарю за отзыв!
Функция P - является сложной зависимостью, охватывающей множество параметров. Вид этой зависимости, а также более подробное описание предлагаемого подхода, Вы сможете найти здесь https://nuclear-power-engineering.ru/authors/podgornysk/
Спасибо, я посмотрел статьи по ссылке, про P нашел что:
Проведенный анализ показал, что для активной зоны реактора с водой под дав
лением может быть использовано соотношение для определения отвода тепла с по
верхности твэла в нестационарном процессе, полученное в работе [10] для актив
ной зоны высокотемпературного газоохлаждаемого реактора:
Ссылку [10] (Кузеванов В.С., Подгорный С.К. Газоохлаждаемые реакторы. Профилирование и интенсификация теплообмена. – Кишинев: Palmarium Academic Publishing, 2019. – 84 с) найти в открытом доступе не могу:(
Меня скорее интересует не сложность самой зависимости, а ее входные аргументы, а точнее чувствительность решения данной задачи к этим аргументам.
Меня скорее интересует не сложность самой зависимости, а ее входные
аргументы, а точнее чувствительность решения данной задачи к этим
аргументам.
Необходимости в такой проверке не возникло, поэтому я не смогу ответить Вам на Ваш вопрос.
Ссылку [10] (Кузеванов В.С., Подгорный С.К. Газоохлаждаемые реакторы. Профилирование и интенсификация теплообмена. – Кишинев: Palmarium Academic Publishing, 2019. – 84 с) найти в открытом доступе не могу:(
Всё верно, монографии в открытом доступе нет. Однако она не содержит какой либо информации, которая была бы неопубликована.
А как получать функцию стока тепла qv.st для применения данного метода для аналогичных задач?
Необходимо было во введении к статье указать ссылку на работу в журнале "Известия вузов", таким образом определенный процент читателей мог бы сразу начать более подробное ознакомление с материалом. Теперь что касается содержания статьи:
Отсутствует обзор использующихся в отечественной и зарубежной практиках программных средств для теплогидравлического анализа активных зон ядерных реакторов. Если говорить про иерархию кодов в соответствии с масштабами расчетных областей в ЯЭУ, то можно выделить системные (system), субканальные (subchannel), CFD, DNS.
Следствием из предыдущего пункта является некорректное сравнение времени счета вашего программного средства с CFD-моделированием, т.к. данные коды находятся на разных ступенях в приведенной градации и предназначены совершенно для различного рода анализа. Ваше программное средство логичнее было сравнивать с классическими одномерными кодами, в которых сложное движение теплоносителя в пределах пучков стержней сводится к рассмотрению одномерного течения внутри эквивалентной трубы.
Т.к. отсутствует сравнение с одномерными кодами, а именно они используются в упомянутых вами тренажерах, то область применения вашего ПС является крайне туманной. В особенности учитывая тот факт, что наиболее интересные динамические процессы в ЯЭУ связаны с кипением теплоносителя, моделирование которого с помощью вашего ПС не было рассмотрено.
В целом основное развитие программных средств сейчас идет в сторону корректного моделирования локальных эффектов, которые раньше либо не учитывались вовсе, либо учитывались благодаря введению соответствующей степени консерватизма. Ваше программное средство идет в абсолютно противоположную сторону...
"Пара слов о реализации данного алгоритма" - пункт специально для Хабра... Вообще, если вы добавляете информацию о выборе ЯП, то он должен быть обоснован. Фраза "по моему личному мнению, Java обладает приятным синтаксисом." не является обоснованием (причем я не уверен, что даже создатели Java считали, что у данного ЯП приятный синтаксис). Для задач вычислительной физики подходят другие языки: Fortran, C++, Matlab, Python, Julia (произвольный порядок перечисления)
Необходимо было во введении к статье указать ссылку на работу в журнале "Известия вузов"
Мне показалось, что это будет рекламой других ресурсов.
Ваше программное средство логичнее было сравнивать с классическими одномерными кодами
Я не являюсь ни аспирантом ни работником какого-либо специального заведения, и, насколько мне известно, такие коды не лежат в открытом доступе. Соответствено, сравнивать мне было больше не с чем.
"Пара слов о реализации данного алгоритма" - пункт специально для Хабра...
Здесь Вы не правы. Данное программное средство было действительно написано на Java, и я не вижу ничего плохо в том, чтобы об этом упомянуть.
Предсказательная вычислительная гидродинамика