Добрый день! Прошу помощи.
Стоит задача написания UDF для случая, когда поток в канале некоторое время движется в одном направлении, затем направление меняется на противоположное.
Единственное, что удалось реализовать - это периодическое изменение знака у значения скорости для inlet. При этом опорное давление на Outlet всё это время равно = 0. Знаю что это неверно!
Но как в UDF прописать еще и периодическое изменение типов зон мне не удалось реализовать. Т.е. мне необходимо периодически менять (inlet_velocity на pressure_outlet) и (pressure_outlet на inlet_velocity) для двух зон.
Вот что у меня есть:
#include "udf.h"
#define time_0 0.0 /* объявление промежутков времени */
#define time_1 2.0
#define time_2 4.0
#define time_3 6.0
#define time_4 8.0
#define time_5 10.0
#define time_6 12.0
DEFINE_PROFILE(inlet_velocity, thread, position)
{
real x[ND_ND]; // added
real time;
real velocity;
face_t f;
begin_f_loop(f, thread)
{
F_CENTROID(x, f, thread); // added
time = RP_Get_Real("flow-time");
{
if ((time > time_0 && time <= time_1) | (time > time_2 && time <= time_3) | (time > time_4 && time <= time_5) )
{
velocity = 2.0; /* поток в "прямом" направлении */
}
else if ((time > time_1 && time <= time_2) | (time > time_3 && time <= time_4) | (time > time_5 && time <= time_6))
{
velocity = -2.0; /* поток в "обратном" направлении */
}
}
end_f_loop(f, thread)
F_PROFILE(f, thread, position) = velocity;
}
}
Заранее благодарен за помощь!
Здравствуйте, можете попробовать команду (или что-то типа того под ваши нужды)
define boundary-conditions modify-zones zone-type inlet pressure-inlet
Может быть вы имели в виду статическое давление на Outlet?
Добрый день! Спасибо большое за напутственное направление!
Буду пробовать.
Т.е. для таких "реверсивных периодических" течений можно обойтись и буз UDF? Это очень удобно
Добавить комментарий