Boundary Integral Equations - The Boundary Element Method with Programming - page 119

dSx= C/k*(ThetB-ThetA)
Fact= CosthB/CosthA
IF (Fact > 0.0) THEN
dSy= -C/k*LOG(Fact)
ELSE
dSy= 0.
END IF
dRx= -C/H*(costhB*SINThB - cosThA*sinThA)
dRy= C/H*(costhB**2 - cosThA**2)
up= up + dU*t0(ns) - dT*u(ns)
qxp= -k*(dSx*t0(ns)-dRx*u(ns)) ! q-x'
qyp= -k*(dSy*t0(ns)-dRy*u(ns)) ! q-y'
qx= qx + qxp*vn(1,ns) - qyp*vn(2,ns)
qy= qy + qxp*vn(2,ns) + qyp*vn(1,ns)
END DO &
Segments5
Up= Up - q/k*xS(2,np) ! superimpose solutions
qy= qy + q
WRITE (11,'(5(A,F10.3))') &
'x=',xS(1,np),', y=',xS(2,np),', T=',up,', q-x=',qx,', q-
y=',qy
END DO &
Interior_points
STOP
END PROGRAM Direct_Method

INPUT DATA for program Direct_method

1.0 Problem specification
q,k, nseg, rq q Heat inflow

k Thermal conductivity
nseg Number of segments
rq Radius of isolator

4.0 Interior point specification
Npoints Number of interior points
5.0 Interior point coordinates (Npoints cards)
x,y x,y coordinates of interior points

5.6.1 Sample input and output

Here we show the input file for the calculation of the problem in Figure 5.11, with 16
segments and interior points along a horizontal and vertical line and the output file
generated by Program 5.2.