PDE de tipo Poisson con solución analítica

Una manera sencilla de tener una ecuación de Poisson en $latex 3D$ de la que conocer su solución analítica es la siguiente. Para empezar, consideramos una función:

$latex u(x,y,z)$

a la que le aplicamos el operador $latex Delta$ y obtendremos otra función:

$latex s(x,y,z)$.

Ya tenemos $latex Delta u = s$, es decir,

$latex frac{partial^2}{partial x^2}u(x,y,z) + frac{partial^2}{partial y^2}u(x,y,z) + frac{partial^2}{partial z^2}u(x,y,z) = s(x,y,z)$

Para las condiciones de contorno es tan sencillo como considerar el domino:

$latex [a,b] times [c,d] times [e,f]$

y ver cuanto vale $latex u$ en cada uno de los extremos, de manera que obtenemos:

$latex u(a,y,z) = g_a(y,z), u(b,y,z) = g_b(y,z)$,

$latex u(x,c,z) = g_c(x,z), u(x,d,z) = g_d(x,z)$,

$latex u(x,y,e) = g_e(x,y), u(x,y,f) = g_f(x,y)$.

Por ejemplo, si consideramos $latex u(x,y,z)=x^2+y^2+z^2$, entonces:

$latex nabla cdot nabla u = (frac{partial}{partial x},frac{partial}{partial y},frac{partial}{partial z}) cdot (u_x,u_y,u_z)$

de manera que:

$latex Delta u = frac{partial}{partial x}2x + frac{partial}{partial y}2y + frac{partial}{partial z}2z = 6$

y tenemos la ecuación de Poisson $latex Delta u = 6$. Las condiciones de contorno, en $latex Omega = [0,1]^3$ quedan:

$latex u(0,y,z) = g_{xm}(y,z) = y^2+z^2$

$latex u(1,y,z) = g_{xM}(y,z) = y^2+z^2 + 1$

$latex u(x,0,z) = g_{ym}(x,z) = x^2+z^2$

$latex u(x,1,z) = g_{yM}(x,z) = x^2+z^2 + 1$

$latex u(x,y,0) = g_{zm}(x,y) = x^2+y^2$

$latex u(x,y,1) = g_{zM}(x,y) = x^2+y^2 + 1$

Resumiendo, la solución de $latex Delta u = 6$ siendo las funciones anteriores los valores de $latex u$ en $latex partial Omega$ es:

$latex u = x^2+y^2 + z^2$.

Otro ejemplo concreto para el caso de $latex u=0$ en $latex partial Omega$ siendo

$latex Omega = { (x,y,z): 0<x<1, 0<y<1,0<z<1}$ el cubo unidad.

Tomamos $latex u(x,y,z) = (x^4-x^2)(y^4-y^2)(z^4-z^2)$. Es sencillo comprobar que $latex u=0$ en $latex partial Omega$ (p.e. $latex u(1,y,z)=(1-1)(y^4-y^2)(z^4-z^2) = 0$). ¿Cuanto vale $latex Delta u$ en este caso?

$latex s(x,y,z) = Delta [(x^4-x^2)(y^4-y^2)(z^4-z^2)] = $

$latex = nabla cdot [(4x^3-2x)(y^4-y^2)(z^4-z^2),$

$latex ,(x^4-x^2)(4y^3-2y)(z^4-z^2),$

$latex (x^4-x^2)(y^4-y^2)(4z^3-2z)] = $

$latex = 2[(6x^2-1)(y^4-y^2)(z^4-z^2) +$

$latex + (x^4-x^2)(6y^2-1)(z^4-z^2) + $

$latex + (x^4-x^2)(y^4-y^2)(6z^2-1)]$.

De manera que la solución en el cubo unidad de la ecuación de Poisson

$latex u_{xx} + u_{yy} + u_{zz} = $

$latex = 2[(6x^2-1)(y^4-y^2)(z^4-z^2) +$

$latex + (x^4-x^2)(6y^2-1)(z^4-z^2) + $

$latex + (x^4-x^2)(y^4-y^2)(6z^2-1)]$

con condiciones homogeneas de tipo Dirichlet en la frontera tiene como solución:

$latex u(x,y,z) = (x^4-x^2)(y^4-y^2)(z^4-z^2)$

Dejar un comentario

Tu dirección de correo electrónico no será publicada. Los campos obligatorios están marcados con *


¡IMPORTANTE! Responde a la pregunta: ¿Cuál es el valor de 9 3 ?