Problema de capa límite

De Laplace

1 Enunciado

Obtener una solución aproximada, con un error máximo de un 1% en todo el intervalo (0,1) para la ecuación diferencial

                

2 Generalizamos el problema

En vez de hacerlo para el caso particular de 1% de error y con un coeficiente 0.01, vamos a suponer un parámetro pequeño eps y estudiamos la ecuación diferencial

                y(1) = 1

y queremos hallar una solución aproximada con un error de orden .

3 La solución ingenua

Como primera aproximación, si toleramos un error de orden eps, despreciamos el primer término y nos quedamos con

cuya solución es de la forma

pero rápidamente tropezamos con un problema: tenemos una sola constante de integración, a, y dos condiciones de contorno. No podemos satisfacer ambas a la vez. Como mucho una.

¿Cuál de las dos? Es fácil ver que si imponemos y(0) = 0 nos quedamos sin nada, así que esa es la problemática. Sí podemos imponer la otra,

      

así que aproximamos la solución por

Pero es claro que esta solución no vale cerca de x = 0, pues ahí vale y debería valer 0. Cerca de x = 0 la solución debe variar muy rápidamente para poder satisfacer la condición de contorno

Estamos en un caso de perturbaciones singulares. Si existe una zona en la que la solución varía muy rápidamente, como aquí cerca de x=0, la segunda derivada será muy grande y el término que hemos despreciado no será pequeño.

Lo que hemos hallado es realmente la solución exterior, que vale fuera de esta zona de cambio rápido (la capa límite), pero no dentro de ella.

4 Solución interior

Veamos como podemos estudiar el comportamiento cerca de x = 0. Necesitamos una lupa, que en este caso es un reescalado. Definamos una nueva variable X, que es finita en la región , como

con delta otro parámetro pequeño que hay que determinar. Si escribimos la ecuación diferencial en términos de X queda

eps/delta^2 Y + Y'/delta + delta X Y = 0

el valor de delta lo sacamos de que al menos dos de los términos sean del mismo orden, lo que en este caso nos da delta = eps y nos deja con

Y + Y' + eps^2 X Y = 0

despreciando ahora el último término, que es mucho más pequeño, nos quedamos con la ecuación diferencial

Y + Y' = 0      ---> Y = A + B exp(-X)

Ahora sí podemos imponer la condición en x=0 y escribir la solución interior como

Y = A(1 - exp(-X))

Pero ¿cuánto vale A?

4) El 'matching'

Para poder hallar A debemos empalmar la solución interior con la exterior. Esto se hace imponiendo que sobre una escala intermedia entre x y X (mucho menor que x pero mucho mayor que X) ambas soluciones coincidan. En su versión más sencilla, se impone que el comportamiento, para x pequeño de la solución exterior coincida con el comportamiento, para X grande, de la solución interior, esto es,

lim_(x->0) y = lim_(X->oo) Y

lo que nos da

rq(e) = A

y por tanto la solución interior es

Y = rq(e)(1-exp(-X)) = rq(e)(1-exp(-X/eps))

5) La solución uniforme

Tenemos entonces dos aproximaciones para la solución de la ecuación diferencial. Una que vale siempre que x no sea muy pequeño

y = rq(e)exp(-x^2/2)

y una que vale cuando x es muy pequeño

Y = rq(e)(1-exp(-x/eps))

¿Cómo podemos conseguir una solución aproximada para todo el intervalo? Recordemos el límite común de ambas funciones. Esto quiere decir que para x pequeño la solución exterior vale L = rq(e) y para x finito la exterior vale lo mismo. Esto nos permite construir una solución uniforme sumando ambas y restando el límite común

y_u(x) = y + Y - L = rq(e)(exp(-x^2/2) - exp(-x/eps))

Esta es la solución que andábamos buscando. No cumple la ecuación diferencial en ningún punto, pero se diferencia de la solución en un término del orden de eps en todo el intervalo (de hecho, el error máximo es menor que el 1% en todo el intervalo.

Algunas gráficas: