From nek5000-users at lists.mcs.anl.gov Fri Jun 8 05:18:23 2012 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 8 Jun 2012 05:18:23 -0500 Subject: [Nek5000-users] Budget of turbulent kinetic energy Message-ID: Hello, I want to calculate the terms in budget of turbulent kinetic energy. Right now I am calculating the avg, rms and rm2 terms using subroutine avg_all(). I was wondering is there a way I can use these to calculate the budget terms. Or do I have to calculate each term separately and find their time averages? I have already spend enough time trying to figure this out. Any help would be greatly appreciated. Regards, Harish. -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Jun 11 13:38:01 2012 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 11 Jun 2012 14:38:01 -0400 Subject: [Nek5000-users] Oscillating Bounday Message-ID: Hi Neks, Is it possible to model a moving boundary which oscillates with time in Nekton? If the answer is yes, I wonder if you could provide me with an example. Regards, Alireza Karimi -- Alireza Karimi PhD Candidate Department of Engineering Science and Mechanics Virginia Polytechnic Institute and State University Blacksburg, VA 24061 From nek5000-users at lists.mcs.anl.gov Mon Jun 11 13:42:02 2012 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 11 Jun 2012 13:42:02 -0500 (CDT) Subject: [Nek5000-users] Oscillating Bounday In-Reply-To: References: Message-ID: Hi Alireza, There is a 3D example in /peris, which is a peristaltic flow example. I also have a 2D one that I set up a few weeks back and will try to get that into the repo. Best regards, Paul On Mon, 11 Jun 2012, nek5000-users at lists.mcs.anl.gov wrote: > Hi Neks, > > Is it possible to model a moving boundary which oscillates with time in > Nekton? If the answer is yes, I wonder if you could provide me with an > example. > > Regards, > Alireza Karimi > > > > -- > Alireza Karimi > PhD Candidate > Department of Engineering Science and Mechanics > Virginia Polytechnic Institute and State University > Blacksburg, VA 24061 > _______________________________________________ > 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 Tue Jun 19 13:18:19 2012 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 19 Jun 2012 13:18:19 -0500 Subject: [Nek5000-users] Setting up solid elements with genbox Message-ID: Hi all, I am attempting to setup a simple 2D conjugate heat transfer problem using genbox. I use the example "turbChannel" as my model. I want to setup a fluid box in between two solid boxes. Therefore, my 2DturbChannel.box file appears: base.rea 2 spatial dimension ( < 0 --> generate .rea/.re2 pair) 2 number of fields Box -8 8 nelx,nely for Box 0 8 1. x0,x1,gain (rescaled in usrdat) -1.0 -0.96 -0.75 -0.45 0.00 0.45 0.75 0.96 1.00 y0,y1,...,yn P ,P ,W ,E bc's (3 chars each!) P ,P ,W ,E bc's (3 chars each!) Box -8 8 nelx,nely for Box 0 8 1. x0,x1,gain (rescaled in usrdat) 1.00 1.04 1.25 1.55 2.00 2.45 2.75 2.96 3.00 y0,y1,...,yn P ,P ,E ,W bc's (3 chars each!) P ,P ,E ,W bc's (3 chars each!) Box -8 8 nelx,nely for Box 0 8 1. x0,x1,gain (rescaled in usrdat) -1.00 -1.04 -1.25 -1.55 -2.00 -2.45 -2.75 -2.96 -3.00 y0,y1,...,yn P ,P ,W ,E bc's (3 chars each!) P ,P ,W ,E bc's (3 chars each!) Running genbox on this file produces a .rea file with 192 fluid elements. I first try to change "nelv" in the .rea file to 64. Running genmap on this file produces a "PERIODIC MISMATCH" error, likely due to the 192 boundary conditions for the fluid field. So I then delete all of the boundary conditions for elements greater than 64 in the fluid (leaving thermal boundary conditions for all 192 elements). Genmap now works, but running nek with this .rea/.map combination gives the error: ERROR READING INITIAL CONDITION/DRIVE FORCE DATA ABORTING IN ROUTINE RDICDF. I thought that this procedure would work, given previous discussion on the problem of conjugate heat transfer in Nek5000 (http://lists.mcs.anl.gov/pipermail/nek5000-users/2010-February/000121.html). Perhaps I generated the .map file from the incorrect form of the .rea file? My question: is there a simple way to deal with solid elements in this configuration without using prenek? Thanks, Kyle From nek5000-users at lists.mcs.anl.gov Tue Jun 19 15:17:21 2012 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 19 Jun 2012 15:17:21 -0500 Subject: [Nek5000-users] Setting up solid elements with genbox In-Reply-To: References: Message-ID: Hi Kyle, Did you take a look at the conj_ht case in the examples folder? It has something pretty close to what you want. Justin On Tue, Jun 19, 2012 at 1:18 PM, wrote: > Hi all, > > I am attempting to setup a simple 2D conjugate heat transfer problem using > genbox. I use the example "turbChannel" as my model. I want to setup a > fluid box in between two solid boxes. Therefore, my 2DturbChannel.box file > appears: > > base.rea > 2 spatial dimension ( < 0 --> generate .rea/.re2 > pair) > 2 number of fields > Box > -8 8 nelx,nely for Box > 0 8 1. x0,x1,gain (rescaled > in usrdat) > -1.0 -0.96 -0.75 -0.45 0.00 0.45 0.75 0.96 1.00 y0,y1,...,yn > P ,P ,W ,E bc's (3 chars each!) > P ,P ,W ,E bc's (3 chars each!) > Box > -8 8 nelx,nely for Box > 0 8 1. x0,x1,gain (rescaled > in usrdat) > 1.00 1.04 1.25 1.55 2.00 2.45 2.75 2.96 3.00 y0,y1,...,yn > P ,P ,E ,W bc's (3 chars each!) > P ,P ,E ,W bc's (3 chars each!) > Box > -8 8 nelx,nely for Box > 0 8 1. x0,x1,gain (rescaled > in usrdat) > -1.00 -1.04 -1.25 -1.55 -2.00 -2.45 -2.75 -2.96 -3.00 y0,y1,...,yn > P ,P ,W ,E bc's (3 chars each!) > P ,P ,W ,E bc's (3 chars each!) > > > Running genbox on this file produces a .rea file with 192 fluid elements. > I first try to change "nelv" in the .rea file to 64. Running genmap on this > file produces a "PERIODIC MISMATCH" error, likely due to the 192 boundary > conditions for the fluid field. So I then delete all of the boundary > conditions for elements greater than 64 in the fluid (leaving thermal > boundary conditions for all 192 elements). Genmap now works, but running > nek with this .rea/.map combination gives the error: > > ERROR READING INITIAL CONDITION/DRIVE FORCE DATA > ABORTING IN ROUTINE RDICDF. > > > I thought that this procedure would work, given previous discussion on the > problem of conjugate heat transfer in Nek5000 ( > http://lists.mcs.anl.gov/pipermail/nek5000-users/2010-February/000121.html). > Perhaps I generated the .map file from the incorrect form of the .rea file? > My question: is there a simple way to deal with solid elements in this > configuration without using prenek? > > Thanks, > > Kyle > _______________________________________________ > 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 Tue Jun 19 18:16:03 2012 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 19 Jun 2012 18:16:03 -0500 (CDT) Subject: [Nek5000-users] Setting up solid elements with genbox In-Reply-To: Message-ID: Thanks Justin. I did look at that example, and it is helpful (especially for setting up the boundary conditions in the .box file, which I have now adjusted to be blank for the fluid bc's in the solid domain). My issue is that I don't understand the general steps necessary to setup an arbitrary domain. Simply changing the dimensions in the .box file, and using the conj_ht.rea as the base .rea file like this: conj_ht.rea 2 spatial dimension ( < 0 --> generate .rea/.re2 pair) 2 number of fields Box fluid -8 8 nelx,nely for Box 0 8 1. x0,x1,gain (rescaled in usrdat) -1.0 -0.96 -0.75 -0.45 0.00 0.45 0.75 0.96 1.00 y0,y1,...,yn P ,P ,W ,W bc's (3 chars each!) P ,P ,E ,E bc's (3 chars each!) Box upper fuel -8 8 nelx,nely for Box 0 8 1. x0,x1,gain (rescaled in usrdat) 1.00 1.04 1.25 1.55 2.00 2.45 2.75 2.96 3.00 y0,y1,...,yn , , , bc's (3 chars each!) I ,I ,E ,I bc's (3 chars each!) Box lower fuel -8 8 nelx,nely for Box 0 8 1. x0,x1,gain (rescaled in usrdat) -1.00 -1.04 -1.25 -1.55 -2.00 -2.45 -2.75 -2.96 -3.00 y0,y1,...,yn , , , bc's (3 chars each!) I ,I ,I ,E bc's (3 chars each!) and then using genbox, changing nelv, deleting the fluid bc's for elements greater than 64, using genmap results in the peculiar error when running nek: setup mesh topology WARNINGa: Detected non-right-handed element. Number 129 C1-4: -0.4000E-01 -0.4000E-01 -0.4000E-01 -0.4000E-01 This is with the latest repo (843). Maybe I am missing some necessary changes to the .rea file? Thanks, Kyle ----- Original Message ----- From: nek5000-users at lists.mcs.anl.gov To: nek5000-users at lists.mcs.anl.gov Sent: Tuesday, June 19, 2012 3:17:21 PM Subject: Re: [Nek5000-users] Setting up solid elements with genbox Hi Kyle, Did you take a look at the conj_ht case in the examples folder? It has something pretty close to what you want. Justin On Tue, Jun 19, 2012 at 1:18 PM, < nek5000-users at lists.mcs.anl.gov > wrote: Hi all, I am attempting to setup a simple 2D conjugate heat transfer problem using genbox. I use the example "turbChannel" as my model. I want to setup a fluid box in between two solid boxes. Therefore, my 2DturbChannel.box file appears: base.rea 2 spatial dimension ( < 0 --> generate .rea/.re2 pair) 2 number of fields Box -8 8 nelx,nely for Box 0 8 1. x0,x1,gain (rescaled in usrdat) -1.0 -0.96 -0.75 -0.45 0.00 0.45 0.75 0.96 1.00 y0,y1,...,yn P ,P ,W ,E bc's (3 chars each!) P ,P ,W ,E bc's (3 chars each!) Box -8 8 nelx,nely for Box 0 8 1. x0,x1,gain (rescaled in usrdat) 1.00 1.04 1.25 1.55 2.00 2.45 2.75 2.96 3.00 y0,y1,...,yn P ,P ,E ,W bc's (3 chars each!) P ,P ,E ,W bc's (3 chars each!) Box -8 8 nelx,nely for Box 0 8 1. x0,x1,gain (rescaled in usrdat) -1.00 -1.04 -1.25 -1.55 -2.00 -2.45 -2.75 -2.96 -3.00 y0,y1,...,yn P ,P ,W ,E bc's (3 chars each!) P ,P ,W ,E bc's (3 chars each!) Running genbox on this file produces a .rea file with 192 fluid elements. I first try to change "nelv" in the .rea file to 64. Running genmap on this file produces a "PERIODIC MISMATCH" error, likely due to the 192 boundary conditions for the fluid field. So I then delete all of the boundary conditions for elements greater than 64 in the fluid (leaving thermal boundary conditions for all 192 elements). Genmap now works, but running nek with this .rea/.map combination gives the error: ERROR READING INITIAL CONDITION/DRIVE FORCE DATA ABORTING IN ROUTINE RDICDF. I thought that this procedure would work, given previous discussion on the problem of conjugate heat transfer in Nek5000 ( http://lists.mcs.anl.gov/pipermail/nek5000-users/2010-February/000121.html ). Perhaps I generated the .map file from the incorrect form of the .rea file? My question: is there a simple way to deal with solid elements in this configuration without using prenek? Thanks, Kyle _______________________________________________ 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 From nek5000-users at lists.mcs.anl.gov Tue Jun 19 18:54:26 2012 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 19 Jun 2012 18:54:26 -0500 Subject: [Nek5000-users] Setting up solid elements with genbox In-Reply-To: References: Message-ID: Hi Kyle, The problem looks like your .box file. You need to list the element boundaries from least to greatest. So for your lower box, reverse the order of the y-points: -3.00 -2.96 etc Justin On Tue, 19 Jun 2012 18:16:03 -0500 (CDT) nek5000-users at lists.mcs.anl.gov wrote: > Thanks Justin. I did look at that example, and it is helpful > (especially for setting up the boundary conditions in the .box file, > which I have now adjusted to be blank for the fluid bc's in the solid > domain). > > My issue is that I don't understand the general steps necessary to > setup an arbitrary domain. Simply changing the dimensions in the .box > file, and using the conj_ht.rea as the base .rea file like this: > > conj_ht.rea > 2 spatial dimension ( < 0 --> > generate .rea/.re2 pair) 2 number of fields > > Box fluid > -8 8 nelx,nely for Box > 0 8 1. x0,x1,gain > (rescaled in usrdat) -1.0 -0.96 -0.75 -0.45 0.00 0.45 0.75 0.96 > 1.00 y0,y1,...,yn P ,P ,W ,W > bc's (3 chars each!) P ,P ,E ,E > bc's (3 chars each!) Box upper fuel > -8 8 nelx,nely for Box > 0 8 1. x0,x1,gain > (rescaled in usrdat) 1.00 1.04 1.25 1.55 2.00 2.45 2.75 2.96 3.00 > y0,y1,...,yn , , , bc's (3 chars > each!) I ,I ,E ,I bc's (3 chars > each!) Box lower fuel > -8 8 nelx,nely for Box > 0 8 1. x0,x1,gain > (rescaled in usrdat) -1.00 -1.04 -1.25 -1.55 -2.00 -2.45 -2.75 -2.96 > -3.00 y0,y1,...,yn , , , > bc's (3 chars each!) I ,I ,I ,E > bc's (3 chars each!) > > and then using genbox, changing nelv, deleting the fluid bc's for > elements greater than 64, using genmap results in the peculiar error > when running nek: > > setup mesh topology > > WARNINGa: Detected non-right-handed element. > Number 129 C1-4: -0.4000E-01 -0.4000E-01 -0.4000E-01 > -0.4000E-01 > > This is with the latest repo (843). Maybe I am missing some necessary > changes to the .rea file? > > Thanks, > > Kyle > > ----- Original Message ----- > From: nek5000-users at lists.mcs.anl.gov > To: nek5000-users at lists.mcs.anl.gov > Sent: Tuesday, June 19, 2012 3:17:21 PM > Subject: Re: [Nek5000-users] Setting up solid elements with genbox > > > Hi Kyle, > > Did you take a look at the conj_ht case in the examples folder? It > has something pretty close to what you want. > > Justin > > > On Tue, Jun 19, 2012 at 1:18 PM, < nek5000-users at lists.mcs.anl.gov > > wrote: > > > Hi all, > > I am attempting to setup a simple 2D conjugate heat transfer problem > using genbox. I use the example "turbChannel" as my model. I want to > setup a fluid box in between two solid boxes. Therefore, my > 2DturbChannel.box file appears: > > base.rea > 2 spatial dimension ( < 0 --> generate .rea/.re2 pair) > 2 number of fields > Box > -8 8 nelx,nely for Box > 0 8 1. x0,x1,gain (rescaled in usrdat) > -1.0 -0.96 -0.75 -0.45 0.00 0.45 0.75 0.96 1.00 y0,y1,...,yn > P ,P ,W ,E bc's (3 chars each!) > P ,P ,W ,E bc's (3 chars each!) > Box > -8 8 nelx,nely for Box > 0 8 1. x0,x1,gain (rescaled in usrdat) > 1.00 1.04 1.25 1.55 2.00 2.45 2.75 2.96 3.00 y0,y1,...,yn > P ,P ,E ,W bc's (3 chars each!) > P ,P ,E ,W bc's (3 chars each!) > Box > -8 8 nelx,nely for Box > 0 8 1. x0,x1,gain (rescaled in usrdat) > -1.00 -1.04 -1.25 -1.55 -2.00 -2.45 -2.75 -2.96 -3.00 y0,y1,...,yn > P ,P ,W ,E bc's (3 chars each!) > P ,P ,W ,E bc's (3 chars each!) > > > Running genbox on this file produces a .rea file with 192 fluid > elements. I first try to change "nelv" in the .rea file to 64. > Running genmap on this file produces a "PERIODIC MISMATCH" error, > likely due to the 192 boundary conditions for the fluid field. So I > then delete all of the boundary conditions for elements greater than > 64 in the fluid (leaving thermal boundary conditions for all 192 > elements). Genmap now works, but running nek with this .rea/.map > combination gives the error: > > ERROR READING INITIAL CONDITION/DRIVE FORCE DATA > ABORTING IN ROUTINE RDICDF. > > > I thought that this procedure would work, given previous discussion > on the problem of conjugate heat transfer in Nek5000 > ( http://lists.mcs.anl.gov/pipermail/nek5000-users/2010-February/000121.html > ). Perhaps I generated the .map file from the incorrect form of > the .rea file? My question: is there a simple way to deal with solid > elements in this configuration without using prenek? > > Thanks, > > Kyle > _______________________________________________ > 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 From nek5000-users at lists.mcs.anl.gov Tue Jun 19 21:02:17 2012 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 19 Jun 2012 21:02:17 -0500 Subject: [Nek5000-users] Setting up solid elements with genbox In-Reply-To: References: Message-ID: Hello Kyle, I do not know much about the set up for your case, however I see an error in the way the .box file is setup. In a .box file, the points should be in an ascending order. Hence for the following lines in the your .box file the points in y-direction should be in the reverse order, i.e. start with -3.0 to -1.0. Box lower fuel -8 8 nelx,nely for Box 0 8 1. x0,x1,gain (rescaled in usrdat) -1.00 -1.04 -1.25 -1.55 -2.00 -2.45 -2.75 -2.96 -3.00 y0,y1,...,yn , , , bc's (3 chars each!) I ,I ,I ,E bc's (3 chars each!) Hope this helps. Harish. On Tue, Jun 19, 2012 at 6:16 PM, wrote: > Thanks Justin. I did look at that example, and it is helpful (especially > for setting up the boundary conditions in the .box file, which I have now > adjusted to be blank for the fluid bc's in the solid domain). > > My issue is that I don't understand the general steps necessary to setup > an arbitrary domain. Simply changing the dimensions in the .box file, and > using the conj_ht.rea as the base .rea file like this: > > conj_ht.rea > 2 spatial dimension ( < 0 --> generate .rea/.re2 > pair) > 2 number of fields > > Box fluid > -8 8 nelx,nely for Box > 0 8 1. x0,x1,gain (rescaled > in usrdat) > -1.0 -0.96 -0.75 -0.45 0.00 0.45 0.75 0.96 1.00 y0,y1,...,yn > P ,P ,W ,W bc's (3 chars each!) > P ,P ,E ,E bc's (3 chars each!) > Box upper fuel > -8 8 nelx,nely for Box > 0 8 1. x0,x1,gain (rescaled > in usrdat) > 1.00 1.04 1.25 1.55 2.00 2.45 2.75 2.96 3.00 y0,y1,...,yn > , , , bc's (3 chars each!) > I ,I ,E ,I bc's (3 chars each!) > Box lower fuel > -8 8 nelx,nely for Box > 0 8 1. x0,x1,gain (rescaled > in usrdat) > -1.00 -1.04 -1.25 -1.55 -2.00 -2.45 -2.75 -2.96 -3.00 y0,y1,...,yn > , , , bc's (3 chars each!) > I ,I ,I ,E bc's (3 chars each!) > > and then using genbox, changing nelv, deleting the fluid bc's for elements > greater than 64, using genmap results in the peculiar error when running > nek: > > setup mesh topology > > WARNINGa: Detected non-right-handed element. > Number 129 C1-4: -0.4000E-01 -0.4000E-01 -0.4000E-01 -0.4000E-01 > > This is with the latest repo (843). Maybe I am missing some necessary > changes to the .rea file? > > Thanks, > > Kyle > > ----- Original Message ----- > From: nek5000-users at lists.mcs.anl.gov > To: nek5000-users at lists.mcs.anl.gov > Sent: Tuesday, June 19, 2012 3:17:21 PM > Subject: Re: [Nek5000-users] Setting up solid elements with genbox > > > Hi Kyle, > > Did you take a look at the conj_ht case in the examples folder? It has > something pretty close to what you want. > > Justin > > > On Tue, Jun 19, 2012 at 1:18 PM, < nek5000-users at lists.mcs.anl.gov > > wrote: > > > Hi all, > > I am attempting to setup a simple 2D conjugate heat transfer problem using > genbox. I use the example "turbChannel" as my model. I want to setup a > fluid box in between two solid boxes. Therefore, my 2DturbChannel.box file > appears: > > base.rea > 2 spatial dimension ( < 0 --> generate .rea/.re2 pair) > 2 number of fields > Box > -8 8 nelx,nely for Box > 0 8 1. x0,x1,gain (rescaled in usrdat) > -1.0 -0.96 -0.75 -0.45 0.00 0.45 0.75 0.96 1.00 y0,y1,...,yn > P ,P ,W ,E bc's (3 chars each!) > P ,P ,W ,E bc's (3 chars each!) > Box > -8 8 nelx,nely for Box > 0 8 1. x0,x1,gain (rescaled in usrdat) > 1.00 1.04 1.25 1.55 2.00 2.45 2.75 2.96 3.00 y0,y1,...,yn > P ,P ,E ,W bc's (3 chars each!) > P ,P ,E ,W bc's (3 chars each!) > Box > -8 8 nelx,nely for Box > 0 8 1. x0,x1,gain (rescaled in usrdat) > -1.00 -1.04 -1.25 -1.55 -2.00 -2.45 -2.75 -2.96 -3.00 y0,y1,...,yn > P ,P ,W ,E bc's (3 chars each!) > P ,P ,W ,E bc's (3 chars each!) > > > Running genbox on this file produces a .rea file with 192 fluid elements. > I first try to change "nelv" in the .rea file to 64. Running genmap on this > file produces a "PERIODIC MISMATCH" error, likely due to the 192 boundary > conditions for the fluid field. So I then delete all of the boundary > conditions for elements greater than 64 in the fluid (leaving thermal > boundary conditions for all 192 elements). Genmap now works, but running > nek with this .rea/.map combination gives the error: > > ERROR READING INITIAL CONDITION/DRIVE FORCE DATA > ABORTING IN ROUTINE RDICDF. > > > I thought that this procedure would work, given previous discussion on the > problem of conjugate heat transfer in Nek5000 ( > http://lists.mcs.anl.gov/pipermail/nek5000-users/2010-February/000121.html). Perhaps I generated the .map file from the incorrect form of the .rea > file? My question: is there a simple way to deal with solid elements in > this configuration without using prenek? > > Thanks, > > Kyle > _______________________________________________ > 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 Wed Jun 20 02:07:31 2012 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 20 Jun 2012 12:37:31 +0530 Subject: [Nek5000-users] Reg: Rayleigh Benard convection Message-ID: Dear Sir, Iam new to nek5000, so far I have ran certain example problems and a simple case where the box file is as follows: ddfree.rea 2 2 Box 1 -10 -8 0.0 2.0 1. 0.0 0.25 1.25 P ,P ,SYM, , P ,P ,t , , Box 2 -10 -10 0.0 2.0 1. 0.25 0.75 1. P ,P , , , P ,P , , , Box 3 -10 -8 0.0 2.0 1. 0.75 1.0 0.8 P ,P , ,SYM, P ,P , ,t , As a new nekuser I have used 8 processors for this problem,(as I dont know how to determine the correct number of processor) . I obtained .f files but these files didn't open in VISIT even when I changed it to .fld files. It was not able to read those files. I have observed that fld files directly obtained can be read by VISIT. Can you please tell me how do we get .fld and not .f files while running nek. I kindly request you to give a benchmark problem where Rayleigh Benard convection is applied. It will be of great help to me. with warm regards, -- *Vinaya Sivanandan* *Pune* ** -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Wed Jun 20 11:51:01 2012 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 20 Jun 2012 11:51:01 -0500 (CDT) Subject: [Nek5000-users] Setting up solid elements with genbox In-Reply-To: Message-ID: Thanks everyone, Switching the order of the mesh points in ascending order fixed my right-handed element error. Best, Kyle ----- Original Message ----- From: nek5000-users at lists.mcs.anl.gov To: nek5000-users at lists.mcs.anl.gov Sent: Tuesday, June 19, 2012 9:02:17 PM Subject: Re: [Nek5000-users] Setting up solid elements with genbox Hello Kyle, I do not know much about the set up for your case, however I see an error in the way the .box file is setup. In a .box file, the points should be in an ascending order. Hence for the following lines in the your .box file the points in y-direction should be in the reverse order, i.e. start with -3.0 to -1.0. Box lower fuel -8 8 nelx,nely for Box 0 8 1. x0,x1,gain (rescaled in usrdat) -1.00 -1.04 -1.25 -1.55 -2.00 -2.45 -2.75 -2.96 -3.00 y0,y1,...,yn , , , bc's (3 chars each!) I ,I ,I ,E bc's (3 chars each!) Hope this helps. Harish. On Tue, Jun 19, 2012 at 6:16 PM, < nek5000-users at lists.mcs.anl.gov > wrote: Thanks Justin. I did look at that example, and it is helpful (especially for setting up the boundary conditions in the .box file, which I have now adjusted to be blank for the fluid bc's in the solid domain). My issue is that I don't understand the general steps necessary to setup an arbitrary domain. Simply changing the dimensions in the .box file, and using the conj_ht.rea as the base .rea file like this: conj_ht.rea 2 spatial dimension ( < 0 --> generate .rea/.re2 pair) 2 number of fields Box fluid -8 8 nelx,nely for Box 0 8 1. x0,x1,gain (rescaled in usrdat) -1.0 -0.96 -0.75 -0.45 0.00 0.45 0.75 0.96 1.00 y0,y1,...,yn P ,P ,W ,W bc's (3 chars each!) P ,P ,E ,E bc's (3 chars each!) Box upper fuel -8 8 nelx,nely for Box 0 8 1. x0,x1,gain (rescaled in usrdat) 1.00 1.04 1.25 1.55 2.00 2.45 2.75 2.96 3.00 y0,y1,...,yn , , , bc's (3 chars each!) I ,I ,E ,I bc's (3 chars each!) Box lower fuel -8 8 nelx,nely for Box 0 8 1. x0,x1,gain (rescaled in usrdat) -1.00 -1.04 -1.25 -1.55 -2.00 -2.45 -2.75 -2.96 -3.00 y0,y1,...,yn , , , bc's (3 chars each!) I ,I ,I ,E bc's (3 chars each!) and then using genbox, changing nelv, deleting the fluid bc's for elements greater than 64, using genmap results in the peculiar error when running nek: setup mesh topology WARNINGa: Detected non-right-handed element. Number 129 C1-4: -0.4000E-01 -0.4000E-01 -0.4000E-01 -0.4000E-01 This is with the latest repo (843). Maybe I am missing some necessary changes to the .rea file? Thanks, Kyle ----- Original Message ----- From: nek5000-users at lists.mcs.anl.gov To: nek5000-users at lists.mcs.anl.gov Sent: Tuesday, June 19, 2012 3:17:21 PM Subject: Re: [Nek5000-users] Setting up solid elements with genbox Hi Kyle, Did you take a look at the conj_ht case in the examples folder? It has something pretty close to what you want. Justin On Tue, Jun 19, 2012 at 1:18 PM, < nek5000-users at lists.mcs.anl.gov > wrote: Hi all, I am attempting to setup a simple 2D conjugate heat transfer problem using genbox. I use the example "turbChannel" as my model. I want to setup a fluid box in between two solid boxes. Therefore, my 2DturbChannel.box file appears: base.rea 2 spatial dimension ( < 0 --> generate .rea/.re2 pair) 2 number of fields Box -8 8 nelx,nely for Box 0 8 1. x0,x1,gain (rescaled in usrdat) -1.0 -0.96 -0.75 -0.45 0.00 0.45 0.75 0.96 1.00 y0,y1,...,yn P ,P ,W ,E bc's (3 chars each!) P ,P ,W ,E bc's (3 chars each!) Box -8 8 nelx,nely for Box 0 8 1. x0,x1,gain (rescaled in usrdat) 1.00 1.04 1.25 1.55 2.00 2.45 2.75 2.96 3.00 y0,y1,...,yn P ,P ,E ,W bc's (3 chars each!) P ,P ,E ,W bc's (3 chars each!) Box -8 8 nelx,nely for Box 0 8 1. x0,x1,gain (rescaled in usrdat) -1.00 -1.04 -1.25 -1.55 -2.00 -2.45 -2.75 -2.96 -3.00 y0,y1,...,yn P ,P ,W ,E bc's (3 chars each!) P ,P ,W ,E bc's (3 chars each!) Running genbox on this file produces a .rea file with 192 fluid elements. I first try to change "nelv" in the .rea file to 64. Running genmap on this file produces a "PERIODIC MISMATCH" error, likely due to the 192 boundary conditions for the fluid field. So I then delete all of the boundary conditions for elements greater than 64 in the fluid (leaving thermal boundary conditions for all 192 elements). Genmap now works, but running nek with this .rea/.map combination gives the error: ERROR READING INITIAL CONDITION/DRIVE FORCE DATA ABORTING IN ROUTINE RDICDF. I thought that this procedure would work, given previous discussion on the problem of conjugate heat transfer in Nek5000 ( http://lists.mcs.anl.gov/pipermail/nek5000-users/2010-February/000121.html ). Perhaps I generated the .map file from the incorrect form of the .rea file? My question: is there a simple way to deal with solid elements in this configuration without using prenek? Thanks, Kyle _______________________________________________ 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 _______________________________________________ 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 Thu Jun 21 06:51:10 2012 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 21 Jun 2012 06:51:10 -0500 Subject: [Nek5000-users] Version 843 has errors in postpro.f Message-ID: Hello Nek Developers, I updated to version 843 and could not compile Nek due to multiple errors in the postpro.f file. There are illegal characters in the postpro.f file and two hpts() subroutines. Just wanted to bring this to your notice. Regards, Harish. -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Thu Jun 21 14:34:10 2012 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 21 Jun 2012 14:34:10 -0500 Subject: [Nek5000-users] Version 843 has errors in postpro.f In-Reply-To: References: Message-ID: Hi Harish, When you updated to revision 843, svn might have merged your postpro.f files. Try moving postpro.f to a different file name, and updating again. If this doesn't solve your problem, let me know. Thanks, Katie On Thu, Jun 21, 2012 at 6:51 AM, wrote: > Hello Nek Developers, > > I updated to version 843 and could not compile Nek due to multiple errors > in the postpro.f file. > > There are illegal characters in the postpro.f file and two hpts() > subroutines. > > Just wanted to bring this to your notice. > > Regards, > > Harish. > > > > _______________________________________________ > 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 Jun 21 14:47:24 2012 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 21 Jun 2012 14:47:24 -0500 Subject: [Nek5000-users] Version 843 has errors in postpro.f In-Reply-To: References: Message-ID: Hi Katie, Thanks for the reply. I deleted my old Nek folder and downloaded a new one. Now the code compiles w/o errors. Regards, Harish. On Thu, Jun 21, 2012 at 2:34 PM, wrote: > Hi Harish, > > When you updated to revision 843, svn might have merged your postpro.f > files. Try moving postpro.f to a different file name, and updating > again. If this doesn't solve your problem, let me know. > > Thanks, > Katie > > On Thu, Jun 21, 2012 at 6:51 AM, wrote: > >> Hello Nek Developers, >> >> I updated to version 843 and could not compile Nek due to multiple errors >> in the postpro.f file. >> >> There are illegal characters in the postpro.f file and two hpts() >> subroutines. >> >> Just wanted to bring this to your notice. >> >> Regards, >> >> Harish. >> >> >> >> _______________________________________________ >> 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 Jun 22 13:26:28 2012 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 22 Jun 2012 14:26:28 -0400 Subject: [Nek5000-users] postnek on a mac Message-ID: Hello all, I have downloaded the most recent version of the software and am trying to compile postnek on a mac running lion with gfortran and gcc. The command: ./maketools postnek Gives the following error: Warning: Deleted feature: ASSIGN statement at (1) gcc -mcmodel=medium -c -O2 xdriver.c xdriver.c: In function ?TypeChar?: xdriver.c:1790: warning: comparison is always true due to limited range of data type xdriver.c: In function ?error?: xdriver.c:1910: warning: format not a string literal and no format arguments xdriver.c:1910: warning: format not a string literal and no format arguments gfortran -mcmodel=medium -c scrdmp.f gfortran -mcmodel=medium -c coef.f gfortran -mcmodel=medium -c postnek7.f gfortran -mcmodel=medium -c speclib.f gfortran -mcmodel=medium -c mxm.f gcc -mcmodel=medium -c -DUNDERSCORE ../../nek/byte.c gfortran -mcmodel=medium -c ssyev.f gfortran -mcmodel=medium -c iolib.f gfortran -mcmodel=medium -o /Users/mrp/bin/postx postnek.o postnek2.o postnek3.o postnek5.o postnek6.o tsort.o postnek8.o postnek9.o plot.o getfld.o legend.o userf.o revert.o trap.o animate.o genxyz.o screen.o g3d.o subs.o xinterface.o locglob.o postnek5a.o blas.o xdriver.o scrdmp.o coef.o postnek7.o speclib.o mxm.o byte.o ssyev.o iolib.o -L/usr/X11R6/lib -lX11 -lm ld: in xdriver.o, in section __TEXT,__text reloc 3: length != 2 and X86_64_RELOC_GOT not supported for architecture x86_64 collect2: ld returned 1 exit status make[1]: *** [postx] Error 1 make: *** [all] Error 1 Any ideas? Thanks, Mark From nek5000-users at lists.mcs.anl.gov Mon Jun 25 03:39:50 2012 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 25 Jun 2012 14:09:50 +0530 Subject: [Nek5000-users] Reminder Reg: Rayleigh Benard convection Message-ID: Dear Sir, Iam new to nek5000, so far I have ran certain example problems and a simple case where the box file is as follows: ddfree.rea 2 2 Box 1 -10 -8 0.0 2.0 1. 0.0 0.25 1.25 P ,P ,SYM, , P ,P ,t , , Box 2 -10 -10 0.0 2.0 1. 0.25 0.75 1. P ,P , , , P ,P , , , Box 3 -10 -8 0.0 2.0 1. 0.75 1.0 0.8 P ,P , ,SYM, P ,P , ,t , As a new nekuser I have used 8 processors for this problem,(as I dont know how to determine the correct number of processor) . I obtained .f files but these files didn't open in VISIT even when I changed it to .fld files. It was not able to read those files. I have observed that fld files directly obtained can be read by VISIT. Can you please tell me how do we get .fld and not .f files while running nek. I kindly request you to give a benchmark problem where Rayleigh Benard convection is applied. It will be of great help to me. with warm regards, -- *Vinaya Sivanandan* *Pune* ** -- *Vinaya Sivanandan* *Pune* ** -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Jun 25 07:05:56 2012 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 25 Jun 2012 07:05:56 -0500 (CDT) Subject: [Nek5000-users] Reminder Reg: Rayleigh Benard convection In-Reply-To: Message-ID: Hi Vinaya, You need to create a metafile with extension nek5000 http://nek5000.mcs.anl.gov/index.php/VISIT which then you open in VisIt. Also you can generate this file with the script http://www.mcs.anl.gov/~obabko/visnek Best. Aleks ----- Original Message ----- From: nek5000-users at lists.mcs.anl.gov To: nek5000-users at lists.mcs.anl.gov Sent: Monday, June 25, 2012 3:39:50 AM Subject: [Nek5000-users] Reminder Reg: Rayleigh Benard convection Dear Sir, Iam new to nek5000, so far I have ran certain example problems and a simple case where the box file is as follows: ddfree.rea 2 2 Box 1 -10 -8 0.0 2.0 1. 0.0 0.25 1.25 P ,P ,SYM, , P ,P ,t , , Box 2 -10 -10 0.0 2.0 1. 0.25 0.75 1. P ,P , , , P ,P , , , Box 3 -10 -8 0.0 2.0 1. 0.75 1.0 0.8 P ,P , ,SYM, P ,P , ,t , As a new nekuser I have used 8 processors for this problem,(as I dont know how to determine the correct number of processor) . I obtained .f files but these files didn't open in VISIT even when I changed it to .fld files. It was not able to read those files. I have observed that fld files directly obtained can be read by VISIT. Can you please tell me how do we get .fld and not .f files while running nek. I kindly request you to give a benchmark problem where Rayleigh Benard convection is applied. It will be of great help to me. with warm regards, -- Vinaya Sivanandan Pune -- Vinaya Sivanandan Pune _______________________________________________ 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 Tue Jun 26 11:00:38 2012 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 26 Jun 2012 12:00:38 -0400 Subject: [Nek5000-users] Anisotropic conductivity Message-ID: Hi Neks, I want to solve for a convective heat transfer problem with a diagonal conductivity tensor which has different constant values in x and y directions. I assume I should use udiff, but I'm not sure how to implement it. I wonder if you have any idea. Regards, Alireza From nek5000-users at lists.mcs.anl.gov Tue Jun 26 14:10:15 2012 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 26 Jun 2012 14:10:15 -0500 (CDT) Subject: [Nek5000-users] Anisotropic conductivity In-Reply-To: References: Message-ID: Hi Alireza, There is a paper on my website describing the approach: http://www.mcs.anl.gov/~fischer/pubhtml/fischer_sci05.pdf Basically, you need to change the geometric factors after nek has set them. I normally do this in userchk, because the first call to userchk implies that the code is ready to start solving the equations and that (most) all of the requisite variables are set. Therefore, if you change them, they won't be reset against your wishes. (This is trickier if you have ifuservp=.true., because in that case nek will keep resetting the diffusivity..., but maybe it will be ok---we can look into that.) I can send an example where I've set up anisotropic diffusion problem... (attached here, in fact, but I've not tested recently). It's a 2D case. Let me know if you need assistance with the 3D case or have questions with this one -- the only tricky part is figuring out where the components of the symmetric 3x3 tensor are stored. Best regards, Paul On Tue, 26 Jun 2012, nek5000-users at lists.mcs.anl.gov wrote: > Hi Neks, > > I want to solve for a convective heat transfer problem with a diagonal > conductivity tensor which has different constant values in x and y > directions. I assume I should use udiff, but I'm not sure how to implement > it. I wonder if you have any idea. > > Regards, > Alireza > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users -------------- next part -------------- A non-text attachment was scrubbed... Name: aniso.tgz Type: application/x-gtar Size: 4514 bytes Desc: URL: From nek5000-users at lists.mcs.anl.gov Fri Jun 29 16:43:46 2012 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Sat, 30 Jun 2012 03:13:46 +0530 Subject: [Nek5000-users] Corrupt mesh Message-ID: Dear Nek devs, I am using genbox to construct a mesh for a box of aspect ratio 1 x 1 x 1. But for some reason, the mesh seems to be corrupted, from what I see in VisIt (screenshot attached). I have attached the .box file used to construct this mesh. Could you help diagnose the reason for this error? Thanks, Mani -------------- next part -------------- A non-text attachment was scrubbed... Name: high_ray.box Type: application/vnd.previewsystems.box Size: 4065 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: Screenshot.png Type: image/png Size: 186662 bytes Desc: not available URL: From nek5000-users at lists.mcs.anl.gov Fri Jun 29 18:08:57 2012 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 29 Jun 2012 18:08:57 -0500 (CDT) Subject: [Nek5000-users] Corrupt mesh In-Reply-To: Message-ID: Hi Mani, I am curious what is a reason for not using a single Box instead of your current strategy of 27 Boxes with uniform element distribution: 1:-8 -8 -8 2:-12 -8 -8 3:-8 -8 -8 4:-8 -12 -8 5:-12 -12 -8 6:-8 -12 -8 7:-8 -8 -8 8:-12 -8 -8 9:-8 -8 -8 10:-8 -8 -12 11:-12 -8 -12 12:-8 -8 -12 13:-8 -12 -12 14:-12 -12 -12 15:-8 -12 -12 16:-8 -8 -12 17:-12 -8 -12 18:-8 -8 -12 19:-8 -8 -8 20:-12 -8 -8 21:-8 -8 -8 22:-8 -12 -8 23:-12 -12 -8 24:-8 -12 -8 25:-8 -8 -8 26:-12 -8 -8 27:-8 -8 -8 ? Thanks. Aleks ----- Original Message ----- From: nek5000-users at lists.mcs.anl.gov To: nek5000-users at lists.mcs.anl.gov Sent: Friday, June 29, 2012 4:43:46 PM Subject: [Nek5000-users] Corrupt mesh Dear Nek devs, I am using genbox to construct a mesh for a box of aspect ratio 1 x 1 x 1. But for some reason, the mesh seems to be corrupted, from what I see in VisIt (screenshot attached). I have attached the .box file used to construct this mesh. Could you help diagnose the reason for this error? Thanks, Mani _______________________________________________ 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 Fri Jun 29 18:54:52 2012 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 29 Jun 2012 18:54:52 -0500 (CDT) Subject: [Nek5000-users] Corrupt mesh In-Reply-To: Message-ID: Hi Mani, I got it -- you want to have boundary layers although it may be simpler just to give element coordinates explicitly in a single Box... Well, I generated .rea file from your .box file with repo version of genbox that looks fine in postx: http://www.mcs.anl.gov/~obabko/m1.tiff http://www.mcs.anl.gov/~obabko/m2.tiff Best. Aleks ----- Forwarded Message ----- From: nek5000-users at lists.mcs.anl.gov To: nek5000-users at lists.mcs.anl.gov Sent: Friday, June 29, 2012 6:08:57 PM Subject: Re: [Nek5000-users] Corrupt mesh Hi Mani, I am curious what is a reason for not using a single Box instead of your current strategy of 27 Boxes with uniform element distribution: 1:-8 -8 -8 2:-12 -8 -8 3:-8 -8 -8 4:-8 -12 -8 5:-12 -12 -8 6:-8 -12 -8 7:-8 -8 -8 8:-12 -8 -8 9:-8 -8 -8 10:-8 -8 -12 11:-12 -8 -12 12:-8 -8 -12 13:-8 -12 -12 14:-12 -12 -12 15:-8 -12 -12 16:-8 -8 -12 17:-12 -8 -12 18:-8 -8 -12 19:-8 -8 -8 20:-12 -8 -8 21:-8 -8 -8 22:-8 -12 -8 23:-12 -12 -8 24:-8 -12 -8 25:-8 -8 -8 26:-12 -8 -8 27:-8 -8 -8 ? Thanks. Aleks ----- Original Message ----- From: nek5000-users at lists.mcs.anl.gov To: nek5000-users at lists.mcs.anl.gov Sent: Friday, June 29, 2012 4:43:46 PM Subject: [Nek5000-users] Corrupt mesh Dear Nek devs, I am using genbox to construct a mesh for a box of aspect ratio 1 x 1 x 1. But for some reason, the mesh seems to be corrupted, from what I see in VisIt (screenshot attached). I have attached the .box file used to construct this mesh. Could you help diagnose the reason for this error? 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 From nek5000-users at lists.mcs.anl.gov Sat Jun 30 17:43:13 2012 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Sun, 01 Jul 2012 04:13:13 +0530 Subject: [Nek5000-users] Corrupt mesh In-Reply-To: References: Message-ID: Hi Alex, Thank you for checking. Regards, Mani On 06/30/2012 05:24 AM, nek5000-users at lists.mcs.anl.gov wrote: > Hi Mani, > > I got it -- you want to have boundary layers although it may be simpler just to give element coordinates explicitly in a single Box... > > Well, I generated .rea file from your .box file with repo version of genbox that looks fine in postx: > > http://www.mcs.anl.gov/~obabko/m1.tiff > http://www.mcs.anl.gov/~obabko/m2.tiff > > Best. > Aleks > > > > ----- Forwarded Message ----- > From: nek5000-users at lists.mcs.anl.gov > To: nek5000-users at lists.mcs.anl.gov > Sent: Friday, June 29, 2012 6:08:57 PM > Subject: Re: [Nek5000-users] Corrupt mesh > > Hi Mani, > > I am curious what is a reason for not using a single Box instead of your current strategy of 27 Boxes with uniform element distribution: > > 1:-8 -8 -8 > 2:-12 -8 -8 > 3:-8 -8 -8 > 4:-8 -12 -8 > 5:-12 -12 -8 > 6:-8 -12 -8 > 7:-8 -8 -8 > 8:-12 -8 -8 > 9:-8 -8 -8 > 10:-8 -8 -12 > 11:-12 -8 -12 > 12:-8 -8 -12 > 13:-8 -12 -12 > 14:-12 -12 -12 > 15:-8 -12 -12 > 16:-8 -8 -12 > 17:-12 -8 -12 > 18:-8 -8 -12 > 19:-8 -8 -8 > 20:-12 -8 -8 > 21:-8 -8 -8 > 22:-8 -12 -8 > 23:-12 -12 -8 > 24:-8 -12 -8 > 25:-8 -8 -8 > 26:-12 -8 -8 > 27:-8 -8 -8 > > ? > > Thanks. > Aleks > > > > ----- Original Message ----- > From: nek5000-users at lists.mcs.anl.gov > To: nek5000-users at lists.mcs.anl.gov > Sent: Friday, June 29, 2012 4:43:46 PM > Subject: [Nek5000-users] Corrupt mesh > > Dear Nek devs, > > I am using genbox to construct a mesh for a box of aspect ratio 1 x 1 x > 1. But for some reason, the mesh seems to be corrupted, from what I see > in VisIt (screenshot attached). I have attached the .box file used to > construct this mesh. Could you help diagnose the reason for this error? > > 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 >