def calculate_q_debye_linear(eps_fluid, lambda_d, zeta):
    """
    Calculate the charge accumulated in the Debye layer
    (Adjari, 2006)
    units: Coulombs
    Notes:

    """
    q = -eps_fluid * zeta / lambda_d
    return q