I. INTRODUCCIÓN
En el Laboratorio de Dieléctricos de la FACET se estudia, desde hace tiempo, la caracterización de propiedades y estados de madurez del tomate a partir de mediciones dieléctricas. Se trabaja con muestras de tomates en sus distintos estados: líquido, pasta, trozo y fruta entera usando una celda Hewlett-Packard 16451B de electrodos planos y paralelos con anillo de guarda para las primeras muestras y una celda de construcción propia para la fruta entera.
En este último caso, para el estudio de la fruta entera, se buscó un método de trabajo no tradicional que permita realizar las mediciones sin dañar el tomate. Se reprodujo la celda propuesta por Varlan-Sansen1 y luego se modificaron los conectores para minimizar los errores de mediciones eléctricas por cables.
En la interpretación de los resultados de las mediciones obtenidas, cuando se pretende contrastar con teorías de comportamiento dieléctrico de medios materiales, aparece la necesidad de conocer cualitativamente la forma de los campos generados por los electrodos de la celda, ya que se supone la interacción con campos generalmente uniformes.
La búsqueda en la bibliografía específica permitió verificar que la situación de interés no estaba resuelta. Por ejemplo, Pallares Muñoz y Rodríguez Calderón2 diseñan y desarrollan una herramienta educativa para estudiar el efecto de bordes en un condensador de placas paralelas usando la ecuación de Laplace en 2D con diferencias finitas, Echeverri y Guarín3 lo hacen usando lenguaje Python considerando simetrías de revolución, Philippe4 ofrece el programa en Maple para una cavidad y Melzani5 sintetiza los temas de estos autores ofreciendo el algoritmo en Phyton para resolver la cavidad o los bordes de un capacitor de placas planas paralelas.
En este trabajo se muestran los resultados de la simulación numérica, usando la ecuación de Laplace en diferencias finitas, para comparar los campos generados por los electrodos planos inclinados de la celda construida y los electrodos planos paralelos. Se discute la utilidad de la representación gráfica obtenida para las consideraciones experimentales posteriores en el Laboratorio de Dieléctricos.
II. MÉTODO DE TRABAJO
La celda de sólidos HP 16451B
Cuando se trabaja con el tomate no entero, ya sea en pasta, trozo o licuado, en el Laboratorio de Dieléctricos se utiliza la celda Hewlett-Packard 16451B para sólidos (Fig. 1). Esta celda comercial tiene dos placas paralelas circulares a la que se le puede variar su separación con un tornillo micrométrico hasta una distancia del orden del centímetro. Las mediciones se realizan colocando entre las placas una caja de Petri con la muestra del tomate. La celda tiene un anillo de guarda que asegura la zona con uniformidad del campo.
La celda de sólidos de Varlan-Sansen
La celda original de Varlan-Sansen considera tres puntos de apoyo, dos son los electrodos y el tercero otorga estabilidad a la fruta (Fig. 2).
El Laboratorio de Dieléctricos utilizó inicialmente esta geometría en mediciones de tomate6 pero, tal como se describe en la Ref.,7 se modificó para quedar actualmente con dos puntos de contacto y apoyo que son los electrodos, inclinados para mejorar el contacto eléctrico con la fruta. Los resultados experimentales fueron presentados por Hemsy et al.8
En la Fig. 3 se muestra la celda que respeta la propuesta de Varlan-Sansen. Los electrodos, dos superficies metálicas circulares de 1 cm de diámetro, se encuentran a potenciales inversos fijados desde el impedanciómetro.
En la modificación realizada se eliminan los cables que unen la celda al equipo como se muestra en la Fig. 4.
Para la medición que se realiza sobre los tomates, interesa que, al separar e inclinar los electrodos, la distorsión del campo generado no introduzca un error grande al comparar los resultados con los obtenidos usando placas planas paralelas. Para realizar este análisis se resolvió la ecuación de Laplace en diferencias finitas, por medio de iteración numérica y se escribió un código en lenguaje Fortran. El programa permite calcular el potencial eléctrico para luego realizar las gráficas de contorno.
Ecuación de Laplace
La ecuación de Laplace se aplica en problemas de física en los cuales el medio es homogéneo y no hay distribuciones volumétricas de carga. Sin embargo, es posible la presencia de cargas puntuales y densidades de carga lineal y superficial. La ecuación de Laplace se escribe:
donde 2 𝑉 es el laplaciano de 𝑉, y 𝑉 el potencial eléctrico.
Para resolver la Ec. (1) por medio de iteración numérica se supone un problema bidimensional en coordenadas 𝑥,𝑦 en el cual se quiere conocer el potencial eléctrico. Se divide el interior de una sección en cuadrados con lados de longitud h. En la Fig. 5 se muestra una porción de esta región. Los valores de potencial 𝑉 0 , 𝑉 1 , 𝑉 2 , 𝑉 3 y 𝑉 4 en los puntos son desconocidos.
La ecuación de Laplace en dos dimensiones en las variables independientes 𝑥,𝑦 es:
Utilizando polinomio de Taylor de grado 3 se puede escribir la siguiente expresión, conocida como ecuación de Laplace en diferencias finitas
Los aspectos teóricos que llevan a plantear la relación (3) entre los potenciales en los puntos indicados en la Fig. 4, se pueden consultar en textos específicos.2,5,9,10 En este trabajo, el interés está centrado en la comparación de los resultados obtenidos con las distintas geometrías.
El potencial obtenido usando la expresión (3) es un resultado aproximado para un caso bidimensional. El resultado mejora cuando h se aproxima a cero, es decir a medida que la región en la que se quiere calcular el potencial se divide en una rejilla con cuadrados de menor tamaño. En este trabajo se ha considerado una grilla uniforme, con puntos separados la misma distancia h.
En forma intuitiva, la expresión (3) indica que el potencial en un punto de la rejilla está representado por un promedio del potencial de los cuatros primeros vecinos. El método iterativo simplemente utiliza la ecuación para determinar el potencial barriendo cada uno de los puntos de la cuadrícula. El proceso del cálculo se repite tantas veces como sea necesario hasta que los valores entre dos cálculos consecutivos de la región no fluctúen dentro de un error predeterminado.
Considerando la Ec. (3) se desarrolló un programa en lenguaje Fortran que permite introducir las condiciones de borde del problema, para que una computadora realice los cálculos de potencial en toda la grilla. Como condiciones de borde del problema se fijaron los potenciales de los electrodos y la ausencia de distribuciones de carga que modifiquen los potenciales en la región.
III. RESULTADOS Y DISCUSIÓN
A continuación se muestran los resultados para la celda adaptada de Varlan-Sansen y para la celda comercial HP 16451B. Luego se compararán los campos a la que es sometida la muestra en el proceso de medición en cada una de las celdas.
Simulación celda modificada de Varlan-Sansen
Con los valores calculados por el programa se graficaron las líneas equipotenciales de la celda utilizada para medir las propiedades eléctricas de fruta entera, con tres valores de separación entre los electrodos.
Si bien las mediciones con la fruta se realizan con potenciales de 0.1 V o 0.5 V, en las simulaciones se consideró que los contactos estaban a 5.0 V y −5.0 V para una mejor representación visual.
En las gráficas de la Fig. 6 se observa la falta de uniformidad de los campos entre electrodos, que se inclinaron a 45º, valor del ángulo tomado como referencia para realizar las mediciones del tomate entero. El comportamiento de estas zonas equipotenciales se deberá incluir en la interpretación circuital de resultados. También se observa que la simetría del sistema se refleja en la simetría de los valores de potencial.
La comparación de las equipotenciales permite definir la mejor zona de trabajo, es decir que se puede relacionar la inclinación y la separación de las placas con el tamaño de la fruta, de modo que se verifiquen las dos condiciones buscadas, buen contacto eléctrico y campo “uniforme” dentro de un error estimable.
Los resultados obtenidos también se pueden representar en 3D. En la Fig. 7 se observa la superficie del potencial eléctrico calculado entre los electrodos, donde el eje vertical indica el potencial en voltios.
La celda de sólidos HP 16451B
La prueba a las que se sometió el programa fue simular el comportamiento de los campos en la celda de electrodos planos paralelos, problema extensamente estudiado. Aquí también se consideraron las placas a 5.0 V y −5.0 V (Fig.8).
Como se observa en la Fig. 8, las líneas equipotenciales son paralelas en la región central (zona de campo uniforme indicada por el fabricante con el anillo de guarda) y luego se empiezan a curvar a medida que se acercan a los bordes. También puede verse una simetría en la gráfica en la recta horizontal que pasa por 0.5 que se corresponde con el potencial 𝑉=0.
IV. CONCLUSIONES
En mediciones dieléctricas por lo general se trabaja con campos uniformes como los generados por la celda Hewlett-Packard 16451B dentro del anillo de guarda.
En este trabajo se simuló el comportamiento del campo generado por la celda modificada de Varlan- Sansen y se observó que los campos son no uniformes. Por consiguiente es que debe definirse la mejor zona de trabajo para introducir el menor error posible por la distorsión del campo en las mediciones. Esto lo que permite es poder comparar los resultados medidos entre una celda y otra, o con resultados publicados por otros grupos donde se trabaja con campos uniformes.