Exemplo n.º 1
0
    def test_n2ls_1500(self):
        comp = CEA_tp()

        this_nguess = np.array([6.58936767e-06,   2.27158983e-02,   6.24254002e-06,   2.27252139e-02])

        this_ch = np.array([[  2.27224876e-02,   4.54383859e-02,   2.27224876e-02],
                     [  4.54383859e-02,   9.08951525e-02,   4.54508710e-02],
                     [  2.27224876e-02,   4.54508710e-02,   3.51630007e-06]])

        this_rhs = np.array([-1.40217711, -2.80452075, -1.40240467])
        
        this_mu = np.array([-43.73879404, -61.71398144, -35.95037481])

        comp.nguess = this_nguess
        comp.T = 1500
        comp.P = 1.034210
        comp.run()

        rel_error(comp.rhs, this_rhs, 0.001)
        rel_error(comp.ch, this_ch, 0.001)
        rel_error(comp.mu, this_mu, 0.001)
Exemplo n.º 2
0
    def test_n2ls_4000(self):
        comp = CEA_tp()

        this_nguess = np.array([ 0.02040748,  0.00231478,  0.01020431,  0.03292581])

        this_ch = np.array([[  2.27222641e-02,   2.50370446e-02,   2.27222641e-02],
         [  2.50370446e-02,   7.04838323e-02,   4.54456580e-02],
         [  2.27222641e-02,   4.54456580e-02,   7.59390390e-07]])

        this_rhs = np.array([-0.81078449, -1.59262881, -1.14346367])
        
        this_mu = np.array([-34.02175706, -50.32264785, -32.60178159])

        comp.nguess = this_nguess
        comp.T = 4000
        comp.P = 1.034210
        comp.run()

        rel_error(comp.rhs, this_rhs, 0.001)
        rel_error(comp.ch, this_ch, 0.001)
        rel_error(comp.mu, this_mu, 0.001)