UNIVERSIDAD DE ZARAGOZA

FACULTAD DE CIENCIAS
DEPARTAMENTO DE FSICA APLICADA
REA DE ELECTROMAGNETISMO

ESTUDIO DE LA RELAJACIN DIELCTRICA

EN MEZCLAS POLAR-NO POLAR POR
REFLECTOMETRA EN DOMINIO DEL
TIEMPO (T.D.R.)

TRABAJO ACADMICAMENTE DIRIGIDO

CURSO 2005-06

ALUMNO: DANIEL DIGN RODRGUEZ

DIRECTOR: JUAN PABLO MARTNEZ JIMNEZ

Agradecimientos.

Quisiera expresar mi agradecimiento:

A mi director de TAD Juan Pablo Martnez y a mi Profesor Jos Mara

Fornis Marquina. No slo por su apoyo y ayuda incondicionales, sino
tambin por la innumerable cantidad de explicaciones y buenos consejos
que me han brindado, y que, sin duda alguna, me han servido para
completar mi formacin como fsico.

Al departamento de Fsica Aplicada, por proporcionarme todas las

herramientas necesarias para la elaboracin del trabajo.

A mis compaeros de 5 de la Licenciatura en Fsica y, por supuesto, a

mi familia.

ndice
Pgina
1.- Introduccin.

2.- Fsica de dielctricos.

2.1.- Comportamiento de la materia bajo la accin de un campo electrosttico.
2.2.- Comportamiento dinmico de la materia.
2.2.1.- En el dominio de la frecuencia.
2.2.2.- En el dominio temporal.
2.3.- Relaciones de dispersin.
2.3.1.- Modelo de Debye para fenmenos de relajacin.
2.3.2.- Correcciones del modelo de Debye.
2.4.- Leyes de mezclas.

9
9
11
12
13
14
15
20
22

3.- Teora de las lneas de transmisin.

3.1.- Modelo circuital de una lnea de transmisin.
3.1.1.- Ecuaciones del telegrafista.
3.1.2.- Lnea de transmisin finita.
3.2.- Modelo electromagntico de una lnea de transmisin.
3.2.1.- Cable coaxial.

25
25
25
28
30
32

4.- Tcnica de medida.

4.1.- Descripcin del montaje experimental
4.2.- Funcin de transferencia del sistema de medida.
4.3.- Tratamiento de las medidas.
4.4.- Clculo de la constante dielctrica.
4.5.- Correccin de la arandela.

35
35
39
42
44
45

5.- Resultados experimentales.

5.1.- La permitividad esttica.
5.2.- Absorcin mxima.
5.3.- El tiempo de relajacin.

47
51
54
54

6.- Conclusiones

59

Bibliografa

61

Apndice: Cdigo del programa en lenguaje C de la aplicacin tdr.exe

63

1.- Introduccin.
A lo largo de la historia de la humanidad se ha generalizado el uso de nuevos
lquidos como pinturas, carburantes, refrigerantes, lubricantes..., tanto ms variados y
complejos cuanto ms desarrollada se encuentra nuestra sociedad. En la actualidad, no
cesa la necesidad de nuevos materiales dielctricos en reas muy diversas. Las mezclas
de lquidos nos proporcionan un amplio abanico de posibilidades para la
implementacin previa de sistemas que presenten una respuesta electromagntica
especfica.
Una de las magnitudes utilizadas para la caracterizacin de mezclas dielctricas,
es la permitividad dielctrica en funcin de la frecuencia. En medios no magnticos, nos
informa totalmente de su respuesta electromagntica. La espectroscopia dielctrica
permite conocer parmetros fundamentales desde el punto de vista de la investigacin
cientfica bsica y desde un punto de vista tecnolgico, conocer la respuesta del material
a cualquier seal electromagntica.
El objetivo de este trabajo es la caracterizacin de mezclas binarias polar-no
polar. Mediante reflectometra en el dominio del tiempo (Time Domain Reflectometry,
TDR) se calcula la permitividad compleja de la mezcla en un rango de frecuencias
DC-5 GHz y se determina el tiempo de relajacin. Las componentes polares son
pentanol, hexanol y heptanol y la no polar es el ciclohexano. Los rangos de
concentracin de alcohol en las muestras a medir son 0.05, 0.10 y 0.15 (fracciones
molares).
En el laboratorio ya se han realizado medidas a concentraciones superiores, de
manera que este Trabajo Acadmicamente Dirigido (T.A.D.) completa su estudio. Para
analizar los resultados se utilizarn leyes dielctricas de mezcla existentes en la
literatura que se basan en las relaciones de concentracin entre sus componentes.

2.- Fsica de dielctricos

En este captulo se van a presentar las magnitudes que caracterizan un medio
dielctrico, y su comportamiento bajo la accin de un campo variable. Tambin se
realiza un estudio de los procesos de relajacin y se introduce el modelo de Debye y sus
correcciones.
Se termina el captulo presentando unas leyes de mezclas dielctricas que
permiten estimar el valor de la permitividad esttica de la mezcla.

2.1.- Comportamiento de la materia bajo la accin de un campo

electrosttico.
La materia est constituida por partculas con carga. En los dielctricos las
cargas estn ligadas entre s, mientras que en los metales, algunos electrones cuyas
energas pertenecen a la banda de conduccin, no estn totalmente ligados a los ncleos
de los tomos, y pueden moverse por el metal. Los campos creados por estas cargas
sufren variaciones espacio-temporales muy grandes con perodos espaciales de 1 , y
temporales del orden de 10-15 s. Estas variaciones no son apreciables desde un punto de
vista macroscpico, as que las nicas variaciones observables, son las producidas por
los campos externos.
El comportamiento elctrico lo podremos caracterizar por su momento dipolar si
suponemos que el medio es neutro y despreciando los momentos de orden superior. A
este tratamiento se le denomina aproximacin dipolar [1]. Es equivalente a considerar el
dielctrico como un conjunto de dipolos.
Se define el vector polarizacin como el momento dipolar por unidad de
volumen:

v dpr
;
P=
dv

(2-1)

En virtud de cules son los procesos que intervienen en la aparicin de la

polarizacin, se hace una clasificacin:
Polarizacin orientacional. Se debe a la reordenacin de los dipolos que componen

un medio sometido a la accin de un campo elctrico. Aparece en sustancias

constituidas por molculas con momento dipolar permanente. En ausencia de
campo, las molculas estn dispuestas de forma aleatoria debido a la agitacin
trmica, por tanto su polarizacin ser nula. Si a continuacin se aplica un campo
esttico, al cabo de un cierto tiempo se alcanza una situacin de equilibrio en la que
habr mayor nmero de dipolos orientados en la direccin del campo y por tanto la
polarizacin ser no nula.
Polarizacin inducida. Los dipolos elctricos aparecen como consecuencia de la

modificacin de la distribucin de cargas al aplicar un campo elctrico al medio. Si

este tipo de polarizacin aparece a nivel atmico como consecuencia del
desplazamiento relativo entre la nube electrnica y el ncleo del tomo, se le
denomina polarizacin atmica electrnica. Cuando ocurre como consecuencia de
la deformacin de una red inica se la denomina polarizacin inica.
Para un medio istropo, la relacin entre el momento dipolar y el campo
elctrico local (es la superposicin del externo y los campos microscpicos debidos las
cargas del medio) puede expresarse como un desarrollo en serie de potencias:
r
r
r
p = E loc + E loc

r 3
+ E loc + K

(2-2)

En la mayora de los casos, podemos quedarnos slo con el primer trmino de

esta expresin. Entonces la polarizacin que es el promedio de los momentos de
polarizacin en el material por unidad de volumen, slo depender del valor promedio
del campo local, es decir, del campo externo:
r
r
P = 0e E ;

10

(2-3)

donde 0 es la permitividad dielctrica del vaco, y e la susceptibilidad elctrica que

es una caracterstica propia del medio.
Teniendo en cuenta la definicin de vector desplazamiento elctrico:
r
r r
D = 0E + P ;

(2-4)

y sustituyendo la expresin de la polarizacin (2-3), se obtiene:

r
r
r
D = 0 (1 + e )E = E

(2-5)

donde es la permitividad dielctrica del medio constante dielctrica absoluta. Si el

medio es el vaco no existirn cargas sobre las que se genere momento dipolar, por lo
que la polarizacin del medio es nula y la permitividad es la del vaco.
Normalmente trabajaremos con la permitividad dielctrica relativa, que se define
como:
r =

= e + 1;
0

(2-6)

en el caso de campo aplicado esttico, denotaremos la permitividad relativa como s .

2. 2.- Comportamiento dinmico de la materia.

Cuando el campo elctrico externo es variable en el tiempo, la respuesta del
medio no ser instantnea, y por tanto el vector polarizacin en un instante de tiempo,
depender del campo en dicho instante y en instantes anteriores.
En el dominio de la frecuencia, esto se traduce en una respuesta del medio
dependiente de la frecuencia del campo.

11

2.2.1.- Dominio de la frecuencia.

Cualquier funcin dependiente del tiempo se puede expresar como superposicin
de funciones armnicas segn el anlisis de Fourier. Por simplicidad, supondremos que
el campo elctrico aplicado es una funcin armnica de frecuencia :
r
r
E(t ) = E 0 e jt ;

(2-9)

Si el comportamiento del medio es lineal, el vector polarizacin ser tambin

una funcin armnica, pero la respuesta no tiene por qu ser instantnea, y es posible la
existencia de un desfase entre ambos:
r
r
P(t ) = P0 e ( jt + ) ;

(2-10)

La relacin entre ambas magnitudes vectoriales es:

r
r
P0 = 0 e E 0 ;

(2-11)

donde ahora e es un nmero complejo dependiente de la frecuencia. En general:

r
r
P() = 0 e ()E () ;

(2-12)

Y el vector desplazamiento en funcin de la frecuencia es:

r
r
r
D() = 0 1 + e () E() = ()E()

(2-13)

En este caso, la permitividad tambin es un nmero complejo que normalmente

se expresa en forma binomial:

() = () j()

12

(2-14)

En el ao 1926, Kramers y Krnig [1] por separado, demostraron que () y

() son dependientes la una de la otra. Tales relaciones conocidas como Kramers-

Krnig son las siguientes:

2
' () = 0 +

2
() =

z (z )
dz
z 2 2

' (z ) 0
z 2 2

dz

(2-15a)

(2-15b)

2.2.2.- Dominio temporal.

Un campo elctrico variable en el tiempo se puede expresar como una sucesin
de pulsos de diferente amplitud. Debido al comportamiento lineal del medio, la
polarizacin originada tambin podr expresarse como una sucesin de pulsos. Como el
medio no siempre responde de forma instantnea, el valor de la polarizacin en un
instante de tiempo t depender de las amplitudes de los pulsos del campo en dicho
instante t y en tiempos anteriores:

r
P (t ) = 0

r
(t ')E(t t ')dt '

(2-16)

El vector desplazamiento en funcin del tiempo que se obtiene utilizando la

expresin anterior es el siguiente:

r
r
r
r
D(t ) = 0 E(t ) + P(t ) = 0 E(t ) + 0

13

r
(t ')E(t t ')dt '

(2-17)

De los dos trminos de (2-17), el primero indica una respuesta instantnea e

independiente del medio, y el segundo es la contribucin de la respuesta del medio,
caracterizado por la susceptibilidad dielctrica temporal (t) al campo elctrico en el
instante t y en los anteriores.

2.3.- Relaciones de dispersin.

La relacin dispersin es la variacin de la permitividad, tanto de ' como de

" , con la frecuencia. En la curva de dispersin de un dielctrico normalmente se

pueden distinguir dos zonas: una en la que la parte imaginaria de la permitividad es nula
y, como predicen las relaciones de Kramers y Krnig, la parte real es constante; y otra
en la que " es positiva y por tanto ' vara. En esta ltima zona se producen prdidas y
hay dos tipos de fenmenos que las producen:
Fenmenos de relajacin debidos a la reorientacin de los dipolos permanentes
que constituyen el material dielctrico al estar sometidos a un campo oscilante.
Los tiempos de respuesta son grandes ( > 15 ps )

Fenmenos de resonancia. Se caracterizan por tener tiempos de respuesta ms

cortos. Se debe al comportamiento dinmico de la polarizacin inducida sobre el
medio.
En la figura 2.1. se muestra la curva de dispersin de un material dielctrico con
dos relajaciones y cuatro resonancias.

14

Figura 2.1. Ejemplo de curva de dispersin de un material dielctrico.

2.3.1.- Modelo de Debye para fenmenos de relajacin.

Supongamos que tenemos una mezcla binaria constituida por molculas polares,
con momento dipolar permanente y por molculas no polares. En este caso tendremos
dos contribuciones a la polarizacin: la debida a los dipolos y la asociada a la
polarizacin electrnica, que consideraremos que tiene una respuesta instantnea, de
r
modo que se denotar por P . Supongamos que en el instante t = 0 se aplica un campo
constante, de manera que el campo elctrico se puede expresar como una funcin
escaln en funcin del tiempo:
r r
E = E 0 (t )

(2-18)

0 t < 0
donde (t ) =
es la funcin escaln normalizada.
1 t 0

Entonces, suponiendo que la polarizacin del medio tiene una respuesta descrita
por una ecuacin diferencial de primer orden:
r r r
P Ps P
=
t

(2-19)

15

La polarizacin total en funcin del tiempo ser:

r
r
r r
P(t ) = P + Ps P

)1 e t / (t )

(2-20)

El vector polarizacin tiene dos trminos: el instantneo, y el retardado, que se

corresponden con las dos contribuciones comentadas anteriormente. En la figura 2.2. se
muestra una grfica donde se representa el mdulo del vector polarizacin frente al
tiempo.

Ps

P
Figura 2.2. Evolucin temporal del mdulo del vector polarizacin bajo la accin de un
campo electrosttico.

es el tiempo de relajacin cuya dependencia con la temperatura, viene dada por la ley

de Arrhenius:
= 0 eW/kT

(2-21)

donde W es la energa de activacin necesaria para alcanzar la orientacin del campo,

k es la constante de Boltzmann y 0 el valor al que tiende el tiempo de relajacin a

temperaturas muy altas.
Los dos trminos de la polarizacin pueden expresarse en funcin del campo
aplicado, de manera que:
r
r
r
r r
P(t ) = ( 0 ) E 0 + ( s )1 e t / E 0 = Pi + Po

16

(2-22)

v
Teniendo en cuenta esa expresin para P , el desplazamiento elctrico se puede

escribir como:
r
r
r
D(t ) = (t ) E 0 + ( s )1 e t / E 0

(2-23)

y derivando con respecto al tiempo, se obtiene la corriente de desplazamiento:

r
r
D
r
1
= E 0 (t ) + ( s ) e -t/ E 0
t

(2-24)

Dado que la expresin de la permitividad en funcin de la frecuencia es [4]:

( ) = + ( s ) (t )e j t dt

(2-25)

podemos identificar la funcin respuesta como (t ) =

e t /
, y se obtiene la siguiente

expresin para la permitividad:

() = +

s
1 + j

(2-26)

que es la expresin de Debye.

Despejando la parte real e imaginaria, se obtiene:

' () = +

" () =

s
1 + 2 2

( s )
1 + 2 2

(2-27a)

(2-27b)

La figura 2.3. muestra la permitividad relativa frente a para un proceso de

relajacin segn el modelo de Debye.

17

Figura 2.3. Relajacin en el modelo de Debye.

Si representamos la permitividad en el plano complejo, es decir, la parte imaginaria en

el eje de ordenadas y la parte real en el de abcisas, obtenemos el llamado diagrama de
Argand de Cole-Cole:

Figura 2.4. Diagrama de Argand en el modelo de Debye.

Que para el modelo de Debye es una semicircunferencia de radio

+
,0 ; El mximo de se obtiene en:
en el punto ( , ) = s
2

18

s
, con centro
2

1
;

0 =

(2-28)

que es la frecuencia de relajacin. Se despeja en las ecuaciones (2-27) el factor

s
1 + 2 2

, y tras igualar los resultados, se obtiene:

(2-29)

Sustituyendo esta expresin en (2-27a), llegamos a la ecuacin de una circunferencia:

( )(1 + 2 2 ) = s

( ) +

2
=
( ) s

2 = ( )( s )

+ s

= s
2
2

(2-30)

que es la que se representa en el Diagrama de Argand de Cole-Cole. As que

representando en el plano complejo la permitividad compleja experimental de un
dielctrico,

podremos

comprobar

geomtricamente

cuanto

se

aproxima

su

comportamiento al modelo de Debye.

Otra posibilidad a la hora de representar la permitividad compleja, son los
diagramas lineales de Cole, que nos permiten calcular la permitividad esttica e infinita,
y el tiempo de relajacin. Despejando de la expresin (2-29):
1
= +

19

(2-31)

1

y representando frente a , se obtiene la ecuacin de una recta de pendiente y

ordenada en el origen . Este es el llamado diagrama lineal de Cole para permitividad

en alta frecuencia Cole infinita.
Si ahora se despeja de la ecuacin (2-27a), y se sustituye en (2-27b), se
obtiene la ecuacin de otra recta, tomando como variable independiente , y como
variable dependiente :

(1 + )
2 2

2 2

= s ()

(2-32)

Este es el llamado diagrama lineal de Cole para permitividad esttica.

Si representamos los valores experimentales de frente a , podemos
ajustar estos puntos a una recta, cuya pendiente, nos dar el valor del tiempo de
relajacin y la ordenada en el origen el valor de la permitividad esttica.
Cuanto ms se ajusten los puntos experimentales a estos diagramas lineales, ms
se aproximar el modelo de Debye al proceso de relajacin medido.

2.3.2.- Correcciones del modelo de Debye.

No todos los medios dielctricos siguen un comportamiento que se pueda
describir con el modelo de Debye. La causa puede ser la existencia de una distribucin
continua discreta de tiempos de relajacin. Esto sucede cuando la frecuencia de
relajacin depende del eje de rotacin en molculas que carecen de simetra esfrica.
A continuacin se muestran dos modelos alternativos basados en ecuaciones
empricas:
20

 Modelo de Cole-Cole. La expresin de la permitividad que propusieron K. S. Cole y

R. H. Cole en 1949, es la siguiente:

() = +

1 + ( j )1

(2-33)

donde es un parmetro de ajuste real cuyos valores estn comprendidos

0 < 1 . Su valor es inversamente proporcional a la temperatura y directamente

proporcional al nmero de grados de libertad internos de la molcula.

En el lmite = 0 , recuperamos la ecuacin de Debye. Este modelo es
apropiado para materiales orgnicos polares formados por cadenas largas. La
relajacin es ms lenta que en el modelo de Debye y el pico de prdidas tiene mayor
anchura y menor altura.

Figura 2.5. Relajacin en el modelo de Cole-Cole.

Modelo de Cole-Davidson. R. H. Cole y D. V. Davidson propusieron en 1951 una

frmula emprica que representara los fenmenos de relajacin asimtricos, como en

el glicerol. Esta es la expresin que propusieron para la permitividad compleja en
funcin de la frecuencia:

21

 () = +

(1 + j )

(2-34)

Donde 0 < 1 . Muchos medios que a temperatura ambiente sufren fenmenos de

relajacin que pueden ser descritos con el modelo de Debye ( = 1 ), se alejan de tal
modelo conforme la temperatura disminuye dado que la asimetra de su relajacin
aumenta, y por tanto el valor de disminuye.

Figura 2.6. Relajacin en el modelo de Cole-Davidson.

2.4.- Leyes de mezclas.

Dada una mezcla dielctrica de dos componentes 1 y 2, cuyas fracciones en
volumen son 1 y 2 = 1 1 , y permitividades estticas relativas s1 y s 2 , podemos
hacer un modelo simple que permita estimar el valor de la permitividad esttica relativa
de la mezcla s .
Un posible modelo consiste en considerar la mezcla compuesta por dos fases
bien diferenciadas a modo de capas. Si colocamos este medio compuesto entre las
placas de un condensador en alguna de estas dos disposiciones,

22

Figura 2.7. Disposiciones perpendicular y paralela del modelo de mezcla entre las placas
de un condensador de placas planoparalelas.

y calculamos la capacidad del condensador, podremos calcular la permitividad esttica

relativa de la mezcla para las dos disposiciones paralela y perpendicular.
Para disposicin paralela, tenemos dos condensadores en serie, de modo que la
capacidad resultante es:
1
1
1
=
+
C C1 C 2

d d
d
= 1 + 2
S s S s1 S s 2

1
= 1 + 2
s s1 s 2

(2-35)

Para disposicin perpendicular, tenemos dos condensadores en paralelo. En este

caso:
C = C1 + C 2

S s 1S s1 2S s 2
=
+
d
d
d

23

 s = 1 s1 + 2 s 2

(2-36)

La ecuacin (2-36) nos proporciona un lmite superior para el valor de s .A este

modelo se le llama de Wiener superior (W.U.). Se define constante dielctrica de exceso
para el modelo W.U. como:

( s )W.U. = ( s )exp ( s )W.U.

(2-37)

La ecuacin (2-35) proporciona un valor lmite inferior para la permitividad

esttica. Este es el modelo de Wiener inferior (W.L.) Se define constante dielctrica de
exceso para el modelo W.L. como:

( s )W.L. = ( s )exp ( s )W.L.

24

(2-38)

3.- Teora de las lneas de transmisin.

Una lnea de transmisin es un dispositivo, constituido por uno varios
conductores, que confina la energa electromagntica, guiando la propagacin de las
ondas electromagnticas. El trmino gua suele reservarse para las lneas de
transmisin constituidas por un solo conductor.

3.1.- Modelo circuital de una lnea de transmisin.

3.1.1.- Ecuaciones del telegrafista.

Si las dimensiones del circuito son muy pequeas en comparacin con la
longitud de onda, la intensidad elctrica en un instante dado puede considerarse
constante en cualquier punto de un cable conductor del circuito.
A continuacin se deducirn las ecuaciones diferenciales que deben satisfacer
los valores del voltaje y corriente en una lnea de transmisin. Asumimos que por la
lnea se propaga una seal electromagntica en el modo TEM puro. Con esto
garantizamos la unicidad en las definiciones de V e I. En caso de propagacin de otros
modos, habra que utilizar un convenio que las definiese.
Dividiendo la lnea en secciones de longitud z << podremos aplicar las leyes
de Kirchhoff (teora de circuitos de baja frecuencia); En la figura 3.1. se muestra el
modelo equivalente de una seccin de lnea de longitud z .

Figura 3.1. Modelo circuital de una seccin de la lnea de transmisin de longitud z.

25

Donde R, G, L y C son parmetros definidos por unidad de longitud. Rz da cuenta de

las prdidas en los conductores, Gz de las prdidas dielctricas, Lz representa el
almacenamiento de energa magntica y Cz el de energa elctrica.
Aplicamos las leyes de Kirchhoff a este modelo incremental y obtenemos:
i(z, t )

+ v(z + z, t )
v(z, t ) = Rz i(z, t ) + Lz

v(z + z, t )
i(z, t ) = Gz v(z + z, t ) + Cz
+ i(z + z, t )

(3-1)

Dividiendo ambas ecuaciones por z , y haciendo el lmite para z 0 ;

obtenemos las llamadas ecuaciones del telegrafista:
i(z, t )
v(z, t )
z = R i(z, t ) L t
i(z, t )
v(z, t )

= G v(z, t ) C
z
t

(3-2)

Para seales armnicas en el tiempo:

i(z, t ) = I(z ) exp( jt )
v(z, t ) = V(z ) exp( jt )

(3-3)

entonces las ecuaciones del telegrafista (3-2) quedan independientes de la variable

temporal:
V(z )
z = (R + jL )I(z )
I(z )

= (G + jC ) V(z )
z

(3-4)

Derivando ambas expresiones con respecto a z y sustituyendo trminos se

obtienen las ecuaciones desacopladas de Helmholtz para V e I:

26

 2 V(z )
2 V(z ) = 0

2
z
2
I(z ) 2 I(z ) = 0
z2

(3-5)

donde es la constante de propagacin:

= + j =

(R + jL )(G + jC)

(3-6)

En ausencia de prdidas R = G = 0; la constante de propagacin se reduce a:

= j = j LC ;

(3-7)

y la velocidad de propagacin de la seal es:

v=

(3-8)

LC

Las soluciones de las ecuaciones de Helmholtz se pueden expresar como

superposicin de una onda incidente y una reflejada:

V(z ) = V0+ e
I(z ) =

I 0+

- z

- z

+ V0 e

I 0

+ z

+ z

(3-9)

Imponiendo la condicin (3-4), se obtiene:

I(z ) =

V + e - z V e + z ;
0

R + jL

(3-10)

Las lneas de transmisin se caracterizan habitualmente por un parmetro con

dimensiones de resistencia que se denomina impedancia caracterstica, y se define como
el cociente entre la tensin y la corriente en una lnea de longitud infinita, y por tanto no
hay onda reflejada:

27

ZC =

V(z )
=
I( z ) V = 0
0

R + j L
G + jC

(3-11)

Cuando la lnea no tiene prdidas, la impedancia caracterstica es un nmero

real:
ZC =

L
C

(3-12)

As que podemos escribir la intensidad como:

I(z ) =

1 + - z
+ z
V0 e
V0 e

ZC

(3-13)

3.1.2.- Lnea de transmisin finita.

Dada una lnea de transmisin de longitud finita L, conectada a un generador
con carga adaptada y cargada al final de la lnea con una impedancia de valor ZL.

Figura 3.2. Lnea de transmisin de longitud finita.

Los valores de V e I en z = L son:

V(L ) = V0+ e
I(L ) =

I 0+

-L

-L

28

+ V0 e
I 0

+ L

+ L

(3-14)

Adems se cumple:
V(L ) = Z L I(L ); V0+ = Z C I 0+ ; V0 = Z C I 0 ;

(3-15)

se despeja V0+ y V0 :
V0+

V0 =

e + L
e

2
L
2

(V(L ) + I(L ) Z C )

(3-16)

(V(L ) I(L ) Z C )

y sustituyendo (3-16) en (3-9):

I(L )
{(Z L + Z C ) exp[(L z )] + (Z L Z C ) exp[ (L z )]}
2
I(L )
{(Z L + Z C ) exp[(L z )] (Z L Z C ) exp[ (L z )]}
I(z ) =
2Z C
V(z ) =

(3-17)

El coeficiente de reflexin se define como el cociente entre la onda reflejada y la

incidente:
(z ) =

V (z )
V + (z )

(3-18)

Por tanto, el coeficiente de reflexin en la carga es:

L =

V (L ) Z L Z C
=
V + (L ) Z L + Z C

29

(3-19)

3. 2.- Modelo electromagntico de una lnea de transmisin.

Para la descripcin de un fenmeno electromagntico, como en nuestro caso, se
utilizarn las ecuaciones de Maxwell, que en un medio dielctrico l.h.i. (lineal,
homogneo e istropo) son las siguientes:
r
E = 0

(3-20)

r
H = 0

(3-21)

r
r
H
E =
t

(3-22)

r
r
E
H =
t

(3-23)

donde y son la permitividad dielctrica y la permeabilidad magntica,

respectivamente.
Se deducir a partir de las ecuaciones de Maxwell la ecuacin de ondas. Para
ello se utiliza la siguiente propiedad vectorial:

) (

r v r
r v r v r
X = X 2X

(3-24)

r
donde X es una funcin vectorial. Aplicando esta propiedad al campo elctrico y

sustituyendo (3-20) se obtiene:

r v r
v r
E = 2 E

Ahora sustituimos (3-22) y (3-23) en (3-25), y se obtiene la ecuacin de

propagacin del campo elctrico:

30

(3-25)

r
v2r
2E
E 2 = 0
t

(3-26)

Y anlogamente para el campo magntico:

r
v2r
2H
H 2 = 0
t

(3-27)

Como cualquier funcin dependiente del tiempo se puede expresar como

superposicin de funciones armnicas, nos centraremos en suponer una solucin
armnica de la ecuacin de ondas:
r r
E = E 0 e z e jt

(3-28)

Sustituyendo esta solucin en la ecuacin de ondas (3-26), se obtiene la ecuacin

de Helmholtz:
r
r
2 E + 2 E = 0

2 + 2 = 0

(3-29)

(3-30)

donde es la constante de propagacin y es un nmero complejo = + j . La

velocidad de fase de la onda es:
v=

( )

Re

(3-31)

Se define la impedancia intrnseca de la onda como:

=

E
H

(3-32)

Sustituyendo la ecuacin de Maxwell (3-22) primero, y (3-30) despus:

j
j

=
=

31

(3-33)

3.2.1.- Cable coaxial.

Un cable coaxial es una lnea de transmisin que consta de dos conductores
cilndricos coaxiales, de radios a y b. Estos conductores estn separados por un medio
dielctrico. Con esta geometra los campos electromagnticos estn totalmente
confinados en su interior. Esto se demuestra con las leyes de Gauss y Ampre para el
campo elctrico y el campo magntico respectivamente. Adems se admite como modo
de propagacin el modo transversal magntico (TEM).

A continuacin se deducir una expresin, que relacione la impedancia

caracterstica del cable coaxial con su geometra y con las magnitudes macroscpicas
del dielctrico que hay entre los conductores.

Supongamos que se propaga por el cable coaxial el modo TEM, de modo que los
campos elctrico y magntico slo tendrn componentes perpendiculares a la direccin
de propagacin:
r E
E = 0 e z e jt

(3-34)

r
E
H = 0 e z e jt

(3-35)

Figura 3.3. Vista de la seccin transversal de un cable coaxial.

32

A partir de estas expresiones se calculan el voltaje y la intensidad:

b

V=

I=

b
E d = E 0 ln e z e jt
a

H ds =

E 0 2 z jt
z jt
e e = 2E 0
e e

(3-36)

(3-37)

De la definicin (3-11) se deduce que la impedancia caracterstica para una lnea

coaxial es:
ZC =

b 1 b
ln =
ln
2 a 2 a

(3-38)

Si el medio entre los conductores es un material no magntico, = 0 , como

ocurre en la mayora de los materiales, se tiene:

ZC =

1 0 b Z0
ln =
2 a
r

(3-39)

donde Z0 es la impedancia caracterstica con el aire como dielctrico entre los

conductores. Esta expresin relaciona parmetros de la lnea de transmisin, como son
la impedancia caracterstica y los radios, con la magnitud fsica que caracteriza al
dielctrico, que es la constante dielctrica.

33

34

4.- Tcnica de medida.

Para determinar la permitividad dielctrica de las muestras a estudio,
utilizaremos como tcnica de medida la reflectometra en el dominio del tiempo (TDR)
por el mtodo de 1 reflexin. Consiste bsicamente en analizar la seal reflejada por la
muestra colocada en una lnea de transmisin coaxial. El generador de seales emite una
seal escaln de voltaje que se propaga por la lnea hasta alcanzar la muestra. Un
osciloscopio digital recoge la seal reflejada. Tambin se recoge la seal reflejada por
un cortocircuito situado en la misma posicin que la muestra.
A partir de las medidas, se obtiene el coeficiente de reflexin (t ) en el dominio
temporal (TDR). Utilizando el mtodo de la curva derivada desarrollado por Samulon,
se obtiene el coeficiente de reflexin en el dominio frecuencial () y la constante
dielctrica compleja () . Para el clculo numrico se ha utilizado la aplicacin
pcmtdr.exe y se ha desarrollado la aplicacin tdr.exe.

4.1.- Descripcin del montaje experimental

A continuacin se muestra un esquema del dispositivo de medida.

Figura 4.1. Esquema del montaje experimental.

35

Est constituido por un TEST SET donde estn el generador de seales y los
sensores, un osciloscopio digital y un ordenador que recoge las medidas a travs del bus
HP-IB y realiza los clculos para obtener la permitividad compleja.
El generador del TEST SET HP-54121A emite seales escaln de 200mV de
amplitud y con un tiempo de subida de 35 ps. De los cuatro canales del TEST SET,
usaremos el canal 1, que es a la vez emisor y detector, y es el ideal para hacer estudios
de reflectometra. El resto de los canales son slo detectores.
El osciloscopio digital HP-54120B tiene una anchura de banda de 12.5 GHz y
adquiere las medidas en tiempo equivalente. Para ello es preciso que exista una
repetitividad de una seal para que pueda ser medida. El osciloscopio recoge unos pocos
puntos de la seal por cada repeticin y se utiliza una base de tiempos equivalente. Esta
es una caracterstica tpica de los osciloscopios con un gran ancho de banda. A
continuacin se muestra un ejemplo esquemtico de medida en tiempo equivalente.

Figura 4.2. Esquema de medida en tiempo equivalente.

El problema ms grande que se tiene utilizando una base de tiempos equivalente

es el jitter, que es el nombre con el que se denomina al ruido en la base de tiempos. Para
poder reconstruir la seal, es necesario que el jitter sea inferior al intervalo de muestreo.
Para este osciloscopio, el intervalo de muestreo mnimo es de 0.2 ps.
36

Para disminuir el ruido en voltaje, el osciloscopio hace promedios de la seal.

Tpicamente, el nmero de promedios es 128.
La lnea coaxial APC-7 es de alta precisin. Los conductores interno y externo
tienen dimetros de 3.04mm y 7mm, respectivamente. Calculando a partir de la
expresin (2-33) la impedancia caracterstica de la lnea es Z C = 50 . Al final de sta,
se coloca otro segmento de lnea que se utiliza como clula de medida. Se han utilizado
dos clulas de medida de longitud diferente: La clula corta de longitud 10 cm y
longitud til 8.5 cm, y la clula larga de longitud 20 cm y longitud til 18.5 cm.

Figura 4.3. Clula de medida.

El cuerpo central es de acero inoxidable que hace el papel de conductor externo

y contiene a la muestra. El conductor central es de cobre y en los extremos posee dos
cavidades donde se enroscan los tornillos que sujetan las arandelas de centrado, que
mantienen la posicin del conductor central y actan como extensiones del conductor
central.

Un esquema de la clula de medida y sus dimensiones pueden verse en las

figuras 4.3. y 4.4.

37

Figura 4.4. Vista del cuerpo de la clula y del conductor central.

En el extremo que se conectar a la lnea coaxial, se coloca una arandela de

centrado y adems una arandela de tefln que evitar el escape de la muestra.
Recordemos que se medirn muestras lquidas. En el otro extremo de la clula se coloca
una arandela de metal para hacer el cortocircuito al final de la lnea.

En la figura 4.5. se muestra un esquema de la arandela de centrado que est

hecha por un material transparente en el rango de las microondas.

Figura 4.5. Vista de la arandela de centrado (HP).

38

Este es un esquema del montaje de la clula de medida:

Figura 4.6. Montaje de la clula.

Y a continuacin se muestra un esquema de la carga de la clula para el mtodo

de 1 reflexin, que es el que hemos utilizado en el laboratorio.

Figura 4.7. Carga de la clula para el mtodo de primera reflexin.

4. 2.- Funcin de transferencia del sistema de medida.

Para conocer la respuesta de nuestro sistema de medida para cualquier tipo de
excitacin, deberemos conocer la relacin entre la seal perturbada por el sistema s(t) y
la de excitacin e(t):
s(t ) = T[e(t )]

39

(4-1)

Como el sistema de medida es un sistema lineal, podemos descomponer la seal

de excitacin en suma de funciones delta desplazadas y pesadas:
+

e(t ) =

e(t)(t t)dt

(4-2)

A partir de la expresin (4-1):

s (t ) = T e(t )(t t )dt

(4-3)

Si el sistema es lineal, la respuesta a una suma de seales ser igual a la suma de

las respuestas a cada una de las seales:
+

s (t ) =

e(t)T [(t t)]dt

(4-4)

As que T [ (t t )] es la funcin respuesta a un impulso, que representaremos

con el smbolo h (t; t ) . Entonces la relacin entre la salida y la entrada se puede
expresar como:
+

s (t ) =

e(t)h(t; t)dt

(4-5)

Si el sistema adems de lineal es invariante temporal, la respuesta impulso depende slo

de la diferencia t-t.
h (t ; t ) = h (t t )

(4-6)

Sustituimos en (4-5):
+

s (t ) =

e(t)h(t t)dt

40

(4-7)

que es la convolucin de la seal de excitacin y la respuesta impulso:

s (t ) = (e h )(t )

(4-8)

Haciendo la transformada de Fourier de esta expresin, se obtiene:

S() = H() E()

(4-9)

donde H() es la funcin de transferencia del sistema. Esto puede escribirse de otra
manera. Si W () representa la seal medida y V() la seal real, conociendo la funcin
de transferencia del sistema, podemos determinar la seal real a partir de la seal
medida:
V() =

W ()
H()

(4-10)

Para conocer la funcin de transferencia, podemos usar una carga cuya respuesta
sea conocida y nos pueda servir de referencia. Este es el caso de colocar un cortocircuito
como carga al final de la lnea. En este caso la onda reflejada real con el cortocircuito
es:

() =
VCC

()
WCC
= V + ()
H()

(4-11)

donde V + es la onda incidente. Y cuando coloquemos una muestra al final de la lnea,

la onda reflejada real ser:

V () =

W ()
= ()V + ()
H()

(4-12)

Haciendo el cociente entre las ecuaciones (4-12) y (4-11), se obtiene este importante
resultado:
V ()
W ()
() = +
=
V ()
WCC ()

41

(4-13)

Midiendo la seal reflejada por una muestra y la obtenida por un cortocircuito,

se puede determinar el valor del factor de reflexin verdadero, sin necesidad de
conocer la respuesta del sistema explcitamente.

4. 3.- Tratamiento de las medidas.

Las medidas registradas por el osciloscopio son los puntos del muestreo de la
seal dentro de una ventana temporal escogida. El paso temporal toma valores del orden
del picosegundo, y el tamao de la ventana toma valores comprendidos entre 300 ps y
900 ps.
Si se realizara la transformada de Fourier discreta sobre estos valores, el
resultado sera la convolucin de la transformada de la seal y la transformada de la
ventana (que es una funcin sinc), y se producira un fenmeno conocido como aliasing.
A continuacin se muestra una figura donde se explica este efecto.

Figura 4.8. Efecto de aliasing.

Para evitar este problema, se utilizar el mtodo de la curva derivada

desarrollado por Samulon [4], que consiste en hacer la transformada de Fourier discreta
de la curva derivada de la seal.

42

La curva derivada se define como:

V(t k ) = V(t k ) V(t k 1 ) = V(kt ) V(t (k 1))

(4-14)

donde tk = kt (k = 0, 1,, N-1) y N es el nmero de muestreos. La curva derivada en

los extremos y fuera de la ventana es nula y por tanto es insensible al efecto de la
ventana.
La transformada de Fourier discreta de la curva derivada es:

1
V( n ) =
N

[V(t ) V(t
k

k 1

)]exp( jn t k )

(4-15)

k =1

donde la frecuencia n = 2nf y f es el paso frecuencial. La expresin (4-15) se

puede expresar como:

1
V( n ) =
N

k =1

exp( j n t )
V(t k ) exp( j n t k )
N

V (t

k 1

) exp( jn t k 1 )

(4-16)

k =1

Teniendo en cuenta que los valores V(t 0 ) = V(t N ) , entonces:

V( n ) = {1 exp( j n t )}V( n )

(4-17)

Esta es la relacin entre la transformada de Fourier discreta de una seal y la

transformada de Fourier de la curva derivada. As que sustituyendo (4-17) en (4-13) se
obtiene la siguiente expresin para el coeficiente de reflexin en el dominio de la
frecuencia:
W () 1 exp( jt CC )
W ()
() =
=

() 1 exp( jt m )
WCC
WCC ()

43

(4-18)

donde t m es el intervalo de muestreo de la seal reflejada por la muestra, y t CC es el

intervalo de muestreo de la seal reflejada por el cortocircuito. En nuestras medidas,
siempre los intervalos de muestreo son iguales. En este caso la expresin (4-18) se
simplifica:
W ()
() =

()
WCC

(4-19)

4. 4.- Clculo de la constante dielctrica.

Se necesita encontrar la expresin que relaciona la permitividad dielctrica y el
coeficiente de reflexin para medios no magnticos. Teniendo en cuenta las ecuaciones
(3-19) y (3-39), se tiene:
1
() =

r ()
1
r ()

1
(4-20)
+1

Se despeja la constante dielctrica compleja en funcin de la frecuencia:

1 ()
= ' () j ()
() =
(
)
1
+

(4-21)

Para hacer todos estos clculos se utiliza la aplicacin pcmtdr.exe [5], que
dispone el laboratorio. Tambin se ha desarrollado en este Trabajo Acadmicamente
Dirigido una aplicacin: tdr.exe; para el mtodo de primera reflexin. El cdigo en
lenguaje C del programa se adjunta en el Apndice de esta Memoria. Las novedades de
este programa, con respecto al pcmtdr.exe son: el clculo del tiempo de relajacin,
presenta los valores medidos normalizados. Permite calcular la permitividad real e
imaginaria con un paso en frecuencias menor que con el pcmtdr.exe. Para una ventana
en frecuencias dada, es capaz de medir * para mil valores de la frecuencia. Por
ejemplo, para una ventana de 10 GHz, el paso mnimo del es de 0.01 GHz, mientras que
para el pcmtdr.exe es de 0.05GHz.
44

4. 5.- Correccin de la arandela.

En todos los clculos que se han realizado anteriormente, no se ha tenido en
cuenta el hecho de que se coloca una arandela de tefln para evitar la fuga de las
muestras lquidas que se van a medir. As que en vez de estar midiendo el factor de
reflexin que se produce en la frontera aire-muestra, lo que medimos en realidad es el
coeficiente de reflexin en las fronteras aire-tefln y tefln-muestra.
La correccin la realizamos en el dominio de frecuencias [6]. Si llamamos Yen a
la admitancia medida para el conjunto tefln-muestra y Ym a la admitancia real de la
muestra:

Ym =

Yca Yen
Z cc Yen 1

(4-22)

Yca es la admitancia de entrada con la salida terminada en circuito abierto y Zcc es la

impedancia de entrada con la salida terminada en cortocircuito. Las relaciones de Yca y
Zcc con el espesor d de la arandela, con la velocidad de propagacin c de la seal en la
lnea y con la permitividad del tefln , se muestran a continuacin:

Yca = Y0 tanh j d / c

Z cc = Z 0

tanh j d / c

45

(4-23)

(4-24 )

46

5.- Resultados experimentales.

En este trabajo se han realizado medidas sobre mezclas dielctricas compuestas
por dos componentes: uno polar y otro no polar. El componente no polar es el
ciclohexano y como componentes polares el pentanol, el hexanol y el heptanol.
Como ilustracin se muestra en la figura 5.1. los voltajes normalizados, frente al
tiempo, de las seales reflejadas por una muestra de 1-pentanol + ciclohexano con una
concentracin de pentanol del 5%.

Voltajes normalizados de las seales reflejadas

por la muestra y por el cortocirc uito.
1.0

corto
muestra

0.8

V/V

MA X

0.6

0.4

0.2

0.0
0

50

100

150

200

250

300

350

400

t (ps)

Figura 5.1. Voltajes normalizados de las seales reflejadas por un cortocircuito y la

muestra frente al tiempo.

Mediante las herramientas de clculo (pcmtdr.exe y tdr.exe), se calculan los

valores de ' y " en funcin de la frecuencia, los valores de la frecuencia y el tiempo
de relajacin, y "max . En las figuras 5.2., 5.3., y 5.4. se muestran las curvas de
dispersin obtenidas para las tres muestras.

47

1-pentanol+ciclohexano
2.50

0.5

2.00

0.4

' (5%)
' (10%)

1.50

0.3

' (15%)

'

''

'' (5%)
1.00

0.2

'' (10%)
'' (15%)

0.50

0.1

0.00

Frec(GHz)

Figura 5.2. Relaciones de dispersin del 1-pentanol + ciclohexano para las tres
concentraciones del componente polar.

48

1-hexanol+ciclohexano
2.5

2.0

0.8
' ( 5%)
' (10%)

1.5

0.6

' (15%)

''

'
'' (5%)

1.0

0.4

'' (10%)
'' (15%)
0.5

0.2

0.0

0
0

Frec (GHz)

Figura 5.3. Relaciones de dispersin del 1-hexanol + ciclohexano para las tres
concentraciones del componente polar.

49

1-heptanol+ciclohexano
2.5

0.5

2.0

0.4
' ( 5%)
' (10%)

1.5

0.3

' (15%)

'

''
0.2

1.0
'' ( 5%)
'' (10%)
0.5

0.1

'' (15%)
0

0.0
0

Frec (GHz)

Figura 5.4. Relaciones de dispersin del 1-heptanol + ciclohexano para las tres
concentraciones del componente polar.

En la tabla 5.1. se muestran los valores de la permitividad esttica, la frecuencia

de relajacin, el tiempo de relajacin, y el valor mximo de la parte imaginaria de la
permitividad

dielctrica;

obtenidos

para

el

1-pentanol+ciclohexano,

1-

hexanol+ciclohexano y 1-heptanol+ciclohexano. Las concentraciones del componente

50

polar son del 5%, 10% y 15%. Todas las medidas se realizaron a temperatura constante
(T = 25C)

5%
PTC
HXC
HPC

10%
PTC
HXC
HPC

15%
PTC
HXC
HPC

S (ps) max
0.03 0.5 (0.01
2.15
2.14
2.12

55
58
67

0.05
0.02
0.018

S (ps) max
0.03 0.5 (0.01
2.25
2.23
2.23

69
71
85

0.07
0.06
0.03

S (ps) max
0.03 0.5 (0.01
2.37
2.35
2.34

85
90
110

0.12
0.10
0.07

Tabla 5.1. Valores de la constante dielctrica esttica, la

frecuencia y el tiempo de relajacin, y el valor mximo
EMBED Equation.3

para: pentanol+ciclohexano

(PTC), hexanol+ciclohexano (HXC) y heptanol+ciclohexano

(HPC).

5. 1.- Permitividad esttica.

En la figura 5.5. se muestran los valores de la constante dielctrica esttica en
funcin de la concentracin del componente polar. Los valores a mayor concentracin,
fueron medidos en estudios previos en el mismo laboratorio con la misma tcnica de
medida [7].

51

16
( ) PTC
s exp

14

( ) HXC

( ) HPC

s exp

s exp

Variacin de la constante dielctrica esttica

Mezclas 1-alcohol + ciclohexano

12
10
8
6
4
2

X-polar

0
0

0.2

0.4

0.6

0.8

Figura 5.5. Representacin de la constante dielctrica esttica en funcin de la

concentracin del componente polar: pentanol (PTC), hexanol (HXC), heptanol (HPC).

Se observa que los valores de la esttica es encuentran comprendidos entre las

estticas del componente no polar (ciclohexano) de valor s,ciclohexano = 2.016 y las de
los componentes polares puros. Se observa tambin que el valor de la esttica disminuye
si el nmero de carbonos del componente polar aumenta. Esto se debe a que cuantos
ms carbonos tiene, la parte no polar es ms amplia que la polar, que es el radical OH.
Calculando los valores de la constante dielctrica esttica con los modelos de
Wiener superior (W.U.) e inferior (W.L.) dados por las expresiones (2-36) y (2-35)

52

respectivamente, para estas concentraciones; y representando la constante dielctrica de

exceso para ambos modelos (2-37) y (2-38), se obtienen las siguientes grficas.

0.5
( )

s W.U.

PTC

( )

s W.U.

HXC

( )

s W.U.

HPC

Variacin de la constante dielctrica de exceso

Modelo Wiener-superior

-0.5
-1
-1.5
-2
-2.5
-3

X-polar

-3.5
0

0.2

0.4

0.6

0.8

Figura 5.6. Representacin de la constante dielctrica de exceso del modelo de Wiener

superior en funcin de la concentracin del componente polar.

53

5
( )

s W.L.

PTC

( )

s W.L.

HXC

( )

s W.L.

HPC

Constante dielctrica de exceso

Modelo Wiener inferior

X-polar

-1

-2
0

0.2

0.4

0.6

0.8

Figura 5.7. Representacin de la constante dielctrica de exceso del modelo de Wiener

inferior en funcin de la concentracin del componente polar.

5. 2.- Absorcin mxima.

En la siguiente grfica se representa " mxima frente a la concentracin de la
componente polar.

54

max

PTC

max

HXC

max

HPC

Variacin de la absorcin mxima

Mezcla 1-alcohol + ciclohexano

1
X-polar

0
0

0.2

Figura 5.8. Representacin de

0.4

0.6

0.8

en funcin de la concentracin del componente polar.

Se observa que la absorcin mxima es creciente con la concentracin. Y al

igual que para el valor de la esttica, " disminuye conforme el nmero de carbonos
aumenta.

5. 3.- El tiempo de relajacin.

En la siguiente grfica se representan los valores del tiempo de relajacin frente
a la concentracin del componente polar para las tres mezclas.

55

1200
PTC
0

1000

HXC

HPC

Variacin del tiempo de relajacin

Mezclas 1-alcohol + ciclohexano

800

600

400

200
X-polar

0
0

0.2

0.4

0.6

0.8

Figura 5.9. Representacin del tiempo de relajacin en funcin de la concentracin del

componente polar para las tres mezclas.

Se observa que el valor del tiempo de relajacin decrece conforme la

concentracin del componente polar es menor. Representando esta misma grfica en
escala logartmica se puede apreciar con mejor detalle su comportamiento:
PTC
0

1000

HXC
0

ln( )-x

polar

HPC
0

100
X
0

0.2

0.4

0.6

0.8

Figura 5.10. Representacin del tiempo de relajacin en funcin de la concentracin del

componente polar para las tres mezclas en escala logartmica.

56

La grfica muestra que existen dos regiones. La zona de transicin aparece para
una concentracin crtica x c 0.17 . Para baja concentracin, se observa una
dependencia lineal entre ln y la concentracin x:
ln = ln 0 + ax

(5-1)

Despejando , se obtiene una ley del tipo Arrhenius:

= 0 e ax

(5-2)

Adems esta dependencia, es prcticamente independiente del nmero de

carbonos del componente polar. Las molculas de ciclohexano inhiben la formacin de
multipolos. Por el contrario, para concentraciones superiores a xc, el tiempo de
relajacin depende mucho ms del nmero de carbonos del componente polar. El
tamao de los multipolos formados depende de la longitud de las cadenas del alcohol.
Si representamos para las 8 concentraciones, el tiempo de relajacin frente al
nmero de carbonos, se obtiene:
1200
Variacin del tiempo de relajacin
( concentracin constante)

(0.05)
(0.10)
(0.15)
(0.2)
(0.4)
(0.6)
(0.8)
(1.0)

1000

800

600

400

200
n C

0
4

Figura 5.11. Representacin del tiempo de relajacin en funcin del nmero de carbonos
para las 8 concentraciones.

57

Si se ajustan a rectas estos puntos para cada concentracin y = mx + n; y se

representa el valor de la pendiente m en funcin de la concentracin, se pone en
evidencia la existencia de dos regiones de comportamiento.

200
m

Variacin de la pendiente de las rectas

"tiempo de relajacin vs n de carbonos"
Mezclas 1-alcohol + ciclohexano

150

100

50

x-polar

0
0

0.2

0.4

0.6

0.8

Figura 5.12. Representacin el valor de la pendiente m en funcin de la concentracin.

58

59

6.- Conclusiones.

El objetivo de este trabajo ha sido el estudio del comportamiento electromagntico

de mezclas dielctricas binarias polar no polar. El componente no polar era el
ciclohexano y los componentes polares estudiados eran los alcoholes primarios
pentanol, hexanol y heptanol.

Se ha utilizado la tcnica de reflectometra en dominio temporal (TDR) por el

mtodo de primera reflexin. Esta tcnica permite la obtencin de la curva de
dispersin en el rango de DC-5GHz donde se presentan las particularidades de la
relajacin dielctrica en estas mezclas, en particular las magnitudes permitividad
esttica, absorcin mxima y tiempo de relajacin.

Se ha desarrollado un programa de clculo que permite obtener curvas de dispersin

con pasos frecuenciales de 0.005 GHz que mejora a otros existentes. Tambin
calcula la absorcin mxima, el tiempo de relajacin de la mezcla y las curvas de los
voltajes normalizados.

El valor de la constante dielctrica esttica disminuye si el nmero de carbonos del

componente polar aumenta debido a que las molculas polares estn constituidas por
una estructura apolar y un radical OH. Conforme el nmero de carbonos aumenta,
la parte apolar crece dificultando la formacin de los multipolos.

A bajas concentraciones del componente polar, la constante dielctrica esttica se

mantiene prcticamente constante y resulta independiente del nmero de carbonos.
Esto se debe a que las molculas de ciclohexano inhiben la formacin de multipolos,
aislando los dipolos de las molculas polares.

Los valores de la constante dielctrica esttica estn comprendidos entre los

estimados por los modelos de Wiener-inferior y Wiener-superior excepto para la

60

mezcla pentanol+ciclohexano donde los valores de la permitividad esttica, toman

valores inferiores a los estimados por dichos modelos. Esto puede deberse al hecho
de que el nmero de carbonos del pentanol sea inferior al del componente no polar,
y esto le permita al ciclohexano, inhibir con mayor facilidad la contribucin a la
permitividad del pentanol.

La absorcin mxima aumenta con la concentracin para las tres mezclas.

El valor del tiempo de relajacin decrece conforme la concentracin del componente

polar disminuye. A menor concentracin el tamao de los multipolos formados son
menores, y por tanto las frecuencias de rotacin sern superiores, que se traduce en
tiempos de relajacin ms cortos.

La variacin del tiempo de relajacin respecto a la concentracin del componente

polar, presenta dos regiones de comportamiento. La zona de transicin est en torno
a x c 0.17 . En la regin de baja concentracin el tiempo de relajacin es
prcticamente independiente del nmero de carbonos, al igual que ocurre con los
valores obtenidos de la permitividad esttica, y los puntos experimentales se pueden
ajustar a una ecuacin tipo Arrhenius.

61

Bibliografa

[1] J.D. Jackson, Classical Electrodynamics, John Wiley & Sons, (1999)

[2] J.M. Albella Martn, J. M. Martnez Duart, Fsica de dielctricos, Ed. Marcombo,
Barcelona (1984)

[3] R.K. Wangsness. Campos electromagnticos, Ed. Limusa, (1996)

[4] H.A. Samulon, Spectrum analysis of transient response curves Proc. IRE, (1951)

[5] J. Cabeza Guilln Puesta a punto de diversos mtodos de caracterizacin de medios

dielctricos y magnticos en microondas, Tesina Universidad de Zaragoza, (2000)

[6] J. Letosa, Caracterizacin de dielctricos por TDR: anlisis de errores aleatorios,

Tesis Universidad de Zaragoza, (1997)

[7] I. A. Ghemes, Estudio de la relajacin dielctrica en mezclas polar-no polar por

TDR. Trabajo estancia Erasmus- Universidad de Zaragoza. (2003).

62

63

Apndice: Cdigo en lenguaje C de la aplicacin tdr.exe

/* Autor: DANIEL DIGON RODRIGUEZ 2006 */
/* tdr.exe*/
/* Reflectometra en el dominio del tiempo. Mtodo de 1 reflexin*/
#include <math.h>
#include <stdio.h>
#include <conio.h>
#include <unistd.h>
#include <stdlib.h>
#define VGA
#include "graficos.h"
#define N 1000
int mue[4], cor[4], i, leido=0;
int dm[N], dc[N],max,min;
double datmue[7], datcor[7];
float vm[N], vc[N],maxdos, mindos;
double t[N], f[N];
float pasof, venf,tamf;/*pasof y venf en GHz*/
float rperm[N],iperm[N], rer[N],imr[N],rtfcdc[N], rtfcdm[N],itfcdc[N],itfcdm[N];
float ipmax, fmax, tr, ipmin,rpmax, rpmin, maxtot,mintot;
char nmue[14], ncor[14], name[14], fr[8]="ini",fw[8];
char inmue,incor, menu, caso;
FILE *esc;
FILE *lemue, *lecor;
void marco(void);
void tfd(void);
void men(void);
void info(void);
int guarda(void);
void dibuja(void);
int leo(void);
int main()
{
inicia_graf(MODO_VGA);
inicio:
men();
eli:
menu=getch();
if(menu=='0'){inicia_graf(MODO_TXT);return(0);}
if(menu!='1'&&menu!='2'&&menu!='3'&&menu!='4'&&menu!='5'){goto eli;}
if(menu=='1'){info();getch();goto inicio;}
if(menu=='2'){leido=leo();goto inicio;}
if(leido==0)goto eli;
if(menu=='3'){
max=dc[0];
min=dc[0];
for(i=0;i<mue[2];i++)
{
if(max<dc[i]){max=dc[i];}
if(min>dc[i]){min=dc[i];}
}
for(i=0;i<mue[2];i++){
vc[i]=((float)dc[i]-(float)min)/(float)(max-min);
vm[i]=((float)dm[i]-(float)min)/(float)(max-min);

64

}}
if(menu=='4'||menu=='5'){
for(i=0;i<mue[2];i++)
{
vc[i]=((float)dc[i]-(float)datcor[5])*(float)datcor[3]+(float)datcor[4];
vm[i]=((float)dm[i]-(float)datmue[5])*(float)datmue[3]+(float)datmue[4];
}
}
if(menu=='4'){
maxdos=vc[0];
mindos=vc[0];
for(i=0;i<mue[2];i++)
{
if(maxdos<vc[i]){maxdos=vc[i];}
if(mindos>vc[i]){mindos=vc[i];}
}}
if(menu=='5')tfd();
fclose(lemue);
fclose(lecor);
optar:
system("CLS");
marco();
gotoxy(20,12);
printf("1) Volver al menu");
gotoxy(20,14);
printf("2) Ver grfica");
gotoxy(20,16);
printf("3) Guardar clculos");
gotoxy(20,18);
printf("0) Salir del programa");
recoge:
caso=getch();
if(caso=='1'){goto inicio;}
if(caso=='3'){guarda();goto optar;}
if(caso=='2'){dibuja();goto optar;}
if(caso=='0'){
inicia_graf(MODO_TXT);
return(0);}
goto recoge;
}
void tfd(void)
{
float esa=3/*espesor arandela en mm*/, pea=2.04;/*permitividad relativa de la arandela*/
float rpa[N],ipa[N];/*permitividades real e imaginarias corregidas de arandela*/
float ryen[N],iyen[N], yca[N], zcc[N];
float reyd[N],imyd[N];/*admitancia corregida*/
char aran;
int j, k;
float deltat=datmue[0]*1000000000.0;
initf:
system("CLS");
marco();
gotoxy(12,4);printf("Introduzca el paso de frecuencia en GHz (ej 0.1)\n");
gotoxy(12,5);scanf("%f",&pasof);
if(pasof==0){gotoxy(12,10);printf("El paso de frecuencia ha de ser mayor que 0");goto initf;}
gotoxy(12,7);printf("Introduzca la ventana deseada en GHz (ej 10)\n");
gotoxy(12,8);scanf("%f",&venf);
tamf=venf/pasof;
gotoxy(12,10);printf("%f puntos",tamf);
if(tamf>N+0.01){gotoxy(18,14);printf(">1000 Demasiados puntos!\n");getch();goto initf;}

65

pregunta:
gotoxy(12,12);printf("Desea correccin de arandela? (s/n)");aran=getch();
if(aran!='s'&& aran!='n') goto pregunta;
if(aran=='s')
{
gotoxy(12,14);printf("Introduzca espesor de la arandela en mm: ");scanf("%f",&esa);
gotoxy(12,16);printf("Introduzca permitividad de la arandela: ");scanf("%f",&pea);
}
for(k=0;k<mue[2];k++)
{
t[k]=deltat*k; /*en nanosegundos*/
}
for(k=0;k<tamf;k++)
{
f[k]=k*pasof+0.01;
}
for(j=0;j<tamf;j++)
{
rtfcdc[j]=0; /*inicializacion*/
itfcdc[j]=0;
rtfcdm[j]=0;
itfcdm[j]=0;
for(k=0;k<mue[2]-1;k++)
{
rtfcdc[j]=rtfcdc[j]+(float)(1.0/((float)mue[2]))*((float)vc[k+1](float)vc[k])*((float)cos(2.0*PI*f[j]*t[k]));
itfcdc[j]=itfcdc[j]-(float)(1.0/((float)mue[2]))*((float)vc[k+1](float)vc[k])*((float)sin(2.0*PI*f[j]*t[k]));
rtfcdm[j]=rtfcdm[j]+(float)(1.0/((float)mue[2]))*((float)vm[k+1](float)vm[k])*((float)cos(2.0*PI*f[j]*t[k]));
itfcdm[j]=itfcdm[j]-(float)(1.0/((float)mue[2]))*((float)vm[k+1](float)vm[k])*((float)sin(2.0*PI*f[j]*t[k]));
}
}
for(j=0;j<tamf;j++)
{
rer[j]=-(rtfcdm[j]*rtfcdc[j]+itfcdm[j]*itfcdc[j])/(rtfcdc[j]*rtfcdc[j]+itfcdc[j]*itfcdc[j]);
imr[j]= (rtfcdm[j]*itfcdc[j]-itfcdm[j]*rtfcdc[j])/(rtfcdc[j]*rtfcdc[j]+itfcdc[j]*itfcdc[j]);
if(aran=='n')
{
rperm[j]=((1-rer[j]*rer[j]-imr[j]*imr[j])*(1-rer[j]*rer[j]-imr[j]*imr[j])4.0*imr[j]*imr[j])/(((1+rer[j])*(1+rer[j])+imr[j]*imr[j])*((1+rer[j])*(1+rer[j])+imr[j]*imr[j]));
iperm[j]=4*imr[j]*(1-rer[j]*rer[j]imr[j]*imr[j])/(((1+rer[j])*(1+rer[j])+imr[j]*imr[j])*((1+rer[j])*(1+rer[j])+imr[j]*imr[j]));
}
}
if(aran=='s')
{
for(j=0;j<tamf;j++)
{
ryen[j]=(1.0/50.0)*(1-rer[j]*rer[j]-imr[j]*imr[j])/((1+rer[j])*(1+rer[j])+imr[j]*imr[j]);
iyen[j]=-(2.0/50.0)*imr[j]/((1+rer[j])*(1+rer[j])+imr[j]*imr[j]);
zcc[j]=50*(float)tan((2.0/299.79)*PI*(double)f[j]*(double)esa*sqrt((double)pea))/(float)sqrt((doubl
e)pea); /*imaginario puro, solo pongo el valor absoluto*/
yca[j]=(float)sqrt((double)pea)*(float)tan((2.0/299.79)*PI*(double)f[j]*(double)esa*sqrt((double)pe
a))/50.0; /*tb imaginario puro*/

66

reyd[j]=ryen[j]*(1+zcc[j]*yca[j])/((float)pow((double)(1+zcc[j]*iyen[j]),2)+(float)pow((double)(zcc[j]*
ryen[j]),2));
imyd[j]=(ryen[j]*ryen[j]*zcc[j]+(iyen[j]yca[j])*(1+zcc[j]*iyen[j]))/((float)pow((double)(1+zcc[j]*iyen[j]),2)+(float)pow((double)(zcc[j]*ryen[j]
),2));
rperm[j]=(50*50)*(reyd[j]*reyd[j]-imyd[j]*imyd[j]);
iperm[j]=-2*50*50*reyd[j]*imyd[j];
}
}
fmax=f[0];
ipmax=iperm[0];
ipmin=iperm[0];
rpmax=rperm[0];
rpmin=rperm[0];
for(j=0;j<tamf;j++)
{
if(ipmax<iperm[j]){ipmax=iperm[j];fmax=f[j];}
if(ipmin>iperm[j]){ipmin=iperm[j];}
if(rpmax<rperm[j]){rpmax=rperm[j];}
if(rpmin>rperm[j]){rpmin=rperm[j];}
}
if(rpmax>ipmax)maxtot=rpmax;else maxtot=ipmax;
if(rpmin<ipmin)mintot=rpmin;else mintot=ipmin;
tr=1000.0/(2.0*PI*fmax);
gotoxy(12,18);printf("M ximo de e\"");
gotoxy(14,21);printf("fmax=%f GHz",fmax);
gotoxy(14,20);printf("e\"max=%f",ipmax);
gotoxy(14,22);printf("tau=%f ps",tr);
getch();
}
void men(void){
system("CLS");
marco();
gotoxy(10,3);
printf("Reflectometra en el dominio del tiempo (T.D.R.)");
gotoxy(18,6);
printf("ELIJA UNA OPCION");
if(leido==1){
gotoxy(10,12);
printf("3) Crear fichero con medidas de MUESTRA y CORTO normalizadas");
gotoxy(10,14);
printf("4) Crear fichero con voltajes de MUESTRA y CORTOCIRCUITO");
gotoxy(10,16);
printf("5) Clculo de la permitividad dielctrica en funcin de");
gotoxy(10,17);
printf(" la frecuencia por el mtodo de primera reflexin");
}
gotoxy(10,8);
printf("1) Informacin sobre el programa");
gotoxy(10,10);
printf("2) Abrir archivo");
gotoxy(10,19);
printf("0) Salir");
gotoxy(10,22); printf("AREA DE ELECTROMAGNETISMO");
gotoxy(10,23); printf("DEPARTAMENTO DE FISICA APLICADA");
gotoxy(10,24); printf("UNIVERSIDAD DE ZARAGOZA CURSO 2005/06");
}
void marco(void)

67

{
rect_lleno(X_M-6,0,X_M-1,Y_M-1,1);
rect_lleno(0,0,6,Y_M-1,1);
rect_lleno(0,0,X_M-1,5,1);
rect_lleno(0,Y_M-6,X_M-1,Y_M-1,1);
}
void info(void)
{
system("CLS");
marco();
gotoxy(18,3);printf("INFORMACION");
gotoxy(8,5);printf("TITULO DEL T.A.D.: \"Estudio de la relajacin dielctrica en");
gotoxy(8,6);printf("
mezclas polar - no polar por T.D.R.\"");
gotoxy(8,8);printf("TUTOR: Juan Pablo Martnez Jimnez");
gotoxy(8,10);printf("ALUMNO: Daniel Dign Rodrguez");
gotoxy(8,12);printf("CURSO 2005/06");
gotoxy(8,15);printf("El programa permite leer los ficheros .mue y .cor donde estn las ");
gotoxy(8,16);printf("medidas, y escribe en un fichero .d las siguientes opciones:");
gotoxy(8,19);printf("-Los valores normalizados del corto y la muestra en funcin del tiempo");
gotoxy(8,21);printf("-Los valores en voltios del corto y la muestra en funcin del tiempo");
gotoxy(8,23);printf("-La permitividad elctrica en funcin de la frecuencia por el mtodo");
gotoxy(9,24);printf("de primera reflexin con correccin de arandela. Para esta opcin, el");
gotoxy(9,25);printf("programa nos da el valor de la frecuencia correspondiente al mximo de");
gotoxy(9,26);printf("la parte imaginaria de la permitividad");
}
int guarda(void)
{
int i;
system("CLS");
marco();
gotoxy(8,4);printf("Los archivos de lectura son %s.mue y %s.cor", fr,fr);
gotoxy(8,6);printf("Escriba el nombre del archivo de escritura");
gotoxy(8,7);scanf("%s",fw);
gotoxy(8,9);printf("El archivo de escritura es %s.d", fw);
sprintf(name,"a:%s.d",fw);
esc=fopen(name,"wt");
if(esc==NULL){gotoxy(8,11);printf("ERROR: No hay disquete en a:");getch();return 0;}
if(menu=='3')
{
fprintf(esc,"t corto muestra\n");
}
if(menu=='4'){fprintf(esc,"t cor(V) mue(V)\n");}
if(menu=='3'||menu=='4')
{
for(i=0;i<mue[2];i++)
{
fprintf(esc,"%d %f %f\n",i, vc[i],vm[i]);
}
}
if(menu=='5')
{
fprintf(esc,"f(GHz) e' e\"\n");
for(i=0;i<tamf;i++){
fprintf(esc,"%.2f %f %f\n",f[i], rperm[i], iperm[i]);
}
}
fclose(esc);
if(esc!=NULL){gotoxy(10,12);printf("Clculos guardados!");}
getch();
return 1;
}

68

void dibuja(void)
{
int i,j=0;
int pas=(int)(mue[2]/512),poscy,posmy;
int posrpx, posrpy, posipx, posipy;
if (pas==0){pas=1;}
system("CLS");
marco();
rectangulo(64,48,576,432,BLANCO);
if(menu=='3')
{
gotoxy(6,4);printf("1");
gotoxy(6,28);printf("0");
gotoxy(6,2);printf("Valores normalizados del corto (AMARILLO) y la muestra (ROJO)");
}
if(menu=='4')
{
gotoxy(3,4);printf("%.3f",maxdos);
gotoxy(3,28);printf("%.3f",mindos);
gotoxy(6,2);printf("Voltajes (en voltios) del corto (AMARILLO) y la muestra (ROJO)");
}
if(menu=='3'||menu=='4')
{
gotoxy(72,29);printf("%d",mue[2]);
for(i=0;i<mue[2];i++)
{
if(menu=='3'){
poscy= (int)((1-vc[i])*384);
posmy= (int)((1-vm[i])*384);
}
if(menu=='4'){
poscy= (int)((maxdos-vc[i])*384/(maxdos-mindos));
posmy= (int)((maxdos-vm[i])*384/(maxdos-mindos));
}
elipse_f((int)(64+(float)i*512.0/((float)mue[2])),poscy+48,2,2,AMARILLO);
elipse_f((int)(64+(float)i*512.0/((float)mue[2])),posmy+48,2,2,ROJO);
}
}
if(menu=='5')
{
gotoxy(3,4);printf("%.3f",maxtot);
gotoxy(3,28);printf("%.3f",mintot);
gotoxy(10,2);printf("Parte real (AMARILLO) y parte imaginaria (ROJO) de");
gotoxy(10,3);printf("la permitividad frente a la frecuencia en GHz");
gotoxy(9,29);printf("0");
gotoxy(72,29);printf("%.1f",venf);
/*linea de cero*/
if(maxtot>0 && mintot<0){
gotoxy(4,(int)(4+24*(maxtot)/(maxtot-mintot)));printf("0");
linea_fast(64,(int)(48+384*(maxtot)/(maxtot-mintot)),576,(int)(48+384*(maxtot)/(maxtotmintot)),BLANCO);
}
for(i=0;i<tamf;i++)
{
elipse_f((int)(64+(f[i]-0.01)*512.0/venf),(int)(432-384*(rperm[i]-mintot)/(maxtotmintot)),2,2,AMARILLO);
elipse_f((int)(64+(f[i]-0.01)*512.0/venf),(int)(432-384*(iperm[i]-mintot)/(maxtotmintot)),2,2,ROJO);
}

69

}
getch();
}
int leo(void)
{
system("CLS");
marco();
gotoxy(8,4);printf("Escriba el nombre del archivo de lectura");
gotoxy(8,6);scanf("%s",fr);
gotoxy(8,8);printf("Los archivos de lectura son %s.mue y %s.cor", fr,fr);
sprintf(nmue,"a:%s.mue",fr);
sprintf(ncor,"a:%s.cor",fr);
lemue=fopen(nmue,"r");
if(lemue == NULL){gotoxy(8,12);printf("El archivo %s.mue no se encuentra en el disquete",fr);
getch();}
lecor=fopen(ncor,"r");
if(lecor == NULL){gotoxy(8,14);printf("El archivo %s.cor no se encuentra en el disquete",fr);
getch();}
if(lecor==NULL&&lemue==NULL){gotoxy(8,16);printf("Es posible que no haya disquete en
a:");
getch();}
if(lecor==NULL||lemue==NULL)return 0;
rewind(lemue);
rewind(lecor);
for(i=0;i<4;i++){
fscanf(lemue,"%d ",&mue[i]);
fscanf(lecor,"%d ",&cor[i]);
}
if(mue[2]>N)
{ gotoxy(8,18);printf("Archivo demasiado grande\n");
getch();
return 0;
}
for(i=0;i<7;i++){
fscanf(lemue,"%E ",&datmue[i]);
fscanf(lecor,"%E ",&datcor[i]);
}
for(i=0;i<mue[2];i++){
fscanf(lemue,"%c",&inmue);
fscanf(lemue,"%d",&dm[i]);
fscanf(lemue,"%c",&inmue);
fscanf(lecor,"%c",&incor);
fscanf(lecor,"%d",&dc[i]);
fscanf(lecor,"%c",&incor);
}
return 1;
}

70