Моделирование порыва ветра.

[4elovek 0]

Уважаемые специалисты, недавно столкнулся со следующей проблемой. Требуется смоделировать порыв ветра, воздействующий на дирижабль. Основная проблема заключается в том, что при моделировании порыва в виде сопла, из которого инжектируется струя, не получается чёткого фронта , то есть ступенчатого порыва.

Вопрос: можно ли смоделировать такой фронт с помощью пользовательской функции, которая описывала бы неравномерный поток? Скорость по оси ординат (порыв поперечный) должна быть функцией абсциссы. Кроме того, если использовать принцип обратимости, то возмущение должно перемещаться навстречу дирижаблю с течением времени.

Видел примеры, где скорость по оси абсцисс является функцией ординаты, но там в основном гладкие функции, а можно ли сделать резкий переход, который бы не размазывался при удалении от границы в центр расчётной области. Кто может подскажите как будет выглядеть такая функция.

/***********************************************************************

vprofile.c

UDF for specifying steady-state velocity profile boundary condition

************************************************************************/

#include "udf.h"

#include "math.h"

DEFINE_PROFILE(inlet_x_velocity, thread, position)

{

real x[ND_ND]; /* this will hold the position vector */

real y,z;

face_t f;

// volume_t f;

real V_an,G,H;

V_an=15.0;

begin_f_loop(f, thread)

{

F_CENTROID(x,f,thread);

y = x[1];

if(y>=0.0)

{

F_PROFILE(f, thread, position) = V_an;

}

}

end_f_loop(f, thread)

}

В первом приближении подошёл бы вот такой вариант, но помогите пожалуйста перевернуть его так, чтобы значение Vy менялось бы от икса. Заранее спасибо.

Изменено 4 октября 2011 пользователем 4elovek

[Disco 1]

Привет. У меня щас похожей проблемой занимаюсь (см. соседнюю тему)

В первом приближении подошёл бы вот такой вариант, но помогите пожалуйста перевернуть его так, чтобы значение Vy менялось бы от икса. Заранее спасибо.

я сделал так:

pressure_far_field

направляющий вектор потока X-компонент - 1 (constant)

а Y-компонент (udf)

/***********************************************************************

profile.c

************************************************************************/

#include "udf.h"

DEFINE_PROFILE(inlet_y_velocity, thread, position)

{

float t, velocity;

face_t f;

t = RP_Get_Real("flow-time");

velocity = (0.025*sin(0.5*t)+0.05);

begin_f_loop(f, thread)

{

F_PROFILE(f, thread, position) = velocity;

}

end_f_loop(f, thread)

}

Вроде дует нормально

Изменено 5 октября 2011 пользователем Disco

[4elovek 0]

Привет. У меня щас похожей проблемой занимаюсь (см. соседнюю тему)

я сделал так:

pressure_far_field

направляющий вектор потока X-компонент - 1 (constant)

а Y-компонент (udf)

А почему вы задаёте модель сжимаемого потока, ведь это глубокий дозвук? Кроме того, интересно, а как ведёт себя поток в середине расчётной области, просто у меня происходит разбухание струи с изменением профиля скорости в ней. Кроме того в невозмущённой части возникают два вихря, с каждой стороны фронта. Проблема в том,что скорости в них достигает 6 метров в секунду при скорости самого порыва 15 метров в секунду. Как у вас с этими вопросами обстоит дело?

Большое спасибо за функцию. Как я понял у вас поток периодический и привязан к во времени. Судя по вашей функции у вас поток меняет скорость разом во всей расчётной области, а мне нужно смоделировать именно фронт, то есть перепад скоростей внутри расчётной области в один момент времени желательно с постоянной зоной смешения, не расширяющейся при удалении от границы.

F_CENTROID(x,f,thread);

y = x[1];

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

Изменено 6 октября 2011 пользователем 4elovek