[schwaral]

Hi,
I figured out, that my Problem is because c_density is wrong adressed. without this line the code works fine. My problem is, that I don't know how to tell the programm, that I need the density of the cell which is attached to the face at the boundary. DEFINE_PROFILE and the begin_f_loop do not provide the cell index I guess. I tried to make a c_loop around the c_density but this didn't help.

Mybe is there a solution with pointers so you can adress the face and the cell at the same time? Unforunatley with my C++ skills I'm a bit lost here. Any Help would be nice.

DEFINE_PROFILE(massfraction_for_const_RH,thread,po sition)
{
cell_t cell;
face_t f;

double a= -5800.2206;
double T,massfraction_new;
double p_sat, c_density;

begin_f_loop(f,thread)
{
T=F_T(f,thread);

p_sat= exp(a*pow(T,-1));

c_density = C_R(cell, thread);

massfraction_new=0.7*p_sat/c_density;
F_PROFILE(f,thread,position)=massfraction_new;
end_f_loop(f,thread)
}
}

[torjo]

I think you need to use

F_C0(f,thread)

and

THREAD_T0(thread)

to get the cell and thread for your call to C_R
try:

c_density = C_R(F_C0(f,thread),THREAD_T0(thread))

BR.
Torjo

[schwaral]

yep that was the solution thanks a lot!

[p.norouzi.gv]

Originally Posted by schwaral

Hi,
I figured out, that my Problem is because c_density is wrong adressed. without this line the code works fine. My problem is, that I don't know how to tell the programm, that I need the density of the cell which is attached to the face at the boundary. DEFINE_PROFILE and the begin_f_loop do not provide the cell index I guess. I tried to make a c_loop around the c_density but this didn't help.

Mybe is there a solution with pointers so you can adress the face and the cell at the same time? Unforunatley with my C++ skills I'm a bit lost here. Any Help would be nice.

DEFINE_PROFILE(massfraction_for_const_RH,thread,po sition)
{
cell_t cell;
face_t f;

double a= -5800.2206;
double T,massfraction_new;
double p_sat, c_density;

begin_f_loop(f,thread)
{
T=F_T(f,thread);

p_sat= exp(a*pow(T,-1));

c_density = C_R(cell, thread);

massfraction_new=0.7*p_sat/c_density;
F_PROFILE(f,thread,position)=massfraction_new;
end_f_loop(f,thread)
}
}

hi , I found out you solved the problem.
Could you please share the correct answer ?

[AlexanderZ]

#include "udf.h"
DEFINE_PROFILE(massfraction_for_const_RH,thread,po sition)
{
cell_t cell,c0;
face_t f;
thread *t0;
real a= -5800.2206;
real T,massfraction_new;
real p_sat, c_density;

begin_f_loop(f,thread)
{
T=F_T(f,thread);
c0=F_C0(f,thread)
t0=THREAD_T0(thread)
p_sat= exp(a*pow(T,-1));
c_density = C_R(c0,t0);
massfraction_new=0.7*p_sat/c_density;
F_PROFILE(f,thread,position)=massfraction_new;
end_f_loop(f,thread)
}
}