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

Нестационарная гидродинамика


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

Добрый день!

 

Ищу совета у опытных ребят, интересует нестационарный гидравлический расчет численным методом.

 

Задача трехмерная, необходимо определить расход и перепад давления на клапане при его открытии. Время открытия 5 сек.

 

Задача решается в нестационарной постановке средствами CFX или Fluent. Насколько я понял, решение можно считать сошедшимся, если на каждом шаге невязки падают более чем на три порядка. Применяю адаптивный шаг по времени с числом Куранта от 1 до 5. При таких параметрах шаг составляет порядка микросекунды, а само решение занимает несколько миллионов итераций.

 

Поток турбулентный, применена K-E модель со скалируемой пристеночной функцией.

 

Хотелось уточнить, такие классы задач действительно так трудоемки или заданные мной настройки не подходят для таких процессов?

 

 Также пытался решить нестационарную задачу с проходом среды через открытый клапан с выходом на установившейся режим и также получал микроскопический временной шаг и миллионы итераций :(

 

Дайте пару советов, если не трудно :)

 

Спасибо!

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


UnPinned posts
21 hours ago, Guterfreund said:

Дайте пару советов, если не трудно :)

Первый и главный: всегда подробно описывайте задачу. Даже не ее постановку и выбранный вами метод решения. А что нужно посчитать (и, желательно, зачем. Да...).

И потом уже: как вы подошли к решению и какие проблемы при этом возникли. Не зная задачи, сложно сказать, правильны ли ваши действия.

 

Открытие клапанов - это обычно экстремальные режимы на начальном этапе (сверхзвук у газа / кавитация у жидкостей), очень малые зазор =>, мелкая сетка. Вместе с тем, такие режимы (когда СЗ или кавитация) обычно длятся сотые / десятые доли секунды, и для типичных инженерных целей ими можно пренебречь при расчете выхода на стационар. 

Вы не упоминаете подвижные сетки и изменяющиеся ГУ. Что в вашей задаче нестационарного? Точно ли нельзя использовать несколько стационарных расчетов разных конфигураций?

 

Если нельзя поделиться геометрией/.cas конкретно вашей задачи, соорудите похожую.

 

21 hours ago, Guterfreund said:

При таких параметрах шаг составляет порядка микросекунды, а само решение занимает несколько миллионов итераций.

Да, возможно, если считать процесс в лоб, со всеми ударными волнами.. Шаг может вырасти, когда процесс станет более установившимся. А может и не вырасти. Может, придется его вручную изменять, или из journal-файла, или удф-ом. Может, в другом солвере (напр. coupled) все будет решаться на больших шагах. Но нужно ли это все и поможет ли в конкретно вашем случае - хз.

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

Я подозреваю, что те же особенности в задаче, что я тут ранее на форуме в FlowVision решал, где также два объема, дикий перепад и щель между ними.

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

Нужно иметь хорошее начальное приближение в области. Т.е. можно начать с полностью открытого клапана и потом его закрывать. Щель будет уменьшаться, скорость при постоянстве расхода расти. Соответственно шаг интегрирования снижаться будет из-за роста скорости, а также из-за уменьшения размера ячеек сетки в зоне зазора. Ну а там и сатурейшн прешша не загорами с кавитацией :)

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

 при постоянстве расхода 

Какое-то физически не реальное условие. Тут вроде что-то вроде лабиринтного уплотнения. Перепад давлений скорее будет постоянным ... 

Ссылка на сообщение
Поделиться на других сайтах
2 часа назад, Fedor сказал:

Какое-то физически не реальное условие. Тут вроде что-то вроде лабиринтного уплотнения. Перепад давлений скорее будет постоянным ... 

Тут, действительно, вопрос. Сопротивление клапана будет расти. Если на входе насос, то расход, по идее, будет падать. Или нет?

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

Если не нравится задавать расход и определять перепад давления, то задайте перепад давления и определите расход

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

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

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

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

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

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

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

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

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

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

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




  • Сообщения

    • Viktor2004
      Туда вроде наоборот заносится значение какое-то, но не считывается Смотрю чужие программы, там вроде только 8549 используют. При чем у них в программах циклы поболее моих, которыми я круги чертил O4198 ; #500=#8549 ; IF[ #8549 EQ 0 ] THEN ; #10409=0 ; ENDIF ; IF[ #8549 EQ 134 && #10409 EQ 0 && #10422 NE 0 ] THEN ; #10409=1 ; #10422=FIX[#10422 / 10] ; #10427=1 ; ENDIF ; IF[ #8549 EQ 48 && #10409 EQ 0 ] THEN ; #10409=1 ; #10422=#10422 * 10 + 0 ; #10427=1 ; ENDIF ; IF[ #8549 EQ 49 && #10409 EQ 0 ] THEN ; #10409=1 ; #10422=#10422 * 10 + 1 ; #10427=1 ; ENDIF ; IF[ #8549 EQ 50 && #10409 EQ 0 ] THEN ; #10409=1 ; #10422=#10422 * 10 + 2 ; #10427=1 ; ENDIF ; IF[ #8549 EQ 51 && #10409 EQ 0 ] THEN ; #10409=1 ; #10422=#10422 * 10 + 3 ; #10427=1 ; ENDIF ; IF[ #8549 EQ 52 && #10409 EQ 0 ] THEN ; #10409=1 ; #10422=#10422 * 10 + 4 ; #10427=1 ; ENDIF ; IF[ #8549 EQ 53 && #10409 EQ 0 ] THEN ; #10409=1 ; #10422=#10422 * 10 + 5 ; #10427=1 ; ENDIF ; IF[ #8549 EQ 54 && #10409 EQ 0 ] THEN ; #10409=1 ; #10422=#10422 * 10 + 6 ; #10427=1 ; ENDIF ; IF[ #8549 EQ 55 && #10409 EQ 0 ] THEN ; #10409=1 ; #10422=#10422 * 10 + 7 ; #10427=1 ; ENDIF ; IF[ #8549 EQ 56 && #10409 EQ 0 ] THEN ; #10409=1 ; #10422=#10422 * 10 + 8 ; #10427=1 ; ENDIF ; IF[ #8549 EQ 57 && #10409 EQ 0 ] THEN ; #10409=1 ; #10422=#10422 * 10 + 9 ; #10427=1 ; ENDIF ; IF[ #10422 GT 9999 ] THEN ; #10422=FIX[#10422 / 10] ; ENDIF ; IF[ #500 EQ 249 && #8500 EQ 4170 ] THEN ; #8500=4601 ; M99 ; ENDIF ; IF[ #500 EQ 240 ] THEN ; IF[ #8500 EQ 4174 || #8500 EQ 4173 || #8500 EQ 4177 ] THEN ; #8500=4170 ; M99 ; ENDIF ; IF[ #8500 EQ 4178 ] THEN ; #8500=4170 ; M99 ; ENDIF ; ENDIF ; IF[ #500 EQ 248 && #8500 EQ 4170 && D34.0 NE 0 ] THEN ; #8500=4174 ; M99 ; ENDIF ; IF[ #500 EQ 247 && #8500 EQ 4170 && D30.4 EQ 1 && D29.6 EQ 0 ] THEN ; #8500=4175 ; M99 ; ENDIF ; IF[ #500 EQ 246 && #8500 EQ 4170 && D35.0 EQ 1 && D35.4 EQ 1 ] THEN ; #8500=4178 ; M99 ; ENDIF ; IF[ #8549 EQ 152 && #10409 EQ 0 && #10422 EQ 5656 ] THEN ; #8500=4177 ; M99 ; ENDIF ; M99 ;  
    • nicomed
      Как победить проблему с блокированием сегментов эскиза не знаю. Поэтому пока полумеры... Ввел активацию чертежного листа после разрыва линии. Для разблокирования сегментов после закрытия программы двойным кликом активировать чертежный вид, где надо двигать сегменты эскиза. Добавил настройки. Можно настроить величину генерируемого зазора в  файле options.xml. По толщине и стилю линии, если кто подскажет как прочитать стиль и толщину выбранной кромки ( силуэтной кромки), то тогда можно будет "напосмотреть".   Пару раз натыкался на такой глюк - остается только один сегмент эскиза, грешу на рассинхронизацию работы солида и приложения. EdgeTrim_v3.3.zip
    • si-64
      да, это работает. спасибо. но как то особо не работал со слоями, нужно прямо заморачиваться(
    • nicomed
    • Stanislav
      Фрезеровка и токарка ЧПУ. Территориально Нижний Новгород, отправка в регионы ТК. Фрезеровка Алюминия и его сплавов, Латунь, Текстолиты, Сталь в небольших объемах. Токарка как цветные так и черный металлов и неметаллов. Габариты деталей для фрезеровки порядка 400мм*600мм. Оплата, пока что, наличные или на карту. Чертежи и модели отправлять на почту 0563506@mail.ru там же все согласования и обсуждения.
    • Александр1979
      Можно попробовать код считывать из переменной #8501. В руководстве написано что значение сохраняется в переменной до тех пор, пока не будет произведено чтение из неё.
    • si-64
      Можно ли сделать так чтобы в чертеже все размеры были черные, а в эскизе детали цвета управляемых и управляющих различались? настройка меняет цвет как в чертеже так и в эскизах: Зачем это нужно? В чертеже есть элементы, отличные от черного, поэтому вариант печати (сохранения в PDF) в ч/б не устраивает.
    • Artas28
      Шота получилось. Всё равно. Жду ответа от гуру NX, как сделать фаску красиво в 5 осях, что бы вектор инструмента менялся.  
    • Viktor2004
      Моя программа считывает его лишь тогда, когда в момент нажатия кнопки программа выполняет команду считывания переменной. Во всех остальных случаях код получается не считан
    • gudstartup
      но кто то же этот код считывает из mdi буфера раз его там не оказывается когда ваша программа его опрашивает? 
×
×
  • Создать...