| 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/05/20 15:40]
aburza [Construction] |
wiki:projets:convection_geosciences [2020/10/05 14:39] (Version actuelle)
|
| ===== Logiciels utilisés ===== | ===== Logiciels utilisés ===== |
| |
| * MatLab R2018b | * 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. | 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. |
| 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 une certaine hauteur d’huile. Pour déterminer le coefficient de dilatation, nous connaissons cette équation : rho = rho0*(1-alphadelta T). On isole alors alpha. L’équation devient : alpha = Sdh(T-T0)/(Sdh+V0). Nous avons besoin de connaître la section du tube, la température et la hauteur à chaque instant, et le volume initial contenu dans le ballon. Pour calculer la pente, nous avons donc calculé la section du tube en verre et retrouver le volume initial par la mesure de la masse de l’huile. Pour connaître la température, nous avons utilisé un thermomètre mis dans le bain-marie. Et pour la hauteur, nous avions placé une règle graduée. L’expérience a été filmée, permettant ainsi de visionner et d’attribuer la hauteur du tube avec sa température correspondante. Nous avons tracé le graphique, et obtenu cette courbe-ci. Au lieu de calculer la pente de la droite en prenant un delta T=15°C, nous avons calculé la pente entre chaque point, puis nous avons fait la moyenne de toute ces droites, obtenant ainsi une valeur de alpha, égale à 1,12.10^-3. Pour l’incertitude, nous avons fait l’écart-type de toutes les valeurs de alpha égale à 7,3.10^-4. | 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 une certaine hauteur d’huile. Pour déterminer le coefficient de dilatation, nous connaissons cette équation : rho = rho0*(1-alphadelta T). On isole alors alpha. L’équation devient : alpha = Sdh(T-T0)/(Sdh+V0). Nous avons besoin de connaître la section du tube, la température et la hauteur à chaque instant, et le volume initial contenu dans le ballon. Pour calculer la pente, nous avons donc calculé la section du tube en verre et retrouver le volume initial par la mesure de la masse de l’huile. Pour connaître la température, nous avons utilisé un thermomètre mis dans le bain-marie. Et pour la hauteur, nous avions placé une règle graduée. L’expérience a été filmée, permettant ainsi de visionner et d’attribuer la hauteur du tube avec sa température correspondante. Nous avons tracé le graphique, et obtenu cette courbe-ci. Au lieu de calculer la pente de la droite en prenant un delta T=15°C, nous avons calculé la pente entre chaque point, puis nous avons fait la moyenne de toute ces droites, obtenant ainsi une valeur de alpha, égale à 1,12.10^-3. Pour l’incertitude, nous avons fait l’écart-type de toutes les valeurs de alpha égale à 7,3.10^-4. |
| 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. | 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. |
| | {{:wiki:projets:experience_1_graphe.png?400|}} {{ :wiki:projets:experience_1_montage.png?200|}} |
| |
| __2ème expérience :__ | __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'Aluminium dans une éprouvette graduée remplie d’huile de moteur. Nous avons calculé la masse volumique de la bille et de l’huile avec la relation classique ⍴=m/V en utilisant une balance et l’éprouvette pour l’huile, son diamètre et sa masse pour la bille. Nous avons rempli l’éprouvette graduée jusqu’à 500 mL, démarré la vidéo et lâché la bille à l’aide d’une pince le plus proche possible de la surface du fluide pour qu’elle parte avec une vitesse initiale la plus proche de 0 possible. Nous avons reproduit l’expérience plusieurs fois afin d’obtenir un grand nombre de valeurs ce qui nous a permis par la suite d’estimer l’incertitude sur notre mesure. à partir de ces vidéo nous avons écrit un programme sous Matlab, nous permettant de récupérer les positions de la bille dans l’éprouvette au cours du temps en traitant chaque vidéo image par image, cela fonctionne avec un système de pointeur, les deux premiers points déterminant la hauteur de fluide totale, les suivants servant à repérer la position de la bille à un intervalle de temps régulier ( 16 images = 16 positions pour un intervalle de temps de 1s soit 1 image toutes les 0,0625s). Finalement on obtient un graphe de la hauteur de la bille (position x verticale) en fonction du temps, on voit que la vitesse de la bille se stabilise et décrit une droite linéaire dont la pente est la vitesse dx/dt, la hauteur de notre éprouvette était donc assez grande pour que la vitesse de la bille soit constante au bout d’un certain temps. On a donc une valeur de la vitesse de chute de la bille on peut donc calculer la viscosité dynamique puis notre viscosité cinématique : nous avons trouvé ʋ=348 ±16 (1σ) mm^2/s, incertitude par écart-type (formule: ..) nous sommes dans le bon ordre de grandeur: v calculée théoriquement =500 | 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'Aluminium dans une éprouvette graduée remplie d’huile de moteur. Nous avons calculé la masse volumique de la bille et de l’huile avec la relation classique ⍴=m/V en utilisant une balance et l’éprouvette pour l’huile, son diamètre et sa masse pour la bille. Nous avons rempli l’éprouvette graduée jusqu’à 500 mL, démarré la vidéo et lâché la bille à l’aide d’une pince le plus proche possible de la surface du fluide pour qu’elle parte avec une vitesse initiale la plus proche de 0 possible. Nous avons reproduit l’expérience plusieurs fois afin d’obtenir un grand nombre de valeurs ce qui nous a permis par la suite d’estimer l’incertitude sur notre mesure. à partir de ces vidéo nous avons écrit un programme sous Matlab, nous permettant de récupérer les positions de la bille dans l’éprouvette au cours du temps en traitant chaque vidéo image par image, cela fonctionne avec un système de pointeur, les deux premiers points déterminant la hauteur de fluide totale, les suivants servant à repérer la position de la bille à un intervalle de temps régulier ( 16 images = 16 positions pour un intervalle de temps de 1s soit 1 image toutes les 0,0625s). Finalement on obtient un graphe de la hauteur de la bille (position x verticale) en fonction du temps, on voit que la vitesse de la bille se stabilise et décrit une droite linéaire dont la pente est la vitesse dx/dt, la hauteur de notre éprouvette était donc assez grande pour que la vitesse de la bille soit constante au bout d’un certain temps. On a donc une valeur de la vitesse de chute de la bille on peut donc calculer la viscosité dynamique puis notre viscosité cinématique : nous avons trouvé ʋ=348 ±16 (1σ) mm^2/s, incertitude par écart-type (formule: ..) nous sommes dans le bon ordre de grandeur: v calculée théoriquement =500\\ |
| | |
| | {{:wiki:projets:experience_2_graphe.png?400|}} {{ :wiki:projets:experience_2_montage.png?200|}} |
| |
| __3ème expérience :__ | __3ème expérience :__ |
| 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. | 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'huile de moteur ayant une hauteur H, celui-ci étant placé sur une plaque chauffante et contenant deux sondes. Celles-ci mesurent une température pendant une durée D et selon un pas de mesure nx. On a alors créé deux vecteurs, x, le vecteur distance, en fonction de dx, et t, le vecteur temps, en fonction de dt. Ensuite, on a formé une grille de valeurs en fonction de x et t. Puis, pour une gamme de kappa donnée, on a rempli un tableau de valeur de T en fonction de la distance x. Connaissant le Δx de notre expérience c’est à dire la distance entre nos deux sondes, il nous a été possible de déterminer les deux valeurs de T. Pour la valeur de la sonde située dans la partie basse, on ne peut pas prendre un x égal à o car l'augmentation de température durant le chauffage est trop importante à cause des effets de bords qui ne permettent pas une diffusion lente de la chaleur. Nous avons pris donc une valeur de x plus grande pour pouvoir respecter les conditions limites. Enfin, nous avons exprimé ΔT en fonction de t ce qui nous permet de comparer nos résultats avec des valeurs théoriques de log (K). Ainsi, par cette méthode nous avons pu déterminer le coefficient de diffusivité K. Celui-ci vaut 2x10-4 m²/s +ou- 1x10-4. Cette incertitude a été obtenue lors de l’encadrement de notre courbe (en noire sur le graphique) par deux valeurs de kappa (en jaune foncé et bleu). | Tout d’abord, nous avons modélisé un bécher rempli d'huile de moteur ayant une hauteur H, celui-ci étant placé sur une plaque chauffante et contenant deux sondes. Celles-ci mesurent une température pendant une durée D et selon un pas de mesure nx. On a alors créé deux vecteurs, x, le vecteur distance, en fonction de dx, et t, le vecteur temps, en fonction de dt. Ensuite, on a formé une grille de valeurs en fonction de x et t. Puis, pour une gamme de kappa donnée, on a rempli un tableau de valeur de T en fonction de la distance x. Connaissant le Δx de notre expérience c’est à dire la distance entre nos deux sondes, il nous a été possible de déterminer les deux valeurs de T. Pour la valeur de la sonde située dans la partie basse, on ne peut pas prendre un x égal à o car l'augmentation de température durant le chauffage est trop importante à cause des effets de bords qui ne permettent pas une diffusion lente de la chaleur. Nous avons pris donc une valeur de x plus grande pour pouvoir respecter les conditions limites. Enfin, nous avons exprimé ΔT en fonction de t ce qui nous permet de comparer nos résultats avec des valeurs théoriques de log (K). Ainsi, par cette méthode nous avons pu déterminer le coefficient de diffusivité K. Celui-ci vaut 2x10-4 m²/s +ou- 1x10-4. Cette incertitude a été obtenue lors de l’encadrement de notre courbe (en noire sur le graphique) par deux valeurs de kappa (en jaune foncé et bleu). |
| | {{:wiki:projets:experience_3_graphe.png?200|}} {{ :wiki:projets:experience_3_montage.png?200|}} |
| | |
| Nous sommes donc en possession de tous les paramètres nécessaires au calcul du nombre de Rayleigh critique. | 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 : | Calcul Ra expérimental et comparaison théorique : |
| Et surtout, pensez à mettre des photos. | Et surtout, pensez à mettre des photos. |
| </WRAP> | </WRAP> |
| | |
