[Steven_Math]

Buenas noches, he estado resolviendo la siguiente ecuación pero no he podido:

 
\(
\begin{cases}{-v\Delta u+\nabla p=f}&\text{en}& \Omega \\
\text{div} u=0 & \text{en}& \Omega\\u=0 & \text{en}& \partial\Omega\end{cases}
 \)

y
                               \( \displaystyle\int_{\Omega}p(r,z)r dr dz=0 \)

Debo encontrar un ejemplo de funciones \( u(r,z) \) y \( p(r,z) \) que satisfagan las ecuaciones anteriores, les agradezco su gran ayuda.

[mathtruco]

Hola Steven_Math. Es mucho más fácil de lo que crees. Evitando el ejemplo trivial (donde todas las funciones son la nula):

Considera cualquier presión con integral nula. Por ejemplo, elige cualquier función (derivable) \( \tilde p(x) \) y define

    \( p(x)=\tilde p(x)-\dfrac{1}{|\Omega|}\displaystyle\int_\Omega \tilde p(x)dx \).

Obviamente \( p \) definida así tendrá integral nula. Puedes elegir un polinomio por simplicidad, aunque da igual.

Ahora elige una función \( u \) con divergencia nula y que sea cero en la frontera de \( \Omega \).  Podrías elegir velocidad nula, pero sería demasiado trivial. Si quieres pensar algo menos trivial, puedes pensar en un polinomio que, que por supuesto, dependerá del dominio  \( \Omega \).

Por último, teniendo \( u \) y \( p \) los reemplazas en tu primera ecuación y obtendrás \( f \).

[mathtruco]

Una solución menos trivial para la velocidad:

    \( \begin{cases}
u_1(x,y,z)=x^2y(x-1)^2(y-1)(2y-1)
\\
u_2(x,y,z):=-u_1(y,x)
\\
u_3(x,y,z)=0
\end{cases}
 \)

Si buscas cualquier paper de análisis numérico que trabaje con la ecuación de Stokes tendrá que darse como primer ejemplo un problema con solución exacta conocida donde velocidad tenga divergencia nula. En algunos de ellos también buscarán que la condición Dirichlet sea nula.

P.D. Olvidé mencionar que el dominio es el cubo unitario: \( \Omega=(0,1)^3 \).

[Steven_Math]

Muchas gracias Mathtruco,
¿podría considerar las funciones:

   \( p(x,y)=(x-1/2)(y-1/2) \  \ \text{y} \   \ u(x,y)=\begin{pmatrix} -256x^{2}(x-1)^2y(y-1)(2y-1)  \\  256y^2(y-1)^2x(x-1)(2x-1) \end{pmatrix} , \ \Omega=(0,1)^2 \)?

[mathtruco]

Claro, es un ejemplo que aparecen en muchos papers (en particular uno mío  )

[Steven_Math]

Me parece interesante que tengas papers en este tema.
Yo estoy afianzando en este tema, sobre todo para implementar el método  de elementos  finitos  en Feenic's.
Si no es mucha molestia, te agradezco que me envíes  documentos  que puedan servir para ir avanzando en estos temas.
Mi correo electrónico  es: [email protected].

Saludos.

[mathtruco]

No creo poder ayudarte con material. Alguna vez intenté usar Fenics, pero tuve líos con versiones de python y documentación no actualizada de su página. Pretendo volver al intentarlo cuando lo requiera. Lo que he hecho es con códigos hechos en casa.

Pero cuando tengas dudas específicas las puedes postear en el foro.

Buena suerte.