helper_functions.py 文件源码

def momentum_thickness(y, v):
    """Compute the momentum thickness using Simpson's method.

    Parameters
    ----------
    y : ndarray
        The independent variable.
    v : ndarray
        The velocity values.

    Returns
    -------
    float
        The value of the momentum thickness.

"""
    return simps(v/v[-1]*(1-v/v[-1]), x=y)