Problema de capa límite

De Laplace

(Diferencias entre revisiones)

última version al 13:02 15 ago 2009

Contenido

1 Enunciado
2 Generalizamos el problema
3 La solución ingenua
4 Solución interior
5 El matching
6 La solución uniforme
7 Algunas gráficas

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

$0.01y'' + y' + x y = 0\,$ $y(0) = 0\,$ $y(1) = 1\,$

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

$\epsilon y'' + y' + x y = 0\,$ $y(0) = 0\,$

y (1) = 1

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

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

$y' + xy = 0\,$

cuya solución es de la forma

$y = a \exp\left(-\frac{x^2}{2}\right)$

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,

$1 = a \exp\left(-\frac{1}{2}\right)$ $\Rightarrow$ $a = \sqrt{\mathrm{e}}$

así que aproximamos la solución por

$y = \sqrt{\mathrm{e}}\exp\left(-\frac{x^2}{2}\right)$

Pero es claro que esta solución no vale cerca de $x = 0$ , pues ahí vale $\sqrt{\mathrm{e}}$ 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 $x \ll 1$ , como

$X = \frac{x}{\delta}$

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

$\frac{\varepsilon}{\delta^2}\frac{\mathrm{d}^2Y}{\mathrm{d}X^2} + \frac{1}{\delta}\frac{\mathrm{d}Y}{\mathrm{d}X} + \delta X Y = 0$

el valor de $δ$ lo sacamos de que al menos dos de los términos sean del mismo orden, lo que en este caso nos da

$\delta = \varepsilon$

y nos deja con

$\frac{\mathrm{d}^2Y}{\mathrm{d}X^2} + \frac{\mathrm{d}Y}{\mathrm{d}X} + \varepsilon^2 X Y = 0$

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

$\frac{\mathrm{d}^2Y}{\mathrm{d}X^2} + \frac{\mathrm{d}Y}{\mathrm{d}X} = 0$ $\Rightarrow$ $Y = A + B \exp\left(-X\right)$

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?

5 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\to 0} y = \lim_{X\to\infty} Y$

lo que nos da

$\sqrt{\mathrm{e}} = A\,$

y por tanto la solución interior es

$Y = \sqrt{\mathrm{e}}\left(1-\exp(-X)\right) = \sqrt{\mathrm{e}}\left(1-\exp\left(-\frac{x}{\varepsilon}\right)\right)$

6 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 = \sqrt{\mathrm{e}}\exp\left(-\frac{x^2}{2}\right)$

y una que vale cuando $x$ es muy pequeño (de orden $\varepsilon$ o menor)

$Y = \sqrt{\mathrm{e}}\left(1-\exp\left(-\frac{x}{\varepsilon}\right)\right)$

¿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 = \sqrt{\mathrm{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 = \sqrt{\mathrm{e}}\left(\exp\left(-\frac{x^2}{2}\right) - \exp\left(-\frac{x}{\varepsilon}\right)\right)$

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 $\varepsilon$ en todo el intervalo (de hecho, el error máximo es menor que el 1% en todo el intervalo.

7 Algunas gráficas

La solución calculada numéricamente

Las soluciones interior y exterior

La exacta y las dos aproximadas

La exacta y la solución uniforme

Error absoluto en el intervalo

Error relativo en todo el intervalo

@@ Línea 24: / Línea 24: @@
 ¿Cuál de las dos? Es fácil ver que si imponemos <math>y(0) = 0</math> nos quedamos sin nada, así que esa es la problemática. Sí podemos imponer la otra,
-<center><math>1 = a \exp\left(-\frac{1}{2}\right)</math>{{tose}}<math>a = \sqrt{e}</math></center>
+<center><math>1 = a \exp\left(-\frac{1}{2}\right)</math>{{tose}}<math>a = \sqrt{\mathrm{e}}</math></center>
 así que aproximamos la solución por
-<center><math>y = \sqrt(e)\exp\left(-\frac{x^2}{2}\right)</math></center>
+<center><math>y = \sqrt{\mathrm{e}}\exp\left(-\frac{x^2}{2}\right)</math></center>
-Pero es claro que esta solución no vale cerca de <math>x=0</math>, pues ahí vale <math>\sqrt{e}</math> y debería valer 0. Cerca de <math>x=0</math> la solución debe variar muy rápidamente para poder satisfacer la condición de contorno
+Pero es claro que esta solución no vale cerca de <math>x=0</math>, pues ahí vale <math>\sqrt{\mathrm{e}}</math> y debería valer 0. Cerca de <math>x=0</math> 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.
@@ Línea 39: / Línea 39: @@
 Veamos como podemos estudiar el comportamiento cerca de <math>x=0</math>. Necesitamos una lupa, que en este caso es un reescalado. Definamos una nueva variable <math>X</math>, que es finita en la región <math>x \ll 1</math>, como
-<math>X = \frac{x}{\delta}</math>
+<center><math>X = \frac{x}{\delta}</math></center>
-con delta otro parámetro pequeño que hay que determinar. Si escribimos la ecuación diferencial en términos de X queda
+con <math>\delta</math> otro parámetro pequeño que hay que determinar. Si escribimos la ecuación diferencial en términos de <math>X</math> queda
- eps/delta^2 Y'' + Y'/delta + delta X Y = 0
+<center><math> \frac{\varepsilon}{\delta^2}\frac{\mathrm{d}^2Y}{\mathrm{d}X^2} + \frac{1}{\delta}\frac{\mathrm{d}Y}{\mathrm{d}X} + \delta X Y = 0</math></center>
-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
+el valor de <math>\delta</math> lo sacamos de que al menos dos de los términos sean del mismo orden, lo que en este caso nos da
- Y'' + Y' + eps^2 X Y = 0
+<center><math>\delta = \varepsilon</math></center>
+y nos deja con
+<center><math>\frac{\mathrm{d}^2Y}{\mathrm{d}X^2} + \frac{\mathrm{d}Y}{\mathrm{d}X} + \varepsilon^2 X Y = 0</math></center>
 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)
+<center><math>\frac{\mathrm{d}^2Y}{\mathrm{d}X^2} + \frac{\mathrm{d}Y}{\mathrm{d}X}  = 0</math>{{tose}}<math>Y = A + B \exp\left(-X\right)</math></center>
-Ahora sí podemos imponer la condición en x=0 y escribir la solución interior como
+Ahora sí podemos imponer la condición en <math>x=0</math> y escribir la solución interior como
- Y = A(1 - exp(-X))
+<center><math>Y = A(1 - \exp(-X))\,</math></center>
 Pero ¿cuánto vale A?
-) El 'matching'
+==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
+<center><math>\lim_{x\to 0} y = \lim_{X\to\infty} Y</math></center>
 lo que nos da
- rq(e) = A
+<center><math>\sqrt{\mathrm{e}} = A\,</math></center>
 y por tanto la solución interior es
- Y = rq(e)(1-exp(-X)) = rq(e)(1-exp(-X/eps))
+<center><math>Y = \sqrt{\mathrm{e}}\left(1-\exp(-X)\right) = \sqrt{\mathrm{e}}\left(1-\exp\left(-\frac{x}{\varepsilon}\right)\right)</math></center>
+==La solución uniforme==
+Tenemos entonces dos aproximaciones para la solución de la ecuación diferencial. Una que vale siempre que <math>x</math> no sea muy pequeño
+<center><math>y = \sqrt{\mathrm{e}}\exp\left(-\frac{x^2}{2}\right)</math></center>
+y una que vale cuando <math>x</math> es muy pequeño (de orden <math>\varepsilon</math> o menor)
+<center><math>Y = \sqrt{\mathrm{e}}\left(1-\exp\left(-\frac{x}{\varepsilon}\right)\right)</math></center>
+¿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 <math>L = \sqrt{\mathrm{e}}</math> y para <math>x</math> finito la exterior vale lo mismo. Esto nos permite construir una solución uniforme sumando ambas y restando el límite común
+<center><math>y_u(x) = y + Y - L = \sqrt{\mathrm{e}}\left(\exp\left(-\frac{x^2}{2}\right) - \exp\left(-\frac{x}{\varepsilon}\right)\right)</math></center>
+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 <math>\varepsilon</math> en todo el intervalo (de hecho, el error máximo es menor que el 1% en todo el intervalo.
+==Algunas gráficas==
+* La solución calculada numéricamente
+<center>[[Imagen:gexacta.gif]]</center>
+* Las soluciones interior y exterior
-) La solución uniforme
+<center>[[Imagen:gio.gif]]</center>
-Tenemos entonces dos aproximaciones para la solución de la ecuación diferencial. Una que vale siempre que x no sea muy pequeño
+* La exacta y las dos aproximadas
- y = rq(e)exp(-x^2/2)
+<center>[[Imagen:geio.gif]]</center>
-y una que vale cuando x es muy pequeño
+* La exacta y la solución uniforme
- Y = rq(e)(1-exp(-x/eps))
+<center>[[Imagen:geu.gif]]</center>
-¿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
+* Error absoluto en el intervalo
- y_u(x) = y + Y - L = rq(e)(exp(-x^2/2) - exp(-x/eps))
+<center>[[Imagen:eabs.gif]]</center>
-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.
+* Error relativo en todo el intervalo
-Algunas gráficas:
+<center>[[Imagen:erel.gif]]</center>

Problema de capa límite

De Laplace

última version al 13:02 15 ago 2009

Contenido

1 Enunciado

2 Generalizamos el problema

3 La solución ingenua

4 Solución interior

5 El matching

6 La solución uniforme

7 Algunas gráficas

Herramientas:

Buscar

Vistas

Herramientas

Herramientas personales

Navegación

SEARCH

TOOLBOX

LANGUAGES