Aplicación del método de los mínimos cuadrados para la obtención de los parámetros de los modelos de Henderson y Chung-Pfost

Application of the least squares method to obtain Henderson and Chung-Pfost's parameters

Ricardo Díaz-Martín^a, Gastón Sanglier-Contreras^a y Adib Guardiola-Mouhaffel^b,*

^a Universidad CEU San Pablo, Campus de Montepríncipe, Madrid, España. *Autor para correspondencia. Correo electrónico: guardiolaadib@gmail.com

Recibido el 26 de mayo de 2014.
Aceptado el 26 de agosto de 2014.

Resumen

Se desarrollan modelos matemáticos que predigan los puntos de equilibrio de humedad de los materiales higroscópicos que habitualmente se transportan, y la explicación de la obtención de los métodos de cálculo de los parámetros mediante regresión lineal de los mínimos cuadrados nos dará una comprensión complementaria.

Palabras clave: Chung-Pfost, Henderson, Mínimos cuadrados, Astillas.

Abstract

Keywords: Chung-Pfost, Henderson, Least squares, Pellets and wood chips.

Introducción

Las astillas son un material higroscópico, lo cual indica que el agua que contiene es constante si se encuentra en equilibrio con el agua en estado de vapor del aire de la atmósfera que la rodea (Prieto García, Gordillo Martínez, Prieto Méndez, Gómez Aldapa y Román Gutérrez, 2006), de modo que las condiciones de los habitáculos o contenedores donde se transporten son determinantes para conservar su buen estado y sus propiedades.

El método de los mínimos cuadrados permite, mediante los pares compuestos por los resultados de la humedad generada por las sales y la humedad adsorbida por la biomasa, calcular los parámetros característicos de los modelos de Henderson y Chung-Pfost utilizados para predicción de estados en el transporte de esta.

Tanto las astillas como otras biomasas derivadas de la madera, como los pellets, son materiales que poseen la característica de adaptarse a la humedad del ambiente en que se encuentren; así, el contenido de humedad se adapta hasta alcanzar el equilibrio (Henderson, 1952). Las curvas de adsorción describen el cambio que sufren las biomasas al verse expuestas a condiciones higrotérmicas determinadas, de aumento progresivo hasta llegar a la humedad de equilibrio; la relación entre el contenido de humedad y la actividad del agua da la información necesaria para modelizar los materiales evaluados (García-Pérez, Cárcel, Clemente y Mulet, 2008).

Las astillas de pino resultan de introducir la madera directamente en una astilladora, la forma básica de producción de biomasa en madera, mientras que el proceso de obtención de los pellets requiere el prensado de la madera; la lignina realiza la función de aglomerante, con lo que aumenta la densidad de la madera, que no tiene que ser necesariamente de la tala de árboles, si no que se puede recuperar de carpinterías y podas.

En la madera del pino, las hemicelulosas que forman parte de la matriz junto con la lignina, donde reside la celulosa, son un 35% y donde se almacena el agua cuando se absorbe; la hemicelulosa y la lignina representan cada una un 24% de la madera.

Dos baños termostáticos (a 15 y 35 °C) sobre diez recipientes con diferentes humedades relativas proporcionan los datos necesarios para la regresión lineal por mínimos cuadrados.

Los datos serán la humedad relativa generada en cada recipiente y la cantidad de agua absorbida por las astillas, según el método gravimétrico estático de las sales ácidas (García Esteban, Guindeo Casasús, de Palacios y García Fernández, 2004).

Modelos matemáticos

La formación de pares entre la predicción y el estado de saturación del agua para cada sal es característico de cada modelo y dan lugar a las rectas de ajuste.

En este caso se ha elegido los modelos de Henderson y Chung-Pfost, aprobados por la Sociedad Americana de Ingenieros Agrícolas para el maíz y el arroz paddy, y los más comunes por su generalidad y relativa precisión (Muñoz et al., 1991), además de ser los únicos que tienen como variable la temperatura, dicho modelos son los únicos que consideran la temperatura en el análisis de isotermas de adsorción (García Fernández, Simón García, Esteban, de Palacios, Fernández-Golfín y Conde, 2013).

A partir de los modelos, se obtienen las respectivas ecuaciones linealizadas.

Ecuación de Henderson modificada

Xeq: contenido de humedad en equilibrio, %. a_w: actividad del agua en el elemento. T: temperatura (°C). C₁: constante de Henderson. C₂: constante de Henderson. (Handbook of Industrial Drying, 2006.)

Fórmula 2, ecuación de Henderson linealizada:

Ecuación de Chung-Pfost modificada

Fórmula 3, ecuación de Chung-Pfost modificada:

]]>

Xeq: contenido de humedad en equilibrio %. a_w: actividad del agua en el elemento. T: temperatura (°C). E: constante de Chung-Pfost. F: constante de Chung-Pfost. (Handbook of Industrial Drying, 2006.)

Fórmula 4, ecuación de Chung-Pfost linealizada:

Parámetros de ajuste de mínimos cuadrados

Fórmula 5, ecuación de ajuste mínimos cuadrados:

Ecuación de Chung-Pfost modificada

Fórmula 6, ecuación de ajuste mínimos cuadrados (método regresivo lineal, cálculo estadístico).

Objetivos

Explicar cómo, a partir de un procedimiento experimental de equilibrios de adsorción de humedad, podemos realizar un comportamiento predictivo de unas muestras de astillas de pino canario mediante los modelos matemáticos de Henderson y Chung-Pfost, desarrollando detalladamente la utilización de la técnica de análisis numérico de los mínimos cuadrados.

Metodología

Se ha utilizado un procedimiento de análisis gravimétrico para evaluar los porcentajes de humedad adsorbidos por las muestras, las cuales fueron previamente expuestas a diez sales ácidas (Li CL, CH₃COOK, MgCl₂, K₂CO₃, Mg (NO₃)₂, Sr Cl₂, Na Cl, K Cl, Ba Cl, K₂S O₄) en recipientes cerrados hasta alcanzar su equilibrio, a las temperaturas de 15 y 35 °C. Cada sal genera una humedad relativa característica en el interior de los recipientes (tablas 1 y 2).

De tal forma utilizamos los resultados finales de los porcentajes de adsorción de humedad de cada una de las diez muestras para aplicarlos en los modelos mencionados anteriormente y confrontarlo con su valor de humedad correspondiente a la sal que haya sido expuesto, formamos pares asociados permitiéndonos desarrollar el método de los mínimos cuadrados para la obtención de los parámetros característicos de cada modelo (tablas 3 y 4).

Resultados

Determinación de parámetros en puntos de equilibrio de adsorción

En el análisis de los resultados de la tabla 5, se obtienen los valores de contenido de humedad en equilibrio a cada temperatura. La construcción de los pares para poder aplicar el método de los mínimos cuadrados se configura de la siguiente manera. Los valores de humedad generados por las sales se le aplican las expresiones.

Aplicación de los datos

La disposición correlativa de los resultados para formar los pares se diferencia según la modelización lineal de los dos modelos; mientras que en el modelo de Henderson se toma el negativo del logaritmo invertido respecto a la recepción de los datos, en el modelo de Chung-Pfost se hace correlativa, tomando el mismo valor directo en valor negativo.

Parámetros Chung-Pfost:

Fórmula 7, resultados de pendiente y coeficiente anexo, para cálculo por regresión en mínimos cuadrados parámetros de la modelización de Chung-Pfost

Parámetros de Henderson:

Una vez calculadas la pendiente y la constante de la recta de ajuste de los parámetros de cada modelo, se recurre a las expresiones matemáticas lineales para determinar los parámetros característicos.

Fórmula 9, fórmulas de Henderson y Chung-Pfost, remarcando el posicionamiento de sus parámetros característicos.

En el caso de la fórmula de Henderson, el valor de C2 tiene como resultado el mismo que la pendiente de ajuste, mientras que en F de Chung-Pfost el valor es la inversa de la pendiente.

Para el segundo parámetro de los modelos, en el caso de Henderson el valor de C3 se eleva al número e; por otro lado, en Chung-Pfost la relación entre la inversa de la pendiente y la inversa del valor auxiliar de la recta da E.

Conclusiones

Mientras que la fórmula de Henderson pertenece a un modelo exponencial inverso, la expresión de Chung-Pfost corresponde a uno lineal, lo que condiciona la capacidad de representar parámetros de sorción de humedad con un error por debajo del 10%.

Mientras el modelo de Henderson está limitado en la modelización de biomasas que tengan sus puntos máximos de humedad de equilibrio reducido, el modelo de Chung-Pfost no es aconsejado para biomasas elevadas en puntos máximos de adsorción.

Como vemos en la figura 1, la pendiente es mayor en la modelización de Henderson, lo que indica que las mediciones son más agudas, pero no por eso más o menos precisas.

Agradecimientos

Este trabajo ha sido posible gracias al apoyo de la empresa Sedam Management, Instituto Tecnológico de Canarias y Tragsa (Empresa de transformación agraria S.A), Oliver González Arias, Miguel Pérez Afonso, Gonzalo Pierna-Vieja Izquierdo, Julián Muñoz de la Fuente.

Referencias

Brunauer, S., Emmett, P.H., y Teller, E. (1938). Adsorption of gases in multimolecular layers. Journal of the American Chemical Society, 60, 308-319.         [ Links ]

García Esteban, A., Guindeo Casasús, P., de Palacios, F., y García Fernández, F. (2004). Histéresis de la madera de Pinus Sylvestris L. para isotermas de 35 °C mediante el método de sales saturadas. Mater construction, 54, 53-56.         [ Links ]

García Fernández, F., Simón García, C., Esteban, L.G., de Palacios, P., Fernández-Golfín, J.I., Conde, M. (2013). Estudio comparativo del método de cámara y de las sales saturadas en la obtención de los puntos de equilibrio higroscópico de limba. Sociedad Española de Ciencias Forestales. VI Congreso Forestal Español.         [ Links ]

García-Pérez, J.V., Cárcel, J.A., Clemente, G., y Mulet, A. (2008). Water sorption isotherms for lemon peel at different temperatures and isosteric heats. Food Science and Technology, 41, 18-25.         [ Links ]

Handbook of Industrial Drying (2007) (3.a ed.). Cleveland: CRC.         [ Links ]

Henderson, S.M. (1952). A basic concept of equilibrium moisture. Agriculture Engeneering, 33, 29-32.         [ Links ]

Kaymak-Ertekin, F., y Gedik, A. (2004) Sorption isotherms andisosteric heat of sorption for grapes, apricots, apples and potatoes. Lebensmittel Wissenschaft und Technologie, 37, 429-438.         [ Links ]

Muñoz, R. (1991 ). Un modelo propuesto para isotermas en equilibrio. Agricultura Técnica (Chile), 51, 87-90.         [ Links ]

Prieto García, F., Gordillo Martínez, A.J., Prieto Méndez, J., Gómez Aldapa, C.A., y Román Gutérrez A D. (2006). Evaluación de las isotermas de sorción en cereales para el desayuno. Avances en Ciencias e lngeniería, 3, 137-150.         [ Links ]