This study is concerned with the computation of turbulent flow and heat transfer in U‐bends of strong curvature. Following the earlier studies within the authors' group on flows through round‐ended U‐bends, here attention is turned to flows through square‐ended U‐bends. Flows at two Reynolds numbers have been computed, one at 100,000 and the other at 36,000. In the heat transfer analysis, the Prandtl number was either 0.72 (air) or, in a further departure from our earlier studies, 5.9 (water). The turbulence modelling approaches examined, include a two‐layer and a low‐Re k‐ε model, a two‐layer and a low‐Re version of the basic differential stress model (DSM) and a more recently developed, realisable version of the differential stress model that is free of wall‐parameters. For the low‐Re effective viscosity model (EVM) and DSMs, an alternative, recently proposed length‐scale correction term, independent of wall distance has also been tested. Even the simplest model employed – two‐layer EVM – reproduces the mean flow development with reasonable accuracy, suggesting that the mean flow development is mainly influenced by mean pressure rather than the turbulence field. The heat transfer parameters, on the other hand, show that only the low‐Re DSMs produce reliable Nusselt number predictions for both Prandtl numbers examined.