The introduction of eddy diffusivity and subsequently the turbulent Prandtl number works as a way to define a simple relationship between the extra shear stress and heat flux that is present in turbulent flow. If the momentum and thermal eddy diffusivities are zero (no apparent turbulent shear stress and heat flux), then the turbulent flow equations reduce to the laminar equations. We can define the eddy diffusivities for momentum transfer ${\displaystyle \varepsilon _{M}}$  and heat transfer ${\displaystyle \varepsilon _{H}}$  as
${\displaystyle -{\overline {u'v'}}=\varepsilon _{M}{\frac {\partial {\bar {u}}}{\partial y}}}$  and ${\displaystyle -{\overline {v'T'}}=\varepsilon _{H}{\frac {\partial {\bar {T}}}{\partial y}}}$ 
where ${\displaystyle -{\overline {u'v'}}}$  is the apparent turbulent shear stress and ${\displaystyle -{\overline {v'T'}}}$  is the apparent turbulent heat flux.
The turbulent Prandtl number is then defined as
${\displaystyle \mathrm {Pr} _{\mathrm {t} }={\frac {\varepsilon _{M}}{\varepsilon _{H}}}.}$ 

The turbulent Prandtl number has been shown to not generally equal unity (e.g. Malhotra and Kang, 1984; Kays, 1994; McEligot and Taylor, 1996; and Churchill, 2002). It is a strong function of the molecular Prandtl number amongst other parameters and the Reynolds Analogy is not applicable when the molecular Prandtl number differs significantly from unity as determined by Malhotra and Kang;^[1] and elaborated by McEligot and Taylor^[2] and Churchill ^[3]

Turbulent momentum boundary layer equation:
${\displaystyle {\bar {u}}{\frac {\partial {\bar {u}}}{\partial x}}+{\bar {v}}{\frac {\partial {\bar {u}}}{\partial y}}=-{\frac {1}{\rho }}{\frac {d{\bar {P}}}{dx}}+{\frac {\partial }{\partial y}}\left[(\nu {\frac {\partial {\bar {u}}}{\partial y}}-{\overline {u'v'}})\right].}$ 
Turbulent thermal boundary layer equation,
${\displaystyle {\bar {u}}{\frac {\partial {\bar {T}}}{\partial x}}+{\bar {v}}{\frac {\partial {\bar {T}}}{\partial y}}={\frac {\partial }{\partial y}}\left(\alpha {\frac {\partial {\bar {T}}}{\partial y}}-{\overline {v'T'}}\right).}$  Substituting the eddy diffusivities into the momentum and thermal equations yields
${\displaystyle {\bar {u}}{\frac {\partial {\bar {u}}}{\partial x}}+{\bar {v}}{\frac {\partial {\bar {u}}}{\partial y}}=-{\frac {1}{\rho }}{\frac {d{\bar {P}}}{dx}}+{\frac {\partial }{\partial y}}\left[(\nu +\varepsilon _{M}){\frac {\partial {\bar {u}}}{\partial y}}\right]}$ 
and
${\displaystyle {\bar {u}}{\frac {\partial {\bar {T}}}{\partial x}}+{\bar {v}}{\frac {\partial {\bar {T}}}{\partial y}}={\frac {\partial }{\partial y}}\left[(\alpha +\varepsilon _{H}){\frac {\partial {\bar {T}}}{\partial y}}\right].}$ 
Substitute into the thermal equation using the definition of the turbulent Prandtl number to get
${\displaystyle {\bar {u}}{\frac {\partial {\bar {T}}}{\partial x}}+{\bar {v}}{\frac {\partial {\bar {T}}}{\partial y}}={\frac {\partial }{\partial y}}\left[(\alpha +{\frac {\varepsilon _{M}}{\mathrm {Pr} _{\mathrm {t} }}}){\frac {\partial {\bar {T}}}{\partial y}}\right].}$ 

In the special case where the Prandtl number and turbulent Prandtl number both equal unity (as in the Reynolds analogy), the velocity profile and temperature profiles are identical. This greatly simplifies the solution of the heat transfer problem. If the Prandtl number and turbulent Prandtl number are different from unity, then a solution is possible by knowing the turbulent Prandtl number so that one can still solve the momentum and thermal equations.

In a general case of three-dimensional turbulence, the concept of eddy viscosity and eddy diffusivity are not valid. Consequently, the turbulent Prandtl number has no meaning. ^[4]

• Kays, William; Crawford, M.; Weigand, B. (2005). Convective Heat and Mass Transfer, Fourth Edition. McGraw-Hill. ISBN 978-0-07-246876-2.