Articulo de referencia

tensor de tensiones de Maxwell

El tensor de tensiones de Maxwell (llamado así en honor a James Clerk Maxwell ) es un tensor simétrico de segundo orden en tres dimensiones que se utiliza en el electromagnetism...

El tensor de tensiones de Maxwell (llamado así en honor a 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 que actúan sobre la carga a partir de la ley de Lorentz . Cuando la situación se complica, este procedimiento habitual puede volverse impracticable, con ecuaciones que abarcan múltiples líneas. Por lo tanto, resulta conveniente agrupar muchos de estos términos en el tensor de tensiones de Maxwell y utilizar la aritmética tensorial para encontrar la solución al problema en cuestión.

En la formulación relativista del electromagnetismo, los nueve componentes del tensor de tensiones de Maxwell aparecen, negados, como componentes del tensor de energía-tensión electromagnético , que es el componente electromagnético del tensor de energía-tensión total . 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 demi{\displaystyle \mathbf {E} }yB{\displaystyle \mathbf {B} }. Utilizando el cálculo vectorial y las ecuaciones de Maxwell , se busca la simetría en los términos que contienenmi{\displaystyle \mathbf {E} }yB{\displaystyle \mathbf {B} }y la introducción del tensor de tensiones de Maxwell simplifica el resultado.

  1. Partiendo de la ley de fuerza de LorentzF=q(mi+v×B)=(mi+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=ρmi+J×B{\displaystyle \mathbf {f} =\rho \mathbf {E} +\mathbf {J} \times \mathbf {B} }
  2. Próximo,ρ{\displaystyle \rho }yJ{\displaystyle \mathbf {J} }pueden ser reemplazados por los camposmi{\displaystyle \mathbf {E} }yB{\displaystyle \mathbf {B} }, utilizando la ley de Gauss y la ley circuital de Ampère : F=ε0(mi)mi+1μ0(×B)×Bε0mit×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 puede interpretar físicamente, a saber, el vector de Poynting . Usando la regla del producto y la ley de inducción de Faraday se obtiene t(mi×B)=mit×B+mi×Bt=mit×Bmi×(×mi){\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} )} y ahora podemos reescribirF{\displaystyle \mathbf {f} }como F=ε0(mi)mi+1μ0(×B)×Bε0t(mi×B)ε0mi×(×mi),{\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} ),} luego recopilando términos conmi{\displaystyle \mathbf {E} }yB{\displaystyle \mathbf {B} }da F=ε0[(mi)mimi×(×mi)]+1μ0[B×(×B)]ε0t(mi×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.mi{\displaystyle \mathbf {E} }yB{\displaystyle \mathbf {B} }, lo cual se puede lograr insertando(B)B{\displaystyle \left({\boldsymbol {\nabla }}\cdot \mathbf {B} \right)\mathbf {B} }debido a la ley de Gauss para el magnetismo : F=ε0[(mi)mimi×(×mi)]+1μ0[(B)BB×(×B)]ε0t(mi×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).} Eliminando los rizos (que son bastante complicados de calcular), usando la identidad del cálculo vectorial12(AA)=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} ,} conduce a: F=ε0[(mi)mi+(mi)mi]+1μ0[(B)B+(B)B]12(ε0mi2+1μ0B2)ε0t(mi×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 tensiones de Maxwell . σijε0(miimij12δijmi2)+1μ0(BiBj12δijB2).{\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).} Todos excepto el último período deF{\displaystyle \mathbf {f} }se puede escribir como la divergencia tensorial del tensor de tensiones de Maxwell, lo que da como resultado: σ=F+ε0μ0St,{\displaystyle {\boldsymbol {\nabla }}\cdot {\boldsymbol {\sigma }}=\mathbf {f} +\varepsilon _{0}\mu _{0}{\frac {\partial \mathbf {S} }{\partial t}}\,,} Como en el teorema de Poynting , el segundo término del lado derecho de la ecuación anterior puede interpretarse como la derivada temporal de la densidad de momento del campo electromagnético, mientras que el primer término es la derivada temporal de la densidad de momento para 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 .S=1μ0mi×B.{\displaystyle \mathbf {S} ={\frac {1}{\mu _{0}}}\mathbf {E} \times \mathbf {B} .}

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

La derivación anterior supone un conocimiento completo de ambos.ρ{\displaystyle \rho }yJ{\displaystyle \mathbf {J} }(tanto cargas y corrientes libres como confinadas). Para el caso de materiales no lineales (como el hierro magnético con una curva B–H), debe utilizarse el tensor de tensiones de Maxwell no lineal. [ 1 ]

Ecuación

En física , el tensor de tensiones de Maxwell es el tensor de tensiones de un campo electromagnético . Como se derivó anteriormente, viene dado por:

σij=ϵ0miimij+1μ0BiBj12(ϵ0mi2+1μ0B2)δij{\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}},

dóndeϵ0{\displaystyle \epsilon _{0}}es la constante eléctrica ,μ0{\displaystyle \mu _{0}}la constante magnética ,mi{\displaystyle \mathbf {E} }el campo eléctrico ,B{\displaystyle \mathbf {B} }el campo magnético yδij{\displaystyle \delta _{ij}}es la delta de Kronecker . En el sistema gaussiano , viene dada por:

σij=14π(miimij+HiHj12(mi2+H2)δij){\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óndeH{\displaystyle \mathbf {H} }es el campo magnetizante .

Una forma alternativa de expresar este tensor es:

σ=14π[mimi+HHmi2+H22I]{\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]}

dónde{\displaystyle \otimes }es el producto diádico , y el último tensor es el diádico unitario:

I(100010001)=incógnita^incógnita^+y^y^+z^z^{\displaystyle \mathbb {I} \equiv {\begin{pmatrix}1&0&0\\0&1&0\\0&0&1\end{pmatrix}}=\mathbf {\hat {x}} \otimes \mathbf {\hat {x}} +\mathbf {\hat {y}} \otimes \mathbf {\hat {y}} +\mathbf {\hat {z}} \otimes \mathbf {\hat {z}} }

El elementoij{\displaystyle ij}del tensor de esfuerzos de Maxwell tiene unidades de momento por unidad de área por unidad de tiempo y da el flujo de momento paralelo a lai{\displaystyle i}el eje que cruza una superficie normal a laj{\displaystyle j}eje (en dirección negativa) por unidad de tiempo.

Estas unidades también pueden verse como unidades de fuerza por unidad de área (presión negativa), y laij{\displaystyle ij}El elemento del tensor también puede interpretarse como la fuerza paralela a lai{\displaystyle i}El eje th sufrió una superficie normal alj{\displaystyle j}eje por unidad de área. De hecho, los elementos diagonales proporcionan la tensión (fuerza de 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 experimenta una fuerza en una dirección que no es normal al elemento. Esta fuerza cortante viene dada por los elementos fuera de la diagonal del tensor de tensiones.

En magnetostática

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

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

En electrostática

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

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

y en forma simbólica por

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

dóndeI{\displaystyle \mathbf {I} }es el tensor identidad apropiado({\displaystyle {\big (}}generalmente3×3){\displaystyle 3\times 3{\big )}}.

Autovalor

Los valores propios del tensor de tensiones de Maxwell vienen dados por:

{λ}={(ϵ02mi2+12μ0B2), ±(ϵ02mi212μ0B2)2+ϵ0μ0(miB)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 autovalores se obtienen aplicando iterativamente el lema del determinante de matrices , junto con la fórmula de Sherman-Morrison .

Observando que la matriz de la ecuación característica,σλI{\displaystyle {\overleftrightarrow {\boldsymbol {\sigma }}}-\lambda \mathbf {\mathbb {I} } }, se puede escribir como

σλI=(λ+V)I+ϵ0mimiT+1μ0BBT{\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=12(ϵ0mi2+1μ0B2){\displaystyle V={\frac {1}{2}}\left(\epsilon _{0}E^{2}+{\frac {1}{\mu _{0}}}B^{2}\right)}

nosotros establecimos

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

Aplicando el lema del determinante de matrices una vez, obtenemos

det(σλI)=(1+1μ0BTU1B)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 de nuevo se obtiene,

det(σλI)=(1+1μ0BTU1B)(1ϵ0miTmiλ+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 del lado derecho, vemos inmediatamente queλ=V{\displaystyle \lambda =-V}es uno de los valores propios.

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

U1=(λ+V)1ϵ0mimiT(λ+V)2(λ+V)ϵ0miTmi{\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} }}}

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

((λ+V)ϵ0(miB)2μ0((λ+V)+ϵ0miTmi))((λ+V)+ϵ0miTmi){\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)+ϵ0mi2)ϵ0μ0(miB)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 autovalores.

Véase también

Referencias

  1. Brauer, John R. (13 de enero de 2014). Actuadores y sensores magnéticos . John Wiley & Sons. ISBN 9781118754979.
  • 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.ª ed.", John Wiley & Sons, Inc., 1999
  • Richard Becker, "Campos electromagnéticos e interacciones", Dover Publications Inc., 1964