Solar radiation impacts diverse aspects of city life, such as harvesting energy with PV panels, passive heating of buildings in winter, cooling the loads of air-conditioning systems in summer, and the urban microclimate. Urban digital twins and 3D city models can support solar studies in the process of urban planning and provide valuable insights for data-driven decision support. This study examines the calculation of solar incident radiation at the city scale in Sofia using remote sensing data for the large shading context in a mountainous region and 3D building data. It aims to explore the methods of geometry optimisation, limitations, and performance issues of a 3D computer-aided design (CAD) tool dedicated to small-scale solar analysis and employed at the city scale. Two cases were considered at the city and district scales, respectively. The total face count of meshes for the simulations constituted approximately 2,000,000 faces. A total of 64,379 roofs for the whole city and 4796 buildings for one district were selected. All calculations were performed in one batch and visualised in a 3D web platform. The use of a 3D CAD environment establishes a seamless process of updating 3D models and simulations, while preprocessing in Geographic Information System (GIS) ensures working with large-scale datasets. The proposed method showed a moderate computation time for both cases and could be extended to include reflected radiation and dense photogrammetric meshes in the future.