[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