L’optimisation de la production d’une huile essentielle passe par la modélisation cinétique de son extraction. Le présent travail se propose de déterminer graphiquement des paramètres macroscopiques des modèles empiriques utilisables dans la modélisation cinétique de l’extraction de l’huile essentielle de C. nardus. Les huiles essentielles étudiées ont été extraites par vapo-hydrodistillation et ont été analysées par chromatographie en phase gazeuse couplée à la spectrométrie de masse. La modélisation a été conduite selon l’approche phénoménologique qui distingue deux étapes : une étape rapide caractérisée par la constante de vitesse k1 et une deuxième étape de diffusion de soluté dans la matrice végétale avant son élution (lavage) par le solvant, caractérisée par k2. L’étude a été menée au laboratoire et en atelier pilote sur quatre modèles empiriques : Sovova-Milojevic, Peleg, Monod et Langmuir. Au laboratoire, les résultats expérimentaux valident le modèle du pseudo premier de Sovova-Milojevic (R2 = 0,981 – 0,998), Peleg (R2= 0,984 – 0.993), Langmuir (R2 = 0,996) et Monod (R2 = 0,996), à l’échelle pilote les modèles de Peleg (R2= 0,9853) Monod et Langmuir (R2 = 0,989). La vapo-hydrodistillation de C. nardus se passe en deux étapes (lavage et diffusion) au laboratoire et en atelier pilote. Au laboratoire les résultats expérimentaux valident les quatre modèles testés ; à l’échelle pilote le modèle de Sovova-Milojevic n’est pas validé. Ce constat suggère des mécanismes légèrement différents pour les deux échelles, différence probablement due à l’effet d’échelle.