# Patent

## Title
Method of on-line monitoring of radial clearances in steam turbines

## Abstract

A method of monitoring radial clearances in a steam turbine during operation of the turbine is provided. The method, in an exemplary embodiment, includes measuring a temperature of the rotor shaft at a time1 and at a time2, measuring a temperature of the rotor blade at time1 and at time2, measuring a temperature of the shell at time1 and at time2, calculating a shaft radial growth between time1 and time2, calculating a blade growth between time1 and time2, calculating a shell radial growth between time1 and time2, and determining a change in a radial gap between the shell and a distal end of the rotor blade from time1 to time2 using the following equation: change in radial gap=shell radial growth−shaft radial growth−blade growth.

## Background
BACKGROUND OF THE INVENTION
The present invention relates generally to rotary machines, such as steam and gas turbines, and, more particularly, relates to a method of monitoring clearance between tips of rotating rotor blades and a stationary outer casing of a reaction design high pressure steam turbine.
Steam and gas turbines are used, among other purposes, to power electric generators. A steam turbine has a steam path which typically includes, in serial-flow relationship, a steam inlet, a turbine, and a steam outlet. A gas turbine has a gas path which typically includes, in serial-flow relationship, an air intake (or inlet), a compressor, a combustor, a turbine, and a gas outlet (or exhaust nozzle). Compressor and turbine sections include at least one circumferential row of rotating blades. The free ends or tips of the rotating blades are surrounded by a stator casing.
The efficiency of the turbine depends in part on the radial clearance or gap between the rotor blade tips and the surrounding casing and the clearance between the rotor and the diaphragm packings. If the clearance is too large, more of the steam or gas flow will leak through the gap between the rotor blade tips and the surrounding casing or between the diaphragm and the rotor, decreasing the turbine's efficiency. If the clearance is too small, the rotor blade tips can strike the surrounding casing during certain turbine operating conditions. Gas or steam leakage, either out of the gas or steam path or into the gas or steam path, from an area of higher pressure to an area of lower pressure, is generally undesirable. For example, gas-path leakage in the turbine or compressor area of a gas turbine, between the rotor of the turbine or compressor and the circumferentially surrounding turbine or compressor casing, will lower the efficiency of the gas turbine leading to increased fuel costs. Also, steam-path leakage in the turbine area of a steam turbine, between the rotor of the turbine and the circumferentially surrounding casing, will lower the efficiency of the steam turbine leading to increased fuel costs.
It is known that the clearance changes during periods of acceleration or deceleration due to changing centrifugal force on the blade tips and due to relative thermal growth between the rotating rotor and stationary casing. During periods of differential centrifugal and thermal growth of the rotor and casing the clearance changes can result in severe rubbing of the moving blade tips against the stationary casing. This increase in blade tip clearance results in efficiency loss.
Clearance control devices, such as rigid abradable shrouds, have been used in the past to accommodate rotor-to-casing clearance change. However, none are believed to represent an optimum design for controlling such clearance. Also, positive pressure packings have been used that include movable packings that permit the packings to be in a retracted position during startup and in an extended position during steady state operation of the turbine. However, the moving parts can stick during operation preventing the packings from moving between the extended and retracted positions.
BRIEF DESCRIPTION OF THE INVENTION
In one aspect, a method of monitoring radial clearances in a steam turbine during operation of the turbine is provided. The turbine includes an outer shell and a rotor including a rotor shaft and a plurality of rotor blades attached to the rotor shaft. The method includes measuring a temperature of the rotor shaft at a time1 and at a time2, measuring a temperature of the rotor blade at time1 and at time2, measuring a temperature of the shell at time1 and at time2, calculating a shaft radial growth between time1 and time2, calculating a blade growth between time1 and time2, calculating a shell radial growth between time1 and time2, and determining a change in a radial gap between the shell and a distal end of the rotor blade from time1 to time2 using the following equation: change in radial gap=shell radial growth−shaft radial growth−blade growth.
In another aspect, a method of monitoring radial clearances in a steam turbine during operation of the turbine is provided. The turbine includes an outer shell and a rotor including a rotor shaft and a plurality of rotor blades attached to the rotor shaft. The method includes measuring a temperature of the rotor shaft continuously during operation, measuring a temperature of the rotor blade continuously during operation, measuring a temperature of the shell continuously during operation, calculating a shaft radial growth as a function of rotor shaft temperature over time, calculating a blade growth as a function of rotor blade temperature over time, calculating a shell radial growth as a function of shell temperature over time, and determining a change in a radial gap between the shell and a distal end of the rotor blade over time using the following equation: change in radial gap=shell radial growth−shaft radial growth−blade growth.
In another aspect, a method of monitoring radial clearances in a steam turbine during operation of the turbine is provided. The turbine includes an outer shell and a rotor including a rotor shaft and a plurality of rotor blades attached to the rotor shaft. The method includes calculating a shaft radial growth as a function of rotor shaft temperature over time, calculating a blade growth as a function of rotor blade temperature over time, calculating a shell radial growth as a function of shell temperature over time, and determining a change in a radial gap between the shell and a distal end of the rotor blade over time using the following equation: change in radial gap=shell radial growth−shaft radial growth−blade growth.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a sectional schematic representation of a reaction design steam turbine.
FIG. 2 is an enlarged sectional schematic representation of a portion of the steam turbine shown in FIG. 1.
FIG. 3 is a flow chart of a method of monitoring radial clearances in a steam turbine during operation of the turbine.
FIG. 4 is a schematic representation of a portion of a steam turbine rotor.



DETAILED DESCRIPTION OF THE INVENTION
A method of monitoring radial clearances in a steam turbine during operation of the turbine is described in more detail below. The method calculates thermal expansions of components in the steam turbine which are proportional to averaged metal temperatures at a given location in the turbine. For example, the averaged temperature for the turbine shell at a given location can be obtained from measurements of shell temperature at one or more points across the thickness of the shell. Also, the temperature distribution in the rotor at a given location can be computed from the measured surface temperature and the rate of change of surface temperature over time. The method is described as used in a reaction design steam turbine; however, the method described below is applicable for other steam turbine designs, such as impulse steam turbines. The method uses the measured data from turbine shells and rotors for on-line computations of radial clearances of turbine components. This real time clearance data can be used by an operator to control turbine transients such that tip clearance changes are within specified limits for high thermal efficiency and to avoid rubbing between rotor tips and the shell.
Referring to the drawings, FIG. 1 is a sectional schematic-representation of a reaction design steam turbine 10. Steam turbine 10 includes a rotor shaft 12 passing through turbine 10 and sealed at each end by packings 14. A plurality of turbine blades 16 are connected to shaft 12. Between turbine blades 16 there is positioned a plurality of non-rotating turbine nozzles 18. Turbine blades or buckets 16 are connected to turbine shaft 12 while turbine nozzles 18 extend from an inner housing or shell 20 surrounding turbine blades 16 and nozzles 18. An outer housing 22 encloses inner housing 20 and rotor shaft 12 Steam is directed through nozzles 18 and through blades 16 causing blades 16 to rotate along with turbine shaft 12.
FIG. 2 is an enlarged sectional schematic representation of inner housing 20 of steam turbine 10. Inner housing 20 includes a plurality of outer ring portions 30. Each outer ring portions 30 include a ring 32 of steam directing nozzles 18 supported within outer ring portion 30, and an inner ring portion 34 contained within nozzle ring 32. Turbine buckets 16 are secured at their inner ends 36 to turbine wheels 38 extending from turbine shaft 12 rotatable about an axis 40. The radial outer ends 42 of buckets 16 include bucket covers 44 which rotate with buckets 16. In one embodiment, a cover 44 is positioned on radial outer end 42 of each bucket 16 and in alternate embodiments on outer ends 42 of two or more buckets 16 in the form of a band so as to permit adjacent buckets 16 to be coupled to a common cover or band.
Inner ring portion 34 of housing 20 includes a packing ring 48. Packing ring 48 is positioned adjacent turbine shaft 12. Turbine shaft 12 includes a sealing means 54 to seal a gap 56 between turbine shaft 12 and inner ring portion 34 of inner housing 20 to prevent the passage of stream through gap 56. Sealing means 54 is positioned adjacent packing ring 48 and includes a plurality of axially spaced brush seals 58 extending from rotor 12. Sealing means 54 can also include axially spaced labyrinth seal teeth (not shown) or a combination of axially spaced labyrinth seal teeth and brush seal seals 58.
Bucket covers 44 include a sealing means 66 to provide a seal in a gap 68 between bucket cover 44 and housing 20 to prevent the passage of steam through gap 68. Sealing means 66 includes a plurality of axially spaced labyrinth seal teeth 70 extending from bucket cover 44. Sealing means 66, in other embodiments include brush seals alone or combined with axially spaced labyrinth seal teeth 70.
Referring to FIG. 3, a method 80 of monitoring radial clearances, for example the size of gap 68 and the size of gap 56, in steam turbine 10 includes measuring 82 a temperature of rotor shaft 12 a time1 and at a time2, measuring 84 a temperature of rotor blade 16 time1 and at time2, measuring 86 a temperature of shell 20 at time1 and at time2, calculating 88 a radial growth of shaft 12 between time1 and time2, calculating 90 a growth of blade 16 between time1 and time2, calculating 92 a radial growth of shell 20 between time1 and time2, and determining 94 a change in radial gap 68 from time1 to time2 using the equation:

Change In Radial Gap=Shell Radial Growth−Shaft Radial Growth−Blade Growth.

Of course, because the radial growth of blade 16 has no effect on gap 56, the term Blade Growth in the above equation is zero for calculations of the change in radial growth of gap 56. Also, the temperatures of shaft 12, shell 20 and blade 16 can be measured at distinct intervals or can be continuously monitored over time.

The shaft radial growth is equal to the coefficient of thermal expansion of the rotor (αR) times an outer radius (RR) of the rotor times an instantaneous volume averaged temperature (TR) of the rotor.

Shaft Radial Growth=αR*RR*TR

The blade radial growth is equal to the coefficient of thermal expansion of the rotor blade (αB) times a length (LB) of the rotor blade times an instantaneous volume averaged temperature (TB) of the rotor blade.

Rotor Blade Growth=αB*LB*TB

In most cases the instantaneous volume averaged temperature (TB) of the rotor blade can be closely approximated by the rotor outer surface temperature or the steam temperature.
The shell radial growth is equal to the coefficient of thermal expansion of the shell (αS) times an inner radius (RS) of the shell at the blade tip times an instantaneous volume averaged temperature (TS) of the shell.

Shell Radial Growth=αS*RS*TS

In double wall shell designs with horizontal flanges, the radial clearance can vary as a function of circumferential location on the shell. To account for these variances shell radial growth is calculated at the top, the bottom and the side of the shell. Particualrly, the instantaneous volume averaged temperature (TS) of the shell is calculated for each location, at the top, the bottom and the side of the shell.
Instantaneous average temperatures TR and TS are computed using a finite difference method employing a finite element model utilizing the finite element of a segment of an infinitely long cylinder. This method is explained below using the rotor as an example, and the same method is applicable for the shell, considering the shell as a hollow cylinder. Referring to FIG. 4, the rotor is divided in a specific number of elements, for example 10 elements.
Controls

    Elements=10 (Elements number)
    Nodes=Elements+2 (Nodes number)
    Nr=Nodes−1 (Last Centroid Node Number) 



  Volume
  =
  
    
      
        (
        
          
            R
            o
            2
          
          -
          
            R
            i
            2
          
        
        )
      
      Elements
    
    ⁡
    
      [
      
        in
        3
      
      ]
    
  




    Volume→Element Volume
    The Temperature & Time Maximum Incremental Changes are set.
    MaxDTemp=5 (Maximum Incremental Temperature Change) 



  
    Max
    ⁢
    
       
    
    ⁢
    DTime
  
  =
  
    
      
        
          (
          
            
              R
              o
            
            -
            
              R
              i
            
          
          )
        
        2
      
      
        8
        ·
        DO
        ·
        
          Elements
          2
        
      
    
    ⁢
    
       
    
    [
    min
    ]
  




    MaxDTime→Maximum Incremental Time Change


Initializing Temperature
This block assigns an initial temperature value to each boundary of the rotor elements.

    
    
        For I=1 To Nodes
        Tr(I)=InitialTemp [° F.]
        Next I
    
    
    Tr(Nodes)→Array of Nodes elements.


Since Tr(Nodes) represents the boundary temperatures of each rotor element, all Tr elements have an initial value of InitialTemp.
Initial Temperature Distribution
Tsurf=InitialTemp: Tavg=InitialTemp: Tbore=InitialTemp
Since the initial temperature distribution of the rotor is uniform, the rotor surface temperature (Tsurf), the rotor bore temperature (Tbore) and the rotor average temperature (Tavg) have an initial value of InitialTemp.
Centrifugal Stress Factor
Due to the large stress gradient existing at the rotor bore, the BoreCntrfStrs (Bore Centrifugal Stress) needs to be evaluated. This block calculates the Centrifugal Stress Factor (SpeedFact) needed to calculate the Bore Centrifugal Stress defined as: 



  SpeedFact
  =
  
    
      InBoreCntrStrs
      
        RatedSpeed
        
          
             
          
          ⁢
          2
        
      
    
    ⁡
    
      [
      
        KSI
        /
        RPM2
      
      ]
    
  




Extrapolation Factor
To calculate the thermal stress at the surfaces at any given time, Tsurf, Tbore, and Tavg must be known. The extrapolation factor (Extrapfact) is found by means of an extrapolation of the ramp rate temperatures to the inner surface.
If we consider that Ri always has a greater value than 0, is possible to eliminate the second option that considers the case for Ri=0. 


  
    R
    =
    
      
        
          
            R
            i
            2
          
          +
          
            Volume
            2
          
        
      
      ⁡
      
        [
        in
        ]
      
    
  




  
    R2
    =
    
      
        
          
            R
            2
          
          +
          Volume
        
      
      ⁡
      
        [
        in
        ]
      
    
  




  
    Extrapfact
    =
    
      
        (
        
          
            
              R
              2
            
            -
            
              
                R
                i
                2
              
              ·
              
                (
                
                  1
                  +
                  
                    2
                    ·
                    
                      log
                      ⁡
                      
                        (
                        
                          R
                          
                            R
                            i
                          
                        
                        )
                      
                    
                  
                
                )
              
            
          
          
            Volume
            -
            
              2
              ·
              
                R
                i
                2
              
              ·
              
                log
                ⁡
                
                  (
                  
                    
                      R
                      2
                    
                    R
                  
                  )
                
              
            
          
        
        )
      
      ⁡
      
        [
        Dimensionless
        ]
      
    
  



ExtrapFact→Extrapolation Factor for Temperature at Bore Surface
Heatflow Factors
This block assigns a heatflow factor for the internal elements of the rotor. The heatflow factor for a specific rotor element is the average area normal to qi divided by the distance from element i to element i+1.

For I=2 To Nr−1

R2=Sqr(R^2+Volume)[in]

A(I)=(R+R2)/(R2−R)[in]

R=R2

Next I

Since the Finite Difference Method uses the conservation of energy for a specific rotor element to get the radial temperature distribution, the values of the heat flow factors vary radially too.
Surface & Bore Heatflow Factors
The previous formula to calculate the heatflow factors applies directly to all internal elements, but must be modified for the boundary elements to meet the boundary conditions. This block assigns the heatflow factor for the boundary elements.

A(1)=0[in]

A(Nr)=(R+Ro)/(Ro−R)[in]

Asurf=2*Ro/A(Nr)[in]

A(Nr)→Array of Nr elements.

Ramp Rates
This block defines the Time, Temperature, HTC and Speed Ramp Rates using the inputs of the previous block. The variable NumSubSteps defines the number of iterations of the new temperature distribution block.

DTime=Time−OldTime[min]

DTemp=Tfluid−OldTfluid[° F.]

TempRamp=DTemp/DTime[° F./min]

HTCRamp=(HTC−OldHTC)/DTime[BTU/hr·ft2·° F./min]

SpeedRamp=(Speed−OldSpeed)/DTime[RPM/min]

 NumSubSteps=Int(Application.WorksheetFunction—

.Max(DTime/MaxDTime, Abs(DTemp/MaxDTemp), 1))[Dimensionless]

DT=DTime/NumSubSteps[min]

The accuracy of the calculation of the new temperature distribution calculation depends on the size of the time step. Sufficient accuracy is obtained if the maximum time step is DT.
New Temperature Distribution
This block calculates the new temperature distribution of the rotor setting new values for Tr elements.

For K=1 To NumSubSteps

TrSum=0

For I=2 To Nr

DF=Diff(D0, DM, Tr(I))[in2/min]

TrNew=Tr(I)−DT*DF/Volume*(A(I)*(Tr(I)−Tr(I+1))+A(I−1)*(Tr(I)−Tr(I−1)))[° F.]

TrSum=TrSum+TrNew[° F.]

Tr(I−1)=PrevNew[° F.]

PrevNew=TrNew [° F.]

Next I

The variable TrNew calculates the new temperature value for a specific rotor element using the corresponding values of thermal diffusivity and heatflow factor. The variable TrSum allows the required storage of information to calculate the average temperature (Tavg).
Time Delta
This block assigns new values for Time, Temperature and Speed using the maximum time step (DT).

Time=OldTime+K*DT[min]

Tfluid=OldTfluid+TempRamp*K*DT[° F.]

Speed=OldSpeed+SpeedRamp*K*DT[RPM]

Tr(Nr)=TrNew [° F.]

Surface Temperature & HTC
This block assigns a new value for HTC and calculates the surface temperature (Tsurf). Since a convection heat transfer process is carried out between the rotor surface and the fluid, the fluid temperature (Tfluid) is required to calculate the temperature at the last rotor element Tr(Nodes).

HTC=OldHTC+HTCRamp*K*DT[BTU/hr·ft2·° F.]

Cond=RhoC*Diff(D0, DM, Tr(Nodes))[BTU·° F./in·min]  (A)

Factor=(HTC/8640)*Asurf/Cond[Dimensionless]

Tr(Nodes)=(Tr(Nr)+Factor*Tfluid)/(1+Factor)[° F.]

Tsurf=Tr(Nodes)[° F.]

Bore & Average Temperatures
This block calculates the bore temperature (Tbore) and the average temperature (Tavg). To calculate the temperature of the first rotor element the extrapolation factor for temperature at bore surface (ExtrapFact) is required.

Tr(1)=Tr(2)−ExtrapFact*(Tr(3)−Tr(2))[° F.]

Tbore=Tr(1)[° F.]

Tavg=TrSum/Elements[° F.]

Surface Stress & Strain
This block calculates the surface stress and strain. The actual coefficient of thermal expansion (Alpha) is required.

ExpnC=Alpha(A0, AM, Tavg)[%/° F.]

SurfStrn=ExpnC*(Tavg−Tsurf)[%]

SurfStrs=Modulus(E0, EM, Tsurf)*SurfStm/0.7[KSI]

PCSurfAllow=100*SurfStrn/AllowSurfStrn[%]

Bore Stress
This block calculates the total bore stress. The actual coefficient of thermal expansion (Alpha) and the actual Young's modulus (Modulus) are required.
 BoreCntrfStrs=SpeedFact*Speed^2[KSI]

BoreThrmStrs=Modulus(E0, EM, Tbore)*ExpnC*(Tavg−Tbore)/0.7[KSI]

TotBoreStrs=BoreThrmStrs+BoreCntrfStrs[KSI]

PCBoreAllow=100*TotBoreStrs/AllowBoreStrs[%]

Next K

While the invention has been described in terms of various specific embodiments, those skilled in the art will recognize that the invention can be practiced with modification within the spirit and scope of the claims.

## Claims

1. A method of monitoring radial clearances in a steam turbine during operation of the turbine, the turbine comprising an outer shell and a rotor, the rotor comprising a rotor shaft and a plurality of rotor blades attached to the rotor shaft, said method comprising:
measuring a temperature of the rotor shaft at a time1 and at a time2; 
measuring a temperature of at least one of the plurality of rotor blades at time1 and at time2; 
measuring a temperature of the shell at time1 and at time2; 
calculating a shaft radial growth between time1 and time2; 
calculating a blade growth between time1 and time2; 
calculating a shell radial growth between time1 and time2; and 
determining a change in a radial gap between the shell and a distal end of the rotor blade from time1 to time2 using the following equation: 


change in radial gap=shell radial growth−shaft radial growth−blade growth, 


where the growth calculations are performed by using a finite difference method to obtain the instantaneous volume averaged temperatures of the rotor and the shell, and wherein at least one of the surrounding gas temperature and the measured blade temperatures is used to approximate the instantaneous averaged temperature of the plurality of rotor blades.

2. A method in accordance with claim 1 wherein calculating a shell radial growth comprises calculating a shell radial growth using the following equation:


shell radial growth=αS*RS*TS 


where
αS is the coefficient of thermal expansion of the shell; 
RS is an inner radius of the shell at the blade tip; 
TS is an instantaneous volume averaged temperature of the shell.

3. A method in accordance with claim 2 wherein TS is an instantaneous volume averaged temperature of the shell at a top location.

4. A method in accordance with claim 2 wherein TS is an instantaneous volume averaged temperature of the shell at a bottom location.

5. A method in accordance with claim 2 wherein TS is an instantaneous volume averaged temperature of the shell at a side location.

6. A method in accordance with claim 1 wherein calculating a shaft radial growth comprises calculating a shaft radial growth using the following equation:


shaft radial growth=αR*RR*TR 


where
αR is the coefficient of thermal expansion of the rotor; 
RR is an outer radius of the rotor; 
TR is an instantaneous volume averaged temperature of the rotor.

7. A method in accordance with claim 1 wherein calculating a rotor blade growth comprises calculating a rotor blade growth using the following equation:


rotor blade growth=αB*RB*TB 


where
αB is the coefficient of thermal expansion of the blade; 
RB is a length of the blade; 
TB is an instantaneous volume averaged temperature of the blade.

8. A method of monitoring radial clearances in a steam turbine during operation of the turbine, the turbine comprising an outer shell and a rotor, the rotor comprising a rotor shaft and a plurality of rotor blades attached to the rotor shaft, said method comprising:
measuring a temperature of the rotor shaft continuously during operation; 
measuring a temperature of at least one of the plurality of rotor blades continuously during operation; 
measuring a temperature of the shell continuously during operation; 
calculating a shaft radial growth as a function of rotor shaft temperature over time; 
calculating a blade growth as a function of rotor blade temperature over time; 
calculating a shell radial growth as a function of shell temperature over time; and 
determining a change in a radial gap between the shell and a distal end of the rotor blade over time using the following equation: 
change in radial gap=shell radial growth−shaft radial growth−blade growth. 
where the growth calculations are performed by using a finite difference method to obtain the instantaneous volume averaged temperatures of the rotor and the shell, and wherein at least one of the surrounding gas temperature and the measured blade temperatures is used to approximate the instantaneous averaged temperature of the plurality of rotor blades.

9. A method in accordance with claim 8 wherein calculating a shell radial growth comprises calculating a shell radial growth using the following equation:


shell radial growth=αS*RS*TS 


where
αS is the coefficient of thermal expansion of the shell; 
RS is an inner radius of the shell at the blade tip; 
TS is an instantaneous volume averaged temperature of the shell.

10. A method in accordance with claim 9 wherein TS is an instantaneous volume averaged temperature of the shell at a top location.

11. A method in accordance with claim 9 wherein TS is an instantaneous volume averaged temperature of the shell at a bottom location.

12. A method in accordance with claim 9 wherein TS is an instantaneous volume averaged temperature of the shell at a side location.

13. A method in accordance with claim 8 wherein calculating a shaft radial growth comprises calculating a shaft radial growth using the following equation:


shaft radial growth=αR*RR*TR 


where
αR is the coefficient of thermal expansion of the rotor; 
RR is an outer radius of the rotor; 
TR is an instantaneous volume averaged temperature of the rotor.

14. A method in accordance with claim 8 wherein said calculating a rotor blade growth comprises calculating a rotor blade growth using the following equation:


rotor blade growth=αB*LB*TB 


where
αB is the coefficient of thermal expansion of the blade; 
LB is a length of the blade; 
TB is an instantaneous volume averaged temperature of the blade.

15. A method of monitoring radial clearances in a steam turbine during operation of the turbine, the turbine comprising an outer shell and a rotor, the rotor comprising a rotor shaft and a plurality of rotor blades attached to the rotor shaft, said method comprising:
calculating a shaft radial growth as a function of rotor shaft temperature over time; 
calculating a blade growth as a function of at least one of the plurality of rotor blade temperature over time; 
calculating a shell radial growth as a function of shell temperature over time; and 
determining a change in a radial gap between the shell and a distal end of the rotor blade over time using the following equation: 
change in radial gap=shell radial growth−shaft radial growth−blade growth; 
where the growth calculations are performed by using a finite difference method to obtain the instantaneous volume averaged temperatures of the rotor and the shell, and wherein at least one of the surrounding gas temperature and the measured blade temperatures is used to approximate the instantaneous averaged temperature of the plurality of rotor blades.

16. A method in accordance with claim 15 wherein calculating a shell radial growth comprises calculating a shell radial growth using the following equation:


shell radial growth=αS*RS*TS 


where
αS is the coefficient of thermal expansion of the shell; 
RS is an inner radius of the shell at the blade tip; 
TS is an instantaneous volume averaged temperature of the shell.

17. A method in accordance with claim 16 wherein TS is an instantaneous volume averaged temperature of the shell at a top location.

18. A method in accordance with claim 16 wherein TS is an instantaneous volume averaged temperature of the shell at a bottom location.

19. A method in accordance with claim 16 wherein TS is an instantaneous volume averaged temperature of the shell at a side location.

20. A method in accordance with claim 15 wherein calculating a shaft radial growth comprises calculating a shaft radial growth using the following equation:


shaft radial growth=αR*RR*TR 


where
αR is the coefficient of thermal expansion of the rotor; 
RR is an outer radius of the rotor; 
TR is an instantaneous volume averaged temperature of the rotor.

21. A method in accordance with claim 15 wherein calculating a rotor blade growth comprises calculating a rotor blade growth using the following equation:


rotor blade growth=αB*LB*TB 


where
αB is the coefficient of thermal expansion of the blade; 
LB is a length of the blade; 
TB is an instantaneous volume averaged temperature of the blade.

22. A system for monitoring radial clearances in a steam turbine during operation of the turbine, the turbine comprising an outer shell and a rotor, the rotor comprising a rotor shaft and a plurality of rotor blades attached to the rotor shaft, said system comprising:
a measurement means configured to: 
measure a temperature of the rotor shaft at a time1 and at a time2; 
measure a temperature of at least one said plurality of rotor blades at time1 and at time2; 
measure a temperature of the shell at time1 and at time2; and 

a calculation means configured to: 
calculate a shaft radial growth between time1 and time2; 
calculate a blade growth between time1 and time2; 
calculate a shell radial growth between time1 and time2; and 
calculate a change in a radial gap between the shell and a distal end of the rotor blade from time1 to time2 using the following equation: 


change in radial gap=shell radial growth−shaft radial growth−blade growth, 



where the growth calculations are performed by using a finite difference method to obtain the instantaneous volume averaged temperatures of the rotor and the shell, and wherein at least one of the surrounding gas temperature and the measured blade temperatures is used to approximate the instantaneous averaged temperature of the plurality of rotor blades.

23. A system in accordance with claim 22 wherein said calculation means is further configured to calculate the shell radial growth using the following equation:


shell radial growth=αS*RS*TS 


where
αS is the coefficient of thermal expansion of the shell; 
RS is an inner radius of the shell at the blade tip; 
TS is an instantaneous volume averaged temperature of the shell.

24. A system in accordance with claim 23 wherein TS is an instantaneous volume averaged temperature of the shell at a top location.

25. A system in accordance with claim 23 wherein TS is an instantaneous volume averaged temperature of the shell at a bottom location.

26. A system in accordance with claim 23 wherein TS is an instantaneous volume averaged temperature of the shell at a side location.

27. A system in accordance with claim 22 wherein said calculation means is further configured to calculate the shaft radial growth using the following equation:


shaft radial growth=αR*RR*TR 


where
αR is the coefficient of thermal expansion of the rotor; 
RR is an outer radius of the rotor; 
TR is an instantaneous volume averaged temperature of the rotor.

28. A system in accordance with claim 22 wherein said calculation means is further configured to calculate the rotor blade growth using the following equation:


rotor blade growth=αB*LB*TB 


where
αB is the coefficient of thermal expansion of the blade; 
LB is a length of the blade; 
TB is an instantaneous volume averaged temperature of the blade.