We have studied the thermal conductivity of confined superfluids on a barli
ke geometry. We use the planar magnet lattice model on a lattice H x H x L
with L much greater thanH. We have applied open boundary conditions on the
bar sides (the confined directions of length H) and periodic along the long
direction. We have adopted a hybrid Monte Carlo algorithm to efficiently d
eal with the critical slowing down and in order to solve the dynamical equa
tions of motion we use a discretization technique which introduces errors o
nly O[(deltat)(6)] in the time step deltat. Our results demonstrate the con
sistency of scaling using known values of the critical exponents and we obt
ained the scaling function of the thermal resistivity. We find that our res
ults for the thermal resistivity scaling function are in very good agreemen
t with the available experimental results for pores using the temperature s
cale and thermal resistivity scale as free fitting parameters.