The aggregation behavior of Asphaltene molecules in solution was studied by molecular dynamic simulation. Three Asphaltene models were studied in Hexane and o-Dichlorobenzene. The simulations were carried out using the GROMACS software, while the GROMOS 54A7 force field was used for the solvent models. The temperature range for the simulation was 300K to 350K, while the system pressure was at 1bar. The Radius of gyration, Cluster size distribution, and radial distribution function g(r), for the three Asphaltenes in solution, was determined from the study. The molecules form aggregates leading to the formation of parallel and non-parallel π- π stacking between the aromatic rings and resulting in a relatively stable aggregate structure. In Hexane, molecules with larger conjugate cores form aggregates more rapidly at lower temperatures while in o-Dichlorobenzene molecules which has a longer lateral chain formed more stable aggregates at lower temperatures. In comparison, aggregates are denser and more stable in Hexane than in o-Dichlorobenzene, while the latter is a better solvent for Asphaltene. It was also demonstrated that at higher temperatures the Asphaltene molecules will dissociate from the bulk mass to form smaller unstable aggregates, causing a decrease in the radius of gyration (R G ) of Asphaltene aggregates. Furthermore, the results demonstrated that the temperature, Asphaltene structures, and the presence of polar groups play a significant role in Asphaltene-solvent interactions.