5 votos

Derivada normal expresada en coordenadas polares

Una derivada normal en un contorno dado $\Gamma$ de una función $u$ expresado en coordenadas cartesianas $u=u(x,y)$ se puede calcular como: $$\frac{\partial u}{\partial n} = \nabla u \cdot n = \frac{\partial u}{\partial x} n_x + \frac{\partial u}{\partial y} n_y $$ donde los símbolos están representados en la figura.

enter image description here

Cuando un segmento de una curva está parametrizado por la variable normalizada $\xi \in [-1,1]$ podemos describir su geometría como:

$x=x(\xi)= \sum_i x_i N_i(\xi)$ y $y=y(\xi)= \sum_i y_i N_i(\xi)$

donde $(x_i, y_i)$ son puntos nodales, y $N_i$ - polinomios de interpolación, por ejemplo, los polinomios de Lagrange. Esto conduce a las derivadas:

$\frac{dx}{d\xi}= \sum_i x_i \frac{dN_i}{d\xi}$ y $\frac{dy}{d\xi}= \sum_i y_i \frac{dN_i}{d\xi}$

El jacobiano de la transformación entre coordenadas locales y globales puede calcularse como

$$d\Gamma = \sqrt{dx^2+dy^2} $$

y las componentes del vector normal son iguales

$n_x = \cos \alpha = \frac{dy}{d\Gamma}$

$n_y = - \sin \alpha = -\frac{dx}{d\Gamma}$

De este modo se pueden describir las componentes del vector normal en términos de coordenadas cartesianas y, en consecuencia, en términos de parámetro $\xi$ .

Ahora bien, ¿qué pasa si la función $u$ y una curva $\Gamma$ se expresan en coordenadas polares $r = r(\theta)$ , $\theta = [0,2\pi]$ , $u=u(r,\theta)$ ? Recientemente, un documento que he encontrado proporciona la siguiente fórmula: $$\frac{\partial u}{\partial n} = \frac{\partial u}{\partial r} - \frac{r'}{r^2} \frac{\partial u}{\partial \theta}.$$

Según la teoría la derivada normal en coordenadas polares viene dada por: $$\frac{\partial u}{\partial n} = \nabla u \cdot n = \frac{\partial u}{\partial r} n_r + \frac{1}{r} \frac{\partial u}{\partial \theta} n_\theta $$ y lo que necesito es expresar las componentes radial y angular del vector normal en términos de coordenadas polares, o de hecho en términos de parámetro $\theta$ porque $r=r(\theta)$ . Teniendo en cuenta que $x=r(\theta) \cos \theta$ , $y=r(\theta) \sin \theta$ y $d\Gamma = \sqrt{dx^2+dy^2} $ se puede obtener que $$d\Gamma = \sqrt{r'^2 + r^2} $$ pero no encuentro ninguna relación trigonométrica entre $dr$ , $d\theta$ y $n_r$ , $n_\theta$ como se puede encontrar para $dx$ , $dy$ , $n_x$ , $n_y$ . Así que mi pregunta es: ¿Podría alguien proporcionar fórmulas para $n_r$ y $n_\theta$ en términos de coordenadas polares?

EDITAR : He encontrado otro documento que contiene la fórmula que estaba buscando. Es diferente de la fórmula del artículo anteriormente citado del mismo autor. Lo que es más, es un poco diferente de lo que Ted Schifrin derivado.

$$\frac{\partial u}{\partial n} = \frac{r}{\sqrt{r'^2 + r^2}} \left[ \frac{\partial u}{\partial r} - \frac{r'}{r^2} \frac{\partial u}{\partial \theta} \right] $$

El apéndice del documento ofrece la derivación de la fórmula. Ahora todo me parece bien.

0 votos

¿Puede proporcionar un enlace a la fuente?

0 votos

Aquí está . Ecuación (3)

0 votos

Entonces, has editado parcialmente. Tu primera ecuación debería estar escrita de la siguiente manera: $$\nabla u = \frac{\partial u}{\partial r}\hat r + \frac1r\frac{\partial u}{\partial\theta}\hat\theta$$ (o, si lo prefiere, puede escribir $e_r,e_\theta$ para los vectores unitarios). La fórmula de la derivada normal es algo totalmente diferente. No se trata de "esto se puede derivar como ...".

2voto

Ted Shifrin Puntos 33487

Se supone que el contorno en cuestión viene dado por $r=r(\theta)$ .

EDITAR : Vale, he intentado descifrar lo que está pasando aquí. Supongo que el artículo que citas se refiere a la derivada normal a lo largo de una curva de nivel (contorno) de $u$ . La fórmula no parece ser correcta con mi suposición de que el autor está parametrizando la curva por $r=r(\theta)$ . Con esta suposición, obtengo la fórmula $$\frac{\partial u}{\partial n} = \frac{r\frac{\partial u}{\partial r}-r'\frac{\partial u}{\partial\theta}}{\sqrt{r^2 + r'^2}}=\frac{\frac{\partial u}{\partial r}-\frac{r'}r\frac{\partial u}{\partial\theta}}{\sqrt{1 + \big(\frac{r'}r\big)^2}}.$$

0 votos

La fórmula que corresponde a (2) en el artículo citado también tiene un signo "-" delante del $(r'/r^2) (\partial u/\partial \theta)$ plazo.

0 votos

Efectivamente, hay este signo menos en el papel y había este menos en mi primer post. He corregido la fórmula. Espero que los editores no la "revisen" de nuevo :-)

0 votos

@ Ted Shifrin Durante los días de Navidad tuve mi tiempo para pensar en el problema una vez más. He editado mi primer post. Desgraciadamente no consigo obtener ninguna de las soluciones dadas. ¿Puedo pedirte que me des más pistas?

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X