Forcing pressure outlet according to pressure prob

Cedric

Hi!

I have to force a pressure outlet with the data coming from a probe that is somewhere else is the fluid domain. I want the pressure in the outlet to be the same as the one in a specific surface somewhere else in the domain.

This is what I have tried, but it doesn't work. I can't even manage to get any average pressure on the surface I have selected according to its ID. Could you please help me?

DEFINE_PROFILE(tube_pressure,thread,nv) {

int ID;

int domain_id;

real p;

real ptot;

ID = 22; /* ID ot the probe*/

domain_id =1.0;

p = 0.0;

ptot = 0.0;

domainp = Get_Domain(domain_id);

/ f_threadp = Lookup_Thread(domainp, ID);

begin_f_loop(fp, f_threadp) /* loops over faces in a face thread */

{ i = (i + 1.0);

p += F_P(fp,f_threadp);

Message("i = %f, p = %f\n",i ,p);

}

end_f_loop(fp, f_threadp);

ptot = (p / i);

Message("Hello world ID= %f\n", ID);

Message ("index = %f, ptot = %f\n", index, ptot);

/* definition of pressure profile */

begin_f_loop (f, thread)

{

F_CENTROID(x, f, thread);

if (x[1] >= top) /*above=never*/

{

F_PROFILE(f, thread, nv) = 0.;

}

else /*below=always*/

{

F_PROFILE(f, thread, nv) = ptot;

}

}

end_f_loop (f,thread) }