Me gustaría aprender cómo funciona el elemento Raviart-Thomas (RT). Con ese fin, me gustaría describir analíticamente cómo se ven las funciones básicas en el cuadrado de referencia. El objetivo aquí no es implementarlo yo mismo, sino simplemente obtener una comprensión intuitiva del elemento.
Estoy basando en gran medida este trabajo en los elementos triangulares discutidos aquí , tal vez extenderlo a cuadriláteros es un error en sí mismo.
Dicho esto, puedo definir las funciones básicas para el primer elemento RK RK0:
parai=1,…,4.
Las condiciones en son que:
donde es la unidad normal que se muestra a continuación, y x j es su coordenada.
Este es el cuadrado de referencia , por lo que esto conduce a un sistema de ecuaciones para cada función base. Para ϕ 1 esto es:
que se puede resolver para dar:
Las otras funciones básicas se pueden encontrar de manera similar.
Suponiendo que esto sea correcto, el siguiente paso es encontrar las funciones básicas para RK1. Aquí es donde me estoy poniendo un poco inseguro de mí mismo. Según el enlace de arriba, el espacio que nos interesa es:
Una base para sería { 1 , x , y }
Creo que esto significa que las funciones básicas de RK1 deberían tomar la forma:
Esto deja 10 incógnitas para cada función básica. Si aplicamos las mismas condiciones que en el caso RK0, a saber:
, donde n j es la unidad normal como se muestra a continuación:
Esto nos da 8 ecuaciones. Los otros 2, creo, se pueden encontrar en algunos momentos. No estoy realmente seguro de cómo exactamente. El enlace de arriba habla sobre la integración contra una base para , pero tengo problemas para entender qué significa eso. ¿Estoy en el camino correcto o me he perdido algo por completo aquí?