[ankit956]

I added one UDF for finding the particle concentration inside the office room. But after initializing some for 9-10 steps, the fluent is crashing down. Could anyone please help me what is wrong with my UDF. I'm posting my UDF here:


#include "udf.h"

/* Define the region boundaries */
#define X_MIN -3.75
#define X_MAX 3.75
#define Y_MIN 0.0
#define Y_MAX 3.5
#define Z_MIN 0.0
#define Z_MAX 5.5

DEFINE_EXECUTE_AT_END(particle_concentration)
{
Domain *domain;
Thread *ct;
cell_t c;

real volume = 0.0;
int particle_count = 0;
real pc; // Change pc to real type for proper division

/* Get the domain pointer */
domain = Get_Domain(1);
if (domain == NULL)
{
printf("Error: Domain is NULL\n");
return;
}

/* Loop over all cell threads in the domain */
thread_loop_c(ct, domain)
{
/* Loop over all cells in the thread */
begin_c_loop(c, ct)
{
/* Get the cell centroid coordinates */
real xc[ND_ND];
C_CENTROID(xc, c, ct);

/* Check if the cell is within the specified region */
if (xc[0] >= X_MIN && xc[0] <= X_MAX &&
xc[1] >= Y_MIN && xc[1] <= Y_MAX &&
xc[2] >= Z_MIN && xc[2] <= Z_MAX)
{
/* Update the volume of the region */
volume += C_VOLUME(c, ct);
}
}
end_c_loop(c, ct)
}

/* Calculate the particle concentration */
if (volume > 0)
{
pc = (real)particle_count / volume;
}
else
{
printf("Warning: Volume is zero\n");
pc = 0.0;
}

/* Output the result to the console */
printf("Particle concentration in the defined volume: %g particles per unit volume\n", pc);
}