[SNeCu]

I wirte UDF to get momentum flux in x direction on a face after each iteration. However, it shows error: receive a fatal signal(Segmentation fault). Error Object: #f and it cannot continue. Could any one tell me where is the mistake(k-e, series)? Thanks.
#include "udf.h"
DEFINE_EXECUTE_AT_END(momentumxflux)
{
real mf=0;
Domain *d;
Thread *t;
face_t f;
thread_loop_f(t,d)
{
begin_f_loop(f,t)
mf+=F_FLUX(f,t)*F_U(f,t);
end_f_loop(f,t)
}
printf("Momentum flux on this face is: %g\n", mf);
}

[Bruno Machado]

Originally Posted by SNeCu

I wirte UDF to get momentum flux in x direction on a face after each iteration. However, it shows error: receive a fatal signal(Segmentation fault). Error Object: #f and it cannot continue. Could any one tell me where is the mistake(k-e, series)? Thanks.
#include "udf.h"
DEFINE_EXECUTE_AT_END(momentumxflux)
{
real mf=0;
Domain *d;
Thread *t;
face_t f;
thread_loop_f(t,d)
{
begin_f_loop(f,t)
mf+=F_FLUX(f,t)*F_U(f,t);
end_f_loop(f,t)
}
printf("Momentum flux on this face is: %g\n", mf);
}

Since you are not hooking it in any boundary, I think you gotta define the value of the thread. Have a look in your case file and change your code in order to match this new face thread value.

[pakk]

Originally Posted by Bruno Machado

Since you are not hooking it in any boundary, I think you gotta define the value of the thread. Have a look in your case file and change your code in order to match this new face thread value.

Yes. This.

Two relevant lines from the code:

Domain *d;

The line above means: reserve variable "d" to contain a pointer a domain.

thread_loop_f(t,d)

This means: loop through all threads in domain "d". But what is domain "d"? You didn't specify it yet! This is causing your program to crash.

And actually, you don't want to loop through all threads. You want to have only one thread: the face that you are interested in. You need to find a way to specify this face, look in the help for this.

[SNeCu]

Yes, I revised the code and it works! Thank you.