Ci-dessous, les différences entre deux révisions de la page.
Les deux révisions précédentes Révision précédente Prochaine révision | Révision précédente | ||
wiki:projets:convection_geosciences [2019/04/17 16:21] aburza [Jour 3: 08/02] |
wiki:projets:convection_geosciences [2020/10/05 14:39] (Version actuelle) |
||
---|---|---|---|
Ligne 9: | Ligne 9: | ||
* Un rapide descriptif du projet: Le projet convection consiste à modéliser la convection mantellique à travers un modèle analogique | * Un rapide descriptif du projet: Le projet convection consiste à modéliser la convection mantellique à travers un modèle analogique | ||
* Le but et les objectifs de ce projet: Le but est de déterminer le nombre de Rayleigh critique (donc le moment où le fluide convecte) pour différentes huiles, en déterminant 3 coefficients à travers 3 expériences : la dilatation thermique, la diffusivité thermique et la viscosité cinématique. On prend pour vérifier nos manipulations expérimentales, | * Le but et les objectifs de ce projet: Le but est de déterminer le nombre de Rayleigh critique (donc le moment où le fluide convecte) pour différentes huiles, en déterminant 3 coefficients à travers 3 expériences : la dilatation thermique, la diffusivité thermique et la viscosité cinématique. On prend pour vérifier nos manipulations expérimentales, | ||
- | * La liste des machines et outils nécessaires: | + | * La liste des machines et outils nécessaires: |
- | * La liste du matériel nécessaire, | + | * La liste du matériel nécessaire, |
===== Logiciels utilisés ===== | ===== Logiciels utilisés ===== | ||
- | * MatLab R2018b | + | |
- | * Arduino | + | * Arduino |
- | * Origin 60 | + | * Origin 60 |
- | ===== Construction ===== | + | ===== Construction ===== |
+ | Notre stage portait sur le phénomène de convection et s’est déroulé de janvier à avril 2019, au sein du Fablab Géosciences et encadré par Loïc Labrousse et Pierre Thery. L’idée générale était d’élaborer des expériences nous permettant de calculer les coefficients intervenants dans le calcul du nombre de Rayleigh critique, nombre qui marque le début de la convection dans un fluide. Nos objectifs étaient donc de tester nos expériences et nos valeurs en prenant comme fluide témoin une huile de moteur dont on connaissait toutes les caractéristiques et d’en déduire un nombre de Rayleigh expérimentale à comparer à une valeur théorique. Le but final étant de pouvoir utiliser ces expériences pour des fluides se rapprochant du manteau terrestre. | ||
+ | |||
+ | __Théorie rapide :__ | ||
+ | On est parti de la formule de Rayleigh pour identifier les coefficients que nous devions retrouver pour calculer le nombre de Rayleigh critique d’un fluide. On peut mesurer h directement, | ||
+ | |||
+ | __1ère expérience :__ | ||
+ | Pour déterminer le coefficient de dilatation, nous avons mis en place une expérience afin de voir la variation de volume. Pour cela, nous avons utilisé un ballon rempli d’huile que l’on a chauffé en faisant un bain-marie. Comme on peut le voir sur la photo, nous avons disposé un tube en verre traversant le bouchon du ballon. Initialement, | ||
+ | Nous avions effectué l’expérience pour l’huile de moteur. On a fait la même expérience pour de l’huile de colza. Seulement, les données obtenues n’ont pas été concluantes : le coefficient de dilatation de l’huile de colza devait être trop grand pour nos ordres de grandeurs. | ||
+ | {{: | ||
+ | |||
+ | __2ème expérience :__ | ||
+ | Pour déterminer la viscosité cinématique nous avons décidé de passer par la vitesse de chute d’une bille dans notre huile, qui est reliée par la Loi de Stokes à la viscosité dynamique, on peut alors déterminer la viscosité cinématique à partir de la relation qui existe entre viscosité dynamique, masse volumique et viscosité cinématique. Pour obtenir cette vitesse nous avons pris en vidéo la chute d’une bille d' | ||
+ | |||
+ | {{: | ||
+ | |||
+ | __3ème expérience :__ | ||
+ | Pour la dernière expérience, | ||
+ | Nous avons pu mener l’expérience suivante : dans un bécher de 3L, nous avons introduit 1L d' | ||
+ | Ainsi donc, nous avons obtenu ce graphique présentant l’évolution de la température de nos sondes en fonction du temps. | ||
+ | Avec ces données nous avons été alors capable de déterminer notre coefficient de diffusivité thermique K. Pour cela, nous avons de nouveau eu recours à l’écriture d'un programme sur Matlab. | ||
+ | Tout d’abord, nous avons modélisé un bécher rempli d' | ||
+ | {{: | ||
+ | |||
+ | Nous sommes donc en possession de tous les paramètres nécessaires au calcul du nombre de Rayleigh critique. | ||
+ | Calcul Ra expérimental et comparaison théorique : | ||
+ | Finalement on a pu calculer un Ra expérimental et un Ra théorique et les comparer, nous avons trouvé pour Raexp= 623 et pour Rath= 1000, en encadrant notre Ra exp avec valeur max et min on trouve un Ra moyen de 584 ce qui n’est pas aberrant pour cette hauteur. Nos valeurs expérimentales sont toujours dans le bon ordre de grandeur, au final nos expériences sont donc à améliorer pour plus de précision mais elles sont tout de même fiables. | ||
+ | |||
+ | |||
+ | __Conclusion + implications aux échelles géologiques :__ | ||
+ | Maintenant que nous avons un protocole fiable pour mesurer les coefficients intervenant dans le calcul du nombre de Rayleigh pour un fluide donné, nous allons pouvoir tester des fluides dont les caractéristiques se rapprochent de celle du manteau terrestre. En effet aux échelles de temps géologiques, | ||
+ | la détermination du nombre de Rayleigh du manteau est plus difficile à mettre en place. En effet, h, α, ʋ et k sont facilement identifiables mais g n'est pas très quantifiable car il varie beaucoup avec la profondeur et on ne connaît pas le gradient de température du manteau, donc le delta T est difficile à évaluer aussi.\\ | ||
+ | === Poster du Stage : === | ||
+ | {{ : | ||
Ligne 47: | Ligne 81: | ||
{{: | {{: | ||
====Jour 4: 15/02==== | ====Jour 4: 15/02==== | ||
- | ====Jour 5: 22/02==== | ||
- | ====Jour 6: 01/03==== | ||
- | ====Jour 7: 08/03==== | ||
- | ====Jour 8: 15/03==== | ||
- | ====Jour 9: 22/03==== | ||
- | ====Jour 10: 29/03==== | ||
- | ====Jour 11: 05/04==== | ||
- | ====Jour 12: 12/04==== | ||
- | |||
- | |||
- | Présentation des projets et choix des binômes\\ | ||
- | |||
- | Audrey/ | ||
- | |||
- | Elsi/ | ||
- | |||
- | ====Jour 2: 01/02==== | ||
- | |||
- | |||
- | Théorie autour de la convection\\ | ||
- | |||
- | Début de recherche pour expériences détermination coefficients α , κ , η\\ | ||
- | |||
- | |||
- | ====Jour 3: 08/02==== | ||
- | |||
- | 1ers Crash Tests: Chloë avec Pierre:\\ | ||
- | expériences de convection cire gel et huile dans bécher 2L\\ | ||
- | |||
- | {{: | ||
- | |||
- | |||
- | ====Jour 4: 15/02==== | ||
- | |||
2ème Crash Tests: Elsi, Antoine et Chloë avec Loïc \\ | 2ème Crash Tests: Elsi, Antoine et Chloë avec Loïc \\ | ||
Ligne 116: | Ligne 116: | ||
{{ : | {{ : | ||
- | ====Jour 5: 22/02=== | + | ====Jour 5: 22/02==== |
- | + | Cf fin du journal | |
- | + | ||
- | Expérience 1 sur la détermination de α coefficient de dilatation thermique au Fablab Chimie avec Ziyad BOUAOUD\\ | + | |
- | + | ||
- | __Détermination du tube en verre adapté parmi 3 proposés: | + | |
- | + | ||
- | On cherche à avoir le tube en verre ayant le plus grand volume, nous permettant ainsi d' | + | |
- | + | ||
- | Caractéristiques des tubes en verre: diamètre (d) et longueur (l)\\ | + | |
- | + | ||
- | -Tube 1: d= 3.5mm ; l= 26.9 cm\\ | + | |
- | -Tube 2: d= 3.8mm ; l= 20.5 cm\\ | + | |
- | -Tube 3: d= 4.5mm ; l= 14.09 cm\\ | + | |
- | + | ||
- | On a déterminé le diamètre des tubes en verre à l'aide d'un pied à coulisse et sa longueur avec une règle\\ | + | |
- | L' | + | |
- | + | ||
- | + | ||
- | Calcul de la surface de la section | + | |
- | + | ||
- | S= πr^2 = π(d/ | + | |
- | + | ||
- | -S(tube 1)= π*(3.5/ | + | |
- | -S(tube 2)= π*(3.8/ | + | |
- | -S(tube 3)= π*(4.5/ | + | |
- | + | ||
- | Calcul | + | |
- | + | ||
- | ΔS(tube 1)= √(3.5(π/ | + | |
- | ΔS(tube 2)= √(3.8(π/ | + | |
- | ΔS(tube 3)= √(4.5(π/ | + | |
- | + | ||
- | Calcul du volume du tube en verre:\\ | + | |
- | plus précisement | + | |
- | V= Sl\\ | + | |
- | + | ||
- | -V(tube 1)= 9.62*269 = 2587.78 mm^3\\ | + | |
- | -V(tube 2)= 11.34*205 = 2324.7 mm^3\\ | + | |
- | -V(tube 3)= 15.90*140.9 = 2240.30 mm^3\\ | + | |
- | + | ||
- | Ainsi, le tube en verre 1 présente le volume le plus grand, c'est donc celui-ci que l'on choisira\\ | + | |
- | + | ||
- | + | ||
- | __Détermnation de la masse volumique de l' | + | |
- | + | ||
- | ρ= m/V\\ | + | |
- | + | ||
- | Placer une éprouvette graduée de 500 mL sur une balance et tarer\\ | + | |
- | Introduire un volume d' | + | |
- | + | ||
- | {{ : | + | |
- | + | ||
- | On a effectué cette manipulations 3 fois pour avoir une meilleur précision de ρ\\ | + | |
- | + | ||
- | -m1= 252.79 g V1= 300 cm^3 ===> ρ= 252.79/300 = 0.8426 g/cm^3\\ | + | |
- | -m2= 421.85 g V2= 500 cm^3 ===> ρ= 421.85/500 = 0.8437 g/cm^3\\ | + | |
- | -m3= 252.7 g V3= 300 cm^3 ===> ρ= 252.7/300 = 0.8457 g/cm^3\\ | + | |
- | + | ||
- | ρ(moyen)= 0.844 g/cm^3 à 22°C\\ | + | |
- | + | ||
- | L' | + | |
- | + | ||
- | __Matériel expérimental: | + | |
- | -Ballon\\ -Plaque-chauffante\\ -Élévateur\\ -Tube en verre de dimensions d= 3,5 mm, l= 26,9 cm pour l' | + | |
- | + | ||
- | __Protocole expérimental :__\\ | + | |
- | -Remplir le ballon jusqu’au niveau du bouchon\\ | + | |
- | -Mettre le bouchon\\ | + | |
- | -Faire attention à ce qu'il n'y ait pas de bulle car on veut calculer la variation du volume du liquide et non pas de l'air \\ | + | |
- | -Introduire le tube dans le bouchon jusqu' | + | |
- | -Introduire le thermomètre dans le bain-marie\\ | + | |
- | -Chauffer la plaque\\ | + | |
- | -Observer la variation de hauteur de l' | + | |
- | -Mettre en marche l' | + | |
- | -Faire des mesures tous les 5°C de la hauteur d' | + | |
- | + | ||
- | {{ : | + | |
- | //__Photo du montage expérimentale : observation de la variation du volume d' | + | |
- | + | ||
- | 1ère manipulation: | + | |
- | Masse ballon+bouchon+tube = 215.49 g\\ | + | |
- | Masse huile de colza = 245.71g\\ ERREUR | + | |
- | Température initiale = 20.2 °C | + | |
- | + | ||
- | + | ||
- | ==== Jour 6: 01/03 ==== \\ | + | |
+ | ====Jour 6: 01/03==== | ||
2ème manipulation: | 2ème manipulation: | ||
Masse ballon+bouchon+tube = 201.74 g\\ | Masse ballon+bouchon+tube = 201.74 g\\ | ||
Ligne 224: | Ligne 140: | ||
On obtient alors la valeur suivante. **alpha = 1.12x10^-3** \\ | On obtient alors la valeur suivante. **alpha = 1.12x10^-3** \\ | ||
Pour l' | Pour l' | ||
- | ==== Jour 7 : 08/03 ==== | ||
- | + | ====Jour 7: 08/03==== | |
Expérience 2 : détermination du coefficient de viscosité dans l' | Expérience 2 : détermination du coefficient de viscosité dans l' | ||
Ligne 240: | Ligne 155: | ||
- Prendre une vidéo \\ | - Prendre une vidéo \\ | ||
- Avec la pince, lacher la bille au plus près de la surface de l' | - Avec la pince, lacher la bille au plus près de la surface de l' | ||
+ | {{: | ||
- | + | ====Jour 8: 15/03==== | |
- | ====Jour 8 : 15/03==== | + | |
- | + | ||
- | + | ||
- Travail sur un programme pour pouvoir analyser les vidéos de la chute d'une bille dans l' | - Travail sur un programme pour pouvoir analyser les vidéos de la chute d'une bille dans l' | ||
Ligne 253: | Ligne 166: | ||
- Capteur thermique \\ | - Capteur thermique \\ | ||
- | + | ====Jour 9: 22/03==== | |
- | ====Jour 9 : 22/03==== | + | |
Ecriture du programme sur Matlab pour analyser les vidéos de la bille afin d' | Ecriture du programme sur Matlab pour analyser les vidéos de la bille afin d' | ||
Ligne 269: | Ligne 180: | ||
ρ = m/V = m(finale)-m(initiale)/ | ρ = m/V = m(finale)-m(initiale)/ | ||
+ | ====Jour 10: 29/03==== | ||
+ | Traitement des vidéos pour la viscosité avec programme écris avec Loïc | ||
+ | programme calcu_visc.: | ||
- | ==== Jour 10 : 29/03 ====\\ | + | ====Jour 11: 05/04==== |
- | + | ||
- | ==== Jour 11 : 05/04 ====\\ | + | |
- | + | ||
- | + | ||
- | \\ | + | |
Ecriture d'un programme permettant de déterminer une température à partir des données digitales fournies par une sonde.\\ | Ecriture d'un programme permettant de déterminer une température à partir des données digitales fournies par une sonde.\\ | ||
- | ==== Jour 12 : 12/04 ====\\ | + | ====Jour 12: 12/04==== |
- | \\ | + | |
Expérience 3 : détermination du coefficient de diffusivité thermique de l' | Expérience 3 : détermination du coefficient de diffusivité thermique de l' | ||
\\ | \\ | ||
Ligne 293: | Ligne 200: | ||
__Protocole expérimental :__\\ | __Protocole expérimental :__\\ | ||
- | Dans un bécher de 3L, introduire 1L d' | + | Dans un bécher de 3L, introduire 1L d' |
Chauffer l' | Chauffer l' | ||
{{: | {{: | ||
Ligne 309: | Ligne 216: | ||
Sur Matlab, on a modélisé un bécher rempli d' | Sur Matlab, on a modélisé un bécher rempli d' | ||
{{: | {{: | ||
+ | |||
+ | =====Programmes===== | ||
+ | {{ : | ||
+ | 1er programme: calcul de Rayleigh de façon théorique Matlab\\ | ||
+ | 2ème programme: traitement des images pour détermination viscosité cinématique Matlab\\ | ||
+ | 3ème programme: récupération des données de 3 capteurs thermiques Arduino\\ | ||
+ | 4ème programme: récupération des Températures Matlab\\ | ||
+ | ====Jour 5: 22/02==== | ||
+ | Expérience 1 sur la détermination de α coefficient de dilatation thermique au Fablab Chimie avec Ziyad BOUAOUD\\ | ||
+ | |||
+ | __Détermination du tube en verre adapté parmi 3 proposés: | ||
+ | |||
+ | On cherche à avoir le tube en verre ayant le plus grand volume, nous permettant ainsi d' | ||
+ | |||
+ | Caractéristiques des tubes en verre: diamètre (d) et longueur (l)\\ | ||
+ | |||
+ | -Tube 1: d= 3.5mm ; l= 26.9 cm\\ | ||
+ | -Tube 2: d= 3.8mm ; l= 20.5 cm\\ | ||
+ | -Tube 3: d= 4.5mm ; l= 14.09 cm\\ | ||
+ | |||
+ | On a déterminé le diamètre des tubes en verre à l'aide d'un pied à coulisse et sa longueur avec une règle\\ | ||
+ | L' | ||
+ | |||
+ | |||
+ | Calcul de la surface de la section du tube en verre:\\ | ||
+ | |||
+ | S= πr^2 = π(d/2)^2\\ | ||
+ | |||
+ | -S(tube 1)= π*(3.5/ | ||
+ | -S(tube 2)= π*(3.8/ | ||
+ | -S(tube 3)= π*(4.5/ | ||
+ | |||
+ | Calcul de ΔS:\\ | ||
+ | |||
+ | ΔS(tube 1)= √(3.5(π/ | ||
+ | ΔS(tube 2)= √(3.8(π/ | ||
+ | ΔS(tube 3)= √(4.5(π/ | ||
+ | |||
+ | Calcul du volume du tube en verre:\\ | ||
+ | plus précisement | ||
+ | V= Sl\\ | ||
+ | |||
+ | -V(tube 1)= 9.62*269 = 2587.78 mm^3\\ | ||
+ | -V(tube 2)= 11.34*205 = 2324.7 mm^3\\ | ||
+ | -V(tube 3)= 15.90*140.9 = 2240.30 mm^3\\ | ||
+ | |||
+ | Ainsi, le tube en verre 1 présente le volume le plus grand, c'est donc celui-ci que l'on choisira\\ | ||
+ | |||
+ | |||
+ | __Détermnation de la masse volumique de l' | ||
+ | |||
+ | ρ= m/V\\ | ||
+ | |||
+ | Placer une éprouvette graduée de 500 mL sur une balance et tarer\\ | ||
+ | Introduire un volume d' | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | On a effectué cette manipulations 3 fois pour avoir une meilleur précision de ρ\\ | ||
+ | |||
+ | -m1= 252.79 g V1= 300 cm^3 ===> ρ= 252.79/300 = 0.8426 g/cm^3\\ | ||
+ | -m2= 421.85 g V2= 500 cm^3 ===> ρ= 421.85/500 = 0.8437 g/cm^3\\ | ||
+ | -m3= 252.7 g V3= 300 cm^3 ===> ρ= 252.7/300 = 0.8457 g/cm^3\\ | ||
+ | |||
+ | ρ(moyen)= 0.844 g/cm^3 à 22°C\\ | ||
+ | |||
+ | L' | ||
+ | |||
+ | __Matériel expérimental: | ||
+ | -Ballon\\ -Plaque-chauffante\\ -Élévateur\\ -Tube en verre de dimensions d= 3,5 mm, l= 26,9 cm pour l' | ||
+ | |||
+ | __Protocole expérimental :__\\ | ||
+ | -Remplir le ballon jusqu’au niveau du bouchon\\ | ||
+ | -Mettre le bouchon\\ | ||
+ | -Faire attention à ce qu'il n'y ait pas de bulle car on veut calculer la variation du volume du liquide et non pas de l'air \\ | ||
+ | -Introduire le tube dans le bouchon jusqu' | ||
+ | -Introduire le thermomètre dans le bain-marie\\ | ||
+ | -Chauffer la plaque\\ | ||
+ | -Observer la variation de hauteur de l' | ||
+ | -Mettre en marche l' | ||
+ | -Faire des mesures tous les 5°C de la hauteur d' | ||
+ | |||
+ | {{ : | ||
+ | //__Photo du montage expérimentale : observation de la variation du volume d' | ||
+ | |||
+ | 1ère manipulation: | ||
+ | Masse ballon+bouchon+tube = 215.49 g\\ | ||
+ | Masse huile de colza = 245.71g\\ | ||
+ | Température initiale = 20.2°C\\ | ||
+ | |||
Ligne 314: | Ligne 311: | ||
Et surtout, pensez à mettre des photos. | Et surtout, pensez à mettre des photos. | ||
</ | </ | ||
+ |