Heavy QCD axion model in light of pulsar timing arrays