Télécharger la plaquette Cast3M


fr en

Télécharger avch3d.eso

Retour à la liste

Numérotation des lignes : 

avch3d

  1. C AVCH3D    SOURCE    PV090527  23/01/27    21:15:06     11574          
  2.       subroutine avch3d(pch,taar,nrjg,tref0,sr1,srsrag,teta1,
  3.      # teta2,dt,aar0,aar1,err1)
  4.  
  5. c      sous programme de calcul de l avancement chimique de rag
  6.  
  7.  
  8. c   ********************************************************************
  9.        implicit real*8 (a-h,o-z)
  10.        implicit integer (i-n)
  11. c   ********************************************************************
  12.  
  13.       real*8 pch,taar,nrjg,tref0,aar0,sr1,srsrag,teta1,teta2,dt     
  14.       real*8 alpharag,Ear,temp1,aar1,vvgel,cgel,rt
  15.       integer err1
  16.  
  17. c     constante cinetique a tref 
  18.       if((taar.ne.0.).and.(pch.ne.0.)) then
  19. c         temperature moyenne pendant la pas de temps
  20.           temp1=0.5d0*(teta1+teta2)
  21. c         la reaction est lente      
  22.           alpharag=taar**(-1)
  23. c         prise en compte de l activation thermique
  24.           Ear=nrjg/8.31d0
  25.           temp1=teta1+273.15d0
  26.           tempr=tref0+273.15d0
  27. c         activation thermique de la reaction      
  28.           AR=dEXP(-EaR*((1.d0/temp1)-(1.d0/tempr)))
  29.           alpharag=alpharag*AR
  30. c         calcul de l avancement chimique
  31.           if(aar0.lt.sr1) then 
  32.              if (sr1.gt.srsrag) then
  33.                  alpharag=alpharag*((sr1-srsrag)/(1.d0-srsrag))**2
  34.                  aar1=sr1-(sr1-aar0)*exp(-alpharag*dt)
  35.              else
  36.                  aar1=aar0
  37.              end if
  38.           else
  39.              aar1=aar0
  40.           end if
  41. c         la reaction ne peut pas rediminuer
  42.           aar1=max(aar0,aar1)          
  43.       else
  44. c         pas de donnee pour la reaction on la laisse a zero   
  45.           aar1=0.d0
  46.       end if
  47.  
  48.       return
  49.       end
  50.  
  51.  
  52.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales