Ce mémoire de maîtrise en milieu industriel présente les travaux réalisés chez Maya Heat Transfer Technologies afin d'implémenter la technique de la simulation des grandes échelles dans le solveur de dynamique des fluides parallèle NIECE. Les modèles de sous-maille [sic] algébriques de Smagorinsky, WALE (Wale Adapting Local Eddy-Viscosity) et le modèle de Vreman ont été retenus et validés.Ce projet de maîtrise se divise [entre] quatre volets. La première partie présente l'implémentation et l'analyse d'un nouveau schéma de discrétisation spatial et temporel dans le contexte d'un algorithme de résolution des équations de Navier-Stokes par une approche mixte [de] volumes éléments finis sur maillage non structuré. Une étude comparative sur la stabilité, la diffusion et [la] dipersion [sic] numérique des nouveaux schémas et ceux déjà existants dans le code permet d'identifier le schéma centré de second ordre en espace, couplé au schéma de Crank-Nicolson en temps comme un compromis adéquat pour la simulation des grandes échelles dans un contexte industriel. Le second volet détaille les résulats [i.e. résultats] de la simulation de la décroissance d'une turbulence homogène isotrope. La décroissance des spectres d'énergie est présentée pour chacun des modèles de sous-maille, puis comparée à la solution sans modèle. L'effet du raffinement du maillage est investigué.Ce volet détaille également l'implémentation et la validation de l'algorithme permettant l'initialisation d'une turbulence isotrope incompressible. Une discussion sur les formulations analytiques des spectres d'énergie de Passot-Pouquet et de von Karman-Pao est finalement présentée. Le problème de la répresentation [i.e. représentation] de l'écoulement à l'entrée du domaine de calcul est ensuite abordé.Ce troisième volet détaille la méthode de génération de turbulence, appelée méthode des tourbillons synthétiques, servant à générer un champ de vitesse cohérent et représentatif d'une turbulence réelle. Une étude de l'influence des paramètres numériques de la méthode permet de dériver un critère automatique de sélection pour l'utilisateur. Pour terminer, une étude comparative de la méthode des tourbillons synthétiques à une méthode d'injection aléatoire est présentée. En dernier lieu, la réalisation du cas de validation du canal plan périodique permet d'investiguer la capacité du code à simuler les propriétés d'une turbulence cisaillée. Les résultats des modèles de sous-maille sont comparés au cas sans modèle de turbulence. L'effet du raffinement du maillage est étudié sur trois différentes grilles. Les statistiques de premier ordre (profil de vitesse moyenne) et les satistiques [i.e. statistiques] de second ordre (profils des composantes du tenseur de Reynolds) sont détaillées et analysées. Finalement, la contribution de l'énergie cinétique résolue et de l'énergie de sous-maille permet de quantifier la contribution des modèles de turbulence sur la solution.

Le développement des hydroliennes se base souvent sur des conditions d'écoulement idéalisées qui ne reflètent pas entièrement le courant présent sur un site réel. Dans cette optique, ce mémoire investigue l'effet d'un écoulement non aligné ou cisaillé sur les performances et le chargement de deux types d'hydrolienne : les hydroliennes à rotor axial (HRA) et à aile oscillante (HAO). Cette étude est réalisée à l'aide de simulations numériques et les résultats démontrent que pour les deux types d'hydrolienne, le désalignement produit une réduction de la puissance et de la traînée, tandis que le cisaillement n'a qu'un effet de second ordre sur ces derniers. Le chargement additionnel sur la structure de support causé par des conditions d'écoulement non idéal est aussi similaire pour les deux technologies. Le désalignement et le cisaillement affectent cependant plus significativement le chargement en fatigue des pales de l'HRA que celui de l'aile de l'HAO. / Le développement des hydroliennes se base souvent sur des conditions d'écoulement idéalisées qui ne reflètent pas entièrement le courant présent sur un site réel. Dans cette optique, ce mémoire investigue l'effet d'un écoulement non aligné ou cisaillé sur les performances et le chargement de deux types d'hydrolienne : les hydroliennes à rotor axial (HRA) et à aile oscillante (HAO). Cette étude est réalisée à l'aide de simulations numériques et les résultats démontrent que pour les deux types d'hydrolienne, le désalignement produit une réduction de la puissance et de la traînée, tandis que le cisaillement n'a qu'un effet de second ordre sur ces derniers. Le chargement additionnel sur la structure de support causé par des conditions d'écoulement non idéal est aussi similaire pour les deux technologies. Le désalignement et le cisaillement affectent cependant plus significativement le chargement en fatigue des pales de l'HRA que celui de l'aile de l'HAO. / The development of hydrokinetic turbines is often based on idealized flow conditions which do not fully refect river or tidal currents. In this regard, this thesis investigates the effect of non-aligned or sheared flows on the performances and loading of two turbine types: the axial rotor (ART) and oscillating foil (OFT) turbines. This study was conducted with unsteady numerical simulations and the obtained results show that for both turbine types, misalignment produces a reduction in power and thrust, while the shear has a limited effect on those quantities. The additional loading on the support structure caused by the non-idealized flow is also similar for both devices. However, misalignment and shear affect more severely the fatigue loading on ART blades than they do for OFT. / The development of hydrokinetic turbines is often based on idealized flow conditions which do not fully refect river or tidal currents. In this regard, this thesis investigates the effect of non-aligned or sheared flows on the performances and loading of two turbine types: the axial rotor (ART) and oscillating foil (OFT) turbines. This study was conducted with unsteady numerical simulations and the obtained results show that for both turbine types, misalignment produces a reduction in power and thrust, while the shear has a limited effect on those quantities. The additional loading on the support structure caused by the non-idealized flow is also similar for both devices. However, misalignment and shear affect more severely the fatigue loading on ART blades than they do for OFT.

Thèse ou mémoire avec insertion d'articles / Titre de l'écran-titre (visionné le 11 janvier 2024) / Cette thèse de doctorat étudie le transport de fluides viscoplastiques à travers des canaux dotés de parois superhydrophobes striées, en utilisant une approche de modélisation complète et des simulations numériques. L'objectif est de comprendre les interactions entre le fluide viscoplastique et la surface superhydrophobe, et d'identifier les facteurs qui influencent le comportement de l'écoulement. On suppose que les rainures emprisonnent des poches d'air à la surface superhydrophobe, maintenant l'interface liquide/air plate et fixée aux bords des rainures, ce qui entraîne une condition de glissement-collage. L'analyse est complète, couvrant à la fois les régimes de fluage et inertiel dans les limites des canaux épais et minces. En particulier, six paramètres, notamment le nombre de Reynolds (R), le nombre de Bingham (B), le nombre de glissement (b), la longueur de périodicité des rainures (ℓ), la fraction de zone de glissement (φ), et l'angle d'orientation des rainures (θ), sont pris en compte. Ensuite, les effets de ces paramètres sont analysés sur les variables d'écoulement d'intérêt, telles que le champ de perturbation, le champ de vitesse, la magnitude du taux de déformation, les zones de bouchon non cédées, la longueur de glissement effective, le facteur de frottement, l'asymétrie de l'écoulement, le mélange, et les classifications de régime. Dans les canaux épais (ℓ ≪ 1) avec des rainures transversales (θ = 90°), l'augmentation de B et de b entraîne une croissance des champs de perturbation et de glissement, qui deviennent asymétriques en augmentant R. Dans certaines conditions d'écoulement, le bouchon central se rompt ou une zone de bouchon non cédée peut se former à l'interface liquide/air de la paroi superhydrophobe, ce qui nous permet de classifier les régimes d'écoulement. Dans les canaux minces (ℓ ≫ 1), l'augmentation de b ou de θ entraîne la déformation du bouchon central non cédé, sauf pour les écoulements avec des rainures longitudinales (θ = 0). Pour les rainures obliques (0 < θ < 90°), un écoulement secondaire apparaît, déclenchant le mélange de l'écoulement. Enfin, l'analyse de stabilité linéaire de la condition de glissement homogène (φ = 1) révèle les effets stabilisants/déstabilisants des conditions de glissement dans le sens du courant/transversales. / This Ph.D. thesis investigates the transport of viscoplastic fluids through channels with grooved superhydrophobic (SH) walls, using a comprehensive modeling approach and numerical simulations. The goal is to understand the interactions between the viscoplastic fluid and the superhydrophobic surface and identify factors that influence the flow behaviour. It is assumed that the grooves trap air pockets on the SH surface, keeping the liquid/air interface flat and pinned at the groove edges, resulting in a slip-stick condition. The analysis is comprehensive, covering both creeping and inertial regimes in thick and thin channel limits. In particular, six parameters, including the Reynolds number (R), Bingham number (B), slip number (b), groove periodicity length (ℓ), slip area fraction (φ), and groove orientation angle (θ), are considered. Subsequently, the effects of these parameters are analyzed on the flow variables of interest, such as the perturbation field, velocity field, strain-rate magnitude, unyielded plug zones, effective slip length, friction factor, flow asymmetry, mixing, and regime classifications. In thick channels (ℓ ≪ 1) with transverse grooves (θ = 90°), increasing B and b leads to growing the perturbation and slip velocity fields, which become asymmetric by increasing R. Under certain flow conditions, the center plug breaks or an unyielded plug zone can form on the SH wall liquid/air interface, allowing us to classify the flow regimes. In thin channels (ℓ ≫ 1), increasing either b or θ causes the the unyielded center plug to deform and eventually break, except for flows with longitudinal grooves (θ = 0). For the oblique grooves (0 < θ < 90°), a secondary flow appears, triggering flow mixing. Finally, the linear stability analysis of the homogeneous slip condition (φ = 1) reveals stabilizing/destabilizing effects of the streamwise/spanwise slip conditions.

Tableau d'honneur de la Faculté des études supérieures et postdoctorales, 2011-2012 / Ce mémoire de maîtrise porte sur la modélisation numérique de la problématique d'ingestion de gaz chauds dans la cavité rotor-stator du premier étage d'une turbine à gaz moderne. Dans le contexte de développement d'une aviation plus écologique, ce sujet mérite une attention marquée car des gains importants en ce qui a trait à l'efficacité globale du moteur peuvent être atteints par de meilleurs designs de joints rotor-stator limitant l'ingestion. La modélisation physique des phénomènes affectant l'ingestion devient alors essentielle afin de développer une meilleure compréhension des écoulements et mécanismes en présence. C'est dans cette optique que Pratt & Whitney Canada a mandaté le Laboratoire de Mécanique des Fluides Numérique (LMFN) de l'Université Laval afin de conduire une étude approfondie du phénomène d'ingestion. Ce mémoire a donc pour objectif principal de décrire le développement, l'application et la validation de certaines techniques de simulation numérique visant à caractériser qualitativement et quantitativement le phénomène d'ingestion observé à travers le joint et dans la cavité du groupe rotor-stator d'une turbine à gaz. Une investigation des limites de fiabilité d'une telle modélisation numérique dans un contexte réel d'opération du moteur est également présentée. La paramétrisation adéquate du problème et la validité d'extrapoler des résultats numériques d'ingestion obtenus aux conditions "banc d'essai" afin de prédire l'ingestion réelle aux conditions "moteur" font également partie des sujets abordés dans ce mémoire. La méthodologie développée considère deux phases de calculs menant à des diagnostics d'ingestion par le biais d'un marqueur passif. Une série de validations numériques et des comparaisons avec des résultats expérimentaux sont également présentées. La méthodologie de modélisation numérique est considérée robuste et fiable, mais seulement partiellement validée.

No description available.

Les suspensivores ont la tâche importante de séparer les particules de l'eau. Bien qu'une grande gamme de morphologies existe pour les structures d'alimentation, elles sont pratiquement toutes constituées de rangées de cylindres qui interagissent avec leur environnement fluide. Le mécanisme de capture des particules utilisé dépend des contraintes morphologiques, des besoins énergétiques et des conditions d'écoulement. Comme nos objectifs étaient de comprendre ces relations, nous avons eu recours à des études de comparaison pour interpréter les tendances en nature et pour comprendre les conditions qui provoquent de nouveaux fonctionnements. Nous avons utilisé la dynamique des fluides numérique (computational fluid dynamics, CFD) pour créer des expériences contrôlées et pour simplifier les analyses. Notre première étude démontre que les coûts énergétiques associés au pompage dans les espaces petits sont élevés. De plus, le CFD suggère que les fentes branchiales des ptérobranches sont des structures rudimentaires, d'un ancêtre plus grande. Ce dernier point confirme l'hypothèse qu'un ver se nourrit par filtration tel que l'ancêtre des deuterostomes. Notre deuxième étude détermine la gamme du nombre de Reynolds number critique où la performance d'un filtre de balane change. Quand le Re est très bas, les différences morphologiques n'ont pas un grand effet sur le fonctionnement. Cependant, une pagaie devient une passoire lorsque le Re se trouve entre 1 et 3,5. Le CFD s’est dévoilé être un outil très utile qui a permis d’obtenir des détails sur les microfluides. Ces études montrent comment la morphologie et les dynamiques des fluides interagissent avec la mécanisme de capture ou de structures utilisées, ainsi que comment des petits changements de taille, de forme, ou de vitesse d'écoulement peuvent conduire à un nouveau fonctionnement. / Suspension feeders have the important task of separating particles from the liquid medium surrounding them. Although a wide range of morphologies exist for feeding structures, essentially all consist of arrays of cylinders interacting with their fluid environment. The particle capture mechanism employed depends on morphological constraints, energy requirements, and flow conditions, therefore our objectives were to understand these relationships through comparison studies to interpret trends in nature and to understand the conditions in which novel functioning arises. Our studies used computational fluid dynamics (CFD) to create controlled experiments and to simplify analyses. The first study demonstrates the unfeasibly high energetic costs of pumping in small and medium pharynx sizes of hemichordates, implying the gill pores of pterobranchs are likely vestigial structures from a larger ancestor. This last point further supports the hypothesis of an acorn worm as the ancestor to all deuterostomes. Our second study determined the critical Reynolds number (Re) range where barnacle filters transition from a paddle to a sieve. At very low Re, morphological differences have no major effect on functioning and filters behave as paddles, however, at Re 1 - 3.5, these differences dictate when a paddle becomes a sieve. CFD proved to be a very useful tool for simplifying studies and providing detailed microfluidics. These studies demonstrate how morphology and fluid dynamics interact to dictate the capture mechanisms or appendages employed and how simple changes in size, shape, or flow speed can lead to novel functioning.

Dans cette dissertation, nous considérons l’écoulement des liquides conducteurs d’électricité dans un champ magnétique externe. De tels écoulements sont décrits par les équations de la magnétohydrodynamique (MHD) quasi-statique, et sont fréquemment rencontrés dans des applications pratiques. Il suit qu’il y a un intérêt fort pour des outils numérques qui peuvent simuler ces écoulements dans des géometries complexes.<p>La première partie de cette thèse (chapitres 2 et 3) est dédiée à la présentation de la machinerie numérique qui a été utilisée et implémentée afin de résoudre les équations de la MHD quasi-statistique (incompressible). Plus précisément, nous avons contribué au développement d’un solveur volumes finis non-structuré parallèle. La discussion sur ces méthodes est accompagnée d’une analyse numérique qui est aussi valable pour des mailles non-structurées. Dans le chapitre 3, nous vérifions notre implémentation par la simulation d’un certain nombre de cas tests avec un accent sur des écoulements dans un champ magnétique intense.<p>Dans la deuxième partie de cette thèse (chapitres 4-6), nous avons utilsé ce solveur pour étudier des écoulements MHD de proche paroi .La première géometrie considérée (chapitre 4) est celle d’une conduite circulaire infini d’axe à haut nombre de Hartmann. Nous avons investitgué la sensitivité des résultats numériques au schéma de discrétisation et à la topologie de la maille. Nos résultats permettent de caractériser in extenso l’écoulement MHD dans une conduite avec des bords bien conducteurs par moyen des lois d’échelle.<p>Le sujet du cinquième chapitre est l’écoulement dans une conduite toroïdale à section carée. Une étude du régime laminaire confirme une analyse asymptotique pour ce qui concerne les couches de cisaillement. Nous avons aussi effectué des simulations des écoulements turbulents afin d’évaluer l’effet d’un champ magnétique externe sur l’état des couches limites limites.<p>Finalement, dans le chapitre 6, nous investiguons l’écoulement MHD et dans un U-bend et dans un coude arrière. Nous expliquons comment générer une maille qui permet de toutes les couches de cisaillement à un coût computationelle acceptable. Nous comparons nos résultats aux solutions asymptotiques. / Doctorat en Sciences / info:eu-repo/semantics/nonPublished

Les suspensivores ont la tâche importante de séparer les particules de l'eau. Bien qu'une grande gamme de morphologies existe pour les structures d'alimentation, elles sont pratiquement toutes constituées de rangées de cylindres qui interagissent avec leur environnement fluide. Le mécanisme de capture des particules utilisé dépend des contraintes morphologiques, des besoins énergétiques et des conditions d'écoulement. Comme nos objectifs étaient de comprendre ces relations, nous avons eu recours à des études de comparaison pour interpréter les tendances en nature et pour comprendre les conditions qui provoquent de nouveaux fonctionnements. Nous avons utilisé la dynamique des fluides numérique (computational fluid dynamics, CFD) pour créer des expériences contrôlées et pour simplifier les analyses. Notre première étude démontre que les coûts énergétiques associés au pompage dans les espaces petits sont élevés. De plus, le CFD suggère que les fentes branchiales des ptérobranches sont des structures rudimentaires, d'un ancêtre plus grande. Ce dernier point confirme l'hypothèse qu'un ver se nourrit par filtration tel que l'ancêtre des deuterostomes. Notre deuxième étude détermine la gamme du nombre de Reynolds number critique où la performance d'un filtre de balane change. Quand le Re est très bas, les différences morphologiques n'ont pas un grand effet sur le fonctionnement. Cependant, une pagaie devient une passoire lorsque le Re se trouve entre 1 et 3,5. Le CFD s’est dévoilé être un outil très utile qui a permis d’obtenir des détails sur les microfluides. Ces études montrent comment la morphologie et les dynamiques des fluides interagissent avec la mécanisme de capture ou de structures utilisées, ainsi que comment des petits changements de taille, de forme, ou de vitesse d'écoulement peuvent conduire à un nouveau fonctionnement. / Suspension feeders have the important task of separating particles from the liquid medium surrounding them. Although a wide range of morphologies exist for feeding structures, essentially all consist of arrays of cylinders interacting with their fluid environment. The particle capture mechanism employed depends on morphological constraints, energy requirements, and flow conditions, therefore our objectives were to understand these relationships through comparison studies to interpret trends in nature and to understand the conditions in which novel functioning arises. Our studies used computational fluid dynamics (CFD) to create controlled experiments and to simplify analyses. The first study demonstrates the unfeasibly high energetic costs of pumping in small and medium pharynx sizes of hemichordates, implying the gill pores of pterobranchs are likely vestigial structures from a larger ancestor. This last point further supports the hypothesis of an acorn worm as the ancestor to all deuterostomes. Our second study determined the critical Reynolds number (Re) range where barnacle filters transition from a paddle to a sieve. At very low Re, morphological differences have no major effect on functioning and filters behave as paddles, however, at Re 1 - 3.5, these differences dictate when a paddle becomes a sieve. CFD proved to be a very useful tool for simplifying studies and providing detailed microfluidics. These studies demonstrate how morphology and fluid dynamics interact to dictate the capture mechanisms or appendages employed and how simple changes in size, shape, or flow speed can lead to novel functioning.

Cette étude développe une plate-forme mathématique allant de la création d’une base de connaissances à la mise en place d’un modèle adaptatif. Une nouvelle approche en ce qui concerne la modélisation du four rotatif du clinker a été nécessaire pour mener à bien notre étude. Ainsi, un modèle d’état du système aux paramètres repartis, à la base des phénomènes physico-chimiques, a été conçu à l’aide d’équations aux dérivées partielles. La structure du modèle est basée sur trois variables d’état qui peuvent décrire correctement les principaux phénomènes (diffusions, déplacements, transferts et échanges thermiques, etc.) qui se déroulent dans le four rotatif de clinker. En outre, ces variables d’état sont accessibles à travers des paramètres mesurables tels que les températures de la surface extérieure du four et de l’air ambiant. Les paramètres du modèle proposé sont des fonctions de trois variables d’état, qui tout en remplaçant les paramètres mal connus et inaccessibles, décrivent les phénomènes physicochimiques qui caractérisent la complexité du système et sont définis comme des fonctions opératoires. Cette approche permet non seulement de minimiser les erreurs (hypothèse de simplification) mais aussi permet un accès facile et aisé des paramètres. Nous avons élaboré une procédure d’identification basée sur des analyses phénoménologique et dimensionnelle. Ainsi, à partir de données issues d’un complément d’instrumentation associé à la base de connaissances, l’identification des fonctions opératoires (paramètres) a été effectuée à partir d’un état stationnaire du système qui représente notre condition initiale. Cette opération est importante dans la mesure où elle donne aux paramètres du modèle les propriétés de la dynamique du système. Ces résultats sont corroborés par les estimations des variables d’état dont les valeurs et évolutions sont acceptables. Après l’évaluation de l’influence des fonctions opératoires sur les variables d’état, une modélisation des paramètres (fonctions opératoires) du modèle a été proposée à la base d’une décomposition physique à la lumière des connaissances fournies dans la littérature. Une procédure d’ajustement a été également développée afin de pouvoir fournir des variables d’état optimales. Ainsi, l’erreur peut être compensée par ajustement des fonctions opératoires. Par conséquent, à travers la rapidité de la procédure d’ajustement et l’adaptation des paramètres du modèle dans un horizon de temps bien déterminé, le modèle renferme les caractéristiques d’un modèle adaptatif qui pourra être un support essentiel dans la commande automatique du système. / This study develops a mathematical platform going from the establishment of a knowledge database to the setting up of an adaptive model. This has required a new approach of modeling of the clinker rotary kiln (CRK). Thus, a state model of distributed parameter systems, based on physico-chemical phenomena, was designed using partial differential equations. The model structure is based on three state variables which are: the gas, clinker temperatures and the clinker mass distributions, and are elaborated with the help of heat, pressure and mass balance equations. The model parameters are defined by the functions of three state variables. Moreover, the resulting state model, decomposed into five phenomenological zones of CRK, is used as a first step to define a set of Operating Functions (OFs). These OFs has also been decomposed into longitudinal distribution of CRK to replace the constant, unknown or unmeasured parameters. We develop an identification procedure based on phenomenological and dimensional analysis where the identification of operational functions (model parameters) was performed from a stationary state of the CRK. Once the restores state variables have been evaluated, the desired input (which is treated as the control of the CRK) can be more easily found by the proposed model than by simple trial and error. Moreover, the fact that the computation time, to estimate-calibrate the OFs above-mentioned, is very short, then this dynamic computation works faster than real-time. In summary, the cooperation and coordination in real-time between industrial computers and the CRK allows for an adaptable model, where each specific set of the OFs must be analyzed by its accuracy.

Le monoxyde de carbone déshydrogénase (CODH) catalyse la réduction réversible du dioxyde de carbone par son site actif. En utilisant une méthode électrochimique appelée voltammétrie de film protéique, nous étudions le mécanisme enzymatique de CODH en immobilisant l'enzyme à une surface d'électrode de graphite où le transfert direct d'électrons est possible. Traditionnellement, pour empêcher la déplétion du substrat à l'électrode, les électrochimistes utilisent des électrodes tournantes (RDE). Cependant, comme la CODH est très active, même la RDE ne peut pas empêcher l'épuisement, ce qui masque les caractéristiques cinétiques importantes de l’enzyme. Nous ne pouvons pas résoudre le problème avec RDE, puisque nous l’utilisons déjà à la vitesse maximum. Par conséquent, nous devons concevoir une nouvelle cellule électrochimique. Pour cela, nous avons utilisé des simulations de dynamique des fluides computationnelles pour explorer diverses géométries afin d'en trouver une appropriée. Nous avons commencé par valider notre méthode numérique avec la solution théorique bien définie de la cellule réelle de RDE. Après la bonne validation, nous avons déterminé les vitesses de transport de masse au sein de plusieurs géométries et à basé sur l'optimisation des paramètres géométriques, nous avons atteint notre conception appropriée. Ce nouveau prototype a une électrode graphite uniformément accessible avec un taux de transport trois fois plus rapide que le RDE à des vitesses de solution acceptables. Nous avons construit, mis en place avec succès le système pour caractériser ses performances de transport, et trouvé un excellent accord entre les résultats numériques et expérimentaux / Carbon monoxide dehydrogenase (CODH) catalyzes the reversible reduction of carbon dioxide by its active site. Thus, CODH participates in the first step of fuel production. Using an electrochemical method called protein film voltammetry, we study the enzymatic mechanism of CODH by immobilizing the enzyme at a graphite electrode surface where direct electron transfer is possible. Traditionally, to prevent depletion of the substrate at the electrode, electrochemists use rotating electrodes (RDE). However, since CODH is very active, even RDE cannot prevent depletion, which masks the important kinetic characteristics of the enzyme and complicates the analysis of the enzymatic response.We cannot solve the problem with RDE, since we already use it at maximum speed. Therefore, we must completely change our approach and design a new electrochemical cell. For this, we used computational fluid dynamics (CFD) simulations to explore various geometries to find a suitable one. We began by validating our numerical method with the well-defined theoretical solution of the real cell of RDE. After good validation, we determined the mass transport velocities within several proposed geometries of the flow cell of hydrodynamic channel and jet electrodes. Based on the optimization of geometric parameters, we have achieved our proper design of jet electrode. This new prototype has a uniformly accessible graphite electrode with a transport rate three times faster than the RDE at acceptable solution speeds. We have successfully built and implemented the system to characterize its transport performance. We found an excellent agreement between the numerical and experimental results