« "JBSIM 1.3 BY J.B." { "Calcul d'un pH" "Dosage" } 1. IF CHOOSE THEN { "Calcul d'un pH" "Dosage" } SWAP POS { Vd VM NVAC V H CALCpH EQPH } PURGE { AcideFort « "Acide Fort" { { "C" "Concentration en acide fort" 0. 6. } } 1. { } { } IF INFORM THEN OBJ→ DROP END -1. 0. » BaseForte « "Base Forte" { { "C" "Concentration en base forte" 0. 6. } } 1. { } { } IF INFORM THEN OBJ→ DROP END 1. 0. » Acidite1 « "1 ACIDITE" { { "CA" "Concentration en [AH]" 0. 6. } { "CB" "Concentration en [A-]" 0. 6. } { "pKa" "-LOG(Constante d'acidite)" 0. 6. } } 1. { } { } IF INFORM THEN OBJ→ DROP END → CA CB pKa « H CB * CA pKa NEG ALOG * - H pKa NEG ALOG + / 1. » 0. » Acidites2 « "2 ACIDITES" { { "C1" "Concentration en AH2" 0. 6. } { "C2" "Concentration en AH-" 0. 6. } { "C3" "Concentration en A2-" 0. 6. } { "pKa1" "-LOG(Constante d'acidite 1)" 0. 6. } { "pKa2" "-LOG(Constante d'acidite 2)" 0. 6. } } 2. { } { } IF INFORM THEN OBJ→ DROP END → C1 C2 C3 pKa1 pKa2 « C2 2. C1 * + H C3 * H C1 * - pKa2 ALOG * - 2. H 2. ^ * C3 * H 2. ^ C2 * + pKa2 pKa1 + ALOG * - NEG 1. H pKa2 ALOG * + H 2. ^ pKa2 pKa1 + ALOG * + / 1. » 0. » Acidites3 « "3 ACIDITES" { { "C1" "Concentration en H3A" 0. 6. } { "C2" "Concentration en H2A-" 0. 6. } { "C3" "Concentration en HA2-" 0. 6. } { "C4" "Concentration en A3-" 0. 6. } { "pKa1" "-LOG(Constante d'acidite 1)" 0. 6. } { "pKa2" "-LOG(Constante d'acidite 2)" 0. 6. } { "pKa3" "-LOG(Constante d'acidite 3)" 0. 6. } } 2. { } { } IF INFORM THEN OBJ→ DROP END → C1 C2 C3 C4 pKa1 pKa2 pKa3 « -1. 10. pKa1 NEG pKa2 - pKa3 - ^ 10. pKa2 NEG pKa1 - ^ H * + 10. pKa1 NEG ^ H 2. ^ * + H 3. ^ + / 3. 10. pKa2 NEG pKa3 - ^ * 2. 10. pKa2 NEG ^ * H * + H 2. ^ + * C1 C2 + C3 + C4 + * 10. pKa1 NEG ^ * C2 + 2. C3 * + 3. C4 * + 1. » 0. » CH1 « { AcideFort BaseForte Acidite1 Acidites2 Acidites3 Continuer } 1. » Calcul1pH « { Vd 10. } | { V 0. } | 'H' { 0. .0000001 1.E-14 } ROOT LOG NEG "pH" →TAG 1. » Courbe « (0.,-2.) VM SORT DUP SIZE GET 14. R→C 'V' # 2h { (0.,0.) { 1. 1. } } { FUNCTION } OBJ→ DROP 'Y' 7. →LIST 'PPAR' STO « 'EQPH' 'H' 0. H 1.E-14 3. →LIST ROOT LOG NEG » 'EQ' STO PICT PURGE DRAX DRAW » } → Ch LISTE « 'H-1.E-14/H' { « DO LISTE DUP CH1 POS 1. + GET "BECHER" SWAP EVAL IF CHOOSE THEN IF DUP Continuer SAME THEN DROP LISTE DUP Calcul1pH POS 1. + GET EVAL ELSE LISTE DUP 3. ROLL POS 1. + GET EVAL 3. ROLL 3. ROLL * 3. ROLL + SWAP END ELSE DROP 1. END UNTIL END » « DO LISTE DUP CH1 POS 1. + GET "BECHER" SWAP EVAL IF CHOOSE THEN IF DUP Continuer SAME THEN DROP DO LISTE DUP CH1 POS 1. + GET "BURETTE" SWAP EVAL IF CHOOSE THEN IF DUP Continuer SAME NOT THEN LISTE DUP 3. ROLL POS 1. + GET EVAL 3. ROLL V * Vd V + / 3. ROLL * 3. ROLL + SWAP ELSE DROP "VOLUME" { { "LIST V" "LISTE DES PH A CALCULER" 5. } { "Vd" "VOLUME DE SOLUTION A DOSER" 0. 6. } } 1. { } { } IF INFORM THEN OBJ→ DROP 'Vd' STO DUP 'VM' STO SIZE 'NVAC' STO EVAL 'EQPH' STO 1. 1. ELSE DROP 1. END END ELSE DROP 1. END UNTIL END ELSE LISTE DUP 3. ROLL POS 1. + GET EVAL 3. ROLL Vd * Vd V + / 3. ROLL * 3. ROLL + SWAP END ELSE DROP 1. END UNTIL END 0. 'H' STO « 'EQPH' 'H' 0. H 1.E-14 3. →LIST ROOT LOG NEG » 'CALCpH' STO " +--------------+" 3. DISP " | Calcul de pH |" 4. DISP " | en cours ... |" 5. DISP " +--------------+" 6. DISP 1. NVAC FOR M VM M GET DUP 'V' STO CALCpH NEXT NVAC 2. 2. →LIST →ARRY "Courbe ?" { "Oui" "Non" } 1. IF CHOOSE THEN IF "Oui" SAME THEN LISTE DUP Courbe POS 1. + GET EVAL END END { Vd VM NVAC } PURGE » } Ch GET EVAL » END »
On télécharge le répertoire JBSIM.DIR (ou les fichiers séparément)
On tape INSTALL.
Le programme demande de choisir le port d'installation (non protégé en écriture).
JBSim se contente de résoudre l'équation en fonction de [H+] (la concentration en ion hydrogène).
Pour une solution d'eau distillée on résout l'équation suivante :
[H+]=[OH-] ⇐⇒ [H+]-10E-14/[H+]=0
Si on rajoute un monoacide (HCl par exemple), il suffit de rajouter [Cl-]
[H+]=[OH-]+[Cl-] ⇐⇒ [H+]-10E-14/[H+]-Co=0
ou Co est la concentration de l'acide.
Chaque fois que l'on rajoute une espèce chimique dans la solution, on rajoute un terme a l'équation d'électroneutralité.
Le programme a été “code” avec l'utilitaire CODE de Levy après avoir été compressé avec PK de Mika.
On lance le programme :
Si on introduit dans le bêcher un acide fort (HCl), on procède ainsi :
Le pH est alors calculé et le résultat apparaît sur la pile.
Si on introduit dans le bêcher un acide faible (CH3COOH), on procède ainsi :
Le pH est alors calcule et le résultat apparaît sur la pile.
Si on introduit dans le bêcher une base faible (NH3), on procède ainsi :
Si on introduit dans le bêcher une solution tampon (CH3COOH 0.1 mol/L et CH3COO- 0.1 mol/L), on procède ainsi :
Le pH est alors calcule et le résultat apparaît sur la pile.
Si on introduit dans le bêcher un mélange d'acide et de base forts (HCl et NaOH), on procède ainsi :
Le pH est alors calcule et le résultat apparaît sur la pile.
On rajoute dans la solution tout ce que l'on veut puis on Continue …
On opère de la même façon que pour le calcul d'un pH :
Ex : {0 5 10 15 20 25 30 35 40}
le pH sera calculé pour
V = 0, 5, 10, 20, 35, 30, 35, 40 mL
et la courbe sera tracée pour un volume total versé de 40mL.
Si on met {20} un seul pH sera calculé a 20mL et la courbe sera tracée de 0 a 20mL.
NB : les volumes doivent être exprimés dans la même unité (mL ou L).