<simulation>
  
<header>
    
<title lang="french">Réaction réversible simple (équilibre)</title>
    
<title lang="english">Reversible reaction(equilibrium)</title>
    
<title lang="spanish">Reacción reversible sencilla (equilibrio)</title>
    
<author>André Pauss, Stéphane Mottelet, Guillaume Barbaud</author>
    
<date>September 1st, 2005</date>
    
<keywords>Simulación, scilab, xml</keywords>
    
<keywords lang="spanish">simulacion, scilab, xml</keywords>
    
<image href="Reversible_reaction.gif" />
  
</header>
  
<notes lang="french">
    
<p>Cette simulation permet de visualiser l'évolution des concentrations
    des réactifs et des produits pour une réaction réversible.
</p>
  
</notes>
  
<notes lang="english">
    
<p>This simulation allows to examine the evolution of the concentrations
    of reactants and products in a reversible reaction.
</p>
  
</notes>
  
<notes lang="spanish">
    
<p>Esta simulación permite visualizar la evolución de las concentraciones de los   reactivos y productos para una reacción reversible</p>
  
</notes>
  
<parameters><section>
      
<title lang="french">Paramètres des constituants A et B</title>
      
<title lang="english">Parameters of reactants A and B</title>
      
<title lang="spanish">Parámetros de los constituyentes A y B</title>
      
<subsection>
        
<scalar label="nuA" scale="1" state="normal"><name lang="french">Coefficient stoechiométrique de A</name><name lang="english">Stoechiometric coefficient of A</name><name lang="spanish">Coeficiente estoiquiometro de A</name><value>1</value></scalar>
        
<scalar increment="1" label="betaA" max="2" min="1" scale="1" state="normal"><name lang="french">Ordre partiel de A (1 ou 2)</name><name lang="english">Kinetic order of A (1 or 2)</name><name lang="spanish">Orden parcial de A (1 o 2)</name><value>1</value></scalar>
      
</subsection>
      
<subsection>
        
<scalar label="nuB" scale="1" state="normal"><name lang="french">Coefficient stoechiométrique de B</name><name lang="english">Stoechiometric coefficient of B</name><name lang="spanish">Coeficiente estoiquiometro de B</name><notes lang="french">
            
<p>Si vous voulez simuler une équation simple avec un seul réactif
            (A donnant R et S), attribuez au réactif B un coefficient
            stoechiométrique, un ordre partiel ET une concentration nulle.
</p>
          
</notes><notes lang="english">
            
<p>If you want simulate a reaction with a single reactant (A
            producing products R and S), please set the stoechiometric
            coefficient, the kinetic order, AND the initial concentration of B
            to zero.
</p>
          
</notes><notes lang="spanish">
            
<p>Si usted desea simular una ecuación sencilla con uno solo reactivo     (A dando R y S), de al reactivo B un coeficiente    estoiquiometro, un orden parcial Y una concentración igual a 0.</p>
          
</notes><value>1</value></scalar>
        
<scalar increment="1" label="betaB" max="2" min="0" scale="1" state="normal"><name lang="french">Ordre partiel de B (0,1 ou 2)</name><name lang="english">Kinetic order of B (0, 1 or 2)</name><name lang="spanish">Orden parcial de B (0,1 o 2)</name><notes lang="french">
            
<p>Si vous voulez simuler une équation simple avec un seul réactif
            (A donnant R et S), attribuez au réactif B un coefficient
            stoechiométrique, un ordre partiel ET une concentration nulle.
</p>
          
</notes><notes lang="english">
            
<p>If you want simulate a reaction with a single reactant (A
            producing products R and S), please set the stoechiometric
            coefficient, the kinetic order, AND the initial concentration of B
            to zero.
</p>
          
</notes><notes lang="spanish">
            
<p>Si usted desea simular una ecuación sencilla con uno solo reactivo     (A dando R y S), de al reactivo B un coeficiente    estoiquiometro, un orden parcial Y una concentración igual a 0.</p>
          
</notes><value>1</value></scalar>
      
</subsection>
    
</section><section>
      
<title lang="french">Paramètres des constituants R et S</title>
      
<title lang="english">Paremters of products R and S</title>
      
<title lang="spanish">Parámetros de los constituyentes R et S</title>
      
<subsection>
        
<scalar label="nuR" scale="1" state="normal"><name lang="french">Coefficient stoechiométrique de R</name><name lang="english">Stoechiometric coefficient of R</name><name lang="spanish">Coeficiente estoiquiometro de R</name><value>1</value></scalar>
        
<scalar increment="1" label="betaR" max="2" min="1" scale="1" state="normal"><name lang="french">Ordre partiel de R (1 ou 2)</name><name lang="english">Kinetic order of R (1 or 2)</name><name lang="spanish">Orden parcial de R (1 o 2)</name><value>1</value></scalar>
      
</subsection>
      
<subsection>
        
<scalar label="nuS" scale="1" state="normal"><name lang="french">Coefficient stoechiométrique de S</name><name lang="english">Stoechiometric coefficient of S</name><name lang="spanish">Coeficiente estoiquiometro de S</name><notes lang="french">
            
<p>Si vous voulez simuler une équation simple avec un seul produit
            (A et B donnant R), attribuez au produit S un coefficient
            stoechiométrique, un ordre partiel ET une concentration nulle.
</p>
          
</notes><notes lang="english">
            
<p>If you want simulate a reaction with a single product (A
            and B producing product R), please set the stoechiometric
            coefficient, the kinetic order, AND the initial concentration of S
            to zero.
</p>
          
</notes><notes lang="spanish">
            
<p>Si usted desea simular una ecuación sencilla con uno solo producto    (A y B dan R), de al producto S un coeficiente    estoiquiometro, un orden parcial Y una concentración igual a 0.</p>
          
</notes><value>1</value></scalar>
        
<scalar increment="1" label="betaS" max="2" min="0" scale="1" state="normal"><name lang="french">Ordre partiel de S (0,1 ou 2)</name><name lang="english">Kinetic order of S (0, 1 or 2)</name><name lang="spanish">Orden parcial de S(0,1 o 2)</name><notes lang="french">
            
<p>Si vous voulez simuler une équation simple avec un seul produit
            (A et B donnant R), attribuez au produit S un coefficient
            stoechiométrique, un ordre partiel ET une concentration nulle.
</p>
          
</notes><notes lang="english">
            
<p>If you want simulate a reaction with a single product (A
            and B producing product R), please set the stoechiometric
            coefficient, the kinetic order, AND the initial concentration of S
            to zero.
</p>
          
</notes><notes lang="spanish">
            
<p>Si usted desea simular una ecuación sencilla con uno solo producto     (A y B dan R), de al producto S un coeficiente    estoiquiometro, un orden parcial Y una concentración igual a 0.</p>
          
</notes><value>1</value></scalar>
      
</subsection>
    
</section><section>
      
<title lang="french">Concentration initiale des constituants</title>
      
<title lang="english">Initial concentration of components</title>
      
<title lang="spanish">Concentración inicial de los constituyentes</title>
      
<subsection>
        
<scalar increment="0.1" label="cA0" max="1" min="0" scale="1" state="normal" unit="mol / L"><name lang="french">Concentration initiale en A (mol/L)</name><name lang="english">Initial concentration of A (mol/L)</name><name lang="spanish">Concentración inicial de A (mol/L)</name><value>1</value></scalar>
        
<scalar increment="0.1" label="cB0" max="1" min="0" scale="1" state="normal" unit="mol / L"><name lang="french">Concentration initiale en B (mol/L)</name><name lang="english">Initial concentration of B (mol/L)</name><name lang="spanish">Concentración inicial de B (mol/L)</name><notes lang="french">
            
<p>Si vous voulez simuler une équation simple avec un seul réactif
            (A donnant R et S), attribuez au réactif B un coefficient
            stoechiométrique, un ordre partiel ET une concentration nulle.
</p>
          
</notes><notes lang="english">
            
<p>If you want simulate a reaction with a single reactant (A
            producing products R and S), please set the stoechiometric
            coefficient, the kinetic order, AND the initial concentration of B
            to zero.
</p>
          
</notes><notes lang="spanish">
            
<p>Si usted desea simular una ecuación sencilla con uno solo reactivo     (A dando R y S), de al reactivo B un coeficiente    estoiquiometro, un orden parcial Y una concentración igual a 0.</p>
          
</notes><value>0.8</value></scalar>
      
</subsection>
      
<subsection>
        
<scalar increment="0.1" label="cR0" max="1" min="0" scale="1" state="normal" unit="mol / L"><name lang="french">Concentration initiale en R (mol/L)</name><name lang="english">Initial concentration of R (mol/L)</name><name lang="spanish">Concentración inicial de R (mol/L)</name><value>0.1</value></scalar>
        
<scalar increment="0.1" label="cS0" max="1" min="0" scale="1" state="normal" unit="mol / L"><name lang="french">Concentration initiale en S (mol/L)</name><name lang="english">Initial concentration of S (mol/L)</name><name lang="spanish">Concentración inicial de S (mol/L)</name><notes lang="french">
            
<p>Si vous voulez simuler une équation simple avec un seul produit
            (A et B donnant R), attribuez au produit S un coefficient
            stoechiométrique, un ordre partiel ET une concentration nulle.
</p>
          
</notes><notes lang="english">
            
<p>If you want simulate a reaction with a single product (A
            and B producing product R), please set the stoechiometric
            coefficient, the kinetic order, AND the initial concentration of S
            to zero.
</p>
          
</notes><notes lang="spanish">
            
<p>Si usted desea simular una ecuación sencilla con uno solo producto    (A y B dan R), de al producto S un coeficiente    estoiquiometro, un orden parcial Y una concentración igual a 0.</p>
          
</notes><value>0</value></scalar>
      
</subsection>
    
</section><section>
      
<title lang="french">Constantes de vitesse</title>
      
<title lang="english">Rate constants</title>
      
<title lang="spanish">Constantes de velocidad</title>
      
<subsection>
        
<scalar increment="1" label="Ea1" max="200" min="1" unit="kJ / mol" state="normal"><name lang="french">Première énergie d'activation (kJ/mol)</name><name lang="english">First activation energy (kJ/mol)</name><name lang="spanish">Primera energía de activación (kJ/mol)</name><notes lang="french">
            
<p> Réaction de A et B vers R et S, de gauche à droite</p>
          
</notes><notes lang="english">
            
<p> Reaction of A and B to R and S, from left to right</p>
          
</notes><notes lang="spanish">
            
<p>Reacción de A y B hacia R y S, de la izquierda a la derecha</p>
          
</notes><value>50</value></scalar>
        
<scalar label="k01" state="normal"><name lang="french">Premier facteur préexponentiel ko </name><name lang="english">First k0 factor </name><name lang="spanish">Primer factor pre exponencial k0</name><notes lang="french">
            
<p> Réaction de A et B vers R et S, de gauche à droite</p>
          
</notes><notes lang="english">
            
<p> Reaction of A and B to R and S, from left to right</p>
          
</notes><notes lang="spanish">
            
<p>Reacción de A y B hacia R y S, de la izquierda a la derecha</p>
          
</notes><value>2E9</value></scalar>
      
</subsection>
      
<subsection>
        
<scalar increment="1" label="Eaminus1" max="200" min="1" unit="kJ / mol" state="normal"><name lang="french">Seconde énergie d'activation (kJ/mol)</name><name lang="english">Second activation energy (kJ/mol)</name><name lang="spanish">Segunda energía de activación (kJ/mol)</name><notes lang="french">
            
<p> Réaction de R et S vers A et B, de droite à gauche</p>
          
</notes><notes lang="english">
            
<p> Reaction of R and S to A and B, from right to left</p>
          
</notes><notes lang="spanish">
            
<p>Reacción de R y S hacia A y B, de la derecha a la izquierda</p>
          
</notes><value>70</value></scalar>
        
<scalar label="k0minus1" state="normal"><name lang="french">Second facteur préexponentiel ko </name><name lang="english">Second k0 factor </name><name lang="spanish">Segundo factor pre exponencial k0</name><notes lang="french">
            
<p> Réaction de R et S vers A et B, de droite à gauche</p>
          
</notes><notes lang="english">
            
<p> Reaction of R and S to A and B, from right to left</p>
          
</notes><notes lang="spanish">
            
<p>Reacción de R y S hacia A y B, de la derecha a la izquierda</p>
          
</notes><value>2E8</value></scalar>
      
</subsection>
    
</section><section>
      
<title lang="french">Température et conditions de simulation</title>
      
<title lang="english">Temperature and simulation parameters</title>
      
<title lang="spanish">Temperatura y condiciones de simulación</title>
      
<scalar label="T" scale="1" state="normal" unit="K"><name lang="french">Température (K)</name><name lang="english">Temperature (K)</name><name lang="spanish">Temperatura (K)</name><value>298</value></scalar>
      
<scalar increment="1" label="tf" max="20" min="0" scale="1" state="normal" unit=""><name lang="french">Temps final</name><name lang="english">Final time</name><name lang="spanish">Tiempo final</name><value>4</value></scalar>
      
<scalar increment="50" label="simulation_step" max="1000" min="100" scale="1" state="normal" unit=""><name lang="french">Nombre de pas de simulation</name><name lang="english">Number of simulation steps</name><name lang="spanish">Numero de pasos de simulación</name><value>200</value></scalar>
    
</section><section>
      
<title lang="french">Points expérimentaux</title>
      
<title lang="english">Experimental data</title>
      
<title lang="spanish">Puntos experimentales</title>
      
<matrix clear="yes" cols="2" label="user_curve_mat" load="yes" rows="64" save="yes" stripcol="no" striprow="yes" widget="normal">
        
<name />
        
<name lang="spanish" />
        
<col>
          
<name lang="french">Abscisses</name>
          
<name lang="english">X-axis data</name>
          
<name lang="spanish">Coordenadas X</name>
          
<value>0</value>
        
</col>
        
<col>
          
<name lang="french">Ordonnées</name>
          
<name lang="english">Y-axis data</name>
          
<name lang="spanish">Coordenadas Y</name>
          
<value>0</value>
        
</col>
      
</matrix>
    
</section></parameters>
  
<compute>
    
<defdomain1d label="t" unit="">
      
<name lang="french">Temps</name>
      
<name lang="english">Time</name>
      
<name lang="spanish">Tiempo</name>
      
<interval discretization="linear" steps="100">
        
<initialvalue>0</initialvalue>
        
<finalvalue>tf</finalvalue>
      
</interval>
    
</defdomain1d>
    
<ode label="reversible">
      
<refdomain1d ref="t" />
      
<states>
        
<variable label="k1">k01*exp(-Ea1*1000/T/8.3144)</variable>
        
<variable label="kminus1">k0minus1*exp(-Eaminus1*1000/T/8.3144)</variable>
        
<state label="cA" unit="mol/l" size="1">
          
<name lang="french">Concentration de A</name>
          
<name lang="english">Concentration of A</name>
          
<name lang="spanish">Concentración de A</name>
          
<derivative>-nuA*((k1*cA^betaA*cB^betaB) -
          (kminus1*cR^betaR*cS^betaS))
</derivative>
          
<initialcondition>cA0</initialcondition>
        
</state>
        
<state label="cB" unit="mol/l" size="1">
          
<name lang="french">Concentration de B</name>
          
<name lang="english">Concentration of B</name>
          
<name lang="spanish">Concentración de B</name>
          
<derivative>-nuB*((k1*cA^betaA*cB^betaB) -
          (kminus1*cR^betaR*cS^betaS))
</derivative>
          
<initialcondition>cB0</initialcondition>
        
</state>
        
<state label="cR" unit="mol/l" size="1">
          
<name lang="french">Concentration de R</name>
          
<name lang="english">Concentration of R</name>
          
<name lang="spanish">Concentración de R</name>
          
<derivative>nuR*((k1*cA^betaA*cB^betaB) -
          (kminus1*cR^betaR*cS^betaS))
</derivative>
          
<initialcondition>cR0</initialcondition>
        
</state>
        
<state label="cS" unit="mol/l" size="1">
          
<name lang="french">Concentration de S</name>
          
<name lang="english">Concentration of S</name>
          
<name lang="spanish">Concentración de S</name>
          
<derivative>nuS*((k1*cA^betaA*cB^betaB) -
          (kminus1*cR^betaR*cS^betaS))
</derivative>
          
<initialcondition>cS0</initialcondition>
        
</state>
      
</states>
    
</ode>
  
</compute>
  
<graphs>
    
<parametriccurve2d label="user_curve">
      
<name lang="french">Points expérimentaux</name>
      
<name lang="english">Experimental data</name>
      
<name lang="spanish">Puntos experimentales</name>
      
<x1>
        
<name />
        
<name lang="spanish" />
        
<value>user_curve_mat(:,1)</value>
      
</x1>
      
<x2>
        
<name />
        
<name lang="spanish" />
        
<value>user_curve_mat(:,2)</value>
      
</x2>
    
</parametriccurve2d>
  
</graphs>
  
<display>
    
<window splitx="1" splity="1">
      
<title lang="french">Résultats</title>
      
<title lang="english">Results</title>
      
<title lang="spanish">Resultados</title>
      
<axis2d xmax="tf" xmin="0" ymax="cA0+cB0+cR0+cS0" ymin="0" type="cartesian" iso="no" position="left" colorbar="off">
        
<drawcurve2d color="red" ref="cA" thickness="1" />
        
<drawcurve2d color="green" ref="cB" thickness="1" />
        
<drawcurve2d color="blue" ref="cR" thickness="1" />
        
<drawcurve2d color="black" ref="cS" thickness="1" />
        
<drawcurve2d color="auto" linetype="solid" marker="plus" ref="user_curve" thickness="1" />
      
</axis2d>
    
</window>
  
</display>
  
<save>
    
<file format="csv" href="essai.csv" labels="nuA nuB nuR nuS betaA betaB betaR betaS Ea1 k01 Eaminus1 k0minus1 T t cA cB cR cS" />
  
</save>
</simulation>