/TITLE, Radial heat conduction in pipe, heatp1.a
  ! 1.9-in pipe insulated with 3-in magnesia
  ! eight elements radially, five axially
  ! Prog Thermo 87, page 7-46 Newnan, EIT
  ! heat /foot of pipe is 5*20*heat = 47.98 BTU/hr-F-foot
  ! axial node numbering
  ! exact answer is 47.78 vs 47.98 below

/PREP7                    ! preprocessor section
ET,1,Plane55,,,1          ! axisymmetric option
MP,KXX,1,0.042            ! thermal conductivity
len = 1/2                 ! pipe length
rad1 = 1.9/2/12           ! pipe radius
rad2 = (1.9/2+3)/12       ! interface radius
 ! elra < 10
elra = 6                  ! radial number of elements
elax = 5                  ! axial number of elements
nax = elax+1              ! axial number of nodes
nrad = elra+1             ! radial number of nodes

N,1,rad1,0                ! lower-left node
N,nax,rad1,len            ! upper-left node
FILL,1,nax                ! axial nodes
NGEN,nrad,10,1,nax,1,((rad2-rad1)/elra),0 ! axial nodes
/PNUM,node,1              ! plot the nodes
NPLOT
E,1,11,12,2        ! lower-left element
EGEN,elax,1,-1     ! axial elements
EGEN,elra,10,-elax ! remaining elements
OUTPR,,1
FINISH             ! preprocessor section

/SOLUTION          ! solution section
antype,static
outpr,basic,1
outpr,nload,1
elfar = (elra*10+1)
D,1,TEMP,338,,nax,1    ! temp of inside nodes
D,elfar,TEMP, 80,,(elfar+elax),1 ! temp of outside nodes
SOLVE
EPLOT
FINISH          ! solution section

/POST1          ! post processor
PRNSOL,TEMP
! /output,Temp,txt,D:\Ansys\   !  Temperature and heat
 /output,Temp,txt,U:\win\Ansys\   !  Temperature and heat

! FINISH          ! post processor

! Number of elem   Heat     Error
!  radial  axial             %
!  exact          47.78
!    4       5    48.5362   1.6
!    6       5    48.1290   0.7
!    8       5    47.9788   0.4

!   NODE  TEMP  NODE  TEMP    NODE  TEMP   NODE  TEMP
!     1  338.00   2  338.00     3  338.00    4  338.00
!     5  338.00   6  338.00    11  262.01   12  262.01
!    13  262.01  14  262.01    15  262.01   16  262.01
!    21  208.37  22  208.37    23  208.37   24  208.37
!    25  208.37  26  208.37    31  166.92   32  166.92
!    33  166.92  34  166.92    35  166.92   36  166.92
!    41  133.14  42  133.14    43  133.14   44  133.14
!    45  133.14  46  133.14    51  104.65   52  104.65
!    53  104.65  54  104.65    55  104.65   56  104.65
!    61  80.000  62  80.000    63  80.000   64  80.000
!    65  80.000  66  80.000

!   NODE      HEAT
!     1     2.40645
!     2     4.81290
!     3     4.81290
!     4     4.81290
!     5     4.81290
!     6     2.40645
!    61    -2.40645
!    62    -4.81290
!    63    -4.81290
!    64    -4.81290
!    65    -4.81290
!    66    -2.40645