Una investigación sobre el efecto combinado del flujo estacionario y pulsátil dentro de un canal ondulado lleno de un medio poroso anisotrópico

Información Básica

• ID del Artículo: 2510.08729
• Título: Una investigación sobre el efecto combinado del flujo estacionario y pulsátil dentro de un canal ondulado lleno de un medio poroso anisotrópico
• Autores: Sanchita Pramanik, Timir Karmakar (Instituto Nacional de Tecnología Meghalaya)
• Clasificación: physics.flu-dyn math-ph math.MP
• Fecha de Publicación: 13 de octubre de 2025
• Enlace del Artículo: https://arxiv.org/abs/2510.08729

Resumen

Este estudio proporciona una solución analítica para investigar los efectos del flujo estacionario y pulsátil en un canal ondulado simétrico lleno de un medio poroso anisotrópico. El flujo dentro del canal se rige por la ecuación de Darcy-Brinkman. El estudio emplea el método de perturbación para resolver el problema, asumiendo que la relación entre el ancho del canal y la longitud de onda es pequeña (δ² ≪ 1). Bajo la suposición de λ²δ² ≪ 1, el problema se resuelve hasta el orden O(δ²), donde λ es la relación de permeabilidad anisotrópica. El objetivo principal es investigar cómo la relación de permeabilidad anisotrópica afecta la inversión del flujo cerca de las crestas y valles de las paredes onduladas, incluyendo condiciones de flujo estacionario y casos donde se superpone una componente pulsátil al flujo estacionario.

Antecedentes de Investigación y Motivación

Contexto del Problema

1. Significancia Científica: El flujo en canales ondulados es un problema clave en mecánica de fluidos, con aplicaciones amplias en transferencia de masa y calor, ingeniería marina, bioingeniería y otros campos
2. Aplicaciones de Ingeniería:
   • Oxigenación de sangre en oxigenadores de membrana
   • Simulación de flujo sanguíneo en el sistema cardiovascular (estenosis arterial y aneurismas)
   • Depósitos geotérmicos y sistemas de captura y almacenamiento de carbono
   • Superficies onduladas en dispositivos de intensificación de transferencia de calor

Motivación de la Investigación

1. Limitaciones Existentes: Investigaciones previas se han enfocado principalmente en medios porosos isotrópicos, mientras que los medios porosos en aplicaciones reales frecuentemente exhiben características anisotrópicas
2. Necesidades Técnicas: Los haces de fibras en oxigenadores de membrana forman redes anisotrópicas, requiriendo optimización de forma para lograr intercambio de gas eficiente
3. Brecha Teórica: Falta investigación profunda sobre fenómenos de separación de flujo causados por los efectos combinados de paredes onduladas y pulsaciones de presión en medios porosos anisotrópicos

Contribuciones Principales

1. Innovación Teórica: Primera investigación sistemática del efecto combinado del flujo estacionario y pulsátil en medios porosos anisotrópicos
2. Método Analítico: Proporciona solución analítica de orden O(δ²) basada en teoría de perturbación, aplicable a canales con relación de aspecto pequeña
3. Mecanismo de Flujo: Revela el mecanismo de cómo la relación de permeabilidad anisotrópica afecta la inversión de flujo cerca de crestas y valles
4. Valor de Aplicación: Proporciona base teórica para el diseño de oxigenadores de membrana y diagnóstico de enfermedades cardiovasculares
5. Análisis de Capa Límite: Establece teoría de aproximación de Darcy en regiones alejadas de la pared ondulada

Detalles Metodológicos

Definición del Problema

Se investigan las características de flujo de fluido incompresible en un canal ondulado bidimensional simétrico lleno de medio poroso anisotrópico, con ecuación de pared: $H^*(x^*) = b(1 + a\sin(2πx^*))$ donde b es el ancho promedio del canal, a es el parámetro de amplitud, y L es la longitud de onda.

Ecuaciones Gobernantes

El flujo se rige por la ecuación 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}^*$

Tensor de permeabilidad anisotrópica: $\mathbf{K} = \begin{pmatrix} K_x & 0 \\ 0 & K_y \end{pmatrix}$

Parámetros Adimensionales

• Relación de aspecto: $δ = b/L$
• Número de Darcy: $Da = K_x/b^2$
• Relación de anisotropía: $λ^2 = K_x/K_y$
• Relación de viscosidad: $M = μ_{eff}/μ$
• Número de Womersley: $Wo = ρ_f Ω_p b^2/μ$

Expansión de Perturbación

Velocidad y presión se expanden en δ²: $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)$

Aproximación Lubricante

Bajo la suposición δ² ≪ 1, se emplea teoría lubricante para resolver problemas de valores en la frontera. Para flujo pulsátil, se introduce el parámetro ξ² = η² + α², donde η = √(iWo) y α² = 1/Da.

Configuración Experimental

Rangos de Parámetros

• Relación de aspecto: δ = 0.1-0.3
• Parámetro de amplitud: a = 0.3
• Número de Darcy: Da = 0.001-1
• Relación de anisotropía: λ = 1-2
• Número de Womersley: Wo = 1-100
• Amplitud de pulsación: χ = 1

Condiciones de Frontera

1. Condición de no deslizamiento en la pared: $u = 0, v = 0$ en $y = H(x)$
2. Condición de simetría: $\partial u/\partial y = 0, v = 0$ en $y = 0$
3. Conservación de flujo volumétrico: $\int_0^{H(x)} u(x,y)dy = 1$

Métodos de Análisis

1. Flujo Estacionario: Se establece η = 0 para resolver la ecuación de Darcy-Brinkman estacionaria
2. Flujo Pulsátil: Se considera flujo de Womersley impulsado por presión pulsátil pura
3. Efecto Combinado: Flujo no estacionario con componente pulsátil superpuesta al flujo estacionario

Resultados Experimentales

Resultados de Flujo Estacionario

1. Características de Líneas de Corriente:
   • Cuando λ ≤ 1.5, las líneas de corriente se distribuyen básicamente siguiendo la forma de la pared
   • Cuando λ ≥ 1.75, aparecen regiones de separación y reflujo cerca de las crestas
   • El tamaño de la región de reflujo se expande con el aumento de λ
2. Gradiente de Presión:
   • Con valores pequeños de λ, el gradiente de presión es negativo (gradiente favorable)
   • Cuando λ aumenta, aparece un gradiente de presión positivo en la cresta (gradiente adverso)
   • El gradiente adverso causa separación del flujo
3. Esfuerzo Cortante en la Pared:
   • Cuando λ = 1.5, el esfuerzo cortante es cercano a cero
   • Cuando λ ≥ 1.75, el esfuerzo cortante se vuelve negativo cerca de x = 0.25, indicando separación del flujo

Resultados de Flujo Pulsátil

1. Evolución de Vórtices:
   • t = 0-Tp/4: El fluido se mueve en la dirección del flujo principal
   • t = 2Tp/4: Se forma un vórtice cerca de la cresta
   • t = 3Tp/4: El vórtice se expande hacia el flujo principal, con flujo completamente invertido
2. Efectos de Anisotropía:
   • λ = 1 (isotrópico): No se forma vórtice en t = 0 y Tp/4
   • λ = 2 (fuertemente anisotrópico): Se observan circulaciones evidentes en todos los tiempos

Resultados del Efecto Combinado

1. Dependencia Temporal: El tiempo de inicio y ubicación de la inversión de flujo dependen de la interacción entre λ y el tiempo t
2. Cambios de Gradiente de Presión: Aparecen gradientes adversos en t = 3Tp/4 y 4Tp/4
3. Espesor de Capa Límite Pulsátil: $δ_{puls} = 1/S$, donde $S = \frac{1}{\sqrt{M}}(α^4 + Wo^2)^{1/4}\cos\left(\frac{1}{2}\tan^{-1}\left(\frac{Wo}{α^2}\right)\right)$

Análisis de Caída de Presión

La caída de presión estacionaria aumenta rápidamente con α² (inverso del número de Darcy), indicando que la baja permeabilidad aumenta significativamente la resistencia al flujo. El aumento de amplitud a también conduce a un aumento en la caída de presión, reflejando los efectos de la rugosidad superficial.

Trabajos Relacionados

Investigación de Flujo en Canales Ondulados

• Tsangaris & Leiter (1984): Método de perturbación para flujo viscoso
• Pozrikidis (1987): Estudio numérico de flujo de Stokes
• Wang & Chen (2002): Transferencia de calor por convección forzada

Flujo en Medios Porosos

• Ng & Wang (2010): Flujo de Darcy-Brinkman
• Wei et al. (2003): Canal ondulado con capa porosa elástica
• Karmakar & Raja Sekhar (2017): Inversión de flujo en medios porosos anisotrópicos

Aplicaciones de Flujo Pulsátil

• Womersley (1955): Base teórica para flujo pulsátil arterial
• Sobey (1980-1983): Investigación de flujo en oxigenadores de membrana
• Griffith et al. (2009): Flujo pulsátil en geometrías estrechadas

Conclusiones y Discusión

Conclusiones Principales

1. Efectos de Anisotropía: El aumento de λ promueve significativamente la separación del flujo, formando regiones de reflujo cerca de crestas y valles
2. Evolución Temporal: La formación y desaparición de vórtices en flujo pulsátil exhibe periodicidad temporal evidente
3. Mecanismo Combinado: El efecto combinado de componentes estacionaria y pulsátil altera el tiempo y ubicación de la separación del flujo
4. Estructura de Capa Límite: Estructura de doble capa con capa de Brinkman cercana a la pared y región de Darcy en campo lejano

Aplicaciones de Ingeniería

1. Oxigenadores de Membrana: La anisotropía y pulsación pueden optimizar la eficiencia de transferencia de oxígeno
2. Diagnóstico Cardiovascular: Proporciona base teórica para análisis de flujo sanguíneo en estenosis arterial y aneurismas
3. Intensificación de Transferencia de Calor: La separación de flujo en superficies onduladas intensifica la mezcla convectiva

Limitaciones

1. Suposición de Parámetro Pequeño: Solo aplicable cuando δ² ≪ 1 y λ²δ² ≪ 1
2. Simplificación Bidimensional: Los efectos tridimensionales en aplicaciones reales no se consideran
3. Linealización: El método de perturbación limita la aplicabilidad a casos de amplitud grande

Direcciones Futuras

1. Extensión Tridimensional: Considerar canales ondulados tridimensionales con variación en dirección transversal
2. Análisis No Lineal: Efectos no lineales en casos de amplitud grande
3. Validación Experimental: Comparación con datos experimentales
4. Diseño Optimizado: Optimización de parámetros geométricos basada en control de flujo

Evaluación Profunda

Fortalezas

1. Rigor Teórico: Método analítico basado en teoría de perturbación con derivación matemática completa
2. Perspectiva Física: Revela el mecanismo profundo de cómo la permeabilidad anisotrópica afecta la separación del flujo
3. Valor de Aplicación: Proporciona herramientas teóricas prácticas para oxigenadores de membrana y flujo cardiovascular
4. Sistematicidad: Análisis completo que abarca flujo estacionario, pulsátil y efectos combinados
5. Teoría de Capa Límite: Establecimiento de la estructura de doble capa Brinkman-Darcy

Deficiencias

1. Restricciones de Parámetros: La suposición de parámetro pequeño limita la universalidad del método
2. Ausencia de Experimentos: Falta de datos experimentales para validar predicciones teóricas
3. Comparación Numérica: Sin comparación con resultados de simulación CFD
4. Características de Materiales: Falta de discusión profunda sobre características anisotrópicas de materiales porosos reales

Impacto

1. Contribución Académica: Proporciona nuevo marco analítico para teoría de flujo en medios porosos
2. Valor de Ingeniería: Tiene perspectivas de aplicación directa en ingeniería biomédica e intensificación de transferencia de calor
3. Innovación Metodológica: El método combinado de perturbación-lubricante puede generalizarse a otros flujos en geometrías complejas

Escenarios de Aplicabilidad

1. Dispositivos Microfluídicos: Control de flujo en canales ondulados a pequeña escala
2. Bioingeniería: Flujo vascular y diseño de órganos artificiales
3. Aplicaciones Industriales: Diseño optimizado de intercambiadores de calor y dispositivos de separación
4. Ingeniería Geológica: Análisis de transporte de fluidos en capas porosas

Referencias Bibliográficas

El artículo cita 62 referencias relacionadas, abarcando logros de investigación importantes en múltiples campos incluyendo mecánica de fluidos, teoría de medios porosos, e ingeniería biomédica, proporcionando una base teórica sólida para esta investigación.