Decía el señor Swett Marden que

“Un guijarro en el lecho de un pobre arroyuelo puede mudar el curso de un río”.

Parece una exageración y sin duda lo es, pero sirve para que nos hagamos la siguiente pregunta:

¿Cómo de importante es el lecho marino para las olas en la superficie?

Ésta es la pregunta que tratamos de contestar en este artículo. El problema que queremos entender es, dados dos fluidos incompresibles en un medio poroso acotado, ¿cómo se comporta la interfase entre ambos y que diferencias presenta con el caso en el que el medio no esta acotado? Bueno, vamos a trasladar ese problema físico a ecuaciones en derivadas parciales. Tenemos una densidad que presenta dos valores según estemos por encima o por debajo de la interfase, que denotamos por ,

Que los fluidos sean incompresibles y se muevan en un medio poroso acotado quiere decir que el dominio espacial de los fluidos es 

y que la velocidad satisface la Ley de Darcy y la condición de incompresibilidad

Estas ecuaciones se puede trasladar a una única ecuación para la interfase:

Ahora que tenemos el problema cabe preguntarnos si el hecho de que el dominio sea y no cambia mucho la situación. Para hacernos una idea podemos hacer unas simulaciones numéricas preliminares. Para ello consideramos un dato inicial y lo hacemos evolucionar en el caso donde el medio tiene profundidad finita (caso acotado) y también en el caso en el que el medio tiene una profundidad infinita (caso no acotado). Por supuesto el resto de los parámetros físicos son los mismos en ambas evoluciones. Así observamos lo siguiente

(Si no ves bien las imágenes pincha en ellas para hacerlas más grandes)

Parece claro a la vista de estos resultados que el hecho de que el medio esté acotado o no es relevante para las olas.

Una vez que tenemos el problema propuesto tenemos que empezar a sacar teoremas :-P. Evitando tecnicismos lo primero que probamos es

1) (Existencia y unicidad) que si el fluido de arriba es más ligero que el que está abajo el problema tiene una solución.

1.b) (Existencia y unicidad 2) que si el fluido de arriba es más pesado que el de abajo pero la interfase inicial es analítica existe una solución.

2) (Efecto regularizante) que dicha solución se vuelve muy regular (analítica) para cualquier (compárese con la ecuación del calor aquí.)

De momento estos 3 teoremas son idénticos en su enunciado a los teoremas cuando la profundidad es infinita. ¿Sorprendido? Bueno, esto sólo quiere decir que para probar matemáticamente las diferencias que hemos visto en los vídeos y las imágenes anteriores tenemos que trabajar un poco más, así que sed pacientes y esperad a la siguiente entrada ;-).

Bueno, si os veis muy impacientes podéis leer (o, en su caso, releer) ésta, ésta y esta entrada.

–Referencias:

D. Córdoba, RGB, R.Orive, The confined Muskat problem: differences with the deep water regime.

–Nota: Con esta entrada participamos en el Carnaval de Matemáticas en su edición 26, organizado esta vez por ZTFNews.

–Nota 2: Con esta entrada participamos también en el Carnaval de Física en su edición XXXIV que organiza Hablando de Ciencia.