Accelerated urbanization increases both the frequency and intensity of heatwaves (HW) and urban heat islands (UHIs). An extreme HW event occurred in 2012 summer that caused temperatures of more than 40°C in Chicago, Illinois, USA, which is a highly urbanized city impacted by UHIs. In this study, multiple numerical models, including the High Resolution Land Data Assimilation System (HRLDAS) and Weather Research and Forecasting (WRF) model, were used to simulate the HW and UHI, and their performance was evaluated. In addition, sensitivity testing of three different WRF configurations was done to determine the impact of increasing model complexity in simulating urban meteorology. Model performances were evaluated based on the statistical performance metrics, the application of a multi-layer urban canopy model (MLUCM) helps WRF to provide the best performance in this study. HW caused rural temperatures to increase by ∼4°C, whereas urban Chicago had lower magnitude increases from the HW (∼2-3°C increases). Nighttime UHI intensity (UHII) ranged from 1.44 to 2.83°C during the study period. Spatiotemporal temperature fields were used to estimate the potential heat-related exposure and to quantify the excessive heat factor (EHF). The EHF during the HW episode provides a risk map indicating that while urban Chicago had higher heat-related stress during this event, the rural area also had high risk, especially during nighttime in central Illinois. This study provides a reliable method to estimate spatiotemporal exposures for future studies of heat-related health impacts.