/* turbulence model: RSM hybrid (fixed value routine - turbulent shear stresses determined via boussinesq hypothesis*/
/* UDF (User Defined Function) Routine for FLUENT                                                             */
/* see also the paper:                                                                                        */
/* Implementation and comparison of different turbulence models for three dimensional wall jets with Fluent   */
/* FLUENT CFD Forum 2005 Bad Nauheim, Germany								      */
/* Heschl Ch*., Sanz W.**, Klanatsky P.*                                                                      */
/* *  Fachhochschulstudiengänge Burgenland GmbH, Steinamangerstraße 21, A-7423 Pinkafeld                      */
/* ** TU Graz, Institut of Thermal Turbomachinery and Machine Dynamics, Inffeldgasse 25, A-8010 Graz          */
/* Contact: christian.heschl@fh-burgenland.at, wolfgang.sanz@tugraz.at, peter.klanatsky@fh-burgenland.at      */

/* Die Autoren übernhemen keine Verantwortung oder Haftung für den Inhalt der unten angeführten UDF-Routinen  */



#include "udf.h"


DEFINE_PROFILE(fixed_uv, t, i)

{
  cell_t c;

  begin_c_loop (c,t)
  {
   F_PROFILE(c,t,i) = - 2.*C_MU_T(c,t)/C_R(c,t)* 0.5 * (C_DUDY(c,t)+C_DVDX(c,t));
  }
  end_c_loop (c,t)
}

DEFINE_PROFILE(fixed_uw, t, i)

{
  cell_t c;

  begin_c_loop (c,t)
  {
   F_PROFILE(c,t,i) = - 2.*C_MU_T(c,t)/C_R(c,t)* 0.5 * (C_DUDZ(c,t)+C_DWDX(c,t));
  }
  end_c_loop (c,t)
}

DEFINE_PROFILE(fixed_vw, t, i)

{
  cell_t c;

  begin_c_loop (c,t)
  {
   F_PROFILE(c,t,i) = - 2.*C_MU_T(c,t)/C_R(c,t)* 0.5 * (C_DWDY(c,t)+C_DVDZ(c,t));
  }
  end_c_loop (c,t)
}