Problemas de contorno · Hoja de problemas

Hoja de problemas del módulo II.2

Seis problemas de dificultad creciente sobre problemas de contorno, con pistas y sin soluciones. La brújula del módulo: toda solución construida aquí vale en una región y en ninguna otra, y el error que más veces se comete es usarla fuera de ella.

Si de este módulo hay que llevarse un solo aviso, es una pregunta que hay que hacerse antes de escribir cada resultado: ¿en qué región vale esto que acabo de obtener? Las cuatro técnicas del módulo construyen una función que resuelve el problema en una región concreta y no dicen absolutamente nada fuera de ella, y todos los errores clásicos del tema consisten en olvidarlo. Debajo de un plano a tierra la fórmula del par carga-imagen predice 449 kV/m justo bajo la carga y el campo real es exactamente cero. Dentro de la esfera de Kelvin no hay solución ninguna: la imagen vive ahí precisamente porque ése no es el sitio donde se está resolviendo nada. Dentro de la esfera más pequeña que contiene la carga, el desarrollo multipolar no aproxima mal — diverge. Y una serie de Fourier converge al dato de la frontera en todos los puntos salvo justo en la discontinuidad, donde se pasa un 17,9 % de V0V_0 por muchos términos que se sumen. Antes de dar un número por bueno, dibuja la región y comprueba que el punto está dentro.

Necesitas: el módulo II.2 completo. Los problemas 1 y 6 se apoyan en los dos teoremas de unicidad del artículo 01; el 2, el 5 y el 6 en el artículo 02; el 3 en el artículo 03; el 4 y el 5 en el artículo 04. Del módulo II.1 hacen falta, sólo para el problema 2, σ=ε0E\sigma = \varepsilon_0E^{\perp} y la presión electrostática σ2/2ε0\sigma^{2}/2\varepsilon_0 del artículo 04.
Problema 1 ●○○ · El dato que no se puede pedir

Una superficie esférica imaginaria de radio R=10R = 10 cm delimita una región de vacío. Un programa de cálculo te deja plantear el problema de contorno dando sobre ella la derivada normal V/n\partial V/\partial n en lugar del potencial: en vez de una condición de contorno de Dirichlet, una condición de contorno de Neumann.

(a) Le das V/n=Acosθ\partial V/\partial n = -A\cos\theta con A=500A = 500 V/m. Monta y resuelve la integral de compatibilidad y di si el problema tiene solución con la región vacía; si la tiene, escribe la VV que la resuelve y di qué campo es. (b) Demuestra que V+CV + C es solución para cualquier constante y que el campo no cambia. ¿Por qué eso no es un defecto del planteamiento? (c) Ahora le das V/n=B\partial V/\partial n = -B constante, con el mismo valor numérico. Repite la integral y di qué le ocurre al problema. Si insistes en ese dato, ¿cuánta carga tiene que haber dentro, y de qué signo? (d) Comprueba tu respuesta a (c) exhibiendo una configuración concreta que produzca exactamente ese dato sobre la esfera, y explica por qué con el potencial prescrito en la frontera esta situación no puede darse nunca, sea cual sea el dato.

Pista

Todo se decide antes de resolver nada: integra V/n\partial V/\partial n sobre la esfera entera y compárala con Qenc/ε0-Q_{\text{enc}}/\varepsilon_0. La trampa es el signo, porque V/n\partial V/\partial n es En-E_n y no EnE_n: con el signo cambiado, la carga que el dato exige te saldrá del revés. Para comprobarte, deriva la VV que propongas y mira si devuelve el dato en r=Rr=R. Anclas: la que resuelve (a) es V=ArcosθV = -A\,r\cos\theta, o sea un campo uniforme de 500 V/m; y el dato constante de (c) exige Qenc=ε0B4πR2=Q_{\text{enc}} = \varepsilon_0B\,4\pi R^{2} = +0,556 nC.

Problema 2 ●●○ · La chapa tira de la carga, y casi toda desde un botón

El artículo 02 calculó la fuerza sobre la carga usando la imagen. Aquí se calcula la fuerza sobre la chapa, que es la misma cambiada de signo, sin usar la imagen para nada más que para conocer σ(r)\sigma(r).

Una carga puntual q=4q = 4 nC está a d=1d = 1 cm de un plano conductor infinito conectado a tierra.

(a) Evalúa σ(0)\sigma(0) y la presión electrostática justo debajo de la carga, en pascales. (b) Monta y resuelve la integral de σ2/2ε0\sigma^{2}/2\varepsilon_0 sobre todo el plano y comprueba que devuelve la fuerza de la imagen; da su valor y la masa que pesaría lo mismo. (c) ¿Dentro de qué radio se genera la mitad de esa fuerza? ¿Y el 90 %? (d) Compara esos dos radios con los que contienen la mitad y el 90 % de la carga inducida, y explica la diferencia con una sola frase sobre los integrandos.

Pista

La σ(r)\sigma(r) ya está deducida en el artículo 02: métela en la presión y monta la integral con da=2πrdrda = 2\pi r\,dr. La trampa es qué campo empuja a cada elemento de la chapa —σ/2ε0\sigma/2\varepsilon_0 y no σ/ε0\sigma/\varepsilon_0, porque una carga no se empuja a sí misma—, y con el campo entero sale exactamente el doble de fuerza. Tu integral tiene que devolver la que da la carga imagen; si no coincide, el 2 perdido está ahí. Anclas: σ(0)=\sigma(0) = −6,37 μC/m², P(0)=P(0) = 2,29 Pa, y la integral tiene que devolver 0,359 mN, el peso de 36,7 mg.

Problema 3 ●●○ · La boca partida en dos

La ranura del artículo 03 —dos placas conductoras a tierra en y=0y=0 e y=ay=a, con a=10a = 10 cm, extendiéndose hacia x>0x>0— pero con la tira de la boca partida por la mitad y las dos partes aisladas entre sí: la de abajo (0<y<a/20<y<a/2) a V0-V_0 y la de arriba (a/2<y<aa/2<y<a) a +V0+V_0, con V0=100V_0 = 100 V.

(a) Monta y resuelve la integral de Fourier que da los CnC_n. ¿Cuáles sobreviven, y cuáles se anulan? (b) Demuestra que el plano y=a/2y=a/2 es una superficie nodal: V=0V=0 en él para toda xx. Deduce de ahí, invocando el teorema de unicidad, que se puede meter una lámina metálica a tierra por y=a/2y=a/2 sin cambiar absolutamente nada, y di qué problema queda entonces a cada lado. (c) Usa lo anterior para construir la forma cerrada a partir de la del artículo, y evalúa VV en (a/4,3a/4)(a/4,\,3a/4) con un solo término y con la forma cerrada; da el error del truncamiento en porcentaje de V0V_0. (d) ¿A qué distancia de la boca la amplitud del potencial baja de 1 V? Compárala con la del escalón del artículo y explica el factor exacto que sale. (e) ¿Sigue habiendo fenómeno de Gibbs, y en cuántos sitios de la boca?

Pista

Rehaz la integral de Fourier con la boca partida en vez de reutilizar los coeficientes del escalón. La trampa: el desarrollo en senos prolonga la función de forma impar respecto a y=0y=0, pero lo que decide qué armónicos aparecen es la simetría de la boca respecto a su propio centro. Suma tu serie en x=0x=0 y mira si reproduce el escalón con el salto en el sitio correcto. Anclas: sobreviven sólo los n2 (mod 4)n \equiv 2 \ (\mathrm{mod}\ 4), y el primero vale C2=8V0/2π=|C_2| = 8V_0/2\pi = 127,3 V, el mismo número que el C1C_1 del escalón.

Problema 4 ●●○ · El anillo que de lejos es una carga puntual, salvo mirándolo de frente

Un anillo circular de radio R=5R = 5 cm lleva Q=10Q = 10 nC repartidos uniformemente. Está en el plano z=0z=0 y centrado en el origen.

(a) Monta el momento multipolar de orden ll y demuestra que vale QRlPl(0)QR^{l}P_l(0). Evalúa los de l=0l = 0 a 44 y di cuáles se anulan y qué simetría lo obliga. (b) Escribe V(r,θ)V(r,\theta) fuera y evalúalo en el eje a 20 cm con uno, dos y tres términos; compáralo con el valor exacto kQ/r2+R2kQ/\sqrt{r^{2}+R^{2}} y da los tres errores. (c) ¿A partir de qué distancia basta el monopolo con un error menor del 1 %, en el eje y en el ecuador? ¿Por qué no es el mismo sitio, si el objeto es el mismo? (d) Encuentra la dirección en la que el anillo es indistinguible de una carga puntual hasta el orden l=4l=4, calcula ahí el error a 20 cm, compáralo con el del eje a esa misma distancia y explica por qué esa dirección no depende del radio del anillo. (e) Escribe el desarrollo válido dentro del anillo (r<Rr<R) y evalúa el centro. ¿Qué rama del desarrollo general de polinomios de Legendre sobrevive ahí, y por qué?

Pista

Todo el anillo está a θ=90\theta' = 90^\circ, así que el momento de orden ll sale de evaluar un solo polinomio en cero. La trampa es confundir los dos ángulos: el de la fuente entra en el momento y el del observador en la forma angular del término. En el eje el potencial exacto es elemental: desarróllalo en potencias de R/rR/r y tus coeficientes tienen que aparecer uno a uno. Anclas: en el eje a 20 cm salen 449,4 / 435,3 / 436,0 V con uno, dos y tres términos, contra 435,96 V exactos; y el centro de (e) vale kQ/R=kQ/R = 1 798 V sin calcular nada.

Problema 5 ●●● · La carga que casi desaparece al poner la chapa a tierra

Una carga q=10q = 10 nC está a d=2d = 2 cm de un plano conductor infinito conectado a tierra, en la misma configuración del artículo 02.

(a) Lejos del sistema y por encima del plano, ¿qué término del desarrollo manda? Da su coeficiente, evalúa VV a 50 cm sobre el eje y compáralo con la suma exacta de la carga y su imagen: da el error y compruébalo a 1 m para deducir su ley con la distancia. (b) La distribución de carga real de la región es qq más la σ\sigma inducida sobre la chapa. Calcula directamente su carga total y su momento dipolar, sin usar la imagen para nada. Te va a salir la mitad del coeficiente de (a), y las dos cosas no pueden ser a la vez el desarrollo multipolar del mismo sistema: encuentra el error. (c) Monta (r)2σda\int (r')^{2}\sigma\,da sobre el plano y mira qué le pasa. ¿Qué hipótesis del teorema del desarrollo multipolar se está violando, y por qué no se viola nunca en los ejemplos de los artículos? (d) ¿A qué distancia produce esta carga 1 V con la chapa puesta? ¿Y sin ella? Da la razón entre las dos distancias y di qué significa para blindar una carga estática.

Pista

Por encima del plano la solución es la del par carga-imagen: desarróllala lejos y quédate con el primer término que no se anule. La trampa es creer que ese término es el momento dipolar de la distribución real, cuando la σ\sigma inducida está toda en z=0z=0 y no aporta nada a pzp_z. Antes de decidir cuál de los dos manda, pregúntate si la distribución real cabe dentro de alguna esfera. Anclas: a 50 cm el dipolo da 14,38 V frente a 14,40 V exactos, y a 1 m el error se divide por cuatro.

Problema 6 ●●● · Dos cargas positivas que se atraen

Una esfera conductora aislada de R=5R = 5 cm lleva una carga total Q=+20Q = +20 nC. Se le acerca radialmente una carga puntual q=+10q = +10 nC, situada a distancia aa del centro. Las dos cargas son positivas.

(a) Construye el sistema de imágenes. ¿Cuál de los dos teoremas de unicidad autoriza este planteamiento, y por qué el otro no sirve aquí? Di dónde va cada imagen y cuánto vale. (b) Escribe F(a)F(a) y evalúala en a=20a = 20 cm y en a=6,5a = 6{,}5 cm. Da el signo de cada una y el potencial de la esfera en los dos casos. (c) Halla la distancia a la que la fuerza cambia de signo y exprésala como distancia a la superficie. (d) Demuestra que para QqQ\gg q el espesor de esa capa atractiva vale aproximadamente (R/2)q/Q(R/2)\sqrt{q/Q}, comprueba con qué error acierta esa fórmula en el caso de arriba y di qué significa físicamente que ese espesor no llegue a cero por mucho que se cargue la esfera.

Pista

El dato es la carga total de la esfera y no su potencial, y eso decide qué teorema de unicidad se usa y por qué hace falta una segunda imagen. La trampa está en dónde va cada una: la de Kelvin sigue en R2/aR^{2}/a y sólo la nueva puede ir al centro, porque es el único sitio que no estropea la equipotencial. Comprueba que de lejos tu fuerza tiende a la repulsión de Coulomb pura entre qq y QQ. Anclas: a 20 cm sale +44,2 μN y a 6,5 cm, −393 μN.

Fin del módulo II.2. Lo que se ha ganado cabe en una frase: la unicidad es lo que convierte adivinar en demostrar. Con ella, este módulo ha puesto sobre la mesa cuatro maneras de resolver un problema de contorno y las cuatro son legales por el mismo motivo: adivinarla cuando la geometría regala una imagen (método de las imágenes), construirla término a término cuando la frontera es una superficie coordenada (separación de variables), aproximarla con la precisión que haga falta cuando sólo interesa el campo lejano (desarrollo multipolar) y calcularla a la fuerza cuando la geometría no coopera (relajación). Sin el teorema, la primera sería un truco, la segunda no demostraría nada, la tercera no tendría cota de error y la cuarta no sabría a qué converge. En términos del texto guía, los problemas de contorno quedan cubiertos salvo dos que conviene nombrar: los problemas sin simetría azimutal, que exigen los armónicos esféricos YlmY_l^{m} en lugar de los Pl(cosθ)P_l(\cos\theta), y la separación en coordenadas cilíndricas, que aquí no se ha montado.

Lo que falta, dicho con precisión. Todo lo de este módulo ocurre en el vacío, con conductores perfectos y con las cargas quietas. Las dos limitaciones se han notado por dentro: la esfera dieléctrica del artículo 03 tuvo que pedir prestada al módulo II.4 la continuidad de la componente normal de D\mathbf{D}, y el teorema de unicidad se demostró para 2V=0\nabla^{2}V=0, una ecuación que deja de valer en cuanto el campo cambia con el tiempo. Con las cargas quietas, además, la electrostática entera es una historia sobre una sola función escalar: eso se acaba en el módulo siguiente.

El módulo II.3, «Magnetostática», pone las cargas en movimiento y con ello aparece un campo que no es el gradiente de nada: la fuerza de Lorentz, Biot-Savart, la ley de Ampère y el potencial vector. Y la maquinaria de este módulo reaparece allí casi entera, lo cual no es una casualidad sino la razón de haberla montado con este cuidado: cada componente cartesiana del potencial vector cumple una ecuación de Laplace donde no hay corriente, el campo lejano de una espira se obtiene con el mismo desarrollo multipolar y los mismos polinomios de Legendre —sólo que allí el primer término no nulo es el dipolo, porque no existen las cargas magnéticas—, y el teorema de unicidad vuelve a ser lo que autoriza a exhibir una solución en vez de deducirla. La ortogonalidad con la que aquí se despejaban coeficientes de uno en uno no reaparece hasta el II.4: hace falta una frontera que imponga condiciones, y en magnetostática no la hay todavía. Después, con el II.4 y el II.5, los campos entran en la materia y empiezan a cambiar con el tiempo, que es cuando esto se convierte en electrodinámica.