Multimodal simulation of large area silicon photomultipliers for time resolution optimization