Ultrasonic investigation of damage detection has been widely used for non-destructive testing of various concrete structures. This study focuses on damage detection analysis with the aid of wave propagation in two-phase composite concrete with aggregate (inclusion) and mortar (matrix). To fabricate a realistic simulation model containing a variety of irregular aggregate shapes, the mesh generation technique using an image processing technique was proposed. Initially, the domains and boundaries of the aggregates were extracted from the digital image of a typical concrete cut-section. This enables two different domains: aggregates and mortar in heterogeneous concrete sections, and applied the grids onto these domains to discretize the model. Subsequently, finite element meshes are generated in terms of spatial and temporal requirements of the model size. For improved analysis results, all meshes are designed to be quadrilateral type, and an additional process is conducted to improve the mesh quality. With this simulation model, wave propagation analyses were conducted with a central frequency of 75 kHz of the Mexican hat incident wave. Several void damages, such as needle-shaped cracks and void-shaped holes, were artificially introduced in the model. Finally, various formats of internal damage were detected by implementing energy mapping based signal processing.