Une investigation sur l'effet combiné de l'écoulement stationnaire et pulsatile dans un canal ondulé rempli d'un milieu poreux anisotrope

Informations fondamentales

• ID de l'article: 2510.08729
• Titre: An investigation into the combined effect of steady and pulsatile flow within a wavy channel filled with an anisotropic porous medium
• Auteurs: Sanchita Pramanik, Timir Karmakar (National Institute of Technology Meghalaya)
• Classification: physics.flu-dyn math-ph math.MP
• Date de publication: 13 octobre 2025
• Lien de l'article: https://arxiv.org/abs/2510.08729

Résumé

Cette étude fournit une solution analytique pour investiguer les effets de l'écoulement stationnaire et pulsatile dans un canal ondulé symétrique rempli d'un milieu poreux anisotrope. L'écoulement dans le canal est gouverné par l'équation de Darcy-Brinkman. L'étude emploie la méthode de perturbation pour résoudre le problème, en supposant que le rapport entre la largeur du canal et la longueur d'onde est très petit (δ² ≪ 1). Sous l'hypothèse λ²δ² ≪ 1, le problème est résolu à l'ordre O(δ²), où λ est le rapport d'anisotropie de la perméabilité. L'objectif principal est d'étudier comment le rapport d'anisotropie de la perméabilité affecte l'inversion d'écoulement près des crêtes et des creux des parois ondulées, incluant les conditions d'écoulement stationnaire et le cas où l'écoulement stationnaire est superposé à une composante pulsatile.

Contexte et motivation de la recherche

Contexte du problème

1. Signification scientifique: L'écoulement dans les canaux ondulés est un problème clé en mécanique des fluides, avec des applications larges dans le transfert de masse et de chaleur, l'ingénierie marine, la biomédicale et autres domaines
2. Applications d'ingénierie:
   • Oxygénation du sang dans les oxygénateurs à membrane
   • Simulation de l'écoulement sanguin dans le système cardiovasculaire (sténose artérielle et anévrisme)
   • Réservoirs géothermiques et systèmes de capture et stockage du carbone
   • Surfaces ondulées dans les appareils d'intensification du transfert de chaleur

Motivation de la recherche

1. Limitations existantes: Les recherches antérieures se sont principalement concentrées sur les milieux poreux isotropes, alors que les milieux poreux dans les applications réelles présentent souvent des caractéristiques anisotropes
2. Besoins technologiques: Les faisceaux de fibres dans les oxygénateurs à membrane forment des réseaux anisotropes, nécessitant une optimisation de la forme pour réaliser un échange gazeux efficace
3. Lacune théorique: L'effet combiné des parois ondulées et des pulsations de pression dans les milieux poreux anisotropes conduisant à la séparation d'écoulement manque d'études approfondies

Contributions principales

1. Innovation théorique: Première étude systématique de l'effet combiné de l'écoulement stationnaire et pulsatile dans un milieu poreux anisotrope
2. Méthode analytique: Fourniture d'une solution analytique à l'ordre O(δ²) basée sur la théorie de perturbation, applicable aux canaux avec petit rapport longueur-largeur
3. Mécanisme d'écoulement: Révélation du mécanisme d'influence du rapport d'anisotropie de la perméabilité sur l'inversion d'écoulement près des crêtes et des creux
4. Valeur applicative: Fourniture d'une base théorique pour la conception des oxygénateurs à membrane et le diagnostic des maladies cardiovasculaires
5. Analyse de la couche limite: Établissement de la théorie d'approximation de Darcy dans les régions éloignées de la paroi ondulée

Détails méthodologiques

Définition du problème

Étude des caractéristiques d'écoulement d'un fluide incompressible dans un canal ondulé bidimensionnel symétrique rempli d'un milieu poreux anisotrope, avec l'équation de la paroi du canal: $H^*(x^*) = b(1 + a\sin(2πx^*))$ où b est la largeur moyenne du canal, a est le paramètre d'amplitude, et L est la longueur d'onde.

Équations gouvernantes

L'écoulement est gouverné par l'équation de Darcy-Brinkman: $\rho_f \frac{\partial \mathbf{V}^*}{\partial t^*} = -\nabla P^* + \mu_{eff}\nabla^2\mathbf{V}^* - \mu\mathbf{K}^{-1}\mathbf{V}^*$

Tenseur de perméabilité anisotrope: $\mathbf{K} = \begin{pmatrix} K_x & 0 \\ 0 & K_y \end{pmatrix}$

Paramètres adimensionnels

• Rapport d'aspect: $δ = b/L$
• Nombre de Darcy: $Da = K_x/b^2$
• Rapport d'anisotropie: $λ^2 = K_x/K_y$
• Rapport de viscosité: $M = μ_{eff}/μ$
• Nombre de Womersley: $Wo = ρ_f Ω_p b^2/μ$

Développement en perturbation

Les vitesses et pressions sont développées selon δ²: $u_{(s,p)} = u_{(s0,p0)} + δ^2u_{(s1,p1)} + O(δ^4)$$v_{(s,p)} = v_{(s0,p0)} + δ^2v_{(s1,p1)} + O(δ^4)$$P_{(s,p)} = P_{(s0,p0)} + δ^2P_{(s1,p1)} + O(δ^4)$

Approximation de lubrification

Sous l'hypothèse δ² ≪ 1, la théorie de lubrification est employée pour résoudre les problèmes aux valeurs limites. Pour l'écoulement pulsatile, le paramètre ξ² = η² + α² est introduit, où η = √(iWo) et α² = 1/Da.

Configuration expérimentale

Plages de paramètres

• Rapport d'aspect: δ = 0,1-0,3
• Paramètre d'amplitude: a = 0,3
• Nombre de Darcy: Da = 0,001-1
• Rapport d'anisotropie: λ = 1-2
• Nombre de Womersley: Wo = 1-100
• Amplitude de pulsation: χ = 1

Conditions aux limites

1. Condition de non-glissement à la paroi: $u = 0, v = 0$ at $y = H(x)$
2. Condition de symétrie: $\partial u/\partial y = 0, v = 0$ at $y = 0$
3. Conservation du débit volumétrique: $\int_0^{H(x)} u(x,y)dy = 1$

Méthodes d'analyse

1. Écoulement stationnaire: Résolution de l'équation de Darcy-Brinkman stationnaire en fixant η = 0
2. Écoulement pulsatile: Considération de l'écoulement de Womersley pur entraîné par gradient de pression pulsatile
3. Effet combiné: Écoulement instationnaire avec superposition d'une composante pulsatile sur l'écoulement stationnaire

Résultats expérimentaux

Résultats d'écoulement stationnaire

1. Caractéristiques des lignes de courant:
   • Pour λ ≤ 1,5, les lignes de courant suivent essentiellement la forme de la paroi
   • Pour λ ≥ 1,75, une séparation et une zone de recirculation apparaissent près de la crête
   • La taille de la zone de recirculation augmente avec λ
2. Gradient de pression:
   • Pour les petites valeurs de λ, le gradient de pression est négatif (gradient favorable)
   • Avec l'augmentation de λ, un gradient de pression positif apparaît à la crête (gradient adverse)
   • Le gradient adverse provoque la séparation d'écoulement
3. Contrainte de cisaillement à la paroi:
   • Pour λ = 1,5, la contrainte de cisaillement est proche de zéro
   • Pour λ ≥ 1,75, la contrainte de cisaillement devient négative près de x = 0,25, indiquant une séparation d'écoulement

Résultats d'écoulement pulsatile

1. Évolution des tourbillons:
   • t = 0-Tp/4: Le fluide s'écoule dans la direction du flux principal
   • t = 2Tp/4: Formation de tourbillons près de la crête
   • t = 3Tp/4: Extension des tourbillons vers le flux principal, écoulement complètement inversé
2. Influence de l'anisotropie:
   • λ = 1 (isotrope): Aucune formation de tourbillon aux temps t = 0 et Tp/4
   • λ = 2 (anisotropie forte): Circulation évidente observable à tous les instants

Résultats d'effet combiné

1. Dépendance temporelle: Le moment et le lieu du début de l'inversion d'écoulement dépendent de l'interaction entre λ et le temps t
2. Variation du gradient de pression: Un gradient adverse apparaît aux temps t = 3Tp/4 et 4Tp/4
3. Épaisseur de la couche limite pulsatile: $δ_{puls} = 1/S$, où $S = \frac{1}{\sqrt{M}}(α^4 + Wo^2)^{1/4}\cos\left(\frac{1}{2}\tan^{-1}\left(\frac{Wo}{α^2}\right)\right)$

Analyse de la chute de pression

La chute de pression stationnaire augmente rapidement avec α² (l'inverse du nombre de Darcy), indiquant qu'une faible perméabilité augmente significativement la résistance à l'écoulement. L'augmentation de l'amplitude a entraîne également une augmentation de la chute de pression, reflétant l'effet de la rugosité de surface.

Travaux connexes

Recherches sur l'écoulement dans les canaux ondulés

• Tsangaris & Leiter (1984): Méthode de perturbation pour écoulement visqueux
• Pozrikidis (1987): Étude numérique de l'écoulement de Stokes
• Wang & Chen (2002): Transfert de chaleur par convection forcée

Écoulement dans les milieux poreux

• Ng & Wang (2010): Écoulement de Darcy-Brinkman
• Wei et al. (2003): Canal ondulé avec couche poreuse élastique
• Karmakar & Raja Sekhar (2017): Inversion d'écoulement dans un milieu poreux anisotrope

Applications d'écoulement pulsatile

• Womersley (1955): Base théorique de l'écoulement pulsatile artériel
• Sobey (1980-1983): Recherche sur l'écoulement dans les oxygénateurs à membrane
• Griffith et al. (2009): Écoulement pulsatile dans les géométries rétrécies

Conclusions et discussion

Conclusions principales

1. Effet d'anisotropie: L'augmentation de λ favorise significativement la séparation d'écoulement, formant des zones de recirculation près des crêtes et des creux
2. Évolution temporelle: La formation et la disparition des tourbillons dans l'écoulement pulsatile présentent une périodicité temporelle évidente
3. Mécanisme combiné: L'effet combiné du stationnaire et du pulsatile modifie le moment et le lieu de la séparation d'écoulement
4. Structure de la couche limite: Structure bicouche avec couche de Brinkman près de la paroi et région de Darcy en champ lointain

Applications d'ingénierie

1. Oxygénateurs à membrane: L'anisotropie et la pulsation peuvent optimiser l'efficacité du transfert d'oxygène
2. Diagnostic cardiovasculaire: Fourniture d'une base théorique pour l'analyse du flux sanguin dans la sténose artérielle et l'anévrisme
3. Intensification du transfert de chaleur: La séparation d'écoulement sur les surfaces ondulées intensifie le mélange par convection

Limitations

1. Hypothèse de petit paramètre: Applicable uniquement pour δ² ≪ 1 et λ²δ² ≪ 1
2. Simplification bidimensionnelle: Les effets tridimensionnels dans les applications réelles ne sont pas considérés
3. Linéarisation: La méthode de perturbation limite l'applicabilité aux cas de grande amplitude

Directions futures

1. Extension tridimensionnelle: Considération des canaux ondulés tridimensionnels avec variation transversale
2. Analyse non-linéaire: Effets non-linéaires dans les cas de grande amplitude
3. Validation expérimentale: Comparaison avec les données expérimentales
4. Conception optimale: Optimisation des paramètres géométriques basée sur le contrôle d'écoulement

Évaluation approfondie

Points forts

1. Rigueur théorique: Dérivation mathématique complète basée sur la théorie de perturbation
2. Intuition physique: Révélation du mécanisme profond de l'influence du rapport d'anisotropie de la perméabilité sur la séparation d'écoulement
3. Valeur applicative: Fourniture d'outils théoriques pratiques pour les oxygénateurs à membrane et l'écoulement cardiovasculaire
4. Systématicité: Analyse complète couvrant l'écoulement stationnaire, pulsatile et l'effet combiné
5. Théorie de la couche limite: Établissement de la structure bicouche Brinkman-Darcy

Insuffisances

1. Restrictions de paramètres: L'hypothèse de petit paramètre limite l'universalité de la méthode
2. Absence d'expériences: Manque de données expérimentales pour valider les prédictions théoriques
3. Comparaison numérique: Absence de comparaison avec les résultats de simulation CFD
4. Caractéristiques des matériaux: Discussion insuffisante des caractéristiques d'anisotropie des matériaux poreux réels

Impact

1. Contribution académique: Fourniture d'un nouveau cadre analytique pour la théorie de l'écoulement dans les milieux poreux
2. Valeur d'ingénierie: Perspectives d'application directe dans l'ingénierie biomédicale et l'intensification du transfert de chaleur
3. Innovation méthodologique: La méthode combinée perturbation-lubrification peut être généralisée à d'autres écoulements en géométries complexes

Scénarios applicables

1. Dispositifs microfluidiques: Contrôle d'écoulement dans les canaux ondulés à petite échelle
2. Biomédicale: Écoulement vasculaire et conception d'organes artificiels
3. Applications industrielles: Conception optimale des échangeurs de chaleur et appareils de séparation
4. Ingénierie géologique: Analyse du transport de fluides dans les couches rocheuses poreuses

Références bibliographiques

L'article cite 62 références pertinentes, couvrant les domaines importants de la mécanique des fluides, de la théorie des milieux poreux et de l'ingénierie biomédicale, fournissant une base théorique solide pour cette recherche.