Instituto Politécnico Nacional

Escuela Superior de Ingeniería Mecánica y Eléctrica

Unidad Zacatenco

Maestría en Ciencias en Ingeniería Eléctrica

Fundamentos de la Teoría Electromagnética.

Práctica 1: Campo electrostático en una línea coaxial.

Profesor: Dr. Fermín Espino Cortes.

Alumno: Julio César Salazar Hernández

Boleta: B171057

Fecha: 09/02/2018
 Práctica No. 1: “Campo electrostático en una línea
coaxial”.
Objetivo.
Calcular el campo eléctrico y el potencial eléctrico en un cable coaxial utilizando las
ecuaciones de maxwell en el caso electrostático de manera analítica y simulada empleando
COMSOL Multiphysics.

Marco teórico.
Un cable coaxial es una línea de transmisión compuesta de dos conductores, uno interno y
otro externo (ver fig. 1.1). El conductor externo es un tubo metálico o blindaje de malla, más
sin en cambio el conductor interno se sostiene por un dieléctrico que puede ser un sólido,
plástico expandido o semisólido (discos de polietileno, cintas helicoidales, tiras de plástico,
etc.).
Esta línea de transmisión se caracteriza por ser una estructura blindada, es decir, el campo
electromagnético está limitado entre el espacio de los conductores interior y exterior.

Ilustración 1.1.- Estructura de un Cable Coaxial.

Un coaxial tiene un diámetro externo b y uno interno a. El espacio entre los conductores está
formado por un dieléctrico con cierta permitividad. Se aplica una diferencia de potencial
entre los conductores donde el externo tiene 0 v y el interno V0 (ver fig. 1.2).
 Ilustración 1.2.- Radio interno y externo del cable coaxial separados por un dieléctrico con cierta permitividad.

Cálculo analítico.
Cálculo del campo eléctrico. La ecuación de campo eléctrico para campos electrostáticos se
utilizan la ecuación de Poisson (1) y la de Laplace (2) en función de campos.
𝜌
∇2 𝐸 = ∇ [ 𝜀 ]………………………………………………………………(1)

∇2 𝐸 = 0…………………………………………………………………..(2)
Ahora en función de potenciales se tiene:
𝜌
∇2 𝑉 = − 𝜀 …………………………………………………………………(1.1)

∇2 𝑉 = 0……………………………………………………………………(2.1)
En la región de a<r<b se encuentra el dieléctrico por lo que no hay carga libre es decir 𝜌=0,
por lo que se emplea la ecuación 2.1 la ecuación de Laplace de potenciales. Por lo que las
condiciones iniciales son:
𝑉(𝜌 = 𝑎) = 𝑉0
𝑉(𝜌 = 𝑏) = 0
Este problema se resuelve directamente porque la estructura es cilíndrica y simétrica. Girando
el coaxial sobre el eje z no cambia la geometría en nada. Como resultado sabemos que el
campo de potencial eléctrico es función de 𝜌 solo:
Al resolver la ecuación 2.1 en coordenadas cilíndricas se tiene la siguiente ecuación:
1 𝜕 𝜕𝑓 1 𝜕 2𝑓 𝜕 2𝑓
∇2 𝑓(𝑟, ∅, 𝑧) = [𝑟 ] + 2 2 + 2
𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝜕∅ 𝜕𝑧
En este caso en la dirección z y ∅ no hay cambios, es decir, siempre son constantes. Por lo
tanto, lo que nos interesa es saber el campo eléctrico en cada punto desde el radio a hasta el
b. cambiando la ecuación con las variables antes usadas se tiene:
∇2 𝑉(𝜌) = 0
1 𝜕 𝜕𝑉(𝜌)
(𝜌 )=0
𝜌 𝜕𝜌 𝜕𝜌
𝜕 𝜕𝑉(𝜌)
(𝜌 )=
𝜕𝜌 𝜕𝜌
0…………………………………………………………………………………..(2.2)
Al integrar los dos lados de la ecuación anterior.
𝜕 𝜕𝑉(𝜌)
∫ (𝜌 ) 𝑑𝜌 = ∫ 0𝑑𝑝
𝜕𝜌 𝜕𝜌
𝜕𝑉(𝜌)
𝜌 = 𝐶1
𝜕𝜌
𝜕𝑉(𝜌) 𝐶1
=
𝜕𝜌 𝜌
Integrando otra vez se tiene lo siguiente:
𝜕𝑉(𝜌) 𝐶1
∫ 𝑑𝜌 = ∫ 𝑑𝑝
𝜕𝜌 𝜌
𝑽(𝝆) = 𝑪𝟏 𝒍𝒏[𝝆] + 𝑪𝟐
Esta ecuación satisface la ecuación 2.2
𝜕 𝜕𝑉(𝜌)
(𝜌 )=0
𝜕𝜌 𝜕𝜌
𝜕 𝜕
(𝜌 (𝐶1 𝑙𝑛[𝜌] + 𝐶2 )) = 0
𝜕𝜌 𝜕𝜌
𝜕 1
(𝜌(𝐶1 )) = 0
𝜕𝜌 𝜌
0=0
Ahora se evalúan las condiciones iniciales para conocer los valores de las dos constantes:
𝑉(𝜌 = 𝑎) = 𝐶1 ln(𝑎) + 𝐶2 = 𝑉0…………(A)
𝑉(𝜌 = 𝑏) = 𝐶1 ln(𝑏) + 𝐶2 = 0…………...(B)
Resolviendo el sistema de ecuaciones de dos incógnitas, despejano B:
 𝐶2 = −𝐶1 ln(𝑏)
Sustituyendo en A
𝐶1 ln(𝑎) − 𝐶1 ln(𝑏) = 𝑉0
𝐶1 ln(𝑏) − 𝐶1 ln(𝑎) = −𝑉0
−𝑽𝟎
𝑪𝟏 =
𝒍𝒏[𝒃/𝒂]
Sustituyendo c1 en c2:
𝑽𝟎 𝐥𝐧(𝒃)
𝑪𝟐 =
𝒍𝒏[𝒃/𝒂]
Por lo tanto, la solución final es:
−𝑽𝟎 𝒍𝒏[𝝆] 𝑽𝟎 𝐥𝐧(𝒃)
𝑽(𝝆) = + ;𝒂 < 𝝆 < 𝒃
𝒍𝒏[𝒃/𝒂] 𝒍𝒏[𝒃/𝒂]
La comprobación se da sustituyendo los valores de frontera.
−𝑉0 𝑙𝑛[𝑎] 𝑉0 ln(𝑏)
𝑉(𝜌 = 𝑎) = +
𝑙𝑛[𝑏/𝑎] 𝑙𝑛[𝑏/𝑎]
𝑉0 (𝑙𝑛[𝑏] − 𝑙𝑛[𝑎])
𝑉(𝜌 = 𝑎) =
𝑙𝑛[𝑏/𝑎]
𝑉0 𝑙𝑛[𝑏/𝑎]
𝑉(𝜌 = 𝑎) =
𝑙𝑛[𝑏/𝑎]
𝑽(𝝆 = 𝒂) = 𝑽𝟎
−𝑉0 𝑙𝑛[𝑏] 𝑉0 ln(𝑏)
𝑉(𝜌 = 𝑏) = +
𝑙𝑛[𝑏/𝑎] 𝑙𝑛[𝑏/𝑎]
𝑽(𝝆 = 𝒃) = 𝟎
Ahora que se tienen el potencial eléctrico podemos encontrar el campo eléctrico como
sigue:
𝐸̅ = −∇𝑉
𝜕 −𝑉0 𝑙𝑛[𝜌] 𝑉0 ln(𝑏)
𝐸(𝜌) = − [ + ]
𝜕𝜌 𝑙𝑛[𝑏/𝑎] 𝑙𝑛[𝑏/𝑎]
𝑽𝟎 𝟏
𝑬(𝝆) = ̂; 𝒂 < 𝝆 < 𝒃
𝒂
𝒍𝒏[𝒃/𝒂] 𝝆 𝝆
También se puede saber la densidad flujo eléctrico:
 ̅ = 𝜀𝐸̅
𝐷
𝐷(𝜌) = 𝜀𝐸(𝜌)
𝜺𝑽𝟎 𝟏
𝑫(𝝆) = ̂; 𝒂 < 𝝆 < 𝒃
𝒂
𝒍𝒏[𝒃/𝒂] 𝝆 𝝆
Una vez resuelto el caso electrostático se evaluaron las ecuaciones con valores entre el radio
interno hasta el externo. Los datos propuestos del cable coaxial son los siguientes:
𝑟𝑖 = 0.022[𝑚]
𝑟𝑒 = 0.045[𝑚]
𝜀𝑟 = 2.3
𝑉0 = 50𝑘[𝑉]
En la tabla 1,2 y 3 se muestran los valores obtenidos de la densidad de campo eléctrico, el
potencial eléctrico y el campo eléctrico respectivamente.
Tabla 1.- Densidad de campo eléctrico en el cable coaxial.

𝝆 𝑫(𝝆)
0.022 6.46x10-5
0.025 5.68x10-5
0.030 4.74x10-5
0.032 4.44x10-5
0.035 4.06x10-5
0.040 3.55x10-5
0.042 3.38x10-5
0.045 3.16x10-5

Tabla 2.-Potencial Eléctrico en el cable coaxial.

𝝆 𝑽(𝝆)
0.022 50000.00
0.025 41068.34
0.030 28329.63
0.032 23820.36
0.035 17559.20
0.040 8229.43
0.042 4820.49
0.045 0
 Tabla 3.-Campo Eléctrico en el cable coaxial.

𝝆 𝑬(𝝆)
0.022 3175885.466
0.025 2794779.21
0.030 2328982.675
0.032 2183421.258
0.035 1996270.864
0.040 1746737.006
0.042 1663559.054
0.045 1552655.117

Simulación utilizando COMSOL Multiphysic.

Al resolver el problema como caso electrostático y dar las dimensiones y características
físicas del problema se obtuvieron las siguientes graficas de campo eléctrico y potencial
eléctrico. En la figura 1.3 se muestra la gráfica de la norma de campo eléctrico que va de un
poco más de 3M0[ V/m] hasta un poco más de 1M [V/m]. Mas sin en cambio en la figura 1.4
se aprecia el potencial eléctrico que va desde 50k [V] hasta 0 [V].

Ilustración 1.3.- Gráfica del campo eléctrico en el cable coaxial.

 Ilustración 1.4.- Gráfica del potencial eléctrico en el cable coaxial.

Cálculo analítico vs simulado.

Para poder comparar los resultados obtenidos se procedió a graficar las ecuaciones analíticas
en Excel. Las gráficas de potencial eléctrico (fig. 1.4 y 1.5) y la de campo eléctrico (fig. 1.3
y 1.6) son muy parecidas ya que se utilizó un mallado extra fino al simularlo.

Tensión (V)
60000

50000

40000

30000

20000

10000

0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
-10000

Ilustración 1.5.- Gráfica de potencial eléctrico analítico en el cable coaxial.

 Campo Eléctrico (V/m)
3500000

3000000

2500000

2000000

1500000

1000000

500000

0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24

Ilustración 1.6.- Gráfica de campo eléctrico analítico del cable coaxial

Para poder visualizar mejor estas dos graficas COMSOL nos permite obtener la distribucion
del potencial electrico y el campo electrico. En la figura 1.7 se observa que el potencial
electrico es mas intenso en la superficie del conductor interno y se va reduciendo hasta cero
cuando llega al conductor externo.

Ilustración 1.7.- Distribución del potencial eléctrico en el cable coaxial.

En la figura 1.8 se aprecia la distribución del campo eléctrico y se comporta semejante al

potencial eléctrico. El campo eléctrico es más intenso en la superficie del conductor.
 Ilustración 8.- Distribución del campo eléctrico en el cable coaxial.

Conclusión.
En la solución de este problema electrostático se tuvo la ventaja que el problema se pudo
reducir a una simple ecuación diferencial debido a que solo variaba el radio del cable coaxial.
Si fuese el caso contrario quizás la solución analítica se volvería más compleja y difícil de
resolver ya que, aunque se obtengan las ecuaciones diferenciales que representen el
comportamiento físico del problema no garantiza que exista una solución analítica por lo que
se opta por emplear Software de simulación que utiliza métodos numéricos para dar una
solución aproximada.
En el caso de COMSOL Multiphysic utiliza un mallado para poder resolver el problema, si
emplea un mallado muy fino los resultados son muy cercanos al analítico, pero se lleva más
tiempo de computo, si se utiliza un mallado pobre los resultados son alejados al analítico y
emplea menor tiempo de computo.
La distribución del potencial eléctrico y el campo eléctrico es muy parecida ya que la mayor
parte se encuentra en la superficie del conductor central y va disminuyendo su intensidad
hasta desaparecer entre más se acerca al conductor externo.

Referencias.