How to compute the integral of a variable along a line using UDF in FLUENT?
Line surfaces in FLUENT have area in 2d and 2daxi but not in 3d, so currently integrals on line surfaces are zero.
This solution provides a sample User Defined Function (UDF) that shows how to calculate the integral of a variable (say temperature) along a line in fluent, then divide it by the length of the line, and output this number at every time step, for a large number of lines in a transient 3d simulation. /*sample UDF code to work with for line integration: you will need to augment it to use it with multiple lines, etc*/ /* UDF to calculate line integral and line lenght weighted average of a variable. by sra and fm 5/05 To have line integral and line average of variable V on a line S u need to Go to Contour panel and do compute of V on S This will fill the values necessary for computation. */ #include "udf.h" #include "surf.h" #include "cxsurf.h" #include "cxiface.h" #define SID 14 DEFINE_ON_DEMAND(line_integrals) { Surface *s; int i, j, k, index; float v[2][3]; float face_val, line_int, line_ave, curve_length, tot_length; /* s is the surface having data of surface id SID s>np no of points. s>nf no of faces. s>points points to array of np points. s>ip pointer to interpolation vector of np points. s>facets facet list of the surface. length of array is nfl. s>cells pointer to cells of each point. s>nfl facet list length. */ s = SurfaceList+SID; tot_length = 0.; line_int = 0.; for (i = 0, k = 0; k < s>np; i +=s>facets[i]+1, k++) { if (s>facets[i] != 2) continue; /* makes sure only line is there */ /*facet value for line surfaces will be the cell value in which the facet lies*/ face_val = SV_Cell_Field(SVFIELD_CURRENT, s>cells+k); for (j = 0; j < s>facets[i]; j++) { index = s>facets[i+j+1]; /* sets point coord x y z*/ SET_POINT_COORD(v[j][0], s>points[index], X_DIM); SET_POINT_COORD(v[j][1], s>points[index], Y_DIM); SET_POINT_COORD(v[j][2], s>points[index], Z_DIM); } /* Length calculation */ curve_length = sqrt(SQR(v[0][0]v[1][0])+SQR(v[0][1]v[1][1])+SQR(v[0][2]v[1][2])); tot_length += curve_length; line_int += curve_length * face_val; } line_ave = line_int/tot_length; Message("Line integral of currently displayed variable on line of ID %d : %6fn", SID, line_int); Message("Line length average of currently displayed variable on line of ID %d : %6fn", SID, line_ave); } 

