Molecular dynamics simulation techniques are used to study the process of aggregation of highly charged stiff polyelectrolytes due to the presence of multivalent salt. The dominant kinetic mode of aggregation is found to be the case of one end of one polyelectrolyte meeting others at right angles, and the kinetic pathway to bundle formation is found to be similar to that of flocculation dynamics of colloids as described by Smoluchowski. The aggregation process is found to favor the formation of finite bundles of 10-11 filaments at long times. Comparing the distribution of the cluster sizes with the Smoluchowski formula suggests that the energy barrier for the aggregation process is negligible. Also, the formation of long-lived metastable structures with similarities to the raft-like structures of actin filaments is observed within a range of salt concentration.