Hybrid Numerical Models For Fast Design Of Terahertz Plasmonic Devices