Ascenso transitorio de varilla en un fluido Oldroyd-B

Información Básica

• ID del Artículo: 2501.01016
• Título: Ascenso transitorio de varilla en un fluido Oldroyd-B
• Autores: Tachin Ruangkriengsin, Rodolfo Brandão, Katie Wu, Jonghyun Hwang, Evgeniy Boyko, Howard A. Stone
• Clasificación: physics.flu-dyn (Mecánica de Fluidos)
• Fecha de Publicación: 2 de enero de 2025 (Enviado a arXiv)
• Revista: Journal of Fluid Mechanics (En revisión)
• Enlace del Artículo: https://arxiv.org/abs/2501.01016

Resumen

Este artículo investiga el efecto Weissenberg (fenómeno de ascenso de varilla) en fluidos no newtonianos, es decir, el ascenso de la interfaz del fluido a lo largo de una varilla rotativa. Los autores emplean el modelo Oldroyd-B, utilizando análisis asintótico y variables adimensionales para estudiar la altura de la interfaz dependiente del tiempo. El estudio revela que bajo pequeñas deformaciones, las ecuaciones de gobierno indican la existencia de una capa límite temporal, donde la interfaz asciende rápidamente en escalas de tiempo cortas y luego se aproxima gradualmente al estado estacionario. Mediante la introducción de variables de tiempo estirado, los autores derivan el campo de velocidad transitorio y el perfil de interfaz correspondiente, reintroduciendo posteriormente efectos de inercia pequeños pero finitos para investigar su interacción con la viscoelasticidad.

Antecedentes y Motivación de la Investigación

Definición del Problema

El fenómeno de ascenso de varilla (efecto Weissenberg) es uno de los comportamientos más característicos de los fluidos no newtonianos. Cuando una varilla rota en un fluido viscoelástico, la superficie del fluido asciende a lo largo de la varilla. Este fenómeno aún carece de una comprensión teórica profunda.

Importancia

1. Significado Teórico: Este fenómeno es un concepto fundamental para comprender el comportamiento de fluidos complejos, con importancia significativa en la enseñanza e investigación
2. Valor Práctico: Tiene aplicaciones importantes en mediciones reológicas, permitiendo determinar parámetros de materiales como el coeficiente de esfuerzo normal y el tiempo de relajación
3. Aplicaciones de Ingeniería: Proporciona orientación para aplicaciones industriales de soluciones poliméricas y otros fluidos complejos

Limitaciones de Métodos Existentes

1. Limitaciones Teóricas: Investigaciones previas se han concentrado principalmente en comportamiento en estado estacionario, utilizando modelos de fluidos de segundo orden
2. Problemas del Modelo: Los modelos de fluidos de segundo orden son válidos solo bajo números de Deborah o Weissenberg pequeños, y presentan alta inestabilidad
3. Dependencia Temporal: Falta análisis teórico de la dinámica transitoria

Motivación de la Investigación

Los autores seleccionan el modelo Oldroyd-B para evitar los problemas del modelo de fluido de segundo orden en flujos no estacionarios, e investigan sistemáticamente por primera vez la dependencia temporal del fenómeno de ascenso de varilla.

Contribuciones Principales

1. Primer Análisis Teórico: Proporciona el primer análisis teórico transitorio completo del fenómeno de ascenso de varilla en fluidos Oldroyd-B
2. Teoría de Capa Límite Temporal: Descubre y analiza la estructura de la capa límite temporal, revelando el proceso de dos etapas de ascenso rápido y aproximación lenta al estado estacionario
3. Marco de Análisis Adimensional: Establece un método sistemático de análisis adimensional basado en el número de Weissenberg
4. Investigación de Efectos de Inercia: Investiga por primera vez la interacción entre efectos de inercia pequeños pero finitos y la viscoelasticidad
5. Criterio de Ascenso: Propone un criterio teórico para determinar las condiciones de ocurrencia del ascenso de varilla: $Re < 4\beta_p Wi$

Explicación Detallada de Métodos

Definición de la Tarea

Investigar la evolución temporal de la altura de la superficie libre del fluido $h(r,t)$ cuando una varilla cilíndrica infinitamente larga (radio $a$) rota con velocidad angular $\Omega$ en una solución polimérica viscoelástica.

Modelo Matemático

Ecuaciones de Gobierno

Se adopta el modelo constitutivo Oldroyd-B:

• Ecuación de Evolución del Tensor de Configuración: $\frac{\partial\mathbf{A}}{\partial t} + \mathbf{u} \cdot \nabla\mathbf{A} - \mathbf{A} \cdot \nabla\mathbf{u} - (\nabla\mathbf{u})^T \cdot \mathbf{A} = -\frac{1}{\lambda}(\mathbf{A} - \mathbf{I})$
• Tensor de Esfuerzo: $\boldsymbol{\sigma} = -p\mathbf{I} + 2\mu_s\mathbf{E} + G(\mathbf{A} - \mathbf{I})$
• Ecuación de Cantidad de Movimiento: $\nabla \cdot \boldsymbol{\sigma} = \rho g\mathbf{e}_z$

Adimensionalización

Se introducen parámetros adimensionales clave:

• Número de Weissenberg: $Wi = \lambda\Omega$
• Parámetro de Gravedad: $G = \frac{\rho ga}{(\mu_s + \mu_p)\Omega}$
• Número de Reynolds: $Re = \frac{\rho a^2\Omega}{\mu_s + \mu_p}$

Escalas Características

Los autores identifican la escala característica clave: $G = Wi \cdot \mathcal{G}$, donde $\mathcal{G} = O(1)$, lo que simplifica la expansión asintótica de dos parámetros a una expansión de un parámetro.

Puntos de Innovación Técnica

1. Separación de Escalas de Tiempo

Se descubren dos escalas de tiempo:

• Tiempo Externo: $T = O(1)$, correspondiente al comportamiento en estado estacionario
• Tiempo Interno: $\tau = T/Wi$, correspondiente a la capa límite transitoria

2. Método de Perturbación de Dominio

Para pequeñas deformaciones de interfaz, se expande en serie de Taylor las condiciones de frontera desde $Z = H(R,T)$ hasta $Z = 0$, simplificando el tratamiento de las condiciones de frontera de superficie libre.

3. Transformación de Weber

Al tratar efectos de inercia, se aplica la transformación de Weber para procesar el operador diferencial de Bessel en dominios semiinfinitos, transformando ecuaciones diferenciales parciales en ecuaciones diferenciales ordinarias.

Configuración Experimental

Marco de Análisis Teórico

Este artículo es una investigación puramente teórica que utiliza métodos de análisis asintótico:

1. Parámetro Pequeño: $Wi \ll 1$ (viscoelasticidad débil)
2. Rango de Gravedad: $G \gg Wi^{-1}$ (efectos gravitacionales fuertes)
3. Suposición de Pequeña Deformación: Deformación de interfaz relativamente pequeña

Estrategia de Análisis

1. Sección 4: Análisis Externo - Perfil de interfaz en estado estacionario
2. Sección 5: Análisis Interno - Perfil de interfaz transitorio
3. Sección 6: Análisis que incluye efectos de inercia pequeños

Resultados Experimentales

Resultados Principales

1. Perfil de Interfaz en Estado Estacionario

$H_{steady}(R) = \frac{2\beta_p Wi^2}{\mathcal{G}R^4}$

Esto es consistente con el resultado obtenido por Joseph & Fosdick (1973) utilizando el modelo de fluido de segundo orden.

2. Perfil de Interfaz Transitorio

$H_{transient}(R,\tau) = \frac{2\beta_p(1-(\tau+1)e^{-\tau})Wi^2}{\mathcal{G}R^4}$

donde $\tau = T/Wi$ es la variable de tiempo estirado.

3. Criterio de Ascenso

Cuando se incluyen efectos de inercia, la condición para que ocurra el ascenso de varilla es: $4El\beta_p > 1$ donde $El = Wi/Re$ es el número de elasticidad.

Características de Evolución Temporal

1. Etapa de Ascenso Rápido: La interfaz asciende rápidamente dentro de $\tau \approx 6$ unidades de tiempo
2. Aproximación al Estado Estacionario: La solución transitoria converge a la solución en estado estacionario en el límite de tiempo largo
3. Comportamiento No Monótono: Cuando se incluye inercia, la altura de la interfaz puede primero ascender y luego descender

Hallazgos Clave

1. Esfuerzo Circunferencial: El esfuerzo circunferencial generado por la extensión del polímero es el mecanismo principal que impulsa el ascenso
2. Competencia de Inercia: La fuerza centrífuga producida por la inercia compite con el esfuerzo circunferencial de la viscoelasticidad
3. Estructura de Capa Límite: La existencia de la capa límite temporal explica el comportamiento de dos etapas de transitorio rápido y estado estacionario lento

Trabajo Relacionado

Desarrollo Histórico

1. Serrin (1959): Primer intento teórico, utilizando modelo de fluido de segundo orden
2. Giesekus (1961): Ecuaciones constitutivas más generales ignorando efectos de inercia
3. Joseph & Fosdick (1973): Marco sistemático del método de perturbación de dominio
4. Oldroyd (1950, 1951): Aplicación temprana del modelo Oldroyd-B

Ventajas de Este Artículo

1. Dependencia Temporal: Primer análisis transitorio completo
2. Estabilidad del Modelo: Evita los problemas de inestabilidad del modelo de fluido de segundo orden
3. Marco Adimensional: Análisis sistemático de parámetros adimensionales
4. Efectos de Inercia: Primer estudio de la interacción entre inercia y viscoelasticidad

Conclusiones y Discusión

Conclusiones Principales

1. Estructura de Capa Límite: El fenómeno de ascenso de varilla presenta una capa límite temporal clara, donde la interfaz asciende rápidamente y luego se aproxima lentamente al estado estacionario
2. Separación de Escalas: La escala de tiempo transitoria es $O(Wi)$, mientras que la escala de tiempo en estado estacionario es $O(1)$
3. Mecanismo de Ascenso: El esfuerzo circunferencial generado por la extensión del polímero es el mecanismo central que impulsa el ascenso
4. Influencia de la Inercia: Los efectos de inercia pequeños pueden influir significativamente en el comportamiento de ascenso, existiendo un número de elasticidad crítico

Limitaciones

1. Suposición de Pequeña Deformación: Solo aplicable a casos de pequeña deformación de interfaz
2. Viscoelasticidad Débil: Limitado al rango $Wi \ll 1$
3. Gravedad Fuerte: Requiere la condición $G \gg Wi^{-1}$
4. Modelo Simplificado: Ignora efectos de tensión superficial y ángulo de contacto

Direcciones Futuras

1. Número de Weissenberg Alto: Investigar el comportamiento cuando $Wi \gg 1$
2. Modelos Constitutivos Complejos: Extender a modelos FENE-P, Giesekus, etc.
3. Deformación Finita: Tratar casos de deformación de interfaz grande
4. Validación Experimental: Comparación y verificación con datos experimentales

Evaluación Profunda

Fortalezas

1. Rigor Teórico: Derivaciones matemáticas rigurosas, aplicación apropiada de métodos de análisis asintótico
2. Perspectiva Física: Revela la estructura de capa límite temporal y el mecanismo de ascenso
3. Innovación Metodológica: Combinación ingeniosa de transformación de Weber y método de perturbación de dominio
4. Valor Práctico: Proporciona base teórica para mediciones reológicas

Deficiencias

1. Rango de Parámetros: Rango de parámetros aplicables relativamente limitado
2. Comparación Experimental: Falta comparación directa con datos experimentales
3. Verificación Numérica: Algunos resultados solo presentan cálculos numéricos, careciendo de expresiones analíticas

Impacto

1. Contribución Académica: Llena el vacío en la teoría transitoria del ascenso de varilla
2. Valor Metodológico: Los métodos de análisis proporcionados pueden generalizarse a otros problemas de fluidos complejos
3. Apoyo de Ingeniería: Proporciona apoyo teórico para el diseño de reómetros y caracterización de materiales

Escenarios Aplicables

1. Soluciones Poliméricas Diluidas: Fluidos viscoelásticos como fluidos de Boger con viscosidad de corte constante
2. Velocidades de Rotación Bajas: Condiciones experimentales con número de Weissenberg pequeño
3. Mediciones Reológicas: Determinación del coeficiente de esfuerzo normal y tiempo de relajación

Referencias

1. Joseph, D.D. & Fosdick, R.L. (1973). The free surface on a liquid between cylinders rotating at different speeds. Arch. Rat. Mech. Anal., 49, 321-380.
2. Oldroyd, J.G. (1950). On the formulation of rheological equations of state. Proc. R. Soc. A, 200, 523-541.
3. Bird, R.B., Armstrong, R.C. & Hassager, O. (1987). Dynamics of Polymeric Liquids, Vol. 1, 2nd edn. John Wiley and Sons.

Evaluación General: Este es un artículo de alta calidad de mecánica de fluidos teórica que analiza sistemáticamente por primera vez el comportamiento transitorio del fenómeno de ascenso de varilla en fluidos Oldroyd-B. El artículo presenta derivaciones matemáticas rigurosas, perspectivas físicas profundas y métodos innovadores, realizando contribuciones importantes a la teoría de fluidos complejos y aplicaciones reológicas. Aunque existen limitaciones en el rango de parámetros, su marco teórico y métodos de análisis poseen excelente valor de generalización.