Taller 1

Nombre: Estefany Carolina Cely Rodriguez

C.C.: 1013659975

Electrodinamica I

Lista de ejercicios

Griffiths: An Introduction to Electrodynamics

2.15

Se tiene Un cascarón esférico con densidad de carga

ρ=kr2(arb)\begin{equation}\rho=\frac{k}{r^2} \quad(a \leq r \leq b)\end{equation}

Encontrar el campo electrico en las tres regiones:(i) r<ar<a, (ii) a<r<ba<r<b, (iii) r>br>b.

Grafique E|\mathbf{E}| en funcion de rr.


Sol:

La ley de Gauss para el campo electrico en el vacio:

E=ρϵ0\nabla \cdot \mathbf{E}=\frac{\rho}{\epsilon_0}

Al ser simetria esferica, el termino de divergencia es suficiente. De este modo al integrar ambos lados sobre el volumen de una superficie gaussiana de radio R dividimos el problema en 3 casos:

x02+y02+z02r2EdV0=x02+y02+z02r2ρϵ0dV0=1ϵ0x02+y02+z02r2ρdV0={1ϵ00π02π0r(0)(r02sinθ0dr0dϕ0dθ0) if r<a1ϵ00π02πarkr02(r02sinθ0dr0dϕ0dθ0) if a<r<b1ϵ00π02πabkr02(r02sinθ0dr0dϕ0dθ0) if r>b\begin{aligned}\iiint_{x_0^2+y_0^2+z_0^2 \leq r^2} \nabla \cdot \mathbf{E} d V_0 &=\iiint \int_{x_0^2+y_0^2+z_0^2 \leq r^2} \frac{\rho}{\epsilon_0} d V_0 \\=& \frac{1}{\epsilon_0} \iiint \int_{x_0^2+y_0^2+z_0^2 \leq r^2} \rho d V_0 \\=& \begin{cases}\frac{1}{\epsilon_0} \int_0^\pi \int_0^{2 \pi} \int_0^r(0)\left(r_0^2 \sin \theta_0 d r_0 d \phi_0 d \theta_0\right) & \text { if } r<a \\\frac{1}{\epsilon_0} \int_0^\pi \int_0^{2 \pi} \int_a^r \frac{k}{r_0^2}\left(r_0^2 \sin \theta_0 d r_0 d \phi_0 d \theta_0\right) & \text { if } a<r<b \\\frac{1}{\epsilon_0} \int_0^\pi \int_0^{2 \pi} \int_a^b \frac{k}{r_0^2}\left(r_0^2 \sin \theta_0 d r_0 d \phi_0 d \theta_0\right) & \text { if } r>b\end{cases}\end{aligned}

Aplicando el teorema de la divergencia y evaluando:

x02+y02+z02=r2EdS0={0 if r<akϵ0(0πsinθ0dθ0)(02πdϕ0)(ardr0) if a<r<bkϵ0(0πsinθ0dθ0)(02πdϕ0)(abdr0) if r>b\oiint_{x_0^2+y_0^2+z_0^2=r^2} \mathbf{E} \cdot d \mathbf{S}_0= \begin{cases}0 & \text { if } r<a \\ \frac{k}{\epsilon_0}\left(\int_0^\pi \sin \theta_0 d \theta_0\right)\left(\int_0^{2 \pi} d \phi_0\right)\left(\int_a^r d r_0\right) & \text { if } a<r<b \\ \frac{k}{\epsilon_0}\left(\int_0^\pi \sin \theta_0 d \theta_0\right)\left(\int_0^{2 \pi} d \phi_0\right)\left(\int_a^b d r_0\right) & \text { if } r>b\end{cases}

Es un sistema con simetria esferica, por lo que el campo electrico debe ser totalmente radial, de la forma E=E(r)r^\mathbf{E}=E(r) \hat{\mathbf{r}}, y se orienta respecto al vector unitario dSd \mathbf{S} que sale de la superficie gaussiana:

r02=r2[E(r0)r^0](r^0dS0)={0 if r<akϵ0(2)(2π)(ra) if a<r<kϵ0(2)(2π)(ba) if r>b\oiint_{r_0^2=r^2}\left[E\left(r_0\right) \hat{\mathbf{r}}_0\right] \cdot\left(\hat{\mathbf{r}}_0 d S_0\right)= \begin{cases}0 & \text { if } r<a \\ \frac{k}{\epsilon_0}(2)(2 \pi)(r-a) & \text { if } a<r< \\ \frac{k}{\epsilon_0}(2)(2 \pi)(b-a) & \text { if } r>b\end{cases}
r02=r2[E(r0)r^0](r^0dS0)={0 if r<akϵ0(2)(2π)(ra) if a<r<kϵ0(2)(2π)(ba) if r>b\oiint_{r_0^2=r^2}\left[E\left(r_0\right) \hat{\mathbf{r}}_0\right] \cdot\left(\hat{\mathbf{r}}_0 d S_0\right)= \begin{cases}0 & \text { if } r<a \\ \frac{k}{\epsilon_0}(2)(2 \pi)(r-a) & \text { if } a<r< \\ \frac{k}{\epsilon_0}(2)(2 \pi)(b-a) & \text { if } r>b\end{cases}

E(r)E(r) es una cantidad constante sobre la superficie por lo que,

E(r)r0=rdS0={0 if r<a4πkϵ0(ra) if a<r<b4πkϵ0(ba) if r>bE(r) \oiint_{r_0=r} d S_0= \begin{cases}0 & \text { if } r<a \\ \frac{4 \pi k}{\epsilon_0}(r-a) & \text { if } a<r<b \\ \frac{4 \pi k}{\epsilon_0}(b-a) & \text { if } r>b\end{cases}

Que al evaluar resulta en:

E(r)(4πr2)={0 if r<a4πkϵ0(ra) if a<r<b4πkϵ0(ba) if r>bE(r)\left(4 \pi r^2\right)= \begin{cases}0 & \text { if } r<a \\ \frac{4 \pi k}{\epsilon_0}(r-a) & \text { if } a<r<b \\ \frac{4 \pi k}{\epsilon_0}(b-a) & \text { if } r>b\end{cases}

2.16

💡
A long coaxial cable (Fig. 2.26) carries a uniform volume charge density ρ\rho on the inner cylinder (radius aa ), and a uniform surface charge density on the outer cylindrical shell (radius bb ). This surface charge is negative and of just the right magnitude so that the cable as a whole is electrically neutral. Find the electric field in each of the three regions: (i) inside the inner cylinder (s<a)(s<a), (ii) between the cylinders (a<s<b)(a<s<b), (iii) outside the cable (s>b)(s>b). Plot E|\mathbf{E}| as a function of ss.

Sabiendo que el campo eléctrico producido por el cilindro en su interior es radial, entonces usaremos la ley de Gauss, usando un volumen cilíndrico, de radio s<as<a y longitud ll, cómo se ilustra en la Figura 2, se integrará en el volumen.

E=ρ(r)εoEd3x=ρ(r)εod3x\begin{aligned}\vec{\nabla} \cdot \vec{E} &=\frac{\rho(\vec{r})}{\varepsilon_o} \\\int \vec{\nabla} \cdot \vec{E} d^3 x &=\int \frac{\rho(\vec{r})}{\varepsilon_o} d^3 x\end{aligned}

Por el teorema de la divergencia se tiene que:

Eda=ρ(r)εod3x\int \vec{E} \cdot d \vec{a}=\int \frac{\rho(\vec{r})}{\varepsilon_o} d^3 x

Ahora teniendo en cuenta que el campo eléctrico es radial, solo tiene flujo en la superficie lateral del cilindro y no en las tapas, además, es constante para un mismo radio y diferente angulo φ\varphi, por lo que:

Eda=ρ(r)εod3x|\vec{E}| \int d a=\int \frac{\rho(\vec{r})}{\varepsilon_o} d^3 x

Recordando que ρ(r)=ρ\rho(\vec{r})=\rho, y que el área de la primera integral es la superficie lateral del cilindro, entonces:

Eda=E2πlsρ(r)εod3x=0l02π0sρεordrdφdz=ρεo2πls22\begin{aligned}|\vec{E}| \int d a &=|\vec{E}| 2 \pi l s \\\int \frac{\rho(\vec{r})}{\varepsilon_o} d^3 x &=\int_0^l \int_0^{2 \pi} \int_0^s \frac{\rho}{\varepsilon_o} r d r d \varphi d z=\frac{\rho}{\varepsilon_o} 2 \pi l \frac{s^2}{2}\end{aligned}

Por lo que:

E2πls=ρ2πls22εoE=ρs2εoE=ρs2εor^\begin{aligned}|\vec{E}| 2 \pi l s &=\rho 2 \pi l \frac{s^2}{2 \varepsilon_o} \\|\vec{E}| &=\frac{\rho s}{2 \varepsilon_o} \\\vec{E} &=\frac{\rho s}{2 \varepsilon_o} \hat{r}\end{aligned}

(iii)(a<s<b)(iii) (a<s<b)

Al igual que en el caso anterior, integraremos en un volumen cilíndrico, que esta vez tendrá un radio de a<s<ba<s<b (Figura 3), esto implica que la integral de volumen de la densidad de carga no irá hasta s, sino que irá hasta a, pues es donde termina. de modo que se tendrá:

Eda=E2πlsρ(r)εod3x=0l02π0aρεordrdφdz=ρεo2πla22\begin{aligned}|\vec{E}| \int d a &=|\vec{E}| 2 \pi l s \\\int \frac{\rho(\vec{r})}{\varepsilon_o} d^3 x &=\int_0^l \int_0^{2 \pi} \int_0^a \frac{\rho}{\varepsilon_o} r d r d \varphi d z=\frac{\rho}{\varepsilon_o} 2 \pi l \frac{a^2}{2}\end{aligned}

Con lo que el campo eléctrico tendrá la forma:

E2πls=ρ2πla22εoE=ρa22sεoE=ρa22sεor^\begin{aligned}|\vec{E}| 2 \pi l s &=\rho 2 \pi l \frac{a^2}{2 \varepsilon_o} \\|\vec{E}| &=\frac{\rho a^2}{2 s \varepsilon_o} \\\vec{E} &=\frac{\rho a^2}{2 s \varepsilon_o} \hat{r}\end{aligned}

(s>b)(s>b)

Ahora la superficie en la que se integrará es la que se muestra en la Figura 4, con lo cual se tendrá ahora en cuenta el tuvo esferico que recubre el cable coaxial; es de recordar que el enunciado dice que la carga de este ultimo es negativa y de tal magnitud que el cable en conjunto es eléctricamente neutro, es decir, que si:

Vρd3x=Qenc\int_V \rho d^3 x=Q_{e n c}

Entonces:

Vσda=Qenc\int_V \sigma d a=-Q_{e n c}

Por lo que:

ρ(r)εod3x=1εo[Vρd3x+Vσda]=1εo[QencQenc]=0\begin{aligned}\int \frac{\rho(\vec{r})}{\varepsilon_o} d^3 x &=\frac{1}{\varepsilon_o}\left[\int_V \rho d^3 x+\int_V \sigma d a\right] \\&=\frac{1}{\varepsilon_o}\left[Q_{e n c}-Q_{e n c}\right]=0\end{aligned}

Por lo que el campo eléctrico es nulo:

E=0\vec{E}=0

Con esto se tiene que el campo eléctrico del cable coaxial es:

E(s)={12εoρsr^s<aρa22sεor^a<s<b0b<s\vec{E}(s)=\left\{\begin{array}{cc}\frac{1}{2 \varepsilon_o} \rho s \hat{r} & s<a \\\frac{\rho a^2}{2 s \varepsilon_o} \hat{r} & a<s<b \\0 & b<s\end{array}\right.

Al graficar La magnitud del campo eléctrico E|\vec{E}| en función de ss, se tiene lo ilustrado en la Figura 5.

2.17

💡
An infinite plane slab, of thickness 2d2 d, carries a uniform volume charge density ρ\rho (Fig. 2.27). Find the electric field, as a function of $$$y$, where y=0y=0 at the center. Plot EE versus yy, calling EE positive when it points in the +y+y direction and negative when it points in the y-y direction.

Una de las ecuaciones que gobiernan el campo electrico en el vacio es la ley de Gauss:

E=ρϵ0\nabla \cdot \mathbf{E}=\frac{\rho}{\epsilon_0}

Usando la simetria del plano xzxz solo necesitamos informacion de la divergencia, ademas tenemos unicamente la componente yy del campo electrico E=E(y)y^\mathbf{E}=E(y) \hat{\mathbf{y}}. Decimos que hay igual cantidad de carga hacia los extremos derecho e izquierdo que define el plano, de modo que el campo electrico en y=0y=0 es nulo.

Integrando ambos lados sobre el volumen de una superficie gaussiana de largo LL ancho yy y altura HH se consideran los siguientes casos:

  1. 0<y<d0<y<d
  1. y>dy>d

La carga encerrada es producto de la densidad de carga y el volumenm de modo que:

0H0y0LE(dx0dy0dz0)={0H0y0Lρϵ0(dx0dy0dz0) if 0<y<d0H0d0Lρϵ0(dx0dy0dz0) if y>d\int_0^H \int_0^y \int_0^L \nabla \cdot \mathbf{E}\left(d x_0 d y_0 d z_0\right)= \begin{cases}\int_0^H \int_0^y \int_0^L \frac{\rho}{\epsilon_0}\left(d x_0 d y_0 d z_0\right) & \text { if } 0<y<d \\ \int_0^H \int_0^d \int_0^L \frac{\rho}{\epsilon_0}\left(d x_0 d y_0 d z_0\right) & \text { if } y>d\end{cases}

Ahora aplicando el teorema de la divergencia a la expresion del lado izuqiedo:

EdS0={ρϵ0(0Hdz0)(0ydy0)(0Ldx0) if 0<y<dρϵ0(0Hdz0)(0ddy0)(0Ldx0) if y>d\oiint \mathbf{E} \cdot d \mathbf{S}_0= \begin{cases}\frac{\rho}{\epsilon_0}\left(\int_0^H d z_0\right)\left(\int_0^y d y_0\right)\left(\int_0^L d x_0\right) & \text { if } 0<y<d \\ \frac{\rho}{\epsilon_0}\left(\int_0^H d z_0\right)\left(\int_0^d d y_0\right)\left(\int_0^L d x_0\right) & \text { if } y>d\end{cases}

Debido a que el campo electrico tiene solo componente yy la integral de superficie puede ser calculada en dos caras de la caja:

0H0L[E(y0)y^0(y^0dx0dz0)]y0=0+0H0L[E(y0)y^0(y^0dx0dz0)]y0=y={ρϵ0(HLy) if 0<y<dρϵ0(HLd) if y>d\left.\int_0^H \int_0^L\left[E\left(y_0\right) \hat{\mathbf{y}}_0 \cdot\left(\hat{\mathbf{y}}_0 d x_0 d z_0\right)\right]\right|_{y_0=0}+\left.\int_0^H \int_0^L\left[E\left(y_0\right) \hat{\mathbf{y}}_0 \cdot\left(\hat{\mathbf{y}}_0 d x_0 d z_0\right)\right]\right|_{y_0=y}= \begin{cases}\frac{\rho}{\epsilon_0}(H L y) & \text { if } 0<y<d \\ \frac{\rho}{\epsilon_0}(H L d) & \text { if } y>d \\ & \end{cases}

Al evaluar los productos punto obtenemos:

0H0LE(0)dx0dz0+0H0LE(y)dx0dz0={ρϵ0(HLy) if 0<y<dρϵ0(HLd) if y>d\int_0^H \int_0^L E(0) d x_0 d z_0+\int_0^H \int_0^L E(y) d x_0 d z_0= \begin{cases}\frac{\rho}{\epsilon_0}(H L y) & \text { if } 0<y<d \\ \frac{\rho}{\epsilon_0}(H L d) & \text { if } y>d\end{cases}

El campo eléctrico es cero en la cara y=0y=0 y constante en y0=yy_0=y

0H0L(0)dx0dz0+E(y)0H0Ldx0dz0={ρϵ0(HLy) if 0<y<dρϵ0(HLd) if y>d\int_0^H \int_0^L(0) d x_0 d z_0+E(y) \int_0^H \int_0^L d x_0 d z_0= \begin{cases}\frac{\rho}{\epsilon_0}(H L y) & \text { if } 0<y<d \\ \frac{\rho}{\epsilon_0}(H L d) & \text { if } y>d\end{cases}

Evaluando las integrales:

E(y)(HL)={ρϵ0(HLy) if 0<y<dρϵ0(HLd) if y>dE(y)(H L)= \begin{cases}\frac{\rho}{\epsilon_0}(H L y) & \text { if } 0<y<d \\ \frac{\rho}{\epsilon_0}(H L d) & \text { if } y>d\end{cases}

Luego, al dividir ambos lados por HLH L para resolver E(y)E(y):

E(y)={ρϵ0y if 0<y<dρϵ0d if y>dE(y)= \begin{cases}\frac{\rho}{\epsilon_0} y & \text { if } 0<y<d \\ \frac{\rho}{\epsilon_0} d & \text { if } y>d\end{cases}

Debido a la simetria sobre el plano y=0y=0 y el campo es negativo para valores y<0y<0, obtenemos el campo electrico para <y<-\infty<y<\infty con los terminos impares de E(y)E(y)

E(y)={ρϵ0d if y<dρϵ0y if d<y<dρϵ0d if y>dE(y)= \begin{cases}-\frac{\rho}{\epsilon_0} d & \text { if } y<-d \\ \frac{\rho}{\epsilon_0} y & \text { if }-d<y<d \\ \frac{\rho}{\epsilon_0} d & \text { if } y>d\end{cases}

De modo que,

E(y)={ρϵ0dy^ if y<dρϵ0yy^ if d<y<d.ρϵ0dy^ if y>d\mathbf{E}(y)= \begin{cases}-\frac{\rho}{\epsilon_0} d \hat{\mathbf{y}} & \text { if } y<-d \\ \frac{\rho}{\epsilon_0} y \hat{\mathbf{y}} & \text { if }-d<y<d . \\ \frac{\rho}{\epsilon_0} d \hat{\mathbf{y}} & \text { if } y>d\end{cases}

Esta es la grafica de E(y)E(y) vs yy.

2.21

💡
Find the potential inside and outside a uniformly charged solid sphere whose radius is RR and whose total charge is qq. Use infinity as your reference point. Compute the gradient of VV in each region, and check that it yields the correct field. Sketch V(r)V(r).

Por definición: V(r)=sE(r)ds V(\vec{r})=\oint_s E(\vec{r}) \cdot d s  Mediante la ley de Gauss vamos a determinar el campo eléctrico en un punto dentro y fuera de la esfera.

ρ=Q43πR3E(r)=14πϵoQR3r \begin{aligned} \rho &=\frac{Q}{\frac{4}{3} \pi R^3} \\ E(\vec{r}) &=\frac{1}{4 \pi \epsilon_o} \frac{Q}{R^3} r \end{aligned} 
E(r)=14πϵoQr2 E(\vec{r})=\frac{1}{4 \pi \epsilon_o} \frac{Q}{r^2} 
Vo=r(14πϵoQr2)dr=14πϵoQrV_o=-\int_{\infty}^r\left(\frac{1}{4 \pi \epsilon_o} \frac{Q}{r^2}\right) d r=\frac{1}{4 \pi \epsilon_o} \frac{Q}{r}
E=VE=Vrr^=14πϵoQr2r^ \begin{gathered} E=-\nabla V \\ E=-\frac{\partial V}{\partial r} \hat{r}=\frac{1}{4 \pi \epsilon_o} \frac{Q}{r^2} \hat{r} \end{gathered} 
E=VE=Vrr^=Q4πϵo12R2rR2r^E=14πϵoQR3rr^\begin{gathered} E=-\nabla V \\ E=-\frac{\partial V}{\partial r} \hat{r}=\frac{Q}{4 \pi \epsilon_o} \frac{1}{2 R} \frac{2 r}{R^2} \hat{r} \\ E=\frac{1}{4 \pi \epsilon_o} \frac{Q}{R^3} r \quad \hat{r} \end{gathered}

2.23

💡
For the charge configuration of Prob. 2.15, find the potential at the center, using infinity as your reference point.

Sol:

Del problema 2.152.15 de Griffiths se obtiene el campo para las siguientes regiones, donde se a utilizado la ley de Gauss debido a la simétria esférica que posee el problema.

E2=k(ra)ϵ0r2r^ \vec{E}_2=\frac{k(r-a)}{\epsilon_0 r^2} \hat{r} 
E3=k(ba)ϵ0r2r^\vec{E}_3=\frac{k(b-a)}{\epsilon_0 r^2} \hat{r}

Este resultado implica que el potencial dentro del cascarón es constante, resultado esperado ya que el campo eléctrico en esta región es 0 , también mediante la relación E=V\vec{E}=-\nabla V se llegá a que el potencial es constante.


2.24

💡
For the configuration of Prob. 2.16, find the potential difference between a point on the axis and a point on the outer cylinder. Note that it is not necessary to commit yourself to a particular reference point if you use Eq. 2.222.22.

Estos campos fueron:

Edentro =ρr2ϵ0r^\mathrm{E}_{\text {dentro }}=\frac{\rho r}{2 \epsilon_0} \hat{r}
Eentre =ρa22ϵ0rr^\mathbf{E}_{\text {entre }} = \frac{\rho a^2}{2 \epsilon_0 r} \hat{r}
Efuera =0\mathbf{E}_{\text {fuera }}=0

Donde r es la distancia a la cual se encuentra en punto en el que se quiere conocer el campo. Tomando en cuenta que la componente que importa en el campo eléctrico es la componente radial, los vectores a y b solo contribuirán con la distancia al eje de los cilindros. Llamaremos d al radio del punto fuera del eje, y el punto en el eje tendrá distancia cero. Por tanto la Eq. 1 queda como:

V(d)V(0)=0dEdl=0dEdentro dlabEentre dl0dEfuera dlV(d)-V(0)=-\int_0^d \mathbf{E} \cdot d \mathbf{l}=-\int_0^d \mathbf{E}_{\text {dentro }} \cdot d \mathbf{l}-\int_a^b \mathbf{E}_{\text {entre }} \cdot d \mathbf{l}-\int_0^d \mathbf{E}_{\text {fuera }} \cdot d \mathbf{l}

Teniendo en cuenta las ecuaciones y que tanto el campo EE como el dldl son radiales:

V(d)V(0)=0aρrΣ0drabρa22r0)dr=ρ2r0a22ρa22r0lnba=μ24c0[1+2lnba]\begin{aligned} V(d)-V(0) &=-\int_0^a \frac{\rho r}{\Sigma_0} d r-\int_a^b \frac{\rho a^2}{\left.2 r_0\right)^{\prime}} d r \\ &=-\frac{\rho}{2 r_0} \frac{a^2}{2}-\frac{\rho a^2}{2 r_0} \ln \frac{b}{a} \\ &=\quad-\frac{\mu^2}{4 c_0}\left[1+2 \ln \frac{b}{a}\right] \end{aligned}

Que corresponde a la diferencia de potencial entre los dos puntos.


2.36


a) σa,σb,σR\sigma_a, \sigma_b, \sigma_R : Tomando como superficie gaussiana a la esfera de radio aa, justo integrando en el conductor (donde el campo es cero) se tiene:

EadS=Qencerrada ϵ0=0qa+qina =0\oint \vec{E}a \cdot d \vec{S}=\frac{Q{\text {encerrada }}}{\epsilon_0}=0 \longrightarrow q_a+q_{\text {ina }}=0

Siendo qina q_{\text {ina }} la carga inducida en la superficie del hueco esférico; por la geometría de la superficie, la carga qina q_{\text {ina }} se distribuye de manera homogénea en el área superficial asi:

σa=qa4πa2\sigma_a=\frac{-q_a}{4 \pi a^2}

Siguiendo el mismo análisis para el hueco esférico de radio b y carga interna qbq_b se tiene:

σb=qb4πb2\sigma_b=\frac{-q_b}{4 \pi b^2}

Tomando en cuenta que la carga total en el conductor es nula y que de haber cargas parciales, deben estar distribuidas en las superficies, tenemos:

0=qinu +qinb +qRqaqb+4πR2σR0=q_{\text {inu }}+q_{\text {inb }}+q_R \longrightarrow-q_a-q_b+4 \pi R^2 \sigma_R

Donde qinb q_{\text {inb }} es la carga inducida en la superficie de la esfera de radio bb; por lo tanto al despejar se tiene:

σR=qa+qb4πR2\sigma_R=\frac{q_{\mathrm{a}}+q_{\mathrm{b}}}{4 \pi R^2}

b) Er\vec{E}r con r>Rr>R : Sea la superficie gaussiana la esfera de radio R+rR+r, se hace la integral justo fuera del conductor conQencerrada =(qa+qb)\operatorname{con} Q{\text {encerrada }}=\left(q_a+q_b\right)  asi:

ErdS=qa+qbϵ0=E4π(R+r)2r^\oint \vec{E}_r \cdot d \vec{S}=\frac{q_a+q_b}{\epsilon_0}=|E| 4 \pi(R+r)^2 \hat{r}

Despejando tenemos:

Er=qa+qb4π(R+r)2ϵ0r˙\vec{E}_r=\frac{q_a+q_b}{4 \pi(R+r)^2 \epsilon_0} \dot{r}

c) Ea,Eb\vec{E}_a, \vec{E}_b : Al tomar la superficie gaussiana esférica de radio aa, justo antes del conductor, tenemos:

EadS=quϵ0Ea=qa4πr2ϵ0r^\oint \vec{E}_a \cdot d \vec{S}=\frac{q_u}{\epsilon_0} \longrightarrow \vec{E}_a=\frac{q_a}{4 \pi r^2 \epsilon_0} \hat{r}

Siguiendo el mismo proceso para la superficie gaussiana esférica de radio bb tenemos:

Eb=qb4πr2ϵ0r~\vec{E}_b=\frac{q_b}{4 \pi r^2 \epsilon_0} \tilde{r}

d) Fa,Fb\vec{F}_a, \vec{F}_b : Como F=qE\vec{F}=q \vec{E}, entonces estudiamos el campo eléctrico en la esfera de radio a. Sabemos que dentro de un cascarón conductor el campo eléctrico es solamente el debido a la carga interna a la superficie.

Luego Fab=Fba=0\vec{F}{a b}=\vec{F}{b a}=0; como en ausencia de la carga $q_a$ el campo dentro de la esfera es nulo, tenemos:

Fa=Fb=0\vec{F}_a=\vec{F}_b=0

e) Si se acerca una carga qcq_c a la superficie externa del conductor, ocurre lo siguiente: (a) Cambia σR\sigma_R. (b) Er\vec{E}_r cambia por la distribución no homogénea de cargas en la superficie externa. (c) $$\vec{E}_a$$ y Eb\vec{E}_b se mantienen constantes. (d) Fa,Fb\vec{F}_a, \vec{F}_b permanecen constantes.


2.35


Una esfera de radio R\mathrm{R} con una densidad de carga ρ(r)=kr\rho(r)=k r.encuentre la energia de esta configuración, Hallelo de dos formas distintas. Solución: Sabiendo que la ley de gauss es

Eda=Qinϵ0=1ϵ0ρ(r)dV\oint E \quad d a=\frac{Q_{i n}}{\epsilon_0}=\frac{1}{\epsilon_0} \int \rho(r) d V

Suponiendo una superficia gaussiana donde el campo eléctrico y el diferencial de área sean paralelos y reemplazando obtenemos que

Eda=4πr2E\oint E \quad d a=4 \pi r^2 E

y por el otro lado

1ϵ0ρ(r)dV=1ϵ0(kr)r2sinθdrdθdφ=4πϵ0r2dr={πc0kr4(r<R)πc0kR4(r>R)\begin{aligned} \frac{1}{\epsilon_0} \int \rho(r) d V &=\frac{1}{\epsilon_0} \int(k r) r^2 \sin \theta d r d \theta d \varphi \\ &=\frac{4 \pi}{\epsilon_0} \int r^2 d r= \begin{cases}\frac{\pi}{c_0} k r^4 & (r<R) \\ \frac{\pi}{c_0} k R^4 & (r>R)\end{cases} \end{aligned}

Igualando ambas partes y despejando el campo eléctrico llegamos a

E(r)={kr24ϵ0r^(r<R)R44ϵ0r2r^(r>R)E(r)= \begin{cases}\frac{k r^2}{4 \epsilon_0} \hat{r} & (r<R) \\ \frac{R^4}{4 \epsilon_0 r^2} \hat{r} & (r>R)\end{cases}

a) Metodo 1 Partiendo de la definición de la energía a partir del campo electrico

W=ϵ02E2dτW=\frac{\epsilon_0}{2} \int E^2 d \tau

Reemplazando tenemos

W=ϵ02E2dτ=ϵ020R(kr24ϵ0)24πr2dτ+ϵ02R(kR44ϵ0r2)24πr2dτ=4πϵ02(k4ϵ0)2{0Rr6dτ+R8R1r2dτ}=πk28ϵ0(R77+R7)=πk2R77ϵ0\begin{aligned} W &=\frac{\epsilon_0}{2} \int E^2 d \tau=\frac{\epsilon_0}{2} \int_0^R\left(\frac{k r^2}{4 \epsilon_0}\right)^2 4 \pi r^2 d \tau+\frac{\epsilon_0}{2} \int_R^{\infty}\left(\frac{k R^4}{4 \epsilon_0 r^2}\right)^2 4 \pi r^2 d \tau \\ &=4 \pi \frac{\epsilon_0}{2}\left(\frac{k}{4 \epsilon_0}\right)^2\left\{\int_0^R r^6 d \tau+R^8 \int_R^{\infty} \frac{1}{r^2} d \tau\right\}=\frac{\pi k^2}{8 \epsilon_0}\left(\frac{R^7}{7}+R^7\right) \\ &=\frac{\pi k^2 R^7}{7 \epsilon_0} \end{aligned}

b) Metodo 2 Sabemos que la energía también se puede hallar a partir del potencial de la siguiente forma

W=12ρVdτW=\frac{1}{2} \int \rho V \quad d \tau

Primero hallamos el potencial de esta configuración para (r<R)(r<R)

V(r)=rEdl=R(kR44ϵ0r2)drRr(kr24ϵ0)dr=k4ϵ0{R4(1r)R+r33Rr}=k3ϵ0(R3r34)\begin{aligned} V(r) &=-\int_{\infty}^r E \quad d l=-\int_{\infty}^R\left(\frac{k R^4}{4 \epsilon_0 r^2}\right) d r-\int_R^r\left(\frac{k r^2}{4 \epsilon_0}\right) d r \\ &=-\frac{k}{4 \epsilon_0}\left\{\left.R^4\left(-\frac{1}{r}\right)\right|_{\infty} ^R+\left.\frac{r^3}{3}\right|_R ^r\right\}=\frac{k}{3 \epsilon_0}\left(R^3-\frac{r^3}{4}\right) \end{aligned}

Reemplazando tenemos

W=12ρVdτ=120R(kr)[k3ϵ0(R3r34)]4πr2dr=2πk23ϵ00R(R3r3r64)dr=2πk23ϵ0{R3R4414R77}=πk2R723ϵ0(67)=πk2R77ϵ0\begin{aligned} W &=\frac{1}{2} \int \rho V \quad d \tau=\frac{1}{2} \int_0^R(k r)\left[\frac{k}{3 \epsilon_0}\left(R^3-\frac{r^3}{4}\right)\right] 4 \pi r^2 d r=\frac{2 \pi k^2}{3 \epsilon_0} \int_0^R\left(R^3 r^3-\frac{r^6}{4}\right) d r \\ &=\frac{2 \pi k^2}{3 \epsilon_0}\left\{R^3 \frac{R^4}{4}-\frac{1}{4} \frac{R^7}{7}\right\}=\frac{\pi k^2 R^7}{2 \cdot 3 \epsilon_0}\left(\frac{6}{7}\right)=\frac{\pi k^2 R^7}{7 \epsilon_0} \end{aligned}

Arfken

Ejercicio 3.10.4

Dado que

V(q1,q2,q3)=1h1h2h3[q1(V1h2h3)+q2(V2h3h1)+q3(V3h1h2)]\begin{equation} \begin{split} \boldsymbol{\nabla} \cdot \mathbf{V}\left(q_1, q_2, q_3\right) &=\frac{1}{h_1 h_2 h_3}\left[\frac{\partial}{\partial q_1}\left(V_1 h_2 h_3\right)+\frac{\partial}{\partial q_2}\left(V_2 h_3 h_1\right)+\frac{\partial}{\partial q_3}\left(V_3 h_1 h_2\right)\right]\\ \end{split} \end{equation}

Entonces si V=e^1\mathbf{V} = \mathbf{\hat e}_1:

V(q1,q2,q3)=1h1h2h3[q1(V1h2h3)+q2(V2h3h1)+q3(V3h1h2)]=1h1h2h3[q1(e^1h2h3)]=1h1h2h3q1(h2h3) \begin{aligned} \boldsymbol{\nabla} \cdot \mathbf{V}\left(q_1, q_2, q_3\right) &=\frac{1}{h_1 h_2 h_3}\left[\frac{\partial}{\partial q_1}\left(V_1 h_2 h_3\right)+\frac{\partial}{\partial q_2}\left(V_2 h_3 h_1\right)+\frac{\partial}{\partial q_3}\left(V_3 h_1 h_2\right)\right]\\ &=\frac{1}{h_1 h_2 h_3}\left[\frac{\partial}{\partial q_1}\left(|\mathbf{\hat e}_1|h_2 h_3\right)\right]\\ &=\frac{1}{h_1 h_2 h_3}\frac{\partial}{\partial q_1}\left( h_2 h_3\right)\\ \end{aligned}

De igual forma como:

×B=1h1h2h3e^1h1e^2h2e^3h3q1q2q3h1B1h2B2h3B3\begin{equation}\boldsymbol{\nabla} \times \mathbf{B}=\frac{1}{h_1 h_2 h_3}\left|\begin{array}{ccc}\hat{\mathbf{e}}_1 h_1 & \hat{\mathbf{e}}_2 h_2 & \hat{\mathbf{e}}_3 h_3 \\\frac{\partial}{\partial q_1} & \frac{\partial}{\partial q_2} & \frac{\partial}{\partial q_3} \\h_1 B_1 & h_2 B_2 & h_3 B_3 \end{array}\right|\end{equation}

Entonces si B=e^1\mathbf{B} = \mathbf{\hat e}_1:

×B=1h1h2h3e^1h1e^2h2e^3h3q1q2q3h100=1h1[e^21h3h1q3e^31h2h1q2]\begin{aligned}\boldsymbol{\nabla} \times \mathbf{B}&=\frac{1}{h_1 h_2 h_3}\left|\begin{array}{ccc}\hat{\mathbf{e}}_1 h_1 & \hat{\mathbf{e}}_2 h_2 & \hat{\mathbf{e}}_3 h_3\\\frac{\partial}{\partial q_1} & \frac{\partial}{\partial q_2} & \frac{\partial}{\partial q_3} \\h_1 & 0 & 0 \end{array}\right|\\ &=\frac{1}{h_1}\left[\hat{\mathbf{e}}_2 \frac{1}{h_3} \frac{\partial h_1}{\partial q_3}-\hat{\mathbf{e}}_3 \frac{1}{h_2} \frac{\partial h_1}{\partial q_2}\right] \end{aligned}

Ejercicio 3.10.8

e^ρ=e^xcosφ+e^ysinφ\hat{e}_{\rho} = \hat{e}_{x} \cos \varphi + \hat{e}_{y} \sin \varphi
e^φ=e^xsinφ+e^ycosφ\hat{e}_{\varphi} = -\hat{e}_{x} \sin \varphi + \hat{e}_{y} \cos \varphi

De las ecuaciones:

e^iqj=e^j1hihjqi,ij\begin{equation}\frac{\partial \hat{\mathbf{e}}_i}{\partial q_j}=\hat{\mathbf{e}}_j \frac{1}{h_i} \frac{\partial h_j}{\partial q_i}, \quad i \neq j\end{equation}
e^iqi=jie^j1hjhiqj\begin{equation}\frac{\partial \hat{\mathbf{e}}_i}{\partial q_i}=-\sum_{j \neq i} \hat{\mathbf{e}}_j \frac{1}{h_j} \frac{\partial h_i}{\partial q_j}\end{equation}

Obtenemos que

e^ρφ=e^φ1hρhφρ=e^φρρ=e^φ\begin{equation*}\frac{\partial \hat{\mathbf{e}}_\rho}{\partial \varphi} = \mathbf{\hat e}_{\varphi}\frac{1}{h_{\rho}}\frac{ \partial h_{\varphi}}{\partial \rho} = \mathbf{\hat e}_{\varphi}\frac{ \partial \rho}{\partial \rho} = \mathbf{\hat e}_{\varphi}\end{equation*}

pues hρ=1,hφ=ρ,hz=1h_\rho=1, \quad h_{\varphi}=\rho, \quad h_z=1, ademas

e^φφ=e^ρ1hρhφρe^z1hzhφz=e^ρρρe^zρz=e^ρ\begin{aligned} \frac{\partial \hat{\mathbf{e}}_{\varphi}}{\partial \varphi} &= -\hat{\mathbf{e}}_{\rho} \frac{1}{h_{\rho}} \frac{\partial h_{\varphi}}{\partial \rho}-\hat{\mathbf{e}}_{z} \frac{1}{h_z} \frac{\partial h_{\varphi}}{\partial z}\\ &= -\hat{\mathbf{e}}_{\rho} \frac{\partial \rho}{\partial \rho}-\hat{\mathbf{e}}_{z} \frac{\partial \rho}{\partial z}= -\hat{\mathbf{e}}_{\rho} \end{aligned}

Ejercicio 3.10.10

r=xe^x+ye^y+ze^zr = x \hat{e}_x+ y\hat{e}_y + z \hat{e}_z

(a) Como x=ρcos(φ)x = \rho \cos(\varphi)  y y=ρsin(φ)y = \rho \sin(\varphi), z=zz=z

r=ρcosφ(e^ρcosφe^φsinφ)+ρsinφ(e^ρsinφe^φcosφ)+ze^z=ecos2φe^ρρcosφsinφe^ϕ+ρsin2φe^ρ+ρcosφsinφe^φ+ze^z=ρe^ρcos2φ+sin2φ)+ze^z=e^ρρ+e^zρr = \rho \cos \varphi (\hat{e}_{\rho}\cos\varphi- \hat{e}_{\varphi}\sin \varphi)+\rho \sin \varphi (\hat{e}_{\rho}\sin\varphi- \hat{e}_{\varphi}\cos \varphi)+z\hat{e}_{z}\\ =e \cos^2 \varphi \hat{e}_{\rho} - \rho \cos \varphi \sin\varphi \hat{e}_{\phi}+ \rho \sin^2 \varphi \hat{e}_\rho+ \rho \cos \varphi \sin\varphi \hat{e}_{\varphi}+z \hat{e}_z \\=\rho \hat{e}_{\rho} \cos^2 \varphi + \sin^2 \varphi) + z \hat{e}z\\ = \hat{e}_\rho \rho + \hat{e}_z \rho

(b)

r=(e^ρρ+e^φ1ρφ+e^zz)(e^ρρ+e^zz)=e^ρρ(e^ρρ+e^zz)+e^φ1ρφ(e^ρρ+e^zz)+e^zz(e^ρρ+e^zz)=e^ρe^ρ+e^φ1ρe^φρ+e^e^z=3\begin{aligned} \mathbf \nabla \cdot \mathbf r &= \left(\hat{\mathbf{e}}_\rho \frac{\partial}{\partial \rho}+\hat{\mathbf{e}}_{\varphi} \frac{1}{\rho} \frac{\partial}{\partial \varphi}+\hat{\mathbf{e}}_z \frac{\partial}{\partial z}\right)\cdot (\hat{\mathbf{e}}_\rho \rho+\hat{\mathbf{e}}_z z)\\ &= \hat{\mathbf{e}}_\rho \frac{\partial}{\partial \rho}\cdot (\hat{\mathbf{e}}_\rho \rho+\hat{\mathbf{e}}_z z)+\hat{\mathbf{e}}_{\varphi} \frac{1}{\rho} \frac{\partial}{\partial \varphi} \cdot (\hat{\mathbf{e}}_\rho \rho+\hat{\mathbf{e}}_z z)+\hat{\mathbf{e}}_z \frac{\partial}{\partial z}\cdot (\hat{\mathbf{e}}_\rho \rho+\hat{\mathbf{e}}_z z)\\ &= \hat{\mathbf{e}}_\rho \cdot \hat{\mathbf{e}}_\rho +\hat{\mathbf{e}}_{\varphi} \cdot \frac{1}{\cancel \rho} \hat{\mathbf{e}}_\varphi \cancel \rho+\hat{\mathbf{e}}\cdot\hat{\mathbf{e}}_z = 3\\ \end{aligned}
×V=1ρe^ρρe^φe^zρφzρ0z\begin{aligned} \mathbf \nabla \times \mathbf V &=\frac{1}{\rho}\left|\begin{array}{ccc} \hat{\mathbf{e}}_\rho & \rho \hat{\mathbf{e}}_{\varphi} & \hat{\mathbf{e}}_z \\\\ \frac{\partial}{\partial \rho} & \frac{\partial}{\partial \varphi} & \frac{\partial}{\partial z} \\\\ \rho & 0 & z \end{array}\right|\\ \end{aligned}

Ahora como Vρ=ρV_{\rho} = \rho, Vφ=0V\varphi=0, Vz=zV_z=z

×r=1ρ((z(z)z(ρ(0)))e^ρ(ρ(ρ(0))z(ρ))ρe^φ+(ρ(ρ(0))φ(ρ))e^z)\nabla \times r = \frac{1}{\rho} \left( \left( \frac{\partial}{\partial z}(z) - \frac{\partial}{\partial z}\left( \rho(0) \right) \right)\hat{e}_\rho - \left( \frac{\partial}{\partial \rho} \left( \rho(0) \right) -\frac{\partial}{\partial z}(\rho) \right) \rho \hat{e}_{\varphi} + \left( \frac{\partial}{\partial \rho} (\rho(0)) - \frac{\partial}{\partial \varphi}(\rho)\right) \hat{e}_z \right)
×r=0\nabla \times r = 0

Ejercicio 3.10.12

(a) Dado que ω=ωe^z\vec\omega = \omega \mathbf{\hat e_z} y r=e^ρρ+e^zz\mathbf{r}=\hat{\mathbf{e}}_\rho \rho+\hat{\mathbf{e}}_z z entonces:

v=ω×r=ωe^z×(e^ρρ+e^zz)=ωe^z×e^ρρ=ρωe^φ\begin{aligned} \mathbf{v} &= \vec \omega \times \mathbf{r}\\ & = \omega \mathbf{\hat e_z} \times ( \hat{\mathbf{e}}_\rho \rho+\hat{\mathbf{e}}_z z)\\ & = \omega \mathbf{\hat e_z} \times \hat{\mathbf{e}}_\rho\rho\\ & = \rho\omega \mathbf{\hat e_\varphi} \end{aligned}

(b)

×v=1ρe^ρρe^φe^zρφz0ρ2ω0=2ρωe^z=2ρω=×V=ρωz^\begin{aligned} \mathbf \nabla \times \mathbf v &= \frac{1}{\rho}\left|\begin{array}{ccc}\hat{\mathbf{e}}_\rho & \rho \hat{\mathbf{e}}_{\varphi} & \hat{\mathbf{e}}_z \\\\ \frac{\partial}{\partial \rho} & \frac{\partial}{\partial \varphi} & \frac{\partial}{\partial z} \\\\ 0 & \rho^2 \omega & 0\end{array}\right|\\ &=2\rho\omega \mathbf{\hat e_z} = 2\rho \vec \omega\\ &= \nabla \times V = \rho \omega \hat{z} \end{aligned}

Ejercicio 3.10.13

v=drdt=ddt[e^ρρ+e^zz]=ddt([e^xcosφ(t)+e^ysinφ(t)]ρ(t)+e^zz(t))=([e^xdφ(t)dtsinφ(t)+dφ(t)dte^ycosφ(t)]ρ(t)+dρ(t)dte^ρ+e^zdz(t)dt)=ρφ˙e^φ+ρ˙e^ρ+z˙e^z\begin{aligned} \mathbf v &= \frac{d \mathbf r}{dt} \\ &= \frac{d }{dt}\left[\hat{\mathbf{e}}_\rho \rho+\hat{\mathbf{e}}_z z\right]\\ &= \frac{d }{dt}\left(\left[\hat{\mathbf{e}}_x \cos \varphi(t)+\hat{\mathbf{e}}_y \sin \varphi(t)\right] \rho(t)+\hat{\mathbf{e}}_z z(t)\right)\\ &= \left(\left[-\hat{\mathbf{e}}_x \frac{d \varphi(t)}{dt}\sin \varphi(t)+\frac{d \varphi(t)}{dt}\hat{\mathbf{e}}_y \cos \varphi(t)\right] \rho(t) + \frac{d \rho(t)}{dt}\mathbf{\hat e_\rho}+\hat{\mathbf{e}}_z \frac{d z(t)}{dt}\right)\\ &=\rho\dot{\varphi} \mathbf{\hat e_\varphi} + \dot \rho \mathbf{\hat e_\rho}+\dot z\hat{\mathbf{e}}_z \end{aligned}

Y ademas:

a=dvdt=ddt[ρφ˙e^φ+ρ˙e^ρ+z˙e^z]=ρ˙φ˙e^φ+ρφ¨e^φρφ˙2e^ρ+ρ¨e^ρ+ρ˙φ˙e^φ+z¨e^z=(ρ¨ρφ˙2)e^ρ+(2ρ˙φ˙+ρφ¨)e^φ+z¨e^z\begin{aligned} \mathbf a &= \frac{d \mathbf v}{dt} \\ &= \frac{d }{dt}\left[\rho\dot{\varphi} \mathbf{\hat e_\varphi} + \dot \rho \mathbf{\hat e_\rho}+\dot z\hat{\mathbf{e}}_z \right]\\ &= \dot \rho \dot{\varphi}\mathbf{\hat e_\varphi} + \rho \ddot \varphi \mathbf{\hat e_\varphi} - \rho \dot \varphi^2 \mathbf{\hat e_\rho} + \ddot \rho \mathbf{\hat e_\rho}+ \dot \rho \dot \varphi\mathbf{\hat e_\varphi}+\ddot z\hat{\mathbf{e}}_z \\ &= (\ddot \rho -\rho \dot \varphi^2) \mathbf{\hat e_\rho}+(2\dot \rho \dot{\varphi} + \rho \ddot \varphi )\mathbf{\hat e_\varphi}+\ddot z\hat{\mathbf{e}}_z \\ \end{aligned}

Ejercicio 3.10.15

×A=1ρe^ρρe^φe^zρφz00μI2πln(1ρ)=ρ(μI2πln(1ρ))e^φ=μI2πρ(ln(1ρ))e^φ=μI2π1ρe^φ\begin{aligned} \mathbf \nabla \times \mathbf A &= \frac{1}{\rho}\left|\begin{array}{ccc}\hat{\mathbf{e}}_\rho & \rho \hat{\mathbf{e}}_{\varphi} & \hat{\mathbf{e}}_z \\\\ \frac{\partial}{\partial \rho} & \frac{\partial}{\partial \varphi} & \frac{\partial}{\partial z} \\\\ 0 & 0 & \frac{\mu I}{2 \pi} \ln \left(\frac{1}{\rho}\right)\end{array}\right|\\ &= \frac{\partial}{\partial \rho}\left(\frac{\mu I}{2 \pi} \ln \left(\frac{1}{\rho}\right)\right)\hat{\mathbf{e}}_{\varphi}\\ &= \frac{\mu I}{2 \pi}\frac{\partial}{\partial \rho}\left( \ln \left(\frac{1}{\rho}\right)\right)\hat{\mathbf{e}}_{\varphi}\\ &= \frac{\mu I}{2 \pi}\frac{1}{\rho}\hat{\mathbf{e}}_{\varphi} \end{aligned}

Ejercicio 3.10.21

Dado que la transformacion de coordenadas cartesianas a cilindricas es:

(Cos(φc)Sin(φc)0Sin(φc)Cos(φc)0001)\begin{pmatrix} Cos(\varphi_c) & Sin(\varphi_c) & 0\\ -Sin(\varphi_c) & Cos(\varphi_c) & 0\\ 0 & 0 & 1 \end{pmatrix}

y de coordenadas esfericas a cartes:

(sin(θe)cos(φe)cos(θe)cos(φe)sin(φe)sin(θe)sin(φe)cos(θe)sin(φe)sin(φe)cos(θe)sin(θe)0)\begin{pmatrix} \sin(\theta_e)\cos(\varphi_e) & \cos(\theta_e)\cos(\varphi_e) & -\sin(\varphi_e)\\ \sin(\theta_e)\sin(\varphi_e) & \cos(\theta_e)\sin(\varphi_e) & \sin(\varphi_e)\\ \cos(\theta_e) & -\sin(\theta_e) & 0\\ \end{pmatrix}

Luego :

(r^eθ^eϕ^e)=(cos(φc)sin(φc)0sin(φc)cos(φc)0001)(sin(θe)cos(φe)cos(θe)cos(φe)sin(φe)sin(θe)sin(φe)cos(θe)sin(φe)sin(φe)cos(θe)sin(θe)0)(ρ^cθ^ck^c)=(sin(θe)cos(φe)cos(φc)+sin(θe)sin(φe)sin(φc)cos(θe)cos(φe)cos(φc)+cos(θe)sin(φe)sin(φc)sin(φe)cos(φc)+cos(φe)sin(φc)sin(θe)cos(φe)sin(φc)+sin(θe)sin(φe)cos(φc)cos(θe)cos(φe)sin(φc)+cos(θe)sin(φe)cos(φc)sin(θe)sin(φe)sin(φc)+cos(φe)cos(φc)cos(θe)sin(θe)0)(ρ^cθ^ck^c)=(sin(θe)cos(φ)cos(φ)+sin(θe)sin(φ)sin(φ)cos(θe)cos(φ)cos(φ)+cos(θe)sin(φ)sin(φ)sin(φ)cos(φ)+cos(φ)sin(φ)sin(θe)cos(φ)sin(φ)+sin(θe)sin(φ)cos(φ)cos(θe)cos(φ)sin(φ)sin(φ)+cos(φ)cos(φ)cos(θe)sin(θe)0)(ρ^cθ^ck^c)\begin{equation*} \begin{split} &\begin{pmatrix} \hat{\bold r}_e\\ \hat{\theta}_e \\ \hat{\phi}_e\\ \end{pmatrix} = \begin{pmatrix} \cos(\varphi_c) & \sin(\varphi_c) & 0\\ -\sin(\varphi_c) & \cos(\varphi_c) & 0\\ 0 & 0 & 1 \end{pmatrix} \begin{pmatrix} \sin(\theta_e)\cos(\varphi_e) & \cos(\theta_e)\cos(\varphi_e) & -\sin(\varphi_e)\\ \sin(\theta_e)\sin(\varphi_e) & \cos(\theta_e)\sin(\varphi_e) & \sin(\varphi_e)\\ \cos(\theta_e) & -\sin(\theta_e) & 0\\ \end{pmatrix} \begin{pmatrix} \hat{\rho}_c\\ \hat{\theta}_c\\ \hat{k}_c \end{pmatrix} \\\\ &=\begin{pmatrix} \sin(\theta_e)\cos(\varphi_e)\cos(\varphi_c) +\sin(\theta_e)\sin(\varphi_e)\sin(\varphi_c) & \cos(\theta_e)\cos(\varphi_e)\cos(\varphi_c) + \cos(\theta_e)\sin(\varphi_e)\sin(\varphi_c) & -\sin(\varphi_e)\cos(\varphi_c) + \cos(\varphi_e)\sin(\varphi_c) \\ -\sin(\theta_e)\cos(\varphi_e)\sin(\varphi_c) + \sin(\theta_e)\sin(\varphi_e)\cos(\varphi_c) & -\cos(\theta_e)\cos(\varphi_e)\sin(\varphi_c)+\cos(\theta_e)\sin(\varphi_e) \cos(\varphi_c) & \sin(\theta_e)\sin(\varphi_e)\sin(\varphi_c) + \cos(\varphi_e) \cos(\varphi_c) \\ \cos(\theta_e) & \sin(\theta_e) & 0 \end{pmatrix} \begin{pmatrix} \hat{\rho}_c\\ \hat{\theta}_c\\ \hat{k}c \end{pmatrix} \\\\ &=\begin{pmatrix} \sin(\theta_e)\cos(\varphi)\cos(\varphi) +\sin(\theta_e)\sin(\varphi)\sin(\varphi) & \cos(\theta_e)\cos(\varphi)\cos(\varphi) + \cos(\theta_e)\sin(\varphi)\sin(\varphi) & -\sin(\varphi)\cos(\varphi) + \cos(\varphi)\sin(\varphi) \\ -\sin(\theta_e)\cos(\varphi)\sin(\varphi) + \sin(\theta_e)\sin(\varphi)\cos(\varphi) & -\cos(\theta_e)\cos(\varphi)\sin(\varphi)\sin(\varphi) + \cos(\varphi) \cos(\varphi) \\ \cos(\theta_e) & \sin(\theta_e) & 0 \end{pmatrix} \begin{pmatrix} \hat{\rho}_c\\ \hat{\theta}_c\\ \hat{k}_c \end{pmatrix} \end{split} \end{equation*}

Donde r^e,θ^e,ϕ^e\hat{\mathbf{r}}_e, \hat{\theta}_e ,\hat{\phi}_e son los vectores unitarios asociados a coordenadas esferica y ρ^c,θ^c,k^c\hat{\rho}_c, \hat{\theta}_c, \hat{k}_c son los vectores unitarios las coordenadas ciclindricas. Ademas como el productos de matrices ortogonales es una matriz ortogonal entonces la matriz inversa sera la traspuesta de esta matriz.

(e^ρe^θe^φ)=(sin(θe)cos(φe)sin(θe)sin(φe)cos(θe)cos(φe)cos(φe)sin(φe)cos(θe)sin(θe)sin(φe)cos(φe)0)(i^j^k^)\begin{pmatrix} \hat{e}_{\rho} \\ \hat{e}_{\theta} \\ \hat{e}_\varphi \end{pmatrix} = \begin{pmatrix} \sin(\theta_e)\cos(\varphi_e) & \sin(\theta_e)\sin(\varphi_e) & -\cos(\theta_e)\\ \cos(\varphi_e)\cos(\varphi_e) & \sin(\varphi_e)\cos(\theta_e) & -\sin(\theta_e)\\ -\sin(\varphi_e) & \cos(\varphi_e) & 0\\ \end{pmatrix} \begin{pmatrix} \hat{i} \\ \hat{j} \\ \hat{k} \end{pmatrix}

En el sistema de coordenadas cilindricas tenemos:

(e^ρe^θe^z)=(cos(θe)sin(θe)0sin(θe)cos(θe)0001)(i^j^k^)\begin{pmatrix} \hat{e}_{\rho} \\ \hat{e}_{\theta} \\ \hat{e}_z \end{pmatrix} = \begin{pmatrix} \cos(\theta_e) & \sin(\theta_e) & 0\\ -\sin(\theta_e) & \cos(\theta_e) & 0\\ 0 & 0 & 1\\ \end{pmatrix} \begin{pmatrix} \hat{i} \\ \hat{j} \\ \hat{k} \end{pmatrix}

La matriz que convierte de esfericas a cilindricas es:

C=ATBC=(sin(θe)cos(φe)cos(φe)cos(θe)sin(φe)sin(θe)sin(φe)sin(φe)cos(θe)cos(φe)cos(θe)sin(θe)0)(cos(θe)sin(θe)0sin(θe)cos(θe)0001)=(sin(θe)cos(θe)0001cos(θe)sin(θe)0)C = A^{T} B \\ C= \begin{pmatrix} \sin(\theta_e)\cos(\varphi_e) & \cos(\varphi_e)\cos(\theta_e) & -\sin(\varphi_e)\\ \sin(\theta_e)\sin(\varphi_e) & \sin(\varphi_e)\cos(\theta_e) & \cos(\varphi_e)\\ \cos(\theta_e) & -\sin(\theta_e) & 0\\ \end{pmatrix} \begin{pmatrix} \cos(\theta_e) & \sin (\theta_e) & 0\\ -\sin (\theta_e) & \cos (\theta_e) & 0\\ 0& 0& 1 \end{pmatrix} \\= \begin{pmatrix} \sin (\theta_e) & \cos (\theta_e) & 0\\ 0 & 0 &1\\ \cos(\theta_e)&- \sin( \theta_e) & 0 \end{pmatrix}

Ejercicio 3.10.35

(a)De la ecuacion (2) tenemos que:

×B=1hrhφhθe^rhre^φhφe^θhθrφθhrFrhφFφhθFθ=1r2sinθe^re^φrsin(θ)e^θrrφθ2Pcosθr30rPr3sinθ=e^φrsin(θ)[r(Pr2sinθ)θ(2Pcosθr3)]=e^φrsin(θ)[2Pr3sinθ(2Psinθr3)]=e^φrsin(θ)[2Pr3sinθ+2Psinθr3]=0\begin{aligned}\boldsymbol{\nabla} \times \mathbf{B} &=\frac{1}{h_r h_\varphi h_\theta}\left|\begin{array}{ccc}\hat{\mathbf{e}}_r h_r & \hat{\mathbf{e}}_\varphi h_\varphi& \hat{\mathbf{e}}_\theta h_\theta \\\\ \frac{\partial}{\partial r} & \frac{\partial}{\partial \varphi} & \frac{\partial}{\partial \theta} \\\\ h_r F_r & h_\varphi F_\varphi & h_\theta F_\theta \end{array}\right|\\\\ &=\frac{1}{r^2 \sin {\theta}}\left|\begin{array}{ccc}\hat{\mathbf{e}}_r & \hat{\mathbf{e}}_\varphi r\sin{(\theta)} & \hat{\mathbf{e}}_\theta r \\\\ \frac{\partial}{\partial r} & \frac{\partial}{\partial \varphi} & \frac{\partial}{\partial \theta} \\\\ \frac{2 P \cos \theta}{r^3} & 0 & r \frac{P}{\cancel r^3} \sin \theta \end{array}\right|\\ &= -\hat{\mathbf{e}}_\varphi r\sin{(\theta)} \left[\frac{\partial}{\partial r}\left(\frac{P}{r^2} \sin \theta\right) - \frac{\partial}{\partial \theta}\left(\frac{2 P \cos \theta}{r^3}\right)\right]\\ &= -\hat{\mathbf{e}}_\varphi r\sin{(\theta)} \left[-\frac{2P}{r^3} \sin \theta - \left(-\frac{2 P \sin \theta}{r^3}\right)\right]\\ &= -\hat{\mathbf{e}}_\varphi r\sin{(\theta)} \left[-\frac{2P}{r^3} \sin \theta +\frac{2 P \sin \theta}{r^3}\right] = 0 \end{aligned}

Por lo cual existe una función escalar ϕ \phi tal que F=ϕ\mathbf F =\mathbf \nabla \phi

(b) Dado que que el trabajo se calcula para para un circulo de r=1r = 1 en plano θ=π2\theta = \frac{\pi}{2}, teniendo la siguiente parametrizacion c:[0,2π]R3\mathbf c:[0,2\pi] \rightarrow \Reals^3 , c(ϕ)=r(1,π/2,ϕ)\mathbf c(\phi) = \mathbf r(1,\pi/2, \phi)

Fdr=02πF(c(ϕ))dc(ϕ)dϕdϕ=02πF(1,π2,φ)dr(1,π/2,ϕ)dϕdϕ=02π(e^r2Pcosπ213+e^θP13sinπ2)ddϕ(e^xsinπ2cosφ+e^ysinπ2sinφ+e^zcosπ2)dϕ=02π(Pe^θ)(e^xsinφ+e^ycosφ)dϕ=02πPe^θe^φdϕ=02π0dϕ=0\begin{aligned} \oint \mathbf{F} \cdot d \mathbf{r} &= \int_0^{2\pi} \mathbf{F}(\mathbf c(\phi)) \cdot \frac{d \mathbf c(\phi)}{d\phi} d \phi\\ &= \int_0^{2\pi} \mathbf{F}(1,\frac{\pi}{2},\varphi) \cdot \frac{d \mathbf r(1,\pi/2, \phi)}{d \phi} d \phi\\ &= \int_0^{2\pi} \left(\hat{\mathbf{e}}_r \frac{2 P \cos \frac{\pi}{2}}{1^3}+\hat{\mathbf{e}}_\theta \frac{P}{1^3} \sin \frac{\pi}{2}\right) \cdot \frac{d}{d \phi}\left(\hat{\mathbf{e}}_x \sin {\frac{\pi}{2}} \cos \varphi+\hat{\mathbf{e}}_y \sin {\frac{\pi}{2}}\sin \varphi+\hat{\mathbf{e}}_z \cos {\frac{\pi}{2}}\right) d \phi\\ &= \int_0^{2\pi} \left( P \hat{\mathbf{e}}_\theta\right) \cdot \left(-\hat{\mathbf{e}}_x\sin \varphi+\hat{\mathbf{e}}_y\cos \varphi\right) d \phi\\ &= \int_0^{2\pi} P\hat{ \mathbf{e}}_\theta \cdot \hat{\mathbf{e}}_\varphi d \phi\\ &= \int_0^{2\pi} 0 d \phi =0 \end{aligned}

De los anterior la fuerza podria ser conservativa.

(c) Como el rotacional de la fuerza es cero entonces la fuerza es una fuerza conservativa, lo cual implica que el trabajo realizado por la fuerza sera la misma por cualquier trayectoria entre dos puntos, en particular en una trayectoria en linea recta a lo largo de la direccion radial e^r\mathbf{\hat e}_r y tomando el potencial como cero en el infinito entonces:

φr=r0Fdr=r0F(r(r))e^rdr=r0(e^r2Pcosθr3+e^θPr3sinθ)e^rdr=r02Pcosθr3dr=(2Pcosθr2)r0=(2Pcosθr02)\begin{aligned} \varphi_r &= -\int_{\infty}^{r_0} \mathbf{F}\cdot \mathbf {dr} \\ &= -\int_{\infty}^{r_0} \mathbf{F}(\mathbf r(r))\cdot \mathbf {\hat e}_r dr \\ &= -\int_{\infty}^{r_0}\left(\hat{\mathbf{e}}_r \frac{2 P \cos \theta}{r^3}+\hat{\mathbf{e}}_\theta \frac{P}{r^3} \sin \theta\right) \cdot \mathbf {\hat e}_r dr \\ &= -\int_{\infty}^{r_0} \frac{2 P \cos \theta}{r^3} dr \\ &= -\left(-\frac{2 P \cos \theta}{r^2}\right)\bigg|_{\infty}^{r_0} \\ &= \left(\frac{2 P \cos \theta}{r_0^2}\right) \\ \end{aligned}

Ejercicio 3.10.37

Sea θ\theta el ángulo tal que pr=prcos(θ) \mathbf{p} \cdot \mathbf{r} = pr\cos{(\theta)} luego:

E=(pr4πε0r3)=(pcos(θ)4πε0r2)=(e^rψr+e^θ1rψθ+e^φ1rsinθψφ)(pcos(θ)4πε0r2)=e^rr(pcos(θ)4πε0r2)+e^θ1rθ(pcos(θ)4πε0r2)+e^φ1rsinθφ(pcos(θ)4πε0r2)=e^r(pcos(θ)2πε0r3)e^θ(psin(θ)4πε0r3)\begin{aligned} \mathbf E &=\mathbf \nabla \left(\frac{\mathbf{p} \cdot \mathbf{r}}{4 \pi \varepsilon_0 r^3}\right)\\ &=\mathbf \nabla \left(\frac{p\cos{(\theta)}}{4 \pi \varepsilon_0 r^2}\right)\\ &= \left(\hat{\mathbf{e}}_r \frac{\partial \psi}{\partial r}+\hat{\mathbf{e}}_\theta \frac{1}{r} \frac{\partial \psi}{\partial \theta}+\hat{\mathbf{e}}_{\varphi} \frac{1}{r \sin \theta} \frac{\partial \psi}{\partial \varphi}\right)\left(\frac{p\cos{(\theta)}}{4 \pi \varepsilon_0 r^2}\right)\\ &= \hat{\mathbf{e}}_r \frac{\partial}{\partial r}\left(\frac{p\cos{(\theta)}}{4 \pi \varepsilon_0 r^2}\right)+\hat{\mathbf{e}}_\theta \frac{1}{r} \frac{\partial}{\partial \theta}\left(\frac{p\cos{(\theta)}}{4 \pi \varepsilon_0 r^2}\right)+\hat{\mathbf{e}}_{\varphi} \frac{1}{r \sin \theta} \frac{\partial }{\partial \varphi}\left(\frac{p\cos{(\theta)}}{4 \pi \varepsilon_0 r^2}\right)\\ &=-\hat{\mathbf{e}}_r \left(\frac{p\cos{(\theta)}}{2 \pi \varepsilon_0 r^3}\right)-\hat{\mathbf{e}}_\theta \left(\frac{p\sin{(\theta)}}{4 \pi \varepsilon_0 r^3}\right)\\ \end{aligned}

Machado Volumen 1

Ejercicio 6.2

Aplicando separación de variables a la ecuación 2V=2Vx2+2Vy2+2Vz2=0\nabla^2 \mathbb{V}=\frac{\partial^2 \mathbb{V}}{\partial x^2}+\frac{\partial^2 \mathbb{V}}{\partial y^2}+\frac{\partial^2 \mathbb{V}}{\partial z^2}=0 se llega a la solución general:

V(x,y,z)=V0,0+Vk,0+V0,+Vk,k+Vk,\begin{equation}\mathbb{V}(x, y, z)=\mathbb{V}_{0,0}+\mathbb{V}_{k, 0}+\mathbb{V}_{0, \ell}+\mathbb{V}_{k, k}+\mathbb{V}_{k, \ell}\end{equation}

Donde:

V0,0=[a0,0+b0,0x][c0,0+d0,0y][f0,0+g0,0z]Vk,0=[ak,0+bk,0x][ck,0cosky+dk,0senky][fk,0ekz+gk,0ekz]V0,=[a0,cosx+b0,senx][c0,ey+d0,ey][f0,+g0,z]Vk,k=[ak,kcoskx+bk,ksenkx][ck,k+dk,ky][fk,kekz+gk,kekz]Vk,=[ak,cosx+bk,senx][ck,cosλy+dk,senλy][fk,ekz+gk,ekz]\begin{aligned} \mathbb{V}_{0,0}&=\left[a_{0,0}+b_{0,0} x\right]\left[c_{0,0}+d_{0,0} y\right]\left[f_{0,0}+g_{0,0} z\right]\\ \mathbb{V}_{k, 0}&=\left[a_{k, 0}+b_{k, 0} x\right]\left[c_{k, 0} \cos k y+d_{k, 0} \operatorname{sen} k y\right]\left[f_{k, 0} e^{k z}+g_{k, 0} e^{-k z}\right]\\ \mathbb{V}_{0, \ell}&=\left[a_{0, \ell} \cos \ell x+b_{0, \ell} \operatorname{sen} \ell x\right]\left[c_{0, \ell} e^{\ell y}+d_{0, \ell} e^{-\ell y}\right]\left[f_{0, \ell}+g_{0, \ell} z\right]\\ \mathbb{V}_{k, k}&=\left[a_{k, k} \cos k x+b_{k, k} \operatorname{sen} k x\right]\left[c_{k, k}+d_{k, k} y\right]\left[f_{k, k} e^{k z}+g_{k, k} e^{-k z}\right]\\ \mathbb{V}_{k, \ell}&=\left[a_{k, \ell} \cos \ell x+b_{k, \ell} \operatorname{sen} \ell x\right]\left[c_{k, \ell} \cos \lambda y+d_{k, \ell} \operatorname{sen} \lambda y\right]\left[f_{k, \ell} e^{k z}+g_{k, \ell} e^{-k z}\right] \end{aligned}

Y aplicando las condiciones de frontera entonces:

V(0,y,z)=V0,0(0,y,z)+Vk,0(0,y,z)+V0,(0,y,z)+Vk,k(0,y,z)+Vk,(0,y,z)=0 \begin{aligned} \mathbb{V}(0, y, z)=\mathbb{V}_{0,0}(0, y, z)+\mathbb{V}_{k, 0}(0, y, z)&+ \mathbb{V}_{0, \ell}(0, y, z) \\ &+\mathbb{V}_{k, k}(0, y, z)+\mathbb{V}_{k, \ell}(0, y, z)=0 \end{aligned} \hspace{0.5cm} \Rightarrow
V(0,y,z)=[a0,0+b0,0.0][c0,0+d0,0y][f0,0+g0,0z]+[ak,0+bk,00][ck,0cosky+dk,0senky][fk,0ekz+gk,0ekz]+[a0,cos0+b0,sen0][c0,ey+d0,ey][f0,+g0,z]+[ak,kcos0+bk,ksen0][ck,k+dk,ky][fk,kekz+gk,kekz]+[ak,cos0+bk,sen0][ck,cosλy+dk,senλy][fk,ekz+gk,ekz]=0\begin{aligned}&\mathbb{V}(0, y, z)=\left[a_{0,0}+b_{0,0} .0\right]\left[c_{0,0}+d_{0,0} y\right]\left[f_{0,0}+g_{0,0} z\right] \\&+\left[a_{k, 0}+b_{k, 0} \cdot 0\right]\left[c_{k, 0} \cos k y+d_{k, 0} \operatorname{sen} k y\right]\left[f_{k, 0} e^{k z}+g_{k, 0} e^{-k z}\right] \\&+\left[a_{0, \ell} \cos 0+b_{0, \ell} \operatorname{sen} 0\right]\left[c_{0, \ell} e^{\ell y}+d_{0, \ell} e^{-\ell y}\right]\left[f_{0, \ell}+g_{0, \ell} z\right] \\&\quad+\left[a_{k, k} \cos 0+b_{k, k} \operatorname{sen} 0\right]\left[c_{k, k}+d_{k, k} y\right]\left[f_{k, k} e^{k z}+g_{k, k} e^{-k z}\right] \\&+\left[a_{k, \ell} \cos 0+b_{k, \ell} \operatorname{sen} 0\right]\left[c_{k, \ell} \cos \lambda y+d_{k, \ell} \operatorname{sen} \lambda y\right]\left[f_{k, \ell} e^{k z}+g_{k, \ell} e^{-k z}\right]=0\end{aligned} \hspace{0.5cm}\Rightarrow
V(0,y,z)=a0,0[c0,0+d0,0y][f0,0+g0,0z]+ak,0[ck,0cosky+dk,0senky][fk,0ekz+gk,0ekz]+a0,[c0,ey+d0,ey][f0,+g0,z]+ak,k[ck,k+dk,ky][fk,kekz+gk,kekz]+ak,[ck,cosλy+dk,senλy][fk,ekz+gk,ekz]=0\begin{aligned}&\mathbb{V}(0, y, z)=a_{0,0}\left[c_{0,0}+d_{0,0} y\right]\left[f_{0,0}+g_{0,0} z\right] \\&+a_{k, 0}\left[c_{k, 0} \cos k y+d_{k, 0} \operatorname{sen} k y\right]\left[f_{k, 0} e^{k z}+g_{k, 0} e^{-k z}\right] \\&+a_{0, \ell}\left[c_{0, \ell} e^{\ell y}+d_{0, \ell} e^{-\ell y}\right]\left[f_{0, \ell}+g_{0, \ell} z\right] \\ &+a_{k, k} \left[c_{k, k}+d_{k, k} y\right]\left[f_{k, k} e^{k z}+g_{k, k} e^{-k z}\right] \\&+a_{k, \ell}\left[c_{k, \ell} \cos \lambda y+d_{k, \ell} \operatorname{sen} \lambda y\right]\left[f_{k, \ell} e^{k z}+g_{k, \ell} e^{-k z}\right]=0\end{aligned}

Lo cual implica que a,k=0a_{\ell,k} = 0  para todo \ell, kk, luego:

V0,0=b0,0x[c0,0+d0,0y][f0,0+g0,0z]Vk,0=bk,0x[ck,0cosky+dk,0senky][fk,0ekz+gk,0ekz]V0,=b0,senx[c0,ey+d0,ey][f0,+g0,z]Vk,k=bk,ksenkx[ck,k+dk,ky][fk,kekz+gk,kekz]Vk,=bk,senx[ck,cosλy+dk,senλy][fk,ekz+gk,ekz]\begin{aligned} \mathbb{V}_{0,0}&=b_{0,0} x\left[c_{0,0}+d_{0,0} y\right]\left[f_{0,0}+g_{0,0} z\right]\\ \mathbb{V}_{k, 0}&=b_{k, 0} x\left[c_{k, 0} \cos k y+d_{k, 0} \operatorname{sen} k y\right]\left[f_{k, 0} e^{k z}+g_{k, 0} e^{-k z}\right]\\ \mathbb{V}_{0, \ell}&=b_{0, \ell} \operatorname{sen} \ell x\left[c_{0, \ell} e^{\ell y}+d_{0, \ell} e^{-\ell y}\right]\left[f_{0, \ell}+g_{0, \ell} z\right]\\ \mathbb{V}_{k, k}&=b_{k, k} \operatorname{sen} k x\left[c_{k, k}+d_{k, k} y\right]\left[f_{k, k} e^{k z}+g_{k, k} e^{-k z}\right]\\ \mathbb{V}_{k, \ell}&=b_{k, \ell} \operatorname{sen} \ell x\left[c_{k, \ell} \cos \lambda y+d_{k, \ell} \operatorname{sen} \lambda y\right]\left[f_{k, \ell} e^{k z}+g_{k, \ell} e^{-k z}\right] \end{aligned}

Además:

V(x,0,z)=b0,0x[c0,0+d0,00][f0,0+g0,0z]+bk,0x[ck,0cos0+dk,0sen0][fk,0ekz+gk,0ekz]+b0,senx[c0,e0+d0,e0][f0,+g0,z]+bk,ksenkx[ck,k+dk,k.0][fk,kekz+gk,kekz]+bk,senx[ck,cos0+dk,sen0][fk,ekz+gk,ekz]=0\begin{aligned}&\mathbb{V}(x, 0, z)=b_{0,0} x\left[c_{0,0}+d_{0,0} \cdot 0\right]\left[f_{0,0}+g_{0,0} z\right] \\&+b_{k, 0} x\left[c_{k, 0} \cos 0+d_{k, 0} \operatorname{sen} 0\right]\left[f_{k, 0} e^{k z}+g_{k, 0} e^{-k z}\right] \\&+b_{0, \ell} \operatorname{sen} \ell x\left[c_{0, \ell} e^0+d_{0, \ell} e^0\right]\left[f_{0, \ell}+g_{0, \ell} z\right] \\&+b_{k, k} \operatorname{sen} k x\left[c_{k, k}+d_{k, k} .0\right]\left[f_{k, k} e^{k z}+g_{k, k} e^{-k z}\right] \\&+b_{k, \ell} \operatorname{sen} \ell x\left[c_{k, \ell} \cos 0+d_{k, \ell} \operatorname{sen} 0\right]\left[f_{k, \ell} e^{k z}+g_{k, \ell} e^{-k z}\right]=0\end{aligned}\hspace{0.5cm} \Rightarrow
V(x,0,z)=b0,0x[c0,0][f0,0+g0,0z]+bk,0xck,0[fk,0ekz+gk,0ekz]+b0,senx[c0,+d0,][f0,+g0,z]+bk,ksenkx[ck,k][fk,kekz+gk,kekz]+bk,senx[ck,][fk,ekz+gk,ekz]=0\begin{aligned}&\mathbb{V}(x, 0, z)=b_{0,0} x\left[c_{0,0}\right]\left[f_{0,0}+g_{0,0} z\right] \\ &+b_{k, 0} xc_{k, 0}\left[f_{k, 0} e^{k z}+g_{k, 0} e^{-k z}\right] \\ &+b_{0, \ell} \operatorname{sen} \ell x\left[c_{0, \ell} +d_{0, \ell} \right]\left[f_{0, \ell}+g_{0, \ell} z\right] \\ &+b_{k, k} \operatorname{sen} k x\left[c_{k, k}\right]\left[f_{k, k} e^{k z}+g_{k, k} e^{-k z}\right] \\ &+b_{k, \ell} \operatorname{sen} \ell x\left[c_{k, \ell}\right]\left[f_{k, \ell} e^{k z}+g_{k, \ell} e^{-k z}\right]=0\end{aligned}\hspace{0.5cm} \Rightarrow

lo cual implica que c0,0=ck,0=ck,k=ck,=0c_{0,0}=c_{k, 0}=c_{k, k}=c_{k, \ell}=0  y c0,=d0,c_{0, \ell}=-d_{0, \ell}, luego

V0,0=b0,0xd0,0y[f0,0+g0,0z]Vk,0=bk,0xdk,0senky[fk,0ekz+gk,0ekz]V0,=b0,senxc0,2senhy[f0,+g0,z]Vk,k=bk,ksenkxdk,ky[fk,kekz+gk,kekz]Vk,=bk,senxdk,senλy[fk,ekz+gk,ekz]\begin{aligned} \mathbb{V}_{0,0}&=b_{0,0} xd_{0,0} y\left[f_{0,0}+g_{0,0} z\right]\\ \mathbb{V}_{k, 0}&=b_{k, 0} xd_{k, 0} \operatorname{sen} k y\left[f_{k, 0} e^{k z}+g_{k, 0} e^{-k z}\right]\\ \mathbb{V}_{0, \ell}&=b_{0, \ell} \operatorname{sen} {\ell x} c_{0, \ell}2\operatorname{senh} \ell y \left[f_{0, \ell}+g_{0, \ell} z\right]\\ \mathbb{V}_{k, k}&=b_{k, k} \operatorname{sen} k x d_{k, k} y\left[f_{k, k} e^{k z}+g_{k, k} e^{-k z}\right]\\ \mathbb{V}_{k, \ell}&=b_{k, \ell} \operatorname{sen} \ell x d_{k, \ell} \operatorname{sen} \lambda y\left[f_{k, \ell} e^{k z}+g_{k, \ell} e^{-k z}\right] \end{aligned}

Y:

V(x,y,0)=b0,0xd0,0yf0,0+bk,0xdk,0senky[fk,0+gk,0]+b0,senxc0,2senhy[f0,]+bk,ksenkxdk,ky[fk,k+gk,k]+bk,senxdk,senλy[fk,+gk,]\begin{aligned} \mathbb{V}(x,y,0)&=b_{0,0} xd_{0,0} yf_{0,0}\\ &+b_{k, 0} xd_{k, 0} \operatorname{sen} k y\left[f_{k, 0}+g_{k, 0}\right]\\ &+b_{0, \ell} \operatorname{sen} {\ell x} c_{0, \ell}2\operatorname{senh} \ell y \left[f_{0, \ell}\right]\\ &+b_{k, k} \operatorname{sen} k x d_{k, k} y\left[f_{k, k} +g_{k, k}\right]\\ &+b_{k, \ell} \operatorname{sen} \ell x d_{k, \ell} \operatorname{sen} \lambda y\left[f_{k, \ell}+g_{k, \ell} \right] \end{aligned}
V0,0=b0,0xd0,0yg0,0zVk,0=bk,0xdk,0senkygk,0senkzV0,=b0,senxd0,2senhyg0,zVk,k=bk,ksenkxdk,kygk,ksenhkzVk,=bk,senxdk,senλygk,senhkz\begin{aligned} \mathbb{V}_{0,0}&=b_{0,0} xd_{0,0} yg_{0,0} z\\ \mathbb{V}_{k, 0}&=-b_{k, 0} xd_{k, 0} \operatorname{sen} k yg_{k, 0}\operatorname{sen} k z\\ \mathbb{V}_{0, \ell}&=-b_{0, \ell} \operatorname{sen} {\ell x} d_{0, \ell}2\operatorname{senh} \ell y g_{0, \ell} z\\ \mathbb{V}_{k, k}&=-b_{k, k} \operatorname{sen} k x d_{k, k} yg_{k, k}\operatorname{senh}{kz}\\ \mathbb{V}_{k, \ell}&=-b_{k, \ell} \operatorname{sen} \ell x d_{k, \ell} \operatorname{sen} \lambda y g_{k, \ell}\operatorname{senh}{kz} \end{aligned}

Ahora haciendo uso del principio de \

superposición, la solución general al problema es:

V=V1+V2+V3\mathbb V= \mathbb V_1 + \mathbb V_2 + \mathbb V_3

Donde V1\mathbb V_1 es la solución de

{V(0,y,z)=0V(a,y,z)=V0V(x,0,z)=0V(x,b,z)=0V(x,y,0)=0V(x,y,c)=0\begin{cases}\mathbb{V}(0, y, z)=0 & \mathbb{V}(a, y, z)=V_0 \\ \mathbb{V}(x, 0, z)=0 & \mathbb{V}(x, b, z)=0 \\ \mathbb{V}(x, y, 0)=0 & \mathbb{V}(x, y, c)=0\end{cases}

Donde V2\mathbb V_2 es la solución de

{V(0,y,z)=0V(a,y,z)=0V(x,0,z)=0V(x,b,z)=V0sinxV(x,y,0)=0V(x,y,c)=0\begin{cases}\mathbb{V}(0, y, z)=0 & \mathbb{V}(a, y, z)=0 \\ \mathbb{V}(x, 0, z)=0 & \mathbb{V}(x, b, z)=V_0 \sin x \\ \mathbb{V}(x, y, 0)=0 & \mathbb{V}(x, y, c)=0\end{cases}

Donde V3\mathbb V_3 es la solución de

{V(0,y,z)=0V(a,y,z)=0V(x,0,z)=0V(x,b,z)=0V(x,y,0)=0V(x,y,c)=V0cosx\begin{cases}\mathbb{V}(0, y, z)=0 & \mathbb{V}(a, y, z)=0 \\ \mathbb{V}(x, 0, z)=0 & \mathbb{V}(x, b, z)=0 \\ \mathbb{V}(x, y, 0)=0 & \mathbb{V}(x, y, c)=V_0 \cos x\end{cases}

Además como el potencial entre las tapas no metálicas debe ser continuo entonces tenemos la siguientes condiciones adicionales:

V1(x,b,z)=b0,0xbz+bk,0xsenkbsenkz+b0,senxsenhbz+bk,ksenkxbsenhkz+bk,senxsenλbsenhkz=0\begin{aligned} \mathbb{V}_1(x,b,z)&=b_{0,0} x b z\\ &+b_{k, 0} x \operatorname{sen} k b\operatorname{sen} k z\\ &+b_{0, \ell} \operatorname{sen} {\ell x} \operatorname{senh} \ell b z\\ &+b_{k, k} \operatorname{sen} kx b\operatorname{senh}{kz}\\ &+b_{k, \ell} \operatorname{sen} \ell x \operatorname{sen} \lambda b \operatorname{senh}{kz} = 0 \end{aligned}

Luego b0,0=b0,=bk,l=bk,k=0b_{0,0}=b_{0, \ell}=b_{k,l}=b_{k, k}=0 y kn=nπbk_n=\frac{n \pi}{b} por lo que

V1(x,y,z)=bk,0xsennπzbsennπzb\begin{aligned} \mathbb{V}_1(x,y,z) =b_{k, 0} x \operatorname{sen} \frac{n \pi z} b\operatorname{sen} \frac{n \pi z}{b} \end{aligned}

Y también:

V1(a,y,z)=bk,0asennπzbsennπzb=V0a=V0[bk,0sennπzbsennπzb]1\begin{aligned} \mathbb{V}_1(a,y,z) =b_{k, 0} a \operatorname{sen} \frac{n \pi z} b\operatorname{sen} \frac{n \pi z}{b} =V_0 \hspace{0.5cm}\Rightarrow \end{aligned}\\ \begin{aligned} a =V_0\left[b_{k, 0} \operatorname{sen} \frac{n \pi z} b\operatorname{sen} \frac{n \pi z}{b}\right]^{-1}\hspace{0.5cm}\Rightarrow \end{aligned}

Por lo que k=0k=0 y

V1(x,y,z)=V0\mathbb V_1(x,y,z) = V_0

Ahora para V2\mathbb V_2

V2(a,y,z)=b0,0ayz+bk,0asenkysenkz+b0,zsenasenhy+bk,kysenkasenhkz+bk,senasenλysenhkz=0\begin{aligned} \mathbb{V}_2(a,y,z)&=b_{0,0} a y z\\ &+b_{k, 0} a \operatorname{sen} k y\operatorname{sen} k z\\ &+b_{0, \ell} z\operatorname{sen} {\ell a} \operatorname{senh} \ell y \\ &+b_{k, k}y \operatorname{sen} ka \operatorname{senh}{kz}\\ &+b_{k, \ell} \operatorname{sen} \ell a \operatorname{sen} \lambda y \operatorname{senh}{kz} = 0 \end{aligned}

Luego b0,0=bk,0=bk,k=0b_{0,0}=b_{k,0}=b_{k, k}=0 y n=nπa\ell_n=\frac{n \pi}{a} por lo que

V2(x,y,c)=sennπxa[b0,seny+bk,senλysenhkc]=0b0,seny+bk,senλysenhkc=0\begin{aligned} \mathbb{V}_2(x,y,c) &= \operatorname{sen} \frac{n\pi x}{a}\left[b_{0, \ell}\operatorname{sen} \ell y+b_{k, \ell}\operatorname{sen} \lambda y \operatorname{senh}{kc}\right]=0\hspace{0.5cm} \Rightarrow \end{aligned}\\ \begin{aligned} b_{0, \ell}\operatorname{sen} \ell y+b_{k, \ell}\operatorname{sen} \lambda y \operatorname{senh}{kc}=0\hspace{0.5cm} \Rightarrow \end{aligned}\\
V2(x,y,c)=b0,0xyc+bk,0xsenkysenkc+b0,csenxsenhy+bk,kysenkxsenhkc+bk,senxsenλysenhkc=0\begin{aligned} \mathbb{V}_2(x,y,c)&=b_{0,0} x y c\\ &+b_{k, 0} x \operatorname{sen} k y\operatorname{sen} k c\\ &+b_{0, \ell} c \operatorname{sen} {\ell x} \operatorname{senh} \ell y \\ &+b_{k, k}y \operatorname{sen} kx \operatorname{senh}{kc}\\ &+b_{k, \ell} \operatorname{sen} \ell x \operatorname{sen} \lambda y \operatorname{senh}{kc} = 0 \end{aligned}

Luego b0,0=bk,0=bk,k=0b_{0,0}=b_{k,0}=b_{k, k}=0 y kn=nπak_n=\frac{n \pi}{a} por lo que

V2(a,y,z)=bk,0asenkysenkz=0b0,seny+bk,senλysenhkc=0\begin{aligned} \mathbb{V}_2(a,y,z) &= b_{k, 0} a \operatorname{sen} k y\operatorname{sen} k z=0\hspace{0.5cm} \Rightarrow \end{aligned}\\ \begin{aligned} b_{0, \ell}\operatorname{sen} \ell y+b_{k, \ell}\operatorname{sen} \lambda y \operatorname{senh}{kc}=0\hspace{0.5cm} \Rightarrow \end{aligned}\\

Y también:

Entonces:

V1(x,y,γ)=bk,senxsennπyβsenhkγ=0\begin{aligned} \mathbb{V}_1(x,y,\gamma) =b_{k, \ell} \operatorname{sen} \ell x \operatorname{sen} \frac{n \pi y}{\beta} \operatorname{senh}{k\gamma} = 0 \end{aligned}
V1(a,y,z)=b0,0ad0,0yg0,0zbk,0adk,0senkygk,0senkzb0,senad0,2senhyg0,zbk,ksenkadk,kygk,ksenhkzbk,senadk,senλygk,senkz=V0\begin{aligned} \mathbb{V}_1(a,y,z)&=b_{0,0} a d_{0,0} yg_{0,0} z\\ &-b_{k, 0} a d_{k, 0} \operatorname{sen} k yg_{k, 0}\operatorname{sen} k z\\ &-b_{0, \ell} \operatorname{sen} {\ell a} d_{0, \ell}2\operatorname{senh} \ell y g_{0, \ell} z\\ &-b_{k, k} \operatorname{sen} k a d_{k, k} yg_{k, k}\operatorname{senh}{kz}\\ &-b_{k, \ell} \operatorname{sen} \ell a d_{k, \ell} \operatorname{sen} \lambda y g_{k, \ell}\operatorname{sen}{kz} = V_0 \end{aligned}

Ejercicio 7.3 Machado

Tenemos que,

V=Q4πε0[1(xX)2+(yY)2+z21(x+X)2+(yY)2+z2+1(x+X)2+(y+Y)2+z21(xX)2+(y+Y)2+z2]V = \frac{Q}{4 \pi \varepsilon_0} \left[ \frac{1}{\sqrt{(x-X)^2+ (y-Y)^2 +z^2}} - \frac{1}{\sqrt{(x+X)^2+ (y-Y)^2 +z^2}}+ \frac{1}{\sqrt{(x+X)^2+ (y+Y)^2 +z^2}} - \frac{1}{\sqrt{(x-X)^2+ (y+Y)^2 +z^2}} \right]

E=φ=φxi^φyj^φzk^\vec{E} = - \nabla \varphi = - \frac{\partial \varphi}{\partial x}\hat{i} - \frac{\partial \varphi}{\partial y}\hat{j} - \frac{\partial \varphi}{\partial z}\hat{k}
φx=Q4πε0((xX)(xX)2+(yY)2+z2+(x+X)(x+X)2(yY2)+z2(x+X)(x+X)2+(y+Y)2+z2+(xX)(xX)2+(y+Y)2+z2)=Q4πε0(2X(x+X)2+(y+Y)2+z2)\frac{\partial \varphi}{\partial x} = \frac{Q}{4 \pi \varepsilon_0} \left( -\frac{(x-X)}{\sqrt{(x-X)^2 + (y-Y)^2 + z^2}} + \frac{(x+X)}{\sqrt{(x+X)^2(y-Y^2)+z^2}} - \frac{(x+X)}{\sqrt{(x+X)^2+ (y+Y)^2 + z^2}}+ \frac{(x-X)}{\sqrt{(x-X)^2 + (y+Y)^2 +z^2}}\right)\\ = \frac{Q}{4 \pi \varepsilon_0} \left( - \frac{2X}{\sqrt{(x+X)^2+ (y+Y)^2 + z^2}}\right)

Luego

E=Q4πϵ0[2Xi^+2Yj^+2zk^(x+X)2+(y+Y)2+z2]\vec{E} = \frac{Q}{4 \pi \epsilon_0} \left[ \frac{2X \hat{i} + 2 Y\hat{j} + 2z \hat{k}}{\sqrt{(x+X)^2 + (y+Y)^2 + z^2}} \right]

Para el caso de σ=ε0En^\sigma = \varepsilon_0 \vec{E} \cdot \hat{n} la densidad de carga se distribuye en i^\hat{i} y j^\hat{j} luego:

σ=Q4π[2Xi^+2Yj^(x+X)2+(y+Y)2+z2]\sigma = \frac{Q}{4 \pi} \left[ \frac{2X \hat{i} + 2Y \hat{j}}{ \sqrt{(x+X)^2 + (y+Y )^2 +z^2}} \right]
F=QσA4πϵ0rrrra3=Q8πε0sXi^+2Yj^(xX)2+(yY)2+z2\partial F = \frac{Q \sigma \partial A}{4 \pi \epsilon_0} \frac{\vec{r}- \vec{r'}}{|\vec{r} - \vec{r_a}|^3} \\ = \frac{Q}{8 \pi \varepsilon_0} \frac{sX \hat{i} + 2 Y \hat{j}}{\sqrt{(x-X)^2+ (y-Y)^2 +z^2}}

Fx=Q8πε002X(xX)2+(yY)2+z2xF_x = \frac{Q}{8 \pi \varepsilon_0} \int_{0}^{\infty} \frac{2X}{\sqrt{(x-X)^2 + (y-Y)^2 +z2}} \partial x

Que resulta ser la fuerza electrica que actua sobre la carga QQ

7.8 Machado

7.10 Machado