<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>