[Nek5000-users] Average solution along lines
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Wed Nov 13 06:09:14 CST 2013
Yes - that is correct.
----- Original Message -----
From: nek5000-users at lists.mcs.anl.gov
To: "nek5000" <nek5000-users at lists.mcs.anl.gov>
Sent: Wednesday, November 13, 2013 6:06:08 AM
Subject: Re: [Nek5000-users] Average solution along lines
Thank you Paul.
The routine should work for non-uniform cartesian mesh also right, if the dx2 is calculated. This is done in y_average function but commented out for dy2 ???
Regards
praveen
On Wed, Nov 13, 2013 at 5:10 PM, < nek5000-users at lists.mcs.anl.gov > wrote:
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
_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
More information about the Nek5000-users
mailing list