[Nek5000-users] Average solution along lines

[nek5000-users at lists.mcs.anl.gov]

Hi Praveen,

Yes, you are correct in assuming that y_average computes

ua(x) = 1/L int_0^L u(x,y) dy

under the assumptions of: 2D, uniform element length in y,
and tensor-product array of elements (e.g., as generated
by genbox).

I wrote an x_average() routine a few weeks back and have
tested it.  I've just now added this to navier5.f --- hopefully it
will work for you.

Be certain you set lelx, lely, and lelz in SIZE, and likewise
nelx,nely,nelz in (say) usrdat(), where you must also add:

     include 'ZPER'

Of course, nelz=lelz=1 in your case.

I

Paul

----- Original Message -----
From: nek5000-users at lists.mcs.anl.gov
To: "nek5000" <nek5000-users at lists.mcs.anl.gov>
Sent: Tuesday, November 12, 2013 10:17:14 PM
Subject: [Nek5000-users] Average solution along lines



Hello 


For a 2d simulation I want to compute averages on y = constant lines 


ua(y) = (1/L) int_0^L u(x,y) dx 


My mesh is cartesian so this should be easy to do. I have seen the y_average function which seems to do just this but I need a x_average function I guess. I can modify this for x averages but would like to check first that y_average computes averages along lines in 2d case ??? 


Thanks 
praveen 
_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users