1. Field of the Invention
The present invention relates generally to differential scanning calorimeters (DSCs), and more specifically to methods of accounting for heat leakage in DSCs.
2. Background Information
A DSC is a symmetrical instrument, comprising a sample and reference calorimeter within a common thermal enclosure, where the two calorimeters are intended to be identical. DSCs have a sensor which measures the temperature difference between the sample and the reference position in the respective calorimeters. A sample to be analyzed is loaded into a pan and placed on the sample position of the sensor and an inert reference material is loaded into a pan and placed on the reference position of the sensor (alternatively, an empty pan is often used as the reference). The sensor is installed in an oven whose temperature is varied dynamically according to a desired temperature program. The temperature program for conventional DSCs typically includes combinations of linear temperature ramps and constant temperature segments. Modulated DSCs use a temperature program in which periodic temperature modulations are superimposed on linear ramps. Modulated DSCs are described in U.S. Pat. No. 5,224,775, which is incorporated by reference herein.
Conventional DSC Heat Flow Rate Measurement
In conventional DSCs, the measured quantities are the difference between the temperatures of the sample and reference sensors and the sample sensor temperature. During the dynamic portion of the DSC experiment, a differential temperature is created between the sample and reference positions on the sensor. The temperature difference is the result of the difference between the heat flow to the sample and the heat flow to the reference. The temperature difference is assumed to be proportional to the difference in heat flow to the sample as compared to the heat flow to the reference, thus the differential temperature may be used to measure the heat flow to the sample using the equation
q=E(T)ΔT
Where, E(T) is a temperature dependent factor that reflects the proportionality of the measured differential temperature to the sample heat flow rate. A fundamental underlying assumption is that the sample and reference sides of the calorimeter are perfectly symmetrical. In reality the balance is less than perfect which is reflected in the observation that when the calorimeter is operated empty, the resulting heat flow rate signal is not zero as it should be if the DSC were symmetrical as assumed. Additionally, it is well known that the heat flow rate signal is smeared, as a result of the heat capacity of the sensor and pans and the differences in heating rate that exist between the sample and reference. The beating rate differences exist in conventional DSC during a transition or generally during a MDSC experiment. However, those shortcomings may be overcome by using the DSC apparatus and heat flow rate measurement methods disclosed below.
DSC Heat Flow Rate Measurement Including Sensor Asymmetry
Performing a heat flow balance on the sample side of the sensor yields the measured sample heat flow rate
which is equal to the heat flow rate through the sensor sample thermal resistance Rs minus the heat stored in Cs. Similarly, a heat balance on the reference side of the sensor gives the measured reference heat flow rate
which is equal to the heat flow rate through sensor reference thermal resistance Rr minus the heat stored in Cr. In the equations herein, τ represents time.
The desired quantity (the differential heat flow to the sample with respect to the reference) is the difference between the sample and reference heat flows:
q=qs−qr
Substituting for qs and qr yields:
Substituting the following expressions for two temperature differences in the differential scanning calorimeter,
ΔT=Ts−Tr
ΔTo=To−Ts
where ΔT is the temperature difference between the sensor sample and the reference positions and ΔTo is the temperature difference between the sample position and a position at the base of the sensor, results in the DSC heat flow equation:
The DSC heat flow equation has four terms. The first term accounts for the effect of the difference between the sensor sample thermal resistance and the sensor reference thermal resistance. The second term is the conventional DSC heat flow. The third term accounts for the effect of the difference between the sensor sample thermal capacitance and the sensor reference thermal capacitance. The fourth term accounts for the effect of the difference between the heating rates of the sample and reference sides of the DSC. Conventionally, when this equation is applied to the DSC heat flow, the first and third terms are zero because Rs and Rr are assumed to be equal and Cs and Cr are also assumed to be equal. The fourth term is zero because the heating rate difference between the sample and reference is ignored.
Heat Leakage Effects
Heat leakage is defined as heat that flows between the sample or reference and the enclosure without being measured. This leakage, or unmeasured heat flow, contributes to the uncertainty of experiments in which DSCs are utilized.
Advantages that are commonly associated with using twin calorimeters in DSCs as described above include cancellation of heat leakage and temperature disturbances common to both calorimeters. If the calorimeters each contain empty sample pans, and are heated at the same rate, there should be no difference between the sensor temperatures of each enclosure. However, in practice, the two calorimeters will always generally measure a different heat flow rate. Even if the calorimeters are perfectly matched, the presence of a sample causes the sample and reference heat flow rates to be imbalanced.
When a DSC is not balanced, the leakage heat flows will not be balanced and will not cancel when the difference is taken between the measured sample and reference heat flow rates, which contributes to a measurement error in the heat flow rate. In general, the leakage heat flows are a small fraction of the measured heat flow. If the imbalance is small, the measurement error is small, as a difference between two relatively small quantities. The magnitude of the errors has been estimated for disk-type heat flux DSCs at between 1 and 5% by numerical simulation.
In principle, the uncertainty can be reduced by calibration of the heat flow rate. However, experimental conditions that are different from the calibration conditions will cause uncertainty due to heat leakage to increase. Differences in sample heat capacity and differences in the transitions within a sample are the primary differences between calibration and experimental conditions, because the experimental sample in general will have a different heat capacity from the calibrant and will have different transitions.
A manifestation of the effect of heat leakage that is well known in differential scanning calorimetry is the difference between heat flow calibration and peak area calibration. In heat flow calibration, a sample having known specific heat capacity, such as sapphire, is heated at a constant heating rate and the heat flow rate is measured. Heat capacity of the sample is obtained by dividing the heat flow rate signal by the heating rate, giving the measured heat capacity. The ratio of the actual value to the measured value gives the heat flow calibration factor for the range of temperatures of the experiment. It is applied to the results of subsequent experiments to improve the accuracy of the heat flow measurement. In peak area calibration, the sample is a material with a transition whose enthalpy is well known and highly repeatable, such as the melt of a pure metal like indium. A calibration experiment is performed and the area of the resulting heat flow rate peak is measured by integrating the heat flow signal versus time over an appropriate baseline. The area of the peak provides the measured enthalpy of the transition. The ratio of the measured enthalpy to the actual enthalpy is the peak area calibration factor.
In principle, at the temperature of the transition, the peak area and the heat flow calibration factors should be identical, or very nearly so, but in general the agreement is rather poor. The reason for this disagreement is the large temperature imbalance created by the transition and the attendant increase in sample heat leakage due to the imbalance. For example, when calibrating with a typical indium sample of 5 mg at 10° C. min, the temperature difference between the sample and the reference before the melt is of the order of only a few tenths of 1° C., but that increases to 1° C. or more at the peak of the transition. Thus, the temperatures of the sample and reference calorimeters are very far from being identical, the twin principle is severely violated, and the leakage heat flows will be much different and will not cancel each other out.
Another example of the effect of heat leakage is the well known effect of thermal conductivity of the purge gas on DSC cell calibration. In either heat flow rate or peak area calibration, it is well known that changing the purge gas causes both the peak area and the heat flow rate calibration to change. In the absence of heat leakage, both the heat flow rate calibration and peak area calibration would depend only upon thermal resistance of the calorimeters and would be independent of the purge gas used.
The present invention is a method and system for calculating a heat flow to a sample in a differential scanning calorimeter (DSC). The DSC has a sensor within an enclosure comprising an absolute temperature measurement detector for measuring the temperature of a base position on the sensor, a first differential temperature detector for measuring the temperature difference between a sample position and the base position, and a second differential temperature detector for measuring the temperature difference between a reference position and a sample position. Thermal resistances and heat capacities of the DSC are calibrated. The DSC is operated, and the heat flow to the sample is calculated using a method that includes the leakage heat flows.
The present invention is also a method and system for calculating heat flow in a differential scanning calorimeter having an independent sample measuring section and an independent reference measuring section. The heat flow rate for the sample being analyzed is determined using the measured quantities in a method that accounts for sample pan heat flow, reference pan heat flow, and heat leakage. The heat flow to a sample in the sample pan is measured in a measuring section and the heat flow to a reference is measured in a reference measuring section.
In a preferred embodiment, the present invention comprises a two step method for calibrating thermal resistances and heat capacities of a differential scanning calorimeter (DSC) having sample and reference calorimeters. in the first step, the DSC is operated empty, at a constant heating rate. A correction factor for leakage heat flow is applied to the time constants of the two calorimeters. Effective time constants are determined for the sample and reference calorimeters, accounting for heat leakage.
In the second calibration step, the DSC is operated with known calibration samples placed without pans on both of the sample and reference calorimeters of the DSC. The heat capacity is determined for each of the sample and reference sensors, accounting for heat leakage, and the sample and reference thermal resistances are calculated using the time constants from the first step and the respective heat capacities.
In performing conventional differential scanning calorimetry experiments, it is assumed that the temperature difference between a sample and reference is proportional to the heat flow to the sample. However, the temperature difference is also influenced by the heat capacity of the pan, the sensing apparatus, and resistances to heat flow within certain components of the DSC. Therefore, as described in the other patents and pending applications that are incorporated by reference, methods can be performed for compensating for these effects through calibration and calculation to increase the accuracy of measurements in DSC experiments. Although these methods and systems significantly improve the instrument baseline heat flow rate, sensitivity and resolution, there has remained a problem that such measurements do not account for heat leakage. The following provides a description of a method and system for further improving the accuracy of the measurement of heat flow rates in a DSC by including the leakage heat flows in calculations of sample heat flow rate from measured signals.
DSC Heat Flow Rate Measurement Including Sensor Asymmetry and the Effects of the Sample and Reference Pans
As noted in the Background with reference to
And the measured reference heat flow rate is just that of the empty reference pan:
Where mps and mpr are the masses of the sample and reference pans, Cpan is the specific heat of the pan material.
Perform a heat balance on the sample pan, substitute for the sample pan heat flow and solve for the sample heat flow, giving:
Solve the reference heat flow equation for the pan specific heat and substitute the resulting equation for pan specific heat into the sample heat flow equation:
This equation gives the actual sample heat flow, accounting for the sample and reference pan heat flow rates. The second term on the right hand side is seen to be the reference heat flow rate multiplied by the ratio of the sample and reference pan masses and by the ratio of the sample and reference pan heating rates. It accounts for the fact that during a transition, or during a modulated DSC experiment, the sample pan heats at a different rate than the reference pan because of the sample heat flow.
When the sample heating rate is different from the reference heating rate, the fraction of reference heat flow subtracted from the sample heat flow is greater or lesser, depending upon whether the sample pan heating rate is greater or lesser than that of the reference pan. Because the reference heat flow is just the reference pan heat flow, this equation accounts for differences between sample and reference pan heating rates. For example, during a melt in a DSC, the heating rate of the sample pan falls below the programmed rate, while the reference pan continues to heat at the programmed rate. In a conventional DSC the reference heat flow rate being subtracted from the sample heat flow rate is for a pan heating at the programmed rate. Thus, during the melt, the reference heat flow rate being subtracted from the sample heat flow rate is too large and the resultant DSC heat flow rate signal is too small. During the baseline return, the sample pan is heating faster than the reference pan and the reference heat flow rate being subtracted from the sample heat flow rate is too small and as a consequence, the DSC heat flow rate signal is too large.
To use the true sample heat flow equation requires that the sample and reference pan temperatures be known so that their derivatives may be taken. Unfortunately, there is no way to measure the pan temperatures directly. However, the pan temperatures can be calculated from the temperature and heat flow signals as follows.
The equations for heat flow between the sensor and the sample and reference pans are:
Solving for the pan temperatures:
Tps=Ts−qs·Rps Tpr=Tr−qr·Rpr
Using these equations, pan temperatures and sample heat flows can be calculated from the measured signals.
The pan thermal resistances Rps and Rpr depend on the pan configuration, the purge gas used in the DSC and the temperature of the DSC. Several techniques for determining Rps and Rpr have been developed in the past and are well known to one of ordinary skill in the art. For example, one well-known technique is to measure the slope of the onset of a metal melt.
A preferred semi-empirical method for determining Rps and Rpr is to use a model equation for heat exchange between the sample and reference pans and the DSC sensor. Because when two nominally flat surfaces are brought into contact, they only make contact over a small fraction of the projected area of contact (because the surfaces are not perfectly flat), heat exchange could occur in principle through three parallel paths: by direct solid heat conduction through the asperities of the two surfaces in contact, by convection or conduction through the interstitial gas between the surfaces and by radiation between the two surfaces. However, except in the case of extremely high beat flux and a large temperature difference between the two surfaces, radiation is insignificant, and the dominant heat exchange through the gas is by conduction. In DSC, the heat flux and the resulting temperature difference between a pan and the DSC is small and the thermal contact resistance between the surfaces can be modeled by heat conduction through two series connected solid conductors (representing each of the surfaces) in parallel with heat conduction through the gas. Using these assumptions to model the DSC pan/sensor heat exchange, the model equation for the DSC pan contact resistance is:
Where, R(T) is the contact resistance as a function of temperature; kp(T), ks(T), kg(T) are thermal conductivities of the pan, sensor and gas; αp, αs and αg are geometric factors for the pan, sensor and gas. The geometric factors have the dimension of length, they may be thought of as the ratio of the area normal to the heat flow to the length of the heat conduction path. Thermal conductivities of the sensor, pan and gas are known.
The geometric factors are determined empirically as follows. The contact resistance between the pan and the sample may be measured using MDSC with average heating rate equal to zero. A sample of known heat capacity is loaded into a pan and placed on the sample position of the DSC, and the temperature is modulated sinusoidally at fixed amplitude about a given temperature. Assuming the sample and pan temperatures are the same, a two-temperature model of the DSC may be solved to find the time constant of the pan and sample.
Where temperature with an over bar means the modulated temperature amplitude, obtained by the MDSC temperature deconvolution algorithm (as described in U.S. Pat. No. 5,224,775, which is incorporated by reference herein), Css is the combined heat capacity of the sample and the pan and ω is the circular frequency of the modulation. The assumption that sample and pan temperatures are the same is justified when the modulation period is sufficiently long (typically 60 seconds or above). A similar set of equations can be written for the reference side of the DSC, and the two contact resistances may be determined in a single experiment.
Using this method for determining the contact resistance at a variety of discrete temperatures, the model equation above may be fitted to the data to determine the geometric parameters αp, αs and αg. Using a multiplicity of different samples and pans, the experiments may be repeated many times and a statistically determined average value of the contact resistance can be computed.
DSC Heat Flow Rate Measurement, Including Sensor Asymmetry, the Effects of Sample and Reference Pans and Heat Leakages
Addition of these two resistance elements requires modifications to the equations used to determine the actual sample heat flow. The measured sample and reference heat flows, qs and qr remain the same as described above with reference to
The heat exchange between the pans, sensor and DSC enclosure are considered separately. For the sample side we find,
For the reference side we find,
Solving the reference heat flow equation for the pan heat capacity and eliminating it from the sample heat flow equation and solving it for qsam:
Thus, it is seen that the sample heat flow equation including the heat flow leakages differs from the equation of the previous invention (i) by factors multiplying the measured heat flows, and (ii) by an additional term for both the sample and the reference sides. The factors and additional terms tend to increase the magnitude of the measured sample heat flow rate. The form of the equation makes sense intuitively, as leakage heat flows tend to decrease the measured heat flow.
The leakage heat flow consists of heat that is exchanged directly between the sample and reference pans and the DSC enclosure by the combination of radiation and heat transfer through the gas within the enclosure. The leakage thermal resistance is much larger than the pan contact thermal resistance, for example one group of researchers found that near 156° C., the leakage resistance was 784° C./watt while the pan contact resistance was 10° C./watt.
Calibration of a DSC that Measures the Heat Flow Rate Including Sensor Asymmetry and the Effects of the Sample and Reference Pans
The following describes a method for calibrating thermal resistances and heat capacities of a DSC that are necessary to measure heat flow rate, for a DSC according to the thermal model of
The first experiment is performed with an empty DSC cell. The DSC cell is first held at an isothermal temperature that is below the temperature range of the calibration, for a time segment sufficient to ensure complete equilibration of the sensor. The DSC cell is then heated at a constant heating rate to a temperature above the temperature range of the calibration, and then held at that temperature for another isothermal segment, for a time segment sufficient to ensure equilibration of the sensor at that temperature. This first experiment is used to calculate the sample and reference time constants as a function of temperature over the calibrated temperature range.
The measured sample and reference heat flow rates should be zero (since the DSC cell is empty). Accordingly, if qs and qr are set equal to zero in the heat balance equations for the sample and reference sides of the sensor, the time constants for the sample and reference are given by:
respectively. These results are stored as a function of temperature.
The second experiment uses a pair of calibration samples. The calibration samples may have the same mass, or may have different masses. Preferably, the calibration samples are monolithic sapphire samples, preferably weighing 25 mg or more. Other reference materials with well-known specific heats and no transitions in the temperature range of the calibration may be used instead of sapphire (in which case Csapph would be replaced in the following equations by Cmat where Cmat is the specific heat of the other reference material).
The sample and reference heat flows from the heat balance equations are set as follows:
where ms , mr are the masses of the sample and reference sapphires, respectively, Csapph is the specific heat of sapphire and Tss and Trs are the temperatures of the sample and reference sapphire.
Assume:
Substituting for qs and τs in the sample heat balance equation and solving for the sensor sample thermal capacitance Cs:
Substituting for qr and τr in the reference heat balance equation and solving for the sensor reference thermal capacitance Cr:
The results from the second experiment using sapphire (or another well-known calibration material) and using the time constants for a DSC cell obtained in the first experiment are used as shown above to calculate the sample and reference sensor heat capacities as a function of temperature. Finally, the sensor sample and reference thermal resistances are calculated from the time constants and the sensor thermal capacitances:
Calibration of a DSC that Measures the Heat Flow Rate Including Sensor Asymmetry, the Effects of the Sample and Reference Pans and Heat Leakages
In accordance with an embodiment of the present invention, the calibration procedure described above for calibrating thermal resistances and heat capacities of a DSC can be modified to include the heat leakages.
The first of the two calibration experiments involves running the DSC empty at a constant heating rate, but because of leakage heat flow, the measured heat flow is not zero. In the model of the empty DSC shown in
The effective time constant of the two calorimeters can be found as before:
The leakage thermal resistances are assumed to be known, while the magnitude of the measuring resistances is generally known to within 5% before calibration. This allows a correction factor to be applied to the time constant as follows. The ratio of the measuring thermal resistance to the effective measuring resistance is:
This time constant correction factor varies with temperature and with the type of purge gas used. Also, in this equation the subscript m is a generic subscript representing either the subscript s or the subscript r. The calorimeter time constants corrected for leakage heat flow are:
τs=Kτ′s τr=Kτ′r
Because the leakage thermal resistance is much larger than the sensor thermal resistances, the time constant correction factor is a number that is close to unity and is fairly insensitive to typical variations in the sensor thermal resistances (Rs and Rr). Thus, the use of a standard value of the correction factor results in relatively small uncertainty. However even that uncertainty may be reduced by performing an iterative analysis of the calculation data where the sensor thermal resistances are calculated and used to calculate a new correction factor which is used to calculate yet another value of the sensor thermal resistances, iteration continues until the changes in sensor thermal resistance diminish sufficiently.
The second calibration experiment is identical to the first, except that calibration sample without pans are placed on both the sample and the reference sides of the calorimeter. The sample heat flow equation given above with reference to
Substituting the sample heat flow measurement equation, assuming that the heating rate of the calibration sample and the sensor are the same, and solving for the sample sensor heat capacity:
And similarly, the reference sensor heat capacity is:
where Css and Csr are the heat capacities of the sample and reference calibration samples respectively, and Rls and Rlr are the leakage thermal resistances and Rss and Rsr are the contact thermal resistances of the sample and reference calibration samples, respectively.
The sample and reference thermal resistances are then calculated from the time constants and thermal resistances:
With the calibration complete, heat flow measurements may be made that includes the leakage heat flows.
Determining Leakage Thermal Resistances
To implement the measurement and calibration methods described above requires that the leakage thermal resistances be known. There are in fact three different leakage thermal resistances required to perform a calibration of the instrument and the subsequent measurement. The first is the leakage resistance of an empty DSC cell as used in the first step of the calibration; it represents the combination of radiation and gas heat transfer between the DSC sensor and the enclosure. The second is the leakage thermal resistance between the calibration samples and the enclosure; it represents the combination of radiation and gas heat transfer between each calibration sample and the enclosure. The third is the leakage thermal resistance between a sample pan and the DSC enclosure; it represents the combination of radiation and gas heat transfer between each pan and the enclosure. The values of leakage thermal resistances can be computed using Finite Element Analysis of the DSC cell including the interior of the enclosure, the sample pans and the sensor.
Finite element modeling is a well-established engineering tool wherein a component of arbitrary shape is considered to be comprised of a series of small elements of finite size. Using commercially or publicly available software, a set of equations is devised that relates the flow of heat through each of the elements to temperatures at the element's nodes. The resulting set of equations is solved to find the temperature distribution within the model from which the heat flow rates between surfaces of interest may be calculated. Simulated experiments are then run to determine the heat flows and to use the temperature differences and the rates of heat flowing to calculate equivalent leakage thermal resistances. Finite element analysis can be used to create values for the DSC leakage thermal resistances at different temperatures, to be coded into the instrument as software. A look-up table can then be used to identify a correct value of resistance for a particular temperature. The values change at different temperatures because the thermal conductivity of all of the materials that the DSC is constructed of changes, such that they are functions of temperature. These results will be used to create tables of leakage resistances as a function of temperature for empty DSC cells for each pan type and using selected purge gases. Radiation and gas conduction will both be included.
Apparatus embodying the present invention are shown schematically in
The foregoing disclosure of embodiments of the present invention has been presented for purposes of illustration and description. It is not intended to be exhaustive or to limit the invention to the precise forms disclosed. Many variations and modifications of the embodiments described herein will be obvious to one of ordinary skill in the art in light of the above disclosure. The scope of the invention is to be defined only by the claims appended hereto, and by their equivalents.
Further, in describing representative embodiments of the present invention, the specification may have presented the method and/or process of the present invention as a particular sequence of steps. However, to the extent that the method or process does not rely on the particular order of steps set forth herein, the method or process should not be limited to the particular sequence of steps described. As one of ordinary skill in the art would appreciate, other sequences of steps may be possible. Therefore, the particular order of the steps set forth in the specification should not be construed as limitations on the claims. In addition, the claims directed to the method and/or process of the present invention should not be limited to the performance of their steps in the order written, and one skilled in the art can readily appreciate that the sequences may be varied and still remain within the spirit and scope of the present invention.
The following symbols are all defined above, but are listed here for convenience when reading and applying the claims:
The present application claims priority to U.S. Provisional Patent Application No. 60/362,546, filed March 8, 2002, and is a continuation of Ser. No. 10/382,580, filed Mar. 7, 2003 now U.S. Pat. No. 6,843,595, which is a continuation-in-part of Ser. No. 09/769,313 now U.S. Pat. No. 6,561,692, filed Jan. 26, 2001. The present application incorporates by reference the following patent application: U.S. patent application Ser. No. 09/643,870, entitled “Heat Flux Differential Scanning Calorimeter,” filed Aug. 23, 2000; U.S. patent application Ser. No. 09/767,903 entitled “Differential Scanning Calorimeter,” filed Jan. 26, 2001; U.S. patent application Ser. No. 09/769,313, filed Jan. 26, 2001, entitled “Modulated Differential Scanning Calorimeter,” U.S. patent application Ser. No. 10/085,054, filed Mar. 1, 2002, entitled “System and Method for Calibrating Contact Thermal Resistances in Differential Scanning Calorimeters;” and U.S. patent application Ser. No. 10/382,580, entitled “Differential Scanning Calorimeter Accounting for Heat Leakage.”
Number | Name | Date | Kind |
---|---|---|---|
6071008 | Hatta et al. | Jun 2000 | A |
6488406 | Danley | Dec 2002 | B2 |
6497509 | Merzliakov et al. | Dec 2002 | B2 |
6530686 | Nakamura | Mar 2003 | B1 |
6572263 | Refalo et al. | Jun 2003 | B1 |
7025497 | Danley | Apr 2006 | B2 |
Number | Date | Country | |
---|---|---|---|
20060187998 A1 | Aug 2006 | US |
Number | Date | Country | |
---|---|---|---|
60362546 | Mar 2002 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 10382580 | Mar 2003 | US |
Child | 11400650 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 09769313 | Jan 2001 | US |
Child | 10382580 | US |