Page 69 - dysertaciyahembara

1 1 1
8

 V
  e
i
m
  
k 2ij   D T H    N j  1  hm    I 11  N q  i  I 12  N s  i  I 13  N  I 11  N q  m  I 12  N s  m  I 13  N    

p 
p 
RT

m
  N  N 1 1 1  N   N  N  N 
  I 21 i  I 22 i  I 23 i  I 21 m  I 22 m  I 23 m   (2.18)
 q  s  p   q  s  p  
  N  N  N   N  N  N  
  I 31 i  I 32 i  I 33 i  I 31 m  I 32 m  I 33 m   det J dqdsdp
 q  s  p   q  s  p  
де
 8  N 8  N 8  N 
 x  y  z    i x i  i y i  i z i 
 q  q  q    i 1 q  i 1 q  i 1 q  
   
8
8
8
  J   x  y  z       N i x i   N i y i   N i z i  (2.19)


 s  s  s    i 1 s  i 1 s  i 1 s  


 x  y  z    8  N 8  N 8  N 
 p  p   p    i x i  i y i  i z i 


 i 1 p  i 1 p  i 1 p  

I – елементи матриці оберненої до (2.19).
ij
Для обчислення розподілу концентрації водню складено програму

для реалізації МСЕ відповідно до розробленого алгоритму [14].

2.3.3. Апробація алгоритму і програми розв’язання задачі дифузії

водню. Апробацію проведено на прикладі розв’язування задачі дифузії



водню (2.7) у паралелепіпеді l   1 , l l   2 , l l   3 , l (рис. 2.4)
z
y
x
1
3
2
за таких початкових і крайових умов

C x,y,z,  0  C ;
0
C  l ,y,z,t 1  C x, l ,z,t   2  C x,y, l ,t   3  C , 1 (2.20)

69