free surface Моделирование в Ansys Fluent гидротехнических сооружений

[sunway 0]

Всем доброго дня!

Собственно вопрос заключается в следующем.

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

Нарыл на просторах интернета следующий UDF

Скрытый текст

#include "udf.h"
DEFINE_PROFILE (unsteady_velocity, thread, position)
 {
      
    face_t f;
  
  real t = CURRENT_TIME; 
  real velocity;
  
    begin_f_loop (f, thread)
    {
     velocity=0.02*t*t;
     F_PROFILE(f, thread, position) = velocity;
    }
    end_f_loop(f, t)
 }

 

по этому коду получилось задать переменную скорость по времени, а как сделать чтобы еще геометрия менялась?

 

И вот пример, того что я пытаюсь замоделировать Тыц!

[piden 2 874]

33 minutes ago, sunway said:

а как сделать чтобы еще геометрия менялась?

 

И вот пример, того что я пытаюсь замоделировать Тыц!

А где там изменение геометрии?

[a_schelyaev 366]

Мне для расширения собственного кругозора: во Флюенте надо кодить на Си, чтобы переменное по времени и пространству ГУ задать?

Изменено 24 ноября 2018 пользователем a_schelyaev

[an_rushin 128]

так точно, на си

 

[sunway 0]

В 23.11.2018 в 19:37, piden сказал:

А где там изменение геометрии?

Пардон, неточно выразился, геометрия не меняется

Мне надо чтобы при задании граничных условий в inlet скорость менялась по времени, от 0-4,5 м/сек например по линейно функции (но можно и квадратичной) и по оси "Y" была переменной.

 

А на моем видео как раз случай, когда скорость постоянная 4,5 м/сек

[sunway 0]

И есть еще вопрос, как дополнить UDF чтобы первые 5 секунд в inlet скорость была переменной, до 4,5 м/сек, а после константой, те же 4,5 м/сек.

[sunway 0]

Кажется частично я разобрался со своими хотелками, в части переменной скорости по времени с константной величиной

 

Скрытый текст

#include "udf.h"
DEFINE_PROFILE (unsteady_velocity, thread, position)
 {
      
    face_t f;
  
  real t = CURRENT_TIME;
  real time; // Время после которого начинается равномерная скорость
  real k; // коэффициент квадратичной функции 
  real velocity;
  time = 5.0;
  k=0.18;
  
    begin_f_loop (f, thread)
    {
        if ( t < time ) {
            velocity=k*t*t; // Переменная величина квадратичная функция
        }
        else {
            velocity=k*time*time; // Константная величина
        }
        F_PROFILE(f, thread, position) = velocity;
    }
    end_f_loop(f, t)
 }

 

[soklakov 1 475]

В 24.11.2018 в 20:02, a_schelyaev сказал:

во Флюенте надо кодить на Си, чтобы переменное по времени и пространству ГУ задать?

или загрузить таблицу из файла. вплоть до версии 19.1

с 19.2 появились expression а-ля cfx.