def thrust_bearing(**kwargs):
    """
    Analysis for the journal bearing
    
    We only take keyword arguments to avoid problems with order of parameters
    Parameters must be specified by name!
    
    Parameters
    ----------
    mu : float
        Friction coefficient [-]
    V : float
        Contact velocity [m/s]        
    N : float
        Thrust force [N]
    
    Notes
    -----
    To derive the speed of contact of the orbiting scroll
    
    .. math::
    
        x = {r_o}\\cos \\left( {{\\phi _{ie}} - \\frac{\\pi }{2} - \\theta } \\right)

    .. math::
    
        y = {r_o}\\sin \\left( {{\\phi _{ie}} - \\frac{\\pi }{2} - \\theta } \\right)
    
    .. math::
    
        \\frac{{dx}}{{d\\theta }} = {r_o}\\sin \\left( {{\\phi _{ie}} - \\frac{\\pi }{2} - \\theta } \\right)\\omega 
    
    .. math::
    
        \\frac{{dy}}{{d\\theta }} =  - {r_o}\\cos \\left( {{\\phi _{ie}} - \\frac{\\pi }{2} - \\theta } \\right)\\omega

    .. math::
    
        \\left| v \\right| = \\sqrt {{{\\left( {\\frac{{dx}}{{d\\theta }}} \\right)}^2} + {{\\left( {\\frac{{dy}}{{d\\theta }}} \\right)}^2}}  = {r_o}\\omega
        
    But it is quite possible that you do not have hydro-dynamic lubrication, and 
    as a result you can get asperity-asperity contact and much higher friction
    coefficients - see for example Kobayashi et al. "Experimental Study on Journal Bearing
    Characteristics in Reciprocating Compressors for HFC-134a" Purdue Compressor
    Conferences 1998.  http://docs.lib.purdue.edu/cgi/viewcontent.cgi?article=2257&context=icec
    
    Can be off by a factor of as much as 10 times at low
    Sommerfeld number
     
    """
    mu = kwargs.pop('mu', None)
    V = kwargs.pop('V', None)
    N = kwargs.pop('N', None)
    Wdot_loss_thrust = mu * V * N
    return dict(mu=mu, V=V, Wdot_loss=Wdot_loss_thrust)
