Electrodinámica - Tarea 3

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 5

Facultad de Física

Universidad Veracruzana
Xalapa, Veracruz

Electrodinámica
Tarea 3

Iván Octavio Argüelles Barra


Miranda Gabriela Martínez Parra
Valente Emmanuel Morgado Morales

Profesor: Juan Efraín Rojas Marcial


Tarea 3 Electrodinámica


1. Al considerar el potencial electrostático, φ(⃗r′ ) = k p⃗|⃗·(⃗
r−⃗
r−⃗
r)
r ′ |3
, producido por un
dipolo p⃗, realiza un gráfico de las trazas de las superficies equipotencia-
les en un plano que contiene al dipolo. Para mayor comodidad, el dipolo
puede estar ubicado en el origen. Utiliza los resultados obtenidos para es-
bozar algunas de las líneas de fuerza. Puedes utilizar Mathematica, Maple
o MATLAB.
Solución.

a) Parámetros del Dipolo: En esta línea, definimos el momento dipolar del dipolo
eléctrico. Establecemos un dipolo con un momento dipolar de (1, 0, 0) en la
dirección x.
1 p = [1 , 0 , 0]; % Momento dipolar en la d i r e c c i n
x

b) Diseño de una Cuadrícula: Aquí, creamos una cuadrícula de puntos en el plano


XY, que abarca desde -4 hasta 4 en pasos de 0.5 en ambas direcciones.
1 [x , y ] = meshgrid ( -4:0.5:4 , -4:0.5:4) ;

c) Definición de la constante de Coulomb: Establecemos la constante de Coulomb


en el valor apropiado (8.99 × 109 N m2 /C 2 ).
1 k = 8.99 e9 ;

d ) Cálculo del Potencial Eléctrico: Aquí, calculamos la magnitud de la posición (r)


para cada punto en la cuadrícula. Luego, inicializamos una matriz de ceros (phi)
para almacenar el potencial eléctrico en cada punto.
1 r = sqrt ( x .^2 + y .^2) ; % Magnitud de la p o s i c i n
2 phi = zeros ( size ( x ) ) ; % Inicializar matriz de
potencial

Utilizamos un bucle for para calcular el potencial eléctrico en cada punto de la


cuadrícula utilizando la fórmula electrostática proporcionada. La fórmula se aplica
para cada punto y se almacena en la matriz phi.
1 for i = 1: numel ( x )
2 phi ( i ) = k * dot (p , [ x ( i ) , y ( i ) , 0]) / r ( i ) ^3;
3 end

e) Creación de Superficies Equipotenciales: Aquí creamos una figura y trazamos las


superficies equipotenciales utilizando la función contour. Estas líneas representan
regiones donde el potencial eléctrico es constante.
1 figure ;
2 contour (x , y , phi , 20) ;

f ) Cálculo del Campo Eléctrico: Calculamos el campo eléctrico en cada punto de la


cuadrícula utilizando la función gradient. Esto nos da las componentes Ex y Ey
del campo eléctrico.
1 figure ;
2 [ Ex , Ey ] = gradient ( - phi , 0.2 , 0.2) ;

1
Tarea 3 Electrodinámica

g) Dibujo de Líneas de Campo Eléctrico: Dibujamos las líneas de campo eléctrico


utilizando la función quiver. Las flechas representan la dirección y magnitud del
campo eléctrico en cada punto. Las flechas se escalan y se dibujan en negro.
1 hold on ;
2 quiver (x , y , scale * Ex , scale * Ey , ’k ’) ;
3 hold off ;

h) Ajustes Gráficos y Configuración del Mapa de Colores: Aseguramos que el gráfico


tenga ejes iguales, activamos la cuadrícula y configuramos un mapa de colores
personalizado para las superficies equipotenciales (azul para el lado negativo y
rojo para el lado positivo). También agregamos una barra de colores para indicar
la dirección del campo eléctrico.
1 axis equal ;
2 grid on ;
3 colormap ([0 0 1; 1 0 0]) ; % Azul ( lado negativo ) y
rojo ( lado positivo )
4 colorbar ;

Así, obtenemos el siguiente script, que ejecutado en Matlab arroja la siguiente imagen:
1 % P a r m e t r o s del dipolo
2 p = [1 , 0 , 0]; % Momento dipolar en la d i r e c c i n x
3

4 % Crear una c u a d r c u l a de puntos en el plano


5 [x , y ] = meshgrid ( -4:0.5:4 , -4:0.5:4) ;
6

7 % Constante de Coulomb
8 k = 8.99 e9 ;
9

10 % Calcular el potencial e l c t r i c o en cada punto de la


cuadr cula
11 r = sqrt ( x .^2 + y .^2) ; % Magnitud de la p o s i c i n
12 phi = zeros ( size ( x ) ) ; % Inicializar matriz de potencial
13

14 % Calcular el potencial para cada punto de la c u a d r c u l a


15 for i = 1: numel ( x )
16 phi ( i ) = k * dot (p , [ x ( i ) , y ( i ) , 0]) / r ( i ) ^3;
17 end
18

19 % Crear g r f i c o de las superficies equipotenciales


20 contour (x , y , phi , 20) ;
21 xlabel ( ’X ’) ;
22 ylabel ( ’Y ’) ;
23 title ( ’ L n e a s de Fuerza de un dipolo en el origen ’) ;
24

25 % Calcular el campo e l c t r i c o en cada punto


26 [ Ex , Ey ] = gradient ( - phi , 0.2 , 0.2) ;
27

28 % Normalizar los vectores del campo e l c t r i c o


29 E_mag = sqrt ( Ex .^2 + Ey .^2) ;

2
Tarea 3 Electrodinámica

30 Ex = Ex ./ E_mag ;
31 Ey = Ey ./ E_mag ;
32

33 % Escalar el t a m a o de las l n e a s de campo e l c t r i c o


34 scale = 0.2;
35

36 % Dibujar l n e a s de campo e l c t r i c o
37 hold on ;
38 quiver (x , y , scale * Ex , scale * Ey , ’k ’) ;
39 hold off ;
40

41 axis equal ;
42 grid on ;

3
Tarea 3 Electrodinámica

2. Suppose a molecule is represented by a charge −2q at the origin of coordi-


nates and charges +q at l⃗1 and l⃗2 with |l⃗1 | = |l⃗2 | = l

a) Find the dipole moment of this molecule


b) For H2 O, consider l = 0.958 × 10−10 m and the angle between l⃗1 and l⃗2 being
θ = 105o . If p = 6.14 × 10−30 C · m, find the effective charge q.

Solución.

Las dos cargas +q pueden colocarse de tal manera


que el ángulo θ entre los vectores de posición de las
cargas, l⃗1 y l⃗2 , quede bisectado por uno de los dos ejes,
en este caso, el eje x, como se puede ver en la figura
1. Así, sean l⃗1 = x1 ı̂ + y1 ȷ̂ y l⃗2 = x2 ı̂ + y2 ȷ̂, donde

θ θ
x1 = l1 cos y1 = −l1 sen
2 2
θ θ
x2 = l2 cos y2 = l2 sen
Figura 1: Configuración de las cargas 2 2

Así, el momento dipolar es

p⃗ = q l⃗1 + q l⃗2 = q(x1 ı̂ + y1 ȷ̂) + q(x2 ı̂ + y2 ȷ̂)


   
θ θ θ θ
= q l1 cos ı̂ − l1 sen ȷ̂ + q l2 cos ı̂ + l2 sen ȷ̂
2 2 2 2
θ
= q(l1 + l2 ) cos ı̂
2
Pero recordando que l1 = l2 = l, entonces, finalmente se obtiene

θ
p⃗ = 2ql cos ı̂ (1)
2

Ahora, tomando únicamente la magnitud de p⃗ y despejando q de la ecuación (1), se obtiene


p
q=
2l cos 2θ

Para la molécula H2 O, usando l = 0.958 × 10−10 m, θ = 105o y p = 6.14 × 10−30 C · m

6.14 × 10−30 C · m −20


q= −10 105 ≈ 5.26 × 10 C
2(0.958 × 10 m) cos 2

En conclusión, la carga efectiva es

q = 5.26 × 10−20 C

También podría gustarte