[Nek5000-users] problem with 2d axisymmetric numerical simulation
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Tue Jan 3 07:12:38 CST 2017
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
More information about the Nek5000-users
mailing list