4.6 Vergleichsrechnung mit Hilfe des SST Modells - Lehrstuhl ...

Anhang A
#include "udf.h"
DEFINE_PROFILE(inlet_x_velocity, thread, index){
real p[ND_ND];
real v; //Geschwindigkeit
real vx; //Geschwindigkeit in x
real x; //Koordinate
real y; //Koordinate
real r; //Radius
real alpha; //Winkel
real pi = 3.141592654; //Pi
real n; //Drehzahl
face_t f;
n = 0.025; //Drehzahl
begin_f_loop(f, thread)
{
F_CENTROID(p,f,thread);
x = p[0];
x = x - 823.0; //Verschiebung auf globale Koordinaten
y = p[1];
y = y - 679.674; //Verschiebung auf globale Koordinaten
r = sqrt(x*x+y*y);
v = r*pi*n*2.0;
alpha = asin(y/r);
vx = v*sin(alpha);
F_PROFILE(f, thread, index) = vx;
}
end_f_loop(f, thread)
}
DEFINE_PROFILE(inlet_y_velocity, thread, index)
{
real p[ND_ND];
real v; //Geschwindigkeit
real vy; //Geschwindigkeit in y
real x; //Koordinate
real y; //Koordinate
real r; //Radius
real alpha; //Winkel
real pi = 3.141592654; //Pi
real n; //Drehzahl
face_t f;
n = 0.025; //Drehzahl
begin_f_loop(f, thread)
{
F_CENTROID(p,f,thread);
x = p[0];
x = x - 823.0; //Verschiebung auf globale Koordinaten
y = p[1];
y = y - 679.674;/ //Verschiebung auf globale Koordinaten
r = sqrt(x*x+y*y);
v = r*pi*n*2.0;
alpha = asin(y/r);
vy = v*cos(alpha); if ( x>= 0) {vy = vy * (-1.0);}
F_PROFILE(f, thread, index) = vy;
}
end_f_loop(f, thread) }
91