[Nek5000-users] total volume for velocity mesh

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

volvm1 is the volume of the velocity mesh, computed

using the velocity-based quadrature weights (m1-->velocity

quadrature points; m2-->pressure quadrature points)


Thus:


                 n1 = nx1*ny1*nz1*nelv

                 n2 = nx2*ny2*nz2*nelv    ! pressure points


                 volvm1 = glsum(bm1,n1)  ! sum over all points, all processors of the mass matrix


                 volvm2 = glsum(bm2,n2)


volvm1 and volvm2 should be very close to each other.


Incidentally:


              nt = nx1*ny1*nz1*nelt

              voltm1 = glsum(bm1,nt)


will give the volume of the fluid + solid mesh in the case

of conjugate heat transfer, where nelt > nelv


Paul


________________________________
From: Nek5000-users <nek5000-users-bounces at lists.mcs.anl.gov> on behalf of nek5000-users at lists.mcs.anl.gov <nek5000-users at lists.mcs.anl.gov>
Sent: Tuesday, June 26, 2018 10:30:00 AM
To: nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] total volume for velocity mesh

Hi nek,

Can anyone explain how  volvm1 ( total volume for velocity mesh) calculates? Can it be used for different geometries?

Regards,

Jian
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180626/93308880/attachment-0001.html>