Clase 08 - Ejemplos de cálculo.

Prof. Juan Mauricio Matera

03 de abril de 2019

Cálculo de potenciales y campos

  • Descomponer el problema como superposición de distribuciones simétricas.
  • Utilizar la Ley de Gauss para calcular la contribución al campo de cada una de esas distribuciones.
  • Calcular el campo total como suma de las contribuciones.

Ejemplo: esfera uniformemente cargada con un agujero esférico.

  • Descomponemos la distribución como una esfera con densidad de carga uniforme, y otra esfera más pequeña con densidad de carga uniforme y opuesta, en la región del agujero.
  • Calculamos los campos de ambas distribuciones. \[ \vec{E}_1 = \frac{\rho \vec{r}}{3 \varepsilon_0}\times \left\{ ^{1}_{\frac{R^3}{|\vec{r}|^3}} \;\;^{|\vec{r}|<R}_{|\vec{r}|>R} \right. \] \[ \vec{E}_2 = \frac{-\rho (\vec{r}-\vec{a})}{3 \varepsilon_0}\times \left\{ ^{1}_{\frac{b^3}{|\vec{r}-\vec{a}|^3}} \;\;^{|\vec{r}-\vec{a}|<b}_{|\vec{r}-\vec{a}|>b} \right. \]

  • Dentro de la cavidad, \[ \vec{E}_1+\vec{E}_2= \frac{\rho \vec{r}}{3 \varepsilon_0}- \frac{\rho \vec{r}-\vec{a}}{3 \varepsilon_0} = \frac{\rho \vec{a}}{3 \varepsilon_0} \]

Ejemplo: campos y potenciales alrededor de un átomo

  • En el modelo atómico, asumimos que el núcleo corresponde a una carga puntual positiva, de magnitud \(Z e\), rodeada por una nube de carga total \(-Z e\), con una densidad volumétrica de carga proporcional a \(\rho(r)\approx e^{-r/a}\) donde \(a=0,5 \times 10^{-10}{\rm m}\).

  • Queremos calcular el campo y el potencial asociado a la distribución.

  • De la simetría del problema, y de la Ley de Gauss \[ \vec{E}(\vec{r})=\frac{Q(r) \vec{r}}{4\pi \varepsilon_0 r^3} \] con \(Q(r)\) la carga encerrada por una esfera gaussiana de radio \(r\). \[ Q(r)= Z e + \int_{r'<r} \rho(r') dV' \]
  • \(\rho(r)=C \exp(-r/a)\). Para determinar \(C\), integramos en todo el espacio e igualamos a la carga total \(-Z e\): \[\begin{eqnarray*} -Z e &=& \int \rho(r)dV = \int_{0}^{\infty} \int_0^\theta\int_{-\pi}^{\pi} C \exp(r/a) r^2\sin( \theta) d\phi d\theta dr\\ &=&4\pi C a^3 \int_0^{\infty}\exp(u)u^2du= 8\pi C a^3\\ C&=& \frac{-Z e}{8 \pi a^3} \end{eqnarray*}\]
  • Integrando de la misma manera, pero a un \(r\) finito, \[ Q(r)= Z e (1-\int_0^{r/a} \frac{\exp(-u)}{2}u^{2} du)= Z e \exp(-r/a) (1+\frac{r}{a}+\frac{r^2}{a^2}) \]
  • Finalmente,

\[ \vec{E}(\vec{r})= \frac{Z e}{4\pi\varepsilon_0}\exp(-r/a)(\frac{1}{r^2}+\frac{1}{a r}+\frac{1}{a^2}) \check{r} \]

  • Para construir el potencial, buscamos una función que
    • Se anule para \(r\rightarrow \infty\)
    • Su derivada respecto a \(r\) sea (menos) la componente radial del campo
  • Propongo entonces
    \[ V(\vec{r})= \frac{Z e}{4 \pi \varepsilon_0}e^{-r/a}\left(\frac{1}{r}+\frac{1}{a}\right) \]

Átomo descentrado

  • Ahora que conocemos el potencial para el caso centrado, podemos resolver el caso en que el centro de la nube electrónica no coincide con el núcleo, sino que está desplazado en un vector \(\vec{d}\)

  • El potencial será la superposición del potencial de la carga puntual en \(\vec{d}\), con el de la nube, centrada en el origen:

\[\begin{eqnarray*} V_{\rm nucleo}(\vec{r}) &=& \frac{Z e}{4\pi\varepsilon_0 |\vec{r}-\vec{d}|}\\ V_{\rm nube}(\vec{r}) &=& \frac{Z e}{4\pi\varepsilon_0}\left(\exp(-r/a)(\frac{1}{r}+\frac{1}{a})-\frac{1}{r}\right)\\ V(\vec{r})&=&V_{\rm nucleo}(\vec{r})+V_{\rm nube}(\vec{r})= \frac{Z e}{4\pi\varepsilon_0}\left(\exp(-r/a)(\frac{1}{|\vec{r}|}+\frac{1}{a}) + \frac{1}{|\vec{r}-\vec{d}|}-\frac{1}{|\vec{r}|}\right) \end{eqnarray*}\]

  • A distancias \(r\gg a\), sólo contribuyen los últimos dos términos, que pueden aproximarse por el potencial de un dipolo: \[ V(\vec{r})\approx \frac{Z e}{4\pi\varepsilon_0}\left(\frac{1}{|\vec{r}-\vec{d}|}-\frac{1}{|\vec{r}|}\right)\approx \frac{Z e\vec{d}\cdot\check{r}}{4\pi\varepsilon_0|\vec{r}|^2} \]

Cálculo de capacidad en sistemas de conductores.

  • Determinar la distribución de cargas en los conductores bajo la condición \(\vec{E}=0\) en el interior de cada uno de ellos y la ley de Gauss.
  • Si el sistema tiene simetría plana, cilíndrica o esférica
    • Proponer la forma del campo eléctrico
    • Elegir superficies de Gauss de manera conveniente.
    • Determinar la componente desconocida vía la ley de Gauss.
    • Elegir un referencial adecuado y calcular el potencial electrostático.
  • Si el sistema no presenta simetrías
    • Proponer una distribución de carga.
    • Descomponer la distribución de cargas en elementos de carga.
    • Calcular las contribuciones al potencial de cada elemento de carga.
    • Sumar/integrar las contribuciones.
    • Calcular el campo eléctrico como el gradiente del potencial.

Capacidad de un conductor coaxial

  • Cargas distribuídas sobre las superficies de los conductores.
  • \(Q_1={2\pi r_1 L} \sigma_1\)
  • \(Q_2=2\pi L( r_{2,i} \sigma_{2,i} + r_{2,e} \sigma_{2,e})\)
  • \(Q_3={2\pi r_3 L} \sigma_1\)
  • \(\sigma_{2,i}=-\sigma_{1}\).

  • Simetría cilíndrica
    • \(z\leftrightarrow -z\)
    • \(x \leftrightarrow -x\), \(y\leftrightarrow -y\)
    • \(z \leftrightarrow z+\Delta z\)
    • Rotaciones alrededor de \(z\): \(x\rightarrow x \cos(\phi)-y\sin(\phi)\), \(y\rightarrow y \cos(\phi) + x\sin(\phi)\),
  • Forma del campo eléctrico: \[ \vec{E}(r \hat{r}+z\hat{z})=\vec{E}(r)\hat{r} \;\;\;\;\;\hat{r}\perp \hat{z}\,\;\;\; r=\sqrt{x^2+y^2} \]

  • Superficie gaussiana: Cilíndro de radio \(r\) y altura \(h\) concéntrico.

  • Ley de Gauss: \(\int_{\cal S}\vec{E}\cdot d\vec{\cal S}=\frac{Q({\cal S})}{\varepsilon_0}=\frac{\lambda(r)h}{\varepsilon_0}\) \[\begin{eqnarray*} \int_{\cal S}\vec{E}\cdot d\vec{S}&=&\int_{{\cal S}_{superior}}\vec{E}\cdot d\vec{\cal S} + \int_{{\cal S}_{inferior}}\vec{E}\cdot d\vec{\cal S} +\int_{{\cal S}_{lateral}}\vec{E}\cdot d\vec{\cal S} \end{eqnarray*}\]

Pero \(\vec{E}\cdot d\vec{\cal S}=E(r)\hat{r}\cdot(\pm \vec{z}) dS =0\) sobre las tapas: \[ \int_{{\cal S}_{superior}}\vec{E}\cdot d\vec{\cal S} = \int_{{\cal S}_{inf}}\vec{E}\cdot d\vec{\cal S}=0 \] y \(\vec{E}\cdot d\vec{\cal S}= E(r)\hat{r} \cdot (\hat{r} dS)=E(r)dS\) sobre la superficie lateral. Luego, \[ \int_{\cal S}\vec{E}\cdot d\vec{S}=\vec{E}(r)\int_{{\cal S}_{lateral}} dS= 2\pi r h E(r) \] Finalmente, \[ E(r)=\frac{\lambda(r)}{2\pi r \varepsilon_0} \]

\[ \lambda(r)=\left\{ \begin{array}{l r} 0 & r<r_1\\ 2\pi r_1 \sigma_1 & r_{2,i}>r>r_1\\ 0 & r_{2,e}>r>r_{2,i}\\ 2\pi r_{2,e} (\sigma_{2,e}) & r_3>r>r_{2,e}\\ 2\pi r_e (\sigma_{2,e}+\sigma_3) & r>r_3 \end{array} \right. \;\;\;\;\; \]

Pero,

  • el conductor \(3\) está conectado a tierra. Luego, \(V(r)-V(\infty)=0\rightarrow E(r)=0\) para \(r>r_3\). Por lo tanto, \(\sigma_3=-\sigma_{2,e}\)
  • \(\sigma_{2,e}= \frac{Q_1 + Q_2}{2\pi L r_{2,e}}\).

Potencial

\[ V(\vec{r})=-\int_{\infty}^{\vec{r}} \vec{E}\cdot d\vec{\ell} \]

  • Elegimos una trayectoria “radial” \(\vec{\ell}(u)=\frac{\vec{r}}{u}\) con \(u\in(0,1)\). Luego, \(d\vec{\ell}= \frac{-\vec{r}}{u^2} du\) y \[ V(\vec{r})=\int_{0}^{1} \vec{E}\left(\frac{\vec{r}}{u}\right)\cdot \frac{\vec{r}}{u^2}du= \int_{0}^{1} \frac{\lambda(|r|/u)}{2\pi \varepsilon_0 (|r|/u)} \frac{|\vec{r}|}{u^2}du \] con \(\lambda((\vec{r}/u))\) constante a trozos.
  • Si \(|\vec{r}|>r_3\), \(\lambda(|\vec{r}|/u)=0\) sobre toda la trayectoria y por lo tanto, \(V(\vec{r})=0\).
  • Si \(r_{3}>|\vec{r}|>r_{e,2}\), la integral se reduce al intervalo \(u\in (\frac{|\vec{r}|}{r_{3}},1)\) y \[ V(r)=\int_{\frac{|\vec{r}|}{r_{3}}}^{1}\frac{\sigma_{2,e}r_{2,e}|\vec{r}|}{\varepsilon_0 |r|)} \frac{1}{u}du= \frac{\sigma_{2,e}r_{2,e}}{\varepsilon_0}\int_{\frac{|\vec{r}|}{r_{3}}}^{1} \frac{1}{u}du= \frac{\sigma_{2,e}r_{2,e}}{\varepsilon_0}\log(\frac{r_{3}}{|\vec{r}|}) \]
  • Si \(r_{2,e}>|\vec{r}|>r_{2,i}\), el potencial es constante y por lo tanto, \(V(r)=\frac{\sigma_{2,e}r_{2,e}}{\varepsilon_0}\log(\frac{r_{3}}{r_{2,e}})\)
  • Si \(r_{2,i}>|\vec{r}|>r_{1}\), \(\lambda(r)=2\pi r_1 \sigma_1\). El potencial tendrá dos contribuciones, una de la integral hasta \(r_{2,i}\), dada por \(V(r_{2,i})=V(r_{2,e})\) y otra entre \(r_{2,i}\) y \(|\vec{r}|\). En términos de la variable de integración, \(u\in (\frac{|r|}{r_{2,i}},1)\). Obtenemos entonces \[\begin{eqnarray*} V(|\vec{r}|)&=& \int_{0}^{\frac{|r|}{r_{2,i}}} \frac{\lambda(|\vec{r}|/u)}{2\pi\varepsilon_0 |r|} \frac{1}{u}du+ \int_{\frac{|r|}{r_{2,i}}}^{1} \frac{\lambda(|\vec{r}|/u)}{2 \pi\varepsilon_0 |r|} \frac{1}{u}du\\ &=&\frac{\sigma_{2,e}r_{2,e}\log\left(\frac{r_3}{r_{2,e}}\right) + \sigma_{1}r_{1}\log\left(\frac{r_{2,i}}{|\vec{r}|}\right)}{\varepsilon_0} \end{eqnarray*}\]

  • Finalmente, si \(|\vec{r}|<r_1\), el potencial debe ser constante, por lo que \[ V(|\vec{r}|)= \frac{\sigma_{2,e}r_{2,e}\log\left(\frac{r_3}{r_{2,e}}\right)+\sigma_{1}r_{1}\log\left(\frac{r_{2,i}}{r_1}\right)}{\varepsilon_0} \]

\[ V(\vec{r})=\left\{ \begin{array}{l r} \frac{\sigma_{2,e}r_{2,e}\log\left(\frac{r_3}{r_{2,e}}\right)+ \sigma_{1}r_{1}\log\left(\frac{r_{2,i}}{r_1}\right)}{\varepsilon_0} & r<r_1\\ \frac{\sigma_{2,e}r_{2,e}\log\left(\frac{r_3}{r_{2,e}}\right) + \sigma_{1}r_{1}\log\left(\frac{r_{2,i}}{|\vec{r}|}\right)}{\varepsilon_0} & r_{2,i}>r>r_1\\ \frac{\sigma_{2,e}r_{2,e}\log\left(\frac{r_3}{r_{2,e}}\right)}{\varepsilon_0} & r_{2,e}>r>r_{2,i}\\ \frac{\sigma_{2,e}r_{2,e}\log\left(\frac{r_3}{|\vec{r}|}\right)}{\varepsilon_0} & r_3>r>r_{2,e}\\ 0 & r>r_3 \end{array} \right. \]

  • Obsérvese que si \(r_3\rightarrow \infty\), el potencial diverge como \(\log(\frac{r_3}{r_{2,e}})\). Sin embargo, siempre es posible definir un referencial en algún \(r_3^{'}\) finito.
  • Si \(r_1 \rightarrow 0\), manteniendo \(Q_1\) finita, \(\lambda(r)\) se mantendrá finita (pero no \(\sigma\)). En tal caso, vemos que tampoco es posible definir el referncial para el potencial en \(r=0\).
  • El modelo de alambre infinito es adecuado si exploramos el sistema cerca de su centro, a distancias \(r\ll L\). A distancias \(r\gg L\), el sistema se comportará como una carga puntual de valor \(Q_1+Q_2+Q_3\).

Capacidad

Como el sistema consiste de tres conductores, tiene sentido evaluar su capacidad respecto a cualquiera de los posibles pares:

  • \(C_{12}\): \(\Delta V_{12}=\frac{\sigma_{1}r_{1}\log\left(\frac{r_{2,i}}{r_1}\right)}{\varepsilon_0}\), \(Q_1=-Q_2=2\pi r_{1}\sigma_{1} L\), luego, \[ C_{12}=\frac{Q_1}{\Delta V_{23}}=\frac{2\pi L \varepsilon_0}{\log\left(\frac{r_{2,i}}{r_1}\right)} \]

  • \(C_{23}\): \(\Delta V_{23}=\frac{\sigma_{2,e}r_{2,e}\log\left(\frac{r_{3}}{r_{2,e}}\right)}{\varepsilon_0}\), \(Q_2=-Q_3=2\pi r_{2,e}\sigma_{2,e} L\), luego, \[ C_{23}=\frac{Q_2}{\Delta V_{23}}=\frac{2\pi L \varepsilon_0}{\log\left(\frac{r_{3}}{r_{2,e}}\right)} \]

  • \(C_{13}\): \(\Delta V_{13}=\Delta V_{12}+\Delta V_{23}\), \(Q_2=0\rightarrow \sigma_{2,e}r_{2,e}=-\sigma_{2,i}r_{2,i}=\sigma_{1}r_{1}\), \(Q_3=-Q_1\). \[ C_{13}=\frac{Q_1}{\Delta V_{13}}= \frac{2\pi L \varepsilon_0}{\log\left(\frac{r_{3}}{r_{2,e}}\right) + \log\left(\frac{r_{2,i}}{r_{1}}\right)} =\frac{1}{\frac{1}{C_{12}}+\frac{1}{C_{23}}} \]

Esto es, el sistema se comporta como una combinación serie, del capacitor \(1-2\) y el capacitor \(2-3\).

Resolución sistemática de circuitos de CC

Combinación de F.E.M. en serie y paralelo

  • Si en una rama de un circuito encontramos varios generadores, de FEM’s \({\cal E}_k\) y resistencias internas \(R_{int,k}\) es posible remplazarlos por un único generador equivalente, cuya F.E.M. \({\cal E}_{eq}=\sum {\cal E}_{k}\) , y de resistencia interna \(R_{int,eq}=\sum_k R_{int,k}\).
  • Por otro lado, en una combinación paralelo, \(R_{int,eq}=\left(\sum_k R^{-1}_{int,k}\right)^{-1}\) y \({\cal E}_{equiv}=R_{int,eq}\left(\sum_k \frac{{\cal E}_{k}}{{\cal R}_{int,k}}\right)\).

(Probar como ejercicio)

Theorema de Thévenin

  • En general, dado un par de nodos en cualquier combinación de resistencias y FEM’s en estado estacionario, es posible encontrar un circuito equivalente compuesto solamente por una resistencia y una FEM en serie. Este resultado se conoce como Theorema de Thévenin
  • Este teorema es una consecuencia de la relación lineal entre tensión y corriente en las resistencias, y deja de ser válido si se introducen elementos no lineales (por ejemplo, diodos).
  • En las condiciones del teorema, la energía entregada/absorbida por el circuito equivalente es igual a la suma algebráica de potencias entregadas/absorbidas por sus componentes.

  • Como herramienta de cálculo, el teorema es muy útil cuando se conoce de antemano una expresión sencilla para los parámetros equivalentes.
  • Desde el punto de vista práctico, permite trabajar con circuitos complicados como cajas negas que se caracterizan por dos parámetros medibles directamente.

Procedimiento para la resolución de circuitos de C.C.

  • Construir un circuito equivalente en el que eliminamos todas las ramas que incluyan capacitores.
  • Aplicar reducciones serie y paralelo para las resistencias y F.E.M.’s presentes en el sistema, en forma iterativa.
  • Plantear las ecuaciones de Kirchhoff para el sistema reducido, y determinar las corrientes en cada rama.
  • Recuperar las corrientes y tensiones correspondientes a cada circuito equivalente reducido.
  • Determinar las diferencias de potencial en los capacitores, y calcular sus cargas.
  • Calcular las potencias disipadas o entregadas de acuerdo con la expresión \(P=\Delta V i\). Con \(\Delta V\) medido en el sentido de la corriente \(|P|\) es potencia entregada al circuito si \(P>0\), y disipada si \(P<0\).

Ejemplo