Dear users,

I am working on a Natural convection problem via freeFEM?

A researcher has plotted the results of Normal(grad(T)) along a circle see the top right figure.

what will be the command of Normal(grad(T)) along a Circle in FreeFEM?

Detail is attached in the image below.

just do

```
macro grad(u) [dx(u),dy(u)]//
real flux = int1d(Th,circle)( grad(T)'*[N.x,N.y]) ;
```

Thanks for your kind guidance. It works fine for entire length but it gives only one values.

How do i plot along the running circle for entire running flux. mean x-axis (0-2pi) and y-axis (flux calculation from 0 to 2pi).

Check this command please .

see the first post where the results are plotted along a circle

for (int ii = 0 ; ii < nn; ii++)

{xx[ii] = ii/nn;

macro grad(u) [dx(u),dy(u)]//

flux[ii] = int1d(Th,circle)( grad(T)â€™*[N.x,N.y])(ii/nn,ii/nn) ;

ofstream gnu ( â€śFlux=.txtâ€ť );

for (int ii = 0; ii < nn; ii++)

gnu << flux[ii] << endl;}

Your question is not clear, so your idea never works.

I think, (but I am not sur) you to compute the flux on arc of size theta

where arc go from 0 to theta.

so the flux on the arc will be something like

int1d on arc(theta) of ( grad(T)â€™*[N.x,N.y])

if ( xc,yc) is the center of the circle the angle of circle is atan2(y-yc,x-xc) in [-pi,p2[ with 0 angle on axis x in -pi, pi

So you can compute make a loop for theta in [0, 2*pi]

int1d(Th,circle)( (pi+atan2(y-yc,x-xc) < theta )*grad(T)â€™*[N.x,N.y]);