[Jenny_Wong]

Hi Everyone,

I've written a UDF code to time-average the mean mixture fraction values computed in FLUENT after each iteration.

However, I received "SEMENTATION VIOLATION" error. I think it is becuase of C_FMEAN. When I use other macro, the averaging code works fine.


thread_loop_c(ct,d)
{
begin_c_loop(c,ct)
{
C_UDMI(c,ct,0) = C_FMEAN(c,ct)+C_UDMI(c,ct,0);
C_UDMI(c,ct,1) = C_UDMI(c,ct,0)/count ;
}
end_c_loop(c,ct)
}

If anyone knows the solution to this problem, please kindly advise me. I will really appreciate it. Thank you very much.
Jenny

[ComputerGuy]

Can you post the rest of the code? You could always test the code by replacing C_FMEAN with a constant (say, 0.0) to see if that's causing the error.

Do you have the premixed combustion model turned on?

ComputerGuy

[Jenny_Wong]

Dear ComputerGuy,
Thank you for your reply.

The premixed combustion model wasn't turned on. I used the non-premixed model.

I tried replacing C_FMEAN with a constant as you suggested and the code worked perfectly.

Thanks again,
Jenny

Here's my code:
#include "udf.h"
int count = 0;
real start_avg = 0.001;
DEFINE_EXECUTE_AT_END(avg)
{
if(CURRENT_TIME > start_avg)
{
Domain *d;
Thread *ct;
cell_t c;
d=Get_Domain(1);
count = count+1;
thread_loop_c(ct,d)
{
begin_c_loop(c,ct)
{
C_UDMI(c,ct,0) = C_FMEAN+C_UDMI(c,ct,0);
C_UDMI(c,ct,1) = C_UDMI(c,ct,udm_offset)/count ;
}
end_c_loop(c,ct)
}
}
}

[ComputerGuy]

The code should be something similar to:

#include "udf.h"
#include "sg_mem.h"

int count = 0;
real start_avg = 0.001;
DEFINE_EXECUTE_AT_END(avg)
{
if(CURRENT_TIME > start_avg)
    {
    Domain *d;
    Thread *ct;
    cell_t c;
    d=Get_Domain(1);
    count = count+1;
    thread_loop_c(ct,d)
        {
        begin_c_loop(c,ct)
            {
            C_UDMI(c,ct,0)  = C_FMEAN(c,ct)+C_UDMI(c,ct,0);                     
            C_UDMI(c,ct,1)  = C_UDMI(c,ct,udm_offset)/count ;
            }
              end_c_loop(c,ct)
        }
    }
}

Try those two changes. You were using C_FMEAN incorrectly and not passing c and t, and I think the first include statement should help. I also don't see any place where you're defining udm_offset. Make sure that if you're using variables, they're defined. I don't know if this is actually the code you're using, but I'm surprised it was able to compile.

Let me know

ComputerGuy

Originally Posted by Jenny_Wong

Dear ComputerGuy,
Thank you for your reply.

The premixed combustion model wasn't turned on. I used the non-premixed model.

I tried replacing C_FMEAN with a constant as you suggested and the code worked perfectly.

Thanks again,
Jenny

Here's my code:
#include "udf.h"
int count = 0;
real start_avg = 0.001;
DEFINE_EXECUTE_AT_END(avg)
{
if(CURRENT_TIME > start_avg)
{
Domain *d;
Thread *ct;
cell_t c;
d=Get_Domain(1);
count = count+1;
thread_loop_c(ct,d)
{
begin_c_loop(c,ct)
{
C_UDMI(c,ct,0) = C_FMEAN+C_UDMI(c,ct,0);
C_UDMI(c,ct,1) = C_UDMI(c,ct,udm_offset)/count ;
}
end_c_loop(c,ct)
}
}
}

[Jenny_Wong]

Dear ComputerGuy,

I'm sorry I made a mistake. I did use C_FMEAN(c,ct) in the code and replace udm_offset with 0.

I have added include "sg_mem.h" as you recommended but still FLUENT gave the same SEGMENTATION VIOLATION error.

I saw in 1 forum, C_FMEAN was used like this
fmean = C_FMEAN(c,ct);
a = fmean + b;
So, I tried this in my code too but it didn't work either.

Now, I'm trying to just print the C_FMEAN values and see if it reads the values alone. Please kindly advise if you know what is wrong. Thank you again for your help.

Jenny

[ComputerGuy]

Try the code below...The user-defined memory is valid everywhere, but the C_FMEAN is likely to be only valid if the cell thread is a fluid thread.

#include "udf.h"
#include "sg_mem.h"

int count = 0;
real start_avg = 0.001;
DEFINE_EXECUTE_AT_END(avg)
{
if(CURRENT_TIME > start_avg)
    {
    Domain *d;
    Thread *ct;
    cell_t c;
    d=Get_Domain(1);
    count = count+1;
    thread_loop_c(ct,d)
        {
        if (FLUID_THREAD_P(ct))
        {
        begin_c_loop(c,ct)
            {
            C_UDMI(c,ct,0)  = C_FMEAN(c,ct)+C_UDMI(c,ct,0);                     
            C_UDMI(c,ct,1)  = C_UDMI(c,ct,udm_offset)/count ;
            }
              end_c_loop(c,ct)
          }
        }
    }
}

ComputerGuy

Originally Posted by Jenny_Wong

Dear ComputerGuy,

I'm sorry I made a mistake. I did use C_FMEAN(c,ct) in the code and replace udm_offset with 0.

I have added include "sg_mem.h" as you recommended but still FLUENT gave the same SEGMENTATION VIOLATION error.

I saw in 1 forum, C_FMEAN was used like this
fmean = C_FMEAN(c,ct);
a = fmean + b;
So, I tried this in my code too but it didn't work either.

Now, I'm trying to just print the C_FMEAN values and see if it reads the values alone. Please kindly advise if you know what is wrong. Thank you again for your help.

Jenny

[Jenny_Wong]

Dear ComputerGuy,
You are right. The code works very well now with no SEGMENTATION VIOLATION error. I have been struggling with this problem for 2 months and now I am very happy. Thank you very much,
Jenny

[ComputerGuy]

Outstanding. Glad to help!

Originally Posted by Jenny_Wong

Dear ComputerGuy,
You are right. The code works very well now with no SEGMENTATION VIOLATION error. I have been struggling with this problem for 2 months and now I am very happy. Thank you very much,
Jenny

[hongfu2233]

could you tell me some usefulness about the head file "sg_mem.h", many thanks!