def reynolds_number(rho, vel, d_m, mu):
    """
    function for determining Reynolds number

    :param rho: density
    :param vel: velocity
    :param d_m: pipe diameter
    :param mu: viscosity
    :return: Reynolds number
    """
    return rho * vel * d_m / mu