The interaction of a crack with perfectly bonded rigid isolated inclusions and clusters of inclusions in a brittle matrix is investigated using numerical simulations. Of particular interest is the role inclusions play on crack paths, stress intensity factors (SIFs) and the energy release rates with potential implications to the fracture behavior of particulate composites. The effects of particle size and eccentricity relative to the initial crack orientation are examined first as a precursor to the study of particle clusters. Simulations are accomplished using a new quasi-static crack-growth prediction tool based on the symmetric-Galerkin boundary element method, a modified quarter-point crack-tip element, the displacement correlation technique for evaluating SIFs, and the maximum principal stress criterion for crack-growth direction prediction. The numerical simulations demonstrate a complex interplay of crack-tip shielding and amplification mechanisms leading to significant toughening of the material.