Measurement in anisotropic heat conduction is difficult due to the presence
of cross-derivative in the governing equation. There are many real world application
involving anisotropic heat conduction, one of them being heat flow
control or manipulation. In this work, numerical method to solve anisotropic
heat conduction phenomena is proposed. The proposed method is also able to
handle nonhomogeneous domain. It is based on meshless Least Square Moving
Particle Semi-Implicit (LSMPS) spatial derivatives discretization. The developed
solver is tested using several benchmark cases, before used to simulate
several designs of thermal cloak and concentrator for heat flow manipulation.
The results of the benchmark cases and heat flow manipulation simulations
show good agreement with the past results from previous works. Therefore,
the proposed method is deemed able to simulate anisotropic heat conduction
in nonhomogeneous domain.