def dP_ESDU_low_fin(m, A_min, A_increase, flow_area_contraction_ratio,
    tube_diameter, fin_height, bare_length, pitch_parallel, pitch_normal,
    tube_rows, rho, mu):
    """Calculates the air-side pressure drop for a low-finned tube bank
    according to the ESDU [1]_ method, as described in [2]_. This includes the
    effects of friction of the fin, and acceleration.
 
    .. math::
        \\Delta P = (K_{acc} + n_{rows} K_{f}) \\frac{1}{2}\\rho v_{max}^2

    .. math::
        K_{f} = 4.71 Re_D^{-0.286} \\left(\\frac{\\text{fin height}}
        {\\text{bare length}}\\right)^{0.51}
        \\left(\\frac{p_1 - D_o}{p_2 - D_o}\\right)^{0.536}
        \\left(\\frac{D_o}{p_1 - D_o}\\right)^{0.36}
        
    .. math::
        K_{acc} = 1 + \\text{(flow area contraction ratio)}^2
            
    Parameters
    ----------
    m : float
        Mass flow rate of air across the tube bank, [kg/s]
    A_min : float
        Minimum air flow area, [m^2]
    A_increase : float
        Ratio of actual surface area to bare tube surface area
        :math:`A_{increase} = \\frac{A_{tube}}{A_{bare, total/tube}}`, [-]
    flow_area_contraction_ratio : float
        Ratio of `A_min` to `A_face`, [-]
    tube_diameter : float
        Diameter of the bare tube, [m]
    fin_height : float
        Height above bare tube of the tube fins, [m]
    bare_length : float
        Length of bare tube between two fins 
        :math:`\\text{bare length} = \\text{fin interval} - t_{fin}`, [m]
    pitch_parallel : float
        Distance between tube center along a line parallel to the flow;
        has been called `longitudinal` pitch, `pp`, `s2`, `SL`, and `p2`, [m]
    pitch_normal : float
        Distance between tube centers in a line 90° to the line of flow;
        has been called the `transverse` pitch, `pn`, `s1`, `ST`, and `p1`, [m]
    tube_rows : int
        Number of tube rows per bundle, [-]
    rho : float
        Average (bulk) density of air across the tube bank, [kg/m^3]
    mu : float
        Average (bulk) viscosity of air across the tube bank, [Pa*s]
        
    Returns
    -------
    dP : float
        Overall pressure drop across the finned tube bank, [Pa]

    Notes
    -----
    Low fins are fins which were formed on the tube outside wall, normally
    by the cold rolling process. The data used by the ESDU covered:
        * fin density 11 to 32/inch
        * tube outer diameters 0.5 to 1.25 inches
        * fin heights 0.03 to 0.1 inches
        * Reynolds numbers 1000 to 80000
        
    [1]_ compared this correlation with 81 results and obtained a standard
    deviation of 7.7%. 
    
    The Reynolds number used in this equation is that based on `V_max`, 
    calculated using the minimum flow area.
    
    Examples
    --------
    >>> AC = AirCooledExchanger(tube_rows=4, tube_passes=4, tubes_per_row=8, tube_length=0.5, 
    ... tube_diameter=0.0164, fin_thickness=0.001, fin_density=1/0.003,
    ... pitch_normal=0.0313, pitch_parallel=0.0271, fin_height=0.0041, corbels=True)

    >>> dP_ESDU_low_fin(m=0.914, A_min=AC.A_min, A_increase=AC.A_increase, 
    ... flow_area_contraction_ratio=AC.flow_area_contraction_ratio,
    ... tube_diameter=AC.tube_diameter, fin_height=AC.fin_height, 
    ... bare_length=AC.bare_length, pitch_parallel=AC.pitch_parallel, 
    ... pitch_normal=AC.pitch_normal, tube_rows=AC.tube_rows, rho=1.217, 
    ... mu=0.000018)
    464.54331418655914
    
    References
    ----------
    .. [1] "High-Fin Staggered Tube Banks: Heat Transfer and Pressure Drop for 
       Turbulent Single Phase Gas Flow." ESDU (October 1, 1986). 
    .. [2] Hewitt, G. L. Shires T. Reg Bott G. F., George L. Shires, and
       T. R. Bott. Process Heat Transfer. 1E. Boca Raton: CRC Press, 1994.
    """
    Vmax = m / (A_min * rho)
    Re = Vmax * tube_diameter * rho / mu
    Kf = 4.72 * Re ** -0.286 * (fin_height / bare_length) ** 0.51 * ((
        pitch_normal - tube_diameter) / (pitch_parallel - tube_diameter)
        ) ** 0.536 * (tube_diameter / (pitch_normal - tube_diameter)) ** 0.36
    Ka = 1.0 + flow_area_contraction_ratio * flow_area_contraction_ratio
    dP = (Ka + tube_rows * Kf) * 0.5 * rho * Vmax * Vmax
    return dP
