# Patent

## Title
Method for monitoring the secondary drying in a freeze-drying process

## Abstract

A method to monitor a secondary drying phase of a freeze-drying process comprises initial steps in which is provided to perform pressure rise tests at different time and to calculate a respective value of experimental desorption rate of product (steps 1 to 3). Subsequently, the method provides to estimate initial conditions and kinetic constants of a kinetic model of the process (step 4) and to calculate at time t=t2 a respective residual moisture content and a respective desorption rate (step 5). The method can be performed in a freeze-dryer apparatus which includes a drying chamber that contains a product to be dried and can be isolated to perform pressure rise tests.

## Background
CROSS REFERENCE TO RELATED APPLICATION
This application is a new U.S. utility application claiming priority benefit of EP 08013243.4, filed Jul. 23, 2008, the entire contents of which are hereby incorporated by reference. The invention relates to methods for monitoring a freeze-drying process in a freeze-dryer; in particular it refers to a method for monitoring secondary drying of a freeze-drying process, for example, of pharmaceutical products arranged in containers.
BACKGROUND OF THE INVENTION
Freeze-drying, also known as lyophilization, is a dehydration process that enables removal by sublimation of water and/or solvents from a substance, such as food, pharmaceutical or biological products. Typically the freeze-drying process is used to preserve a perishable product since the greatly reduced water content that results inhibits the action of microorganisms and enzymes that would normally spoil or degrade the product. Furthermore, the process makes the product more convenient for transport. Freeze-dried products can be sealed in containers to prevent the reabsorption of moisture and can be easily rehydrated or reconstituted by addition of removed water and/or solvents. In this way the product may be stored at room temperature without refrigeration, and be protected against spoilage for many years.
Since freeze-drying is a low temperature process in which the temperature of product does not exceed typically 30° C. during the operating phases, it causes less damage or degradation to the product than other dehydration processes using higher temperatures. Freeze-drying does not usually cause significant shrinkage or toughening of the product being dried. Freeze-dried products can be rehydrated much more quickly and easily because of the porous structure created during the sublimation of ice.
In the pharmaceutical field, freeze-drying process is widely used in the production of pharmaceuticals, mainly for parenteral and oral administration, also because freeze-drying process can be carried out in sterile conditions.
A known freeze-dryer apparatus for performing a freeze-drying process usually comprises a drying chamber and a condenser chamber interconnected by a duct that is provided with a valve that allows isolating the drying chamber when required during the process.
FIG. 1 shows the drying chamber which comprises a plurality of temperature-controlled shelves arranged for receiving containers of product to the dried. The condenser chamber includes condenser plates or coils having surfaces maintained at very low temperature, e.g. −50° C., by means of a refrigerant of freezing device. The condenser chamber is also connected to one or more vacuum pumps so as to achieve high vacuum values inside both chambers.
Freeze-drying process typically comprises three phases: a freezing phase, a primary drying phase and a secondary drying phase.
During the freezing phase the shelf temperature is reduced up to typically −30/−40° C. in order to convert into ice most of the water and/or solvents contained in the product.
In the primary drying phase the shelf temperature is increased, while the pressure inside the drying chamber is lowered below 1-5 mbar so as to allow the frozen water and/or solvents in the product to sublime directly from solid phase to gas phase. The application of high vacuum makes possible the water sublimation at low temperatures.
Heat is supplied to the product and the vapour generated by sublimation of frozen water and/or solvents is removed from the drying chamber by means of condenser plates or coils of condenser chamber wherein the vapour can be re-solidified.
Secondary drying phase is provided for removing by desorption the residual moisture of the product, namely the amount of unfrozen water and/or solvents that cannot be removed during primary drying when sublimation of ice takes place. During this phase the shelf temperature is further increased up to a maximum of 30-60° C. to heat the product, while the pressure inside the drying chamber is set typically below 0.1 mbar.
At the end of secondary drying phase the product is sufficiently dried with residual moisture content typically of 1-3%.
Secondary drying has to be carefully monitored in order to point out when the drying process is completed, i.e. when the desired amount of residual moisture in the product has been achieved.
There are known methods for monitoring secondary drying phase.
According to a known method the residual moisture of the product can be determined by extracting samples from the freeze-dryer without interrupting the freeze-drying (e.g. using a “sample thief”) and measuring off-line their moisture content by means of Karl Fischer titration, thermal gravimetric analysis, or near Infra-Red spectroscopy.
U.S. Pat. No. 6,971,187 proposes another method wherein the estimation of the drying rate of the product during the secondary drying is obtained by performing a Pressure Rise Test (PRT).
During a PRT the drying chamber is isolated from the condenser chamber by closing the valve positioned in the duct connecting the two chambers. As the heating is not stopped, the ice sublimation continues, thus increasing in the drying chamber the pressure that can be measured.
Given the curve of pressure vs. time, the slope at the beginning of this curve allows estimating the flow rate of water and/or solvent from the product by the equation:




  
    
      
        
          
            
              d
              ⁢
              
                
              
              ⁢
              P
            
            
              d
              ⁢
              
                
              
              ⁢
              t
            
          
          ⁢
          
            |
            
              t
              =
              
                t
                0
              
            
          
        
        =
        
          
            RT
            V
          
          ⁢
          
            j
            
              w
              ,
              n
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          1
        
        )
      
    
  




where:
P: measured pressure, [Pa]
t: time, [s]
t0: time instant at the beginning of the PRT, [s]
R: gas constant [8.314 J mol−1 K−1]
T: temperature of the vapour, [K]
V: (free) volume of the chamber, [m3]
jw,n: flow rate of water and/or solvent from the product, [mol s−1]
Thus, the mass flow of water and/or solvent can be calculated:




  
    
      
        
          j
          
            w
            ,
            m
          
        
        =
        
          
            
              M
              w
            
            ⁢
            
              V
              RT
            
            ⁢
            
              
                d
                ⁢
                
                  
                
                ⁢
                P
              
              
                d
                ⁢
                
                  
                
                ⁢
                t
              
            
          
          ⁢
          
            |
            
              t
              =
              
                t
                0
              
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          2
        
        )
      
    
  




where:
jw,m: mass flow of water and/or solvent from the product, [kg s−1]
Mw: molecular weight of water and/or solvent, [kg mol−1]
From this value, the loss in water and/or solvent during the measurement period elapsed between two consecutive PRTs can be estimated by:

Δwm,j=jw,m,j−1Δtj  (eq. 3)

where:

    
    
        tj=tj−tj−1 time elapsed between j-th PRT and (j−1)th PRT, [s]
        wm,j: loss in water during the time interval tj, [kg]
    
    


jw,m,j−1: mass flow of water and/or solvent from the product calculated from the (j−1)-th PRT, [kg s−1].
The total amount of water and/or solvent removed between a reference time t0 (e.g. the start of the secondary drying) and any given time of interest tj is simply the summation of all the wm,j occurring in the various intervals between PRTs. Exploiting one independent experimental value for detecting the residual water content at a reference time, e.g. at the end of primary drying, the real time actual moisture content vs. time can be calculated. This requires extracting a sample from the drying chamber or using expensive sensors (e.g. NIR-based sensors) to get this value in-line.
Given this experimental value, some empirical or common sense indications are given to calculate the “optimal” temperature to minimize the time required to complete the secondary drying.
A disadvantage of the above known methods consists in that they require extracting samples from the drying chamber and using expensive sensors for measuring the experimental values of residual water and/or solvent. Samples extraction is an invasive operation that perturbs the freeze-drying process and thus it is not suitable in sterile and/or aseptic processes and/or when automatic loading/unloading of the containers is used. Furthermore, sample extraction is time consuming and requires skilled operators.
Another disadvantage of the method disclosed in U.S. Pat. No. 6,971,187 is that the empirical and common sense indications used for calculating the “optimal” temperature do not allow to optimize the process.
A different approach is disclosed in U.S. Pat. No. 6,176,121 wherein using two successive measurements of desorption rate (DR), i.e. the mass flow rate of the water and/or solvent vapour due to desorption, calculated from jw,m, it is possible to extrapolate the point in time at which a given small value of DR is obtained. In order to do this, the valve placed between the drying chamber and the condenser chamber should be regularly closed for a certain time and the pressure rise curve (PRC), caused by the desorbing water vapour, has to be acquired. Thus, the mass of desorbed water and/or solvent over the time, or rather the desorption rate, can be calculated from the initial slope of the PRC as follows:




  
    
      
        
          DR
          exp
        
        =
        
          
            
              VM
              w
            
            RT
          
          ⁢
          
            
              (
              
                
                  d
                  ⁢
                  
                    
                  
                  ⁢
                  P
                
                
                  d
                  ⁢
                  
                    
                  
                  ⁢
                  t
                
              
              )
            
            
              t
              =
              
                t
                0
              
            
          
          ⁢
          
            100
            
              m
              dried
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          4
        
        )
      
    
  




where:
mdried: mass of the dried product, [kg]
DRexp: experimental desorption rate, [% of water and/or solvent over dried product s−1]
A disadvantage of this method consists in that, due to the very simplified approach, it is shown to fail in correspondence of the end of secondary drying. Moreover, it does not allow to estimate the absolute residual moisture, but only the difference with respect to the equilibrium moisture, which depends on the operating conditions (shelf temperature and drying chamber pressure), and therefore no target about this value can be set.
An object of the invention is to improve the methods for monitoring a freeze-drying process in a freeze-dryer, particularly for monitoring a secondary drying phase of said freeze-drying process.
A further object is to provide a method for calculating process parameters, such as residual moisture content and/or desorption rate of a dried product, that is non-invasive and not-perturbing the freeze-drying process and thus is suitable for being used in sterile and/or aseptic processes and/or when automatic loading/unloading of the containers is used.
Another object is to provide a method capable to precisely estimate initial conditions and kinetic constants of a kinetic model of the drying process, suitable for calculating the process parameters.
Still another object is to provide a method for estimating in a reliable and precise way a residual moisture concentration and/or desorption rate of the dried product during secondary drying phase and a time required for terminating said secondary drying phase.
Another further object is to provide a method wherein estimation of process parameters is progressively improved and refined during progress of secondary drying phase, said estimation being nevertheless good with respect to known methods even at the beginning of secondary drying phase.



BRIEF DESCRIPTION OF THE DRAWINGS
The invention can be better understood and carried into effect with reference to the enclosed drawings, that show an embodiment of the invention by way of non-limitative example, in which
FIG. 1 is a freeze-dryer apparatus for performing a freeze-drying process;
FIG. 2 is a flowchart schematically showing the method of the invention for monitoring a secondary drying phase in a freeze-drying process;
FIG. 3 is a graph showing a sequence of experimental measured values of desorption rate vs time during secondary drying;
FIG. 4 are graphs showing estimation of time evolution respectively of residual moisture concentration and desorption rate of dried product at a further defined time;
FIG. 5 are graphs showing estimation of time evolution respectively of residual moisture concentration and desorption rate of dried product at a further defined time;
FIG. 6 is a graph showing a time evolution sequence of estimations of time required to complete secondary drying;
FIG. 7 illustrates a comparison between estimations of time required to complete secondary drying obtained using the method of the invention and using the method according to U.S. Pat. No. 6,176,121;
FIGS. 8 and 9 show a comparison between experimental values and values predicted by the method of the invention respectively of the desorption rate and of the residual water content.



DESCRIPTION OF THE INVENTION
According to the invention, a method is provided for monitoring a secondary drying phase of a freeze-drying process in a freeze-dryer apparatus including a drying chamber that contains a product to be dried and can be isolated for performing pressure rise tests, said method comprising the steps of:

    
    
        performing a first pressure rise test at time t=t0 and calculating a first value of experimental desorption rate of said product (step 1);
        performing a second pressure rise test at time t=t1 and calculating a second value of experimental desorption rate of said product (step 2);
        performing a third pressure rise test at time t=t2 and calculate a third value of experimental desorption rate of said product (step 3);
        estimating initial conditions and kinetic constants of a kinetic model of the drying process, said kinetic model being suitable for calculating a residual moisture content and/or a desorption rate of said product (step 4);
        calculating at time t=t2 a respective residual moisture content and a respective desorption rate (step 5).
    
    


The method further comprises, after step 5, the step of:

    
    
        comparing said residual moisture content and/or said desorption rate calculated at time t=t2 respectively with a desired final residual moisture concentration and/or a desired final desorption rate (step 6); if said residual moisture content is lower than, or equal to, said final residual moisture concentration or said desorption rate is lower than, or equal to, said final desorption rate, then the secondary drying phase is considered ended; if not the method further comprising the steps of:
        estimating a final time at which said final residual moisture concentration or said final desorption rate is obtained (step 7);
        performing a further pressure rise at time t=tj and calculating at said time t=tj a respective residual moisture content and a respective desorption rate (step 8);
        estimating initial conditions and kinetic constants of said kinetic model (step 9);
        calculating at said time t=tj said respective residual moisture content and/or said respective desorption rate (step 10);
        comparing said residual moisture content and/or said desorption rate calculated at said time t=tj respectively with said final residual moisture concentration and/or said final desorption rate (step 11); if said residual moisture content is lower than, or equal to, said final residual moisture concentration or said desorption rate is lower than, or equal to, said final desorption rate then the secondary drying phase is considered ended; if not step 7 to 11 are repeated.
    
    


Owing to the invention it is possible to obtain a method for calculating in a reliable and precise way the residual moisture concentration and/or desorption rate of a dried product during a secondary drying phase of a freeze-drying process. The method is also capable to precisely estimate initial conditions and kinetic constants of a kinetic model of the drying process, which calculates the residual moisture concentration and/or desorption rate process, without extracting any samples from the drying chamber and without using expensive sensors to get this value in-line. Thus, the monitoring method of the invention is non-invasive and non-perturbing the freeze-drying process and is suitable for being used in sterile and/or aseptic processes and/or when automatic loading/unloading of the containers is used.
Furthermore, the method allows calculating the time required for terminating said secondary drying phase, wherein the stop requirement can be that the residual moisture concentration, or the desorption rate, has a respective desired final value.
Since the steps of the method are iterated till the end of secondary drying phase is reached, estimation of process parameters is progressively improved and refined during progress of secondary drying phase, said estimation being nevertheless good with respect to known methods even at the beginning of secondary drying phase.
The method of the invention monitors a secondary drying phase of a freeze-drying process in a freeze-dryer. In particular, the method calculates the residual moisture content of a dried product and provides a reliable estimation of the time that is necessary to complete this phase, according to the desired target (final moisture content and/or final value of desorption rate).
The method requires performing periodically a Pressure Rise Test (PRT) and thus can be applied to those freeze-drying processes that are carried out in freeze-dryers comprising a drying chamber, where the product to be dried is placed, and a separate condenser chamber, where the vapour generated by drying process flow and can be re-solidified or frozen.
The PRT is carried out by closing for a short time interval (from few tens of seconds, e.g. 30 s, to few minutes) a valve that is placed on the duct that connects drying chamber to condenser chamber and measuring (and recording) the time evolution of the total pressure in the chamber.
From the slope of the curve at the beginning of the test the current water and/or solvent desorption rate (DR, % s−1) can be calculated. The PRT is repeated every pre-specified time interval (e.g. 30 minutes) in order to know the time evolution of the water and/or solvent desorption rate. The time interval can be constant or can be changed during the operation.
All the methods based on the PRT for monitoring the primary drying step of a freeze-drying process take advantage from the fact that, during the test, the pressure in the drying chamber increases until equilibrium is reached. As this is not the case for secondary drying (due to the low values of the flow rate of water and/or solvent), the only information that can be exploited from PRT is the estimation of the water and/or solvent flow rate, that can thus be integrated in order to evaluate the water and/or solvent loss in time. The estimation of the moisture content requires knowing the initial moisture concentration, which is calculated according to the method of the invention, as described in detail in the following, without extracting any samples from the drying chamber and without using expensive sensors to get this value in-line. In other words, the monitoring method is non-invasive and non-perturbing the freeze-drying process and thus is suitable for being used in sterile and/or aseptic processes and/or when automatic loading/unloading of the containers is used.
The method of the invention requires modelling the dependence of the Desorption Rate (DR) on the residual moisture content (CS) in the dried product. Various known mathematical equations can be used to this purpose. The method comprises an algorithm able to work efficiently whichever correlation is used.
Various kinetic models have been proposed to model the desorption rate of water and/or solvent. The desorption rate can be assumed to depend on the residual moisture content, or on the difference between the residual moisture content and the equilibrium value.
Both types of models have been demonstrated to perform more or less in the same way; moreover, there is uncertainty about the real physical mechanism of water and/or solvent desorption that may depend on the product considered.
In a first version of the method, the desorption rate DR is assumed to depend on the residual moisture CS in the solid matrix of the dried product, according to equation:

DR=−kCS  (eq. 5)

the time evolution of the residual moisture Cs, given in % of water and/or solvent per dried mass, can be calculated by the integration of the following differential equation:




  
    
      
        
          
            d
            ⁢
            
              
            
            ⁢
            
              C
              S
            
          
          
            d
            ⁢
            
              
            
            ⁢
            t
          
        
        =
        
          DR
          =
          
            -
            
              kC
              S
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          6
        
        )
      
    
  




where t is the time [s] and k is the kinetic constant of the process [s−1].
The kinetic constant can be a function of the temperature and, thus, it can change with time as the temperature of the product can change with time, in particular at the beginning of the secondary drying when the temperature is risen from the value used during primary drying to that of the secondary drying.
If a PRT is made at time t=tj−1 and the successive PRT is made at time t=tj and the product temperature, that is slightly varying in the interval [tj−tj−1], is assumed to be constant and equal to a mean value, the variation of the moisture concentration in the solid can be described by the equation:




  
    
      
        
          
            d
            ⁢
            
              
            
            ⁢
            
              C
              S
            
          
          
            d
            ⁢
            
              
            
            ⁢
            t
          
        
        =
        
          
            DR
            j
          
          =
          
            
              -
              
                k
                j
              
            
            ⁢
            
              C
              S
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          7
        
        )
      
    
  




The solution of eq. 7 requires the initial condition, i.e. the value of the residual moisture CS at time t=tj−1:

CS=CSj−1e−kj(t−tj−1)  (eq. 8)

The value of CS,j−1 can be calculated from the time integration of eq. 6 in the previous time interval:

CSj−1=CS,j−2e−kj−1(tj−1−tj−2)  (eq. 9)

and thus:

CS=CS,j−2e−kj−1(tj−1−tj−2)e−kj(t−tj−1)  (eq. 10)

This procedure can be iterated until the value CS,0 of the residual moisture at the beginning of the secondary drying phase (t=t0) appears. Thus, in the time interval between tj and tj−1 the evolution of the residual moisture concentration is given by:




  
    
      
        
          C
          S
        
        =
        
          
            C
            
              S
              ,
              0
            
          
          ⁢
          
            
              ∏
              
                i
                =
                1
              
              
                j
                -
                1
              
            
            ⁢
            
              
            
            ⁢
            
              
                e
                
                  -
                  
                    
                      k
                      i
                    
                    ⁡
                    
                      (
                      
                        
                          t
                          i
                        
                        -
                        
                          t
                          
                            i
                            -
                            1
                          
                        
                      
                      )
                    
                  
                
              
              ⁢
              
                e
                
                  -
                  
                    
                      k
                      j
                    
                    ⁡
                    
                      (
                      
                        t
                        -
                        
                          t
                          
                            j
                            -
                            1
                          
                        
                      
                      )
                    
                  
                
              
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          11
        
        )
      
    
  




The solution of eq. 11 requires the value of initial moisture concentration CS,0.
The evolution of the theoretical value of the desorption rate in the time interval between tj and tj−1 is thus given by:




  
    
      
        
          DR
          theor
        
        =
        
          
            -
            
              k
              j
            
          
          ⁢
          
            C
            
              S
              ,
              0
            
          
          ⁢
          
            
              ∏
              
                i
                =
                1
              
              
                j
                -
                1
              
            
            ⁢
            
              
            
            ⁢
            
              
                e
                
                  -
                  
                    
                      k
                      i
                    
                    ⁡
                    
                      (
                      
                        
                          t
                          i
                        
                        -
                        
                          t
                          
                            i
                            -
                            1
                          
                        
                      
                      )
                    
                  
                
              
              ⁢
              
                e
                
                  -
                  
                    
                      k
                      j
                    
                    ⁡
                    
                      (
                      
                        t
                        -
                        
                          t
                          
                            j
                            -
                            1
                          
                        
                      
                      )
                    
                  
                
              
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          12
        
        )
      
    
  




If CS,0 and the values of the various kj are perfectly known and the model given by eq. 6 is adequate to describe the dynamics of the system, eq. 11 can be used to know the time evolution of the residual moisture content and thus the time that is required to fulfill the requirements on the final value of the moisture content in the product. If the requirement is on the value of the desorption rate, eq. 12 can be used to this purpose.
The above situation is quite rare, since the value of initial moisture concentration has to be measured by extracting samples and the various kinetic constants are never known a priori.
The method according to the invention provides calculating initial condition Cs, 0 and kinetic constants performing the following steps as shown in the flowchart of FIG. 2.
Step 1
At time t=t0 a PRT is performed and a respective desorption rate DR (indicated in the following as DRexp,0) is calculated, i.e. using eq. 4.
From eq. 12 it is:

DRexp,0=DRtheor,0=−k0CS,0  (eq. 13)

Step 2
At time t=t1 a PRT is performed and a respective desorption rate DR (indicated in the following as DRexp,0) is calculated, i.e. using eq. 4.
From eq. 12 it is:

DRexp,1=DRtheor,1=−k1CS,0e−k1(t1−t0)  (eq. 14)

Step 3
At time t=t2 a PRT is performed and the desorption rate DR (indicated in the following ad DRexp,2) is calculated, i.e. using eq. 4.
From eq. 12 it is:

DRexp,2=DRtheor,2=−k2Cs,0e−k1(t1−t0)e−k2(t2−t1)  (eq. 15)

Step 4
Values of CS,0, k0, k1 and k2 are estimated so that the calculated values of the desorption rate matches with all the experimental values available (DRexp,0, DRexp,1 and DRexp,2). This can be done using a minimization algorithm to solve the following non-linear least-square problem:




  
    
      
        
          min
          
            
              C
              
                S
                ,
                0
              
            
            ,
            
              k
              i
            
          
        
        ⁢
        
          
            ∑
            
              i
              =
              0
            
            2
          
          ⁢
          
            
              (
              
                
                  DR
                  
                    exp
                    ,
                    i
                  
                
                -
                
                  DR
                  
                    theor
                    ,
                    i
                  
                
              
              )
            
            2
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          16
        
        )
      
    
  




and assuming, for example, that k2 is equal to k1, due to fact that the time interval between two PRTs is generally small, e.g. 30 minutes, and to fact that the temperature of the product is almost constant during secondary drying (only at the beginning of secondary drying the temperature of the product varies, from that of primary drying to that required by secondary drying, but this variation is generally slow, due to the thermal inertia of the system).
As starting values for eq. 16 it is possible to use the rough approximations of k0, k1 and CS,0 that can be calculated from eq. 13 and eq. 14 after the first two PRTs:




  
    
      
        
          
            k
            0
          
          =
          
            
              k
              1
            
            =
            
              
                -
                
                  1
                  
                    
                      t
                      1
                    
                    -
                    
                      t
                      0
                    
                  
                
              
              ⁢
              ln
              ⁢
              
                
                  DR
                  
                    exp
                    ,
                    1
                  
                
                
                  DR
                  
                    exp
                    ,
                    0
                  
                
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            C
            
              S
              ,
              0
            
          
          =
          
            
              DR
              
                exp
                ,
                0
              
            
            
              
                1
                
                  
                    t
                    1
                  
                  -
                  
                    t
                    0
                  
                
              
              ⁢
              ln
              ⁢
              
                
                  DR
                  
                    exp
                    ,
                    1
                  
                
                
                  DR
                  
                    exp
                    ,
                    0
                  
                
              
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          17
        
        )
      
    
  




These values are just a first approximation of the kinetic constants k0 and k1 and of residual moisture content CS,0; these estimations will be refined after each PRT.
Step 5
Once estimated the values of CS,0, k0, k1 and k2, at time t=t2 it is calculated the residual moisture concentration CS,2, using eq. 11, or the desorption rate DRtheor,2, using eq. 12.
Step 6
The calculated residual moisture concentration CS,2, or desorption rate DRtheor,2, is compared with a desired value of final or target residual moisture concentration Cs,f, or a desired value of final or target desorption rate DRf.
If the calculated residual moisture concentration CS,2, or desorption rate DRtheor,2, is lower than, or equal to, the final residual moisture concentration CS,f, or final desorption rate DRf, then the secondary drying phase is completed.
Step 7
If the calculated residual moisture concentration CS,2 is higher than the final residual moisture concentration CS,f, or the calculated desorption rate DRtheor,2 is higher than the final desorption rate DRf, then using the calculated values of CS,0 and of kinetic constants k0, k1 and k2, it is possible to estimate the final time tf at which the desired residual moisture concentration CS,f, or final desorption rate DRf, is obtained, assuming that the temperature of the product does not change. This can be done by using eq. 11 where CS is replaced by CS,f and, thus, t corresponds to tf:




  
    
      
        
          t
          f
        
        =
        
          
            t
            2
          
          -
          
            
              1
              
                k
                2
              
            
            ⁢
            
              ln
              ⁡
              
                (
                
                  
                    C
                    
                      S
                      ,
                      f
                    
                  
                  
                    C
                    
                      S
                      ,
                      2
                    
                  
                
                )
              
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          18
        
        )
      
    
  




A different stop criterion can be assumed, e.g. the requirement that the desorption rate has a certain low value. For this purpose eq. 12 can be used where DR is replaced by the target value and, thus, t corresponds to tf.
Step 8
A new PRT is performed at time t=tj and a respective desorption rate DRexp,j is calculated; from eq. 12:




  
    
      
        
          
            
              
                DR
                
                  exp
                  ,
                  j
                
              
              =
              
              ⁢
              
                DR
                
                  theor
                  ,
                  j
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  -
                  
                    k
                    j
                  
                
                ⁢
                
                  C
                  
                    S
                    ,
                    0
                  
                
                ⁢
                
                  
                    ∏
                    
                      i
                      =
                      1
                    
                    
                      j
                      -
                      1
                    
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      e
                      
                        -
                        
                          
                            k
                            i
                          
                          ⁡
                          
                            (
                            
                              
                                t
                                i
                              
                              -
                              
                                t
                                
                                  i
                                  -
                                  1
                                
                              
                            
                            )
                          
                        
                      
                    
                    ⁢
                    
                      e
                      
                        -
                        
                          
                            k
                            j
                          
                          ⁡
                          
                            (
                            
                              
                                t
                                j
                              
                              -
                              
                                t
                                
                                  j
                                  -
                                  1
                                
                              
                            
                            )
                          
                        
                      
                    
                  
                
              
            
          
        
      
    
    
      
        (
        
          
            eq
            .
            
              
            
            ⁢
            15
          
          ⁢
          bis
        
        )
      
    
  




This step can be repeated several times, as better explained in the following, and after each PRT a new value of DR is available and a better estimation of the values of CS,0, k0, k1, . . . , kj and tf is obtained, until the end of the secondary drying phase.
For example, at time t=t3 the PRT gives DRexp,3 and from eq. 12 it follows:




  
    
      
        
          
            
              
                DR
                
                  exp
                  ,
                  3
                
              
              =
              
              ⁢
              
                
                  DR
                  
                    theor
                    ,
                    3
                  
                
                =
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  -
                  
                    k
                    3
                  
                
                ⁢
                
                  C
                  
                    s
                    ,
                    0
                  
                
                ⁢
                
                  
                
                ⁢
                
                  e
                  
                    -
                    
                      
                        k
                        1
                      
                      ⁡
                      
                        (
                        
                          
                            t
                            1
                          
                          -
                          
                            t
                            0
                          
                        
                        )
                      
                    
                  
                
                ⁢
                
                  e
                  
                    -
                    
                      
                        k
                        2
                      
                      ⁡
                      
                        (
                        
                          
                            t
                            2
                          
                          -
                          
                            t
                            1
                          
                        
                        )
                      
                    
                  
                
                ⁢
                
                  e
                  
                    -
                    
                      
                        k
                        3
                      
                      ⁡
                      
                        (
                        
                          
                            t
                            3
                          
                          -
                          
                            t
                            2
                          
                        
                        )
                      
                    
                  
                
              
            
          
        
      
    
    
      
        (
        
          
            eq
            .
            
              
            
            ⁢
            15
          
          ⁢
          ter
        
        )
      
    
  




Step 9
Values of constants CS,0, k0, k1, . . . , kj are estimated by solving the non-linear least-square problem:




  
    
      
        
          min
          
            
              C
              
                S
                ,
                0
              
            
            ,
            
              k
              i
            
          
        
        ⁢
        
          
            ∑
            
              i
              =
              0
            
            j
          
          ⁢
          
            
              (
              
                
                  DR
                  
                    exp
                    ,
                    i
                  
                
                -
                
                  DR
                  
                    theor
                    ,
                    i
                  
                
              
              )
            
            2
          
        
      
    
    
      
        (
        
          
            eq
            .
            
              
            
            ⁢
            16
          
          ⁢
          bis
        
        )
      
    
  




assuming, for example, that kj is equal to kj−1, as previously stated.
For example, at time t=t3, the values of constants CS,0, k0, k1, k2 and k3 are calculated by solving the non-linear least-square problem:




  
    
      
        
          min
          
            
              C
              
                S
                ,
                0
              
            
            ,
            
              k
              i
            
          
        
        ⁢
        
          
            ∑
            
              i
              =
              0
            
            3
          
          ⁢
          
            
              (
              
                
                  DR
                  
                    exp
                    ,
                    i
                  
                
                -
                
                  DR
                  
                    theor
                    ,
                    i
                  
                
              
              )
            
            2
          
        
      
    
    
      
        (
        
          
            eq
            .
            
              
            
            ⁢
            16
          
          ⁢
          ter
        
        )
      
    
  




Step 10
Once estimated the values of CS,0, k0, k1, . . . , kj, it is possible to calculate at time t=tj the residual moisture concentration CS,j, using eq. 11, or the desorption rate DRtheor,j, using eq. 12. For example, at time t=t3, once estimated the values of CS,0, k0, k1, k2 and k3, it is possible to calculate the residual moisture concentration CS,3 or the desorption rate DRtheor,3.
Step 11
The calculated value of residual moisture concentration CS,3, or desorption rate DRtheor,3, are compared at time t=tj with the final residual moisture concentration CS,f, or the final desorption rate DRf.
If the calculated residual moisture concentration CS,j, or desorption rate DRtheor,j, is lower than, or equal to, the final residual moisture concentration CS,f, or the final desorption rate DRf, then the secondary drying phase is terminated.
If the calculated residual moisture concentration CS,j, or desorption rate DRtheor,j, is higher than the final residual moisture concentration CS,f, or the final desorption rate DRf, than step 7 is repeated with t=tj for estimating the final time tf at which the final residual moisture concentration CS,f, or final desorption rate DRf, is obtained:




  
    
      
        
          t
          f
        
        =
        
          
            t
            j
          
          -
          
            
              1
              
                k
                j
              
            
            ⁢
            
              ln
              ⁡
              
                (
                
                  
                    C
                    
                      S
                      ,
                      f
                    
                  
                  
                    C
                    
                      S
                      ,
                      j
                    
                  
                
                )
              
            
          
        
      
    
    
      
        (
        
          
            eq
            .
            
              
            
            ⁢
            18
          
          ⁢
          bis
        
        )
      
    
  




For example, at time t=t3 it is possible to estimate the final time tf at which the desired residual moisture concentration CS,f, or final desorption rate DRf, is obtained, assuming that the temperature of the product does not change. This can be done by using eq. 11 where CS is replaced by CS,f and, thus, t corresponds to tf:




  
    
      
        
          t
          f
        
        =
        
          
            t
            3
          
          -
          
            
              1
              
                k
                3
              
            
            ⁢
            
              ln
              ⁡
              
                (
                
                  
                    C
                    
                      S
                      ,
                      f
                    
                  
                  
                    C
                    
                      S
                      ,
                      3
                    
                  
                
                )
              
            
          
        
      
    
    
      
        (
        
          
            eq
            .
            
              
            
            ⁢
            18
          
          ⁢
          ter
        
        )
      
    
  




A different stop criterion can be assumed, i.e. the requirement that the desorption rate has a certain final low value. For this purpose eq. 12 can be used where DR is replaced by the target value and, thus, t corresponds to tf.
Steps 7 to 11 are repeated till the end of secondary drying phase is reached, i.e. till the estimated value of residual moisture concentration Cs,j, or desorption rate DRtheor,j at time tj, is lower than, or equal to, the desired value of residual moisture concentration CS,f, or desorption rate DRf.
In a second version of the method, the desorption rate DR is assumed to depend on the difference between the residual moisture content CS in the solid matrix of the dried product and the equilibrium moisture concentration CS,eq:

DR=−k(CS−CS,eq)  (eq. 19)

The equilibrium moisture concentration Cs,eq is an additional parameter, the value of which can be known (it must be determined experimentally).
Starting from this different expression of desorption rate and repeating the same procedure above described, it is possible to achieve similar results.
The kinetic constant k can be a function of the temperature and can change with time; also the equilibrium moisture concentration Cs,eq changes with temperature, and thus, with time. Again, even if the temperature of the product can change with time, this variation is assumed to be negligible during the time interval between one PRT and the successive, thus allowing the analytical solution of the mass balance equation.
If one PRT is made at t=tj−1 and the successive PRT is made at t=tj, the evolution of the residual moisture concentration, given in % of water and/or solvent per dried mass, is given in the interval [tj−tj−1] by the integration of the following differential equation:




  
    
      
        
          
            d
            ⁢
            
              
            
            ⁢
            
              C
              S
            
          
          
            d
            ⁢
            
              
            
            ⁢
            t
          
        
        =
        
          
            D
            ⁢
            
              
            
            ⁢
            
              R
              j
            
          
          =
          
            -
            
              
                k
                j
              
              ⁡
              
                (
                
                  
                    C
                    S
                  
                  -
                  
                    C
                    
                      S
                      ,
                      eq
                      ,
                      j
                    
                  
                
                )
              
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          20
        
        )
      
    
  




The solution of eq. 20 requires the initial condition, i.e. the value of the residual moisture CS at t=j−1:

CS=CS,j−1e−kj(t−tj−1)+kjCS,eq,j[t−tj−1e−kj(t−tj−1)]  (eq. 21)

The value of CS,j−1 can be calculated from the time integration of eq. 20 in the previous time interval:

CS,j−1=CS,j−2e−kj−1(tj−1−tj−2)++kj−1CS,eq,j−1[tj−1−tj−2e−kj−1(tj−1−tj−2)]  (eq. 22)

and thus:

CS={CS,j−2e−kj−1(tj−1−tj−2)++kj−1CS,eq,j−1[tj−1−tj−2e−kj−1(tj−1−tj−2)]}e−kj(t−tj−1)++kjCS,eq,j[t−tj−1e−kj(t−tj−1)]  (eq. 23)

Similarly, CS,j−2, that is required to get CS,j−1, can be calculated as follow:

CS,j−2=CS,j−3e−kj−2(tj−2−tj−3)++kj−2CS,eq,j−2[tj−2−tj−3e−kj−2(tj−2−tj−3)]  (eq. 24)

This procedure can be iterated until the value of the residual moisture CS,0 at the beginning of the secondary drying stage (t=t0) appears:

CS,1=CS,0e−k1(t1−t0)+k1CS,eq,1[t1−t0e−k1(t1−t0)]  (eq. 25)

Thus, in the time interval between tj and tj−1 the evolution of the residual moisture concentration can be obtained as a function of CS,0, CS,eq,r (with r=1, . . . ,j) and kr (with r=1, . . . ,j).
The evolution of the theoretical value of the desorption rate in the time interval between tj and tj−1 is given by:

DRtheor=−kj{CSj−1e−kj(t−tj−1)++kjCS,eq,j[t−tj−1e−kj(t−tj−1)]−CS,eq,j}  (eq. 26)

and thus it is a function of CS,0, CS,eq,r (with r=1, . . . ,j) and of kr (with r=1, . . . ,j).
If CS,0 and the values of the various kinetic constants kj are perfectly known and the model given by eq. 20 is adequate to describe the dynamics of the system, eq. 21 can be used to know the time evolution of the residual moisture content and thus the time that is required to fulfill the requirements on the final value of the residual moisture content in the product. If the requirement is on the value of the desorption rate, eq. 26 can be used to this purpose.
The above situation is quite rare, since the value of initial moisture concentration has to be calculated by extracting samples and the various kinetic constants are never known a priori.
The method according to the invention provides calculating initial condition Cs, 0 and kinetic constants performing the following steps as shown in the flowchart of FIG. 2.
Step 1
At time t=t0 a PRT is performed and the desorption rate DR (indicated in the following as DRexp,0) is calculated, e.g. using eq. 4.
From eq. 20 it is:

DRexp,0=DRtheor,0=−k0(CS,0−CS,eq,0)  (eq. 27)

Step 2
At time t=t1 a PRT is performed and the desorption rate DR (indicated in the following as DRexp,1) is calculated, e.g. using eq. 4.
From eq. 26 it is:

DRexp,1=DRtheor,1=−k1{CS,0e−k1(t1−t0)++k1CS,eq,1[t1−t0e−k1(t1−t0)]−CS,eq,1}  (eq. 28)

Step 3
At time t=t2 a PRT is performed and the desorption rate DR (indicated in the following ad DRexp,2) is calculated, e.g. using eq. 4.
From eq. 26 it is:

DRexp,2=DRtheor,2=−k2{CS,1e−k2(t2−t1)++k2CS,eq,2[t2−t1e−k2(t2−t1)]−CS,eq,2}  (eq. 29)

Step 4
Values of CS,0, k0, k1 and k2 are estimated so that the calculated values of the desorption rates matches with all the experimental values available (DRexp,0, DRexp,1, DRexp,2). This can be done using a minimization algorithm to solve the following non-linear least-square problem:




  
    
      
        
          min
          
            
              C
              
                S
                ,
                0
              
            
            ,
            
              k
              i
            
          
        
        ⁢
        
          
            ∑
            
              i
              =
              0
            
            2
          
          ⁢
          
            
              (
              
                
                  D
                  ⁢
                  
                    
                  
                  ⁢
                  
                    R
                    
                      exp
                      ,
                      i
                    
                  
                
                -
                
                  D
                  ⁢
                  
                    
                  
                  ⁢
                  
                    R
                    
                      theor
                      ,
                      i
                    
                  
                
              
              )
            
            2
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          30
        
        )
      
    
  




assuming, for example, that k2 is equal to k1, as previously stated.
The values of CS,eq,0, CS,eq,1 and CS,eq,2 must be known (from experimentation).
Step 5
Once estimated the values of CS,0, k0, k1 and k2, it is possible to calculate at time t=t2 the residual moisture concentration CS,2 (or the desorption rate), using eq. 20.
Step 6
The calculated value of residual moisture concentration CS,2 is compared with a desired value of a final residual moisture concentration CS,f.
If the calculated value of residual moisture concentration CS,2 is lower than, or equal to, the final residual moisture concentration CS,f, then the secondary drying phase is completed.
Step 7
If the calculated value of residual moisture concentration CS,2 is higher than the desired final residual moisture concentration CS,f, then using the calculated values of CS,0 and of the kinetic constants it is possible to estimate the time tf at which the desired value of residual moisture concentration CS,f is obtained, assuming that the temperature of the product does not change. This can be done by using eq. 21 where CS is replaced by CS,f and thus t corresponds to tf. In this case the following non-linear equation must be solved:

CS,f=CS,2e−k2(tf−t2)+k2CS,eq,2[tf−t2e−k2(tf−t2)]  (eq. 31)

A different stop criterion can be assumed, e.g. the requirement that the desorption rate DR has a certain final low value DRf. For this purpose eq. 26 can be used wherein DR is replaced by final desorption rate DRf.
Step 8
A new PRT is performed at time t=tj and a respective desorption rate DRexp,j is calculated; from eq. 26:

DRexp,j=DRtheor,j=−kj{CSj−1e−kj(tj−tj−1)++kjCS,eq,j[tj−tj−1e−kj(tj−tj−1)]−CS,eq,j}  (eq. 29bis)

This step can be repeated several times and after each PRT a new value of DR is available and a better estimation of the values of CS,0, k0, k1, . . . , kj and tf is obtained, until the end of the secondary drying phase.
For example, at time t=t3 the PRT gives DRexp,3 and from eq. 26 it is:

DRexp,3=DRtheor,3=−k3{CS,2e−k3(t3−t2)++k3CS,eq,3[t3−t2e−k3(t3−t2)]−CS,eq,3}(eq. 29ter)

Step 9
Values of CS,0, k0, k1, . . . , kj are estimated by solving the non-linear least-square problem:




  
    
      
        
          min
          
            
              C
              
                S
                ,
                0
              
            
            ,
            
              k
              i
            
          
        
        ⁢
        
          
            ∑
            
              i
              =
              0
            
            j
          
          ⁢
          
            
              (
              
                
                  D
                  ⁢
                  
                    
                  
                  ⁢
                  
                    R
                    
                      exp
                      ,
                      i
                    
                  
                
                -
                
                  D
                  ⁢
                  
                    
                  
                  ⁢
                  
                    R
                    
                      theor
                      ,
                      i
                    
                  
                
              
              )
            
            2
          
        
      
    
    
      
        (
        
          
            eq
            .
            
              
            
            ⁢
            30
          
          ⁢
          bis
        
        )
      
    
  




assuming, for example, that kj is equal to kj−1, as previously stated.
For example, at time t=t3, the values CS,0, k0, k1, k2 and k3 are calculated by solving the non-linear least-square problem:




  
    
      
        
          min
          
            
              C
              
                S
                ,
                0
              
            
            ,
            
              k
              i
            
          
        
        ⁢
        
          
            ∑
            
              i
              =
              0
            
            3
          
          ⁢
          
            
              (
              
                
                  D
                  ⁢
                  
                    
                  
                  ⁢
                  
                    R
                    
                      exp
                      ,
                      i
                    
                  
                
                -
                
                  D
                  ⁢
                  
                    
                  
                  ⁢
                  
                    R
                    
                      theor
                      ,
                      i
                    
                  
                
              
              )
            
            2
          
        
      
    
    
      
        (
        
          
            eq
            .
            
              
            
            ⁢
            30
          
          ⁢
          ter
        
        )
      
    
  




Step 10
Once estimated the values of CS,0, k0, k1, . . . , kj, it is possible to calculate at time t=tj the residual moisture concentration Cs,j using eq. 20, or the desorption rate DRtheor,j.
Step 11
The calculated value of residual moisture concentration CS,j, or desorption rate DRtheor,j, is compared with the final residual moisture concentration CS,f, or the final desorption rate DRf.
If the estimated value of residual moisture concentration CS,j, or desorption rate DRtheor,j, is lower than, or equal to, the final residual moisture concentration CS,f, or the final desorption rate DRf, then secondary drying phase is completed.
If the estimated value of residual moisture concentration CS,j, or desorption rate DRtheor,j, is higher than final residual moisture concentration CS,f, or final desorption rate DRf, than step 7 is repeated with t=tj for estimating the final time tf at which the final residual moisture concentration CS,f (or final desorption rate DRf) is obtained:

CS,f=CS,je−kj(tf−tj)+kjCS,eq,j[tf−tje−kj(tf−tj)]  (eq. 31bis)

For example, at time t=t3, using the calculated values of CS,0 and of the kinetic constants it is possible to estimate the time instant tf at which the final residual moisture concentration CS,f is obtained, assuming that the temperature of the product does not change; the following non-linear equation must be solved:

CS,f=CS,3e−k3(tf−t3)+k3CS,eq,3[tf−t3e−k3(tf−t3)]  (eq. 31ter)

A different stop criterion can be assumed, e.g. the requirement that the desorption rate has a certain low value.
In the following and with reference to FIGS. 3 to 7, it is provided an example of the application of the method of the invention for monitoring a secondary drying phase of a drying process.
FIG. 3 shows an experimental campaign with provides values of desorption rate vs. time during the secondary drying.
The first version of the method is used.
Step 1
At time t=t0=0 s from the PRT (and eq. 4) it comes that DRexp,0=0.00056% water over dried product s−1.
Step 2
At time t=t1=1296 s from the PRT (and eq. 4) it comes that DRexp,1=0.00049% water s−1.
Step 3
At time t=t2=2592 s from PRT (and eq. 4) it comes that DRexp,2=0.00035% water s−1.
Step 4
Using the preliminary estimation of the kinetic constants k0 and k1 and of CS,0 from eq. 17 (k0=k1=1.03·10−4 s−1, CS,0=5.48% water over dried product), eq. 16 is used to calculate CS,0 and the kinetic constants (CS,0=4.13% water over dried product).
Steps 5, 7
Using the calculated values of CS,0 and of the kinetic constants and eq. 18 it is possible to estimate the time instant tf at which the desired value of final moisture concentration CS,f (e.g. 0.2% water over dried product) is obtained. In this case it is calculated that 25056 s are still required.
FIG. 4 shows an estimation of the time evolution of the concentration CS and of the desorption rate DR obtained using the estimation of CS, 0 and of the kinetic constants.
At this point, the above described procedure can be iterated (steps 7 to 11).
At time t=t3=3888 s from PRT (and eq. 4) it comes that DRexp,3=0.00028% water s−1.
Using eq. 16 it calculated CS,0=4.06% water over dried product and that 26352 s are still required.
FIG. 5 shows the estimation of the time evolution of the concentration CS and of the desorption rate DR obtained using the new estimation of CS, 0 and of the kinetic constants.
It is possible to see that at each iteration the estimation of the values of CS,0 is improved, as well as the estimation of the time tf required to complete the secondary drying phase.
FIG. 6 shows how the estimate of the final time tf required to complete the secondary drying phase changes with time.
FIG. 7 illustrates a comparison between estimations of final time tf required to complete the secondary drying phase (end-points of secondary drying phase) using the method of the invention (broken line with round dots) and using the method according to the U.S. Pat. No. 6,176,121 (broken line with square dots).
It is possible to see that the estimations of the time required to get the end of the secondary drying using the method of the invention is quite good even at the beginning of the phase and is refined as the secondary drying goes on. On the contrary, using the method disclosed in U.S. Pat. No. 6,176,121 the prediction of the time required to complete secondary drying is not reliable at the beginning and after each PRT the prediction is updated until the end of the drying is obtained.
The method of the invention was also validated by means of a series of experiments carried out in laboratory.
FIGS. 8 and 9 are an example of the results that can be obtained with the algorithm of the method are used.
In particular, FIGS. 8 and 9 are a comparison between the experimental values (symbols) and those predicted by the algorithm of the invention (solid line) respectively of the desorption rate (FIG. 8) and of the residual water content (FIG. 9). The time evolution of a shelf temperature is also shown (FIG. 8, dotted line). Time is equal to zero at the beginning of the secondary drying.
The example refers to a freeze-drying cycle of an aqueous solution of sucrose at 20% by weight (155 vials having a diameter of 20.85·10−3 m, filled with 3·10−3 1 of solution). The freezing phase was carried out at −50° C. for 17 h, primary drying phase was carried out at −15° C. and 10 Pa for 25 h and secondary drying phase was carried out at 20° C.
The experimental values of desorption rate have been obtained by means of the Pressure Rise Test (see eq. (4)), while the residual water content was determined by weighing some vials taken from the drying chamber using a sample thief.
The kinetic model for the desorption of water that was used by the algorithm is the same of the first version of the method (eq. 5-18), i.e. the desorption rate was assumed to be proportional to the residual water content.
The time evolution of the desorption rate is a consequence of the fact that when secondary drying is started the shelf temperature is increased and, during this time interval, the product temperature, and thus the desorption rate, increases. After this, the temperature remains constant and, due to the lowering of the residual water content, the desorption rate decreases.

## Claims

1. Method for calculating a residual moisture concentration and/or desorption rate of a product during a secondary drying phase of a freeze-drying process in a freeze-dryer apparatus including a drying chamber and a condenser chamber where vapour generated by drying process flows, said apparatus being provided with a pressure sensor that can be isolated for performing pressure rise tests and measuring a total pressure inside said drying chamber, said method comprising the steps of:
a) setting a desired final residue moisture concentration and/or a desired final desorption rate of said product;
b) measuring initial residual moisture concentration and/or desorption rate and estimating kinetic constants of a kinetic model of the drying process, said kinetic model being suitable for calculating the residual moisture concentration and/or desorption rate of said product;
c) closing a valve placed on a duct connecting said drying chamber to said condenser chamber for a preselected period of time;
d) measuring a pressure change in said drying chamber;
e) calculating a desorption rate from the closing time period of step c) and the pressure change of step d);
f) repeating steps c)-e) to calculate residual moisture concentrations and/or desorption rates at pre-specified time intervals;
g) integrating the calculated residual moisture concentrations and/or desorption rate from said initial conditions to the calculated residual moisture concentrations and/or desorption rates of the pre-specified time intervals of step e) to determine a current residual moisture concentration of the drying product;
wherein if said current residual moisture concentration and/or said current final desorption rate is lower than or equal to said desired residual moisture concentration and/or to said desired desorption rate then said secondary drying phase ends, and wherein
a first pressure rise test at time t=t0 at the beginning of the secondary drying phase is performed and a first value of experimental desorption rate (DRexp, 0) of said product is calculated using the equation:






  
    
      
        
          D
          ⁢
          
            
          
          ⁢
          
            R
            exp
          
        
        =
        
          
            
              VM
              w
            
            RT
          
          ⁢
          
            
              (
              
                
                  d
                  ⁢
                  
                    
                  
                  ⁢
                  P
                
                
                  d
                  ⁢
                  
                    
                  
                  ⁢
                  t
                
              
              )
            
            
              t
              =
              
                t
                0
              
            
          
          ⁢
          
            100
            
              m
              dried
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          4
        
        )
      
    
  




where:
DRexp: experimental desorption rate, [% water and/or solvent s−1]
P: measured pressure, [Pa]
t: time, [s]
t0: time instant at the beginning of the pressure rise test, [s]
R: gas constant [8,314 J mol−1 K−1]
T: temperature of the vapour, [K]
V: (free) volume of drying chamber, [m3]
Mw: molecular weight of water and/or solvent, [kg mol−1]
mdried: mass of the dried product, [kg]
a second pressure rise test at a successive time t=t1 is performed and a second value of experimental desorption rate (DRexp, 1) of said product is calculated using the equation:






  
    
      
        
          
            DR
            exp
          
          =
          
            
              
                VM
                W
              
              RT
            
            ⁢
            
              
                (
                
                  dP
                  dt
                
                )
              
              
                t
                =
                
                  t
                  0
                
              
            
            ⁢
            
              100
              
                m
                dried
              
            
          
        
        ;
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          4
        
        )
      
    
  




a third pressure rise test at time t=t2 is performed and a third value of experimental desorption rate (DRexp, 2) of said product is calculated using the equation eq. 4






  
    
      
        
          
            DR
            exp
          
          =
          
            
              
                VM
                W
              
              RT
            
            ⁢
            
              
                (
                
                  dP
                  dt
                
                )
              
              
                t
                =
                
                  t
                  0
                
              
            
            ⁢
            
              100
              
                m
                dried
              
            
          
        
        ;
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          4
        
        )
      
    
  




and wherein
said initial conditions (CS,0) and kinetic constants (k0, k1, k2) of a kinetic model of the drying process, said kinetic model being suitable for calculating a residual moisture content (CS) and/or desorption rate (DRtheor) of said product are calculated using the equations:


DRexp,0=DRtheor,0=−k0CS,0  (eq. 13)



DRexp,1=DRtheor,1=−k1CS,0e−k1(t1−t0)  (eq. 14)



DRexp,2=DRtheor,2=−k2Cs,0e−k1(t1−t0)e−k2(t2−t1)  (eq. 15)


where:
DRexp,j: experimental desorption rate at time tj, [% water and/or solvent s−1]
DRtheor,j: desorption rate (theoretical value) at time tj [% water and/or solvent s−1]
CS,0: value of the residual moisture [% water and/or solvent over dried product] at the beginning of the secondary drying phase (t=t0);
kj: kinetic constant of the process at time t=tj (with j=0, 1, 2), [s−1];
and a minimisation algorithm to solve the minimum least square problem described by equation:






  
    
      
        
          min
          
            
              C
              
                S
                ,
                0
              
            
            ,
            
              k
              i
            
          
        
        ⁢
        
          
            ∑
            
              i
              =
              0
            
            2
          
          ⁢
          
            
              (
              
                
                  DR
                  
                    exp
                    ,
                    i
                  
                
                -
                
                  DR
                  
                    theor
                    ,
                    i
                  
                
              
              )
            
            2
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          16
        
        )
      
    
  




calculating at time t=t2 a respective residual moisture content (CS,2) and a respective desorption rate (DRtheor, 2) (step 5) respectively using the equations:






  
    
      
        
          C
          S
        
        =
        
          
            C
            
              S
              ,
              0
            
          
          ⁢
          
            
              ∏
              
                i
                =
                1
              
              
                j
                -
                1
              
            
            ⁢
            
              
            
            ⁢
            
              
                e
                
                  -
                  
                    
                      k
                      i
                    
                    ⁡
                    
                      (
                      
                        
                          t
                          i
                        
                        -
                        
                          t
                          
                            i
                            -
                            1
                          
                        
                      
                      )
                    
                  
                
              
              ⁢
              
                e
                
                  -
                  
                    
                      k
                      j
                    
                    ⁡
                    
                      (
                      
                        t
                        -
                        
                          t
                          
                            j
                            -
                            1
                          
                        
                      
                      )
                    
                  
                
              
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          11
        
        )
      
    
  
  
    
      
        
          DR
          theor
        
        =
        
          
            -
            
              k
              j
            
          
          ⁢
          
            C
            
              S
              ,
              0
            
          
          ⁢
          
            
              ∏
              
                i
                =
                1
              
              
                j
                -
                1
              
            
            ⁢
            
              
            
            ⁢
            
              
                e
                
                  -
                  
                    
                      k
                      i
                    
                    ⁡
                    
                      (
                      
                        
                          t
                          i
                        
                        -
                        
                          t
                          
                            i
                            -
                            1
                          
                        
                      
                      )
                    
                  
                
              
              ⁢
              
                
                  e
                  
                    -
                    
                      
                        k
                        j
                      
                      ⁡
                      
                        (
                        
                          t
                          -
                          
                            t
                            
                              j
                              -
                              1
                            
                          
                        
                        )
                      
                    
                  
                
                .
              
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          12
        
        )

2. Method according to claim 1, further comprising after step 5 the steps of:
comparing said residual moisture content (CS,2) and/or said desorption rate (DRtheor,2) calculated at time t=t2 respectively with a desired final residual moisture concentration (CS,f) and/or a desired final desorption rate (DRf) (step 6); if said residual moisture content (CS,2) is lower than, or equal to, said final residual moisture concentration (CS,f) or said desorption rate (DRtheor,2) is lower than, or equal to, said final desorption rate (DRf), then the secondary drying phase is considered ended; if not the method further comprising the steps of:
estimating a final time (tf) at which said final residual moisture concentration (CS,f) or said final desorption rate (DRf) is obtained (step 7);
performing a further pressure rise test at time t=tj and calculating at said time t=tj a respective residual moisture content (CS,j) and a respective desorption rate (DRtheor,j) (step 8);
estimating initial conditions (CS,0) and kinetic constants (k0, k1, k2, . . . , kj) of said kinetic model (step 9);
calculating at said time t=tj a respective residual moisture content (CS,j) and/or a respective desorption rate (DRtheor,j) (step 10);
comparing said residual moisture content (CS,j) and/or said desorption rate (DRtheor,j) calculated at said time t=tj respectively with said final residual moisture concentration (CS,f) and/or said final desorption rate (DRf) (step 11); if said residual moisture content (CS,j) is lower than, or equal to, said final residual moisture concentration (CS,f) or said desorption rate (DRtheor,j) is lower than, or equal to, said final desorption rate (DRf) then the secondary drying phase is considered ended; if not steps 7 to 11 are repeated.

3. Method according to claim 2, wherein said experimental desorption rates (DRexp,O, DRexp,1, DRexp,2) are calculated using the equation:




  
    
      
        
          DR
          exp
        
        =
        
          
            
              VM
              W
            
            RT
          
          ⁢
          
            
              (
              
                dP
                dt
              
              )
            
            
              t
              =
              
                t
                0
              
            
          
          ⁢
          
            100
            
              m
              dried
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          4
        
        )
      
    
  



where:
DRexp: experimental desorption rate, [% water and/or solvent s−1]
P: measured pressure, [Pa]
t: time, [s]
t0: time instant at the beginning of the pressure rise test, [s]
R: gas constant [8,314 J mol−1 K−1]
T: temperature of the vapour, [K]
V: (free) volume of drying chamber, [m3]
Mw: molecular weight of water and/or solvent, [kg mol−1]
mdried: mass of the dried product, [kg].

4. Method according to claim 3, wherein said kinetic model comprises mathematical equations suitable for modeling the dependence of the desorption rate (DR) on the residual moisture content (CS) in the product.

5. Method according to claim 3, wherein said desorption rate is assumed to depend on said residual moisture content in said product according to the equation:


DR=−kCS  (eq. 5)

where:
DR: desorption rate, [% water and/or solvent s−1]
k: kinetic constant of the process, [s−1]
CS: residual moisture content, [% water/solvent over dried product].

6. Method according to claim 5, wherein a time evolution of said residual moisture concentration (CS) at time t=tj is given by the integration of the following differential equation:




  
    
      
        
          
            d
            ⁢
            
              
            
            ⁢
            
              C
              S
            
          
          
            d
            ⁢
            
              
            
            ⁢
            t
          
        
        =
        
          
            D
            ⁢
            
              
            
            ⁢
            
              R
              j
            
          
          =
          
            
              -
              
                k
                j
              
            
            ⁢
            
              C
              S
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          7
        
        )
      
    
  



where:
DRj: desorption rate at time t=tj, [% water and/or solvent s−1]
t: time, [s]
kj: kinetic constant of the process at time t=tj, [s−1].

7. Method according to claim 6, wherein said calculating a residual moisture content (Cs) is made by means of the equation:




  
    
      
        
          C
          S
        
        =
        
          
            C
            
              S
              ,
              0
            
          
          ⁢
          
            
              ∏
              
                i
                =
                1
              
              
                j
                -
                1
              
            
            ⁢
            
              
            
            ⁢
            
              
                e
                
                  -
                  
                    
                      k
                      i
                    
                    ⁡
                    
                      (
                      
                        
                          t
                          i
                        
                        -
                        
                          t
                          
                            i
                            -
                            1
                          
                        
                      
                      )
                    
                  
                
              
              ⁢
              
                e
                
                  -
                  
                    
                      k
                      j
                    
                    ⁡
                    
                      (
                      
                        t
                        -
                        
                          t
                          
                            j
                            -
                            1
                          
                        
                      
                      )
                    
                  
                
              
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          11
        
        )
      
    
  



where:
CS,0: value of the residual moisture [% water and/or solvent over dried product] at the beginning of the secondary drying phase (t=t0);
kr: kinetic constant of the process at time t=tr (with r=1, 2, . . . , j), [s−1].

8. Method according to claim 7, wherein said calculating a desorption rate (DRtheor) is made by means of the equation:




  
    
      
        
          D
          ⁢
          
            
          
          ⁢
          
            R
            theor
          
        
        =
        
          
            -
            
              k
              j
            
          
          ⁢
          
            C
            
              S
              ,
              0
            
          
          ⁢
          
            
              ∏
              
                i
                =
                1
              
              
                j
                -
                1
              
            
            ⁢
            
              
            
            ⁢
            
              
                e
                
                  -
                  
                    
                      k
                      i
                    
                    ⁡
                    
                      (
                      
                        
                          t
                          i
                        
                        -
                        
                          t
                          
                            i
                            -
                            1
                          
                        
                      
                      )
                    
                  
                
              
              ⁢
              
                
                  e
                  
                    -
                    
                      
                        k
                        j
                      
                      ⁡
                      
                        (
                        
                          t
                          -
                          
                            t
                            
                              j
                              -
                              1
                            
                          
                        
                        )
                      
                    
                  
                
                .
              
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          12
        
        )

9. Method according to claim 8, wherein said estimating initial conditions (CS,0) and kinetic constants (k0, k1, k2, . . . , kj), at time t=tj, is made by means of the following equations:




  
    
      
        
          D
          ⁢
          
            
          
          ⁢
          
            R
            
              exp
              ,
              0
            
          
        
        =
        
          
            D
            ⁢
            
              
            
            ⁢
            
              R
              
                theor
                ,
                0
              
            
          
          =
          
            
              -
              
                k
                0
              
            
            ⁢
            
              C
              
                S
                ,
                0
              
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          13
        
        )
      
    
  
  
    
      
        
          D
          ⁢
          
            
          
          ⁢
          
            R
            
              exp
              ,
              1
            
          
        
        =
        
          
            D
            ⁢
            
              
            
            ⁢
            
              R
              
                theor
                ,
                1
              
            
          
          =
          
            
              -
              
                k
                1
              
            
            ⁢
            
              C
              
                S
                ,
                0
              
            
            ⁢
            
              e
              
                -
                
                  
                    k
                    1
                  
                  ⁡
                  
                    (
                    
                      
                        t
                        1
                      
                      -
                      
                        t
                        0
                      
                    
                    )
                  
                
              
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          14
        
        )
      
    
  
  
    
      
        
          D
          ⁢
          
            
          
          ⁢
          
            R
            
              exp
              ,
              2
            
          
        
        =
        
          
            D
            ⁢
            
              
            
            ⁢
            
              R
              
                theor
                ,
                2
              
            
          
          =
          
            
              -
              
                k
                2
              
            
            ⁢
            
              C
              
                s
                ,
                0
              
            
            ⁢
            
              e
              
                -
                
                  
                    k
                    1
                  
                  ⁡
                  
                    (
                    
                      
                        t
                        1
                      
                      -
                      
                        t
                        0
                      
                    
                    )
                  
                
              
            
            ⁢
            
              e
              
                -
                
                  
                    k
                    2
                  
                  ⁡
                  
                    (
                    
                      
                        t
                        2
                      
                      -
                      
                        t
                        1
                      
                    
                    )
                  
                
              
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          15
        
        )
      
    
  
  
    
      
        
          
            
              
                D
                ⁢
                
                  
                
                ⁢
                
                  R
                  
                    exp
                    ,
                    j
                  
                
              
              =
              
              ⁢
              
                
                  D
                  ⁢
                  
                    
                  
                  ⁢
                  
                    R
                    
                      theor
                      ,
                      j
                    
                  
                
                =
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  -
                  
                    k
                    j
                  
                
                ⁢
                
                  C
                  
                    S
                    ,
                    0
                  
                
                ⁢
                
                  
                    ∏
                    
                      i
                      =
                      1
                    
                    
                      j
                      -
                      1
                    
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      e
                      
                        -
                        
                          
                            k
                            i
                          
                          ⁡
                          
                            (
                            
                              
                                t
                                i
                              
                              -
                              
                                t
                                
                                  i
                                  -
                                  1
                                
                              
                            
                            )
                          
                        
                      
                    
                    ⁢
                    
                      e
                      
                        -
                        
                          
                            k
                            j
                          
                          ⁡
                          
                            (
                            
                              
                                t
                                j
                              
                              -
                              
                                t
                                
                                  j
                                  -
                                  1
                                
                              
                            
                            )
                          
                        
                      
                    
                  
                
              
            
          
        
      
    
    
      
        (
        
          
            eq
            .
            
              
            
            ⁢
            15
          
          ⁢
          
            
          
          ⁢
          bis
        
        )
      
    
  



and solving the following non-linear least square problem:





  
    
      
        
          min
          
            
              C
              
                S
                ,
                0
              
            
            ,
            
              k
              i
            
          
        
        ⁢
        
          
            ∑
            
              i
              =
              0
            
            j
          
          ⁢
          
            
              
                (
                
                  
                    D
                    ⁢
                    
                      
                    
                    ⁢
                    
                      R
                      
                        exp
                        ,
                        i
                      
                    
                  
                  -
                  
                    D
                    ⁢
                    
                      
                    
                    ⁢
                    
                      R
                      
                        theor
                        ,
                        i
                      
                    
                  
                
                )
              
              2
            
            .
          
        
      
    
    
      
        (
        
          
            eq
            .
            
              
            
            ⁢
            16
          
          ⁢
          bis
        
        )

10. Method according to claim 7, wherein said final time (tf) is calculated, assuming that temperature of said product does not change, by means of the following equation, resulted from (eq. 11):




  
    
      
        
          t
          f
        
        =
        
          
            t
            j
          
          -
          
            
              1
              
                k
                j
              
            
            ⁢
            
              ln
              ⁡
              
                (
                
                  
                    C
                    
                      S
                      ,
                      f
                    
                  
                  
                    C
                    
                      S
                      ,
                      j
                    
                  
                
                )
              
            
          
        
      
    
    
      
        (
        
          
            eq
            .
            
              
            
            ⁢
            18
          
          ⁢
          bis
        
        )
      
    
  



where:
CS,f: final residual moisture concentration [% water and/or solvent over dried product];
CS,j: residual moisture concentration at time t=tj [% water and/or solvent over dried product].

11. Method for calculating a residual moisture concentration and/or desorption rate of a product during a secondary drying phase of a freeze-drying process in a freeze-dryer apparatus including a drying chamber and a condenser chamber where vapour generated by drying process flows, said apparatus being provided with a pressure sensor that can be isolated for performing pressure rise tests and measuring a total pressure inside said drying chamber, said method comprising the steps of:
a) setting a desired final residue moisture concentration and/or a desired final desorption rate of said product;
b) measuring initial residual moisture concentration and/or desorption rate and estimating kinetic constants of a kinetic model of the drying process, said kinetic model being suitable for calculating the residual moisture concentration and/or desorption rate of said product;
c) closing a valve placed on a duct connecting said drying chamber to said condenser chamber for a preselected period of time;
d) measuring a pressure change in said drying chamber;
e) calculating a desorption rate from the closing time period of step c) and the pressure change of step d);
f) repeating steps c)-e) to calculate residual moisture concentration and/or desorption rate at pre-specified time intervals;
g) integrating the calculated residual moisture concentration and/or desorption rate from said initial conditions to the calculated residual moisture concentration and/or desorption rates of the pre-specified time intervals of step e) to determine a current residual moisture of the drying product;
wherein if said current residual moisture concentration and/or said current final desorption rate is lower than or equal to said desired residual moisture concentration and/or said desired desorption rate then said secondary drying phase ends, and wherein
a first pressure rise test at time t=t0 at the beginning of the secondary drying phase is performed and a first value of experimental desorption rate (DRexp, 0) of said product is calculated using the equation:






  
    
      
        
          DR
          exp
        
        =
        
          
            
              VM
              W
            
            RT
          
          ⁢
          
            
              (
              
                dP
                dt
              
              )
            
            
              t
              =
              
                t
                0
              
            
          
          ⁢
          
            100
            
              m
              dried
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          4
        
        )
      
    
  




where:
DRexp: experimental desorption rate, [% water and/or solvent s−1]
P: measured pressure, [Pa]
t: time, [s]
t0: time instant at the beginning of the pressure rise test, [s]
R: gas constant [8,314 J mol−1 K−1]
T: temperature of the vapour, [K]
V: (free) volume of drying chamber, [m3]
Mw: molecular weight of water and/or solvent, [kg mol−1]
mdried: mass of the dried product, [kg]
a second pressure rise test at a successive time t=t1 is performed and a second value of experimental desorption rate (DRexp, 1) of said product is calculated using the equation:






  
    
      
        
          
            DR
            exp
          
          =
          
            
              
                VM
                W
              
              RT
            
            ⁢
            
              
                (
                
                  dP
                  dt
                
                )
              
              
                t
                =
                
                  t
                  0
                
              
            
            ⁢
            
              100
              
                m
                dried
              
            
          
        
        ;
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          4
        
        )
      
    
  




a third pressure rise test at time t=t2 is performed and a third value of experimental desorption rate (DRexp, 2) of said product is calculated using the equation eq. 4






  
    
      
        
          
            DR
            exp
          
          =
          
            
              
                VM
                W
              
              RT
            
            ⁢
            
              
                (
                
                  dP
                  dt
                
                )
              
              
                t
                =
                
                  t
                  0
                
              
            
            ⁢
            
              100
              
                m
                dried
              
            
          
        
        ;
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          4
        
        )
      
    
  




and wherein
said initial conditions (CS,0) and kinetic constants (k0, k1, k2) of a kinetic model of the drying process, said kinetic model being suitable for calculating a residual moisture content (CS) and/or desorption rate (DRtheor) of said product are calculated using the equations:


DRexp,0=DRtheor,0=−k0(CS,0−CS,eq,0)  (eq. 27)



DRexp,1=DRtheor,1=−k1{CS,0e−k1(t1−t0)++k1CS,eq,1[t1−t0e−k1(t1−t0)]−CS,eq,1}  (eq. 28)



DRexp,2=DRtheor,2=−k2{CS,1e−k2(t2−t1)++k2CS,eq,2[t2−t1e−k2(t2−t1)]−CS,eq,2}  (eq. 29)


where
DRexp,j: experimental desorption rate at time tj, [% water and/or solvent s−1]
Cs,eq,j equilibrium moisture concentration at time tj [% water and/or solvent s−1]:
CS,0: value of the residual moisture [% water and/or solvent over dried product] at the beginning of the secondary drying phase (t=t0);
kj: kinetic constant of the process at time t=tj (with j=0, 1, 2), [s−1]. and a minimisation algorithm to solve the minimum least square problem described by equation:






  
    
      
        
          min
          
            
              C
              
                S
                ,
                0
              
            
            ,
            
              k
              i
            
          
        
        ⁢
        
          
            ∑
            
              i
              =
              0
            
            2
          
          ⁢
          
            
              (
              
                
                  DR
                  
                    exp
                    ,
                    i
                  
                
                -
                
                  DR
                  
                    theor
                    ,
                    i
                  
                
              
              )
            
            2
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          30
        
        )
      
    
  




calculating at time t=t2 a respective residual moisture content (CS,2) and a respective desorption rate (DRtheor, 2) (step 5) using for calculating the residual moisture content (CS,2) the equations:


CS=CS,j−1e−kj(t−tj−1)++kjCS,eq,j[t−tj−1e−kj(t−tj−1)]  (eq. 21)



CS,j−1=CS,j−2e−kj−1(tj−1−tj−2)++kj−1CS,eq,j−1[tj−1−tj−2e−kj−1(tj−1−tj−2)]  (eq. 22)



CS,j−2=CS,j−3e−kj−2(tj−2−tj−3)++kj−2CS,eq,j−2[tj−2−tj−3e−kj−2(tj−2−tj−3)]  (eq. 24)



CS,1=CS,0e−k1(t1−t0)+k1CS,eq,1[t1−t0e−k1(t1−t0)]  (eq. 25)

where:
CS,0: value of the residual moisture [% water and/or solvent over dried product] at the beginning of the secondary drying phase (t=t0);
kr: kinetic constant of the process at time t=tr (with r=1, 2, . . . , j), [s−1];
Cs,eq,r: equilibrium moisture concentration at time t=tr with r=1, 2, . . . , j), [% water and/or solvent over dried product];

and using for calculating the respective desorption rate (DRtheor, 2) the equation:


DRtheor=−kj{CSj−1e−kj(t−tj−1)++kjCS,eq,j[t−tj−1e−kj(t−tj−1)]−CS,eq,j}  (eq. 26).

12. Method according to claim 11, wherein said desorption rate (DRtheor) is assumed to depend on said residual moisture content (CS) in said product according to the equation:


DR=−k(CS−CS,eq)  (eq. 19)

where:
DR: desorption rate, [% water and/or solvent s−1]
k: kinetic constant of the process, [s−1]
CS: residual moisture concentration, [% water and/or solvent over dried product]
Cs,eq: equilibrium moisture concentration, [% water and/or solvent over dried product].

13. Method according to claim 12, wherein a time evolution of said residual moisture concentration (CS) at time t=tj is given by the integration of the following differential equation:




  
    
      
        
          
            d
            ⁢
            
              
            
            ⁢
            
              C
              S
            
          
          
            d
            ⁢
            
              
            
            ⁢
            t
          
        
        =
        
          
            D
            ⁢
            
              
            
            ⁢
            
              R
              j
            
          
          =
          
            -
            
              
                k
                j
              
              ⁡
              
                (
                
                  
                    C
                    S
                  
                  -
                  
                    C
                    
                      S
                      ,
                      eq
                      ,
                      j
                    
                  
                
                )
              
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          20
        
        )
      
    
  



where:
DRj: desorption rate at time t=tj, [% water and/or solvent s−1]
t: time, [s]
k1: kinetic constant of the process, [s−1],
Cs,eq,j: equilibrium moisture concentration at time t=tj, [% water and/or solvent over dried product].

14. Method according to claim 13, wherein said calculating a residual moisture content (CS) at time t=tj is made by means of the following equations:


CS=CS,j−1e−kj(t−tj−1)+—kjCS,eq,j[t−tj−1e−kj(t−tj−1)]  (eq. 21)



and



CS,j−1=CS,j−2e−kj−1(tj−1−tj−2)++kj−1CS,eq,j−1[tj−1−tj−2e−kj−1(tj−1−tj−2)]  (eq. 22)



CS,j−2=CS,j−3e−kj−2(tj−2−tj−3)++kj−2CS,eq,j−2[tj−2−tj−3e−kj−2(tj−2−tj−3)]  (eq. 24)



CS,1=CS,0e−k1(t1−t0)+k1CS,eq,1[t1−t0e−k1(t1−t0)]  (eq. 25)

where:
CS,0: value of the residual moisture [% water and/or solvent over dried product] at the beginning of the secondary drying phase (t=t0);
kr: kinetic constant of the process at time t=tr (with r=1, 2, . . . , j), [s−1];
Cs,eq,r: equilibrium moisture concentration at time t=tr with r=1, 2, . . . , j), [% water and/or solvent over dried product].

15. Method according to claim 14, wherein said calculating a desorption rate (DRtheor) is made by means of the equation:


DRtheor=−kj{CS,j−1e−kj(t−tj−1)++kjCS,eq,j[t−tj−1e−kj(t−tj−1)]−CS,eq,j}  (eq. 26)

16. Method according to claim 15, wherein said experimental desorption rates (DRexp,0, DRexp,1, DRexp,2) are calculated using the equation:




  
    
      
        
          D
          ⁢
          
            
          
          ⁢
          
            R
            exp
          
        
        =
        
          
            
              VM
              w
            
            RT
          
          ⁢
          
            
              (
              
                
                  d
                  ⁢
                  
                    
                  
                  ⁢
                  P
                
                
                  d
                  ⁢
                  
                    
                  
                  ⁢
                  t
                
              
              )
            
            
              t
              =
              
                t
                0
              
            
          
          ⁢
          
            100
            
              m
              dried
            
          
        
      
    
    
      
        (
        
          eq
          .
          
            
          
          ⁢
          4
        
        )
      
    
  



where:
DRexp: experimental desorption rate, [% water and/or solvent s−1]
P: measured pressure, [Pa]
t: time, [s]
t0: time instant at the beginning of the pressure rise test, [s]
R: gas constant [8,314 J mol−1 K−1]
T: temperature of the vapour, [K]
V: (free) volume of drying chamber, [m3]
Mw: molecular weight of water and/or solvent, [kg mol−1]
mdried: mass of the dried product, [kg]
and wherein said estimating initial conditions (CS,0) and kinetic constants (k0, k1, k2, . . . , kj), at time t=tj, is made by means of the following equations:


DRexp,0=DRtheor,0=−k0(CS,0−CS,eq,0)  (eq. 27)



DRexp,1=DRtheor,1=−k1{CS,0e−k1(t1−t0)++k1CS,eq,1[t1−t0e−k1(t1−t0)]−CS,eq,1}  (eq. 28)



DRexp,2=DRtheor,2=−k2{CS,1e−k2(t2−t1)++k2CS,eq,2[t2−t1e−k2(t2−t1)]−CS,eq,2}  (eq. 29)



DRexp,j=DRtheor,j=−kj{CSj−1e−kj(tj−tj−1)++kjCS,eq,j[tj−tj−1e−kj(tj−tj−1)]−CS,eq,j}(eq. 29ter)


and solving the following non-linear least square problem:





  
    
      
        
          min
          
            
              C
              
                S
                ,
                0
              
            
            ,
            
              k
              i
            
          
        
        ⁢
        
          
            ∑
            
              i
              =
              0
            
            j
          
          ⁢
          
            
              
                (
                
                  
                    D
                    ⁢
                    
                      
                    
                    ⁢
                    
                      R
                      
                        exp
                        ,
                        i
                      
                    
                  
                  -
                  
                    D
                    ⁢
                    
                      
                    
                    ⁢
                    
                      R
                      
                        theor
                        ,
                        i
                      
                    
                  
                
                )
              
              2
            
            .
          
        
      
    
    
      
        (
        
          
            eq
            .
            
              
            
            ⁢
            30
          
          ⁢
          bis
        
        )

17. Method according to claim 16, wherein said final time (tf) is calculated, assuming that a temperature of said product does not change, by means of the following equation, resulted from (eq. 21):


CS,f=CS,je−kj(tf−tj)++kjCS,eq,j[tf−tje−kj(tf−tj)]  (eq. 31bis)

where:
CS,f: final residual moisture concentration [% water and/or solvent over dried product];
CS,j: residual moisture concentration at time t=tj [% water and/or solvent over dried product].