Operadores diferenciales - OCW UPM
Operadores diferenciales - OCW UPM
Operadores diferenciales - OCW UPM
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Dpto. Física y Mecánica<br />
<strong>Operadores</strong> <strong>diferenciales</strong>
Se denominan líneas coordenadas de un espacio<br />
euclídeo tridimensional a aquellas que se obtienen<br />
partiendo un punto dado P de coordenadas (q (q1, 1, q q2, 2, q q3), 3),<br />
variando una de ellas y manteniendo fijas las otras dos.<br />
Si las líneas coordenadas son ortogonales en un punto<br />
las coordenadas se denominan coordenadas<br />
ortogonales. Por ejemplo, las coordenadas cartesianas<br />
(x,y,z) son ortogonales, al igual que los son las<br />
coordenadas d d cilíndricas ilí d i ( (r, ϕ,z) ) y llas coordenadas d d<br />
esféricas (r, θ,ϕ).<br />
<strong>OCW</strong>‐<strong>UPM</strong>
En matemáticas un operador diferencial es un operador lineal<br />
definido como una función del operado de diferenciación.<br />
El uso más común del operador diferencial es realizar la<br />
derivada en sí misma.<br />
Considerando que la variable es x, la notación más común de<br />
éste operador es<br />
d<br />
dx<br />
Las sucesivas derivadas se expresan por<br />
2 3 n<br />
d d d<br />
, ,...<br />
2 3 n<br />
dx dx dx<br />
<strong>OCW</strong>‐<strong>UPM</strong>
Cuando se consideran derivadas respecto a variables<br />
diferentes las derivadas parciales pueden escribirse como:<br />
∂ ∂ ∂<br />
, ,<br />
∂x ∂y ∂z<br />
Las sucesivas derivadas se expresan por<br />
2 2 2<br />
∂ ∂ ∂<br />
, ,<br />
∂x ∂y ∂z<br />
2 2 2<br />
......................<br />
∂ ∂ ∂<br />
, ,<br />
∂x ∂y ∂z<br />
n n n<br />
n n n<br />
<strong>OCW</strong>‐<strong>UPM</strong>
Y las derivadas parciales cruzadas<br />
2 2 2 2 2 2<br />
∂ ∂ ∂ ∂ ∂ ∂<br />
, , , , ,<br />
∂∂ x y ∂∂ y x ∂∂ x z ∂∂ z x ∂∂ y z ∂∂ z y<br />
<strong>OCW</strong>‐<strong>UPM</strong>
Líneas í coordenadas<br />
Se denominan líneas coordenadas de un espacio euclídeo<br />
tridimensional a aquellas que se obtienen partiendo un punto<br />
dado P de coordenadas (q 1,q 2,q 3), variando una de ellas y<br />
manteniendo fijas las otras dos.<br />
Si las líneas coordenadas son ortogonales en un punto las<br />
coordenadas se denominan coordenadas ortogonales ortogonales. Por<br />
ejemplo, las coordenadas cartesianas (x,y,z) son ortogonales,<br />
al igual que los son las coordenadas cilíndricas (r, ϕ,z) y las<br />
coordenadas esféricas (r, θ,ϕ).
Líneas í coordenadas<br />
TTomando d tres vectores tangentes a cada d lí línea en un punto,<br />
se obtienen tres vectores ortogonales entre sí, pero no<br />
unitarios unitarios.<br />
�<br />
e<br />
i<br />
�<br />
∂r<br />
=<br />
∂ q<br />
Para obtener un sistema ortonormal, se divide cada vector<br />
por su módulo ód l<br />
hi( q1, q2, q3) =<br />
�<br />
ei<br />
�<br />
∂r<br />
=<br />
∂<br />
q<br />
enominado factor de escala<br />
i<br />
i
Líneas í coordenadas<br />
Dado un conjunto de<br />
coordenadas (q1, 1 q2, 2 q3) 3 sobre<br />
un espacio euclídeo, cuyas<br />
líneas son ortogonales entre sí,<br />
se puede d construir t i una bbase<br />
ortonormal en cada punto a<br />
partir de los vectores tangentes<br />
en cada una de las líneas<br />
coordenadas.<br />
q3<br />
u1 �<br />
u3 �<br />
u �<br />
u<br />
2<br />
q1<br />
q2
Factor de escala<br />
En el cálculo vectorial se definen unas cantidades cantidades,<br />
denominadas factor<br />
proporción entre lo<br />
de escala<br />
que varía<br />
hi, que representan<br />
una coordenada y<br />
la<br />
el<br />
desplazamiento que produce esta variación.<br />
∂r �<br />
� ∂r<br />
hi( q1, q2, q3) = ei<br />
=<br />
∂q<br />
A partir ti dde lla expresión ió siendo i d<br />
�<br />
⎛ ⎞ ⎛ ⎞ ⎛ ⎞⎛ ⎞<br />
⎜ ⎟ ⎜ ⎟ ⎜ ⎟⎜ ⎟<br />
⎝ ⎠ ⎝ ⎠ ⎝ ⎠⎝ ⎠<br />
2 2<br />
2 ∂r<br />
� dx � dy � dz � dx � dy � dz � dx � dy � dz<br />
( h i ) = = e x + e y + e z = e x + e y + e z · e x + e y + e z<br />
∂ ∂qi<br />
ddqi ddqi ddqi dq d i dq d i dq d i dq d i dq d i dq d i<br />
i
Factor de escala<br />
( h )<br />
i<br />
2<br />
2 2 2<br />
⎛ ∂x ⎞ ⎛ ∂y ⎞ ⎛ ∂z<br />
⎞<br />
= ⎜ ⎟ + ⎜ ⎟ +<br />
⎜ ⎟<br />
⎝∂qi ⎠ ⎝∂qi ⎠ ⎝∂qi ⎠<br />
dde donde d d el l factor f t de d escala l se expresa mediante di t la l expresión ió<br />
h<br />
i<br />
2 2 2<br />
⎛ ∂x ⎞ ⎛ ∂y ⎞ ⎛ ∂z<br />
⎞<br />
= ⎜ ⎟ + ⎜ ⎟ + ⎜ ⎟<br />
⎝∂qi ⎠ ⎝∂qi ⎠ ⎝∂qi ⎠
Factor de escala<br />
Considerando un sistema de coordenadas cartesianas (x (x,y,z) y z) y<br />
un sistema de coordenadas curvilíneas (q1, q2, q3). Las coordenadas cartesianas se pueden expresar en función<br />
de las coordenadas curvilíneas mediante relaciones de la<br />
forma<br />
x=x(q 1 q 2 q 3)<br />
x x(q 1, q 2, q 3)<br />
y=y(q 1, q 2, q 3)<br />
z=z(q 1, q 2, q 3)
Factor de escala<br />
Conociendo estas relaciones los factores de escala (h ( 1, 1, h2, 2, h3) 3)<br />
se determinan mediante las expresiones<br />
h<br />
h<br />
2<br />
1<br />
h<br />
2 2 2<br />
⎛ ∂x ⎞ ⎛ ∂y ⎞ ⎛ ∂z<br />
⎞<br />
= ⎜ ⎟ + ⎜ ⎟ + ⎜ ⎟<br />
⎝∂q1⎠ ⎝∂q1⎠ ⎝∂q1⎠ 2 2 2<br />
⎛ ∂x ⎞ ⎛ ∂y y ⎞ ⎛ ∂z⎞<br />
= ⎜ ⎟ + ⎜ ⎟ + ⎜ ⎟<br />
⎝∂q2 ⎠ ⎝∂q2 ⎠ ⎝∂q2 ⎠<br />
3<br />
2 2 2<br />
⎛ ∂x⎞ ⎛ ∂y ⎞ ⎛ ∂z<br />
⎞<br />
= ⎜ ⎟ + ⎜ ⎟ + ⎜ ⎟<br />
⎝ ∂ q q3 ⎠ ⎝ ∂ q q3 ⎠ ⎝ ∂ q q3<br />
⎠
Factor de escala en coordenadas cartesianas<br />
Cuando el sistema de coordenadas es el de coordenadas<br />
cartesianas ( q1, q2, q3) = ( x, y, z)<br />
, los vectores unitarios definidos<br />
� � � � � �<br />
sobre cada una de las direcciones son ( u , u , u ) = ( i, j, k)<br />
1 2 3<br />
La relación entre coordenadas es (x=x, y=y, z=z), y los<br />
ffactores t dde escala l son<br />
2 2 2<br />
⎛∂x⎞ ⎛∂y⎞ ⎛∂z⎞ h 1 = h x = ⎜ ⎟ + ⎜ ⎟ + ⎜ ⎟ = 1<br />
⎝ ⎝∂x ⎠ ⎝ ⎝∂x ⎠ ⎝ ⎝∂x ⎠<br />
2 2 2<br />
⎛∂x⎞ ⎛∂y⎞ ⎛∂z⎞ h h2 = h hy<br />
= ⎜ ⎟ + ⎜ ⎟ + ⎜ ⎟ = 1<br />
⎝∂y⎠ ⎝∂y⎠ ⎝∂y⎠ 2 2 2<br />
⎛∂x⎞ ⎛∂ y ⎞ ⎛∂z⎞ h h3= h hz=<br />
⎜ ⎟ + ⎜ ⎟ + ⎜ ⎟ = 1<br />
⎝∂z⎠ ⎝∂z ⎠ ⎝∂z⎠
Factor de escala en coordenadas cilíndricas<br />
Cuando el sistema de coordenadas es de coordenadas<br />
cilíndricas ( q 1 1, q 2 2, q 3 3)<br />
) = ( r, ϕϕ<br />
, z)<br />
) ,los vectores unitarios definidos sobre<br />
cada una de las direcciones son � � � � � �<br />
( u ,la relación<br />
1, u2, u3) = ( ur, uϕ, k)<br />
entre coordenadas es<br />
x= rcosϕ y = rsenϕ<br />
z =<br />
z
Factor de escala en coordenadas cilíndricas<br />
LLos ffactores t dde escala l son<br />
2 2 2<br />
⎛ ∂ x ⎞ ⎛ ∂ y ⎞ ⎛ ∂ z ⎞<br />
2 2<br />
h1= hr= ⎜ ⎟ + ⎜ ⎟ + ⎜ ⎟ = ( cosϕ) + ( senϕ)<br />
= 1<br />
⎝∂r ⎠ ⎝∂r ⎠ ⎝∂r⎠ 2 2 2<br />
⎛ ∂x ⎞ ⎛ ∂y ⎞ ⎛ ∂z<br />
⎞<br />
2 2<br />
h2 = hϕ = ⎜ ⎟ + ⎜ ⎟ + ⎜ ⎟ = ( − rsenϕ) + ( r cosϕ)<br />
= r<br />
⎝ ⎝∂ϕϕ ⎠ ⎝ ⎝∂ϕϕ ⎠ ⎝ ⎝∂ϕϕ ⎠<br />
2 2 2<br />
⎛ ∂ x ⎞ ⎛ ∂ y ⎞ ⎛ ∂ z ⎞<br />
h3= hz=<br />
⎜ ⎟ + ⎜ ⎟ + ⎜ ⎟ = 1<br />
⎝∂z⎠ ⎝∂z ⎠ ⎝∂z⎠
Factor de escala en coordenadas esféricas<br />
Cuando el sistema de coordenadas es de coordenadas<br />
esféricas ( q1, q2, q3) = ( r, ϕ,<br />
z)<br />
los vectores unitarios definidos sobre<br />
� � � � � �<br />
cada una de las direcciones son ( u1, u2, u3) = ( ur, uθ, u,<br />
ϕ ) la relación entre<br />
coordenadas cartesianas y esféricas es<br />
x = rsenθ rsenθ cos cosϕϕ<br />
y = rsenθ senϕ<br />
z =<br />
r<br />
cosθ
Factor de escala en coordenadas esféricas<br />
Los factores de escala son<br />
2 2 2<br />
⎛∂x⎞ ⎛∂y⎞ ⎛∂z⎞ 2 2 2<br />
h h1 = h hr = ⎜ ⎟ + ⎜ ⎟ + ⎜ ⎟ = ( sen senθθ cos ϕϕ ) + ( sen senθθ sen senϕϕ<br />
) + (cos ϕϕ<br />
) = 1<br />
⎝∂r ⎠ ⎝∂r ⎠ ⎝∂r⎠ 2 2 2<br />
⎛ ∂x ⎞ ⎛ ∂y ⎞ ⎛ ∂z<br />
⎞<br />
2 2 2<br />
h2 = hθ = ⎜ ⎟ + ⎜ ⎟ + ⎜ ⎟ = ( r cosθcosϕ) + ( r cos θsenϕ) + ( − rsenθ) = r<br />
⎝∂θ ⎠ ⎝∂θ ⎠ ⎝∂θ ⎠<br />
2 2 2<br />
⎛ ∂x ⎞ ⎛ ∂y ⎞ ⎛ ∂z<br />
⎞<br />
2 2<br />
h3 = hϕ = ⎜ ⎟ + ⎜ ⎟ + ⎜ ⎟ = ( − rsenθ senϕ) + ( rsenθcosϕ) = rsenθ<br />
⎝ ∂ ∂ϕϕ ⎠ ⎝ ∂ ∂ϕϕ ⎠ ⎝ ∂ ∂ϕϕ<br />
⎠
Un operador diferencial vectorial es un operador lineal<br />
que actúa sobre campos vectoriales definidos sobre una<br />
variedad diferenciables<br />
El operador vectorial nabla es un operador matemático que<br />
tiene carácter vectorial; sin embargo carece de algunas de las<br />
propiedades de las que gozan las magnitudes vectoriales,<br />
como por ejemplo el módulo módulo.<br />
En coordenadas coo de adas cartesianas ca tes a as se expresar epesapor po<br />
� ∂ � ∂ � ∂ �<br />
∇ = i + j + k<br />
∂ ∂x ∂ ∂y ∂ ∂z
Si el operador nabla se aplica escalarmente a un escalar se<br />
obtiene el gradiente de dicha función escalar; si se aplica<br />
escalarmente esca a e te a un u vector ecto se obtiene obt e e laa divergencia d e ge ca de laa<br />
función vectorial.<br />
Si el operador nabla se aplica vectorialmente a una función<br />
vectorial se obtiene el rotacional de la función vectorial
Gradiente de un campo escalar f(q 1,q 2q 3)<br />
Dada un campo escalar, escalar el gradiente de dicha función es un<br />
campo vectorial, que en coordenadas curvilíneas se expresa<br />
mediante ed a te la a ecuac ecuación ó<br />
� 1 ∂f � 1 ∂f � 1 ∂f � ⎛ 1 ∂ � 1 ∂ � 1 ∂ � ⎞<br />
∇ f = u1+ u2 + u3 = ⎜ u1+ u2 + u3 ⎟ f<br />
h1 ∂q1 h2 ∂q2 h3 ∂q3 ⎝h1 ∂q1 h2 ∂q2 h3 ∂q3<br />
⎠<br />
Siendo<br />
� ⎛ 1 ∂ � 1 ∂ � 1 ∂ �<br />
⎞<br />
∇= ⎜ u1+ u2 + u3<br />
⎟<br />
⎝h1 ∂q1 h2 ∂q2 h3 ∂q3<br />
⎠
Gradiente de un campo escalar f(q 1,q 2q 3)<br />
EEn un sistema i t de d coordenadas d d cartesianas t i<br />
� ∂f � ∂f � ∂f<br />
�<br />
∇ f = i + j + k<br />
∂x ∂y ∂z<br />
En coordenadas cilíndricas<br />
� ∂f � 1 ∂f � ∂f<br />
�<br />
∇ f = u + uϕ+ u<br />
∂ r r ∂ ∂ϕϕ ∂ z<br />
Y en coordenadas esféricas<br />
r z<br />
� ∂f � 1 ∂f � 1 ∂f<br />
�<br />
∇ f = ur+ u + u<br />
∂ r r ∂ ∂θθ rsen rsenθθ<br />
∂ ∂ϕϕ<br />
θ ϕ
Divergencia de un campo vectorial<br />
�<br />
Dado un campo vectorial vq ( 1, q q2, q q3<br />
) , la divergencia de dicho<br />
campo es un campo escalar, que se puede expresar como la<br />
aplicación ap cac ó de del ope operador ado nabla ab a escalarmente esca a e te po por eel ca campo po<br />
vectorial<br />
� � ⎛ 1 ∂ � 1 ∂ � 1 ∂ � ⎞ � � �<br />
∇ v = ⎜ u + u + u ⎟·<br />
vu + v u + v u<br />
⎝ h q h q h q ⎠<br />
( )<br />
1 2 3 1 1 2 2 3 3<br />
1 ∂ 1 2 ∂ 2 3 ∂ 3<br />
� �<br />
⎛ 1 ∂ ∂v 1 ∂ 1 ∂ ⎞ 1 ⎛∂( hh) ∂( h h) ∂(<br />
hh ) ⎞<br />
1 1 ∂v2<br />
1 ∂v ⎞ 3 1 ⎛∂( vhh 1 2 3) ∂( hv 1 2h3) ∂(<br />
hh 1 2v3) ⎞<br />
∇ v = ⎜ + + ⎟= ⎜ + + ⎟<br />
⎝h1 ∂q1 h2 ∂q2 h3 ∂q3 ⎠ hh 1 2h3 ⎝ ∂q1 ∂q2 ∂q3<br />
⎠
En un sistema de coordenadas cartesianas<br />
� � ∂v ∂v<br />
∂v<br />
∇ v = + +<br />
∂ ∂x ∂ ∂y ∂ ∂z<br />
En coordenadas cilíndricas<br />
Y en coordenadas d d esféricas fé i<br />
x y z<br />
� � 1 ∂( rv ) 1 ∂ ( vϕ<br />
) ∂v<br />
∇ v = + +<br />
r ∂r r ∂ϕ∂z r z<br />
� 2<br />
� 1 ∂(<br />
rvr)<br />
1 ∂(<br />
senθ v ) 1 ∂v<br />
θ<br />
ϕ<br />
∇ v = + +<br />
2<br />
r ∂ ∂rrsenθθ ∂ ∂θθ rsenθθ<br />
∂ ∂ϕϕ
Rotacional de un campo vectorial<br />
�<br />
Dado un campo vectorial vq ( ) , el rotacional de es el<br />
1, q2, q3)<br />
producto vectorial, del operador nabla por el vector, lo que<br />
resulta<br />
� �<br />
∇∧ v =<br />
� � �<br />
hu hu hu<br />
1 1 2 2 3 3<br />
1 ∂ ∂ ∂<br />
hhh ∂ q ∂ q ∂<br />
q<br />
hv hv hv<br />
1 2 3 1 2 3<br />
1 1 2 2 3 3
En un sistema de coordenadas cartesianas<br />
En coordenadas cilíndricas í<br />
Y en coordenadas esféricas<br />
�<br />
ur �<br />
ruϕ �<br />
uz<br />
� � 1 ∂<br />
∇∧ v =<br />
r ∂r ∂<br />
∂ϕ ∂<br />
∂z<br />
� �<br />
∇∧ v =<br />
v rv v<br />
r ϕ z<br />
1<br />
i j k � � �<br />
j<br />
� � ∂<br />
∇∧ v =<br />
∂x ∂<br />
∂y ∂<br />
∂z<br />
r<br />
2<br />
rsen r<br />
vxvy vz<br />
� � �<br />
u ru rsenθu θ ϕϕ<br />
∂ ∂ ∂<br />
θ ∂ ∂θ ∂ϕ<br />
v rv rsenθθv r<br />
θ ϕ
Laplaciano de un campo escalar f(q 1,q 2q 3)<br />
La expresión del laplaciano de un campo escalar f es un<br />
campo escalar, que en coordenadas curvilíneas se expresa<br />
mediante la ecuación<br />
⎡ ∂ ⎛hh∂f ⎞ ∂ ⎛hh ∂f ⎞ ∂ ⎛hh ∂f<br />
⎞⎤<br />
2 1 ⎡ ∂ ⎛hh 2 2∂f ⎞ ∂ ⎛hh 1 1∂f ⎞ ∂ ⎛hh 1 2∂f<br />
⎞⎤<br />
Δ f =∇ f = 3 ⎢ ⎜ ⎟+ ⎜ ⎟+<br />
⎜ ⎟⎥<br />
hhh 1 2 ⎢⎣ ∂q1⎝ h1 ∂q1⎠ ∂q2⎝ h2 ∂q2⎠ ∂q3⎝ h3 ∂q3⎠⎥⎦
En un sistema de coordenadas cartesianas<br />
En coordenadas cilíndricas<br />
2 2 2<br />
∂ f ∂ f ∂ f<br />
Δ f = + +<br />
2 2 2<br />
∂ ∂x ∂ ∂y ∂ ∂z<br />
1 ∂ ⎛ ∂ ⎞ 1 ⎛∂ ⎞ ∂<br />
Δ f = ⎜ r ⎟ + ⎜ ⎟ +<br />
r ∂r⎝ ∂r ⎠ r ⎝∂ϕ⎠ ∂z<br />
2 2<br />
f f f<br />
2 2 2<br />
Y en coordenadas d d esféricas fé i<br />
2<br />
1 ∂ ⎛ 2 ∂f⎞ 1 ∂ ⎛ ∂f ⎞ 1 ⎛∂ f ⎞<br />
Δ f = ⎜r ⎟+ senθ<br />
2 ⎜ ⎟+<br />
2 2 ⎜ 2 ⎟<br />
r ∂ r ⎝ ∂ r ⎠ r sen senθθ ∂ ∂θθ ⎝ ∂ ∂θθ ⎠ r sen θθ ⎝ ∂ ∂ϕϕ<br />
⎠