In this work, we present an efficient numerical tool for the prediction of the final microstructure, mechanical properties, and distortions of automotive steel spindles subjected to quenching processes by immersion in liquid tanks. The complete model, which consists of a two-way coupled thermal–metallurgical model and a subsequent (one-way coupled) mechanical model, was numerically implemented using finite element methods. The thermal model includes a novel generalized solid-to-liquid heat transfer model that depends explicitly on the piece’s characteristic size, the physical properties of the quenching fluid, and quenching process parameters. The resulting numerical tool is experimentally validated by comparison with the final microstructure and hardness distributions obtained on automotive spindles subjected to two different industrial quenching processes: (i) a batch-type quenching process with a soaking air-furnace stage prior to the quenching, and (ii) a direct quenching process where the pieces are submerged directly in the liquid just after forging. The complete model retains accurately, at a reduced computational cost, the main features of the different heat transfer mechanisms, with deviations in the temperature evolution and final microstructure lower than 7.5% and 12%, respectively. In the framework of the increasing relevance of digital twins in industry, this model is a useful tool not only to predict the final properties of quenched industrial pieces but also to redesign and optimize the quenching process.