#NEWLINE# def test_super_pore_conductance( self ) : #NEWLINE# #INDENT# g_super = [ ] #NEWLINE# BC1_pores = np . arange ( 20 , 30 ) #NEWLINE# self . alg . set_boundary_conditions ( bctype = 'Dirichlet' , bcvalue = 0.4 , pores = BC1_pores ) #NEWLINE# BC2_pores = np . arange ( 45 , 66 ) #NEWLINE# self . alg . set_boundary_conditions ( bctype = 'Neumann_group' , bcvalue = 1.4e-10 , pores = BC2_pores ) #NEWLINE# g_super . append ( 2e-12 ) #NEWLINE# BC3_pores = np . arange ( 87 , 94 ) #NEWLINE# self . alg . set_boundary_conditions ( bctype = 'Neumann_group' , bcvalue = ( - 9e-11 ) , pores = BC3_pores ) #NEWLINE# g_super . append ( ( np . ones ( len ( BC3_pores ) ) * 1.5e-12 ) ) #NEWLINE# BC4_pores = np . arange ( 3 , 7 ) #NEWLINE# self . alg . set_boundary_conditions ( bctype = 'Neumann_group' , bcvalue = 1e-11 , pores = BC4_pores ) #NEWLINE# g_super . append ( np . array ( [ 6.42e-13 ] ) ) #NEWLINE# self . alg . run ( conductance = 'throat.cond' , quantity = 'pore.mole_fraction' , super_pore_conductance = g_super ) #NEWLINE# self . alg . return_results ( ) #NEWLINE# r1 = self . alg . rate ( BC1_pores ) [ 0 ] #NEWLINE# r2 = self . alg . rate ( BC2_pores ) [ 0 ] #NEWLINE# r3 = self . alg . rate ( BC3_pores ) [ 0 ] #NEWLINE# r4 = self . alg . rate ( BC4_pores ) [ 0 ] #NEWLINE# assert ( np . absolute ( ( ( ( r1 + r2 ) + r3 ) + r4 ) ) < 1e-20 ) #NEWLINE# assert ( np . size ( self . alg . super_pore_conductance [ 0 ] ) == 1 ) #NEWLINE# assert ( np . size ( self . alg . super_pore_conductance [ 1 ] ) == 7 ) #NEWLINE# assert ( np . size ( self . alg . super_pore_conductance [ 2 ] ) == 1 )