We performed hierarchical multi-scale simulations to study the adhesion properties of various epoxy–aluminium interfaces in the absence and presence of water. The epoxies studied differ from each other in their hexagonal ring structures where one contains aromatic and the other aliphatic rings. As aluminium is unavoidably covered with alumina, a cross-linked epoxy structure near an alumina substrate is created and relaxed by performing coarse-grained simulations. To that purpose, we employ a recently developed parameterization method for variable bead sizes. For polymer–metal interactions, a multi-scale parameterization scheme is applied where the relative adsorption of each bead type is quantified. At the mesoscopic scale, the adhesion pro...