Modern audio systems and musical effects feature multi-core processing units. The development of parallel audio processing algorithms capable of exploiting the architecture of such hardware is thus in order. In this paper, we present a parallel version of the Hierarchical Scattering Iterative Method (HSIM), a technique based on Wave Digital Filter principles recently proposed for the emulation of multiphysics audio circuits containing multiple nonlinear oneports and nonlinear transformers. HSIM operates in a modular fashion, and it is characterized by a high number of embarrassingly parallelizable operations, making it a good candidate for parallel execution. After analyzing HSIM from the parallel computing perspective, we propose three different strategies for the distribution of HSIM workload among threads of execution, and we show how to compute the maximum achievable speedup. The emulation of a possible output stage of a vacuum-tube guitar amplifier is considered, and a performance comparison between parallel and serial implementations of HSIM is presented, pointing out a speedup of nearly 30%. The proposed method proves thus to be promising for Virtual Analog modeling applications, leading the way towards the parallel digital emulation of increasingly complex audio circuits.