[Nek5000-users] [*] Re: Function to calculate volume intergral

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

Mani,


A variable temp can be used only in the subroutines which are called with the arguments of collocation point coordinates and global element number:

ix,iy,iz,ieg


These routines are userf, userq, userbc and useric where temp has the following meaning (updated through NEKUSE)

temp = t(ix,iy,iz,gllel(ieg))

Note that gllel(ieg) returns the local element number on a processor for a global element number ieg.


Similar meaning there have ux,uy,uz and ps(1),ps(2)...,ps(ldimt), namely, velocity components and passive scalars at the point ix,iy,iz,ieg


Best,
Aleks







----- Original Message -----
From: nek5000-users at lists.mcs.anl.gov
To: nek5000-users at lists.mcs.anl.gov
Sent: Thursday, September 2, 2010 7:51:55 AM GMT -06:00 US/Canada Central
Subject: Re: [Nek5000-users] [*] Re: Function to calculate volume intergral

  Thanks! I tried that earlier but it was seg faulting cause I put in 
temperature as 'temp' and not 't'.

Mani
On 09/02/2010 05:56 PM, nek5000-users at lists.mcs.anl.gov wrote:
> Hi Mani,
>
> You can use
>
>         n      = nx1*ny1*nz*nelv
>         vytvol = glsc3(vy,t,bm1,n)
>
> for the volume integral of Vy and temperature
>
> Best,
> Aleks
>
>
> ----- Original Message -----
> From: nek5000-users at lists.mcs.anl.gov
> To: nek5000-users at lists.mcs.anl.gov
> Sent: Thursday, September 2, 2010 12:38:10 AM GMT -06:00 US/Canada Central
> Subject: [Nek5000-users] Function to calculate volume intergral
>
>    Hi,
>
>       Is there a function one can use in Nek that computes volume
> integrals. I'd like to compute the Nusselt number for Rayleigh-Benard
> convection (=1 +<Vy*Temp>), where Vy is the vertical velocity and<>  is
> the volume integral.
>
> Thanks,
> Mani
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>

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