openLaval.py 文件源码

def get_C(self):
        """ get reduction in maximum weight flow due to two-dimensional flow
        See also:
            NASA TN D-4421 eq. 34b
        """
        Kstar_max = self.get_Kstar_max()
        a = 1 - np.sqrt((self.gamma + 1)/(self.gamma - 1)) * ((self.gamma + 1) / 2)**(1 / (self.gamma - 1)) * (self.Mustar / (self.Mustar - self.Mlstar))
        f = lambda Mstar: Kstar_max * (1 - (Kstar_max / self.Mlstar)**2 * Mstar**2)**(1 / (self.gamma - 1))
        b = integrate.quad(f, self.Mlstar, self.Mustar)
        C = a * b[0]
        return C