[behest]

Hello all,
I would like to calculate this k.w (k is gradient of the kinetic energy and w is the gradient of the specific dissipation rate) for each cell. Can you tell me how should write it in UDF?
I used several ways to compute the dot product multipler like as:
#include "udf.h"
#include "math.h"

DEFINE_ON_DEMAND(grad)
{
Domain *domain=Get_Domain(ROOT_DOMAIN_ID);

cell_t c;
Thread *t;
real s7;

thread_loop_c(t,domain) /*loops over all cell threads in domain*/
{
begin_c_loop(c,t) /* loops over cells in a cell thread */
{
s7=C_K_G(c,t)[0]*C_O_G(c,t)[0]+C_K_G(c,t)[1]*C_O_G(c,t)[1]+C_K_G(c,t)[2]*C_O_G(c,t)[2];
Message("Gradient of k and w in a cell are: (%f) \n", s7);
}
end_c_loop(c,t)
}
}

[behest]

The problem is solved now. Actually, it needs to active solve/se/expert and let to software to save all gradient data.
Thanks all

Originally Posted by behest

Hello all,
I would like to calculate this k.w (k is gradient of the kinetic energy and w is the gradient of the specific dissipation rate) for each cell. Can you tell me how should write it in UDF?
I used several ways to compute the dot product multipler like as:
#include "udf.h"
#include "math.h"

DEFINE_ON_DEMAND(grad)
{
Domain *domain=Get_Domain(ROOT_DOMAIN_ID);

cell_t c;
Thread *t;
real s7;

thread_loop_c(t,domain) /*loops over all cell threads in domain*/
{
begin_c_loop(c,t) /* loops over cells in a cell thread */
{
s7=C_K_G(c,t)[0]*C_O_G(c,t)[0]+C_K_G(c,t)[1]*C_O_G(c,t)[1]+C_K_G(c,t)[2]*C_O_G(c,t)[2];
Message("Gradient of k and w in a cell are: (%f) \n", s7);
}
end_c_loop(c,t)
}
}