[ymd]

Hi, I am performing a heat loss simulation of skin tissue. I have to apply convection, radiation, and evaporation boundary conditions to the outer skin surface. For evaporation, I have heat flux value. I was able to apply convection and radiation boundary values through mixed boundary condition. How can I add heat flux to the same surface? Is there a UDF macro for such an additional boundary condition. Help will be much appreciated.

[vinerm]

You cannot do that directly. However, you can implement the known heat flux as a heat source.

[ymd]

Thank you Vinerm. Should the heat flux be provided by giving wall thickness and heat generation? Or, is it through UDF using DEFINE_SOURCE?

[vinerm]

Both should work fine. Try first with thin wall. If it doesn't work as expected, then you can apply cell based source terms

[ymd]

Thank you. I got it working.

[GreyMat]

ymd, I am new in C language, can you please guide me how to add convection and heat flux boundary condition on same wall. thank you

Can you please provide the C code for heat flux.

[AlexanderZ]

put heat generation rate and wall thickness in convection tab, which will represent your heat flux

[GreyMat]

Thanks AlexanderZ, I am aware about wall thickness and heat generation.

But can you provide the C language (UDF) for adding heat flux at a particular wall?

[AlexanderZ]

I don't know the way, how you can apply heat flux and convection directly
UDF allows you to modifiy heat source according to your needs, but you still should hook it, and ther is no other option accept heat source and wall thickness (at least I don't know)

[GreyMat]

Ok, thanks mate

[ymd]

One way to do it is by using DEFINE_PROFILE macro. Here is an example code:

#include "udf.h"

#define heat_flux 40 //Unit W/m2
#define convection 10 //Unit W/m2 K
#define Temp 298.15 //Unit K (25 deg C)
DEFINE_PROFILE(wallbc,thread,i)
{

face_t f;

begin_f_loop(f,thread)
F_PROFILE(f,thread,i) = -heat_flux + convection*(Temp-F_T(f,thread));
end_f_loop(f,thread)
}
Hope it works.

[GreyMat]

Ok thanks ymd, ill try and reply

[GreyMat]

(AlexanderZ and ymd)
-I have attached a UDF which define heat flux and convection boundary condition at same wall
- But ANSYS give me SIGSEGV error and fluent suddenly closed.
- Can you please provide me any solution or changes within the UDF, thank you
-please check the attached files

[AlexanderZ]

change the name of the function
DEFINE_PROFILE(name_of_the_function,thread,i)

it must be different from names of variables, you are using

[guti]

I want to apply heat flux on rotating tube with convection and radiation. Already, I have created my first UDF for stationary tube and faced big problem with rotating one. One thing, I do not understand how to read the face wall temperature and other I could not able to locate point where the heat flux gonna to be applied. Can I combine DEFINE_ZONE_MOTION and DEFINE_PROFILE macros or which one can I use for the problem? Can anyone help me? Please tell me how can I proceed it? Thank you for your cooperation.

[AlexanderZ]

EFINE_ZONE_MOTION and DEFINE_PROFILE macros could be used inside one UDF
they are independent from each other

[guti]

Ok thanks for your reply.

So, Can I do like this?

#include"udf.h"
#include"mem.h"
#define sigma 5.6697*10^-8
#define epsilon 0.1
#define H_air 8.52 /* convective heat transfer coefficient*/
#define T_air 298
#define T_sky 278
#define R 0.0147 /*internal radius of the tube*/
DEFINE_ZONE_MOTION(rotation,axis,origin,vel,omega, time,dtime)
{

vel[0]=0.0;/* reset velocities */
vel[1]=0.0;
vel[2]=0.0025;

omega[0]=0.0;
omega[1]=0.0;
omega[2]=0.52;

N3V_D (axis, =, 3.0,3.0,3.0);
N3V_S (origin, =,2.0);
}

DEFINE_PROFILE(heat_flux,thread,i)
real x[ND_ND];/* this will hold the position vector*/
real y;
real z;
real r;

begin_f_loop (f,t)
{
c0=F_C0(f,t);
t0=THREAD_T0(t);
T_abs=F_T(f,t);

F_CENTROID(x,f,t);
y=x[1];
z=x[2];
r=z/sqrt((z*z)+(y*y));

F_PROFILE(f,t,i)= (1408*r)-(F_T(f,t)-T_air)-(0.1*0.000000056697)*(pow(F_T(f,t),4)-pow(T_sky,4));
}

end_f_loop(f,t)
}

or how can I just proceed? please show me the way.

[AlexanderZ]

compile code, read log, fix error if any

then you need to apply functions to corresponding faces
select surface of your interest -> boundary conditions -> heat flux -> select your udf

for more information read ansys fluent customization manual