An unfitted HDG method for a semi-linear elliptic problem arising from nuclear fusion reaction