Tensor de tensión de Maxwell

Descripción matemática en electromagnetismo

El tensor de tensión de Maxwell (llamado así por James Clerk Maxwell ) es un tensor simétrico de segundo orden en tres dimensiones que se utiliza en el electromagnetismo clásico para representar la interacción entre las fuerzas electromagnéticas y el momento mecánico . En situaciones sencillas, 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 fuerza de Lorentz . Cuando la situación se vuelve más complicada, este procedimiento ordinario puede volverse poco práctico, con ecuaciones que abarcan varias líneas. Por lo tanto, es conveniente recopilar muchos de estos términos en el tensor de tensión de Maxwell y utilizar la aritmética tensorial para encontrar la respuesta al problema en cuestión.

En la formulación relativista del electromagnetismo, los nueve componentes del tensor de tensión de Maxwell aparecen, negados, como componentes del tensor electromagnético de tensión-energía , que es el componente electromagnético del tensor total de tensión-energía . Este último describe la densidad y el flujo de energía y momento en el espacio-tiempo .

Motivación

Como se describe a continuación, la fuerza electromagnética se escribe en términos de y . Utilizando el cálculo vectorial y las ecuaciones de Maxwell , se busca la simetría en los términos que contienen y , y la introducción del tensor de tensión de Maxwell simplifica el resultado. E {\displaystyle \mathbf {E} } B {\displaystyle \mathbf {B} } E {\displaystyle \mathbf {E} } B {\displaystyle \mathbf {B} }

Ecuaciones de Maxwell en unidades del SI en el vacío
(como referencia)
NombreForma diferencial
Ley de Gauss (en el vacío) E = ρ ε 0 {\displaystyle {\boldsymbol {\nabla }}\cdot \mathbf {E} ={\frac {\rho }{\varepsilon _{0}}}}
Ley de Gauss para el magnetismo B = 0 {\displaystyle {\boldsymbol {\nabla }}\cdot \mathbf {B} =0}
Ecuación de Maxwell-Faraday
(ley de inducción de Faraday)
× E = B t {\displaystyle {\boldsymbol {\nabla }}\times \mathbf {E} =-{\frac {\partial \mathbf {B} }{\partial t}}}
Ley circuital de Ampère (en el vacío)
(con corrección de Maxwell)
× B = μ 0 J + μ 0 ε 0 E t   {\displaystyle {\boldsymbol {\nabla }}\times \mathbf {B} =\mu _{0}\mathbf {J} +\mu _{0}\varepsilon _{0}{\frac {\partial \mathbf {E} }{\partial t}}\ }
  1. Comenzando con la ley de fuerza de Lorentz

    F = q ( E + v × B ) = ( E + v × B ) ρ d τ {\displaystyle {\begin{aligned}\mathbf {F} &=q(\mathbf {E} +\mathbf {v} \times \mathbf {B} )\\[3pt]&=\int (\mathbf {E} +\mathbf {v} \times \mathbf {B} )\rho \mathrm {d} \tau \end{aligned}}} La fuerza por unidad de volumen es

    f = ρ E + J × B {\displaystyle \mathbf {f} =\rho \mathbf {E} +\mathbf {J} \times \mathbf {B} }
  2. A continuación, y se pueden reemplazar por los campos y , utilizando la ley de Gauss y la ley circuital de Ampère : ρ {\displaystyle \rho } J {\displaystyle \mathbf {J} } E {\displaystyle \mathbf {E} } B {\displaystyle \mathbf {B} } f = ε 0 ( E ) E + 1 μ 0 ( × 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} }
  3. La derivada temporal se puede reescribir como algo que se pueda interpretar físicamente, es decir, el vector de Poynting . Si utilizamos la regla del producto y la ley de inducción de Faraday obtenemos y ahora podemos reescribir como y , entonces, al reunir los términos con y obtenemos t ( E × B ) = E t × B + E × B t = E t × B E × ( × E ) {\displaystyle {\frac {\partial }{\partial t}}(\mathbf {E} \times \mathbf {B} )={\frac {\partial \mathbf {E} }{\partial t}}\times \mathbf {B} +\mathbf {E} \times {\frac {\partial \mathbf {B} }{\partial t}}={\frac {\partial \mathbf {E} }{\partial t}}\times \mathbf {B} -\mathbf {E} \times ({\boldsymbol {\nabla }}\times \mathbf {E} )} f {\displaystyle \mathbf {f} } f = ε 0 ( E ) E + 1 μ 0 ( × B ) × B ε 0 t ( E × B ) ε 0 E × ( × 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 \mathbf {E} ),} E {\displaystyle \mathbf {E} } B {\displaystyle \mathbf {B} } f = ε 0 [ ( E ) E E × ( × E ) ] + 1 μ 0 [ B × ( × B ) ] ε 0 t ( E × B ) . {\displaystyle \mathbf {f} =\varepsilon _{0}\left[({\boldsymbol {\nabla }}\cdot \mathbf {E} )\mathbf {E} -\mathbf {E} \times ({\boldsymbol {\nabla }}\times \mathbf {E} )\right]+{\frac {1}{\mu _{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).}
  4. Parece que falta un término en la simetría en y , lo que se puede lograr insertando debido a la ley de Gauss para el magnetismo : Al eliminar los rizos (que son bastante complicados de calcular), el uso de la identidad del cálculo vectorial conduce a: E {\displaystyle \mathbf {E} } B {\displaystyle \mathbf {B} } ( B ) B {\displaystyle \left({\boldsymbol {\nabla }}\cdot \mathbf {B} \right)\mathbf {B} } f = ε 0 [ ( E ) E E × ( × E ) ] + 1 μ 0 [ ( B ) B B × ( × B ) ] ε 0 t ( E × B ) . {\displaystyle \mathbf {f} =\varepsilon _{0}\left[({\boldsymbol {\nabla }}\cdot \mathbf {E} )\mathbf {E} -\mathbf {E} \times ({\boldsymbol {\nabla }}\times \mathbf {E} )\right]+{\frac {1}{\mu _{0}}}\left[({\boldsymbol {\nabla }}\cdot \mathbf {B} )\mathbf {B} -\mathbf {B} \times \left({\boldsymbol {\nabla }}\times \mathbf {B} \right)\right]-\varepsilon _{0}{\frac {\partial }{\partial t}}\left(\mathbf {E} \times \mathbf {B} \right).} 1 2 ( A A ) = A × ( × A ) + ( A ) A , {\displaystyle {\frac {1}{2}}{\boldsymbol {\nabla }}(\mathbf {A} \cdot \mathbf {A} )=\mathbf {A} \times ({\boldsymbol {\nabla }}\times \mathbf {A} )+(\mathbf {A} \cdot {\boldsymbol {\nabla }})\mathbf {A} ,} f = ε 0 [ ( E ) E + ( E ) E ] + 1 μ 0 [ ( B ) B + ( B ) B ] 1 2 ( ε 0 E 2 + 1 μ 0 B 2 ) ε 0 t ( E × B ) . {\displaystyle \mathbf {f} =\varepsilon _{0}\left[({\boldsymbol {\nabla }}\cdot \mathbf {E} )\mathbf {E} +(\mathbf {E} \cdot {\boldsymbol {\nabla }})\mathbf {E} \right]+{\frac {1}{\mu _{0}}}\left[({\boldsymbol {\nabla }}\cdot \mathbf {B} )\mathbf {B} +(\mathbf {B} \cdot {\boldsymbol {\nabla }})\mathbf {B} \right]-{\frac {1}{2}}{\boldsymbol {\nabla }}\left(\varepsilon _{0}E^{2}+{\frac {1}{\mu _{0}}}B^{2}\right)-\varepsilon _{0}{\frac {\partial }{\partial t}}\left(\mathbf {E} \times \mathbf {B} \right).}
  5. Esta expresión contiene todos los aspectos del electromagnetismo y el momento y es relativamente fácil de calcular. Se puede escribir de forma más compacta introduciendo el tensor de tensión de Maxwell . Todos los términos de, excepto el último, se pueden escribir como la divergencia del tensor del tensor de tensión de Maxwell, lo que da: Como en el teorema de Poynting , el segundo término del lado derecho de la ecuación anterior se puede interpretar como la derivada temporal de la densidad de momento del campo EM, mientras que el primer término es la derivada temporal de la densidad de momento de las partículas masivas. De esta manera, la ecuación anterior será la ley de conservación del momento en la electrodinámica clásica; donde se ha introducido el vector de Poynting. σ i j ε 0 ( E i E j 1 2 δ i j E 2 ) + 1 μ 0 ( B i B j 1 2 δ i j B 2 ) . {\displaystyle \sigma _{ij}\equiv \varepsilon _{0}\left(E_{i}E_{j}-{\frac {1}{2}}\delta _{ij}E^{2}\right)+{\frac {1}{\mu _{0}}}\left(B_{i}B_{j}-{\frac {1}{2}}\delta _{ij}B^{2}\right).} f {\displaystyle \mathbf {f} } σ = f + ε 0 μ 0 S t , {\displaystyle {\boldsymbol {\nabla }}\cdot {\boldsymbol {\sigma }}=\mathbf {f} +\varepsilon _{0}\mu _{0}{\frac {\partial \mathbf {S} }{\partial t}}\,,} S = 1 μ 0 E × B . {\displaystyle \mathbf {S} ={\frac {1}{\mu _{0}}}\mathbf {E} \times \mathbf {B} .}

En la relación anterior para la conservación del momento, es la densidad de flujo de momento y juega un papel similar al del teorema de Poynting . σ {\displaystyle {\boldsymbol {\nabla }}\cdot {\boldsymbol {\sigma }}} S {\displaystyle \mathbf {S} }

La derivación anterior supone un conocimiento completo de y (tanto de cargas y corrientes libres como acotadas). Para el caso de materiales no lineales (como el hierro magnético con una curva BH), se debe utilizar el tensor de tensión de Maxwell no lineal. [1] ρ {\displaystyle \rho } J {\displaystyle \mathbf {J} }

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 dedujo anteriormente, viene dado por:

σ i j = ϵ 0 E i E j + 1 μ 0 B i B j 1 2 ( ϵ 0 E 2 + 1 μ 0 B 2 ) δ i j {\displaystyle \sigma _{ij}=\epsilon _{0}E_{i}E_{j}+{\frac {1}{\mu _{0}}}B_{i}B_{j}-{\frac {1}{2}}\left(\epsilon _{0}E^{2}+{\frac {1}{\mu _{0}}}B^{2}\right)\delta _{ij}} ,

donde es la constante eléctrica y es la constante magnética , es el campo eléctrico , es el campo magnético y es la delta de Kronecker . Con magnitudes gaussianas , viene dada por: ϵ 0 {\displaystyle \epsilon _{0}} μ 0 {\displaystyle \mu _{0}} E {\displaystyle \mathbf {E} } B {\displaystyle \mathbf {B} } δ i j {\displaystyle \delta _{ij}}

σ i j = 1 4 π ( E i E j + H i H j 1 2 ( E 2 + H 2 ) δ i j ) {\displaystyle \sigma _{ij}={\frac {1}{4\pi }}\left(E_{i}E_{j}+H_{i}H_{j}-{\frac {1}{2}}\left(E^{2}+H^{2}\right)\delta _{ij}\right)} ,

¿Dónde está el campo magnetizante ? H {\displaystyle \mathbf {H} }

Una forma alternativa de expresar este tensor es:

σ = 1 4 π [ E E + H H E 2 + H 2 2 I ] {\displaystyle {\overset {\leftrightarrow }{\boldsymbol {\sigma }}}={\frac {1}{4\pi }}\left[\mathbf {E} \otimes \mathbf {E} +\mathbf {H} \otimes \mathbf {H} -{\frac {E^{2}+H^{2}}{2}}\mathbb {I} \right]}

donde es el producto diádico , y el último tensor es la díada unitaria: {\displaystyle \otimes }

I ( 1 0 0 0 1 0 0 0 1 ) = ( x ^ x ^ + y ^ y ^ + z ^ z ^ ) {\displaystyle \mathbb {I} \equiv {\begin{pmatrix}1&0&0\\0&1&0\\0&0&1\end{pmatrix}}=\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 del tensor de tensión de Maxwell tiene unidades de momento por unidad de área por unidad de tiempo y da el flujo de momento paralelo al eje ésimo que cruza una superficie normal al eje ésimo (en la dirección negativa) por unidad de tiempo. i j {\displaystyle ij} i {\displaystyle i} j {\displaystyle j}

Estas unidades también pueden verse como unidades de fuerza por unidad de área (presión negativa), y el elemento del tensor también puede interpretarse como la fuerza paralela al eje th que sufre una superficie normal al eje th por unidad de área. De hecho, los elementos diagonales dan la tensión (tracción) que actúa sobre un elemento de área diferencial normal al eje correspondiente. A diferencia de las fuerzas debidas 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 al elemento. Esta cizalladura está dada por los elementos no diagonales del tensor de tensión. i j {\displaystyle ij} i {\displaystyle i} j {\displaystyle j}

Recientemente se ha demostrado que el tensor de tensión de Maxwell es la parte real de un tensor de tensión electromagnético complejo más general cuya parte imaginaria explica las fuerzas electrodinámicas reactivas. [2]

En magnetostática

Si el campo es solo magnético (lo que es en gran medida cierto en los motores, por ejemplo), algunos de los términos desaparecen y la ecuación en unidades del SI se convierte en:

σ i j = 1 μ 0 B i B j 1 2 μ 0 B 2 δ i j . {\displaystyle \sigma _{ij}={\frac {1}{\mu _{0}}}B_{i}B_{j}-{\frac {1}{2\mu _{0}}}B^{2}\delta _{ij}\,.}

Para objetos cilíndricos, como el rotor de un motor, esto se simplifica aún más:

σ r t = 1 μ 0 B r B t 1 2 μ 0 B 2 δ r t . {\displaystyle \sigma _{rt}={\frac {1}{\mu _{0}}}B_{r}B_{t}-{\frac {1}{2\mu _{0}}}B^{2}\delta _{rt}\,.}

donde es la fuerza cortante en la dirección radial (hacia afuera del cilindro) y es la fuerza cortante en la dirección tangencial (alrededor del cilindro). Es la fuerza tangencial que hace girar el motor. es la densidad de flujo en la dirección radial y es la densidad de flujo en la dirección tangencial. r {\displaystyle r} t {\displaystyle t} B r {\displaystyle B_{r}} B t {\displaystyle B_{t}}

En electrostática

En electrostática no están presentes los efectos del magnetismo. En este caso el campo magnético se anula, es decir , y obtenemos el tensor de tensión de Maxwell electrostático . Se da en forma de componentes por B = 0 {\displaystyle \mathbf {B} =\mathbf {0} }

σ i j = ε 0 E i E j 1 2 ε 0 E 2 δ i j {\displaystyle \sigma _{ij}=\varepsilon _{0}E_{i}E_{j}-{\frac {1}{2}}\varepsilon _{0}E^{2}\delta _{ij}}

y en forma simbólica por

σ = ε 0 E E 1 2 ε 0 ( E E ) I {\displaystyle {\boldsymbol {\sigma }}=\varepsilon _{0}\mathbf {E} \otimes \mathbf {E} -{\frac {1}{2}}\varepsilon _{0}(\mathbf {E} \cdot \mathbf {E} )\mathbf {I} }

donde está el tensor identidad apropiado generalmente . I {\displaystyle \mathbf {I} } ( {\displaystyle {\big (}} 3 × 3 ) {\displaystyle 3\times 3{\big )}}

Valor propio

Los valores propios del tensor de tensión de Maxwell vienen dados por:

{ λ } = { ( ϵ 0 2 E 2 + 1 2 μ 0 B 2 ) ,   ± ( ϵ 0 2 E 2 1 2 μ 0 B 2 ) 2 + ϵ 0 μ 0 ( E B ) 2 } {\displaystyle \{\lambda \}=\left\{-\left({\frac {\epsilon _{0}}{2}}E^{2}+{\frac {1}{2\mu _{0}}}B^{2}\right),~\pm {\sqrt {\left({\frac {\epsilon _{0}}{2}}E^{2}-{\frac {1}{2\mu _{0}}}B^{2}\right)^{2}+{\frac {\epsilon _{0}}{\mu _{0}}}\left({\boldsymbol {E}}\cdot {\boldsymbol {B}}\right)^{2}}}\right\}}

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, , se puede escribir como σ λ I {\displaystyle {\overleftrightarrow {\boldsymbol {\sigma }}}-\lambda \mathbf {\mathbb {I} } }

σ λ I = ( λ + V ) I + ϵ 0 E E T + 1 μ 0 B B T {\displaystyle {\overleftrightarrow {\boldsymbol {\sigma }}}-\lambda \mathbf {\mathbb {I} } =-\left(\lambda +V\right)\mathbf {\mathbb {I} } +\epsilon _{0}\mathbf {E} \mathbf {E} ^{\textsf {T}}+{\frac {1}{\mu _{0}}}\mathbf {B} \mathbf {B} ^{\textsf {T}}}

dónde

V = 1 2 ( ϵ 0 E 2 + 1 μ 0 B 2 ) {\displaystyle V={\frac {1}{2}}\left(\epsilon _{0}E^{2}+{\frac {1}{\mu _{0}}}B^{2}\right)}

Nosotros nos pusimos

U = ( λ + V ) I + ϵ 0 E E T {\displaystyle \mathbf {U} =-\left(\lambda +V\right)\mathbf {\mathbb {I} } +\epsilon _{0}\mathbf {E} \mathbf {E} ^{\textsf {T}}}

Aplicando el lema del determinante matricial una vez, esto nos da

det ( σ λ I ) = ( 1 + 1 μ 0 B T U 1 B ) det ( U ) {\displaystyle \det {\left({\overleftrightarrow {\boldsymbol {\sigma }}}-\lambda \mathbf {\mathbb {I} } \right)}=\left(1+{\frac {1}{\mu _{0}}}\mathbf {B} ^{\textsf {T}}\mathbf {U} ^{-1}\mathbf {B} \right)\det {\left(\mathbf {U} \right)}}

Al aplicarlo nuevamente se obtiene:

det ( σ λ I ) = ( 1 + 1 μ 0 B T U 1 B ) ( 1 ϵ 0 E T E λ + V ) ( λ V ) 3 {\displaystyle \det {\left({\overleftrightarrow {\boldsymbol {\sigma }}}-\lambda \mathbf {\mathbb {I} } \right)}=\left(1+{\frac {1}{\mu _{0}}}\mathbf {B} ^{\textsf {T}}\mathbf {U} ^{-1}\mathbf {B} \right)\left(1-{\frac {\epsilon _{0}\mathbf {E} ^{\textsf {T}}\mathbf {E} }{\lambda +V}}\right)\left(-\lambda -V\right)^{3}}

Del último multiplicando en el lado derecho, vemos inmediatamente que es uno de los valores propios. λ = V {\displaystyle \lambda =-V}

Para encontrar la inversa de , utilizamos la fórmula de Sherman-Morrison: U {\displaystyle \mathbf {U} }

U 1 = ( λ + V ) 1 ϵ 0 E E T ( λ + V ) 2 ( λ + V ) ϵ 0 E T E {\displaystyle \mathbf {U} ^{-1}=-\left(\lambda +V\right)^{-1}-{\frac {\epsilon _{0}\mathbf {E} \mathbf {E} ^{\textsf {T}}}{\left(\lambda +V\right)^{2}-\left(\lambda +V\right)\epsilon _{0}\mathbf {E} ^{\textsf {T}}\mathbf {E} }}}

Factorizando un término en el determinante, nos queda encontrar los ceros de la función racional: ( λ V ) {\displaystyle \left(-\lambda -V\right)}

( ( λ + V ) ϵ 0 ( E B ) 2 μ 0 ( ( λ + V ) + ϵ 0 E T E ) ) ( ( λ + V ) + ϵ 0 E T E ) {\displaystyle \left(-\left(\lambda +V\right)-{\frac {\epsilon _{0}\left(\mathbf {E} \cdot \mathbf {B} \right)^{2}}{\mu _{0}\left(-\left(\lambda +V\right)+\epsilon _{0}\mathbf {E} ^{\textsf {T}}\mathbf {E} \right)}}\right)\left(-\left(\lambda +V\right)+\epsilon _{0}\mathbf {E} ^{\textsf {T}}\mathbf {E} \right)}

Por lo tanto, una vez que resolvamos

( λ + V ) ( ( λ + V ) + ϵ 0 E 2 ) ϵ 0 μ 0 ( E B ) 2 = 0 {\displaystyle -\left(\lambda +V\right)\left(-\left(\lambda +V\right)+\epsilon _{0}E^{2}\right)-{\frac {\epsilon _{0}}{\mu _{0}}}\left(\mathbf {E} \cdot \mathbf {B} \right)^{2}=0}

obtenemos los otros dos valores propios.

Véase también

Referencias

  1. ^ Brauer, John R. (13 de enero de 2014). Actuadores y sensores magnéticos. John Wiley & Sons. ISBN 9781118754979.
  2. ^ Nieto-Vesperinas, Manuel; Xu, Xiaohao (12 de octubre de 2022). "El teorema complejo del tensor de tensión de Maxwell: el tensor de tensión imaginario y la fuerza reactiva del momento orbital. Un nuevo escenario subyacente a las fuerzas ópticas electromagnéticas". Light: Science & Applications . 11 (1): 297. doi :10.1038/s41377-022-00979-2. PMC 9556612 . PMID  36224170. 
  • David J. Griffiths , "Introducción a la electrodinámica", págs. 351-352, Benjamin Cummings Inc., 2008
  • John David Jackson, "Electrodinámica clásica, 3.ª edición", John Wiley & Sons, Inc., 1999
  • Richard Becker, "Campos electromagnéticos e interacciones", Dover Publications Inc., 1964
Retrieved from "https://en.wikipedia.org/w/index.php?title=Maxwell_stress_tensor&oldid=1244494697"