[mohammadreza_hj]

Hi everybody,

I have a question that I'm stuck in it; I asked this question at the end of the last thread that I made, but because I'm in a hurry (deadlines) and I thought that because the discussion's getting long, nobody answered. So I ask it again here.

I want to catch the VOF of the phases on the wall by a defined profile. Is there any macro for this purpose, and with what data accessing and looping macro?

I send my code here, so if you can help me, I will be appreciated. but when I apply my code, the fluent crashes out as soon as I initialize the solution. (image attache)


#include "udf.h"
#include "mem.h"
#include "sg.h"
#include "sg_mphase.h"
#include "flow.h"
#include "metric.h"
DEFINE_PROFILE(heatflux_slip_shoulder_W, t, i)
{
/*Domain *d=Get_Domain(1);*/
int phase_domain_index = 0; /* primary phase index is 0 */
Thread *tm = THREAD_SUPER_THREAD(t);
Thread *subthread = THREAD_SUB_THREAD(tm,phase_domain_index);
double p[ND_ND]; /* this will hold the position vector */
double x, y, r, qslip, sigma, tav_stick, tav_fric, temp, Vx, Vy, Vz, V_W, V_T, VF;
double w = 1120;
double pi = 3.1415927;
double U = 0.002;
double press = 50000000;
double delta = 0.65;
double etta = 0.7;
double heat_ratio = 0.6383;
double Rshol = 0.0075;
double fric ;
double Rpin = 0.0025;
face_t f;
/*thread_loop_f(f, d)*/
{
begin_f_loop(f,t)
{
VF =F_VOF(f,subthread);
F_CENTROID(p, f, t);
x = p[0];
y = p[1];
r = sqrt((x*x) + (y*y));
temp = F_T(f,t);if ((297.<=temp)&&(temp<=273.))
{
sigma_6 = 325.80263157895 + (-0.171052631578952 * temp);
sigma_5 = 282.355263157897 + (-0.0921052631578974*temp);
}
...(and some conditions like above)...

if ((297.<=temp)&&(temp<=855.))
fric = 0.383752244165171 + (-0.000448833034111311 * temp);
else
fric = 0;

tav_6 = sigma_6*1000000/1.732;
tav_5 = sigma_5*1000000/1.732;
tav = VF * tav_6 + (1 - VF ) * tav_5;
qslip = ((pi*w*r/30) - (U*y/r)) * ((etta * (1 - delta) * tav) + (delta * fric * press));

F_PROFILE(f, t, i) = heat_ratio * qslip;
}
end_f_loop(f,t)
}

[AlexanderZ]

compile this code

#include "udf.h"
#include "mem.h"
#include "sg.h"
#include "sg_mphase.h"
#include "flow.h"
#include "metric.h"
DEFINE_PROFILE(heatflux_slip_shoulder_W, t, i)
{
/*Domain *d=Get_Domain(1);*/
int phase_domain_index = 0; /* primary phase index is 0 */
Thread *subthread = THREAD_SUB_THREAD(t,phase_domain_index);
double p_coord[ND_ND]; /* this will hold the position vector */
double x, y, r, qslip, sigma, tav_stick, tav_fric, temp, Vx, Vy, Vz, V_W, V_T, VF;
double w = 1120;
double pi = 3.1415927;
double U = 0.002;
double press = 50000000;
double delta = 0.65;
double etta = 0.7;
double heat_ratio = 0.6383;
double Rshol = 0.0075;
double fric ;
double Rpin = 0.0025;
face_t f;

begin_f_loop(f,t)
{
VF =F_VOF(f,subthread);
F_CENTROID(p_coord, f, t);
x = p_coord[0];
y = p_coord[1];
r = sqrt((x*x) + (y*y));
temp = F_T(f,t);if ((297.<=temp)&&(temp<=273.))
{
sigma_6 = 325.80263157895 + (-0.171052631578952 * temp);
sigma_5 = 282.355263157897 + (-0.0921052631578974*temp);
}
/*...(and some conditions like above)...*/

if ((297.<=temp)&&(temp<=855.))
fric = 0.383752244165171 + (-0.000448833034111311 * temp);
else
fric = 0;

tav_6 = sigma_6*1000000/1.732;
tav_5 = sigma_5*1000000/1.732;
tav = VF * tav_6 + (1 - VF ) * tav_5;
qslip = ((pi*w*r/30) - (U*y/r)) * ((etta * (1 - delta) * tav) + (delta * fric * press));

F_PROFILE(f, t, i) = heat_ratio * qslip;
}
end_f_loop(f,t)
}

[mohammadreza_hj]

Originally Posted by AlexanderZ

compile this code

#include "udf.h"
#include "mem.h"
#include "sg.h"
#include "sg_mphase.h"
#include "flow.h"
#include "metric.h"
DEFINE_PROFILE(heatflux_slip_shoulder_W, t, i)
{
/*Domain *d=Get_Domain(1);*/
int phase_domain_index = 0; /* primary phase index is 0 */
Thread *subthread = THREAD_SUB_THREAD(t,phase_domain_index);
double p_coord[ND_ND]; /* this will hold the position vector */
double x, y, r, qslip, sigma, tav_stick, tav_fric, temp, Vx, Vy, Vz, V_W, V_T, VF;
double w = 1120;
double pi = 3.1415927;
double U = 0.002;
double press = 50000000;
double delta = 0.65;
double etta = 0.7;
double heat_ratio = 0.6383;
double Rshol = 0.0075;
double fric ;
double Rpin = 0.0025;
face_t f;

begin_f_loop(f,t)
{
VF =F_VOF(f,subthread);
F_CENTROID(p_coord, f, t);
x = p_coord[0];
y = p_coord[1];
r = sqrt((x*x) + (y*y));
temp = F_T(f,t);if ((297.<=temp)&&(temp<=273.))
{
sigma_6 = 325.80263157895 + (-0.171052631578952 * temp);
sigma_5 = 282.355263157897 + (-0.0921052631578974*temp);
}
/*...(and some conditions like above)...*/

if ((297.<=temp)&&(temp<=855.))
fric = 0.383752244165171 + (-0.000448833034111311 * temp);
else
fric = 0;

tav_6 = sigma_6*1000000/1.732;
tav_5 = sigma_5*1000000/1.732;
tav = VF * tav_6 + (1 - VF ) * tav_5;
qslip = ((pi*w*r/30) - (U*y/r)) * ((etta * (1 - delta) * tav) + (delta * fric * press));

F_PROFILE(f, t, i) = heat_ratio * qslip;
}
end_f_loop(f,t)
}

Hello AlexanderZ,
Thank you so much,

I do your command, but nothing changed; I attached the error log file also here.
and I should mention that I put my simulation on the volume of fluid multiphase. (not mixture or Elureian)

So, what do you think? Where is the problem? It's from F_VOF macro, or it's from its domain and thread?

[mohammadreza_hj]

Originally Posted by mohammadreza_hj

Hello AlexanderZ,
Thank you so much,

I do your command, but nothing changed; I attached the error log file also here.
and I should mention that I put my simulation on the volume of fluid multiphase. (not mixture or Elureian)

So, what do you think? Where is the problem? It's from F_VOF macro, or it's from its domain and thread?

I'm so sorry; I forgot to attach the error file.
so here is :

[AlexanderZ]

try this code

#include "udf.h"

DEFINE_PROFILE(heatflux_slip_shoulder_W, t, i)
{
/*Domain *d=Get_Domain(1);*/
int phase_domain_index = 0; /* primary phase index is 0 */
Thread *subthread = THREAD_SUB_THREAD(t,phase_domain_index);
face_t f;

begin_f_loop(f,t)
{
VF =F_VOF(f,subthread);
F_UDMI(f,t,0) = VF;
F_PROFILE(f, t, i) = 1.0;
}
end_f_loop(f,t)
}

allocate 1 user-defined memory location in fluent GUI.
now you may plot VOF distribution

[mohammadreza_hj]

Originally Posted by AlexanderZ

try this code

#include "udf.h"

DEFINE_PROFILE(heatflux_slip_shoulder_W, t, i)
{
/*Domain *d=Get_Domain(1);*/
int phase_domain_index = 0; /* primary phase index is 0 */
Thread *subthread = THREAD_SUB_THREAD(t,phase_domain_index);
face_t f;

begin_f_loop(f,t)
{
VF =F_VOF(f,subthread);
F_UDMI(f,t,0) = VF;
F_PROFILE(f, t, i) = 1.0;
}
end_f_loop(f,t)
}

Allocate 1 user-defined memory location in fluent GUI.
now you may plot VOF distribution

Thanks for your attention, Mr. alexander

Well, I try this one, and again fluent crashes out.
I'm confused, and I don't know what's happening?
I am suspicious that maybe I do something wrong.
So I want to say what I did; first of all, I compiled the UDF, it was ok.
Then I allocated one udm.
Then I defined my model, multiphase vof (implicit), and energy and laminar flow.
Then I defined my material from fluent material, air, and water vapor.
Then in the list of phases, I made air primary and the other secondary phase.
And in boundary, I define velocity inlet of phases and outflow and apply udf on wall shoulder.
Then initialize and crash.
Here my image of geometry. This image is when I initialize and patch it first and then apply the UDF on the wall.

Is there any problem with my steps?

[mohammadreza_hj]

Originally Posted by AlexanderZ

try this code

#include "udf.h"

DEFINE_PROFILE(heatflux_slip_shoulder_W, t, i)
{
/*Domain *d=Get_Domain(1);*/
int phase_domain_index = 0; /* primary phase index is 0 */
Thread *subthread = THREAD_SUB_THREAD(t,phase_domain_index);
face_t f;

begin_f_loop(f,t)
{
VF =F_VOF(f,subthread);
F_UDMI(f,t,0) = VF;
F_PROFILE(f, t, i) = 1.0;
}
end_f_loop(f,t)
}

Allocate 1 user-defined memory location in fluent GUI.
now you may plot VOF distribution

Hi, again Mr alexanderz,
I want to share some result with you,
well, I explorer in the UDF manual, and I find one example that was near my project, and I apply that to my work. The example was: "Example UDF that Utilizes UDM and UDS Variables" on page 332 in UDF manual 19 R3.
When I apply to adjust, fluent crashes out again.
What do you think about that? Where is the problem?
and I should mention that I allocated udm and uds in GUI.
Thanks a lot.

[AlexanderZ]

most likely, if you use code from manual and get error - there is a problem with your case, not with a code

you've told

This image is when I initialize and patch it first and then apply the UDF on the wall.

try this approach, actually, initialize first, then apply UDF and run the case.

[mohammadreza_hj]

Originally Posted by AlexanderZ

most likely, if you use code from manual and get error - there is a problem with your case, not with a code

you've told


try this approach, actually, initialize first, then apply UDF and run the case.

Thank you alexanderz

Can you tell me more about the case problem, because when I don't apply the VOF code, it runs very well?

And about the initialize first and then apply UDF, I tried that. When I want to apply the macro on the wall, it crashes out.

[mohammadreza_hj]

Originally Posted by mohammadreza_hj

Thank you alexanderz

Can you tell me more about the case problem, because when I don't apply the VOF code, it runs very well?

And about the initialize first and then apply UDF, I tried that. When I want to apply the macro on the wall, it crashes out.

Oh, I'm sorry I did not get your point; I think I was busy; I will do that.
Thanks.

[mohammadreza_hj]

Originally Posted by AlexanderZ

most likely, if you use code from manual and get error - there is a problem with your case, not with a code

you've told


try this approach, actually, initialize first, then apply UDF and run the case.

Mr. Alexanderz,

I try what you said. I initialize first and then apply the UDF and put them on the boundary; up to now, there was no problem.
Then I click on the calculation, and I face this error.