BASED ON RQTQR VELOCITY

Background of the Invention

Field of the Invention The present invention relates to a centrifuge instrument having a system for automatically identifying a rotor introduced thereinto.

Description of the Prior Art A centrifuge instrument is a device adapted to expose a liquid sample carried in a rotating member, called a rotor, to a centrifugal force field. The centrifuge instrument includes a drive shaft, or spindle, adapted to receive any one of a predetermined plurality of rotors. It is important to correctly ascertain the identity of a particular rotor being used in the instrument at any given time. Such rotor identity information is important, among other reasons, for automatically controlling acceleration and deceleration times and for controlling the temperature or other parameters related to the particular separation being effected. Perhaps more importantly rotor identification is vital to insure that the particular rotor being used is not rotated to a speed that would cause rotor disintegration at an energy level high enough to breach the containment system of the instrument.

Presently rotor identification may be performed manually by requiring the operator of the instrument to introduce information via the control panel regarding the identity of the particular rotor being utilized. This system is open to inadvertent error or deliberate misrepresentation by the operator and thus cannot be relied upon for providing rotor identification information if the same is being used in connection with any safety- related consideration.

Automatic systems for rotor identification are available. Exemplary of such systems are those shown in United States Patent 4,551,715 (Durbin) and United States Patent 4,601,697 (Kamm). These systems utilize some form of coding elements usually disposed on the undersurface of the rotor. The coding elements are read by an appropriate optical or magnetic detector mounted in an operative location in the instrument. These systems share the disadvantage that the detector element, due to its location within the instrument, may be subject to corrosion which would vitiate its ability to accurately detect the coding elements provided on the rotor. Moreover such a system would be inapplicable in ascertaining the identity of rotors not equipped with the appropriate coding elements. Thus these identification systems would be unable to identify a significant population of rotors unless those rotors were retrofit with the appropriate coding elements. Furthermore retrofitting carries with it the risk of accidental or deliberate mismarking of the rotor and for this reason shares the same disadvantages as discussed above.

A rotor identification system relying on the interruption of a beam of light from a source to a detector is disclosed in United States Patent 4,450,391 (Hara).

United States Patent 4,827,197 (Giebeler) discloses a rotor identification system based on the inertia of the rotor when the rotor is used in what is believed to be an evacuated chamber. Such a system may be applicable for use in a nonevacuated or partially evacuated chamber so long as the inertia measurement is made at an angular velocity which is sufficiently low so that windage effects are negligible. This system would appear to become unreliable when windage effects become dominant.

An ultrasonic rotor recognition system is disclosed and claimed in copending application Serial Number 07/363,907, filed May 18, 1989, based on international application PCT/US87/03221 (Romanauskas), assigned to the assignee of the present invention.

Summary of the Invention

The present invention relates to an apparatus and to a method for identifying which one of a plurality of rotors is mounted within a centrifuge instrument. Each rotor has a predetermined velocity versus time profile associated therewith. The instrument includes a motive source having a shaft adapted to receive one of the plurality of rotors thereon.

In accordance with a first embodiment of the invention, one or more signal(s) is(are) generated representative of the actual velocity ω _a of a rotor disposed on the shaft at one or more measurement times t _m following initiation of rotation of the rotor. The predetermined measurement times t _m is(are) selected such that windage effects imposed on the rotor cause the velocity of the rotor on the shaft to differ by a measurable amount from the velocity of each of the others of the plurality of rotors. A rotor identity signal based upon the windage of the rotor is generated in response to the signal(s) representative of the velocity ω _a . The rotor identity signal may be generated using a look-up table or a comparator for comparing the actual velocity signal ω _a to a velocity reference signal ω _r ef- The velocity reference signal ω _r ef may be derived from a first rotor identification system.

In accordance with a second embodiment of the invention one or more signal(s) is(are) generated representative of the time t _a that a rotor disposed on the shaft first reaches one or more predetermined measurement velocities ω _m - The

predetermined measurement velocity ω _m is(are) selected so that windage effects imposed on the rotor cause the time(s) required by the rotor on the shaft to reach the measurement velocity differ by a measurable amount from the time required by each of the others of the plurality of rotors to reach the measurement velocity ω - A rotor identity signal based upon the windage of the rotor is generated in response the signal(s) representative of the time t _a . The rotor identity signal may be generated using a look-up table or a comparator for comparing the time signal ta to a time reference signal tret- The time reference signal t _re f may be derived from a first rotor identification system.

In accordance with a third embodiment of the present invention a selector for selectively applying either the velocity signal ω _a or the time signal ta to the rotor identity signal generator. The selector applies the selected signal in accordance with the relationship between the actual velocity of the rotor as measured at a predetermined time td with respect to a predetermined velocity cod-

Brief Description of the Drawings

The invention may be more fully understood from the following detailed description thereof, taken in connection with the accompanying drawings, which form a part of this application and in which:

Figure 1 is a highly stylized pictorial representation of a centrifuge instrument with which a control system in accordance with the present invention may be used, and includes a functional block diagram of the control system of the present invention; and

Figure 2 is a graph of the relationship between rotor speed and time for a hypothetical family of centrifuge rotors.

The Appendix (two pages) attached hereto following the description and preceding the claims is a C language source code listing of a program for implementing the present invention.

Detailed Description of the Invention

Throughout the following detailed description, similar reference numerals refer to similar elements in all Figures of the drawings.

Shown in Figure 1 is a stylized pictorial representation of a centrifuge instrument generally indicated by reference character 10 with which a rotor identification arrangement generally indicated by the reference character 50 embodying the teachings of the present invention may be used. The instrument 10 includes a framework schematically indicated at 12. The framework 12 supports a bowl 14. The interior of the bowl 14 defines a generally enclosed chamber 16 in which a rotating element, or rotor, 18 may be received. Access to the chamber 16 is afforded through a door 20. The bowl 14 may be provided with suitable evaporator coils (not shown) in the event that it is desired to refrigerate the bowl 14, the rotor 18 and its contents.

One or more energy containment members, or guard rings, 22 is(are) carried by the framework 12. The guard ring 22 is arranged concentrically with respect to the bowl 14 and serves to absorb the kinetic energy of the rotor 18 or fragments thereof should a catastrophic failure of the rotor 18 occur. The guard ring 22 is movably mounted within the framework 12, as schematically indicated by the rollers 24, to

permit free rotation of the ring 22 to absorb any rotational component of the energy of the rotor fragments. It is important to absorb the energy of the rotor and to contain the possible fragments which if permitted to exit the instrument may cause injury to an operator.

A motive source 30 is mounted within the framework 12. The motive source 30 may be any one of a well-known variety of sources, such a brushless DC electric motor, an induction motor, or an oil turbine drive. The motive source 30 is connected to or includes as an element thereof a drive shaft 34. The drive shaft 34 projects into the chamber 16. The upper end of the shaft 34 is provided with a mounting spud 36 which receives the rotor 18. Any one of a predetermined number of rotor elements may be received on the spud 36.

Whatever form is used the source 30 exhibits a predetermined output torque versus angular velocity profile. When asserted the source 30 is operative to accelerate a rotor 18 mounted on the shaft 34 to a predetermined operational angular velocity. A tachometer generally indicated by the reference character 38 is arranged to monitor the rotational speed (i. e., angular velocity) of the shaft 34 and thereby the rotational speed (i. e., angular velocity) of the rotor 18 received thereon. Any convenient form of tachometer arrangement may be utilized and remain within the contemplation of the present invention. An electrical signal representative of the actual angular velocity of the shaft 34 and of a rotor mounted thereon is carried by an output line 38L from the tachometer 38.

As mentioned earlier it is vitally important to accurately identify the rotor 18 mounted to the shaft 34. To this end the instrument 10 may include a first rotor identification system 42. The first rotor identification

system 42 includes a sensor 42S disposed within the chamber 16. The system 42 is operative to provide an identification signal on a line 44 representative of the identity of the particular rotor 18 mounted within the chamber 16. The ultrasonic rotor recognition system disclosed and claimed in copending application Serial Number 07/363,907, assigned to the assignee of the present invention, is preferred.

For reasons that will become more apparent herein the identification signal produced by the first rotor identification system 42 on the line 44 is utilized as an entry into a suitable reference table 46. Output lines 46V, 46T extend from the reference table 46. The signal on the line 46V represents the angular velocity ω _re f able to be achieved by the particular rotor as identified by the first rotor identification system 42 within a predetermined time following the initiation of a centrifugation run. Similarly, the signal on the line 46T represents the time t _re f required following the initiation of a centrifugation run for the particular rotor identified by the first rotor identification system 42 to achieve a predetermined angular velocity.

When a motive force is applied a body disposed in a nonevacuated or a partially evacuated environment, such as a rotor 18 mounted to the shaft 34 within the chamber 16, the body manifests two forms of resistance to motion in response to the application of the motive force. The first form of resistance is functionally related to the mass of the body and to its radially distribution. This form of resistance is termed inertia. Inertial resistance to acceleration is dominant in a nonevacuated or a partially evacuated environment at relatively low rotational speeds. The second form of resistance to motion is a fluid frictional effect functionally related to the configuration of the body. This effect is termed

windage. Windage is dominant in a nonevacuated or a partially evacuated environment at relatively high rotational speeds.

With reference to Figure 2 shown is a graphical depiction of the angular velocity ω versus time t for a family of four centrifuge rotors. Rotors 1 and 2 are regarded as low windage rotors, while rotors 3 and 4 may be viewed as high windage rotors. The windage of rotor 1 is less than that of rotor 2. Similarly the windage of rotor 3 is less than that of rotor 4. Every centrifuge rotor usable within a given centrifuge instrument exhibits a predetermined angular velocity versus time profile such as is indicated in Figure 2.

As may be viewed from inspection of Figure 2 a low windage rotor such as rotor 1 (or rotor 2) undergoes a relatively substantial increase in angular velocity Δωι_ for a relatively small time increment Δ_L. Conversely, a high windage rotor such as rotor 4 (or rotor 3) undergoes a relatively small increase in angular velocity Δωμ over a relatively substantial time increment ΔtH-

Thus, there may be defined a demarcating curve, shown in Figure 2 as a line Ld, which may be used to separate rotors that exhibit low windage effects from those that exhibit high windage effects. This circumstance is utilized in one aspect of the present invention, as will be described presently. For later reference there is defined a first predetermined decision point Pd along the curve of demarcation Ld- The point decision Pd is defined by the decision time td and the decision velocity ωd. The point Pd thus has the coordinates (td, ωd). A second predetermined point Pd 2 along the curve of demarcation l_d defined by the decision time t 2 and the decision velocity ωd 2 is also shown in Figure 2. The decision point Pd2 has the coordinates (td2, ωd2).

The rotor identification arrangement 50 in accordance with the present invention includes a timer 52 for providing a signal on a line 52L representative of elapsed time following initiation of a centrifugation run. Typically the timer 52 is initiated upon energization of the motive source 30.

In accordance with a first general embodiment of the present invention the rotor identification arrangement 50 includes means 54 responsive to the tachometer signal on the line 38L and to the timer signal on the line 52L for generating a signal on a line 54L representative of the actual measured angular velocity ω _a exhibited by a rotor 18 mounted on the shaft 34 at at least a first predetermined measurement time tm following initiation of rotation of the rotor 18. A signal representative of the measurement time t _m is applied to the means 54 on a line 58. The predetermined measurement time t _m is selected to correspond to a time when windage effects imposed on the rotor cause the velocity of the rotor on the shaft to differ by a measurable amount from the velocity of each of the others of the plurality of rotors. That is, the measurement time is selected at a point in the centrifugation run where windage effects will be significant and can be used to discern the identity of the rotor.

The signal on the line 54L representative of actual measured angular velocity ω _a at the measurement time tm is applied to means generally indicated by the reference character 60. The means 60 is responsive to the signal representative of the actual measured angular velocity ω _a for generating a rotor identity signal based upon the windage of the rotor 18. In accordance with the present invention the means 60 may take one of several forms.

In a first form the means 60 comprises a look-up table 62. Using the signal on the line 54L as an address the table 62

produces a signal on an output line 64 representative of the identity of the rotor 18 on the shaft 34. The identity signal on the line 64 may serve as the primary rotor identification signal. Alternatively, if the first rotor identification system 42 is provided, the signal on the line 64 may be used as a verification of the rotor identity derived by that means. For example, the identity signal on the line 64 may be compared with the identity signal on the line 44 to determine if an identification mismatch has occurred.

In a second aspect the means 60 may be implemented in the form of a comparator 66. The actual measured angular velocity ω _a on the line 54L is applied to one side of the comparator 66 while a reference angular velocity value ω _re f corresponding to a known rotor is applied to the comparator 66 over a line 68. The truth of the comparison determines the identity of the rotor 18 which is carried on an output line 70.

This arrangement is also believed to be useful as a verification of the first rotor recognition system 42. The reference angular velocity value ω _re f is derived from the reference table 46 responsive to the identity determined by the first rotor recognition system 42. A true comparison between the actual angular velocity ω _a and the reference angular velocity ω _re f verifies the identity determination made by the first rotor recognition system 42.

Alternatively the reference angular velocity value ω _r βf may be applied to the comparator 66 in accordance with a predetermined sequence, as by stepping through a table of angular velocity values corresponding to particular rotors stored in a suitable table 72.

In accordance with a second general embodiment the rotor identification arrangement 50 includes means 74 also

responsive to the tachometer signal on the line 38L and to the timer signal on the line 52L for generating a signal on a line 74L representative of the actual time t _a following initiation of rotation at which the rotor first reaches a predetermined measurement angular velocity ω _m - A signal representative of the measurement velocity ω _m is applied to the means 74 on a line 76. The predetermined measurement velocity ω _m is selected to correspond to a velocity when windage effects imposed on the rotor causes the time needed by each of the rotors able to be used on the shaft to differ by a measurable amount from the time required by the others of the plurality of rotors. That is, the measurement velocity is selected at a point in the centrifugation run where windage effects will be significant and can be used to discern the identity of the rotor.

The signal on the line 74L representative of actual measured elapsed time t _a needed to reach the measurement velocity ω _m is also applied to the means 60. In this instance the means 60 is responsive to the signal representative of the actual measured measured elapsed time t _a for generating a rotor identity signal based upon the windage of the rotor.

The signal on the line 74L may be used as an address to access an identity signal from the table 62. The resultant rotor identity signal is again presented on the line 64.

Alternatively the actual measured time t _a is applied to one side of the comparator 66 with a reference time value t _re f corresponding to a known rotor being again applied to the comparator 66 over the line 68. The identity of the rotor signal is again presented on the line 70 based on the truth of the comparison effected by the comparator 66. As before the reference time value tref may be derived from the reference table 46 responsive to the identity determined by the first rotor recognition system 42. The reference time value tref

may also be again applied to the comparator 66 in a predetermined sequence from the table 72.

Reverting to Figure 2 it may be observed that for a high windage rotor at the measurement time t _m , the actual measured velocities are relatively close to each other. This may be observed from the velocities ω _a -3 and ω _a -4 for the high windage rotor 3 and the high windage rotor 4, respectively. Conversely at the measurement time tm the respective actual measured velocities for low windage rotors are spaced relatively widely. The velocities ω _a -ι and ω _a -2 for the low windage rotor 1 and the low windage rotor 2 bear witness to this statement. It may similarly be observed that the actual times t _a- ι and t _a -2 required for the low windage rotor 1 and the low windage rotor 2, respectively, to attain the measurement velocity ω _m are relatively close to each other. However the actual times t _a- 3 and ta-4 required for the high windage rotor 3 and the high windage rotor 4, respectively, to reach the measurement velocity ωm are spaced more widely apart.

These observations make it clear that a low windage rotor may be more accurately identified using the means 54 which measures the actual speed ω _a of a rotor at a predetermined measurement time tm- Conversely, for a high windage rotor, more accurate identification may be made using the means 74 to measure the actual time fe required by the rotor to reach a predetermined measurement speed ω _m .

In accordance with yet another embodiment of the present invention a selector 78 responsive to both the tachometer signal on the line 38L and the timer output on the line 52L utilizes the coordinates (td, ωd) of a predetermined decision point Pd on the curve of demarcation l_d to determine whether an unknown rotor 18 on the shaft 34 lies in either the

high windage or the low windage regime. Based on the results of this determination either the means 54 or the means 74 is selected. If the actual velocity ω _a of the rotor on the shaft at the time td is greater than the velocity ωd the rotor lies in the low windage regime. In this event the output line 78L is asserted. Alternatively, if the actual velocity ω _a of the rotor on the shaft at the time td is less than the velocity ω the rotor lies in the high windage regime. This causes the line 78H to be asserted.

One convenient implementation uses the output on a line 78H (high windage) or 78L (low windage) to assert a switch 80H or 80L thereby to connect the output of either the means 54 or the means 74 to the means 60. In addition, if the means 60 is implemented using the comparator 66 the output of the selector 78 may be used to close a switch 82 which applies either the reference time value or a reference velocity value ω _re f from the table 46 to the line 68 to the comparator 66.

If only one decision point Pd on the curve of demarcation Ld is used it should be judiciously chosen so that a decision as to the regime in which the rotor falls (i. e., low windage or high windage) is made as early as practicable in the centrifugation run. This circumstance permits the identity determination to be made at a time when windage effects are significant yet before a potential safety hazard may develop. The decision point Pd should be selected to properly categorize a low inertia, high windage rotor, which may undergo an initial rapid acceleration due to its relatively low inertia before windage effects become dominant.

It should be understood that if the second decision point Pd2 on the curve of demarcation Ld is used the slope of the curve of demarcation between the decision points Pd and Pd 2

may serve as a useful indicator of the appropriate regime to which the rotor on the shaft belongs.

Since the present invention relies on windage effects produced by a given rotor to identify the same it may be appreciated that in a practical application the control arrangement 50 should include a calibration scheme to compensate for the effects of atmospheric pressure at the locality where the instrument is being used and to compensate for idiosyncrasies (as in drive torque versus velocity profile, for example) between centrifuge instruments.

To this end means generally indicated by the reference character 86, 88 are respectively connected into the output lines from the means 54 and 74 for scaling the signals on the respective lines 54L, 74L by a predetermined scaling factor. The scaling factor serves to adjust the value of the signal on the line in which it is connected to compensate for any locality and/or individual instrument effects.

In practice the calibration is done using a reference rotor of precisely known windage and having a precisely known velocity versus time profile in a standardized instrument at a standardized pressure (e. g., atmospheric pressure at sea level). The reference rotor is used in the instrument and the compensating means 86, 88 appropriately adjusted to bring the actual signal values on the lines 54L, 74L into predetermined close tolerance with the reference values known to be produced by the reference rotor under the standardized conditions.

-o-O-o-

As noted earlier each centrifuge rotor usable within a given centrifuge instrument exhibits a predetermined angular

velocity versus time profile. Figure 2 illustrates such hypothetical profiles for each of four rotors. Heretofore the description of the present invention made clear the manner in which a rotor may be identified on the basis of a single point along the profile. However, accuracy of identification may be enhanced if a plurality of points (i. e., two or more points) along the velocity versus time curve are used to identify an unknown rotor.

For example, reference to the velocity versus time curves for the rotor 1 or the rotor 3 as shown in Figure 2 is invited. In accordance with this further aspect of the present invention the means 54 may be asserted to generate on a line 54L a signal representative of the actual measured angular velocity ω _a exhibited by a rotor 18 mounted on the shaft 34 at at least a second predetermined measurement time t _m 2 following the first predetermined measurement time t _m . In this manner a velocity versus time profile of the unknown rotor may be constructed.

If, for example, the unknown rotor is in fact the rotor 1 the second signal on the line 54L is thus representative of actual measured angular velocity ω _a- r for that rotor at the measurement time t _m 2- From the profile generated using the information representative of actual measured angular velocities ω _a- ι and ω _a .r at the respective measurement times t _m and t _m 2 it is believed that a more accurate identity signal of the unknown rotor can be generated.

The actual angular velocities ω _a -ι and ω _a- r may be used in a variety of ways to generate the identity signal. For example, each velocity measurement signal ω _a- ι and ω _a -r may be used as an address to the table 62 and a consensus (or unanimity) of identity outputs from the table 62 may be required before an identity signal is presented on the line 64.

Alternatively, a point-by-point comparison may be made using the comparator 66. Each of the actual angular velocities ω _a -ι and ω _a -r is compared to a respective reference velocity ωref and ω _re f corresponding to each respective reference time tm and tm2- The reference velocities ω _f ef and ω _r ef are derived from the table 46 (responsive to the first identification system 42) or from the store 72.

As yet another alternative the set of actual angular velocities may be used to generate the slope of a velocity versus time curve of the unknown rotor. The slope of the curve may be compared to a reference slope (e.g., as derived from the first identification system) or to the slopes of a family of rotors to determine the rotor identity. If more than two actual velocities are measured an equation may be fit to the set of angular velocities. The coefficients of the terms of the equation may be compared to a reference set of coefficients (e. g., as derived from the first identification system or from the coefficients of the equations of a family of rotors stored in the store 72) to determine the rotor identity.

The measurement of each of the angular velocities may be made with reference to zero velocity. However, especially when dealing with a plurality of actual velocities comprising a velocity versus time profile, it is believed more advantageous to use the incremental difference between the angular velocity ω _a -ι and the angular velocity ω _a- r to identify the unknown rotor. The change in velocity over the time increment t _m to tm2 (i- e., the acceleration) is thus used to identify the rotor.

If the means 74 is used the signals on the line 74L represent the actual measured elapsed times ta and ta' needed for the rotor to reach the respective measurement velocities ω _m and ω _m 2- In the context of Figure 2, for the instance of the

rotor 3, the elapsed times t _a -3 and ta-3' and the respective measurement velocities ω _m and ω 2 are shown. As analagously described immediately above these time signals may be applied to the means 62 or to the comparator 66 (deriving its references from the store 72 or from the table 46). In addition to using the time values on a point by point basis the value of the difference (slope) between the times t _a -3 and t _a -3' may also be applied to the table 62 or to the comparator 66.

-o-o-o-

It should be understood that the invention as hereinabove described is preferably implemented using a microprocessor based programmable device operating in accordance with a suitable program.

As an example of a possible microprocessor implementation a listing in C source language of a program implementing the present invention is set forth in the Appendix hereto. The subroutine "primary id" performs the functions 54, 74, and the selector function 78. Based on the determination of whether the rotor lies in the low windage or the high windage regime the "low windage loop" or "the high windage loop" implement the functions of the comparator 66 and the store 72.

-o-0-o-

Although, in general, the invention is believed to find its greatest utility in a nonevacuated centrifuge instrument it should be understood that the present invention may also be used with advantage in a partially evacuated centrifuge instrument. Such an instrument is one that operates at a chamber pressure that is less than atmospheric but still

sufficiently high to exert windage effects on a rotor being spun therein.

It should be noted that the various embodiments of the present invention are all set forth on the diagram of Figure 1 for economy of illustration. A skilled artisan may find it convenient to select only portions of the arrangement as shown in Figure 1 and still remain within the contemplation of the present invention. For example, any discussed form of the means 60 as shown in Figure 1 may be used to implement the present invention. In addition the exact time or velocity values defining the points Pd, the velocity values ω _m and ω _m 2. or the time values t and t 2, could vary based on the torque output of the motive source. However any appropriate time or angular velocity values may be chosen so long the identity determination can be made when windage is dominant but before a safety hazard develops.

Moreover, those skilled in the art, having the benefits of the teachings of the present invention as hereinabove set forth, may effect numerous modifications thereto. Such modifications are to be construed as lying within the contemplation of the present invention, as defined by the appended claims.

/* Thi* progran in c demon*trateβ th* prianry rotor windage identification * /* system using the βelecto' .

/* Global Variable* are a* followsi */ int rotor; /* rotor id index used to indentify rotor */ int ta - 30; /* measurement tiae in second* at which to measure rotor speed */ int W_B ■ 128; /* measurement speed in rpm at which to measure tiae of run */ int t; /* tiae index */ int w; /* rotational velocity measured a tiae t */ int t_d 15; /* decision time to determine which windage loop */ int w~d 64; /* decision speed to determine which windage loop */ int hTgh_windage _# _tbl; /* table containing high windage times for each rotor */ int low_windag*_tbl; /* table containing low windage times for each rotor */

/* Low Windage Loop - This loop determines the rotor that most */ /* closely aatches the windage at the given minimum' tiae */ int ^■ > ow_windage_loop { ) { for ( rotor_id-θ; rotor_id< nuB_rotors_in_tablβ; rotor_id++)

_ if (low_windagβ_tbltrotor_id] v) return(rotό ^" r_id); ~

70

)

/* High indage Loop - This loop determines the rotor that most */ /* closely aatches the windage at the given minimum speed */ int high _e windage_loo ( J { for ( rotor_id«θ; rotor_id< num_rotorβ__in_tabl* ; rotor_id++)

{ if (high .indaga .bl [rotor id) t) return (rotor_id) ; *. ^"" l

)

I* Thiβ id system will return a rotor given the velocity and tiae */ int ρriaary_id() { int total run_time; /* total run tiae of a rotor tracked by t */ int low_wTndage_flag-0; /* this flag selects a low windage rotor */ int high_windage_ lag»0; /* this flag selects a high windage rotor */

/* Loop for the entire run and aonitor the time of the run */

7β

APPENDIX

while (t < total run time) { w « mβasure velocityO; /* get velocity at time t

/* At the decision time, decide whether we have a low windage or */ /* high windage rotor */

>

/* For low windage rotors, determine the rotor id at the measured time */ if ( ( t — t_m ) fifi ( low windage_flag) ) rotor ■ low windage_loop(); return(rotor) ;

* For high windage rotors, determine the rotor id at the measured speed * if (( v — w m) &f (high windage lag)) _u rotor - high windage_loop<)/ return(rotor] ^" ;

^}

•t++;

_. } /* end of run */

\ j* end of primary.J.d routine */ ainO

rotor ■ primary id() ; print (" The rotor was identified to be id φ %i . ", rotor) ;

}

APPENDIX