From nek5000-users at lists.mcs.anl.gov Tue Jan 5 08:32:17 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 5 Jan 2016 15:32:17 +0100 Subject: [Nek5000-users] Stress formulation In-Reply-To: References: Message-ID: Hi Paul, No, I have run my own example with my grid, but now I'm running "axi" case in examples and after 20000 steps I don't get the parabolic profile. In addition, the solution does not change after 1000 steps. Best, Andrew 2015-12-18 17:19 GMT+01:00 : > > > Hi Andrew > > Are you running the "axi" case in the examples? > > This case started with u=0 as an initial condition. > > A quick glance at the output it produces in f00004, after 200 steps, shows > plots > similar to yours. > > If you run longer (change nsteps in the .rea file), the viscous boundary > layer > will diffuse inward from the walls and you should see a parabolic > profile. This > is what I observed just now. > > Please let me know if this resolves the issue --- if not, I'm happy to > check > further. > > Best, Paul > > ------------------------------ > *From:* nek5000-users-bounces at lists.mcs.anl.gov [ > 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:* Thursday, December 17, 2015 7:07 PM > > *To:* nek5000-users at lists.mcs.anl.gov > *Subject:* Re: [Nek5000-users] Stress formulation > > Hi Paul, > > I think I have misspoken. The thing is that I haven't got the profile corresponding > to the inlet one anywhere in pipe. > I add some clarifications. You can find attached file with velocity > curves for my simulaltions. > I have done the simulations of Poiseuille flow in the cylindrical pipe. > With IFAXIS = true, IFSTRS = false, param (30) = 0 the input profile > doesn't chage. With IFAXIS = true, IFSTSR = true, param (30) = 0 the input > velocity profile changes considerably along the pipe and on the axis there > is strange pressure peak, the nature of which I can not understand. Can I > have clarification on this question? > > Thanks, > Andrew > > > 2015-12-14 0:55 GMT+01:00 : > >> >> Hi Andrew, >> >> In my view, you are already getting the correct flow physics. >> >> What happens at outflow, especially at moderate Reynolds number, is not >> necessarily parallel flow. >> >> The stress formulation is giving the stress-free result. >> >> Paul >> >> >> >> ------------------------------ >> *From:* nek5000-users-bounces at lists.mcs.anl.gov [ >> 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:* Sunday, December 13, 2015 5:19 PM >> *To:* nek5000-users at lists.mcs.anl.gov >> *Subject:* Re: [Nek5000-users] Stress formulation >> >> Hi Paul, >> >> Thank you for your reply. >> Yes, if i increase Reynolds number, the outlet profile become much close >> to the inlet one. But in my simulations the Reynolds number is between 300 >> and 1200. >> Can you recommend me what kind of boundary conditions i can utilise for >> ensure for ensure the flow physics ? I would implement these one in >> code. >> >> Thanks, >> Andrew >> >> >> 2015-12-05 3:37 GMT+01:00 : >> >>> >>> Hi Andrew, >>> >>> I checked into this... I think what is happening is that the 'O ' bc >>> for the stress formulation >>> means stress-free, which is not guaranteed to yield a parabolic profile >>> at the outlet. >>> >>> If you increase your Reynolds number I'm guessing that you'll recover >>> the parabolic >>> profile because the viscous stresses will diminish -- this is what I >>> observed. >>> >>> From my perspective, the outlet boundary is not a region where I would >>> count on accurate >>> physics --- it is, after all, a truncated domain, so I don't generally >>> worry too much about >>> the behavior there. >>> >>> Paul >>> >>> ------------------------------ >>> *From:* nek5000-users-bounces at lists.mcs.anl.gov [ >>> 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, December 01, 2015 2:44 AM >>> *To:* nek5000-users at lists.mcs.anl.gov >>> *Subject:* [Nek5000-users] Stress formulation >>> >>> Hi Neks, >>> >>> I have one question about stress formulation in NEK5000. >>> I tried to make a simulation for Poiseuil flow with stress >>> formulation. In inlet boundary condition use was made of parabolic profile. >>> If IFSTRS=false the outlet profile is exactly the same as inlet one but in >>> case of IFSTRS=true the outlet profile changes significantly (both the >>> maximum value and profile shape). Could you tell me how I can get the >>> parabolic profile at outlet, please? >>> The same thing occurs with pressure values. The maximum pressure >>> value with IFSTRS=true is three times higher then the one with IFSTRS=false. >>> >>> I found that the difference comes from subroutines where the >>> stiffness matrix is calculated (axehlm for IFSTRS = false and axhmsf for >>> IFSTRS = true) but I don't undestand what is happening exactly. >>> >>> Thanks in advance for your help, >>> Best regards, >>> Andrew >>> >>> _______________________________________________ >>> 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 > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Jan 8 13:20:23 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 8 Jan 2016 11:20:23 -0800 Subject: [Nek5000-users] AMG Dump Message-ID: Hi All, I'm trying to set up the AMG solver in Nek5000 to solve for the pressures. I tried following the instructions on Ch5 of the NekDoc available in github, which contains instructions for running the AMG_Dump. However, for some reason I'm unable to generate the *.dat files in my case folder. I have uncommented the relevant lines in the makenek document and tried it. However, it doesn't seem to be working. Could anyone help out? Sincerely, *Amitvikram Dutta* *MASc Student* *Graduate Research Assistant | Okanagan CFD Laboratory* *Solidification** Processing and Simulation Lab* *School of Engineering* *University of British Columbia | Okanagan Campus* -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Jan 8 16:59:12 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 8 Jan 2016 22:59:12 +0000 Subject: [Nek5000-users] AMG Dump In-Reply-To: References: Message-ID: Hi Amitvikram, After you edited makenek, did you "makenek clean" followed by your standard makenek invocation? That's important.... HTH, Paul ________________________________ From: nek5000-users-bounces at lists.mcs.anl.gov [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: Friday, January 08, 2016 1:20 PM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] AMG Dump Hi All, I'm trying to set up the AMG solver in Nek5000 to solve for the pressures. I tried following the instructions on Ch5 of the NekDoc available in github, which contains instructions for running the AMG_Dump. However, for some reason I'm unable to generate the *.dat files in my case folder. I have uncommented the relevant lines in the makenek document and tried it. However, it doesn't seem to be working. Could anyone help out? Sincerely, Amitvikram Dutta MASc Student Graduate Research Assistant | Okanagan CFD Laboratory Solidification Processing and Simulation Lab School of Engineering University of British Columbia | Okanagan Campus -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Sat Jan 9 13:27:48 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Sat, 9 Jan 2016 11:27:48 -0800 Subject: [Nek5000-users] Nek5000-users Digest, Vol 83, Issue 2 In-Reply-To: References: Message-ID: ?Hi Paul, Yes I did do so. The compilation of the case goes fine. After that however, on running nek, instead of? generating .dat file, it simply crashes. Note, the case itself runs just fine without the AMG solver, so I know that's not a problem. *Amitvikram Dutta* *MASc Student* *Graduate Research Assistant | Okanagan CFD Laboratory* *Solidification** Processing and Simulation Lab* *School of Engineering* *University of British Columbia | Okanagan Campus* On Sat, Jan 9, 2016 at 10:00 AM, wrote: > Send Nek5000-users mailing list submissions to > nek5000-users at lists.mcs.anl.gov > > To subscribe or unsubscribe via the World Wide Web, visit > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > or, via email, send a message with subject or body 'help' to > nek5000-users-request at lists.mcs.anl.gov > > You can reach the person managing the list at > nek5000-users-owner at lists.mcs.anl.gov > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of Nek5000-users digest..." > > > Today's Topics: > > 1. AMG Dump (nek5000-users at lists.mcs.anl.gov) > 2. Re: AMG Dump (nek5000-users at lists.mcs.anl.gov) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Fri, 8 Jan 2016 11:20:23 -0800 > From: nek5000-users at lists.mcs.anl.gov > To: nek5000-users at lists.mcs.anl.gov > Subject: [Nek5000-users] AMG Dump > Message-ID: > > Content-Type: text/plain; charset="utf-8" > > Hi All, > > I'm trying to set up the AMG solver in Nek5000 to solve for the pressures. > I tried following the instructions on Ch5 of the NekDoc available in > github, which contains instructions for running the AMG_Dump. > > However, for some reason I'm unable to generate the *.dat files in my case > folder. I have uncommented the relevant lines in the makenek document and > tried it. However, it doesn't seem to be working. > > Could anyone help out? > > Sincerely, > > *Amitvikram Dutta* > > *MASc Student* > > *Graduate Research Assistant | Okanagan CFD Laboratory* > > *Solidification** Processing and Simulation Lab* > > *School of Engineering* > > *University of British Columbia | Okanagan Campus* > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: < > http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20160108/11e28d57/attachment-0001.html > > > > ------------------------------ > > Message: 2 > Date: Fri, 8 Jan 2016 22:59:12 +0000 > From: nek5000-users at lists.mcs.anl.gov > To: "nek5000-users at lists.mcs.anl.gov" > > Subject: Re: [Nek5000-users] AMG Dump > Message-ID: > > Content-Type: text/plain; charset="iso-8859-1" > > > Hi Amitvikram, > > After you edited makenek, did you "makenek clean" followed by your > standard makenek invocation? > > That's important.... > > HTH, > > Paul > > ________________________________ > From: nek5000-users-bounces at lists.mcs.anl.gov [ > 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: Friday, January 08, 2016 1:20 PM > To: nek5000-users at lists.mcs.anl.gov > Subject: [Nek5000-users] AMG Dump > > Hi All, > > I'm trying to set up the AMG solver in Nek5000 to solve for the pressures. > I tried following the instructions on Ch5 of the NekDoc available in > github, which contains instructions for running the AMG_Dump. > > However, for some reason I'm unable to generate the *.dat files in my case > folder. I have uncommented the relevant lines in the makenek document and > tried it. However, it doesn't seem to be working. > > Could anyone help out? > > Sincerely, > > Amitvikram Dutta > > MASc Student > > Graduate Research Assistant | Okanagan CFD Laboratory > > Solidification Processing and Simulation Lab > > School of Engineering > > University of British Columbia | Okanagan Campus > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: < > http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20160108/e92f70e4/attachment-0001.html > > > > ------------------------------ > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > > End of Nek5000-users Digest, Vol 83, Issue 2 > ******************************************** > -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Sat Jan 9 19:47:14 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Sun, 10 Jan 2016 01:47:14 +0000 Subject: [Nek5000-users] Nek5000-users Digest, Vol 83, Issue 2 In-Reply-To: References: , Message-ID: How large is your problem (number of elements and number of processors) ? Paul ________________________________ From: nek5000-users-bounces at lists.mcs.anl.gov [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: Saturday, January 09, 2016 1:27 PM To: nek5000-users at lists.mcs.anl.gov Subject: Re: [Nek5000-users] Nek5000-users Digest, Vol 83, Issue 2 ?Hi Paul, Yes I did do so. The compilation of the case goes fine. After that however, on running nek, instead of? generating .dat file, it simply crashes. Note, the case itself runs just fine without the AMG solver, so I know that's not a problem. Amitvikram Dutta MASc Student Graduate Research Assistant | Okanagan CFD Laboratory Solidification Processing and Simulation Lab School of Engineering University of British Columbia | Okanagan Campus On Sat, Jan 9, 2016 at 10:00 AM, > wrote: Send Nek5000-users mailing list submissions to nek5000-users at lists.mcs.anl.gov To subscribe or unsubscribe via the World Wide Web, visit https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users or, via email, send a message with subject or body 'help' to nek5000-users-request at lists.mcs.anl.gov You can reach the person managing the list at nek5000-users-owner at lists.mcs.anl.gov When replying, please edit your Subject line so it is more specific than "Re: Contents of Nek5000-users digest..." Today's Topics: 1. AMG Dump (nek5000-users at lists.mcs.anl.gov) 2. Re: AMG Dump (nek5000-users at lists.mcs.anl.gov) ---------------------------------------------------------------------- Message: 1 Date: Fri, 8 Jan 2016 11:20:23 -0800 From: nek5000-users at lists.mcs.anl.gov To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] AMG Dump Message-ID: > Content-Type: text/plain; charset="utf-8" Hi All, I'm trying to set up the AMG solver in Nek5000 to solve for the pressures. I tried following the instructions on Ch5 of the NekDoc available in github, which contains instructions for running the AMG_Dump. However, for some reason I'm unable to generate the *.dat files in my case folder. I have uncommented the relevant lines in the makenek document and tried it. However, it doesn't seem to be working. Could anyone help out? Sincerely, *Amitvikram Dutta* *MASc Student* *Graduate Research Assistant | Okanagan CFD Laboratory* *Solidification** Processing and Simulation Lab* *School of Engineering* *University of British Columbia | Okanagan Campus* -------------- next part -------------- An HTML attachment was scrubbed... URL: ------------------------------ Message: 2 Date: Fri, 8 Jan 2016 22:59:12 +0000 From: nek5000-users at lists.mcs.anl.gov To: "nek5000-users at lists.mcs.anl.gov" > Subject: Re: [Nek5000-users] AMG Dump Message-ID: > Content-Type: text/plain; charset="iso-8859-1" Hi Amitvikram, After you edited makenek, did you "makenek clean" followed by your standard makenek invocation? That's important.... HTH, Paul ________________________________ From: nek5000-users-bounces at lists.mcs.anl.gov [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: Friday, January 08, 2016 1:20 PM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] AMG Dump Hi All, I'm trying to set up the AMG solver in Nek5000 to solve for the pressures. I tried following the instructions on Ch5 of the NekDoc available in github, which contains instructions for running the AMG_Dump. However, for some reason I'm unable to generate the *.dat files in my case folder. I have uncommented the relevant lines in the makenek document and tried it. However, it doesn't seem to be working. Could anyone help out? Sincerely, Amitvikram Dutta MASc Student Graduate Research Assistant | Okanagan CFD Laboratory Solidification Processing and Simulation Lab School of Engineering University of British Columbia | Okanagan Campus -------------- next part -------------- An HTML attachment was scrubbed... URL: ------------------------------ _______________________________________________ Nek5000-users mailing list Nek5000-users at lists.mcs.anl.gov https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users End of Nek5000-users Digest, Vol 83, Issue 2 ******************************************** -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Jan 11 08:21:24 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 11 Jan 2016 15:21:24 +0100 Subject: [Nek5000-users] VisIt Message-ID: Hi to everyone, this is no an help Nek5000 request. I need the ubuntu 14.04 VisIt installer, can someone send me it? I tried to download it from their web site, but all links seems to be offline. Thanks a lot. From nek5000-users at lists.mcs.anl.gov Tue Jan 12 07:43:58 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 12 Jan 2016 13:43:58 +0000 Subject: [Nek5000-users] WARNING: DIV(V)-QTL too large Message-ID: dear developers , I'm doing a rotator-stator type of simulation using Pn-Pn method (lx1=8 lxd=12 ). Annular mesh has been created using genbox (box file is attached). Inner radius is considered as the hub and outer radius is the case (or the sidewalls of the disk ). I use following settings in rea file 3 p27 TORDER 20.00000 p93 Number of prev pressure solns saved 5.000000 p94 start projecting vel. after p94 step 5.000000 p95 start projecting pr after p95 step 4.000000 p99 dealiasing:if <0 disable 1.000000 p101 No. of additional filter modes 5.000000E-02 p103 weight of stabilizing filter (.01) Hub and rotor (z=0) rotating in a constant angular velocity. Stator (z=0.2) and side walls are stationary. singularities among rotor and stationary components are resolved using a method which is very similar to vortex example (using an exponential ramp instead of a polynomial ramp) . After 3000-4000 iterations , CFL number reached around 0.48 and I got following warning L1/L2 DIV(V) : -1.6903E-16 1.0198E-01 L1/L2 QTL : 0.0000E+00 0.0000E+00 L1/L2 DIV(V)-QTL: -1.6903E-16 1.0198E-01 WARNING: DIV(V)-QTL too large! I'm not using any inlets or outlets in my simulation. I'm not sure the reason for this warning. can anybody give me a little advice regarding this situation. I can provide more details on request. Desanga -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: 3da.box Type: application/octet-stream Size: 1363 bytes Desc: not available URL: From nek5000-users at lists.mcs.anl.gov Wed Jan 13 14:23:14 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 13 Jan 2016 15:23:14 -0500 Subject: [Nek5000-users] LES equations In-Reply-To: References: Message-ID: Hi all, Is the non-linear term in the LES model decomposed in spectral space? Thanks. David -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Jan 18 09:02:36 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 18 Jan 2016 16:02:36 +0100 Subject: [Nek5000-users] Volume average (local) Message-ID: Hi Neks, I am interested in the scalar average inside a specified volume (similar to the glsc2 routine, but only in a small volume). The volume is given by means of the x/y/z-coordinate. V=(x_2-x_1)*(y_2-y_1)*(z_2-z_1), where ?1? indicates the lower and ?2? the upper coordinate. My first idea was (e.g. for the velocity in streamwise direction): ========================================================= real x1, x2, y1, y2, z1, z2 real lval, lvol, gval, gvol real vxavg integer i,j,k,e [...] ! Specify the volume x1=1.0 x2=2.0 y1=0.0 y2=0.5 z1=0.0 z2=1.0 ! Init lvol=0.0 lval=0.0 gval=0.0 gvol=0.0 ! Loop do e=1,nelt do k=1,nz1 do j=1,ny1 do i=1,nx1 if (xm1(i,j,k,e)>=x1 .and. xm1(i,j,k,e)<=x2 & .and. & ym1(i,j,k,e)>=y1 .and. ym1(i,j,k,e)<=y2 & .and. & zm1(i,j,k,e)>=z1 .and. zm1(i,j,k,e)<=z2 & ) then lvol=lvol+bm1(i,j,k,e) ! Volume lval=lval+vx(i,j,k,e)*bm1(i,j,k,e) ! Weighted scalar value end if end do end do end do end do ! MPI call gop(gval,lval,'+ ',1) call gop(gvol,lvol,'+ ',1) ! Volume averaged streamwise velocity vxavg=gval/gvol ========================================================= This unfortunately doesn't work. Do you see where the problem is? Thanks in advance! Jan F. -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Jan 18 19:07:26 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 18 Jan 2016 20:07:26 -0500 Subject: [Nek5000-users] lid driven cavity Message-ID: Hi Neks, I'm trying to solve a simple 2D cavity problem with Nek and I've got the following error: EXIT: E INVALID BC FOUND in genfast 64 * For your info here is the genbox script I used for the mesh: test5.rea 2 spatial dimension 1 number of fields # #======================================================================= # box_1 1st character .ne. "c" or "C" -8 -8 1 nelx,nely,nelz 0 1 1 x_0 x_Nelx ratio 0 1 1 W ,W ,W ,w V bc's ! west,east,south,north,bottom,top * Also, this is the userbc part in user file: c----------------------------------------------------------------------- subroutine userbc (ix,iy,iz,iside,ieg) include 'SIZE' include 'TSTEP' include 'INPUT' include 'NEKUSE' ux=1. uy=0. return end This should be a straightforward problem but am not sure which part I'm getting wrong. Any advice is appreciated. Hopefully, I can make a sample tutorial for this well-known problem and share with new comers. Thanks, Saleh -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Jan 18 19:23:01 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 18 Jan 2016 20:23:01 -0500 Subject: [Nek5000-users] lid driven cavity In-Reply-To: References: Message-ID: Sorry to spam. I could solve the issue. W ,W ,W ,w V bc's ! west,east,south,north,bottom,top should be changed to W ,W ,W ,v V bc's ! west,east,south,north,bottom,top So just a typo in boundary conditions. On 1/18/2016 8:07 PM, nek5000-users at lists.mcs.anl.gov wrote: > Hi Neks, > > I'm trying to solve a simple 2D cavity problem with Nek and I've got > the following error: > > EXIT: E INVALID BC FOUND in genfast 64 > > * For your info here is the genbox script I used for the mesh: > > > test5.rea > 2 spatial dimension > 1 number of fields > # > #======================================================================= > # > box_1 1st character .ne. "c" or "C" > -8 -8 1 nelx,nely,nelz > 0 1 1 x_0 x_Nelx ratio > 0 1 1 > W ,W ,W ,w V bc's ! west,east,south,north,bottom,top > > * Also, this is the userbc part in user file: > > c----------------------------------------------------------------------- > subroutine userbc (ix,iy,iz,iside,ieg) > include 'SIZE' > include 'TSTEP' > include 'INPUT' > include 'NEKUSE' > > ux=1. > uy=0. > > return > end > > This should be a straightforward problem but am not sure which part > I'm getting wrong. Any advice is appreciated. Hopefully, I can make a > sample tutorial for this well-known problem and share with new comers. > > Thanks, > Saleh > > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Thu Jan 21 23:33:31 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 22 Jan 2016 11:03:31 +0530 Subject: [Nek5000-users] Need for clarification on the SOURCE TERM Message-ID: Dear Sir, How to find out whether the source term which I wish to incorporate in the Navier Stokes equation, will be added to the momentum or energy equation. Thanks & Regards Prashanth Mohanasundaram -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Jan 22 00:56:27 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 22 Jan 2016 14:56:27 +0800 Subject: [Nek5000-users] error in generalev Message-ID: Hello Nek, I just met a trouble in running the case ?high_ray? In ?benard?. I changed the original 2X1 calculation field into a 1x1 cavity, and when I ran the case, the error message ?error in generalev? appeared on the screen. All elements of the matrix ?Aeig? were ?NaN?. I guess that I made some mistakes in the .usr file, maybe wrong B.C., but I?m not sure. I hope I can get some help from you. Thank you! Kind regards Xu Bolun -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Jan 22 06:55:09 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 22 Jan 2016 12:55:09 +0000 Subject: [Nek5000-users] error in generalev In-Reply-To: References: Message-ID: Hi Xu, Did you use genbox to build this case? You need BCs on all 4 sides ... Paul ________________________________ From: nek5000-users-bounces at lists.mcs.anl.gov [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: Friday, January 22, 2016 12:56 AM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] error in generalev Hello Nek, I just met a trouble in running the case ?high_ray? In ?benard?. I changed the original 2X1 calculation field into a 1x1 cavity, and when I ran the case, the error message ?error in generalev? appeared on the screen. All elements of the matrix ?Aeig? were ?NaN?. I guess that I made some mistakes in the .usr file, maybe wrong B.C., but I?m not sure. I hope I can get some help from you. Thank you! Kind regards Xu Bolun -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Jan 22 15:02:38 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 22 Jan 2016 21:02:38 +0000 Subject: [Nek5000-users] Need for clarification on the SOURCE TERM In-Reply-To: References: Message-ID: In userf you specify the forcing for the flow (i.e. Navier-Stokes equation) and in userq the source for the energy. In case you have additional passive scalars you can distinguish between them via the variable 'ifield'. For more info see the online documentation. oana ________________________________________ From: nek5000-users-bounces at lists.mcs.anl.gov [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: Thursday, January 21, 2016 11:33 PM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] Need for clarification on the SOURCE TERM Dear Sir, How to find out whether the source term which I wish to incorporate in the Navier Stokes equation, will be added to the momentum or energy equation. Thanks & Regards Prashanth Mohanasundaram From nek5000-users at lists.mcs.anl.gov Sat Jan 23 18:43:58 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Sun, 24 Jan 2016 00:43:58 +0000 Subject: [Nek5000-users] WARNING: DIV(V)-QTL too large In-Reply-To: References: Message-ID: dear neks, nobody is replying me??? Thanks > On 12 Jan 2016, at 13:43, Desanga Fernando wrote: > > > dear developers , > > I'm doing a rotator-stator type of simulation using Pn-Pn method (lx1=8 lxd=12 ). Annular mesh has been created using genbox (box file is attached). Inner radius is considered as the hub and outer radius is the case (or the sidewalls of the disk ). I use following settings in rea file > > 3 p27 TORDER > 20.00000 p93 Number of prev pressure solns saved > 5.000000 p94 start projecting vel. after p94 step > 5.000000 p95 start projecting pr after p95 step > 4.000000 p99 dealiasing:if <0 disable > 1.000000 p101 No. of additional filter modes > 5.000000E-02 p103 weight of stabilizing filter (.01) > > Hub and rotor (z=0) rotating in a constant angular velocity. Stator (z=0.2) and side walls are stationary. singularities among rotor and stationary components are resolved using a method which is very similar to vortex example (using an exponential ramp instead of a polynomial ramp) . After 3000-4000 iterations , CFL number reached around 0.48 and I got following warning > > > L1/L2 DIV(V) : -1.6903E-16 1.0198E-01 > L1/L2 QTL : 0.0000E+00 0.0000E+00 > L1/L2 DIV(V)-QTL: -1.6903E-16 1.0198E-01 > WARNING: DIV(V)-QTL too large! > > I'm not using any inlets or outlets in my simulation. I'm not sure the reason for this warning. can anybody give me a little advice regarding this situation. I can provide more details on request. > > > Desanga > > <3da.box> From nek5000-users at lists.mcs.anl.gov Sun Jan 24 04:00:20 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Sun, 24 Jan 2016 11:00:20 +0100 Subject: [Nek5000-users] WARNING: DIV(V)-QTL too large In-Reply-To: References: Message-ID: You?ll get a WARNING because the difference in divergence is quite large. Is your simulation blowing up? I guess in your case it is related to temporal/spatial resolution. Try to use a smaller time step and/or increased resolution. Cheers, Stefan Am 1/24/16, 1:43 AM schrieb "nek5000-users-bounces at lists.mcs.anl.gov on behalf of nek5000-users at lists.mcs.anl.gov" unter : >dear neks, > >nobody is replying me??? > >Thanks >> On 12 Jan 2016, at 13:43, Desanga Fernando wrote: >> >> >> dear developers , >> >> I'm doing a rotator-stator type of simulation using Pn-Pn method (lx1=8 >>lxd=12 ). Annular mesh has been created using genbox (box file is >>attached). Inner radius is considered as the hub and outer radius is the >>case (or the sidewalls of the disk ). I use following settings in rea >>file >> >> 3 p27 TORDER >> 20.00000 p93 Number of prev pressure solns saved >> 5.000000 p94 start projecting vel. after p94 step >> 5.000000 p95 start projecting pr after p95 step >> 4.000000 p99 dealiasing:if <0 disable >> 1.000000 p101 No. of additional filter modes >> 5.000000E-02 p103 weight of stabilizing filter (.01) >> >> Hub and rotor (z=0) rotating in a constant angular velocity. Stator >>(z=0.2) and side walls are stationary. singularities among rotor and >>stationary components are resolved using a method which is very similar >>to vortex example (using an exponential ramp instead of a polynomial >>ramp) . After 3000-4000 iterations , CFL number reached around 0.48 and >>I got following warning >> >> >> L1/L2 DIV(V) : -1.6903E-16 1.0198E-01 >> L1/L2 QTL : 0.0000E+00 0.0000E+00 >> L1/L2 DIV(V)-QTL: -1.6903E-16 1.0198E-01 >> WARNING: DIV(V)-QTL too large! >> >> I'm not using any inlets or outlets in my simulation. I'm not sure the >>reason for this warning. can anybody give me a little advice regarding >>this situation. I can provide more details on request. >> >> >> Desanga >> >> <3da.box> > >_______________________________________________ >Nek5000-users mailing list >Nek5000-users at lists.mcs.anl.gov >https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users From nek5000-users at lists.mcs.anl.gov Mon Jan 25 14:08:09 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 25 Jan 2016 20:08:09 +0000 Subject: [Nek5000-users] Nekbone correctness testing In-Reply-To: References: Message-ID: -----BEGIN PGP SIGNED MESSAGE----- Hash: SHA512 Dear Neks, I was wondering if someone could help us to identify correctness tests for Nekbone. How can we make sure that after modifying the code we still get a valid output? Is there anything we can use from the printed output, for instance |cg: ||0| |3||.9174E+||02| |comm time ||0||.2370E-||01| |cg: ||101| |2||.9280E-||05| |4||.3327E-||01| |7||.6135E-||01| |2||.5748E-||09| |cg: ||0| |4||.0900E+||02| |comm time ||0||.2232E-||01| |cg: ||101| |1||.0800E-||05| |4||.4507E-||01| |7||.4233E-||01| |3||.4828E-||10| | | |nelt = ||64||, np = ||8||, nx1 = ||10||, elements = ||512| |Tot MFlops = ||4||.0983E+||04||, MFlops = ||5||.1228E+||03| |Setup Flop = ||8||.8960E+||08||, Solver Flop = ||9||.6960E+||07| |Solve Time = ||0||.1926E+||00| |Avg MFlops = ||4||.0983E+||04| Thanks for your help! -----BEGIN PGP SIGNATURE----- Comment: GPGTools - https://gpgtools.org iQEcBAEBCgAGBQJWpoCpAAoJEIhcgrtekvd+sZMH/RZgY0BlJyYZfvZ5gassbirI Z9lbKaoxMRlrxrHREutBFeo2ABRub2YjlVrulcJy95M503eMIzLp+/AbxLgDib7H FZmgQnXYZUOaM/oGjhVkU7Li9nMb2N0dvIxefvOkvNr8FL7Mnz/aGvdwATCx/RTf L5fqx13+ZoqdDAUAoCyslkXuQKgn1ntj8y4lH+kEve1mHQzuvskRBqQa4NgcIc+u M1F3OCCGv7Ik/h/U1BJE1fade+XhCclfXqG8KeYLopYh/D/zyDkc3aGIDceczknG iQ742rKX2Xd7GIZueAhSk1TS1LiDbA/0QExXYKNb1ilLoTxc/3g4enx4k2cIVAE= =fwgy -----END PGP SIGNATURE----- -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: not available URL: From nek5000-users at lists.mcs.anl.gov Mon Jan 25 20:02:30 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 26 Jan 2016 10:02:30 +0800 (GMT+08:00) Subject: [Nek5000-users] BC type Message-ID: Hi,Nek I want to set Neumann boundary conditions of temperature on the right and left sides of a cavity. I wonder how to set it in .box file.Thank you for your help. Kind regards Xu Bolun -- ??? -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Jan 25 20:32:56 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 26 Jan 2016 02:32:56 +0000 Subject: [Nek5000-users] BC type In-Reply-To: References: Message-ID: Dear Xu, You either want "I " for insulated (zero flux), or "f " for inhomogeneous flux specified in userbc() as: flux = 5*y say... (or whatever functional form you want). If you want a nonzero flux that is a constant, you can have "F " in the .rea file, with VAL as the first real in the boundary condition line. Some more information can be found here: http://www.mcs.anl.gov/~oanam/old_page/index.php/Boundary_Definition Paul ________________________________ From: nek5000-users-bounces at lists.mcs.anl.gov [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: Monday, January 25, 2016 8:02 PM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] BC type Hi,Nek I want to set Neumann boundary conditions of temperature on the right and left sides of a cavity. I wonder how to set it in .box file.Thank you for your help. Kind regards Xu Bolun -- ??? -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Tue Jan 26 15:21:11 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 26 Jan 2016 16:21:11 -0500 Subject: [Nek5000-users] Local Refinement Message-ID: Hi All, I have simulated an eccentric stenosis problem at Re=1000 to verify against Varghese et al., JFM 2007 but my RMS and TKE plots are relatively damped compared to their study. The mesh used in their study appears to show local refinement near the jet breakdown region, which would indicated non-conforming elements. On the other hand, I have no local refinement where the jet breakdown occurs, which probably explains the discrepancy. Without local refinement, my mesh size would be on the order of ~20M grid points, compared to 5M grid points used by Varghese et al. to achieve DNS. I have searched the documentations but couldn?t find a clear answer on this. I have two questions and I would greatly appreciate if someone could answer them, or point me in the correct direction: 1) Is local refinement already implemented in the current version of NEK5000? If yes, is there any documentation available to which I can refer? 2) If local refinement is implemented, what is the best method to refine my mesh? Can I somehow loop over the elements and tag the ones I want to refine in .usr file? Or is there a graphical interface to do so? Thank you very much in advance and any help is greatly appreciated :) Regards, Owais Khan From nek5000-users at lists.mcs.anl.gov Fri Jan 29 03:12:29 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 29 Jan 2016 17:12:29 +0800 Subject: [Nek5000-users] Point votices at inlet Message-ID: Hi neks, Is there a way to prescribe point vortices at the inlet? I am trying to apply the vortex method (see attached link for description) at the inlet to generate turbulence. In that formulation, we have | x_i - x |. I am not sure how to incorporate this part to the usr file. Has anybody done this before? http://aerojet.engr.ucdavis.edu/fluenthelp/html/ug/node509.htm Many thanks Rif -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Jan 29 06:43:17 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 29 Jan 2016 12:43:17 +0000 Subject: [Nek5000-users] Two-way coupled LPT In-Reply-To: References: Message-ID: Hi Neks, I'm currently working on improving my two-way coupling set up for my Lagrangian Particle Tracker which I'm developing as part of my PhD. To do this I require to know the coordinates of the corners of each cell (the space around a mesh gridpoint), in order to determine whether a particle is contained in that cell. I am having trouble doing this for the entire mesh, and I think it's something to do with the parallel processing. My approach is as follows: --------------------------------------------------- if(TWOWAYON.EQ.1) then do e = 1, nelv ieg = lglel(e) write(*,*) 'IEG:',e write(*,*) 'nelv:',nelv write(*,*) 'nelt:',nelt do i = 1, nx1*nelx do j = 1, ny1*nely do k = 1, nz1*nelz c write(*,*) 'X:',i,'Y:',j,'Z:',k if(i.GE.2) then xmc(i,j,k,e) = xm1(i,j,k,e) - $ (xm1(i,j,k,e)-xm1(i-1,j,k,e))*0.5D0 else xmc(i,j,k,e) = xm1(i,j,k,e) endif if(i.LE.7) then xxc(i,j,k,e) = xm1(i,j,k,e) + $ (xm1(i+1,j,k,e)-xm1(i,j,k,e))*0.5D0 else xxc(i,j,k,e) = xm1(i,j,k,e) endif if(j.GE.2) then ymc(i,j,k,e) = ym1(i,j,k,e) - $ (ym1(i,j,k,e)-ym1(i,j-1,k,e))*0.5D0 else ymc(i,j,k,e) = ym1(i,j,k,e) endif if(j.LE.7) then yxc(i,j,k,e) = ym1(i,j,k,e) + $ (ym1(i,j+1,k,e)-ym1(i,j,k,e))*0.5D0 else yxc(i,j,k,e) = ym1(i,j,k,e) endif if(k.GE.2) then zmc(i,j,k,e) = zm1(i,j,k,e) - $ (zm1(i,j,k,e)-zm1(i,j,k-1,e))*0.5D0 else zmc(i,j,k,e) = zm1(i,j,k,e) endif if(k.LE.7) then zxc(i,j,k,e) = zm1(i,j,k,e) + $ (zm1(i,j,k+1,e)-zm1(i,j,k,e))*0.5D0 else zxc(i,j,k,e) = zm1(i,j,k,e) endif vol(i,j,k,e) = $ (xxc(i,j,k,e) - xmc(i,j,k,e)) $ *(yxc(i,j,k,e) - ymc(i,j,k,e)) $ *(zxc(i,j,k,e) - zmc(i,j,k,e)) c write(*,*) 'Xmin: ', xmc(i,j,k,e) c write(*,*) 'Xmax: ', xxc(i,j,k,e) c write(*,*) 'Ymin: ', ymc(i,j,k,e) c write(*,*) 'Ymax: ', yxc(i,j,k,e) c write(*,*) 'Zmin: ', zmc(i,j,k,e) c write(*,*) 'Zmax: ', zxc(i,j,k,e) c write(*,*) 'VOL : ', vol(i,j,k,e) enddo enddo enddo enddo endif ----------------------------------- Now if I sum up the volume using: tvol = glsum(vol,n), I get a weird volume that only seems to represent one element (ish). If you can see where I'm going wrong I'd greatly appreciate it. Thanks, Lee Mortimer University of Leeds -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Jan 29 06:49:02 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 29 Jan 2016 12:49:02 +0000 Subject: [Nek5000-users] Low-Mach equations with source term for mass Message-ID: Hi Neks, We want to solve the Low-Mach equations with a source term S in the continuity equation, i.e. d(rho)/dt + div(rho*u) = S. Is there a way to define S via the Fortan77 interface? Best, Philipp Weiss -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Jan 29 09:38:39 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 29 Jan 2016 15:38:39 +0000 Subject: [Nek5000-users] Low-Mach equations with source term for mass In-Reply-To: References: Message-ID: Hi Philipp, Thew source term for the Low Mach equations is available in qthermal.f stored as QTL. This is however hard coded based on the ideal gas equation of state. According to how your equations turn out you may want to change this term entirely. Oana ________________________________________ From: nek5000-users-bounces at lists.mcs.anl.gov [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: Friday, January 29, 2016 6:49 AM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] Low-Mach equations with source term for mass Hi Neks, We want to solve the Low-Mach equations with a source term S in the continuity equation, i.e. d(rho)/dt + div(rho*u) = S. Is there a way to define S via the Fortan77 interface? Best, Philipp Weiss From nek5000-users at lists.mcs.anl.gov Sat Jan 30 07:50:19 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Sat, 30 Jan 2016 19:20:19 +0530 Subject: [Nek5000-users] Need for clarification on boundary condition (inlet velocity) Message-ID: Dear Sir I am basically trying to simulate a plunging airfoil using NEK5000. What I thought was, keeping the NACA profile stationary while making the flow to oscillate. I wish to give the inlet velocity boundary condition as oscillating ( some sine or cos function dependent on time). How shall I incorporate this in my subroutine userbc()? Thanks/Regards Prashanth Mohanasundaram -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Sat Jan 30 16:35:26 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Sat, 30 Jan 2016 15:35:26 -0700 Subject: [Nek5000-users] Need for clarification on boundary condition (inlet velocity) In-Reply-To: References: Message-ID: Hi Prasanth, In userbc you can implement something like this. A simple example in userbc() ux = sin(omega*time) uy = cos(omega*time) uz = 0 On Sat, Jan 30, 2016 at 6:50 AM, wrote: > Dear Sir > > I am basically trying to simulate a plunging airfoil using NEK5000. What I > thought was, keeping the NACA profile stationary while making the flow to > oscillate. > I wish to give the inlet velocity boundary condition as oscillating ( some > sine or cos function dependent on time). How shall I incorporate this in my > subroutine userbc()? > > Thanks/Regards > Prashanth Mohanasundaram > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Sat Jan 30 14:48:34 2016 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Sat, 30 Jan 2016 21:48:34 +0100 Subject: [Nek5000-users] Does Nek5000 use any kind of parallelism inside one element? Message-ID: Dear Nek5000 developpers and users, my name is Lars Haupt and I'am writing my phd-thesis at the moment, which has something to do with scalable p-multigrid based approaches. Therefore i have to understand the state of the art multigrid solver Nek5000 and thats why i'am trying to get some insides into Nek5000 parallelism. I read the article about the scalability of the Nek5000 Code which was achieved on the Juelich Supercomputer in 2010. I thought that the impressive strong scaling speedup was a result of an optimised hybrid parallelization strategy (MPI and OpenMP). But the general info on the Nek5000 website says that only MPI does the parallel job. Which means, domain decomposition is the only kind of parallelism? What happens with the huge amount of nodes inside one element? The FAQ tells something about high efficient matrix-matrix multiplications. Which are implemented using MKL-BLAS 3 dgemm routines, i guess? If that is the case, does Nek5000 use any OpenMP based parallelism inside MKL or why not? Many thanks in advance Lars Haupt Dipl.-Math. Lars Haupt Wissenschaftlicher Mitarbeiter Technische Universit?t Dresden Fakult?t Maschinenwesen Institut f?r Energietechnik Professur f?r Geb?udeenergietechnik und W?rmeversorgung 01062 Dresden Tel.: +49(351) 463-37619 Fax.: +49(351) 463-37076 E-Mail: lars.haupt at tu-dresden.de -------------- next part -------------- A non-text attachment was scrubbed... Name: smime.p7s Type: application/pkcs7-signature Size: 4755 bytes Desc: not available URL: