[homay95]

so I have written a code for " Tf = Ts + ks(dT/dn)" in which Tf is fluid temperature, Ts is the wall temperature and ks is temperature jump coefficient.

DEFINE_PROFILE(temperaturejump, thread, position)
{
face_t f;
cell_t c;
Thread *tc;
double b = 0.000001;
double dTdy ,tem ,walltemp;
begin_f_loop(f, thread)
{
c = F_C0(f,thread);
tc = THREAD_T0(thread);
dTdy = C_T_G(c,tc)[1];
walltemp = F_T(f, thread);
tem = (b*dTdy) + walltemp;
F_PROFILE(f, thread, position) = tem;
}
end_f_loop(f, thread)
}

is it correct or does it have any problem??

[vinerm]

The code is correct, but missing the commands that I shared earlier for storing temperature gradient.

[homay95]

so should i add them within the loop?
and another question is that does the expression "F_T(f, thread)" gives wall temperature? or do i have to use "C_T(c, thread)" instead?

[vinerm]

You need to add those three lines outside the loop.

F_T is good. C_T gives cell temperature.

[homay95]

okay. but im having a trouble with hooking the UDFs. where do i have to hook this temperature jump? for the face (adjacent to solid) or for the shadow(adjacent to fluid) ??

[vinerm]

The profiles are same for a wall and its shadow. If you hook it to one, it automatically works for the other. And that's where you have trouble since the jump should exist only on the fluid side. So, DEFINE_PROFILE may not be a good idea. Use DEFINE_ADJUST, access the thread on the fluid size, and then change the temperature. Let the coupled condition be as it is so that thermal energy conservation could be maintained.

[homay95]

aha I think I understand what you are saying, do you mean if I use DEFINE_ADJUST, it would implement the heat transfer(conduction) between the face and its shadow by itself too, right?

[vinerm]

That's what a default coupled condition does. In reality, there is no conduction between a boundary and its shadow. They are same. When they are coupled, the heat flux entering one is equal to heat flux exiting other, i.e., thermal energy is conserved.

[homay95]

but the problem Im having here is that I need to apply both conditions to the liquid-solid interface, I mean I need to apply temperature jump to shadow(adjacent to fluid) temperature and also I need to have a coupled condition between face and shadow. but in fluent I can only either apply temperature jump or just have coupled condition, not both of them at a time. what should I do?

[vinerm]

That's why you need to use DEFINE_ADJUST. You can apply jump using UDF and coupled condition will maintain the conservation across shadow pairs.

[homay95]

okay. now I understand. thank you very much for your help

[homay95]

so I have read about DEFINE_ADJUST macro. but I dont know how does it work!? what are the differences between this and DEFINE_PROFILE? which arguments should I use ? and to which part of FLUENT should I apply that?

[homay95]

You mean first I have to write a DEFINE_ADJUST which includes temperature jump function in it and hook it to function hooks and then leave the thermal boundary conditions for faces and their shadows as coupled. right?

[homay95]

Hi again
I have written a code by DEFINE_ADJUST for temperature jump as follows :

DEFINE_ADJUST(temperaturejump3, domain)
{
Domain *d;
int zone_ID = 16;
Thread *f_thread = Lookup_Thread(domain, zone_ID);
Thread *c_thread;
cell_t cell;
face_t face;
double dTdy, temp, walltemp;
double b = 0.000001;
thread_loop_f (f_thread,domain)
{
begin_f_loop(face,f_thread)
{
cell = F_C0(face,f_thread);
c_thread = THREAD_T0(f_thread);
walltemp = F_T(face,f_thread);
dTdy = C_T_G(cell,c_thread)[1];
F_PROFILE(face, f_thread, zone_ID) = (b*dTdy) + walltemp;
}
end_f_loop(face,f_thread)
}
}

is it correct or does it need any modifying??

Thanks in advance