This function calculates the Nusselt number for a rectangular duct according to Muzychka and Yovanovich [1].
The dimensionless length zStern as input has to be calculated as follows.
zStern = (sduct ⁄ √Across-section ) ⁄ (Re Pr )
[1]: Muzychka, Y. S.; Yovanovich, M. M. : Laminar Forced Convection Heat Transfer in the Combined Entry Region of Non-Circular Ducts ; Transactions of the ASME; Vol. 126; February 2004
function NusseltNumberMuzychka input Real Re "Reynolds number"; input Real Pr "Prandtl number"; input Real aspRat "aspect ratio between duct height and width"; input Real zStern "dimensionless length"; input Boolean uniWalTem "true if uniform wall temperature boundary conditions"; input Boolean local "true if local nusslet number or false if average shall be calculated"; input Real gamma "shape parameter (rectangular duct: 0.1)"; output Real Nu "Nusselt number"; end NusseltNumberMuzychka;