In this paper, a mesoscale model of concrete is presented, which considers particles, matrix material and the interfacial transition zone (ITZ) as separate constituents. Particles are represented as ellipsoides, generated according to a prescribed grading curve and placed randomly into the specimen. Algorithms are proposed to generate realistic particle configurations efficiently. The nonlinear behavior is simulated with a cohesive interface model for the ITZ. For the matrix material, different damage/plasticity models are investigated. The simulation of localization requires to regularize the solution, which is performed by using integral type nonlocal models with strain or displacement averaging. Due to the complexity of a mesoscale model for a realistic structure, a multiscale method to couple the homogeneous macroscale with the heterogeneous mesoscale model in a concurrent embedded approach is proposed. This allows an adaptive transition from a full macroscale model to a multiscale model, where only the relevant parts are resolved on a finer scale. Special emphasis is placed on the investigation of different coupling schemes between the different scales, such as the mortar method and the arlequin method, and a discussion of their advantages and disadvantages within the current context. The applicability of the proposed methodology is illustrated for a variety of examples in tension and compression.