ES2376675T3 - Método de control del secado secundario en un proceso de secado por congelación. - Google Patents
Método de control del secado secundario en un proceso de secado por congelación. Download PDFInfo
- Publication number
- ES2376675T3 ES2376675T3 ES08013243T ES08013243T ES2376675T3 ES 2376675 T3 ES2376675 T3 ES 2376675T3 ES 08013243 T ES08013243 T ES 08013243T ES 08013243 T ES08013243 T ES 08013243T ES 2376675 T3 ES2376675 T3 ES 2376675T3
- Authority
- ES
- Spain
- Prior art keywords
- time
- residual moisture
- desorption
- product
- water
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F26—DRYING
- F26B—DRYING SOLID MATERIALS OR OBJECTS BY REMOVING LIQUID THEREFROM
- F26B5/00—Drying solid materials or objects by processes not involving the application of heat
- F26B5/04—Drying solid materials or objects by processes not involving the application of heat by evaporation or sublimation of moisture under reduced pressure, e.g. in a vacuum
- F26B5/06—Drying solid materials or objects by processes not involving the application of heat by evaporation or sublimation of moisture under reduced pressure, e.g. in a vacuum the process involving freezing
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Molecular Biology (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Drying Of Solid Materials (AREA)
Abstract
Método de control de una fase de secado secundario de un proceso de secado por congelación en un aparato secador por congelación que incluye una cámara de secado que contiene un producto a secar y que puede ser aislada para realizar ensayos de aumento de presión, comprendiendo dicho método las etapas de: - realizar un primer ensayo de aumento de presión en un tiempo t=t0 y calcular un primer valor de tasa (DRexp,0) de desorción experimental de dicho producto (etapa 1); - realizar un segundo ensayo de aumento de presión en un tiempo t=t1 y calcular un segundo valor de tasa (DRexp,1) de desorción experimental de dicho producto (etapa 2); - realizar un tercer ensayo de aumento de presión en un tiempo t=t2 y calcular un tercer valor de tasa (DRexp,2) de desorción experimental de dicho producto (etapa 3); - estimar las condiciones iniciales (CS,0) y las constantes cinéticas (k0, k1, k2) de un modelo cinético del proceso de secado, siendo adecuado dicho modelo cinético para calcular un contenido (CS) de humedad residual y/o una tasa (DRteor) de desorción de dicho producto (etapa 4); - calcular en un tiempo t=t2 un contenido (CS,2) de humedad residual respectivo y una tasa (DRteor,2) de desorción respectiva (etapa 5).
Description
Método de control del secado secundario en un proceso de secado por congelación
La invención se refiere a métodos de control de un proceso de secado por congelación en un secador por congelación; de forma específica, la misma se refiere a un método de control del secado secundario de un proceso de secado por congelación, por ejemplo, de productos farmacéuticos dispuestos en recipientes.
El secado por congelación, también conocido como liofilización, es un proceso de deshidratación que permite la eliminación por sublimación de agua y/o disolventes de una sustancia, tal como productos alimentarios, farmacéuticos o biológicos. De forma típica, el proceso de secado por congelación se usa para conservar un producto perecedero, ya que el contenido en agua reducido en gran medida resultante evita la acción de microorganismos y enzimas que normalmente estropearían o degradarían el producto. Además, el proceso hace que el producto sea más fácil de transportar. Los productos secados por congelación pueden ser precintados en recipientes para evitar la reabsorción de humedad y pueden ser rehidratados o reconstituidos fácilmente mediante la aplicación del agua o los disolventes eliminados. De esta manera, el producto puede ser almacenado a temperatura ambiente sin refrigeración y quedar protegido contra su deterioro durante muchos años.
Debido a que el secado por congelación es un proceso de baja temperatura en el que la temperatura del producto no excede de forma típica los 30ºC durante las fases operativas, el mismo provoca menos daños o degradación en el producto que otros procesos de deshidratación que usan temperaturas más altas. Normalmente, el secado por congelación no provoca una contracción o endurecimiento significativos del producto secado. Los productos secados por congelación pueden ser rehidratados de forma mucho más rápida y fácil debido a la estructura porosa creada durante la sublimación del hielo.
En el campo farmacéutico, el proceso de secado por congelación es usado ampliamente en la producción de fármacos, principalmente de administración parenteral y oral, también porque el proceso de secado por congelación puede ser llevado a cabo en condiciones estériles.
Un aparato secador por congelación conocido para realizar un proceso de secado por congelación comprende normalmente una cámara de secado y una cámara condensadora interconectadas por un conducto que está dotado de una válvula que permite aislar la cámara de secado durante el proceso en caso necesario.
La cámara de secado comprende una pluralidad de estantes de temperatura controlada dispuestos para alojar recipientes del producto a secar. La cámara condensadora incluye placas o bobinas condensadoras que tienen unas superficies que se mantienen a una temperatura muy baja, p. ej., -50ºC, mediante un dispositivo refrigerante o de congelación. La cámara condensadora también está conectada a una o más bombas de vacío para conseguir valores de vacío elevados en el interior de ambas cámaras.
De forma típica, el proceso de secado por congelación comprende tres fases: una fase de congelación, una fase de secado principal y una fase de secado secundario.
Durante la fase de congelación, la temperatura de almacenamiento se reduce de forma típica hasta -30/-40ºC para convertir en hielo la mayor parte de agua y/o disolventes contenidos en el producto.
En la fase de secado principal, la temperatura de almacenamiento aumenta, mientras la presión en el interior de la cámara de secado disminuye por debajo de 1-5 mbar, a efectos de permitir que el agua y/o disolventes congelados en el producto se sublimen directamente de fase sólida a fase gaseosa. La aplicación de un vacío elevado hace posible la sublimación del agua a temperaturas bajas.
Se suministra calor al producto y el vapor generado por la sublimación del agua y/o los disolventes congelados es eliminado de la cámara de secado mediante las placas o bobinas condensadoras de la cámara condensadora, donde el vapor puede resolidificarse. La fase de secado secundario se realiza para eliminar por desorción la humedad residual del producto, es decir, la cantidad de agua y/o disolventes no congelados que no puede ser eliminada durante el secado principal, cuando tiene lugar la sublimación del hielo. Durante esta fase, la temperatura de almacenamiento aumenta adicionalmente hasta un máximo de 30-60ºC para calentar el producto, mientras la presión en el interior de la cámara de secado se establece de forma típica por debajo de 0,1 mbar.
Al final de la fase de secado secundario, el producto está suficientemente seco, de forma típica, con un contenido de humedad residual del 1-3%.
El secado secundario debe ser controlado de forma precisa para señalar el momento en el que el proceso de secado ha finalizado, es decir, el momento en el que se ha alcanzado la cantidad deseada de humedad residual en el producto.
Existen métodos de control conocidos de una fase de secado secundario.
Según un método conocido, es posible determinar la humedad residual del producto extrayendo muestras del secador por congelación sin interrumpir el secado por congelación (p. ej., usando un “captador de muestras”) y midiendo fuera de línea su contenido en humedad mediante valoración Karl Fischer, análisis gravimétrico térmico o espectroscopia próxima al infrarrojo.
US 6971187 propone otro método, en el que la estimación de la tasa de secado del producto durante el secado secundario se obtiene realizando un ensayo de aumento de presión (EAP).
Durante un EAP, la cámara de secado es aislada de la cámara condensadora cerrando la válvula colocada en el conducto que conecta las dos cámaras. Debido a que el calentamiento no se detiene, la sublimación del hielo continúa, aumentando por lo tanto en la cámara de secado la presión medida.
Teniendo en cuenta la curva de presión-tiempo, la pendiente al inicio de esta curva permite estimar el caudal de agua y/o disolvente procedente del producto mediante la ecuación:
donde:
P: presión medida, [Pa]
t: tiempo, [s]
t0: instante de tiempo al inicio del EAP, [s]
R: constante de gas [8,314 J mol-1 K-1]
T: temperatura del vapor, [K]
V: volumen (libre) de la cámara, [m3]
jw,n: caudal de agua y/o disolvente procedente del producto, [mol s-1]
Por lo tanto, es posible calcular el flujo de masa de agua y/o disolvente:
donde:
jw,m: flujo de masa de agua y/o disolvente procedente del producto, [kg s-1]
Mw: peso molecular de agua y/o disolvente, [kg mol-1]
A partir de este valor, es posible estimar la pérdida en agua y/o disolvente durante el periodo de medición
transcurrido entre dos EAP consecutivos, mediante:
donde:
ftj = tj-tj-1 tiempo transcurrido entre EAP j y EAP (j-1), [s]
fwm,j: pérdida en agua durante el intervalo de tiempo ftj, [kg] jw,m, j-1: flujo de masa de agua y/o disolvente procedente del producto calculado a partir de EAP (j-1), [kg s-1].
La cantidad total de agua y/o disolvente eliminada entre un tiempo de referencia (p. ej., al inicio del secado secundario) y cualquier tiempo tj de interés determinado es simplemente la suma de toda la fwm,j que se produce en los diversos intervalos entre EAP. Utilizando un valor experimental independiente para detectar el contenido de agua residual en un tiempo de referencia, p. ej., al final del secado principal, es posible calcular el contenido real de humedad con respecto al tiempo en tiempo real. Esto requiere extraer una muestra de la cámara de secado o usar detectores caros (p. ej., detectores basados en NIR) para obtener este valor en línea.
Teniendo en cuenta este valor experimental, se usan algunas indicaciones empíricas o de sentido común para calcular la temperatura “óptima” para minimizar el tiempo necesario para finalizar el secado secundario.
Un inconveniente de los métodos anteriores consiste en que los mismos requieren extraer muestras de la cámara de secado y usar detectores caros para medir los valores experimentales de agua residual y/o disolvente. La extracción de muestras es una operación invasiva que perturba el proceso de secado por congelación y, por lo tanto, no resulta adecuada en procesos estériles y/o asépticos y/o cuando se usa una carga/descarga automática de los recipientes. Además, la extracción de muestras consume tiempo y requiere operarios experimentados.
Otro inconveniente del método descrito en US 6971187 consiste en que las indicaciones empíricas y de sentido común usadas para calcular la temperatura “óptima” no permiten optimizar el proceso.
En US 6176121 se da a conocer una realización diferente, en la que, mediante el uso de dos mediciones sucesivas de la tasa (DR) de desorción, es decir, del caudal de masa de vapor de agua y/o disolvente debido a la desorción, calculada a partir de jw,m, es posible extrapolar el punto en el tiempo en el que se obtiene un valor pequeño determinado de DR. Para conseguirlo, la válvula situada entre la cámara de secado y la cámara condensadora debería cerrarse regularmente durante cierto tiempo, siendo necesario obtener la curva de aumento de presión (CAP) provocada por el vapor de agua de la desorción. Por lo tanto, es posible calcular la masa de agua y/o disolvente desorbidos con respecto al tiempo o la tasa de desorción a partir de la pendiente inicial de la CAP tal como sigue:
donde:
mseca: masa del producto seco, [kg]
DRexp: tasa de desorción experimental, [% de agua y/o disolvente con respecto a producto seco s-1]
Un inconveniente de este método consiste en que, debido a la gran simplificación de la realización, se ha comprobado que el mismo falla con respecto al final del secado secundario. Además, el mismo no permite estimar la humedad residual absoluta, sino solamente la diferencia con respecto a la humedad de equilibrio, que depende de las condiciones de funcionamiento (temperatura de almacenamiento y presión de la cámara de secado), no siendo posible por lo tanto establecer un objetivo sobre este valor.
Un objetivo de la invención consiste en mejorar los métodos de control de un proceso de secado por congelación en un secador por congelación, especialmente para controlar una fase de secado secundario de dicho proceso de secado por congelación.
Otro objetivo consiste en dar a conocer un método para calcular parámetros de proceso, tales como el contenido de humedad residual y/o la tasa de desorción de un producto seco, que no es invasivo y no perturba el proceso de secado por congelación y, por lo tanto, es adecuado para su uso en procesos estériles y/o asépticos y/o cuando se usa una carga/descarga automática de los recipientes. Otro objetivo consiste en dar a conocer un método que permite estimar de forma precisa las condiciones iniciales y las constantes cinéticas de un modelo cinético del proceso de secado, adecuado para calcular los parámetros de proceso.
Otro objetivo consiste en dar a conocer un método para estimar de manera fiable y precisa una concentración de humedad residual y/o una tasa de desorción del producto seco durante una fase de secado secundario y el tiempo necesario para finalizar dicha fase de secado secundario.
Otro objetivo consiste en dar a conocer un método en el que la estimación de los parámetros de proceso se mejora y refina progresivamente durante el progreso de la fase de secado secundario, resultando no obstante dicha estimación buena con respecto a los métodos conocidos, incluso al inicio de la fase de secado secundario.
Según la invención, se da a conocer un método de control de una fase de secado secundario de un proceso de secado por congelación en un aparato secador por congelación que incluye una cámara de secado que contiene un producto a secar y que puede ser aislada para realizar ensayos de aumento de presión, comprendiendo dicho método las etapas de:
- -
- realizar un primer ensayo de aumento de presión en un tiempo t=t0 y calcular un primer valor de tasa de desorción experimental de dicho producto (etapa 1);
- -
- realizar un segundo ensayo de aumento de presión en un tiempo t=t1 y calcular un segundo valor de tasa de desorción experimental de dicho producto (etapa 2);
- -
- realizar un tercer ensayo de aumento de presión en un tiempo t=t2 y calcular un tercer valor de tasa de desorción experimental de dicho producto (etapa 3);
- -
- estimar las condiciones iniciales y las constantes cinéticas de un modelo cinético del proceso de secado, siendo adecuado dicho modelo cinético para calcular un contenido de humedad residual y/o una tasa de desorción de dicho producto (etapa 4);
- -
- calcular en un tiempo t=t2 un contenido de humedad residual respectivo y una tasa de desorción respectiva (etapa 5).
El método comprende además, después de la etapa 5, la etapa de:
- -
- comparar dicho contenido de humedad residual y/o dicha tasa de desorción calculados en un tiempo t=t2 con una concentración de humedad residual final deseada y/o una tasa de desorción final deseada, respectivamente, (etapa 6); si dicho contenido de humedad residual es inferior o igual con respecto a dicha concentración de humedad residual final o dicha tasa de desorción es inferior o igual con respecto a dicha tasa de desorción final, considerar la fase de secado secundario finalizada; en caso contrario, el método comprende además las etapas de:
- -
- estimar un tiempo final en el que se obtiene dicha concentración de humedad residual final o dicha tasa de desorción final (etapa 7);
- -
- realizar un ensayo de aumento de presión adicional en un tiempo t=tj y calcular en dicho tiempo t=tj un contenido de humedad residual respectivo y una tasa de desorción respectiva (etapa 8);
- -
- estimar las condiciones iniciales y las constantes cinéticas de dicho modelo cinético (etapa 9);
- -
- calcular en dicho tiempo t=tj dicho contenido de humedad residual respectivo y/o dicha tasa de desorción respectiva (etapa 10);
- -
- comparar dicho contenido de humedad residual y/o dicha tasa de desorción calculados en dicho tiempo t=tj con dicha concentración de humedad residual final y/o dicha tasa de desorción final, respectivamente, (etapa 11); si dicho contenido de humedad residual es inferior o igual con respecto a dicha concentración de humedad residual final o dicha tasa de desorción es inferior o igual con respecto a dicha tasa de desorción final, considerar la fase de secado secundario finalizada; en caso contrario, repetir las etapas 7 a 11.
Gracias a la invención, es posible obtener un método para calcular de manera fiable y precisa la concentración de humedad residual y/o la tasa de desorción de un producto seco durante una fase de secado secundario de un proceso de secado por congelación. El método también permite estimar de forma precisa las condiciones iniciales y las constantes cinéticas de un modelo cinético del proceso de secado, que calcula la concentración de humedad residual y/o la tasa de desorción del proceso, sin extraer ninguna muestra de la cámara de secado y sin usar detectores caros para obtener este valor en línea. Por lo tanto, el método de control de la invención no es invasivo y no perturba el proceso de secado por congelación, y es adecuado para su uso en procesos estériles y/o asépticos y/o cuando se usa una carga/descarga automática de los recipientes. Además, el método permite calcular el tiempo necesario para finalizar dicha fase de secado secundario, pudiendo consistir el requisito de detención en que la concentración de humedad residual o la tasa de desorción tenga un valor final deseado respectivo. Debido a que las etapas del método se repiten hasta que se alcanza el final de la fase de secado secundario, la estimación de los parámetros de proceso se mejora y refina progresivamente durante el progreso de la fase de secado secundario, resultando no obstante dicha estimación buena con respecto a los métodos conocidos, incluso al inicio de la fase de secado secundario.
Es posible mejorar la comprensión de la invención y su puesta en práctica haciendo referencia a los dibujos adjuntos, que muestran una realización de la invención a título de ejemplo no limitativo, en los que:
La Figura 1 es un diagrama de flujo que muestra esquemáticamente el método de control de una fase de secado secundario en un proceso de secado por congelación de la invención;
La Figura 2 es un gráfico que muestra una secuencia de valores medidos experimentales de la tasa de desorción con respecto al tiempo durante el secado secundario;
La Figura 3 son gráficos que muestran una estimación de evolución temporal de la concentración de humedad residual y la tasa de desorción de un producto seco en un tiempo definido, respectivamente;
La Figura 4 son gráficos que muestran una estimación de evolución temporal de la concentración de humedad residual y la tasa de desorción de un producto seco en un tiempo definido adicional, respectivamente;
La Figura 5 es un gráfico que muestra una secuencia de evolución temporal de estimaciones de tiempo necesario para finalizar el secado secundario;
La Figura 6 muestra una comparación entre estimaciones de tiempo necesario para finalizar el secado secundario obtenidas usando el método de la invención y usando el método según US 6176121;
Las Figuras 7 y 8 muestran una comparación entre valores experimentales y valores previstos por el método de la invención de la tasa de desorción y del contenido de agua residual, respectivamente.
El método de la invención controla una fase de secado secundario de un proceso de secado por congelación en un secador por congelación. De forma específica, el método calcula el contenido de humedad residual de un producto seco y permite obtener una estimación fiable del tiempo que es necesario para finalizar esta fase según el objetivo deseado (contenido de humedad final y/o valor final de la tasa de desorción).
El método requiere realizar periódicamente un ensayo de aumento de presión (EAP) y, por lo tanto, puede aplicarse en procesos de secado por congelación que se llevan a cabo en secadores por congelación que comprenden una cámara de secado, donde se dispone el producto a secar, y una cámara condensadora separada, donde circula el vapor generado por el proceso de secado y donde es posible resolidificarlo o congelarlo.
El EAP se lleva a cabo cerrando durante un intervalo de tiempo corto (de unas pocas decenas de segundos, p. ej., 30 s, a unos pocos minutos) una válvula que está dispuesta en el conducto que conecta la cámara de secado a la cámara condensadora y midiendo (y registrando) la evolución temporal de la presión total en la cámara.
A partir de la pendiente de la curva al inicio del ensayo, es posible calcular la tasa (DR, % s-1) de desorción de agua y/o disolvente actual. El EAP se repite cada intervalo de tiempo especificado previamente (p. ej., 30 minutos) para conocer la evolución temporal de la tasa de desorción de agua y/o disolvente. El intervalo de tiempo puede ser constante o puede cambiar durante la operación.
Todos los métodos basados en el EAP para controlar la etapa de secado principal de un proceso de secado por congelación se aprovechan de que, durante el ensayo, la presión en la cámara de secado aumenta hasta alcanzar el equilibrio. Como este no es el caso del secado secundario (debido a los valores reducidos de caudal de agua y/o disolvente), la única información que es posible utilizar a partir del EAP es la estimación del caudal de agua y/o disolvente, que por lo tanto puede ser integrada para evaluar la pérdida de agua y/o disolvente con el tiempo. La estimación del contenido de humedad requiere conocer la concentración de humedad inicial, que se calcula según el método de la invención, tal como se describe de forma detallada a continuación, sin extraer ninguna muestra de la cámara de secado y sin usar detectores caros para obtener este valor en línea. En otras palabras, el método de control no es invasivo y no perturba el proceso de secado por congelación y, por lo tanto, es adecuado para su uso en procesos estériles y/o asépticos y/o cuando se usa una carga/descarga automática de los recipientes.
El método de la invención requiere modelizar la dependencia de la tasa (DR) de desorción con respecto al contenido (CS) de humedad residual en el producto seco. Es posible usar varias ecuaciones matemáticas para tal fin. El método comprende un algoritmo que puede funcionar de forma eficaz independientemente de la correlación usada.
Se han propuesto varios modelos cinéticos para modelizar la tasa de desorción de agua y/o disolvente. Es posible asumir que la tasa de desorción depende del contenido de humedad residual o de la diferencia entre el contenido de humedad residual y el valor de equilibrio.
Se ha demostrado que ambos tipos de modelo se comportan más o menos de la misma manera; además, existe incertidumbre sobre el mecanismo físico real de la desorción del agua y/o disolvente, que puede depender del producto considerado.
En una primera versión del método, se asume que la tasa DR de desorción depende de la humedad residual CS en la matriz sólida del producto seco según la ecuación:
siendo posible calcular la evolución temporal de la humedad residual CS, dada en % de agua y/o disolvente con respecto a masa seca, mediante la integración de la siguiente ecuación diferencial:
donde t es el tiempo [s] y k es la constante cinética del proceso [s-1].
La constante cinética puede ser una función de la temperatura y, por lo tanto, puede cambiar con el tiempo, ya que la temperatura del producto puede cambiar con el tiempo, especialmente al inicio del secado secundario, cuando la temperatura aumenta desde el valor usado durante el secado principal hasta el del secado secundario.
Si un EAP se realiza en un tiempo t= tj-1 y el EAP sucesivo se realiza en un tiempo t=tj, y se asume que la temperatura del producto, que varía ligeramente en el intervalo [tj-tj-1], es constante e igual a un valor promedio, es posible describir la variación de la concentración de humedad en el sólido mediante la ecuación:
La solución de la ec. 7 requiere la condición inicial, es decir, el valor de la humedad residual CS en un tiempo t=tj-1:
Es posible calcular el valor de CS,j-1 a partir de la integración de tiempo de la ec. 6 en el intervalo de tiempo anterior:
10 y, por lo tanto:
Es posible repetir este procedimiento hasta que aparece el valor CS,0 de la humedad residual al inicio de la fase de secado secundario (t=t0). Por lo tanto, en el intervalo de tiempo entre tj y tj-1, la evolución de la concentración de humedad residual viene dada por:
La solución de la ec. 11 requiere el valor de la concentración CS,0 de humedad inicial.
Por lo tanto, la evolución del valor teórico de la tasa de desorción en el intervalo de tiempo entre tj y tj-1 viene dada por:
20 Si CS,0 y los valores de las diversas kj son perfectamente conocidos y el modelo dado por la ec. 6 es adecuado para describir la dinámica del sistema, es posible usar la ec. 11 para conocer la evolución temporal del contenido de humedad residual y, por lo tanto, el tiempo que es necesario para cumplir los requisitos con respecto al valor final del contenido de humedad del producto. Si el requisito es con respecto al valor de la tasa de desorción, es posible usar la ec. 12 para tal fin.
25 La situación anterior es bastante rara, ya que el valor de la concentración de humedad inicial debe medirse extrayendo muestras y las diversas constantes cinéticas nunca son conocidas a priori.
El método según la invención permite calcular la condición inicial CS,0 y las constantes cinéticas realizando las siguientes etapas, tal como se muestra en el diagrama de flujo de la Figura 1.
Etapa 1
En un tiempo t=t0, se realiza un EAP y se calcula una tasa DR de desorción respectiva (indicada a continuación
como DRexp,0), es decir, usando la ec. 4.
A partir de la ec. 12, la misma es:
Etapa 2
En un tiempo t=t1, se realiza un EAP y se calcula una tasa DR de desorción respectiva (indicada a continuación 10 como DRexp,1), es decir, usando la ec. 4.
A partir de la ec. 12, la misma es:
Etapa 3
En un tiempo t=t2, se realiza un EAP y se calcula la tasa DR de desorción (indicada a continuación como DRexp,2), es15 decir, usando la ec. 4.
A partir de la ec. 12, la misma es:
Etapa 4
Los valores de CS,0, k0, k1y k2 se estiman de modo que los valores calculados de la tasa de desorción concuerdan 20 con todos los valores experimentales disponibles (DRexp,0, DRexp,1 y DRexp,,2). Esto puede hacerse usando un algoritmo de minimización para resolver el siguiente problema de mínimos cuadrados no lineales:
25 y asumiendo, por ejemplo, que k2 es igual a k1, debido al hecho de que el intervalo de tiempo entre dos EAP es generalmente pequeño, p. ej., 30 minutos, y al hecho de que la temperatura del producto es casi constante durante el secado secundario (la temperatura del producto solamente varía al inicio del secado secundario, desde la del secado principal hasta la necesaria por el secado secundario, aunque esta variación es generalmente lenta, debido a la inercia térmica del sistema).
30 Como valores de inicio en la ec. 16 es posible usar las aproximaciones no exactas a k0, k1y CS,0, que pueden ser calculadas a partir de la ec. 13 y la ec 14 después de los dos primeros EAP:
Estos valores son simplemente una aproximación a las constantes cinéticas k0y k1 y al contenido CS,0 de humedad residual; estas aproximaciones se refinarán después de cada EAP.
Etapa 5
Una vez estimados los valores de CS,0, k0, k1y k2, en un tiempo t=t2, se calcula la concentración CS,2 de humedad residual, usando la ec. 11, o la tasa DRteor,2 de desorción, usando la ec. 12.
Etapa 6
La concentración CS,2 de humedad residual calculada o la tasa DRteor,2 de desorción se compara con un valor deseado de la concentración Cs,f de humedad residual final u objetivo o con un valor deseado de la tasa DRf de desorción final u objetivo.
Si la concentración CS,2 de humedad residual calculada o la tasa DRteor,2 de desorción es inferior o igual con respecto a la concentración CS,f de humedad residual final o la tasa DRf de desorción final, la fase de secado secundario finaliza.
Etapa 7
Si la concentración CS,2 de humedad residual calculada es superior a la concentración CS,f de humedad residual final, o si la tasa DRteor,2 de desorción calculada es superior a la tasa DRf de desorción final, usando los valores calculados de CS,0 y de las constantes cinéticas k0, k1y k2 es posible estimar el tiempo final tf en el que se obtiene la concentración CS,f de humedad residual deseada o la tasa DRf de desorción final, asumiendo que la temperatura del producto no cambia. Esto puede hacerse usando la ec. 11, donde CS se sustituye por CS,f y, por lo tanto, t se corresponde con tf:
Es posible asumir un criterio de detención diferente, p. ej., el requisito de que la tasa de desorción tenga un valor bajo determinado. Para tal fin, es posible usar la ec. 12, donde DR se sustituye por el valor objetivo y, por lo tanto, t se corresponde con tf.
Etapa 8
Se realiza un nuevo EAP en un tiempo t=tj y se calcula una tasa DRexp,j de desorción respectiva; a partir de la ec. 12:
Esta etapa puede repetirse varias veces, tal como se explica mejor a continuación, y después de cada EAP un nuevo valor de DR está disponible y se obtiene una mejor estimación de los valores de CS,0, k0, k1, ..., kj y tf, hasta el final de la fase de secado secundario.
Por ejemplo, en un tiempo t=t3, el EAP permite obtener DRexp,3, y a partir de la ec. 12: Etapa 9
Los valores de las constantes CS,0, k0, k1, ..., kj se estiman resolviendo el problema de mínimos cuadrados no lineales:
asumiendo, por ejemplo, que kj es igual a kj-1, tal como se ha mencionado anteriormente.
Por ejemplo, en un tiempo t=t3, los valores de las constantes CS,0, k0, k1, k2y k3 se calculan resolviendo el problema de mínimos cuadrados no lineales:
Etapa 10
Una vez estimados los valores de CS,0, k0, k1, ..., kj, es posible calcular en un tiempo t=tj la concentración CS,j de humedad residual, usando la ec. 11, o la tasa DRteor,j de desorción, usando la ec. 12.
Por ejemplo, en un tiempo t=t3, una vez estimados los valores de CS,0, k0, k1, k2y k3, es posible calcular la concentración CS,3 de humedad residual o la tasa DRteor,3 de desorción.
Etapa 11
El valor calculado de la concentración CS,3 de humedad residual o la tasa DRteor,3 de desorción se comparan en un tiempo t=tj con la concentración CS,f de humedad residual final o la tasa DRf de desorción final.
Si la concentración CS,j de humedad residual calculada o la tasa DRteor,j de desorción es inferior o igual con respecto a la concentración CS,f de humedad residual final o la tasa DRf de desorción final, la fase de secado secundario finaliza.
Si la concentración CS,j de humedad residual calculada o la tasa DRteor,j de desorción es superior a la concentración CS,f de humedad residual final o la tasa DRf de desorción final, la etapa 7 se repite con t=tj para estimar el tiempo final tf en el que se obtiene la concentración CS,f de humedad residual final o la tasa DRf de desorción final:
Por ejemplo, en un tiempo t=t3 es posible estimar el tiempo final tf en el que se obtiene la concentración CS,f de humedad residual deseada o la tasa DRf, de desorción final, asumiendo que la temperatura del producto no cambia. Esto puede hacerse usando la ec. 11, donde CS se sustituye por CS,f y, por lo tanto, t se corresponde con tf:
Es posible asumir un criterio de detención diferente, es decir, el requisito de que la tasa de desorción tenga un valor bajo final determinado. Para tal fin, es posible usar la ec. 12, donde DR se sustituye por el valor objetivo y, por lo 5 tanto, t se corresponde con tf. Las etapas 7 a 11 se repiten hasta el final de la fase de secado secundario, es decir, hasta que el valor estimado de la concentración CS,j de humedad residual o la tasa DRteor,j de desorción en un tiempo tj es inferior o igual con respecto al valor deseado de concentración CS,f de humedad residual o de tasa DRf de desorción. En una segunda versión del método, se asume que la tasa DR de desorción depende de la diferencia entre el contenido CS de humedad residual en la matriz sólida del producto seco y la concentración CS,eq de
10 humedad de equilibrio:
15 La concentración Cs,eq de humedad de equilibrio es un parámetro adicional cuyo valor puede ser conocido (debe ser determinado experimentalmente).
Empezando por esta expresión diferente de la tasa de desorción y repitiendo el mismo procedimiento descrito anteriormente, es posible obtener resultados similares.
La constante cinética k puede ser una función de la temperatura y puede cambiar con el tiempo; también la
20 concentración Cs,eq de humedad de equilibrio puede cambiar con el tiempo, asumiéndose que esta variación es despreciable durante el intervalo de tiempo entre un EAP y el sucesivo, permitiendo por lo tanto una solución analítica de la ecuación de balance de masa.
Si un EAP se realiza en t=tj-1 y el EAP sucesivo se realiza en t=tj, la evolución de la concentración de humedad residual, en % de agua y/o disolvente con respecto a masa seca, viene dada en el intervalo [tj-tj-1] por la integración 25 de la siguiente ecuación diferencial:
La solución de la ec. 20 requiere la condición inicial, es decir, el valor de la humedad residual CS en t= tj-1:
Es posible calcular el valor de CS,j-1 a partir de la integración de tiempo de la ec. 20 en el intervalo de tiempo anterior:
y, por lo tanto: De forma similar, CS,j-2, que es necesaria para conseguir CS,j-1, puede ser calculada tal como sigue:
Este procedimiento puede repetirse hasta que aparece el valor de la humedad residual CS,0 al inicio de la etapa de secado secundario (t=t0):
Por lo tanto, en el intervalo de tiempo entre tj y tj-1, es posible obtener la evolución de la concentración de humedad residual en función de CS,0, CS,eq,r (con r = 1,..,j)y kr (con r=1,..,j).
La evolución del valor teórico de la tasa de desorción en el intervalo de tiempo entre tj y tj-1 viene dada por:
y, por lo tanto, la misma es una función de CS,0, CS,eq,r (con r=1,..,j)y de kr (con r=1,..,j).
Si CS,0 y los valores de las diversas constantes cinéticas kj son perfectamente conocidos y el modelo dado por la ec. 20 es adecuado para describir la dinámica del sistema, es posible usar la ec. 21 para conocer la evolución temporal del contenido de humedad residual y, por lo tanto, el tiempo que es necesario para cumplir los requisitos con
15 respecto al valor final del contenido de humedad residual del producto. Si el requisito es con respecto al valor de la tasa de desorción, es posible usar la ec. 26 para tal fin.
La situación anterior es bastante rara, ya que el valor de la concentración de humedad inicial debe calcularse extrayendo muestras y las diversas constantes cinéticas nunca son conocidas a priori.
El método según la invención permite calcular la condición inicial CS,0 y las constantes cinéticas realizando las 20 siguientes etapas, tal como se muestra en el diagrama de flujo de la Figura 1.
Etapa 1 En un tiempo t=t0, se realiza un EAP y se calcula la tasa DR de desorción (indicada a continuación como DRexp,0), p. ej., usando la ec. 4.
A partir de la ec. 20, la misma es:
Etapa 2
En un tiempo t=t1, se realiza un EAP y se calcula la tasa DR de desorción (indicada a continuación como DRexp,1), p.
ej., usando la ec. 4.
A partir de la ec. 26, la misma es:
Etapa 3
En un tiempo t=t2, se realiza un EAP y se calcula la tasa DR de desorción (indicada a continuación como DRexp,2), p.
ej., usando la ec. 4.
A partir de la ec. 26, la misma es:
Etapa 4
Los valores de CS,0, k0, k1y k2 se estiman de modo que los valores calculados de las tasas de desorción concuerdan con todos los valores experimentales disponibles (DRexp,0, DRexp,1 y DRexp,2). Esto puede hacerse usando un algoritmo de minimización para resolver el siguiente problema de mínimos cuadrados no lineales:
asumiendo, por ejemplo, que k2 es igual a k1, tal como se ha mencionado anteriormente. Los valores de CS,eq,0, CS,eq,1 y CS,eq,2 deben ser conocidos (a partir de experimentación). Etapa 5
Una vez estimados los valores de CS,0, k0, k1y k2, es posible calcular en un tiempo t=t2 la concentración CS,2 de humedad residual (o la tasa de desorción), usando la ec. 20.
Etapa 6
El valor calculado de la concentración CS,2 de humedad residual se compara con un valor deseado de la concentración CS,f de humedad residual final.
Si el valor calculado de la concentración CS,2 de humedad residual es inferior o igual con respecto a la concentración CS,f de humedad residual final, la fase de secado secundario finaliza.
Etapa 7
Si el valor calculado de la concentración CS,2 de humedad residual es superior a la concentración CS,f de humedad residual final deseada, usando los valores calculados de CS,0 y de las constantes cinéticas es posible estimar el tiempo tf en el que se obtiene el valor deseado de la concentración CS,f de humedad residual, asumiendo que la temperatura del producto no cambia. Esto puede hacerse usando la ec. 21, donde CS se sustituye por CS,f y, por lo tanto, t se corresponde con tf. En este caso, debe resolverse la siguiente ecuación no lineal:
Es posible asumir un criterio de detención diferente, p. ej., el requisito de que la tasa DR de desorción tenga un valor DRf bajo final determinado. Para tal fin, es posible usar la ec. 26, donde DR se sustituye por la tasa DRf de desorción final.
Etapa 8
Se realiza un nuevo EAP en un tiempo t=tj y se calcula una tasa DRexp,j de desorción respectiva; a partir de la ec. 26:
Esta etapa puede repetirse varias veces, y después de cada EAP un nuevo valor de DR está disponible y se obtiene una mejor estimación de los valores de CS,0, k0, k1, ..., kj y tf, hasta el final de la fase de secado secundario.
Por ejemplo, en un tiempo t=t3, el EAP permite obtener DRexp,3, y a partir de la ec. 26:
Etapa 9
Los valores de CS,0, k0, k1, ..., kj se estiman resolviendo el problema de mínimos cuadrados no lineales:
asumiendo, por ejemplo, que kj es igual a kj-1, tal como se ha mencionado anteriormente.
Por ejemplo, en un tiempo t=t3, los valores CS,0, k0, k1, k2y k3 se calculan resolviendo el problema de mínimos cuadrados no lineales:
Etapa 10
Una vez estimados los valores de CS,0, k0, k1, ..., kj, es posible calcular en un tiempo t=tj la concentración CS,j de humedad residual, usando la ec. 20, o la tasa DRteor,j de desorción.
Etapa 11
El valor calculado de la concentración CS,j de humedad residual o la tasa DRteor,j de desorción se compara con la concentración CS,f de humedad residual final o la tasa DRf de desorción final.
Si el valor estimado de la concentración CS,j de humedad residual o la tasa DRteor,j de desorción es inferior o igual con respecto a la concentración CS,f de humedad residual final o la tasa DRf de desorción final, la fase de secado secundario finaliza.
Si el valor estimado de la concentración CS,j de humedad residual o la tasa DRteor,j de desorción es superior a la concentración CS,f de humedad residual final o la tasa DRf de desorción final, la etapa 7 se repite con t=tj para estimar el tiempo final tf en el que se obtiene la concentración CS,f de humedad residual final (o la tasa DRf de desorción final):
Por ejemplo, en un tiempo t=t3, usando los valores calculados de CS,0 y de las constantes cinéticas, es posible estimar el instante tf de tiempo en el que se obtiene la concentración CS,f de humedad residual final, asumiendo que la temperatura del producto no cambia; debe resolverse la siguiente ecuación no lineal:
Es posible asumir un criterio de detención diferente, p. ej., el requisito de que la tasa de desorción tenga un valor bajo determinado. A continuación, y haciendo referencia a las Figuras 2 a 6, se muestra un ejemplo de la aplicación del método de control de una fase de secado secundario de un proceso de secado de la invención.
La Figura 2 muestra un ensayo experimental que proporciona valores de tasa de desorción con respecto al tiempo durante el secado secundario. Se usa la primera versión del método. Etapa 1 En un tiempo t=t0=0 s, a partir del EAP (y la ec. 4) se obtiene que DRexp,0=0,00056 % agua con respecto a producto
- -
- 1
seco s.
Etapa 2
En un tiempo t=t1=1296 s, a partir del EAP (y la ec. 4) se obtiene que DRexp,1 = 0,00049 % agua s-1.
Etapa 3
En un tiempo t=t2=2592 s, a partir del EAP (y la ec. 4) se obtiene que DRexp,2 = 0,00035 % agua s-1.
Etapa 4
Usando la estimación preliminar de las constantes cinéticas k0y k1 y de CS,0 a partir de la ec. 17 (k0=k1=1,03·10-4 s-1, CS,0=5,48 % agua con respecto a producto seco), se usa la ec. 16 para calcular CS,0 y las constantes cinéticas (CS,0=4,13 % agua con respecto a producto seco).
Etapas 5, 7
Usando los valores calculados de CS,0 y de las constantes cinéticas y la ec 18, es posible estimar el instante tf de tiempo en el que se obtiene el valor deseado de la concentración CS,f de humedad final (p. ej., 0,2 % agua con respecto a producto seco). En este caso, se calcula que siguen siendo necesarios 25056 s.
La Figura 3 muestra una estimación de la evolución temporal de la concentración CS y de la tasa DR de desorción obtenida usando la estimación de CS,0 y de las constantes cinéticas.
En este punto, es posible repetir el procedimiento descrito anteriormente (etapas 7 a 11).
En un tiempo t=t3=3888 s, a partir del EAP (y la ec. 4) se obtiene que DRexp,3 = 0,00028 % agua s-1.
Usando la ec, 16, se calcula que CS,0=4,06 % agua con respecto a producto seco y que siguen siendo necesarios 26352 s.
La Figura 4 muestra la estimación de la evolución temporal de la concentración CS y de la tasa DR de desorción obtenida usando la nueva estimación de CS,0 y de las constantes cinéticas.
Es posible observar en cada repetición que la estimación de los valores de CS,0 mejora, así como la estimación del tiempo tf necesario para finalizar la fase de secado secundario.
La Figura 5 muestra cómo la estimación del tiempo final tf necesario para finalizar la fase de secado secundario cambia con el tiempo. La Figura 6 muestra una comparación entre estimaciones del tiempo final tf necesario para finalizar la fase de secado secundario (puntos extremos de la fase de secado secundario) usando el método de la invención (línea discontinua con puntos redondos) y usando el método según US 6176121 (línea discontinua con puntos cuadrados).
Es posible observar que las estimaciones del tiempo necesario para llegar al final del secado secundario usando el método de la invención son bastante buenas incluso al inicio de la fase y que las mismas se refinan a medida que el secado secundario continúa. En cambio, usando el método descrito en US 6176121, la predicción del tiempo necesario para finalizar el secado secundario no es fiable al inicio y, después de cada EAP, la predicción se actualiza hasta llegar al final del secado.
El método de la invención también ha sido validado mediante una serie de experimentos llevados a cabo en laboratorio.
Las Figuras 7 y 8 son un ejemplo de los resultados que es posible obtener cuando se usa el algoritmo del método.
De forma específica, las Figuras 7 y 8 son una comparación entre los valores experimentales (símbolos) y los previstos por el algoritmo de la invención (línea continua) de la tasa de desorción (Figura 7) y del contenido de agua residual (Figura 8), respectivamente. También se muestra la evolución temporal de la temperatura de almacenamiento (Figura 7, línea en puntos). El tiempo es igual a cero al inicio del secado secundario.
El ejemplo se refiere a un ciclo de secado por congelación de una solución acuosa de sacarosa al 20% en peso (155 viales con un diámetro de 20,85·10-3 m, llenos con 3·10-3 1 de solución). La fase de congelación se llevó a cabo a 50°C durante 17 h, la fase de secado principal se llevó a cabo a -15°C y 10 Pa durante 25 h y la fase de secado secundario se llevó a cabo a 20°C.
5 Los valores experimentales de la tasa de desorción se obtuvieron mediante el ensayo de aumento de presión (ver ec. 4), mientras que el contenido de agua residual se determinó pesando algunos viales extraídos de la cámara de secado usando un captador de muestras.
El modelo cinético para la desorción de agua que se usó con el algoritmo es el mismo de la primera versión del método (ec. 5-18), es decir, se asumió que la tasa de desorción era proporcional al contenido de agua residual.
10 La evolución temporal de la tasa de desorción es una consecuencia del hecho de que, cuando se inicia el secado secundario, la temperatura de almacenamiento aumenta y, durante este intervalo de tiempo, la temperatura del producto y, por lo tanto, la tasa de desorción, aumentan. Después de esto, la temperatura permanece constante y, debido a la disminución del contenido de agua residual, la tasa de desorción disminuye.
Claims (16)
- REIVINDICACIONES1. Método de control de una fase de secado secundario de un proceso de secado por congelación en un aparato secador por congelación que incluye una cámara de secado que contiene un producto a secar y que puede ser aislada para realizar ensayos de aumento de presión, comprendiendo dicho método las etapas de:
- -
- realizar un primer ensayo de aumento de presión en un tiempo t=t0 y calcular un primer valor de tasa (DRexp,0) de desorción experimental de dicho producto (etapa 1);
- -
- realizar un segundo ensayo de aumento de presión en un tiempo t=t1 y calcular un segundo valor de tasa (DRexp,1) de desorción experimental de dicho producto (etapa 2);
- -
- realizar un tercer ensayo de aumento de presión en un tiempo t=t2 y calcular un tercer valor de tasa (DRexp,2) de desorción experimental de dicho producto (etapa 3);
- -
- estimar las condiciones iniciales (CS,0) y las constantes cinéticas (k0, k1, k2) de un modelo cinético del proceso de secado, siendo adecuado dicho modelo cinético para calcular un contenido (CS) de humedad residual y/o una tasa (DRteor) de desorción de dicho producto (etapa 4);
- -
- calcular en un tiempo t=t2 un contenido (CS,2) de humedad residual respectivo y una tasa (DRteor,2) de desorción respectiva (etapa 5).
- 2. Método según la reivindicación 1, que comprende además, después de la etapa 5, las etapas de:
- -
- comparar dicho contenido (CS,2) de humedad residual y/o dicha tasa (DRteor,2) de desorción calculados en un tiempo t=t2 con una concentración (CS,f) de humedad residual final deseada y/o una tasa (DRf) de desorción final deseada, respectivamente, (etapa 6); si dicho contenido (CS,2) de humedad residual es inferior o igual con respecto a dicha concentración (Cs,f) de humedad residual final o dicha tasa (DRteor,2) de desorción es inferior o igual con respecto a dicha tasa (DRf) de desorción final, considerar la fase de secado secundario finalizada; en caso contrario, el método comprende además las etapas de:
- -
- estimar un tiempo final (tf) en el que se obtiene dicha concentración (CS,f) de humedad residual final o dicha tasa (DRf) de desorción final (etapa 7);
- -
- realizar un ensayo de aumento de presión adicional en un tiempo t=tj y calcular en dicho tiempo t=tj un contenido (CS,j) de humedad residual respectivo y una tasa (DRteor,j) de desorción respectiva (etapa 8);
- -
- estimar las condiciones iniciales (CS,0) y las constantes cinéticas (k0, k1, k2, ..., kj) de dicho modelo cinético (etapa 9);
- -
- calcular en dicho tiempo t=tj un contenido (CS,j) de humedad residual respectivo y/o una tasa (DRteor,j)de desorción respectiva (etapa 10);
- -
- comparar dicho contenido (CS,j) de humedad residual y/o dicha tasa (DRteor,j) de desorción calculados en dicho tiempo t=tj con dicha concentración (CS,f) de humedad residual final y/o dicha tasa (DRf) de desorción final, respectivamente, (etapa 11); si dicho contenido (CS,j) de humedad residual es inferior o igual con respecto a dicha concentración (CS,f) de humedad residual final o dicha tasa (DRteor,j) de desorción es inferior o igual con respecto a dicha tasa (DRf) de desorción final, considerar la fase de secado secundario finalizada; en caso contrario, repetir las etapas 7 a 11.
- 3. Método según la reivindicación 1 o 2, en el que dicha tasas (DRexp,0, DRexp,1, DRexp,2) de desorción experimentales se calculan usando la ecuación:donde: DRexp: tasa de desorción experimental, [% agua y/o disolvente s-1]P: presión medida, [Pa]t: tiempo, [s]t0: instante de tiempo al inicio del ensayo de aumento de presión, [s]R: constante de gas [8, 314 J mol-1 K-1]T: temperatura del vapor, [K]V: volumen (libre) de la cámara de secado, [m3] Mw: peso molecular de agua y/o disolvente, [kg mol-1] mseca: masa del producto seco, [kg]
-
- 4.
- Método según cualquiera de las reivindicaciones anteriores, en el que dicho modelo cinético comprende ecuaciones matemáticas adecuadas para modelizar la dependencia de la tasa (DR) de desorción con respecto al contenido (CS) de humedad residual en el producto.
-
- 5.
- Método según cualquiera de las reivindicaciones anteriores, en el que se asume que dicha tasa (DRteor)de desorción depende de dicho contenido (CS) de humedad residual en dicho producto según la ecuación:
donde:DR: tasa de desorción, [% agua y/o disolvente s-1] K : constante cinética del proceso, [s-1] CS: contenido de humedad residual, [% agua/disolvente con respecto a producto seco] - 6. Método según la reivindicación 5, en el que una evolución temporal de dicha concentración (CS) de humedad residual en un tiempo t=tj viene dada por la integración de la siguiente ecuación diferencial:donde: DRj: tasa de desorción en un tiempo t=tj, [% agua y/o disolvente s-1]t : tiempo, [s] kj : constante cinética del proceso en un tiempo t=tj, [s-1].
- 7. Método según la reivindicación 6, en el que dicho cálculo de un contenido (Cs) de humedad residual se realiza mediante la ecuación:donde:CS,0 : valor de la humedad residual [% agua y/o disolvente con respecto a producto seco] al inicio de la fase de secado secundario (t=t0); kr : constante cinética del proceso en un tiempo t=tr (con r=1, 2, ..., j), [s-1].
-
- 8.
- Método según la reivindicación 7, en el que dicho cálculo de una tasa (DRteor) de desorción se realiza mediante la ecuación:
-
- 9.
- Método según la reivindicación 8, dependiendo la reivindicación 5 de la reivindicación 3, en el que dicha estimación de las condiciones iniciales (CS,0) y las constantes cinéticas (k0, k1, k2, ..., kj) en un tiempo t=tj se realiza mediante las siguientes ecuaciones:
y resolviendo el siguiente problema de mínimos cuadrados no lineales: - 10. Método según cualquiera de las reivindicaciones 7 a 9, dependiendo la reivindicación 5 de cualquiera de las reivindicaciones 2 a 4 y dependiendo las reivindicaciones 3 y 4 de la reivindicación 2, en el que dicho tiempo final (tf) se calcula, asumiendo que la temperatura de dicho producto no cambia, mediante la siguiente ecuación, resultante de (ec. 11):donde:CS,f: concentración de humedad residual final [% agua y/o disolvente con respecto a producto seco];CS,j: concentración de humedad residual en un tiempo t=tj [% agua y/o disolvente con respecto a productoseco];
- 11. Método según cualquiera de las reivindicaciones 1 a 4, en el que se asume que dicha tasa (DRteor) de desorción depende de dicho contenido (CS) de humedad residual en dicho producto según la ecuación:DR: tasa de desorción, [% agua y/o disolvente s-1] k : constante cinética del proceso, [s-1] CS: concentración de humedad residual, [% agua y/o disolvente con respecto a producto seco]Cs,eq: concentración de humedad de equilibrio, [% agua y/o disolvente con respecto a producto seco]
- 12. Método según la reivindicación 11, en el que una evolución temporal de dicha concentración (CS) de humedad residual en un tiempo t=tj viene dada por la integración de la siguiente ecuación diferencial:donde:DRj: tasa de desorción en un tiempo t=tj, [% agua y/o disolvente s-1]t : tiempo, [s] kj : constante cinética del proceso, [s-1]. Cs,eq,j: concentración de humedad de equilibrio en un tiempo t=tj, [% agua y/o disolvente con respecto aproducto seco]
- 13. Método según la reivindicación 12, en el que dicho cálculo de un contenido (CS) de humedad residual en un tiempo t=tj se realiza mediante las siguientes ecuaciones:donde:CS,0 : valor de la humedad residual [% agua y/o disolvente con respecto a producto seco] al inicio de la fase de secado secundario (t=t0); kr : constante cinética del proceso en un tiempo t=tr (con r=1, 2, ..., j), [s-1]; Cs,eq,r: concentración de humedad de equilibrio en un tiempo t=tr con r=1, 2, ..., j), [% agua y/o disolventecon respecto a producto seco].
-
- 14.
- Método según la reivindicación 13, en el que dicho cálculo de una tasa (DRteor) de desorción se realiza mediante la ecuación:
-
- 15.
- Método según la reivindicación 14, dependiendo la reivindicación 11 de la reivindicación 3, en el que dicha estimación de las condiciones iniciales (CS,0) y las constantes cinéticas (k0, k1, k2, ..., kj) en un tiempo t=tj se realiza mediante las siguientes ecuaciones:
-
- 16.
- Método según cualquiera de las reivindicaciones 11 a 15, dependiendo la reivindicación 11 de cualquiera de las reivindicaciones 2 a 4 y dependiendo las reivindicaciones 3 y 4 de la reivindicación 2, en el que dicho tiempo final (tf) se calcula, asumiendo que la temperatura de dicho producto no cambia, mediante la siguiente ecuación, resultante de (ec. 21):
y resolviendo el siguiente problema de mínimos cuadrados no lineales:donde:CS,f: concentración de humedad residual final [% agua y/o disolvente con respecto a producto seco];CS,j: concentración de humedad residual en un tiempo t=tj [% agua y/o disolvente con respecto a productoseco].
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| EP08013243A EP2148158B1 (en) | 2008-07-23 | 2008-07-23 | Method for monitoring the secondary drying in a freeze-drying process |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| ES2376675T3 true ES2376675T3 (es) | 2012-03-15 |
Family
ID=40232972
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| ES08013243T Active ES2376675T3 (es) | 2008-07-23 | 2008-07-23 | Método de control del secado secundario en un proceso de secado por congelación. |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US9879909B2 (es) |
| EP (1) | EP2148158B1 (es) |
| CN (1) | CN101634845B (es) |
| AT (1) | ATE532016T1 (es) |
| ES (1) | ES2376675T3 (es) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6503894B1 (en) * | 2000-08-30 | 2003-01-07 | Unimed Pharmaceuticals, Inc. | Pharmaceutical composition and method for treating hypogonadism |
| EP1870649A1 (en) * | 2006-06-20 | 2007-12-26 | Octapharma AG | Lyophilisation targetting defined residual moisture by limited desorption energy levels |
| WO2008042408A2 (en) * | 2006-10-03 | 2008-04-10 | Wyeth | Lyophilization methods and apparatuses |
| IT1397930B1 (it) * | 2009-12-23 | 2013-02-04 | Telstar Technologies S L | Metodo per monitorare l'essiccamento primario di un processo di liofilizzazione. |
| US9459044B1 (en) | 2013-03-15 | 2016-10-04 | Harvest Right, LLC | Freeze drying methods and apparatuses |
| US11209391B2 (en) * | 2016-09-08 | 2021-12-28 | Atonarp Inc. | System having a pre-separation unit |
| CN106853417B (zh) * | 2016-11-18 | 2019-02-26 | 中核兰州铀浓缩有限公司 | 离心级联小量离心机装架真空干燥方法 |
| US11359861B2 (en) | 2018-04-10 | 2022-06-14 | Ima Life North America Inc. | Freeze drying process and equipment health monitoring |
| US11744257B1 (en) | 2018-10-19 | 2023-09-05 | Harvest Right, LLC | Freeze-drying methods including vacuum freezing |
| US11287185B1 (en) | 2020-09-09 | 2022-03-29 | Stay Fresh Technology, LLC | Freeze drying with constant-pressure and constant-temperature phases |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| ES2125652T3 (es) | 1995-02-14 | 1999-03-01 | Georg Wilhelm Oetjen | Procedimiento para la determinacion de la humedad residual durante la desecacion secundaria en un proceso de liofilizacion. |
| US6971187B1 (en) | 2002-07-18 | 2005-12-06 | University Of Connecticut | Automated process control using manometric temperature measurement |
| EP1903291A1 (en) * | 2006-09-19 | 2008-03-26 | Ima-Telstar S.L. | Method and system for controlling a freeze drying process |
-
2008
- 2008-07-23 ES ES08013243T patent/ES2376675T3/es active Active
- 2008-07-23 EP EP08013243A patent/EP2148158B1/en active Active
- 2008-07-23 AT AT08013243T patent/ATE532016T1/de active
-
2009
- 2009-07-14 US US12/502,863 patent/US9879909B2/en active Active
- 2009-07-23 CN CN200910165151.1A patent/CN101634845B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN101634845B (zh) | 2014-05-14 |
| US9879909B2 (en) | 2018-01-30 |
| EP2148158B1 (en) | 2011-11-02 |
| EP2148158A1 (en) | 2010-01-27 |
| US20100018073A1 (en) | 2010-01-28 |
| ATE532016T1 (de) | 2011-11-15 |
| CN101634845A (zh) | 2010-01-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| ES2376675T3 (es) | Método de control del secado secundario en un proceso de secado por congelación. | |
| ES2387071T3 (es) | Método para controlar un proceso de secado por congelación | |
| ES2300947T3 (es) | Dispositivo y procedimiento de guiado de la operacion de deshidratacion durante un tratamiento de liofilizacion. | |
| Patel et al. | Determination of end point of primary drying in freeze-drying process control | |
| Gatlin et al. | Freeze drying: a practical overview | |
| Fissore et al. | Monitoring of the secondary drying in freeze-drying of pharmaceuticals | |
| US6971187B1 (en) | Automated process control using manometric temperature measurement | |
| ES2471016T3 (es) | Método de control de secado principal de un proceso de secado por congelación | |
| KR20190035808A (ko) | 냉동-건조를 위한 방법 및 장치 및 컨테이너 | |
| EP2069699A2 (en) | Lyophilization methods and apparatuses | |
| US11243029B2 (en) | Process monitoring and control for lyophilization using a wireless sensor network | |
| ES2833378T3 (es) | Uso estéril repetido de un medidor en un sistema de secado por congelación esterilizable con vapor | |
| Fissore et al. | On the use of temperature measurement to monitor a freeze-drying process for pharmaceuticals | |
| ES2779023T3 (es) | Procedimiento para la determinación a base de presión de un parámetro de producto en un liofilizador, liofilizador y producto de software | |
| Kuu et al. | Rapid freeze-drying cycle optimization using computer programs developed based on heat and mass transfer models and facilitated by tunable diode laser absorption spectroscopy (TDLAS) | |
| JP3850442B2 (ja) | 凍結乾燥工程での後乾燥中の残留湿分の測定方法 | |
| ES2626001T3 (es) | Optimización de la nucleación y cristalización en la liofilización utilizando una congelación con intersticios | |
| ES2734148T3 (es) | Dispositivo para regular la temperatura de los productos que se van a liofilizar | |
| JP7232097B2 (ja) | 凍結乾燥方法及び凍結乾燥装置 | |
| EP4105585A1 (en) | Freeze-drying method and apparatus | |
| CN114353440B (zh) | 用于冻干机的控制方法及装置、冻干机 | |
| JP2016125682A (ja) | 凍結乾燥機に適用される被乾燥材料の乾燥状態監視装置及び乾燥状態監視方法 | |
| Thompson et al. | Use of a micro freeze-dryer for developing a freeze-drying process | |
| Heldner | Pharmaceutical Freeze‐Drying Systems | |
| ES2787923T3 (es) | Un vehículo que comprende varios compartimentos de carga refrigerados |