[Mohinder Suresh]

Hello Fluent users!
I need some help regarding the application of an x-component momentum source term.

I want to add the source term:
S = -1/x*C*F

to the x-momentum equation. "x" represents the x coordinate, "C" is a constant and F is the inviscid Flux in the x-direction.




In the FLUENT manual I found the example:

#include "udf.h"

#define CON 20.0

DEFINE_SOURCE(cell_x_source, cell, thread, dS, eqn)
{
  real source;

  if (C_T(cell,thread) <= 288.)
  {
    source = -CON*C_U(cell,thread);
    dS[eqn] = -CON;
  }
  else
  {
    source = dS[eqn] = 0.;
  }

  return source;
}

and I understand how it works, but how can I access the inviscid flux in FLUENT? And how would the derivate dS[eqn] look like for this case?

About any help I would be very grateful!


Many greets

Mohinder