Parallel and distributed computing systems are foundational to the success of cloud computing and big data analytics. These systems process computational workflows in a way that can be mathematically modeled by a fork-and-join queueing network with blocking (FJQN/B). While engineering solutions have long been made to build and scale such systems, it is challenging to rigorously characterize their throughput performance at scale theoretically. What further complicates the study is the presence of heavy-tailed delays that have been widely documented therein. To this end, we introduce two fundamental concepts for networks of arbitrary topology (scaling dimension and extended metric dimension) and utilize an infinite sequence of growing FJQN/Bs to study the throughput limit. The throughput is said to be scalable if the throughput limit infimum of the sequence is strictly positive as the network size grows to infinity. We investigate throughput scalability by focusing on heavy-tailed service times that are regularly varying (with index $α>1$) and featuring the network topology described by the two aforementioned dimensions. In particular, we show that an infinite sequence of FJQN/Bs is throughput scalable if the extended metric dimension $<α-1$ and only if the scaling dimension $\leα-1$. These theoretical results provide new insights on the scalability of a rich class of FJQN/Bs with various structures, including tandem, lattice, hexagon, pyramid, tree, and fractals.