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

yp= rp*SIN(Thetp) ! y-coordinate of source point
dxr(1)= xp-xi
dxr(2)= yp-yi
r= SQRT(dxr(1)**2 + dxr(2)**2)
uq= uq + U(r,k,2)*F(npp)
Thetp= Thetp + Delth ! angle to next source point P
END DO &
Source_points2
uq=uq-q/k*yi
WRITE(11,*) 'Temperature at x=',xi,', y=',yi,' =',uq
END DO &
Int_points
STOP
END PROGRAM Trefftz

INPUT DATA for program Trefftz

1.0 Problem specification
q,k, npnts, rq,rp

q … Heat inflow
k … Thermal conductivity
npnts … Number of points P,Q
rq … Radius of isolator
rp … Radius of sources

2.0 Interior point specification
Npoints Number of interior points
3.0 Interior point coordinates (Npoints cards)
x,y x,y coordinates of interior points

5.3.1 Sample input and output

Here we show an example of the input for an isolator of radius 1.0 with 32 points P and
Q, where the source points P are situated along a circle with a radius 0.7.

File INPUT.DAT
1.0 1.0 32 1.0 0.7
18
0. -5.
0. -4.5
0. -4.
0. -3.5
0. -3.
0. -2.5
0. -2.
0. -1.5