openLaval.py 文件源码

def get_Q(self):
        """ get vortex flow parameter Q
        See also:
            NASA TN D-4421 eq. 34a
        """
        f = lambda Mstar: ((self.gamma + 1) / 2 - (self.gamma - 1) / 2 * Mstar**2)**(1 / (self.gamma - 1))
        a = integrate.quad(f, self.Mlstar, self.Mustar)
        Q = (self.Mlstar * self.Mustar) / (self.Mustar - self.Mlstar) * a[0]
        return Q