Inestabilidad espaciotemporal, advección y panorama de turbulencia

Año 2 · Módulo 15: Dinámica espacial y sistemas continuos

148.1. Cuando la perturbación cambia de lugar mientras cambia de intensidad

Una perturbación puede crecer sin permanecer donde nació. Pensemos en una pequeña mancha de colorante introducida en un canal de agua. La difusión ensancha la mancha; la corriente transporta su centro. Si además ocurre una reacción que aumenta localmente la concentración de una sustancia, el máximo puede crecer mientras viaja corriente abajo. Un sensor instalado en un punto fijo quizá registre una perturbación breve que después desaparece, aunque la nube siga aumentando su amplitud a medida que se aleja. Preguntar solo si la perturbación crece o decrece deja fuera una parte esencial de la dinámica.

Hasta ahora hemos estudiado campos, difusión, modos espaciales y sistemas capaces de formar patrones o propagar frentes. En los casos más simples, la posición de una estructura es una referencia relativamente estable. Muchos sistemas físicos, sin embargo, presentan un movimiento de transporte superpuesto a la producción local de estructuras. Otros desarrollan irregularidades que cambian continuamente tanto en el espacio como en el tiempo. Estas situaciones plantean preguntas nuevas. ¿Crece la perturbación para un observador que se mueve con ella? ¿Crece para alguien que espera en un lugar? ¿Cuánto depende la respuesta del tamaño del sistema y de las fronteras por donde entra o sale material?

La palabra inestabilidad precisa siempre de un objeto y un criterio. Una solución uniforme puede ser linealmente inestable si pequeñas perturbaciones aumentan según una norma especificada. Un perfil estacionario puede sufrir una perturbación que se amplifique transitoriamente y luego desaparezca. Un régimen de flujo puede mantener fluctuaciones pese a que el sistema global conserva ciertas cantidades. Estas alternativas no se ordenan necesariamente en una escala simple que conduzca de estabilidad a turbulencia mediante un único parámetro. La geometría y los mecanismos físicos deciden el camino concreto.

El término turbulencia también requiere disciplina conceptual. En mecánica de fluidos designa un conjunto de comportamientos irregulares de los campos de velocidad, caracterizados por fluctuaciones espaciales y temporales, interacciones no lineales y transporte entre escalas, entre otros rasgos. En algunos modelos de reacción–difusión y ecuaciones de amplitud se habla de turbulencia de fase o de defectos en regímenes complejos que no implican necesariamente movimiento turbulento de un líquido. El caos espaciotemporal, por su parte, es una dinámica determinista extendida con irregularidad y sensibilidad a perturbaciones distribuidas. Hay conexiones entre estos conceptos, pero sus criterios empíricos y matemáticos no son intercambiables.

Examinaremos primero el transporte dirigido y la diferencia entre marcos de observación. Después estudiaremos un modelo que permite calcular exactamente cuándo una amplificación viaja y cuándo permanece activa en un punto fijo. Finalmente compararemos la irregularidad de patrones extendidos con la turbulencia hidrodinámica y veremos qué podemos inferir realmente de una serie temporal, un campo registrado o una simulación numérica.

148.2. La advección según se siga un campo o una partícula

Sea $u\left( \mathbf{x},t \right)$ la concentración de una sustancia y $\mathbf{v}\left( \mathbf{x},t \right)$ el campo de velocidades del medio que la transporta. El cambio temporal de $u$ medido por un observador fijo se representa con $\partial_{t}u$. Una pequeña porción material que se mueve según $d\mathbf{x}/dt = \mathbf{v}$ experimenta, en cambio, la derivada material

El segundo término es la advección de la propiedad medida. La derivación aplica la regla de la cadena a lo largo de la trayectoria $\mathbf{x}(t)$, de modo que la concentración varía porque el campo cambia con el tiempo y porque la partícula visita posiciones diferentes. La expresión no presupone turbulencia; una corriente constante y perfectamente uniforme también advecta una mancha.

Para una concentración conservada, una ecuación general de balance que incorpora el flujo advectivo $u\mathbf{v}$, el flujo difusivo $- D\nabla u$ y, si corresponde, una fuente $R$ es

Cuando $\nabla \cdot \mathbf{v} = 0$, el transporte advectivo del lado izquierdo se reduce a $\mathbf{v} \cdot \nabla u$. Si el fluido es compresible, esa sustitución omite el término $u\nabla \cdot \mathbf{v}$, relevante para la concentración por unidad de volumen. La forma conservativa permite expresar correctamente lo que entra y sale de un volumen de control y suele ser la apropiada para una discretización por volúmenes finitos.

Imaginemos una corriente unidimensional constante $v > 0$, sin difusión ni reacciones. La ecuación $\partial_{t}u + v\partial_{x}u = 0$ tiene solución $u(x,t) = u_{0}(x - vt)$, de modo que el perfil se desplaza una distancia $vt$ sin cambiar su forma. Un sensor fijo observa que la señal llega y se va, mientras que un observador que se mueve con velocidad $v$ ve un campo estacionario; ambos describen la misma evolución física desde marcos diferentes. La cuestión se complica cuando la velocidad depende de la posición, porque la deformación y el estiramiento de una mancha pueden intensificar gradientes, que a su vez facilitan la difusión.

Advección y difusión son mecanismos distintos. La primera puede redistribuir una sustancia mediante el movimiento del medio; la segunda es un flujo respecto del medio debido a gradientes de concentración. En un fluido incompresible ideal y sin difusión, la concentración de cada partícula material permanece constante. Aun así, el patrón visto desde un laboratorio puede estirarse y plegarse de manera muy compleja. Cuando existe difusión, los gradientes producidos por el estiramiento determinan a qué velocidad se borran las diferencias más finas.

Esta distinción también evita confundir una onda con el simple desplazamiento de un marcador. Una onda química puede propagarse respecto del propio fluido gracias a la reacción y a la difusión, mientras el fluido transporta todo el patrón. Medir la velocidad en el laboratorio sin conocer la velocidad local del medio puede mezclar ambas contribuciones. Si el flujo es homogéneo y el modelo permite superponerlas, la velocidad observada del frente cambia aproximadamente en la dirección de la advección; con flujo espacialmente variable o reacción compleja hace falta estudiar la ecuación completa.

148.3. Escalas de transporte y números de Péclet y de Reynolds

La advección y la difusión no compiten siempre con la misma intensidad. Elegida una longitud característica $L$, una velocidad $U$ y un coeficiente de difusión molecular $D$, podemos comparar un tiempo advectivo $t_{a} = L/U$ con un tiempo difusivo $t_{d} = L^{2}/D$. Su cociente es el número de Péclet,

Si $Pe \ll 1$, la difusión puede homogeneizar variaciones de escala $L$ más rápidamente que el transporte advectivo sobre esa misma distancia. Si $Pe \gg 1$, el campo suele recorrer esa escala antes de difundirse apreciablemente. El juicio depende de la longitud escogida; una nube de gran tamaño puede estar dominada por transporte dirigido mientras que sus filamentos microscópicos se alisan por difusión.

Por ejemplo, con $U = 2\ cm/s$, $L = 5\ cm$ y $D = 0,1\ cm^{2}/s$, obtenemos $Pe = 100$. El traslado de cinco centímetros ocurre en unos $2,5$ segundos, mientras que la escala difusiva correspondiente es $250$ segundos. Es de esperar que la mancha recorra esa longitud antes de dispersarse uniformemente en ella. La inferencia sería distinta para una microescala o si el valor de $D$ correspondiera a una dispersión efectiva que ya incluye mezcla inducida por el flujo. Confundir difusión molecular con un coeficiente empírico ajustado puede hacer que dos modelos aparentemente idénticos representen procesos diferentes.

En hidrodinámica aparece además el número de Reynolds, $Re = UL/\nu$, donde $\nu$ es la viscosidad cinemática. Este número compara los términos inerciales y viscosos en la ecuación de movimiento del fluido. Aunque ambos números tienen forma matemática parecida, miden relaciones distintas. Un fluido puede transportar muy eficazmente un soluto con $Pe$ elevado sin que su velocidad sea turbulenta. Recíprocamente, un flujo con fluctuaciones intensas puede transportar sustancias con difusividades muy distintas.

La relación entre ambos números se expresa mediante el número de Schmidt, $Sc = \nu/D$, con $Pe = Re\, Sc$ cuando usamos las mismas escalas. Esta igualdad ayuda a identificar cuándo el transporte de una propiedad y el transporte de cantidad de movimiento presentan balances diferentes. De ella no se obtiene un umbral de turbulencia, porque los valores críticos de Reynolds, cuando existen, dependen de geometría, perturbaciones, perfil de velocidad y condiciones de frontera.

La reacción agrega un tercer reloj. Para una tasa característica $r$, podemos formar un tiempo de reacción $t_{r} = 1/|r|$ si $r \neq 0$. Comparar $t_{r}$ con $t_{a}$ y $t_{d}$ permite anticipar si una especie se transforma antes de ser transportada o difundida. En sistemas con varias reacciones y especies, una tasa única puede ser una simplificación insuficiente. El análisis dimensional reduce la variedad de parámetros y ayuda a seleccionar experimentos pertinentes; las transiciones efectivas deben determinarse mediante la dinámica.

148.4. Inestabilidad convectiva y absoluta en un modelo resoluble

Consideremos una perturbación escalar $a(x,t)$ que crece localmente a tasa $r$, se difunde con coeficiente $D > 0$ y viaja en una corriente uniforme de velocidad $c$,

Aunque este modelo lineal para una amplitud pequeña no describe por completo un fluido turbulento, permite separar exactamente crecimiento y transporte. Si buscamos un modo $a \propto \exp(\sigma t + ikx)$, obtenemos la relación de dispersión

Su parte real, $r - Dk^{2}$, describe amplificación o amortiguación de modos periódicos; la parte imaginaria codifica su transporte. Para $r > 0$, los modos de longitud de onda suficientemente grande crecen. Aun así, esta afirmación modal no responde por completo a lo que ocurre con una perturbación localizada ante un sensor inmóvil.

Si el estado inicial es una perturbación idealmente puntual de masa unidad en la recta infinita, la solución fundamental es

El centro de la perturbación viaja con velocidad $c$, su ancho crece como $\sqrt{Dt}$ y su amplitud sobre la trayectoria del centro se multiplica por $e^{rt}/\sqrt{t}$. Si fijamos ahora una posición $x$ y dejamos que $t$ aumente, al desarrollar el exponente se obtiene como tasa temporal dominante $r - c^{2}/(4D)$. Si $0 < r < c^{2}/(4D)$, el paquete se amplifica mientras viaja corriente abajo, pero la respuesta a un impulso localizado decae en cualquier lugar fijo. Se habla entonces de inestabilidad convectiva para este problema. Si $r > c^{2}/(4D)$, la señal producida por el impulso aumenta también en un punto fijo y la inestabilidad es absoluta. En la igualdad hay un caso umbral en el que debe considerarse el factor algebraico.

Podemos calcularlo sin simulación. Si $r = 0,5\ s^{- 1}$, $c = 1\ cm/s$ y $D = 0,25\ cm^{2}/s$, el valor crítico es $c^{2}/(4D) = 1\ s^{- 1}$. En un sensor fijo, la envolvente asintótica decae con tasa exponencial de $- 0,5\ s^{- 1}$, aunque el centro móvil crezca con tasa $0,5\ s^{- 1}$. Si mantenemos los demás parámetros y cambiamos $r$ a $1,5\ s^{- 1}$, la tasa en el lugar fijo pasa a $+ 0,5\ s^{- 1}$.

Estas conclusiones valen para un impulso localizado en un dominio infinito con coeficientes constantes. Una fuente que perturba continuamente la entrada de un canal puede mantener fluctuaciones visibles incluso en régimen convectivo. Las reflexiones, las fronteras, la saturación no lineal y los perfiles espaciales de velocidad pueden producir un comportamiento global distinto. Para sistemas generales, la distinción exige estudiar la evolución de paquetes de ondas y, a menudo, técnicas de análisis espectral más elaboradas.

148.5. Crecimiento transitorio, selección de ondas y papel de las fronteras

Un modo que crece exponencialmente no es la única explicación posible de un aumento temporal de amplitud. En sistemas lineales con variables acopladas, los modos propios pueden ser no ortogonales. Dos contribuciones que inicialmente casi se cancelan pueden separarse y producir un incremento importante antes de que ambas terminen decayendo. Este fenómeno recibe el nombre de crecimiento transitorio y puede existir aunque todos los autovalores relevantes tengan parte real negativa. El estudio de Trefethen y colaboradores sobre estabilidad hidrodinámica mostró por qué confiar únicamente en autovalores puede ser insuficiente para entender la sensibilidad de ciertos flujos.

Un ejemplo algebraico sencillo es el sistema $\dot{x} = - x + Ky$, $\dot{y} = - 2y$, con $K > 0$. Sus autovalores son $- 1$ y $- 2$; todas las soluciones tienden finalmente a cero. Si comenzamos, sin embargo, con $x(0) = 0$, $y(0) = 1$, obtenemos $y(t) = e^{- 2t}$ y

La variable $x$ aumenta desde cero, alcanza en $t = \ln 2$ un máximo de valor $K/4$ y luego decae. Con $K$ grande puede haber una amplificación importante de un componente sin inestabilidad exponencial asintótica. El ejemplo no es espacial, pero exhibe el mecanismo de no normalidad que también puede intervenir al linearizar ecuaciones de flujo. En sistemas no lineales, una excursión transitoria suficientemente grande puede incluso desencadenar una transición hacia otro régimen, aunque esa posibilidad queda fuera del alcance de las ecuaciones lineales.

Las perturbaciones propagantes requieren además distinguir velocidades. La velocidad de fase describe el desplazamiento de una cresta de fase constante; la velocidad de grupo describe la propagación de la envolvente de un paquete bajo aproximaciones de dispersión débil y espectro concentrado. En un medio dispersivo pueden ser distintas. La velocidad con que un pulso invade una región inicialmente estable o inestable es otra magnitud, que puede depender de la selección no lineal del frente. Usar la velocidad de una cresta para inferir cómo se desplaza la energía o la información es incorrecto sin hipótesis adicionales.

En un canal abierto, la entrada puede suministrar perturbaciones y la salida evacuarlas. Una inestabilidad convectiva puede comportarse como un amplificador de ruido; el campo aumenta corriente abajo mientras la región inicial recupera su estado. Una inestabilidad absoluta favorece, bajo las condiciones apropiadas, actividad que crece en un lugar y puede participar en la formación de un modo global. En dominios finitos importan también las reflexiones y la variación espacial del estado de base. Las fronteras pueden convertir una respuesta aparentemente local en una oscilación autoalimentada o suprimir una onda que en un dominio infinito sería posible.

Por ello, además de registrar el crecimiento de un máximo, hay que medir su posición y evaluar señales en lugares fijos. Dos perturbaciones con igual máxima amplitud pueden tener consecuencias distintas si una cruza la región de interés en segundos y la otra permanece amplificándose junto a una estructura sensible. La elección del marco de observación forma parte de la pregunta física.

148.6. De una onda inestable al caos espaciotemporal

Una inestabilidad lineal informa qué perturbaciones empiezan a crecer, pero no resuelve su evolución a amplitud finita. Cuando la perturbación altera los propios términos que la hacen crecer, la dinámica no lineal puede limitar su amplitud, desplazar energía entre modos o crear nuevas estructuras. En un sistema espacialmente extendido, la interacción alcanza muchos grados de libertad (una región puede oscilar mientras otra desarrolla un defecto, varias ondas pueden competir y la geometría local puede cambiar después de cada colisión). El resultado no tiene por qué ser una textura fija como la de ciertos patrones de Turing.

La ecuación de Kuramoto–Sivashinsky es un ejemplo estilizado de esta competencia. En una convención adimensional para un campo real $h(x,t)$, puede escribirse

Al linearizar alrededor del estado uniforme $h = 0$, una perturbación proporcional a $e^{\sigma t + ikx}$ tiene tasa $\sigma(k) = k^{2} - k^{4}$. Las longitudes de onda grandes, correspondientes a $0 < |k| < 1$, crecen; las de número de onda elevado se amortiguan. El término no lineal $h\partial_{x}h$ redistribuye la estructura espacial y puede impedir que el crecimiento de los modos inestables continúe indefinidamente. En dominios suficientemente extensos y para ciertos regímenes se observan fluctuaciones espaciotemporales complejas.

Incluso en este ejemplo interviene el tamaño del dominio. Si imponemos periodicidad sobre un intervalo de longitud $L$, los números de onda permitidos son $k_{n} = 2\pi n/L$, con $n$ entero. Si $L \leq 2\pi$, no hay modos espaciales no constantes dentro de la banda linealmente inestable. En un dominio más largo pueden caber modos que sí crezcan. De esta observación no se sigue que cualquier dominio apenas mayor que $2\pi$ sea necesariamente caótico, porque entre la inestabilidad inicial y el caos pueden existir soluciones estacionarias, periódicas u otros regímenes seleccionados por las condiciones iniciales y los parámetros.

Otra familia central es la ecuación compleja de Ginzburg–Landau, usada como modelo de amplitud cerca de ciertas inestabilidades oscilatorias. Una de sus formas normalizadas es $\partial_{t}A = A + (1 + ib)\nabla^{2}A - (1 + ic)|A|^{2}A$, donde $A$ es complejo y $b,c$ son parámetros reales. El módulo y la fase de $A$ permiten diferenciar regímenes en los que la fase cambia irregularmente sin que la amplitud se anule de otros donde surgen ceros de amplitud, que son defectos. Esa variedad de comportamientos explica su uso en el estudio de ondas y patrones, aunque cada aplicación requiere comprobar que la reducción a una amplitud compleja está justificada.

Estas ecuaciones permiten investigar caos espaciotemporal sin simular todas las moléculas ni resolver un flujo hidrodinámico tridimensional. Su utilidad está en aislar mecanismos dinámicos. Sería un abuso convertir cualquier fluctuación irregular en solución de un modelo universal. El adjetivo caótico debe sostenerse con análisis temporal y espacial; una fotografía de textura desordenada también puede proceder de ruido, heterogeneidad fija o errores numéricos.

Este caos debe separarse de la turbulencia hidrodinámica. Ottino (1990) mostró que un campo de velocidad puede mezclar partículas mediante advección caótica incluso cuando no existe turbulencia plenamente desarrollada. De la sola irregularidad espacial tampoco pueden deducirse la cascada de energía ni las propiedades estadísticas de un flujo turbulento.

148.7. Turbulencia hidrodinámica, interacción de escalas y disipación

En un fluido newtoniano incompresible de densidad constante, el campo de velocidad $\mathbf{v}$ se describe, bajo las hipótesis apropiadas, mediante las ecuaciones de Navier–Stokes

La aceleración convectiva $\left( \mathbf{v} \cdot \nabla \right)\mathbf{v}$ es no lineal, porque el campo transporta su propia cantidad de movimiento. La presión impone la restricción de incompresibilidad y redistribuye fuerzas; la viscosidad amortigua gradientes de velocidad. En un flujo turbulento tridimensional, estructuras de tamaños diferentes interactúan y pueden transferir energía entre escalas. El mecanismo de transferencia y la disipación son propiedades del campo de velocidad, por lo que no deben confundirse con la difusión pasiva de un trazador, aunque ambos procesos puedan darse a la vez.

El balance energético muestra una restricción elemental. En un dominio periódico y sin fuerzas externas, el producto escalar de la ecuación con la velocidad, integrado espacialmente, lleva a

suponiendo suficiente regularidad y condiciones para que desaparezcan los flujos de frontera. La viscosidad disipa energía cinética; el término no lineal la redistribuye espacialmente y entre modos, pero no la crea en el balance total idealizado. Si existe una fuerza externa o trabajo sobre las fronteras, aparece una fuente de energía que puede mantener las fluctuaciones. Una turbulencia sostenida en régimen estacionario necesita, por tanto, un aporte compatible con la disipación promedio.

En ciertos flujos tridimensionales a Reynolds elevado, las ideas de Kolmogórov permiten describir un intervalo de escalas donde las estadísticas dependen principalmente de la tasa media de disipación $\epsilon$ y del número de onda $k$. El espectro de energía puede aproximarse allí por $E(k) \propto \epsilon^{2/3}k^{- 5/3}$. La relación requiere hipótesis sobre el régimen, el rango de escalas y la definición del espectro. No es una firma universal que deba aparecer en cualquier fluido turbulento ni, aislada, una prueba suficiente de turbulencia. Las condiciones de forzamiento, la anisotropía, la estratificación y la dimensionalidad pueden alterar radicalmente la transferencia energética.

En la dinámica bidimensional, la energía y la enstrofía están sometidas a restricciones de conservación diferentes en el límite inviscido, lo que favorece escenarios de transferencia espectral distintos de los habituales en tres dimensiones. Esta diferencia impide extrapolar sin reservas la intuición de una cascada de energía hacia escalas pequeñas a cualquier ecuación de flujo. También ayuda a explicar por qué una simulación bidimensional visualmente atractiva puede representar de manera insuficiente un problema hidrodinámico esencialmente tridimensional.

148.8. Qué medir en un campo irregular

La irregularidad espacial y temporal debe describirse mediante observables definidos. Una primera posibilidad es medir la media y la varianza de un campo después de separar, cuando tenga sentido, componentes estacionarias y fluctuantes. Otra consiste en calcular la correlación entre valores registrados en dos puntos y tiempos diferentes. Para un campo estacionario y aproximadamente homogéneo, una función de correlación puede escribirse como $C\left( \mathbf{r},\tau \right) = \langle\delta u\left( \mathbf{x},t \right)\,\delta u\left( \mathbf{x} + \mathbf{r},t + \tau \right)\rangle$, donde $\delta u = u - \langle u\rangle$. Hay que especificar si los promedios son temporales, espaciales, entre realizaciones o alguna combinación justificada.

Una correlación que disminuye con la separación puede definir una escala característica de organización espacial. Una correlación temporal informa cuánto persiste la memoria estadística de las fluctuaciones. Con solo estas correlaciones no puede garantizarse que la dinámica sea caótica, porque un proceso aleatorio con correlación corta puede imitar estadísticas de un sistema determinista complejo. En el sentido inverso, un sistema caótico puede mantener estructuras coherentes que generan correlaciones a distancias o tiempos relativamente grandes.

La advección modifica la lectura de esas mediciones. Una estructura trasladada sin deformación pierde correlación rápidamente para un sensor fijo, aunque sea casi idéntica cuando la comparamos en un marco que se mueve con ella. En un diagrama espacio-tiempo, una perturbación transportada forma trazos inclinados cuya pendiente permite estimar su velocidad aparente. Por eso puede ser útil examinar $C\left( \mathbf{r},\tau \right)$ cerca de desplazamientos $\mathbf{r} \approx \mathbf{U}\tau$. La conocida hipótesis de turbulencia congelada de Taylor aproxima ciertas fluctuaciones temporales como estructuras transportadas por una velocidad media, pero exige que evolucionen lentamente durante el transporte; falla si se deforman intensamente o si la velocidad advectiva fluctúa de manera comparable a la media.

Otra herramienta es el espectro espacial o espaciotemporal. Un espectro concentrado en pocos números de onda sugiere estructura periódica predominante; uno más distribuido puede reflejar múltiples escalas. Sin embargo, espectros amplios también surgen por ruido de medición, discontinuidades o artefactos de ventanas finitas. Los exponentes de Lyapunov y otros indicadores de sensibilidad a condiciones iniciales pueden ayudar a estudiar caos determinista, pero estimarlos a partir de pocos sensores y registros cortos es especialmente difícil en sistemas extendidos.

Antes de interpretar cualquier indicador conviene preguntar qué escala espacial resuelve el instrumento, cuánta duración efectiva tiene la observación y si el muestreo evita aliasing. También debemos verificar si la señal es suficientemente estacionaria para calcular promedios. Cambios lentos de parámetros pueden generar fluctuaciones aparentes que se atribuirían erróneamente a una dinámica autónoma. Una comparación sólida suele usar más de un observable y comprueba que las conclusiones sobrevivan a decisiones razonables de filtrado y segmentación.

148.9. Simular transporte e irregularidad sin fabricar turbulencia numérica

Las ecuaciones de advección–difusión son particularmente sensibles a cómo representamos el transporte en una malla. Un esquema centrado sencillo puede producir oscilaciones artificiales cerca de gradientes abruptos cuando la advección domina. Un esquema demasiado disipativo puede borrar frentes genuinos y hacer que el sistema parezca más estable de lo que es. En otros casos, la discretización crea patrones que desaparecen al refinar la resolución. Si la conclusión científica depende de la presencia de irregularidades, estas diferencias son un problema de interpretación física.

Hay restricciones temporales conocidas. Para la advección unidimensional con velocidad constante positiva, un esquema explícito de primer orden a contracorriente tiene una condición de estabilidad característica $c\Delta t/\Delta x \leq 1$. Para difusión pura en una dimensión, Euler explícito con diferencias centrales exige $D\Delta t/(\Delta x)^{2} \leq 1/2$. Cuando se combinan procesos, fronteras y varios componentes, deben examinarse las restricciones del esquema completo. Un algoritmo estable en advección pura y otro estable en difusión pura no garantizan automáticamente que su combinación sea estable.

La conservación da pruebas independientes del aspecto visual. En un sistema sin fuentes ni flujo a través de la frontera, la integral de una concentración conservada debe mantenerse. Si el modelo incluye una reacción, la variación de esa integral debe coincidir con la producción o consumo neto correspondiente. La violación persistente de un balance puede señalar un error en el código o en la implementación de una condición de frontera. También es necesario comprobar la convergencia con tamaño de malla y paso temporal. Una aparente ley de potencia observada sobre menos de una década de escalas exige especial cautela.

En turbulencia hidrodinámica se suman otras dificultades, como la resolución de escalas disipativas, el tratamiento de presión e incompresibilidad, las condiciones de pared y la representación del forzamiento. La simulación directa de todas las escalas activas puede ser computacionalmente prohibitiva. En ese caso se usan métodos que modelan parte de las escalas no resueltas, como la simulación de grandes remolinos; el resultado depende de esas decisiones de cierre. Un campo obtenido con un modelo de submalla no debe presentarse como si resolviera exhaustivamente todos los movimientos microscópicos.

La validación requiere problemas de referencia pertinentes. Una gaussiana advectada y difundida, cuyo centro, ancho y amplitud pueden calcularse, prueba aspectos del transporte. Un flujo laminar con solución conocida examina otros elementos del código hidrodinámico. Para una dinámica irregular, las comparaciones con estadísticas experimentales o resultados numéricos de referencia son necesarias, pero deben conservar los mismos observables y condiciones. Reproducir imágenes parecidas de remolinos demuestra mucho menos que recuperar balances, espectros, correlaciones y respuestas a cambios paramétricos dentro de sus incertidumbres.

148.10. Una guía para distinguir transporte, amplificación y turbulencia

Supongamos que un equipo registra una banda de fluctuaciones que avanza por un canal. La primera pregunta es si hay un transporte medio capaz de explicar su movimiento. Se puede medir la velocidad de la corriente o del material portador y compararla con la velocidad observada de la banda. Si ambas coinciden aproximadamente, la advección es una hipótesis plausible, aunque todavía falta explicar su formación. Si la estructura se mueve respecto del medio, quizá intervengan propagación de ondas, reacción o inestabilidades que seleccionen una velocidad propia.

La segunda pregunta es qué ocurre con las perturbaciones pequeñas. Se pueden introducir estímulos localizados de amplitud controlada, seguir su centro y medir la evolución de su envolvente tanto en un marco móvil como en lugares fijos. Un crecimiento sobre una trayectoria desplazada, acompañado por decaimiento en sensores inmóviles, orienta hacia amplificación convectiva. Un aumento sostenido en una ubicación fija puede sugerir inestabilidad absoluta o un mecanismo de realimentación global, pero la distinción requiere excluir forzamiento continuo desde la entrada y efectos de frontera. Cuando las perturbaciones son grandes, el comportamiento puede dejar de seguir la predicción lineal.

La tercera pregunta trata de la organización de las fluctuaciones. Si se mantienen pocos modos oscilatorios, conviene buscar periodicidad y relaciones de fase. Si aparecen múltiples defectos y escalas, es necesario medir su creación, sus interacciones y las correlaciones espaciales. Si se investiga turbulencia de un fluido, también hay que estudiar el campo de velocidad, las escalas dinámicas, el aporte energético y la disipación. Una concentración irregular transportada por un flujo laminar no prueba turbulencia; un patrón químico caótico tampoco demuestra que el fluido que lo contiene tenga un movimiento turbulento.

Las mediciones deben contrastarse con la incertidumbre instrumental y con la sensibilidad al dominio. Cambiar la posición de una entrada, el tamaño de la preparación o la condición de salida puede transformar una amplificación pasajera en una oscilación persistente. Del mismo modo, la reducción de la malla numérica puede hacer desaparecer una irregularidad que parecía físicamente convincente. Atribuir una transición al mecanismo correcto exige estudiar esas dependencias antes de anunciarla como propiedad general del medio.

Las herramientas desarrolladas permiten ordenar los mecanismos en juego. La difusión elimina gradientes finos, la advección desplaza y deforma estructuras, las reacciones pueden amplificar perturbaciones y los acoplamientos no lineales pueden generar dinámicas prolongadas de muchos modos. En fluidos, la transferencia de energía entre escalas agrega mecanismos específicos de turbulencia. Ninguna de estas piezas permite inferir las otras a partir de la apariencia de una imagen. El siguiente problema práctico será cómo discretizar el espacio sin introducir artificialmente inestabilidades o disipación que cambien la respuesta que queríamos estudiar.

Referencias

• Aranson, I. S., y Kramer, L. (2002). The world of the complex Ginzburg–Landau equation. Reviews of Modern Physics, 74, 99–143. https://doi.org/10.1103/RevModPhys.74.99
• Cross, M. C., y Hohenberg, P. C. (1993). Pattern formation outside of equilibrium. Reviews of Modern Physics, 65, 851–1112. https://doi.org/10.1103/RevModPhys.65.851
• Frisch, U. (1995). Turbulence: The Legacy of A. N. Kolmogorov. Cambridge University Press.
• Huerre, P., y Monkewitz, P. A. (1990). Local and global instabilities in spatially developing flows. Annual Review of Fluid Mechanics, 22, 473–537. https://doi.org/10.1146/annurev.fl.22.010190.002353
• Pope, S. B. (2000). Turbulent Flows. Cambridge University Press. https://doi.org/10.1017/CBO9781316179475
• Trefethen, L. N., Trefethen, A. E., Reddy, S. C., y Driscoll, T. A. (1993). Hydrodynamic stability without eigenvalues. Science, 261, 578–584. https://doi.org/10.1126/science.261.5121.578
• Ottino, J. M. (1990). Mixing, chaotic advection, and turbulence. Annual Review of Fluid Mechanics, 22, 207–254. https://doi.org/10.1146/annurev.fl.22.010190.001231