From nek5000-users at lists.mcs.anl.gov Mon Jan 2 07:17:23 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 02 Jan 2017 13:17:23 +0000 Subject: [Nek5000-users] double mesh ? In-Reply-To: References: Message-ID: Hi all, Would anyone like to volunteer a response? Our user kindly requested a response this morning. Many thanks, Ron On Thu, Dec 22, 2016 at 7:40 AM wrote: > Hi > > i am quite new, so i?m sorry if i?m getting things wrong. > > I am trying to simulate a turbulent channel flow, with a passive scalar. > in the streamwise direction, i would like to impose periodic boundary > conditions for the flow, and non-periodic for the scalar. > in fact the scalar has a point source near the inlet, and i want to see > how that develops downstream. > > as far as i understand this is not supported by nek5000, because the last > mesh point is identified with the first mesh point when periodic BC are > imposed. > correct ? > > my question is : could i define two meshes, one for the scalar and one for > the velocity ? > > thanks a lot for your help !!! > agnese > > > Agnese Seminara > -------------------------------- > CNRS > Laboratoire de physique de la mati?re condens?e > Parc Valrose > avenue J Vallot > 06108 Nice, France > +33 (0) 492 076 775 > http://sites.unice.fr/site/aseminara/ > > _______________________________________________ > 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 Jan 3 05:12:14 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 3 Jan 2017 11:12:14 +0000 Subject: [Nek5000-users] convert openfoam mesh to nek5000 Message-ID: Hi all, I want to convert openfoam mesh to nek5000. As the readme file says, the following scripts should be changed to the other one. May I know why the inlet and outlet types changed to be v and O, respectively? wall { type wall; nFaces 2400; startFace 88300; } inlet { type *patch;* nFaces 500; startFace 90700; } outlet { type *patch;* nFaces 500; startFace 91200; } Because foamToNek will read the boundary condition type in this file. Please change the boundary file to: wall { type W; nFaces 2400; startFace 88300 } inlet { type *v;* nFaces 500; startFace 90700; } outlet { type *O;* nFaces 500; startFace 91200; } Kind regards, Jian -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Tue Jan 3 07:36:21 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 3 Jan 2017 14:36:21 +0100 Subject: [Nek5000-users] convert openfoam mesh to nek5000 In-Reply-To: References: Message-ID: You may want to ask the author of the tool. This converter is not part of the Nek5000 project. On 3 Jan 2017, at 14:15, "nek5000-users at lists.mcs.anl.gov " > wrote: Hi all, I want to convert openfoam mesh to nek5000. As the readme file says, the following scripts should be changed to the other one. May I know why the inlet and outlet types changed to be v and O, respectively? wall { type wall; nFaces 2400; startFace 88300; } inlet { type patch; nFaces 500; startFace 90700; } outlet { type patch; nFaces 500; startFace 91200; } Because foamToNek will read the boundary condition type in this file. Please change the boundary file to: wall { type W; nFaces 2400; startFace 88300 } inlet { type v; nFaces 500; startFace 90700; } outlet { type O; nFaces 500; startFace 91200; } Kind regards, Jian _______________________________________________ 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 Jan 3 05:39:51 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 3 Jan 2017 11:39:51 +0000 Subject: [Nek5000-users] convert openfoam mesh to nek5000 In-Reply-To: References: Message-ID: Jian, ?v? and ?O? are letters to designate boundary types in Nek5000. ?v? is a user defined Dirichlet velocity and ?O? is an outflow boundary condition. On Jan 3, 2017, at 6:12 AM, nek5000-users at lists.mcs.anl.gov wrote: Hi all, I want to convert openfoam mesh to nek5000. As the readme file says, the following scripts should be changed to the other one. May I know why the inlet and outlet types changed to be v and O, respectively? wall { type wall; nFaces 2400; startFace 88300; } inlet { type patch; nFaces 500; startFace 90700; } outlet { type patch; nFaces 500; startFace 91200; } Because foamToNek will read the boundary condition type in this file. Please change the boundary file to: wall { type W; nFaces 2400; startFace 88300 } inlet { type v; nFaces 500; startFace 90700; } outlet { type O; nFaces 500; startFace 91200; } Kind regards, Jian _______________________________________________ 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 Jan 3 05:57:20 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 3 Jan 2017 11:57:20 +0000 Subject: [Nek5000-users] convert openfoam mesh to nek5000 Message-ID: Hi, thanks for your reply. What confuses me a lot is why both types of inlet and outlet in openfoam file were changed to be v and O, since they are same patch type. wall { type wall; nFaces 2400; startFace 88300; } inlet { type *patch;* nFaces 500; startFace 90700; } outlet { type *patch;* nFaces 500; startFace 91200; } Because foamToNek will read the boundary condition type in this file. Please change the boundary file to: wall { type W; nFaces 2400; startFace 88300 } inlet { type *v;* nFaces 500; startFace 90700; } outlet { type *O;* nFaces 500; startFace 91200; } Kind regards, Jian -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Tue Jan 3 07:12:38 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 03 Jan 2017 14:12:38 +0100 Subject: [Nek5000-users] problem with 2d axisymmetric numerical simulation Message-ID: Hi Neks, I am trying to simulate a circular pipe. To start with, I am trying to perform the 2d axisymmetric simulation. I have generated my mesh using gmsh, and then converted it to nek5000 format. The 2d numerical simulation is properly working, however, I am finding some troubles with the axisymmetric simulation. When I try to run my simulation, I find the following message, and the code stop running: read .rea file nelgt/nelgv/lelt: 1192 1192 200 lx1 /lx2 /lx3 : 6 4 4 mapping elements to processors 0 149 149 1192 1192 NELV 1 149 149 1192 1192 NELV 2 149 149 1192 1192 NELV 3 149 149 1192 1192 NELV 4 149 149 1192 1192 NELV 5 149 149 1192 1192 NELV 6 149 149 1192 1192 NELV 7 149 149 1192 1192 NELV RANK 0 IEG 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 602 603 604 605 606 607 608 609 610 611 612 620 621 622 623 624 625 626 627 628 639 640 641 642 643 644 657 658 659 660 676 677 element load imbalance: 0 149 149 done :: mapping elements to processors 0 objects found done :: read .rea file 0.25710E-01 sec setup mesh topology Element face on the axis of symmetry must be FACE 1 Element 6 face 2 is on the axis. Element face on the axis of symmetry must be FACE 1 Element 1 face 2 is on the axis. Element face on the axis of symmetry must be FACE 1 Element 1 face 3 is on the axis. Right-handed check complete for 1192 elements. OK. Element face on the axis of symmetry must be FACE 1 Element 2 face 2 is on the axis. Element face on the axis of symmetry must be FACE 1 Element 1 face 2 is on the axis. Following this post: http://lists.mcs.anl.gov/pipermail/nek5000-users/2012-August/001906.html 1) I have set T IFAXIS T IFAZIV 2) I have set 'A ' BC in the region of y=0. 3) I have copied the FLUID BC to the region of THERMAL BOUNDARY CONDITIONS, modifying W, v and O by t 4) In the .usr file I have defined: c----------------------------------------------------------------------- > subroutine userf (ix,iy,iz,eg) > include 'SIZE' > include 'TOTAL' > include 'NEKUSE' > integer e,f,eg > if(y.gt.0) ffy = temp*temp/y > return > end > c----------------------------------------------------------------------- > subroutine userq (ix,iy,iz,eg) > include 'SIZE' > include 'TOTAL' > include 'NEKUSE' > integer e,f,eg > if(y.gt.0) then > visc = param(2) > qvol = -uy*temp/y > endif > return > end > c----------------------------------------------------------------------- > subroutine useric (ix,iy,iz,ieg) > include 'SIZE' > include 'TOTAL' > include 'NEKUSE' > temp = y > return > end > c----------------------------------------------------------------------- > Please, could you give me a piece of advise with this numerical simulation? Thank you very much in advance. SL From nek5000-users at lists.mcs.anl.gov Tue Jan 3 07:38:41 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 3 Jan 2017 13:38:41 +0000 Subject: [Nek5000-users] problem with 2d axisymmetric numerical simulation In-Reply-To: References: Message-ID: Dear SL, For axisymmetric, if the domain contains the axis, then only points along Edge 1 are allowed to be on the axis. That means that the reference domain, shown here in the preprocessor notation, 3 +-------------+ | | 4 | | 2 | | +-------------+ 1 must have Edge 1 on the axis and that the BC for edge 1 must be "A " If you generate your domain with genbox, this is easily accomplished with something like the following in a case.box file: your.rea 2 spatial dimension -1 number of fields # # comments # # #======================================================== # Box -6 -4 nelx,nely,nelz for Box) 0 6 1.0 x0 x1 ratio 0 1 0.7 y0 y1 ratio P ,P ,A ,f , bc's 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: Tuesday, January 03, 2017 7:12 AM To: Nek5000 Users Subject: [Nek5000-users] problem with 2d axisymmetric numerical simulation Hi Neks, I am trying to simulate a circular pipe. To start with, I am trying to perform the 2d axisymmetric simulation. I have generated my mesh using gmsh, and then converted it to nek5000 format. The 2d numerical simulation is properly working, however, I am finding some troubles with the axisymmetric simulation. When I try to run my simulation, I find the following message, and the code stop running: read .rea file nelgt/nelgv/lelt: 1192 1192 200 lx1 /lx2 /lx3 : 6 4 4 mapping elements to processors 0 149 149 1192 1192 NELV 1 149 149 1192 1192 NELV 2 149 149 1192 1192 NELV 3 149 149 1192 1192 NELV 4 149 149 1192 1192 NELV 5 149 149 1192 1192 NELV 6 149 149 1192 1192 NELV 7 149 149 1192 1192 NELV RANK 0 IEG 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 602 603 604 605 606 607 608 609 610 611 612 620 621 622 623 624 625 626 627 628 639 640 641 642 643 644 657 658 659 660 676 677 element load imbalance: 0 149 149 done :: mapping elements to processors 0 objects found done :: read .rea file 0.25710E-01 sec setup mesh topology Element face on the axis of symmetry must be FACE 1 Element 6 face 2 is on the axis. Element face on the axis of symmetry must be FACE 1 Element 1 face 2 is on the axis. Element face on the axis of symmetry must be FACE 1 Element 1 face 3 is on the axis. Right-handed check complete for 1192 elements. OK. Element face on the axis of symmetry must be FACE 1 Element 2 face 2 is on the axis. Element face on the axis of symmetry must be FACE 1 Element 1 face 2 is on the axis. Following this post: http://lists.mcs.anl.gov/pipermail/nek5000-users/2012-August/001906.html 1) I have set T IFAXIS T IFAZIV 2) I have set 'A ' BC in the region of y=0. 3) I have copied the FLUID BC to the region of THERMAL BOUNDARY CONDITIONS, modifying W, v and O by t 4) In the .usr file I have defined: c----------------------------------------------------------------------- > subroutine userf (ix,iy,iz,eg) > include 'SIZE' > include 'TOTAL' > include 'NEKUSE' > integer e,f,eg > if(y.gt.0) ffy = temp*temp/y > return > end > c----------------------------------------------------------------------- > subroutine userq (ix,iy,iz,eg) > include 'SIZE' > include 'TOTAL' > include 'NEKUSE' > integer e,f,eg > if(y.gt.0) then > visc = param(2) > qvol = -uy*temp/y > endif > return > end > c----------------------------------------------------------------------- > subroutine useric (ix,iy,iz,ieg) > include 'SIZE' > include 'TOTAL' > include 'NEKUSE' > temp = y > return > end > c----------------------------------------------------------------------- > Please, could you give me a piece of advise with this numerical simulation? Thank you very much in advance. SL _______________________________________________ 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 Jan 3 10:07:06 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 03 Jan 2017 17:07:06 +0100 Subject: [Nek5000-users] problem with 2d axisymmetric numerical simulation In-Reply-To: References: Message-ID: Dear Paul, thank you very much for your rapid response. I have used 'A ' BC in the edge 1, however the problem persists... I am wondering if the problem could be in how I have defined my geometry: Is it necessary to define edge 1 as x \in [0,10]? Because I have defined x \in [-2,8]... Is it possible to generate edge 1 as a combination of several subdomains? Because I have divided the geometry in 3 subdomains... Thank you in advance. SL El 03-01-2017 14:38, nek5000-users at lists.mcs.anl.gov escribi?: > Dear SL, > > For axisymmetric, if the domain contains the axis, then only points > along Edge 1 > are allowed to be on the axis. That means that the reference domain, > shown > here in the preprocessor notation, > > 3 > +-------------+ > | | > 4 | | 2 > | | > +-------------+ > 1 > > must have Edge 1 on the axis and that the BC for edge 1 must be "A " > > If you generate your domain with genbox, this is easily accomplished > with something like > the following in a case.box file: > > your.rea > 2 spatial dimension > -1 number of fields > # > # comments > # > # > #======================================================== > # > Box > -6 -4 nelx,nely,nelz for Box) > 0 6 1.0 x0 x1 ratio > 0 1 0.7 y0 y1 ratio > P ,P ,A ,f , bc's > > > 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: Tuesday, January 03, 2017 7:12 AM > To: Nek5000 Users > Subject: [Nek5000-users] problem with 2d axisymmetric numerical > simulation > > Hi Neks, > > I am trying to simulate a circular pipe. To start with, I am trying to > perform the 2d axisymmetric simulation. I have generated my mesh using > gmsh, and then converted it to nek5000 format. The 2d numerical > simulation is properly working, however, I am finding some troubles > with > the axisymmetric simulation. > > When I try to run my simulation, I find the following message, and the > code stop running: > read .rea file > nelgt/nelgv/lelt: 1192 1192 200 > lx1 /lx2 /lx3 : 6 4 4 > > mapping elements to processors > 0 149 149 1192 1192 NELV > 1 149 149 1192 1192 NELV > 2 149 149 1192 1192 NELV > 3 149 149 1192 1192 NELV > 4 149 149 1192 1192 NELV > 5 149 149 1192 1192 NELV > 6 149 149 1192 1192 NELV > 7 149 149 1192 1192 NELV > RANK 0 IEG 474 475 476 477 478 479 > 480 > 481 > 482 483 484 485 486 487 > 488 > 489 > 490 492 493 494 495 496 > 497 > 498 > 499 500 501 502 503 504 > 505 > 506 > 507 508 510 511 512 513 > 514 > 515 > 516 517 518 519 520 521 > 522 > 523 > 524 525 526 528 529 530 > 531 > 532 > 533 534 535 536 537 538 > 539 > 540 > 541 542 543 544 546 547 > 548 > 549 > 550 551 552 553 554 555 > 556 > 557 > 558 559 560 561 562 565 > 566 > 567 > 568 569 570 571 572 573 > 574 > 575 > 576 577 578 579 580 583 > 584 > 585 > 586 587 588 589 590 591 > 592 > 593 > 594 595 596 597 598 602 > 603 > 604 > 605 606 607 608 609 610 > 611 > 612 > 620 621 622 623 624 625 > 626 > 627 > 628 639 640 641 642 643 > 644 > 657 > 658 659 660 676 677 > element load imbalance: 0 149 149 > done :: mapping elements to processors > > 0 objects found > done :: read .rea file 0.25710E-01 sec > > setup mesh topology > Element face on the axis of symmetry must be FACE 1 > Element 6 face 2 is on the axis. > Element face on the axis of symmetry must be FACE 1 > Element 1 face 2 is on the axis. > Element face on the axis of symmetry must be FACE 1 > Element 1 face 3 is on the axis. > Right-handed check complete for 1192 elements. OK. > Element face on the axis of symmetry must be FACE 1 > Element 2 face 2 is on the axis. > Element face on the axis of symmetry must be FACE 1 > Element 1 face 2 is on the axis. > > > > Following this post: > http://lists.mcs.anl.gov/pipermail/nek5000-users/2012-August/001906.html > > 1) I have set > T IFAXIS > T IFAZIV > 2) I have set 'A ' BC in the region of y=0. > 3) I have copied the FLUID BC to the region of THERMAL BOUNDARY > CONDITIONS, modifying > W, v and O by t > 4) In the .usr file I have defined: > c----------------------------------------------------------------------- >> subroutine userf (ix,iy,iz,eg) >> include 'SIZE' >> include 'TOTAL' >> include 'NEKUSE' >> integer e,f,eg >> if(y.gt.0) ffy = temp*temp/y >> return >> end >> c----------------------------------------------------------------------- >> subroutine userq (ix,iy,iz,eg) >> include 'SIZE' >> include 'TOTAL' >> include 'NEKUSE' >> integer e,f,eg >> if(y.gt.0) then >> visc = param(2) >> qvol = -uy*temp/y >> endif >> return >> end >> c----------------------------------------------------------------------- >> subroutine useric (ix,iy,iz,ieg) >> include 'SIZE' >> include 'TOTAL' >> include 'NEKUSE' >> temp = y >> return >> end >> c----------------------------------------------------------------------- >> > > Please, could you give me a piece of advise with this numerical > simulation? > > Thank you very much in advance. > > SL > > _______________________________________________ > 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 Jan 3 11:26:14 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 3 Jan 2017 17:26:14 +0000 Subject: [Nek5000-users] problem with 2d axisymmetric numerical simulation In-Reply-To: References: , Message-ID: Hi SL, Attached is an axisymmetric example that I've just set up. Several subdomains are also possible, but you must pay attention to the "A " bc. Please let me know if the attached works and gives the right answer. It looked ok, but I confess I didn't get out a paper and pencil to verify that the answer is indeed correct.... 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, January 03, 2017 10:07 AM To: nek5000-users at lists.mcs.anl.gov Subject: Re: [Nek5000-users] problem with 2d axisymmetric numerical simulation Dear Paul, thank you very much for your rapid response. I have used 'A ' BC in the edge 1, however the problem persists... I am wondering if the problem could be in how I have defined my geometry: Is it necessary to define edge 1 as x \in [0,10]? Because I have defined x \in [-2,8]... Is it possible to generate edge 1 as a combination of several subdomains? Because I have divided the geometry in 3 subdomains... Thank you in advance. SL El 03-01-2017 14:38, nek5000-users at lists.mcs.anl.gov escribi?: > Dear SL, > > For axisymmetric, if the domain contains the axis, then only points > along Edge 1 > are allowed to be on the axis. That means that the reference domain, > shown > here in the preprocessor notation, > > 3 > +-------------+ > | | > 4 | | 2 > | | > +-------------+ > 1 > > must have Edge 1 on the axis and that the BC for edge 1 must be "A " > > If you generate your domain with genbox, this is easily accomplished > with something like > the following in a case.box file: > > your.rea > 2 spatial dimension > -1 number of fields > # > # comments > # > # > #======================================================== > # > Box > -6 -4 nelx,nely,nelz for Box) > 0 6 1.0 x0 x1 ratio > 0 1 0.7 y0 y1 ratio > P ,P ,A ,f , bc's > > > 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: Tuesday, January 03, 2017 7:12 AM > To: Nek5000 Users > Subject: [Nek5000-users] problem with 2d axisymmetric numerical > simulation > > Hi Neks, > > I am trying to simulate a circular pipe. To start with, I am trying to > perform the 2d axisymmetric simulation. I have generated my mesh using > gmsh, and then converted it to nek5000 format. The 2d numerical > simulation is properly working, however, I am finding some troubles > with > the axisymmetric simulation. > > When I try to run my simulation, I find the following message, and the > code stop running: > read .rea file > nelgt/nelgv/lelt: 1192 1192 200 > lx1 /lx2 /lx3 : 6 4 4 > > mapping elements to processors > 0 149 149 1192 1192 NELV > 1 149 149 1192 1192 NELV > 2 149 149 1192 1192 NELV > 3 149 149 1192 1192 NELV > 4 149 149 1192 1192 NELV > 5 149 149 1192 1192 NELV > 6 149 149 1192 1192 NELV > 7 149 149 1192 1192 NELV > RANK 0 IEG 474 475 476 477 478 479 > 480 > 481 > 482 483 484 485 486 487 > 488 > 489 > 490 492 493 494 495 496 > 497 > 498 > 499 500 501 502 503 504 > 505 > 506 > 507 508 510 511 512 513 > 514 > 515 > 516 517 518 519 520 521 > 522 > 523 > 524 525 526 528 529 530 > 531 > 532 > 533 534 535 536 537 538 > 539 > 540 > 541 542 543 544 546 547 > 548 > 549 > 550 551 552 553 554 555 > 556 > 557 > 558 559 560 561 562 565 > 566 > 567 > 568 569 570 571 572 573 > 574 > 575 > 576 577 578 579 580 583 > 584 > 585 > 586 587 588 589 590 591 > 592 > 593 > 594 595 596 597 598 602 > 603 > 604 > 605 606 607 608 609 610 > 611 > 612 > 620 621 622 623 624 625 > 626 > 627 > 628 639 640 641 642 643 > 644 > 657 > 658 659 660 676 677 > element load imbalance: 0 149 149 > done :: mapping elements to processors > > 0 objects found > done :: read .rea file 0.25710E-01 sec > > setup mesh topology > Element face on the axis of symmetry must be FACE 1 > Element 6 face 2 is on the axis. > Element face on the axis of symmetry must be FACE 1 > Element 1 face 2 is on the axis. > Element face on the axis of symmetry must be FACE 1 > Element 1 face 3 is on the axis. > Right-handed check complete for 1192 elements. OK. > Element face on the axis of symmetry must be FACE 1 > Element 2 face 2 is on the axis. > Element face on the axis of symmetry must be FACE 1 > Element 1 face 2 is on the axis. > > > > Following this post: > http://lists.mcs.anl.gov/pipermail/nek5000-users/2012-August/001906.html > > 1) I have set > T IFAXIS > T IFAZIV > 2) I have set 'A ' BC in the region of y=0. > 3) I have copied the FLUID BC to the region of THERMAL BOUNDARY > CONDITIONS, modifying > W, v and O by t > 4) In the .usr file I have defined: > c----------------------------------------------------------------------- >> subroutine userf (ix,iy,iz,eg) >> include 'SIZE' >> include 'TOTAL' >> include 'NEKUSE' >> integer e,f,eg >> if(y.gt.0) ffy = temp*temp/y >> return >> end >> c----------------------------------------------------------------------- >> subroutine userq (ix,iy,iz,eg) >> include 'SIZE' >> include 'TOTAL' >> include 'NEKUSE' >> integer e,f,eg >> if(y.gt.0) then >> visc = param(2) >> qvol = -uy*temp/y >> endif >> return >> end >> c----------------------------------------------------------------------- >> subroutine useric (ix,iy,iz,ieg) >> include 'SIZE' >> include 'TOTAL' >> include 'NEKUSE' >> temp = y >> return >> end >> c----------------------------------------------------------------------- >> > > Please, could you give me a piece of advise with this numerical > simulation? > > Thank you very much in advance. > > SL > > _______________________________________________ > 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 -------------- A non-text attachment was scrubbed... Name: axipipe.tar Type: application/x-tar Size: 21504 bytes Desc: axipipe.tar URL: From nek5000-users at lists.mcs.anl.gov Tue Jan 3 12:04:13 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 03 Jan 2017 19:04:13 +0100 Subject: [Nek5000-users] problem with 2d axisymmetric numerical simulation In-Reply-To: References: , Message-ID: Hi Paul, thank you very much, it works! However, I am trying to re-define the geometry generated by genbox: # Box -12 -4 nelx,nely,nelz for Box) -2.0 4.0 1.0 x0 x1 ratio 0.0 1.0 0.6 y0 y1 ratio P ,P ,A ,W , bc's The idea is to generate 3 different squared domains: ----------- | | -------- y1 | |----| y2 | | ------------------------| x1 x2 So, I have included in usrdat2 subroutine: do i=1,ntot if (xm1(i,1,1,1).lt.x1) then ym1(i,1,1,1)=y1 elseif (xm1(i,1,1,1).ge.x1.and.xm1(i,1,1,1).lt.x2) then ym2(i,1,1,1)=y2 end if end do However, when I run Nek5000 I obtain the following error... call usrdat2 done :: usrdat2 regenerate geometry data 1 0 ERROR: Vanishing Jacobian near 5th node of element 4. -2.0182169769993447E-015 0.18749999999999875 0 xyz: -1.20930E+00 5.00000E-01 0 xyz: -7.90701E-01 6.25000E-01 Is there any different way to create these several domains using genbox? Thank you SL El 03-01-2017 18:26, nek5000-users at lists.mcs.anl.gov escribi?: > Hi SL, > > Attached is an axisymmetric example that I've just set up. > > Several subdomains are also possible, but you must pay attention > to the "A " bc. > > Please let me know if the attached works and gives the right answer. > It looked ok, but I confess I didn't get out a paper and pencil to > verify > that the answer is indeed correct.... > > 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, January 03, 2017 10:07 AM > To: nek5000-users at lists.mcs.anl.gov > Subject: Re: [Nek5000-users] problem with 2d axisymmetric numerical > simulation > > Dear Paul, > > thank you very much for your rapid response. > > I have used 'A ' BC in the edge 1, however the problem persists... I am > wondering if the problem could be in how I have defined my geometry: > > Is it necessary to define edge 1 as x \in [0,10]? Because I have > defined > x \in [-2,8]... > Is it possible to generate edge 1 as a combination of several > subdomains? Because I have divided the geometry in 3 subdomains... > > Thank you in advance. > SL > > El 03-01-2017 14:38, nek5000-users at lists.mcs.anl.gov escribi?: >> Dear SL, >> >> For axisymmetric, if the domain contains the axis, then only points >> along Edge 1 >> are allowed to be on the axis. That means that the reference domain, >> shown >> here in the preprocessor notation, >> >> 3 >> +-------------+ >> | | >> 4 | | 2 >> | | >> +-------------+ >> 1 >> >> must have Edge 1 on the axis and that the BC for edge 1 must be "A " >> >> If you generate your domain with genbox, this is easily accomplished >> with something like >> the following in a case.box file: >> >> your.rea >> 2 spatial dimension >> -1 number of fields >> # >> # comments >> # >> # >> #======================================================== >> # >> Box >> -6 -4 nelx,nely,nelz for Box) >> 0 6 1.0 x0 x1 ratio >> 0 1 0.7 y0 y1 ratio >> P ,P ,A ,f , bc's >> >> >> 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: Tuesday, January 03, 2017 7:12 AM >> To: Nek5000 Users >> Subject: [Nek5000-users] problem with 2d axisymmetric numerical >> simulation >> >> Hi Neks, >> >> I am trying to simulate a circular pipe. To start with, I am trying to >> perform the 2d axisymmetric simulation. I have generated my mesh using >> gmsh, and then converted it to nek5000 format. The 2d numerical >> simulation is properly working, however, I am finding some troubles >> with >> the axisymmetric simulation. >> >> When I try to run my simulation, I find the following message, and the >> code stop running: >> read .rea file >> nelgt/nelgv/lelt: 1192 1192 200 >> lx1 /lx2 /lx3 : 6 4 4 >> >> mapping elements to processors >> 0 149 149 1192 1192 NELV >> 1 149 149 1192 1192 NELV >> 2 149 149 1192 1192 NELV >> 3 149 149 1192 1192 NELV >> 4 149 149 1192 1192 NELV >> 5 149 149 1192 1192 NELV >> 6 149 149 1192 1192 NELV >> 7 149 149 1192 1192 NELV >> RANK 0 IEG 474 475 476 477 478 479 >> 480 >> 481 >> 482 483 484 485 486 487 >> 488 >> 489 >> 490 492 493 494 495 496 >> 497 >> 498 >> 499 500 501 502 503 504 >> 505 >> 506 >> 507 508 510 511 512 513 >> 514 >> 515 >> 516 517 518 519 520 521 >> 522 >> 523 >> 524 525 526 528 529 530 >> 531 >> 532 >> 533 534 535 536 537 538 >> 539 >> 540 >> 541 542 543 544 546 547 >> 548 >> 549 >> 550 551 552 553 554 555 >> 556 >> 557 >> 558 559 560 561 562 565 >> 566 >> 567 >> 568 569 570 571 572 573 >> 574 >> 575 >> 576 577 578 579 580 583 >> 584 >> 585 >> 586 587 588 589 590 591 >> 592 >> 593 >> 594 595 596 597 598 602 >> 603 >> 604 >> 605 606 607 608 609 610 >> 611 >> 612 >> 620 621 622 623 624 625 >> 626 >> 627 >> 628 639 640 641 642 643 >> 644 >> 657 >> 658 659 660 676 677 >> element load imbalance: 0 149 149 >> done :: mapping elements to processors >> >> 0 objects found >> done :: read .rea file 0.25710E-01 sec >> >> setup mesh topology >> Element face on the axis of symmetry must be FACE 1 >> Element 6 face 2 is on the axis. >> Element face on the axis of symmetry must be FACE 1 >> Element 1 face 2 is on the axis. >> Element face on the axis of symmetry must be FACE 1 >> Element 1 face 3 is on the axis. >> Right-handed check complete for 1192 elements. OK. >> Element face on the axis of symmetry must be FACE 1 >> Element 2 face 2 is on the axis. >> Element face on the axis of symmetry must be FACE 1 >> Element 1 face 2 is on the axis. >> >> >> >> Following this post: >> http://lists.mcs.anl.gov/pipermail/nek5000-users/2012-August/001906.html >> >> 1) I have set >> T IFAXIS >> T IFAZIV >> 2) I have set 'A ' BC in the region of y=0. >> 3) I have copied the FLUID BC to the region of THERMAL BOUNDARY >> CONDITIONS, modifying >> W, v and O by t >> 4) In the .usr file I have defined: >> c----------------------------------------------------------------------- >>> subroutine userf (ix,iy,iz,eg) >>> include 'SIZE' >>> include 'TOTAL' >>> include 'NEKUSE' >>> integer e,f,eg >>> if(y.gt.0) ffy = temp*temp/y >>> return >>> end >>> c----------------------------------------------------------------------- >>> subroutine userq (ix,iy,iz,eg) >>> include 'SIZE' >>> include 'TOTAL' >>> include 'NEKUSE' >>> integer e,f,eg >>> if(y.gt.0) then >>> visc = param(2) >>> qvol = -uy*temp/y >>> endif >>> return >>> end >>> c----------------------------------------------------------------------- >>> subroutine useric (ix,iy,iz,ieg) >>> include 'SIZE' >>> include 'TOTAL' >>> include 'NEKUSE' >>> temp = y >>> return >>> end >>> c----------------------------------------------------------------------- >>> >> >> Please, could you give me a piece of advise with this numerical >> simulation? >> >> Thank you very much in advance. >> >> SL >> >> _______________________________________________ >> 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 Jan 5 10:40:11 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 5 Jan 2017 17:40:11 +0100 Subject: [Nek5000-users] gradient of temperature in conj_ht Message-ID: Hello all, I computed the gradient of temperature in conjugate heat transfer case. The gradient of temperature is being computed correctly for the fluid domain (I validated analytically). But, in the solid domain grad T is showing to be zero, which is incorrect. Can anyone please tell me how to compute gradient of a scalar in solid domain for conj_ht case. Is it possible to compute the gradient of a scalar in this case? The scalar is C0 and not differentiable at the interface between solid and fluid domain. I added the following in the userchk subroutine parameter (lt=lx1*ly1*lz1*lelt) common /mystuff/ dtdx(lt),dtdy(lt),dtdz(lt) common /myvort/ cgtx(lt),cgty(lt),cgtz3(lt) call gradm1(dtdx,dtdy,dtdz,t) call outpost(dtdx,dtdy,dtdz,pr,t,'gdt') ! write to gdtblah.f000... call comp_vort3(cgtx,cgty,cgtz,dtdx,dtdy,dtdz) I attach usr, rea and box files for reference. Thank you in advance. Yours sincerely, Sandeep -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: conj_ht.box Type: application/octet-stream Size: 640 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: conj_ht.usr Type: application/octet-stream Size: 4067 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: conj_ht.rea Type: application/octet-stream Size: 130535 bytes Desc: not available URL: From nek5000-users at lists.mcs.anl.gov Thu Jan 5 12:40:54 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 5 Jan 2017 18:40:54 +0000 Subject: [Nek5000-users] gradient of temperature in conj_ht In-Reply-To: References: Message-ID: Hi Sandeep, My guess is that prepost/outpost is not writing the vector field in the solid elements. Using outpost in this manner, which in fact is what I would suggest, the code does not realize that you are writing a field that should be nonzero in the solid. I'll try to think of a work around. One quick one that might work: nelf = nelv ! Save number of fluid elements nelfld(1) = nelt nelv = nelt call outpost(tx,ty,tz,pr,t,'gdt') nelv = nelf nelfld(1) = nelf ... this might be sufficient to trigger all the flags that prepost may look for. 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, January 05, 2017 10:40 AM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] gradient of temperature in conj_ht Hello all, I computed the gradient of temperature in conjugate heat transfer case. The gradient of temperature is being computed correctly for the fluid domain (I validated analytically). But, in the solid domain grad T is showing to be zero, which is incorrect. Can anyone please tell me how to compute gradient of a scalar in solid domain for conj_ht case. Is it possible to compute the gradient of a scalar in this case? The scalar is C0 and not differentiable at the interface between solid and fluid domain. I added the following in the userchk subroutine parameter (lt=lx1*ly1*lz1*lelt) common /mystuff/ dtdx(lt),dtdy(lt),dtdz(lt) common /myvort/ cgtx(lt),cgty(lt),cgtz3(lt) call gradm1(dtdx,dtdy,dtdz,t) call outpost(dtdx,dtdy,dtdz,pr,t,'gdt') ! write to gdtblah.f000... call comp_vort3(cgtx,cgty,cgtz,dtdx,dtdy,dtdz) I attach usr, rea and box files for reference. Thank you in advance. Yours sincerely, Sandeep -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Thu Jan 5 14:09:08 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 5 Jan 2017 20:09:08 +0000 Subject: [Nek5000-users] problem with 2d axisymmetric numerical simulation In-Reply-To: References: , , Message-ID: Hi SL, Attached is a tar file with an axisymmetric expansion case and a couple of README files. Please give this a try -- hope this helps you to understand how genbox can do this. 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, January 03, 2017 12:04 PM To: nek5000-users at lists.mcs.anl.gov Subject: Re: [Nek5000-users] problem with 2d axisymmetric numerical simulation Hi Paul, thank you very much, it works! However, I am trying to re-define the geometry generated by genbox: # Box -12 -4 nelx,nely,nelz for Box) -2.0 4.0 1.0 x0 x1 ratio 0.0 1.0 0.6 y0 y1 ratio P ,P ,A ,W , bc's The idea is to generate 3 different squared domains: ----------- | | -------- y1 | |----| y2 | | ------------------------| x1 x2 So, I have included in usrdat2 subroutine: do i=1,ntot if (xm1(i,1,1,1).lt.x1) then ym1(i,1,1,1)=y1 elseif (xm1(i,1,1,1).ge.x1.and.xm1(i,1,1,1).lt.x2) then ym2(i,1,1,1)=y2 end if end do However, when I run Nek5000 I obtain the following error... call usrdat2 done :: usrdat2 regenerate geometry data 1 0 ERROR: Vanishing Jacobian near 5th node of element 4. -2.0182169769993447E-015 0.18749999999999875 0 xyz: -1.20930E+00 5.00000E-01 0 xyz: -7.90701E-01 6.25000E-01 Is there any different way to create these several domains using genbox? Thank you SL El 03-01-2017 18:26, nek5000-users at lists.mcs.anl.gov escribi?: > Hi SL, > > Attached is an axisymmetric example that I've just set up. > > Several subdomains are also possible, but you must pay attention > to the "A " bc. > > Please let me know if the attached works and gives the right answer. > It looked ok, but I confess I didn't get out a paper and pencil to > verify > that the answer is indeed correct.... > > 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, January 03, 2017 10:07 AM > To: nek5000-users at lists.mcs.anl.gov > Subject: Re: [Nek5000-users] problem with 2d axisymmetric numerical > simulation > > Dear Paul, > > thank you very much for your rapid response. > > I have used 'A ' BC in the edge 1, however the problem persists... I am > wondering if the problem could be in how I have defined my geometry: > > Is it necessary to define edge 1 as x \in [0,10]? Because I have > defined > x \in [-2,8]... > Is it possible to generate edge 1 as a combination of several > subdomains? Because I have divided the geometry in 3 subdomains... > > Thank you in advance. > SL > > El 03-01-2017 14:38, nek5000-users at lists.mcs.anl.gov escribi?: >> Dear SL, >> >> For axisymmetric, if the domain contains the axis, then only points >> along Edge 1 >> are allowed to be on the axis. That means that the reference domain, >> shown >> here in the preprocessor notation, >> >> 3 >> +-------------+ >> | | >> 4 | | 2 >> | | >> +-------------+ >> 1 >> >> must have Edge 1 on the axis and that the BC for edge 1 must be "A " >> >> If you generate your domain with genbox, this is easily accomplished >> with something like >> the following in a case.box file: >> >> your.rea >> 2 spatial dimension >> -1 number of fields >> # >> # comments >> # >> # >> #======================================================== >> # >> Box >> -6 -4 nelx,nely,nelz for Box) >> 0 6 1.0 x0 x1 ratio >> 0 1 0.7 y0 y1 ratio >> P ,P ,A ,f , bc's >> >> >> 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: Tuesday, January 03, 2017 7:12 AM >> To: Nek5000 Users >> Subject: [Nek5000-users] problem with 2d axisymmetric numerical >> simulation >> >> Hi Neks, >> >> I am trying to simulate a circular pipe. To start with, I am trying to >> perform the 2d axisymmetric simulation. I have generated my mesh using >> gmsh, and then converted it to nek5000 format. The 2d numerical >> simulation is properly working, however, I am finding some troubles >> with >> the axisymmetric simulation. >> >> When I try to run my simulation, I find the following message, and the >> code stop running: >> read .rea file >> nelgt/nelgv/lelt: 1192 1192 200 >> lx1 /lx2 /lx3 : 6 4 4 >> >> mapping elements to processors >> 0 149 149 1192 1192 NELV >> 1 149 149 1192 1192 NELV >> 2 149 149 1192 1192 NELV >> 3 149 149 1192 1192 NELV >> 4 149 149 1192 1192 NELV >> 5 149 149 1192 1192 NELV >> 6 149 149 1192 1192 NELV >> 7 149 149 1192 1192 NELV >> RANK 0 IEG 474 475 476 477 478 479 >> 480 >> 481 >> 482 483 484 485 486 487 >> 488 >> 489 >> 490 492 493 494 495 496 >> 497 >> 498 >> 499 500 501 502 503 504 >> 505 >> 506 >> 507 508 510 511 512 513 >> 514 >> 515 >> 516 517 518 519 520 521 >> 522 >> 523 >> 524 525 526 528 529 530 >> 531 >> 532 >> 533 534 535 536 537 538 >> 539 >> 540 >> 541 542 543 544 546 547 >> 548 >> 549 >> 550 551 552 553 554 555 >> 556 >> 557 >> 558 559 560 561 562 565 >> 566 >> 567 >> 568 569 570 571 572 573 >> 574 >> 575 >> 576 577 578 579 580 583 >> 584 >> 585 >> 586 587 588 589 590 591 >> 592 >> 593 >> 594 595 596 597 598 602 >> 603 >> 604 >> 605 606 607 608 609 610 >> 611 >> 612 >> 620 621 622 623 624 625 >> 626 >> 627 >> 628 639 640 641 642 643 >> 644 >> 657 >> 658 659 660 676 677 >> element load imbalance: 0 149 149 >> done :: mapping elements to processors >> >> 0 objects found >> done :: read .rea file 0.25710E-01 sec >> >> setup mesh topology >> Element face on the axis of symmetry must be FACE 1 >> Element 6 face 2 is on the axis. >> Element face on the axis of symmetry must be FACE 1 >> Element 1 face 2 is on the axis. >> Element face on the axis of symmetry must be FACE 1 >> Element 1 face 3 is on the axis. >> Right-handed check complete for 1192 elements. OK. >> Element face on the axis of symmetry must be FACE 1 >> Element 2 face 2 is on the axis. >> Element face on the axis of symmetry must be FACE 1 >> Element 1 face 2 is on the axis. >> >> >> >> Following this post: >> http://lists.mcs.anl.gov/pipermail/nek5000-users/2012-August/001906.html >> >> 1) I have set >> T IFAXIS >> T IFAZIV >> 2) I have set 'A ' BC in the region of y=0. >> 3) I have copied the FLUID BC to the region of THERMAL BOUNDARY >> CONDITIONS, modifying >> W, v and O by t >> 4) In the .usr file I have defined: >> c----------------------------------------------------------------------- >>> subroutine userf (ix,iy,iz,eg) >>> include 'SIZE' >>> include 'TOTAL' >>> include 'NEKUSE' >>> integer e,f,eg >>> if(y.gt.0) ffy = temp*temp/y >>> return >>> end >>> c----------------------------------------------------------------------- >>> subroutine userq (ix,iy,iz,eg) >>> include 'SIZE' >>> include 'TOTAL' >>> include 'NEKUSE' >>> integer e,f,eg >>> if(y.gt.0) then >>> visc = param(2) >>> qvol = -uy*temp/y >>> endif >>> return >>> end >>> c----------------------------------------------------------------------- >>> subroutine useric (ix,iy,iz,ieg) >>> include 'SIZE' >>> include 'TOTAL' >>> include 'NEKUSE' >>> temp = y >>> return >>> end >>> c----------------------------------------------------------------------- >>> >> >> Please, could you give me a piece of advise with this numerical >> simulation? >> >> Thank you very much in advance. >> >> SL >> >> _______________________________________________ >> 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 _______________________________________________ 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: t.tar Type: application/x-tar Size: 88064 bytes Desc: t.tar URL: From nek5000-users at lists.mcs.anl.gov Fri Jan 6 05:09:16 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 6 Jan 2017 12:09:16 +0100 Subject: [Nek5000-users] gradient of temperature in conj_ht Message-ID: Hello Paul, Thank you for your reply. It works now. I am able to recover the gradient of temperature to be same as in the analytical solution. Yours sincerely, Sandeep -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Wed Jan 11 15:41:54 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 11 Jan 2017 21:41:54 +0000 Subject: [Nek5000-users] genmap Message-ID: Hi nek, Would you please tell me what's going on with the following errors. The rea file I used comes from the OpenFOAM. reading .rea file data ... start locglob_lexico: 8 67584 540672 0.10000000000000001 locglob: 1 1 540672 locglob: 2 31 540672 locglob: 3 347 540672 locglob: 1 4687 540672 locglob: 2 49135 540672 locglob: 3 70703 540672 locglob: 1 70868 540672 locglob: 2 70913 540672 locglob: 3 70913 540672 done locglob_lexico: 70913 70913 540672 32 MISSING FACE CONNECTION a 2 20155 26123 26 36 19947 19948 19963 19963 19964 19964 20139 20139 20140 20140 20156 20156 20156 20156 20331 20332 20347 20347 20348 20348 25931 25932 26123 26123 26123 26124 26124 26124 26139 26140 26315 26315 26316 26316 26331 26332 MISSING FACE CONNECTION a 2 20155 26123 27 37 19947 19948 19963 19963 19964 19964 20139 20139 20140 20140 20156 20156 20156 20156 20331 20332 20347 20347 20348 20348 25931 25932 26122 26123 26123 26123 26124 26124 26124 26139 26140 26315 26315 26316 26316 26331 26332 MISSING FACE CONNECTION a 6 19963 26123 30 30 19946 19947 19947 19948 19962 19962 19964 19964 20138 20139 20139 20140 20154 20154 20155 20155 20155 20155 20156 20156 25930 25931 25931 25932 26122 26122 26123 26123 26123 26124 MISSING FACE CONNECTION a 2 19962 25931 28 29 19754 19755 19770 19770 19771 19771 19946 19946 19947 19947 19963 19963 19963 19963 20138 20139 20154 20154 20155 20155 25738 25739 25930 25930 25931 25931 25931 26122 26123 MISSING FACE CONNECTION a 2 19754 19755 16 29 19546 19547 19562 19562 19563 19563 19578 19579 19738 19738 19739 19739 19755 19755 19755 19770 19771 19930 19931 19946 19947 25738 25739 25754 25755 25930 25931 25946 25947 WARNING:Missing Face Connection Not Resolved Kind regards, Jian -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Fri Jan 13 05:48:46 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 13 Jan 2017 12:48:46 +0100 Subject: [Nek5000-users] Convergence problem in 2D square cylinder Message-ID: Dear Nek users, I am a PhD student at School of Aeronautics and Astronautics, Technical University of Madrid, Spain. I am really interested in using Nek5000 in order to generate a flow field database. After that, I will apply reduced order modelling techniques. The final goal is to study the reduced order models and to investigate possible ways to speed up the calculations. My first run in Nek5000 is bidimensional flow over a square cylinder for low Reynolds numbers (Re = 100). Sadly, I am getting some convergence problems, Courant number grows up and explodes in a few iterations. Plotting the velocity fields, I have noticed that the velocity near the square cylinder corners is very big compared to what should be expected. Anyone has any recommendation? Is there any way to introduce artificial viscosity or to increment dissipation? Regards, VBeltranM -- V?ctor Beltr?n Mart?nez Ph.D. Student in Applied Mathematics Escuela T?cnica Superior de Ingenier?a Aeron?utica y del Espacio Universidad Polit?cnica de Madrid From nek5000-users at lists.mcs.anl.gov Fri Jan 13 16:10:16 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 13 Jan 2017 15:10:16 -0700 Subject: [Nek5000-users] Convergence problem in 2D square cylinder In-Reply-To: References: Message-ID: Hi Victor, Do you get any improvement by decreasing dt (.i.e. timstep) in the .rea file or playing around with the order of temporal scheme (TORDER) ?. You can introduce some filtering, through parameters 101 and 103, which would mimic the effects of artificial viscosity. You can read more about the filtering by carefully reading the subroutine q_filter(wght) from the file navier5.f. If you really want to go through the artificial viscosity procedure in its literal sense, you can modify the subroutine uservp() (You can go through the example in turbChannel to see how to make modifications) and set param(30) = 1. Best Regards, Tanmoy On Fri, Jan 13, 2017 at 4:48 AM, wrote: > Dear Nek users, > > I am a PhD student at School of Aeronautics and Astronautics, Technical > University of Madrid, Spain. > I am really interested in using Nek5000 in order to generate a flow field > database. After that, I will apply reduced order modelling techniques. The > final goal is to study the reduced order models and to investigate possible > ways to speed up the calculations. > > My first run in Nek5000 is bidimensional flow over a square cylinder for > low Reynolds numbers (Re = 100). Sadly, I am getting some convergence > problems, Courant number grows up and explodes in a few iterations. > Plotting the velocity fields, I have noticed that the velocity near the > square cylinder corners is very big compared to what should be expected. > > Anyone has any recommendation? Is there any way to introduce artificial > viscosity or to increment dissipation? > > Regards, > > VBeltranM > > -- > V?ctor Beltr?n Mart?nez > Ph.D. Student in Applied Mathematics > Escuela T?cnica Superior de Ingenier?a Aeron?utica y del Espacio > Universidad Polit?cnica de Madrid > > _______________________________________________ > 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 Mon Jan 16 05:33:02 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 16 Jan 2017 12:33:02 +0100 Subject: [Nek5000-users] Rank mismatch in array reference in qvol Message-ID: Hello All, I would like to solve dA/dt = Laplacian A + U cross (curl A) 'A' is a combination of three scalar arrays, T, PS1 and PS2. I modified the following to switch off the advection for 'A' in the .rea file 2.0000000E+00 p23 NPSCAL T IFFLOW T IFHEAT T IFTRAN T F F F F F F F F F F IFNAV & IFADVC (convection in P.S. fields) F F T T T F F F F F F F IFTMSH (IF mesh for this field is T mesh) Followed by fluid, T, PS1, and PS2 bcs. In userchk: parameter (lt=lx1*ly1*lz1*lelt) common /myvort/ b1(lt,3),work1(lt),work2(lt) common /vpot/ ucbx(lt),ucby(lt),ucbz(lt) m = nx1*ny1*nz1*nelv call comp_vort3(b1,work1,work2,t(1,1,1,1,1), t(1,1,1,1,2),t(1,1,1,1,3)) call vcross(ucbx,ucby,ucbz,vx,vy,vz,b1(1,1),b1(1,2),b1(1,3),m) In userq(ix,iy,iz,ieg) parameter (lt=lx1*ly1*lz1*lelt) common /vpot/ ucbx(lt),ucby(lt),ucbz(lt) integer e,f,eg e = gllel(eg) if (ifield.eq.2) then qvol = ucbx(ix,iy,iz,ieg) elseif (ifield.eq.3) then qvol = ucby(ix,iy,iz,ieg) elseif (ifield.eq.4) then qvol = ucbz(ix,iy,iz,ieg) endif When I compile i get and error qvol = ucbx(ix,iy,iz,ieg) 1 Error: Rank mismatch in array reference at (1) (4/1) /home/sandeep_reddy/work/Dynamo/GP/vectpot2/conj_ht.f:187.20: qvol = ucby(ix,iy,iz,ieg) 1 Error: Rank mismatch in array reference at (1) (4/1) /home/sandeep_reddy/work/Dynamo/GP/vectpot2/conj_ht.f:190.20: qvol = ucbz(ix,iy,iz,ieg) 1 Error: Rank mismatch in array reference at (1) (4/1) Could someone kindly help me with this? Thanking you. Your sincerely, Sandeep -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Jan 16 08:53:12 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 16 Jan 2017 14:53:12 +0000 Subject: [Nek5000-users] Rank mismatch in array reference in qvol In-Reply-To: References: Message-ID: Dear Sandeep, You need to declare ucbx(lx1,ly1,lz1,lelt), etc. 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 16, 2017 5:33 AM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] Rank mismatch in array reference in qvol Hello All, I would like to solve dA/dt = Laplacian A + U cross (curl A) 'A' is a combination of three scalar arrays, T, PS1 and PS2. I modified the following to switch off the advection for 'A' in the .rea file 2.0000000E+00 p23 NPSCAL T IFFLOW T IFHEAT T IFTRAN T F F F F F F F F F F IFNAV & IFADVC (convection in P.S. fields) F F T T T F F F F F F F IFTMSH (IF mesh for this field is T mesh) Followed by fluid, T, PS1, and PS2 bcs. In userchk: parameter (lt=lx1*ly1*lz1*lelt) common /myvort/ b1(lt,3),work1(lt),work2(lt) common /vpot/ ucbx(lt),ucby(lt),ucbz(lt) m = nx1*ny1*nz1*nelv call comp_vort3(b1,work1,work2,t(1,1,1,1,1), t(1,1,1,1,2),t(1,1,1,1,3)) call vcross(ucbx,ucby,ucbz,vx,vy,vz,b1(1,1),b1(1,2),b1(1,3),m) In userq(ix,iy,iz,ieg) parameter (lt=lx1*ly1*lz1*lelt) common /vpot/ ucbx(lt),ucby(lt),ucbz(lt) integer e,f,eg e = gllel(eg) if (ifield.eq.2) then qvol = ucbx(ix,iy,iz,ieg) elseif (ifield.eq.3) then qvol = ucby(ix,iy,iz,ieg) elseif (ifield.eq.4) then qvol = ucbz(ix,iy,iz,ieg) endif When I compile i get and error qvol = ucbx(ix,iy,iz,ieg) 1 Error: Rank mismatch in array reference at (1) (4/1) /home/sandeep_reddy/work/Dynamo/GP/vectpot2/conj_ht.f:187.20: qvol = ucby(ix,iy,iz,ieg) 1 Error: Rank mismatch in array reference at (1) (4/1) /home/sandeep_reddy/work/Dynamo/GP/vectpot2/conj_ht.f:190.20: qvol = ucbz(ix,iy,iz,ieg) 1 Error: Rank mismatch in array reference at (1) (4/1) Could someone kindly help me with this? Thanking you. Your sincerely, Sandeep -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Jan 16 08:54:27 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 16 Jan 2017 14:54:27 +0000 Subject: [Nek5000-users] Convergence problem in 2D square cylinder In-Reply-To: References: Message-ID: Dear Victor, Do you have sufficient boundary layer resolution near your cylinder? You'll want to have a relatively thin layer of elements surrounding the cylinder. 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 13, 2017 5:48 AM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] Convergence problem in 2D square cylinder Dear Nek users, I am a PhD student at School of Aeronautics and Astronautics, Technical University of Madrid, Spain. I am really interested in using Nek5000 in order to generate a flow field database. After that, I will apply reduced order modelling techniques. The final goal is to study the reduced order models and to investigate possible ways to speed up the calculations. My first run in Nek5000 is bidimensional flow over a square cylinder for low Reynolds numbers (Re = 100). Sadly, I am getting some convergence problems, Courant number grows up and explodes in a few iterations. Plotting the velocity fields, I have noticed that the velocity near the square cylinder corners is very big compared to what should be expected. Anyone has any recommendation? Is there any way to introduce artificial viscosity or to increment dissipation? Regards, VBeltranM -- V?ctor Beltr?n Mart?nez Ph.D. Student in Applied Mathematics Escuela T?cnica Superior de Ingenier?a Aeron?utica y del Espacio Universidad Polit?cnica de Madrid _______________________________________________ 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 16 09:05:21 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 16 Jan 2017 15:05:21 +0000 Subject: [Nek5000-users] double mesh ? In-Reply-To: References: Message-ID: Dear Agnese, Sorry for the long delay in answering this question. The issue is that in fact it is possible (and used to be a supported but rarely used feature). The challenge, however, is how to make it relatively easy in the context of the current setup for a variety of applications. Nek supports a different gather-scatter handle for each field (velocity, temperature, passive scalar, etc.) We have a handle that allows the user access to the connectivity list, which would give you an opportunity to "disconnect" the periodic faces for the temperature array. I'll try to code up an example for the case of channel flow with a mesh generated by genbox and post it shortly. Incidentally, an alternative approach that we often use for more complex domains would be to use a recycling boundary condition. I think it would not be as efficient, however, as periodicity for the case that you propose. Once again, apologies for the delayed response. 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 22, 2016 4:38 AM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] double mesh ? Hi i am quite new, so i?m sorry if i?m getting things wrong. I am trying to simulate a turbulent channel flow, with a passive scalar. in the streamwise direction, i would like to impose periodic boundary conditions for the flow, and non-periodic for the scalar. in fact the scalar has a point source near the inlet, and i want to see how that develops downstream. as far as i understand this is not supported by nek5000, because the last mesh point is identified with the first mesh point when periodic BC are imposed. correct ? my question is : could i define two meshes, one for the scalar and one for the velocity ? thanks a lot for your help !!! agnese Agnese Seminara -------------------------------- CNRS Laboratoire de physique de la mati?re condens?e Parc Valrose avenue J Vallot 06108 Nice, France +33 (0) 492 076 775 http://sites.unice.fr/site/aseminara/ -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Jan 16 09:31:35 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 16 Jan 2017 16:31:35 +0100 Subject: [Nek5000-users] double mesh ? In-Reply-To: References: Message-ID: Dear Paul thank you so much for your help ! I am looking forward to see your example, it will be of great help ! could you tell me something more about the recycling BC as well ? Thank you again for your help Agnese > > From: nek5000-users at lists.mcs.anl.gov > Subject: Re: [Nek5000-users] double mesh ? > Date: January 16, 2017 at 4:05:21 PM GMT+1 > To: "nek5000-users at lists.mcs.anl.gov " > > Reply-To: nek5000-users at lists.mcs.anl.gov > > Dear Agnese, > > Sorry for the long delay in answering this question. The issue is that in fact it is possible > (and used to be a supported but rarely used feature). The challenge, however, is how to > make it relatively easy in the context of the current setup for a variety of applications. > > Nek supports a different gather-scatter handle for each field (velocity, temperature, passive > scalar, etc.) > > We have a handle that allows the user access to the connectivity list, which would give > you an opportunity to "disconnect" the periodic faces for the temperature array. > > I'll try to code up an example for the case of channel flow with a mesh generated by > genbox and post it shortly. > > Incidentally, an alternative approach that we often use for more complex domains > would be to use a recycling boundary condition. I think it would not be as efficient, > however, as periodicity for the case that you propose. > > Once again, apologies for the delayed response. > > 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 22, 2016 4:38 AM > To: nek5000-users at lists.mcs.anl.gov > Subject: [Nek5000-users] double mesh ? > > Hi > > i am quite new, so i?m sorry if i?m getting things wrong. > > I am trying to simulate a turbulent channel flow, with a passive scalar. > in the streamwise direction, i would like to impose periodic boundary conditions for the flow, and non-periodic for the scalar. > in fact the scalar has a point source near the inlet, and i want to see how that develops downstream. > > as far as i understand this is not supported by nek5000, because the last mesh point is identified with the first mesh point when periodic BC are imposed. > correct ? > > my question is : could i define two meshes, one for the scalar and one for the velocity ? > > thanks a lot for your help !!! > agnese > > > Agnese Seminara > -------------------------------- > CNRS > Laboratoire de physique de la mati?re condens?e > Parc Valrose > avenue J Vallot > 06108 Nice, France > +33 (0) 492 076 775 > http://sites.unice.fr/site/aseminara/ > > > _______________________________________________ > 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 Jan 18 08:07:05 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 18 Jan 2017 15:07:05 +0100 Subject: [Nek5000-users] param(2) Message-ID: Hi all, There seems to be some inconsistencies regarding dynamic (\mu)/kinematic (\nu) viscosity in nek: - The manual says that P002 is kinematic viscosity (\nu) - param(2) is stored in VDIFF, which is split by nu_star in split_vis - In the turbChannel example, param(2) is added to the Dynamic Smagorinsky term, which has dimension m^2/s, as has \nu On the other hand: - MAKEVIS is stated to calculate 2*DEL*[mue*(S ....)], using VDIFF_E - The manual says about P008: mavier5.f: param(8) = param(2) ! conduct = dyn.visc. - \mu is listed as a parameter in Chapter 3.3.1 of the manual, regarding the .rea file. In all the examples, param(1) = 1, so \nu = \mu, but for some physical cases it is convenient to keep the dimensional variables. So what is the correct use of param(2)? Thanks, Carl Erik, FFI, Norway From nek5000-users at lists.mcs.anl.gov Wed Jan 18 09:34:15 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 18 Jan 2017 15:34:15 +0000 Subject: [Nek5000-users] param(2) In-Reply-To: References: Message-ID: Thanks Carl Erik. It should be p2 = \mu. (I hadn't noticed that error in the manual.) Examples are stand-alone entities and thus as long as they are self-consistent, that suffices. Thus, in turbChannel (originally set up by me), I know that rho==1, and therefore ignore it always... such that \mu=p2==\nu, and I never bother to divide by rho. So, I think the answer to your question, if you want physical variables is: p1 = rho p2 = mu p7 = rhoCp p8 = conduct. You should also then set p21=0 p22=0 and p24=.01 p25=.01 These latter changes will turn off the non dimensional tolerances and tell the iterative solvers to set the tolerances in a relative sense according to the anticipated scale of your variables. 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: Wednesday, January 18, 2017 8:07 AM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] param(2) Hi all, There seems to be some inconsistencies regarding dynamic (\mu)/kinematic (\nu) viscosity in nek: - The manual says that P002 is kinematic viscosity (\nu) - param(2) is stored in VDIFF, which is split by nu_star in split_vis - In the turbChannel example, param(2) is added to the Dynamic Smagorinsky term, which has dimension m^2/s, as has \nu On the other hand: - MAKEVIS is stated to calculate 2*DEL*[mue*(S ....)], using VDIFF_E - The manual says about P008: mavier5.f: param(8) = param(2) ! conduct = dyn.visc. - \mu is listed as a parameter in Chapter 3.3.1 of the manual, regarding the .rea file. In all the examples, param(1) = 1, so \nu = \mu, but for some physical cases it is convenient to keep the dimensional variables. So what is the correct use of param(2)? Thanks, Carl Erik, FFI, Norway _______________________________________________ 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 Jan 19 06:45:12 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 19 Jan 2017 13:45:12 +0100 Subject: [Nek5000-users] param(2) In-Reply-To: References: Message-ID: Thank you Paul! p2 = \mu is indeed what we have used so far, and it is consistent with changes in p1 and p2 (same \nu = p2/p1 gives same solution). The only "problem" then is the LES model from the turbChannel example, where UDIFF = EDIFF in uservp, and EDIFF = param(2) + dynamic Smagorinsky term in eddy_visc, Here, param(2) is \mu and the dyn.Smag term has dimension as \nu. If we multiply the dyn.Smag.-term by param(1) to obtain a "\mu_T", the simulation crashes after a couple of time-steps (with \rho=1000, for water). But on closer inspection, for PN-PN, the dyn.Smag. term is split off as VDIFF_E (if nu_star = param(2)) and the correponding stress is added to BFX/BFY/BFZ in makevis, and then BFX (etc) is multiplied by \rho in makeabf, so it gets the correct dimension. For PN-PN-2, it's not clear to me how this is handled, but the effect seems to be the same. EDIFF should be set to param(2) + \nu_T in eddy_visc (as it is turbChannel) to get the correct result for different \rho. Carl Erik On 01/18/2017 04:34 PM, nek5000-users at lists.mcs.anl.gov wrote: > > Thanks Carl Erik. > > It should be p2 = \mu. (I hadn't noticed that error in the manual.) > > Examples are stand-alone entities and thus as long as they are self-consistent, > that suffices. Thus, in turbChannel (originally set up by me), I know that rho==1, > and therefore ignore it always... such that \mu=p2==\nu, and I never bother to > divide by rho. > > So, I think the answer to your question, if you want physical variables is: > > p1 = rho > p2 = mu > p7 = rhoCp > p8 = conduct. > > You should also then set > > p21=0 > p22=0 > > and > > p24=.01 > p25=.01 > > These latter changes will turn off the non dimensional tolerances and tell > the iterative solvers to set the tolerances in a relative sense according to the > anticipated scale of your variables. > > 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: Wednesday, January 18, 2017 8:07 AM > To: nek5000-users at lists.mcs.anl.gov > Subject: [Nek5000-users] param(2) > > Hi all, > > There seems to be some inconsistencies regarding dynamic (\mu)/kinematic > (\nu) viscosity in nek: > > - The manual says that P002 is kinematic viscosity (\nu) > - param(2) is stored in VDIFF, which is split by nu_star in split_vis > - In the turbChannel example, param(2) is added to the Dynamic > Smagorinsky term, which has dimension m^2/s, as has \nu > > On the other hand: > - MAKEVIS is stated to calculate 2*DEL*[mue*(S ....)], using VDIFF_E > - The manual says about P008: mavier5.f: param(8) = param(2) ! conduct = > dyn.visc. > - \mu is listed as a parameter in Chapter 3.3.1 of the manual, regarding > the .rea file. > > In all the examples, param(1) = 1, so \nu = \mu, but for some physical > cases it is convenient to keep the dimensional variables. > > So what is the correct use of param(2)? > > Thanks, > Carl Erik, > FFI, Norway > > > _______________________________________________ > 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 Jan 19 10:58:44 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 19 Jan 2017 16:58:44 +0000 Subject: [Nek5000-users] mesh smoothing Message-ID: Greetings Nekians, Rumor on the street is there are some mesh smoothing algorithms available. We've been experiencing some problems with high pressure iteration counts and one of our current suspects is a poor quality mesh. Any help would be appreciated! Thank you, Dillon -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Thu Jan 19 22:49:00 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 20 Jan 2017 04:49:00 +0000 Subject: [Nek5000-users] mesh smoothing Message-ID: Hi Dillon, I have a mesh smoother. Can you email me on kmittal2 at illinois.edu and we can discuss the smoother over there. Ketan From nek5000-users at lists.mcs.anl.gov Sat Jan 21 06:18:41 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Sat, 21 Jan 2017 12:18:41 +0000 Subject: [Nek5000-users] double mesh ? In-Reply-To: References: , Message-ID: Dear Agnese, Attached is an example that shows how to modify the connectivity for the thermal field to be inflow/outflow. If your domain is generated with genbox and your flow is in the x direction then this should also work for 3D. Note that no changes are required to the .rea file (you should have "P " for both velocity and temperature in the .rea/.box file). The changes occur in the .usr file. 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 16, 2017 9:31 AM To: nek5000-users at lists.mcs.anl.gov Subject: Re: [Nek5000-users] double mesh ? Dear Paul thank you so much for your help ! I am looking forward to see your example, it will be of great help ! could you tell me something more about the recycling BC as well ? Thank you again for your help Agnese From: nek5000-users at lists.mcs.anl.gov Subject: Re: [Nek5000-users] double mesh ? Date: January 16, 2017 at 4:05:21 PM GMT+1 To: "nek5000-users at lists.mcs.anl.gov" > Reply-To: nek5000-users at lists.mcs.anl.gov Dear Agnese, Sorry for the long delay in answering this question. The issue is that in fact it is possible (and used to be a supported but rarely used feature). The challenge, however, is how to make it relatively easy in the context of the current setup for a variety of applications. Nek supports a different gather-scatter handle for each field (velocity, temperature, passive scalar, etc.) We have a handle that allows the user access to the connectivity list, which would give you an opportunity to "disconnect" the periodic faces for the temperature array. I'll try to code up an example for the case of channel flow with a mesh generated by genbox and post it shortly. Incidentally, an alternative approach that we often use for more complex domains would be to use a recycling boundary condition. I think it would not be as efficient, however, as periodicity for the case that you propose. Once again, apologies for the delayed response. 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 22, 2016 4:38 AM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] double mesh ? Hi i am quite new, so i?m sorry if i?m getting things wrong. I am trying to simulate a turbulent channel flow, with a passive scalar. in the streamwise direction, i would like to impose periodic boundary conditions for the flow, and non-periodic for the scalar. in fact the scalar has a point source near the inlet, and i want to see how that develops downstream. as far as i understand this is not supported by nek5000, because the last mesh point is identified with the first mesh point when periodic BC are imposed. correct ? my question is : could i define two meshes, one for the scalar and one for the velocity ? thanks a lot for your help !!! agnese Agnese Seminara -------------------------------- CNRS Laboratoire de physique de la mati?re condens?e Parc Valrose avenue J Vallot 06108 Nice, France +33 (0) 492 076 775 http://sites.unice.fr/site/aseminara/ _______________________________________________ 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: -------------- next part -------------- A non-text attachment was scrubbed... Name: t.tgz Type: application/x-gzip Size: 5006 bytes Desc: t.tgz URL: From nek5000-users at lists.mcs.anl.gov Tue Jan 24 03:36:29 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 24 Jan 2017 10:36:29 +0100 Subject: [Nek5000-users] Convergence problem in 2D square cylinder Message-ID: Thank you Tanmoy and Paul. I have redone the mesh. I do not know why, but less layers in the boundary helped to stabilize the convergence. I have also played around with parameters 101 and 103. However, I should check wether the dynamics have changed or not. VBeltranM -- V?ctor Beltr?n Mart?nez Ph.D. Student in Applied Mathematics Escuela T?cnica Superior de Ingenier?a Aeron?utica y del Espacio Universidad Polit?cnica de Madrid From nek5000-users at lists.mcs.anl.gov Tue Jan 24 08:27:57 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 24 Jan 2017 15:27:57 +0100 Subject: [Nek5000-users] nelgv in uservp of conj_ht Message-ID: Hello All, I have a question in uservp of conjugate heat transfer. I am trying to modify the conjugate heat transfer problem (with additional passive scalars) and use for solving MHD problem with vector potential approach. I tested the conjugate heat transfer case changing the IF condition in uservp in following ways. Case 1: if (ieg .gt. nelgv) then ! properties in the solid udiff = 0.1*param(8) ! conductivity utrans = 1.0 endif Case 2: if (y .ge. 1.5D0) then udiff= 0.1*param(8) elseif (y. le. 0.5D0) then udiff= 0.1*param(8) else udiff = param(8) endif I got a small difference in the temperature profile at the interface between solid and fluid. Could you please let me know, how do the above two cases affect the conditions imposed at the interface (what is the mathematical treatment and how does it change). Thank you in advance. Yours sincerely, Sandeep -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Tue Jan 24 09:33:48 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 24 Jan 2017 15:33:48 +0000 Subject: [Nek5000-users] nelgv in uservp of conj_ht In-Reply-To: References: Message-ID: Hi Sandeep, The temperature should be continuous if everything is working correctly. One way things can go wrong is if the initial temperature is not continuous, because subsequent fields are computed as updates of the original. All corrections are continuous: T^n = T^n-1 + dT dT is continuous. T^n will be continuous only if T^n-1 is continuous. useric() generally ensures continuity of temperature and velocity via the dsavg operation before the start of the calculation. Even if the user prescribes a discontinuous initial condition, useric() invokes dsavg to force continuity. How do you know that the function is discontinuous ? 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, January 24, 2017 8:27 AM To: nek5000-users at lists.mcs.anl.gov Subject: [Nek5000-users] nelgv in uservp of conj_ht Hello All, I have a question in uservp of conjugate heat transfer. I am trying to modify the conjugate heat transfer problem (with additional passive scalars) and use for solving MHD problem with vector potential approach. I tested the conjugate heat transfer case changing the IF condition in uservp in following ways. Case 1: if (ieg .gt. nelgv) then ! properties in the solid udiff = 0.1*param(8) ! conductivity utrans = 1.0 endif Case 2: if (y .ge. 1.5D0) then udiff= 0.1*param(8) elseif (y. le. 0.5D0) then udiff= 0.1*param(8) else udiff = param(8) endif I got a small difference in the temperature profile at the interface between solid and fluid. Could you please let me know, how do the above two cases affect the conditions imposed at the interface (what is the mathematical treatment and how does it change). Thank you in advance. Yours sincerely, Sandeep -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Tue Jan 24 11:41:31 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 24 Jan 2017 18:41:31 +0100 Subject: [Nek5000-users] nelgv in uservp of conj_ht Message-ID: Case (1) is correct way to impose properties for the solid. The interface condition will match the fluxes (in a weak sense) while heat capacity can discontinuous. -----Original message----- > From:nek5000-users at lists.mcs.anl.gov > Sent: Tuesday 24th January 2017 17:28 > To: nek5000-users at lists.mcs.anl.gov > Subject: [Nek5000-users] nelgv in uservp of conj_ht > > Hello All, > I have a question in uservp of conjugate heat transfer. I am trying to modify the conjugate heat transfer problem (with additional passive scalars) and use for solving MHD problem with vector potential approach. > > I tested the conjugate heat transfer case changing the IF condition in uservp in following ways. ? ? ? > > Case 1: > ? ? ? ? if (ieg .gt. nelgv) then ?! properties in the solid > ? ? ? ? ? ? udiff ? = 0.1*param(8) ! conductivity > ? ? ? ? ? ? utrans ?= 1.0 > ? ? ? ? ?endif > > Case 2: > if (y .ge. 1.5D0) then > udiff= 0.1*param(8) > elseif (y. le. 0.5D0) then > udiff= 0.1*param(8) > else > udiff ? = param(8) > endif > > I got a small difference in the temperature profile at the interface between solid and fluid.? > Could you please let me know, how do the above two cases affect the conditions imposed at the interface (what is the mathematical treatment and how does it change). > Thank you in advance. > > Yours sincerely, > Sandeep > _______________________________________________ > 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 Jan 24 10:35:40 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 24 Jan 2017 16:35:40 +0000 Subject: [Nek5000-users] nelgv in uservp of conj_ht In-Reply-To: References: Message-ID: Correct -- Case 1 is the right way.. (Sorry, my previous answer was off the mark as I misinterpreted the question..) The point is that properties do not have to be continuous. So in one element (solid), you have a property all the way up to and including the interface. In another (fluid), you have a different property all the way up to an including the interface. In this way, the integrand in the variational formulation sees piecewise smooth functions and you get an excellent (sometimes exact) approximation to the integral of the step function. 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, January 24, 2017 11:41 AM To: nek5000-users at lists.mcs.anl.gov Subject: Re: [Nek5000-users] nelgv in uservp of conj_ht Case (1) is correct way to impose properties for the solid. The interface condition will match the fluxes (in a weak sense) while heat capacity can discontinuous. -----Original message----- > From:nek5000-users at lists.mcs.anl.gov > Sent: Tuesday 24th January 2017 17:28 > To: nek5000-users at lists.mcs.anl.gov > Subject: [Nek5000-users] nelgv in uservp of conj_ht > > Hello All, > I have a question in uservp of conjugate heat transfer. I am trying to modify the conjugate heat transfer problem (with additional passive scalars) and use for solving MHD problem with vector potential approach. > > I tested the conjugate heat transfer case changing the IF condition in uservp in following ways. > > Case 1: > if (ieg .gt. nelgv) then ! properties in the solid > udiff = 0.1*param(8) ! conductivity > utrans = 1.0 > endif > > Case 2: > if (y .ge. 1.5D0) then > udiff= 0.1*param(8) > elseif (y. le. 0.5D0) then > udiff= 0.1*param(8) > else > udiff = param(8) > endif > > I got a small difference in the temperature profile at the interface between solid and fluid. > Could you please let me know, how do the above two cases affect the conditions imposed at the interface (what is the mathematical treatment and how does it change). > Thank you in advance. > > Yours sincerely, > Sandeep > _______________________________________________ > 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 Jan 26 02:39:59 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 26 Jan 2017 09:39:59 +0100 Subject: [Nek5000-users] Nek compiles on IBM, but not on Intel Message-ID: Hi all, I am using Nek5000 for simulations on a Intel Xeon Westmere Cluster and on an IBM BlueGene/Q system. Both systesm have worked fine so far, but as of this morning, Nek suddenly stopped to compile on the Intel system. On the IBM system, I use the mpixlf77, mpixlc compilers and both my own cases and the examples compile. On the Intel system, I use the mpiifort, mpiicc (version 16.0.1) compilers and compiling stops with the message drive1.f: fatal: can't fgets from file: drive1.f drive2.f: fatal: can't fgets from file: drive2.f Both systems are using the current version from github and exactly the same .usr files. Has there been any update to Nek or the intel compilers that might cause this problem? I just tried the 16.0.3 version (the newest installed from Intel) and the error persisted. Thanks in advance for any thougts on this, Jan From nek5000-users at lists.mcs.anl.gov Thu Jan 26 05:02:12 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 26 Jan 2017 12:02:12 +0100 Subject: [Nek5000-users] Nek compiles on IBM, but not on Intel In-Reply-To: References: Message-ID: Ok, I found the error. For some reason, git pull created some empty files (drive?.f and others). A new git clone solved this. A somewhat related question: Is it correct that the interpolation subroutines (intpts_ etc) have been removed from the latest release? I haven?t used them for some time but still had them in the usr file and noticed an error that they are missing. Jan > Am 26.01.2017 um 09:39 schrieb nek5000-users at lists.mcs.anl.gov: > > Hi all, > I am using Nek5000 for simulations on a Intel Xeon Westmere Cluster and on an IBM BlueGene/Q system. Both systesm have worked fine so far, but as of this morning, Nek suddenly stopped to compile on the Intel system. > On the IBM system, I use the mpixlf77, mpixlc compilers and both my own cases and the examples compile. > On the Intel system, I use the mpiifort, mpiicc (version 16.0.1) compilers and compiling stops with the message > drive1.f: fatal: can't fgets from file: drive1.f > drive2.f: fatal: can't fgets from file: drive2.f > Both systems are using the current version from github and exactly the same .usr files. > Has there been any update to Nek or the intel compilers that might cause this problem? I just tried the 16.0.3 version (the newest installed from Intel) and the error persisted. > Thanks in advance for any thougts on this, > Jan > _______________________________________________ > 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 Jan 26 09:28:02 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 26 Jan 2017 16:28:02 +0100 Subject: [Nek5000-users] Nek compiles on IBM, but not on Intel Message-ID: Sorry for the inconvenience, we try hard not to break user space. The interpolation wrapper has been refactored and is no longer compiled by default. Add the following two pre-processor commands to the top of your .usr file #include "intp.f" #define INTP_MAXPTS 10000 Then, just call intp_setup(), intp_do(), intp_free(). Please note, that the API has changed. Same for g2gi() - grid to grid interpolation - which was replaced by our new generic fld reader gfldr() Cheers, Stefan -----Original message----- > From:nek5000-users at lists.mcs.anl.gov > Sent: Thursday 26th January 2017 14:03 > To: nek5000-users at lists.mcs.anl.gov > Subject: Re: [Nek5000-users] Nek compiles on IBM, but not on Intel > > Ok, I found the error. For some reason, git pull created some empty files (drive?.f and others). A new git clone solved this. > A somewhat related question: Is it correct that the interpolation subroutines (intpts_ etc) have been removed from the latest release? I haven?t used them for some time but still had them in the usr file and noticed an error that they are missing. > > Jan > > Am 26.01.2017 um 09:39 schrieb nek5000-users at lists.mcs.anl.gov: > > > > Hi all, > > I am using Nek5000 for simulations on a Intel Xeon Westmere Cluster and on an IBM BlueGene/Q system. Both systesm have worked fine so far, but as of this morning, Nek suddenly stopped to compile on the Intel system. > > On the IBM system, I use the mpixlf77, mpixlc compilers and both my own cases and the examples compile. > > On the Intel system, I use the mpiifort, mpiicc (version 16.0.1) compilers and compiling stops with the message > > drive1.f: fatal: can't fgets from file: drive1.f > > drive2.f: fatal: can't fgets from file: drive2.f > > Both systems are using the current version from github and exactly the same .usr files. > > Has there been any update to Nek or the intel compilers that might cause this problem? I just tried the 16.0.3 version (the newest installed from Intel) and the error persisted. > > Thanks in advance for any thougts on this, > > Jan > > _______________________________________________ > > 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 Mon Jan 30 07:04:31 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 30 Jan 2017 14:04:31 +0100 Subject: [Nek5000-users] Setics subroutine in load_fld Message-ID: Dear Neks, I am relatively new to NEK5000. I am going through some the .usr files that have been used previously in my research group. I notice that in the /userchk/ subroutine, /load_fld /has been called a few times. Upon checking the /load_fld/ function in in the source code, I found that /setics /subroutine is called at the end of /load_fld/ function. Does this mean that each time /load_fld /function is called, the /useric/ subroutine in the .usr file is executed? I would like to mention that no restart file is given in the .rea file. Many thanks. Best regards, Navrose. -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Jan 30 13:15:00 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 30 Jan 2017 12:15:00 -0700 Subject: [Nek5000-users] Setics subroutine in load_fld In-Reply-To: References: Message-ID: Hi Navrose, You are right. Two things are important in the load_fld() subroutine. The include file 'RESTART' and the routine setics. So everytime a load_fld subroutine is being called, the variable values read from the current file would be stored in vx(i,j,k,e), vy(i,j,k,e), vz(i,j,k,e), pr, t and so on... Best Regards, Tanmoy On Mon, Jan 30, 2017 at 6:04 AM, wrote: > Dear Neks, > > I am relatively new to NEK5000. I am going through some > the .usr files that have been used previously in my research > group. I notice that in the *userchk* subroutine, *load_fld *has been > called a few times. Upon checking the *load_fld* function in > in the source code, I found that *setics *subroutine is called at the > end of *load_fld* function. > > Does this mean that each time *load_fld *function is called, the *useric* > subroutine in the .usr file is executed? I would like to mention that > no restart file is given in the .rea file. > > Many thanks. > > Best regards, > > Navrose. > > > > > _______________________________________________ > 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 Jan 31 01:25:49 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 31 Jan 2017 08:25:49 +0100 Subject: [Nek5000-users] Setics subroutine in load_fld In-Reply-To: References: Message-ID: Many thanks Tanmoy. So each time load_fld is called it will load the base flow variables (I am carrying out computations for incompressible N-S equations). What if I want to load the perturbations also from a different file? One way that I see it is to use load_fld for the restart file for perturbations, copy it in the perturbation variables (vxp, vyp, vzp, prp), and then again use load_fld for the base flow file. Is there a direct or better method for this? Thank you. Best wishes, Navrose. On Monday 30 January 2017 08:15 PM, nek5000-users at lists.mcs.anl.gov wrote: > Hi Navrose, > > You are right. Two things are important in the load_fld() subroutine. > The include file 'RESTART' and the routine setics. So everytime a > load_fld subroutine is being called, the variable values read from the > current file would be stored in vx(i,j,k,e), vy(i,j,k,e), vz(i,j,k,e), > pr, t and so on... > > Best Regards, > Tanmoy > > On Mon, Jan 30, 2017 at 6:04 AM, > wrote: > > Dear Neks, > > I am relatively new to NEK5000. I am going through some > the .usr files that have been used previously in my research > group. I notice that in the /userchk/ subroutine, /load_fld /has been > called a few times. Upon checking the /load_fld/ function in > in the source code, I found that /setics /subroutine is called at the > end of /load_fld/ function. > > Does this mean that each time /load_fld /function is called, the > /useric/ > subroutine in the .usr file is executed? I would like to mention that > no restart file is given in the .rea file. > > Many thanks. > > Best regards, > > Navrose. > > > > > _______________________________________________ > 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 Tue Jan 31 11:11:54 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 31 Jan 2017 17:11:54 +0000 Subject: [Nek5000-users] Multi-fluids simulation Message-ID: Dear Nek users, Can Nek5000 simulate multi fluids? e.g. air blowing over water. It would be better if an example of the codes could be provided. Thank you, Peng From nek5000-users at lists.mcs.anl.gov Tue Jan 31 11:49:29 2017 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 31 Jan 2017 10:49:29 -0700 Subject: [Nek5000-users] Setics subroutine in load_fld In-Reply-To: References: Message-ID: Hi Navrose, I do not know if there is any better method, because I would do it probably the same way. Another way-- you can probably define another set of routine say load_fld_pert(), with modified setics and RESTART file, but probably it is as good as using opcopy(). On Tue, Jan 31, 2017 at 12:25 AM, wrote: > Many thanks Tanmoy. > > So each time load_fld is called it will load the base > flow variables (I am carrying out computations for incompressible N-S > equations). What if I want to load the perturbations also from a different > file? One way that I see it is to use load_fld for the restart file for > perturbations, copy it in the perturbation variables (vxp, vyp, vzp, prp), > and then again use load_fld for the base flow file. Is there a direct or > better method for this? > > Thank you. > > Best wishes, > Navrose. > On Monday 30 January 2017 08:15 PM, nek5000-users at lists.mcs.anl.gov wrote: > > Hi Navrose, > > You are right. Two things are important in the load_fld() subroutine. The > include file 'RESTART' and the routine setics. So everytime a load_fld > subroutine is being called, the variable values read from the current file > would be stored in vx(i,j,k,e), vy(i,j,k,e), vz(i,j,k,e), pr, t and so on... > > Best Regards, > Tanmoy > > On Mon, Jan 30, 2017 at 6:04 AM, wrote: > >> Dear Neks, >> >> I am relatively new to NEK5000. I am going through some >> the .usr files that have been used previously in my research >> group. I notice that in the *userchk* subroutine, *load_fld *has been >> called a few times. Upon checking the *load_fld* function in >> in the source code, I found that *setics *subroutine is called at the >> end of *load_fld* function. >> >> Does this mean that each time *load_fld *function is called, the *useric* >> subroutine in the .usr file is executed? I would like to mention that >> no restart file is given in the .rea file. >> >> Many thanks. >> >> Best regards, >> >> Navrose. >> >> >> >> >> _______________________________________________ >> Nek5000-users mailing list >> Nek5000-users at lists.mcs.anl.gov >> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users >> >> > > > _______________________________________________ > Nek5000-users mailing listNek5000-users at lists.mcs.anl.govhttps://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: