Tensor de tensión de Maxwell
El tensor de tensión de Maxwell (llamado así en honor a James Clerk Maxwell) es un tensor simétrico de segundo orden utilizado en el electromagnetismo clásico para representar la interacción entre las fuerzas electromagnéticas y el momento mecánico. En situaciones simples, como una carga puntual que se mueve libremente en un campo magnético homogéneo, es fácil calcular las fuerzas sobre la carga a partir de la ley de fuerzas de Lorentz. Cuando la situación se vuelve más complicada, este procedimiento ordinario puede volverse imprácticamente difícil, con ecuaciones que abarcan varias líneas. Por tanto, es conveniente reunir muchos de estos términos en el tensor de tensión de Maxwell y utilizar la aritmética de tensor para encontrar la respuesta al problema que nos ocupa.
En la formulación relativista del electromagnetismo, el tensor de Maxwell aparece como parte del tensor electromagnético-energético que es el componente electromagnético del tensor total del estrés-energía. Este último describe la densidad y el flujo de energía y el impulso en tiempo espacial.
Motivación
Como se indica a continuación, la fuerza electromagnética está escrita en términos de E{displaystyle mathbf} y B{displaystyle mathbf {B}. Utilizando cálculos vectoriales y ecuaciones de Maxwell, la simetría se busca en los términos que contienen E{displaystyle mathbf} y B{displaystyle mathbf {B}, e introduciendo el tensor de estrés Maxwell simplifica el resultado.
Nombre | Forma diferencial |
---|---|
Ley de Gauss (en vacío) | Silencio Silencio ⋅ ⋅ E=*** *** ε ε 0{displaystyle {boldsymbol {nabla}cdot mathbf {E} ={frac {rho }{varepsilon ♪♪ |
Ley de Gauss para el magnetismo | Silencio Silencio ⋅ ⋅ B=0{displaystyle {boldsymbol {nabla}cdot mathbf {B} =0} |
Ecuación de Maxwell-Faraday (La ley de inducción de Faraday) | Silencio Silencio × × E=− − ∂ ∂ B∂ ∂ t{displaystyle {boldsymbol {nabla}times mathbf {E} =-{frac {partial mathbf {B}{partial t}} |
Ampère's circuital law (en vacío) (con la corrección de Maxwell) | Silencio Silencio × × B=μ μ 0J+μ μ 0ε ε 0∂ ∂ E∂ ∂ t {displaystyle {boldsymbol {nabla}times mathbf {B} =mu ~ {0}Mathbf {J} +mu} _{0}varepsilon ¿Por qué? |
- Empezando con la ley de la fuerza
Lorentz
el volumen de fuerza porF=q()E+v× × B)=∫ ∫ ()E+v× × B)*** *** dτ τ {displaystyle {begin{aligned}mathbf {F}=q(mathbf {E} +mathbf {v} times mathbf {B})[3pt] âTMa âTMa {E} +mathbf {v} times mathbf {B})mhomatedunidad
f=*** *** E+J× × B{displaystyle mathbf {f} =rho mathbf {E} +mathbf {J} times mathbf {B} - Siguiente, *** *** {displaystyle rho } y J{displaystyle mathbf {J} puede ser reemplazado por los campos E{displaystyle mathbf} y B{displaystyle mathbf {B}, usando la ley de Gauss y la ley de Ampère:
f=ε ε 0()Silencio Silencio ⋅ ⋅ E)E+1μ μ 0()Silencio Silencio × × B)× × B− − ε ε 0∂ ∂ E∂ ∂ t× × B{displaystyle mathbf {f} =varepsilon _{0}left({boldsymbol {nabla }cdot mathbf {E} right)mathbf {E} +{frac {1}{mu} _{0}}left({boldsymbol {nabla }times mathbf {B} right)times mathbf {B} -varepsilon _{0}{frac {partial mathbf {E}{partial t}times mathbf {B}
- El derivado del tiempo puede ser reescrito a algo que se puede interpretar físicamente, a saber, el vector Poynting. Usando la regla del producto y la ley de inducción de Faraday da
y ahora podemos reescribir f{displaystyle mathbf {f} como∂ ∂ ∂ ∂ t()E× × B)=∂ ∂ E∂ ∂ t× × B+E× × ∂ ∂ B∂ ∂ t=∂ ∂ E∂ ∂ t× × B− − E× × ()Silencio Silencio × × E){fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif}entonces coleccionar términos con E{displaystyle mathbf} y B{displaystyle mathbf {B} daf=ε ε 0()Silencio Silencio ⋅ ⋅ E)E+1μ μ 0()Silencio Silencio × × B)× × B− − ε ε 0∂ ∂ ∂ ∂ t()E× × B)− − ε ε 0E× × ()Silencio Silencio × × E),{displaystyle mathbf {f} =varepsilon _{0}left({boldsymbol {nabla }cdot mathbf {E} right)mathbf {E} +{frac {1}{mu} _{0}}left({boldsymbol {nabla }times mathbf {B} right)times mathbf {B} -varepsilon _{0}{frac {partial }{partial t}left(mathbf {E} times mathbf {B}right)-varepsilon _{0}mathbf {E} times ({boldsymbol {nabla }times {i} {f}f}f}f}f=ε ε 0[()Silencio Silencio ⋅ ⋅ E)E− − E× × ()Silencio Silencio × × E)]+1μ μ 0[− − B× × ()Silencio Silencio × × B)]− − ε ε 0∂ ∂ ∂ ∂ t()E× × B).{fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif}cdot mathbf {})mathbf {E} -mathbf {E} times ({boldsymbol {nabla }}timesmathbf} {i} {i}} {i}}}} {i}i}}}}}f} {f}}f}}}}}}f}}}}}}}}f}f}f}f} {f}f} {f}f} _{0}}left[-mathbf {B} times left({boldsymbol {nabla }times mathbf {B}right)right]-varepsilon _{0}{frac {partial }{partial t}left(mathbf {E} times mathbf {B} right). }
- Un término parece estar "perdiendo" de la simetría en E{displaystyle mathbf} y B{displaystyle mathbf {B}, que se puede lograr mediante la inserción ()Silencio Silencio ⋅ ⋅ B)B{displaystyle left({boldsymbol {nabla }cdot mathbf {B} right)mathbf {B} por la ley de Gauss para el magnetismo:
Eliminar los rizos (que son bastante complicados para calcular), utilizando la identidad del cálculo vectorialf=ε ε 0[()Silencio Silencio ⋅ ⋅ E)E− − E× × ()Silencio Silencio × × E)]+1μ μ 0[()Silencio Silencio ⋅ ⋅ B)B− − B× × ()Silencio Silencio × × B)]− − ε ε 0∂ ∂ ∂ ∂ t()E× × B).{fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif}cdot mathbf {})mathbf {E} -mathbf {E} times ({boldsymbol {nabla }}timesmathbf} {i} {i}} {i}}}} {i}i}}}}}f} {f}}f}}}}}}f}}}}}}}}f}f}f}f} {f}f} {f}f} ¿Por qué? _{0}{frac {partial }{partial t}left(mathbf {E} times mathbf {B} right). }conduce a:12Silencio Silencio ()A⋅ ⋅ A)=A× × ()Silencio Silencio × × A)+()A⋅ ⋅ Silencio Silencio )A,{fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} {cdot mathbf {} {fnMicrosoft Sans Serif} {fnMicrosoft} {fnMicrosoft} {f}fnMicrosoft {fnMicrosoft} {f}f}f}fnMicrosoft}s}s}f}fnMicrosoft {fnMicrosoft}fnMicrosoft}fnMicrosoft}s}fnMicrosoft}s}fnMicrosoft}fnMis}fnMinMicrosoft {fnMis}fnMis {fnMicrosoft {s {fnMis {f}f}fnMicrosoft {fnMicrosoft {f}f}fnMicrosoft}fnMicrosoft}fnMicrof=ε ε 0[()Silencio Silencio ⋅ ⋅ E)E+()E⋅ ⋅ Silencio Silencio )E]+1μ μ 0[()Silencio Silencio ⋅ ⋅ B)B+()B⋅ ⋅ Silencio Silencio )B]− − 12Silencio Silencio ()ε ε 0E2+1μ μ 0B2)− − ε ε 0∂ ∂ ∂ ∂ t()E× × B).{fnMitbf {f}cdot mathbf {})mathbf {f}cdot mathbf {})mathbf {E} +(mathbf {E} {i} {i} {i} {i} {i}i} {i}i}i}f}] {fnK}cdot mathbf {B}cdot {B})mathbf {B} +(mathbf {B}cdot {boldsymbol {nabla })mathbf {B} {right]-{bold}{2}{0}}}} {}}}}} {m}} {m} {b} {}}}}}}}}b}b}}} {b} {}}}}}}}}}} {b} {b}b}}}}}}}b}}b}}}}}}}}b}}}}}b}b}}}b}b}}}b}b}b}bh}bh}bh} {bh}}}bh}}bh}bh} {bh}}}}bh}}bh}bh} ¿Qué? ¿Por qué? _{0}{frac {partial }{partial t}left(mathbf {E} times mathbf {B} right). }
- Esta expresión contiene todos los aspectos del electromagnetismo y el impulso y es relativamente fácil de calcular. Se puede escribir más compactamente introduciendo el Tensor de estrés Maxwell,
Todo menos el último mandato f{displaystyle mathbf {f} se puede escribir como la divergencia tensora del tensor de estrés Maxwell, dando:σ σ ij↑ ↑ ε ε 0()EiEj− − 12δ δ ijE2)+1μ μ 0()BiBj− − 12δ δ ijB2).{displaystyle sigma _{ij}equiv varepsilon ¿Qué? {1}{2}delta _{2}right)+{frac {1}{mu} - ¿Qué? (B_{i}B_{j}-{frac {1}{2}delta _{ij}B^{2}right). }Como en el teorema de Poynting, el segundo término en el lado derecho de la ecuación anterior se puede interpretar como el derivado del tiempo de la densidad de impulso del campo EM, mientras que el primer término es el derivado del tiempo de la densidad de impulso para las partículas masivas. De esta manera, la ecuación anterior será la ley de conservación del impulso en la electrodinámica clásica; donde se ha introducido el vector PoyntingSilencio Silencio ⋅ ⋅ σ σ =f+ε ε 0μ μ 0∂ ∂ S∂ ∂ t,{displaystyle {boldsymbol {nabla}cdot {boldsymbol {sigma} }=Mathbf {f} +varepsilon _{0}mu _{0}{frac {partial mathbf {S} {fncipal t},}S=1μ μ 0E× × B.{displaystyle mathbf {S} ={frac {1}{mu} - Mathbf {E} times mathbf {B}
en la relación anterior para la conservación del impulso, Silencio Silencio ⋅ ⋅ σ σ {displaystyle {boldsymbol {nabla}cdot {boldsymbol {sigma} } es densidad de flujo y juega un papel similar a S{displaystyle mathbf {S} en el teorema de Poynting.
La derivación anterior asume conocimiento completo de ambos *** *** {displaystyle rho } y J{displaystyle mathbf {J} (ambos cargos y corrientes libres y obligados). Para el caso de materiales no lineales (como el hierro magnético con una curva BH), el tensor de tensión Maxwell no lineal debe ser utilizado.
Ecuación
En física, el tensor de tensión de Maxwell es el tensor de tensión de un campo electromagnético. Como se derivó anteriormente en unidades SI, viene dado por:
- σ σ ij=ε ε 0EiEj+1μ μ 0BiBj− − 12()ε ε 0E2+1μ μ 0B2)δ δ ij{displaystyle sigma _{ij}=epsilon ¿Qué? ¿Qué? {1}{2}left(epsilon) ¿Qué? ¿Por qué? ¿Qué?,
Donde ε ε 0{displaystyle epsilon _{0} es la constante eléctrica y μ μ 0{displaystyle mu _{0}} es la constante magnética, E{displaystyle mathbf} es el campo eléctrico, B{displaystyle mathbf {B} es el campo magnético y δ δ ij{displaystyle delta _{ij} Es el delta de Kronecker. En Gaussian cgs unit, se da por:
- σ σ ij=14π π ()EiEj+HiHj− − 12()E2+H2)δ δ ij){displaystyle sigma {fnMicroc} {fnMicroc} {fnK} {fnMicrosoft} {fnMicrosoft} {1}{2}left(E^{2}+H^{2}right)delta _{ij}right)},
Donde H{displaystyle mathbf {H} es el campo de magnetización.
Una forma alternativa de expresar este tensor es:
- σ σ Administración Administración =14π π [E⊗ ⊗ E+H⊗ ⊗ H− − E2+H22I]{\displaystyle {\\ssets {s\\s\\\\\\\\\\\\\\\\\\\\\\\\\cH00\\\cH00\\\cH\\\\\\cH\\\\\\\\\cH\\\\\\\\cH\\\\\\cH\\\\\\\\\\cH\\\\\\cH\\cH\\cH\\\\\\\\cH\\\\\\cH\\\\\ ##{boldsymbol {sigma }={frac {1}{4pi}left[mathbf] {E} otimes mathbf {E} +mathbf {H} otimes mathbf {fnMicroc {fnMicroc}}mhbb {I} {I} {fnh}}
Donde ⊗ ⊗ {displaystyle otimes } es el producto dyadic, y el último tensor es el dyad unidad:
- I↑ ↑ ()100010001)=()x^ ^ ⊗ ⊗ x^ ^ +Sí.^ ^ ⊗ ⊗ Sí.^ ^ +z^ ^ ⊗ ⊗ z^ ^ ){displaystyle mathbb {I} equiv {begin{pmatrix}1 tendría un doble0}otimes mathbf {hat {x}}=left(mathbf {hat {x} otimes mathbf {hat {x}} # Mathbf {hat {y} otimes mathbf {hat {y} # Mathbf {hat {z} otimes mathbf {hat {z} right)}
El elemento ij{displaystyle ij} del tensor de estrés Maxwell tiene unidades de impulso por unidad de tiempo de área por unidad y da el flujo de impulso paralelo al i{displaystyle i}eje cruzando una superficie normal al j{displaystyle j}eje (en la dirección negativa) por unidad de tiempo.
Estas unidades también se pueden ver como unidades de fuerza por unidad de área (presión negativa) y la ij{displaystyle ij} elemento del tensor también se puede interpretar como la fuerza paralela al i{displaystyle i}el eje sufrido por una superficie normal al j{displaystyle j}eje por unidad de área. De hecho, los elementos diagonales dan la tensión (pulido) actuando en un elemento de área diferencial normal al eje correspondiente. A diferencia de las fuerzas debido a la presión de un gas ideal, un elemento de área en el campo electromagnético también siente una fuerza en una dirección que no es normal para el elemento. Este tirón es dado por los elementos fuera de la diagonal del tensor de estrés.
En magnetostática
Si el campo es sólo magnético (lo que es cierto en gran medida en los motores, por ejemplo), algunos de los términos desaparecen y la ecuación en unidades SI se convierte en:
- σ σ ij=1μ μ 0BiBj− − 12μ μ 0B2δ δ ij.{displaystyle sigma {fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} ¿Qué? {1}{2mu ¿Por qué? - ¿Qué?
Para objetos cilíndricos, como el rotor de un motor, esto se simplifica aún más a:
- σ σ rt=1μ μ 0BrBt− − 12μ μ 0B2δ δ rt.{displaystyle sigma ¿Qué? ¿Qué? {1}{2mu ¿Por qué? - ¿Qué?
Donde r{displaystyle r} es el cobertizo en la dirección radial (hacia fuera del cilindro) y t{displaystyle t} es el cobertizo en la dirección tangencial (alrededor del cilindro). Es la fuerza tangencial que gira el motor. Br{displaystyle B_{r} es la densidad del flujo en la dirección radial, y Bt{displaystyle B_{t} es la densidad del flujo en la dirección tangencial.
En electrostática
En electrostáticos los efectos del magnetismo no están presentes. En este caso el campo magnético desaparece, es decir, B=0{displaystyle mathbf {B} =mathbf {0}, y obtenemos electrostática Tensor de estrés Maxwell. Se da en forma de componente por
- σ σ ij=ε ε 0EiEj− − 12ε ε 0E2δ δ ij{displaystyle sigma _{ij}=varepsilon ¿Qué? {1}{2}varepsilon ¿Qué? ¿Qué?
y en forma simbólica por
- σ σ =ε ε 0E⊗ ⊗ E− − 12ε ε 0()E⋅ ⋅ E)I{displaystyle {boldsymbol {sigma }= 'varepsilon ♪♪Mathbf {E} otimes mathbf Varepsilon ¿Por qué?
Donde I{displaystyle mathbf} es el tensor de identidad adecuado (){displaystyle {big}generalmente 3× × 3){displaystyle 3times 3{big]}.
Valor propio
Los valores propios del tensor de tensión de Maxwell vienen dados por:
- {}λ λ }={}− − ()ε ε 02E2+12μ μ 0B2), ± ± ()ε ε 02E2− − 12μ μ 0B2)2+ε ε 0μ μ 0()E⋅ ⋅ B)2}{displaystyle {lambda #=left{frac {epsilon ¿Qué? - Sí. {fnMicroc {epsilon} {0}{2}}E^{2}-{frac {1}{2mu} ¿Por qué? ¿Qué? ¿Qué? {E}cdot {boldsymbol Bien.
Estos valores propios se obtienen aplicando iterativamente el lema del determinante matricial, junto con la fórmula de Sherman-Morrison.
Observando que la matriz de ecuación característica, σ σ Administración Administración − − λ λ I{\displaystyle {\\deleftrightarrow {boldsymbol {sigma # Lambda mathbf {mathbb {I}, puede ser escrito como
- σ σ Administración Administración − − λ λ I=− − ()λ λ +V)I+ε ε 0EET+1μ μ 0BBT{\displaystyle {\\deleftrightarrow {boldsymbol {sigma # Lambda mathbf {mathbb {I} =-left(lambda) +Vright)mathbf {mathbb {I} +epsilon _{0}mathbf {E} mathbf {E} {textosf {T}+{frac {1}{mu} ¿Por qué?
dónde
- V=12()ε ε 0E2+1μ μ 0B2){displaystyle V={frac}{2}left(epsilon) ¿Qué? - Sí.
nosotros configuramos
- U=− − ()λ λ +V)I+ε ε 0EET{displaystyle mathbf {U} =-left(lambda) +Vright)mathbf {mathbb {I} } +epsilon _{0}mathbf {E} ^{textsf {T}}
Aplicando el lema del determinante de la matriz una vez, esto nos da
- Det()σ σ Administración Administración − − λ λ I)=()1+1μ μ 0BTU− − 1B)Det()U){displaystyle det {left({overleftrightarrow {boldsymbol {sigma # Lambda mathbf {mathbb {I} } right)}=left(1+{frac {1}{mu} Mathbf {U} ^{-1}mathbf {B}det {left(mathbf {U}right)}
Aplicarlo nuevamente produce,
- Det()σ σ Administración Administración − − λ λ I)=()1+1μ μ 0BTU− − 1B)()1− − ε ε 0ETEλ λ +V)()− − λ λ − − V)3{displaystyle det {left({overleftrightarrow {boldsymbol {sigma # Lambda mathbf {mathbb {I} } right)}=left(1+{frac {1}{mu} Mathbf {U} ^{-1}mathbf {B} right)left(1-{frac {epsilon _{0}mathbf {E} {textosf {T}Mathbf Está bien. -V 'right'
De la última multiplicación en el RHS, vemos inmediatamente que λ λ =− − V{displaystyle lambda = V} es uno de los eigenvalues.
Para encontrar el inverso de U{displaystyle mathbf {U}, usamos la fórmula Sherman-Morrison:
- U− − 1=− − ()λ λ +V)− − 1− − ε ε 0EET()λ λ +V)2− − ()λ λ +V)ε ε 0ETE{displaystyle mathbf {fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} {fnMicroc {epsilon _{0}mathbf {E} mathbf {E} {textosf {T}}{left(lambda +Vright)}left(lambda +Vright)epsilon ♪♪Mathbf {E} {textosf {T} 'Mathbf {E}
Factoring out a ()− − λ λ − − V){displaystyle left(-lambda -Vright)} término en el determinante, nos quedamos con encontrar los ceros de la función racional:
- ()− − ()λ λ +V)− − ε ε 0()E⋅ ⋅ B)2μ μ 0()− − ()λ λ +V)+ε ε 0ETE))()− − ()λ λ +V)+ε ε 0ETE){displaystyle left(-left(lambda +Vright)-{epsilon _{0}left(mathbf {E} cdot mathbf {B} right)}{2}{mu _{0}left(left)lambda +Vright)+epsilon _{0}mathbf {E} {textosf Está bien. +Vright)+epsilon _{0}mathbf {E} {textsf {T}mathbf {E}right)}
Así, una vez que resolvamos
- − − ()λ λ +V)()− − ()λ λ +V)+ε ε 0E2)− − ε ε 0μ μ 0()E⋅ ⋅ B)2=0{displaystyle -left(lambda +Vright)left(left(lambda) +Vright)+epsilon ¿Por qué? ¿Qué? ¿Por qué? {B} right)} {2}=0}
obtenemos los otros dos valores propios.
Contenido relacionado
Julio (unidad)
Ley de Fick
Pascal (unidad)
Newton (unidad)
Trabajo (física)