Inverse patchy particles are promising colloids to develop new architectures in ceramic materials based on their self-assembly. Nonetheless, a good understanding of their aggregation is required. Several previous studies have shown that the behavior of ceramic colloids can be well described by the DLVO interaction potential. In the present paper, we develop new coarse-grained Brownian dynamics simulations, where particles are represented by an assembly of beads interacting via DLVO interactions, whose parameters can be directly linked to experimental characterization. First, the validity of the simulations is proved by studying the heteroaggregation of homogeneously charged particles. Then, simulations are applied to one-patch inverse patchy particles to study the effect of the patch size. They show that the smaller the patch, the more elongated the aggregates. Simulations are also performed to understand the role of the Debye screening length in the particular case of large patches and they show that aggregation leads always to compact aggregates.