Mathematics">
[go: up one dir, main page]

0% ont trouvé ce document utile (0 vote)
194 vues23 pages

Excel Pour TP

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1/ 23

La résolution des équations en génie chimique

En utilisant L’EXCEL
- Les équations linéaires:

- Les équation différentielles

- Modélisation des résultats expérimentaux


TP1: La résolution des équations algébriques

Les équations algébrique peuvent inclure:


- Polynômes: x2, x3, x4, x5, ,,,,,,,
- Log(x) or ln (x)
-Expenstiel exp(x)
-Fonctions trigonométriques: sin(x), cos(s), sinh(x), cosh(x)

Solution par « itérations »


La résolution des équation polynominales

La résolution de ces types d’équation s’effectue dans Microsoft Excel par itération en utilisant deux commandes:
- Valeur cible
- Solver

Exemple: équation de Van der Walls dans le cas de Benzène


1- en utilisant: la valeur cible
Dans le menu on va choisir : Données, Analyse scénario, valeur cible

L’ouverture de la petite fenêtre suivante:

La cellule à changer

La valeur cible
1- en utilisant: Solver
Dans le menu on va choisir : Données, solver

L’ouverture de la petite fenêtre suivante:

La cellule à changer

La valeur cible

Les contraintes pour les variables


La résolution d’un système d’équation

1- donner des valeurs aléatoires à x, y et z


2- calculer les résultats de ces valeurs pour chaque
équation Dans des cellules séparées,
3- calculer R2= (valeur voulue-valeur trouvée)2 dans une
colonne puis la somme de ces R2
4- dans la barre des menus on clique sur Données puis sur solver et une petite fenêtres s’ouvre

- Dans objet à définir on sélectionne la cellule de la somme


- Dans Valeur on donne la valeur de 0
- Dans cellules variables on sélectionne les cellules ou il y a
les valeur de x, y, et z
- Cliquer sur Résoudre
TP 2: Modélisation des résultats expérimentaux: fitting
Exemple1: fitting linéaire
Trouver des constantes d’une équation à partir des résultats expérimentaux

Trouver la valeur de K et de Ea:


1- calculer Ln(k) et 1/T
2- dessiner le graphe de ln(k)=f(1/T)

- Pour faire le fitting:


Dans la barre des menus:
Création,
Ajouter un élément de graphiques,
Courbe de tangence
Exemple2: fitting non linéaire (méthode des moindres carrés)
Trouver des constantes d’une équation à partir des résultats expérimentaux

Pour déterminer a, b, et c on va procéder comme suit:


1- donner des valeurs aléatoire à a, b et c
2- calculer les valeurs de ycal avec l’équation en utilisant
Les valeurs aléatoire,
N,B: pour fixer une cellule lors des calculer on ajoute
$ dans l’adresse de la cellule: $F$3
3- calculer R2
4- dans la barre du menu : Données puis Solver

l’objectif cible: R2
La valeur cible: 1
Les cellules variable: a, b, c
TP 3: opération sur les matrices et résolution des équations différentielles
Dégrée de liberté
- Indice la solvabilité du système

DOF=nombre de variables – nombre d’équations, relations et variables données

le système est solvable si seulement si DOF=0


Exemple:

DOF>0 le système ne peut être résolu

DOF<0 le système ne peut être résolu


Exemple 2:

Bilan de matière
Après substitution

De même pour chaque unité


opération sur les matrices pour résoudre un ensemble
d'équations algébriques linéaires
Avec X, Y et Z sont des variables à trouver
1- construire une matrice A et chaque colonne
va contenir les coefficient du même variable et
une matrice B avec les résultats des équations

2- calculer le déterminant de la matrice qui doit être diffèrent de 0

3- calculer l’inverse de la matrice A:


- sélectionner une plage de cellule qui doit être de même dimension que la matrice A
- Sur la première cellule on insère : puis appuyer sur
Ctrl+Maj+Entrée
4- calculer le produit de l’inverse de la matrice A avec la matrice B,
- Sélectionner une plage de cellule pour le produit: la matrice produit est une matrice
d’une seule colonne
- Sur la première cellule de la matrice produit on écrit: =produitmat(…….;………) puis
on appuie sur Ctrl+maj+entrée
La résolution de la matrice sur EXCEL
La résolution des équations différentielles
- donner une solution pour une équation différentielle spécifique pour des conditions initiales spécifiques et pour
une plage spécifique.
- célèbres techniques de résolution numérique:
méthode Euler
La méthode Runge-Kutta
méthode d'Euler modifiée
Qu’on va nommer f(x)
Runge-Kutta

C-à-dire que
On va poser: x=y et h=t dans la méthode
Exemple 1:
de Runge-Hutta

- On va calculer les ki et le H pour chaque


steps
- Nous avons une colonne pour t , une autre
pour H et 4 colonne pour les ki
- Dans chaque colonne on va introduire la
fonction correspondante
- On répète le calcul jusqu’à la valeur désiré
de H
Exemple 2:
TP5: les réacteurs idéaux homogènes
RAC et réacteur piston
I- les réacteurs continus homogènes parfaitement agités

Soit la réaction suivante mis en évidence dans un RAC d’ordre 1, H2 est en excès, le taux de conversion et
60%, le mélange sortant du réacteur va subir une distillation Flash comme indique le schéma,

A la sortie du Flash en 4, la composition du mélange est la suivante:

Dans l’Excel:
• Calculer le dégrée de liberté pour le réacteur, le Flash, le processus et le total et
monter que le système est solvable,
• Trouver la composition du mélange à la sortie du réacteur et dans la phase vapeur à
la sortie du Flash,
III- le réacteur piston
Soit la réaction en phase gaz, irréversible et du 2ere ordre :, A  2B
mise en œuvre dans un réacteur piston (S= 0,2m2 , L=10m). À la température T et pression p. la
constante de vitesse est égale à 0.69 min-1. le débit est 0,3m3/min

En utilisant le bilan de matière, monter que: la méthode de Runge-Kutta

𝑑𝐶𝐴 𝑘𝑆𝐶𝐴2
=−
𝑑𝐿 𝑄
* Posons C=Y et L=X
* En utilisant la méthode de Runge-Kutta,
determiner le concentration final de A ainsi
que le taux de conversion
I- les réacteurs continus homogènes parfaitement agités
III- le réacteur piston
TP6: les réacteurs idéaux homogènes
Bilan énergétique
Détermination de l'ordre de la vitesse d’une réaction:

L’action de l’iodure de méthyle sur la triméthylamine en solution dans le tétrachlorure de carbone est schématisée
par la réaction :
CH 3 I   C2 H 5 3 N 
 CH 3  N    C2 H 5 3   I

L’évolution l’amine en fonction du temps a été suivis pour deux cas :


1er cas : alimentation équimolaire : [amine] = CA0 = CB0 = 0,224 mol.L-1
2ème cas : [amine] = CA0 = 0,403 mol.L-1, CB0 = 0,046 mol.L-1

Cas 1:
• Donner l’expression de le concentration de A dans le cas ou l’ordre global est 1 et le cas ou n=2
• Dans l’Excel: calculer Ln({A] et 1/[A]
• Tracer les courbes ln([A])=f(t) et 1/{A]= f(t)
• Trouver l’équation de chaque courbes (fitting)
• Comparer les R2 obtenus, conclure
Cas 2
• Calculer les concentration de B et refaire les mêmes étapes précédentes
• Déterminer les ordres partiels de A et de B,
Détermination du point de fonctionnement d’un réacteurs adiabatique

Soit la réaction en phase liquide : A 


k1
B
mise en œuvre dans RAC adiabatique de volume 0.5L est alimenté avec un débit volumique 0.1L.min-1. La
concentration d’entrée et CA0=7mol.l-1, la température est de 75°C.
La capacité calorifique de cette solution est de 0.65 kJ.kg-1.K-1 et sa masse volumique est de 1.07 kg.L-1.

H r   2.4 KJ .mol 1 10000


k  33.10 exp(
12
) min 1
T
Montrer :
- En utilisant le bilan de matière que:

- En utilisant le bilan énergétique que:

Dans l’Excel:
- Tracer les courbes de XAS1=f(T) à déférentes temps de passage, conclure
- Tracer les deux courbes XAS1=f(T) et XAS2=f(T) et déterminer le point de fonctionnement du réacteur adiabatique

Vous aimerez peut-être aussi