rsfbpic/rswake/lbn_wake.py
def calc_Ez_on_axis_no_beam(n_pe, rb, rb_max):
    """
    Calculate the longitudinal electric field in a plasma bubble.

    Valid in "strong bubble regime", where rb_max*k_pe >> 1.
    Valid only in back of bubble, where there is no drive beam.
    The calculation is local; it depends on bubble radius and slope.
    Args:
        n_pe:    number density of the electron plasma
        rb:      the local bubble radius
        rb_max: maximum value of the bubble radius
    Returns:
        Ez: longitudinal electric field along the axis
    """
    drb_dxi = calc_drb_dxi_no_beam(rb, rb_max)

    # from Eq. (3) of LBN2017
    # there is ambiguity in the sign, which needs to be resolved
    Ez = -2.*math.pi*n_pe*np.abs(scipy.constants.e*rsmath.const.MKS_factor)*rb*drb_dxi
    return Ez

    Calculate the longitudinal electric field in a plasma bubble.

    Valid in "strong bubble regime", where rb_max*k_pe >> 1.
    Valid only in back of bubble, where there is no drive beam.
    The calculation is local; it depends on bubble radius and slope.
    Args:
        n_pe:    number density of the electron plasma
        rb:      the local bubble radius
        rb_max: maximum value of the bubble radius
    Returns:
        Ez: longitudinal electric field along the axis
    
Calculate the longitudinal electric field in a plasma bubble.
Valid in "strong bubble regime", where rb_max*k_pe >> 1.
Valid only in back of bubble, where there is no drive beam.
The calculation is local; it depends on bubble radius and slope.
Calculate the longitudinal electric field in a plasma bubble . Valid in " strong bubble regime " where rb_max * k_pe >> 1 . Valid only in back of bubble where there is no drive beam . The calculation is local ; it depends on bubble radius and slope .
def calc_Ez_on_axis_no_beam(n_pe, rb, rb_max):
    drb_dxi = calc_drb_dxi_no_beam(rb, rb_max)
    Ez = -2.*math.pi*n_pe*np.abs(scipy.constants.e*rsmath.const.MKS_factor)*rb*drb_dxi
    return Ez
def calc_Ez_on_axis_no_beam ( n_pe , rb , rb_max ) : drb_dxi = calc_drb_dxi_no_beam ( rb , rb_max ) Ez = - 2. * math . pi * n_pe * np . abs ( scipy . constants . e * rsmath . const . MKS_factor ) * rb * drb_dxi return Ez
Calculate the longitudinal electric field in a plasma bubble.
Calculate the longitudinal electric field in a plasma bubble .
"""
    Calculate the longitudinal electric field in a plasma bubble.

    Valid in "strong bubble regime", where rb_max*k_pe >> 1.
    Valid only in back of bubble, where there is no drive beam.
    The calculation is local; it depends on bubble radius and slope.
    Args:
        n_pe:    number density of the electron plasma
        rb:      the local bubble radius
        rb_max: maximum value of the bubble radius
    Returns:
        Ez: longitudinal electric field along the axis
    """ # from Eq. (3) of LBN2017 # there is ambiguity in the sign, which needs to be resolved
longitudinal electric field along the axis
longitudinal electric field along the axis
Ez
n_pe
number density of the electron plasma
number density of the electron plasma
rb
the local bubble radius
the local bubble radius
rb_max
maximum value of the bubble radius
maximum value of the bubble radius