: 12~050'7 1 50,~72 METHOD AND APPARATUS FOR MEASURING TEMPERATURE
PROFILE WITH A SINGLE OPTICAL FIBER
BACKGROUND OF THE INVENTION
The inv~ntion is generally related to the use of fiber optic technology for thermal radiation detection.
More particularly, the invention provides a method and an apparatus for determining a temperature profile along a single optical fiber.
It is becoming widely recognized that various types of sensor devices can be advantageously made by using optical fibers either as the medium for data trans-mission, or as the sensor transducer, or both. Fiberoptics for sensor systems may be useful in the presence of a high electromagnetic noise background or in environments where electrical signals cannot be used, such as in the presence of explosive atmospheres. Fiber optic sensor systems can be implemented in the measurement of a variety of parameters, such as current, pressure, moisture and temperature.
Temperature may be measured with optical fihers in ~everal Xnown ways. For example, optical fibers have been developed in which the light propagation character~
istics of the fiber are dependent on temperature. In U.S.
Patent 4,151,747, an optical fiber cooperates with a light ~ource and a detector for sensing changes in the tempera-ture being monitored. The amount of light which passes through the fiber varies with changes in the temperature ` ; 1;i~1~50~7
2 50,972 of the fiber. The disadvantage of this sensor is that the optical fiber is very difficult to manufacture because extremely accurate control of the optical and thermal properties of the fiber core and cladding is required.
Another type of optical fiber temperature sensor measures the internally generated, black body radiation emitted from the fiber when it is heated. This sensor is taught in "Fiber-optic Temperature Sensor Based on Intèrnally Generated Thermal Radiation", M. Gottlieb and G. B.
~randt, APplied OPtics~ Vol. 20, No. 19, October 1981.
Both of the temperature sensors described therein only provide an indication of temperature along the hottest region of the fiber and thus function only as a "hot-spot"
probe. Although such fiber optic "hot spot" probes are of potential value in a variety of applications, the exis-tence of multiple hot spots along the fiber can, in some cases, significantly degrade the accuracy of the sensor.
It is an object of this invention to provide a method and an apparatus for determining a temperature profile along a single optical fiber by evaluating the spec~ral power density of the internally generated thermal radiation of an optical fiber.
Infrared radiation is emitted by all solid objects which are not perfectly transparent. This emis-sion is often referred to as "black body radiation".According to the Stefan-Boltzmann law, the total rate per unit area of emission of energy of all wavelengths is directly proportional to the fourth power of the absolute temperature. And, according to Wien's law, the wavelength of maximum intensity is inversely proportional to the absolute temperature of the emitting body. All solid objects have an emissivity between zero and one; a per-fectly transparent object which absorbs nothing and emits nothing, has an emissivity of zero, while an ideal, per-fectly black object has an emissivity of one. F~ber opticmaterial is a semi-transparent material with an emissivity value between zero and one~
~2~0soi
3 50,972 The concept of monitoring temperature with an optical fiber by measuring the internally generated radia-tion of the fiber itself is dislosed in the article "Fiber-optic Temperature Sensor Based on Internally Gener-ated Thermal Radiation" by M. Gottlieb and G. B. Brandt;
Applied Optics; Vol. 20, No. 19; 1 October 1981. In this article, a method of determining the hottest spot along a fiber is described. This method assumes that only one hot spot exists along the fiber and that this single hot spot generates nearly all the thermal radiation in the fiber. This radiation propagates through the fiber and is measured. A single number, the output voltage of the detector at the end of the fiber is converted into a temperature.
There are several assumptions made in the above-identified article. One hot spot exists along the fiber which has a constant temperature 'T' across its length 'Q'. The hot spot is a distance 'L' away from the detect-or. The length of the fiber, not including the 'hot spo~' is cold and contributes no thermal radiation. The dis-tance 'L' of cold fiber, attenuates the radiation emanat-ing from the hot spot by a total amount of e dL, where is the absorptivity of the fiber and e is a constant (e=2.732), the base of the natural log. It is also as-sumed that the fiber is not perfectly clear and that someradiation passing through the fiber is lost because the fiber absorbs it. The clearness of the fiber, its absorp-tivity, is indicated by ~ . The lower the ~ the clearer or more perfect the fiber with ~-0 being a perfect fiber.
Accordingly, in a perfect fiber, if one watt of light enters the fiber (Pin) and travels the length of the fiber, one watt will come out of the fiber (PoUt): PoUt =
p in e-dL. Thus, if ~=: Pout Pin However, fibers are not perfectly clear and~>0. It is also assumed by the Gottlieb et al. disclosure that the clearness, or the ability of the fiber to pass light, ~, ` ~2~0S07
4 50,972 is not dependent upon the color of light selected to measure fiber clearness. Accordingly, the absorptivity of the fiber is assumed to be a constant. Every part of the fiber is exactly the same as any other part, i.e., it has the same diameter, absorption constant, etc. It is also assumed that the fiber is in good, uniform, thermal con-tact with the object being monitored. Finally, the de-tector at the end of the fiber has a flat response and measures integrated power between two selected wave-lengths.
The single hot spot sensor of Gottlieb et al.only provides an estimate of the hottest spot along a fiber's entire length. If in fact, two or more hot spots of similar temperature exist along the fiber, only one estimate of the hottest spot along the fiber is obtained, and in this situation, that value is obtained through a weighted sum of both hot spots. Accordingly, the weighted sum could yield a totally erroneous result.
In order to determine the temperature of an object at every point along the length of the optical fiber in contact therewith, a temperature distribution monitor was developed and comprises a multifiber bundle which monitors an entire temperature distribution. Each fiber of the bundle has a predetermined absorption constant ~ which is measurably distinct from that of every other fiber. This disclosure assumes that 'N' hot spots exist along the fiber which has a temperature Tl, T2.,,,TN across its length and each hot spot temperature is a constant value across its length.
It is further assumed that all parts of each fiber gener-ate radiation which is transmitted to the detector as a ` ; lXlOS07 50,972 weighted sum. Thus for a bundle containing N fibers, N
values of integrated output power are obtained. Because each fiber in the bundle has a distinct absorptive con-stant ~, the attenuation and generating characteristics of each fiber are different, thus N values of integrated power are measured by N detectors. It is also assumed that absorbtivity ~ is constant with wavelength and that the total length and diameter of each fiber are homogene-ous .
This multifiber bundle optical device provides an estimation of the actual thermal distribution, the accuracy of which is reflected in the number of fibers in the bundle. However, the device does utilize an inte-grated power output which is applied to a linear system of 15 equations.
SUMMARY OF THE INVENTION
A method and apparatus for analyzing the thermal radiation output of an optical fiber capable of the self-generation of a thermal radiation output in response to 20 the temperature of the environment proximate the optical fiber. The optical fiber is characterized by a predeter-mined discrete absorption value for each wavelength. The 3 method includes the steps of separating the spectral power i density of the thermal radiation output into a deter-25 minable number of bandwidths and measuring the output spectral power of each bandwidth to provide an estimation of the temperature profile along at least a portion of the optical fiber.
The invention also provides a temperature moni-30 toring apparatus which utilizes a single optical fiber and means for measuring the spectral distribution of the self generated thermal radiation. The apparatus also includes means for selectively identifying a plurality of predeter-mined wavelengths which are analyzed for spectral output 35 power content.
~loSaq 6 50,972 BRIEF DESCRIPTION OF THE DRAWINGS
The above as well as other features and advant-ages of this invention will become apparent through con-sideration of the detailed description in connection with the accompanying drawing which schematically illustrates a single optical fiber temperature profile sensor system, all in accordance with the teachings of this invention.
DETAILED DESCRIPTION OF THE INVENTION
-This invention provides a method and an appar-atus for determining a temperature profile along a single optical fiber. Temperature profiling with a single fiber is possible by assuming that the absorption value in a single fiber is different for a given number of discrete wavelengths and therefore, a given number of equations of spectral power distribution can be derived. As will be more fully described below, the spectral power distribu-tion measurement is analyzed by a given number of non-linear equations.
Considering the schematic illustration of the apparatus o~ this invention, a single optical .iber tem-perature sensor is generally indicated by the reference character 11. The system 11 includes a single optical fiber 13 which is capable of the self-generating of thermal radiation in response to the temperature of the environ-ment or thermal source 15 proximate thereto. The opticalfiber 13 is characterized by having a discrete absorption value for each wavelength of a spectral distribution of the self-generated thermal radiation. A means for measur-ing the spectral distribution of the thermal radiation output is operably associated with the optical fiber 13 and includes a monochromator or similar device which isolates narrow portions of the spectrum through the dispersion of light. The accompanying illustration depicts a lens 17 which focuses the spectral output of the fiber 13 onto a diffraction grating schematically indicated at 19. A detector array means 21 measures the spectral output power of N bandwidth as diffracted by the grating 19 or, for example, isolated by ~ - ~2~0S(~7 7 50,972 a monochromator. Alternatively, an acousto-optic tunable filter can be employed to select predetermined bandwidths in the spectral output of the optical fiber 13 for analy-sis. The analysis of the spectral power output requires an evaluation through a somewhat cumbersome series of simultaneous equations. Accordin~ly, the evaluation can be affected in the spectral distribution measuring means by connecting the detector array to a c~mputer system 23 which can include a multiplexer 25, a high gain amplifier 27, an analog to digital converter 29 and a computer 31.
The spectral power analysis which consists of n number of thermal radiation measurements, can be sent to the com-puter from the detector array and mathematically asso-ciated with the absorption constant, a for the specific bandwidth being analyzed. Following the mathematical solution of the n number of simultaneous equations, as will be outlined below, the results can be automatically displayed on a readout device 33.
The use of a single optical fiber to measure a temperature distribution is based on the assumption that the value of the absorption constant in one fiber is - different at n wavelengths and that therefore N equations of spectral power distribution can be derived. However, as indicated above, these equations are nonlinear. More-over, because spectral power distribution is being measured, a diffraction grating or similarly functioning device is used to define the radiation output into N
wavelength bands. The strength of each of the N wave-length bands is then measured to provide an estimate of ; 30 the spectral power density.
Certain assumptions on the behavioral character-istics of an optical fiber are relied upon in the method and apparatus of this invention. There exists along the length of the fiber N hotspots, each of which has a tem-perature Tl, T2...,TN. All parts of a single fiber gener-ate radiation which is transmitted to a diffraction grat-ing at the end of the fiber and which is broken into N
1;210S07 8 50,972 wavelength bands. The strength of each wavelength band is measured with N detectors and N values if spectral power is obtained. The absorption constant of a single optical fiber has N values at N different wavelengths. The atten- -uation and spectral output power generation character-istics of the fiber at N wavelengths is different, thus N
~alues of spectral power distribution are measured. The absorptivity of the fiber varies with wavelengths but does not vary throughout the fiber. In other words, the char-acteristics of the fiber described above remain constantalong the entire length of the fiber.
In order to fully appreciate the derivation of the equations by which the spectral power output of a single optical fiber sensor of this invention is measured, it is necessary to build on some of the basic concepts described in the background portion of this specification.
The per unit length ability of a material to absorb external thermal radiation is given by its absorp~
tivity ~. To measure a, which herein is examined as a function of wavelength, ~. If a light of wavelength ~1 provides a known input power, Pin, into a fiber of known length, L, the power output, PoUt, is defined by the equation:
P t = Pi e ~(~)L (1) 8ased on this equation, the absorptivity a of a fiber at N
wavelength ~ can be calculated.
~ n analysis of a temperature sensing fiber can be made on the basis of the well-known theory of blackbody radiation and the light guiding properties of optical fibers. The total power density radiated from the surface of a heated body, between the wavelengths ~S and ~L~ is:
f L Cld~
~S E (~T) ~5leC2/~T-1] (2) where E is the emissivity or absorptance of the solid 9 50,972 object and is in general a function of temperature and wavelength and P is the total integrated power emanating from the surface. Building on equation 2 above, the emissivity of a slab having a thickness dx on the surface of the object is given by the statement:
E - absorptance = dx (3) Combining equations 2 and 3 yields:
f ~L C1d~
P = J~ a dx ~51ec2/~T-l]
where T is the temperature of the slab and ~ is absorptiv-ity. Now, radiation is attenuated as it travels through ! the semitransparent object to the surface thereof by:
-ax (5) where x is the distance to the surfa~~e and a i s the ab-sorptivity. Summing the radiation emanating from within the semitransparent object renders an output radiation according to this statement:
L
x=0 ~5 Se-~x dx ~S~eC2/~T(X)-l] (6) where T(x) is the temperature of the object as a function of x. In an optical fiber, only free radiation along the ¦ 20 axis of the fiber is going to reach the detector at one end thereof. Thus, a geometric constant must be placed in front of this equation. For an optical fiber with dia-~eter D, cladding having an index of refraction nClad, and a core having an index of refraction nCOre~ the geometric constant is:
a = ~4 ( 1 ~ n ) (7) -` ~210S07 50,972 The final expression for integrated power emanating from a fiber of absorption a(~), with a temperature distribution T(x) across its length is:
~b x=L
~b r ( Cla(~)e a~)Xdxd~
~ ~ = a J x=O ~51eC2/~T(X)-1] (8) where a is equal to equation 7, D is diameter of the fiber in centimeters, Cl is 3.74xlO 12 (watts/cm2~; C2 is 1.43 (cm-K), a(~) is absorption as a function of wavelength (cm 1), T(x) is temperature as a function of distance along the fiber (cm), L is the total length of the fiber, and ~a and ~b are detector cutoff wavelengths (~m).
Hot spot sensors described above, merely utilize a detector with cutoff wavelengths ~a and ~b to evaluate equation (8). Then by assuming a temperature profile consisting of one hot spot with a typical width (deter-mined by the application) an accurate estimate of the ! ~ hottest point along the profile can be determined. Incontrast, the present invention examines the spectral content of the output radiation. Using the fundamental theorm of calculus aP/a~ can be shown to be:
a~ ~x=o Cla(~) (e) ( ) dx (9) The temperature distribution T(x) is trea~ed as a step-like function along the ~iber where each step is desig-nated by Tl, T2,...TN. Additionally, a(~) is assumed to be known at N points in wavelength and is designated by 1' a2,...aN. P(~) is then measured at each of the wave-lengths yielding:
aP aP aP
a~l; a~2; a~N ' This formulation yields N non-linear equations in N un-knowns as follows:
` lZ~OS~7 11 50,972 (10) aP a clol ( ~ e 1 d3 eC2/~Tl ~J e-llx eC2/~lT2~1 (Xn-l e ~lXd~ ec~
an ~ ec2/~nTl~ e n (xn--L e-~nxc~) c2/~Tn-l~
h I b e~anX dx = e~aa _ e-ab These equations can be solved by using various numerical techniques, for example, Newtons Method or the method of successive iterations. They would be most appropriately used knowing P(~).
Using these calculations, an estimate of the actual p~wer distribution can be obtained, with the qual-ity o~ the estimate being enhanced according to the number of wavelength bands evaluated. This system has the advan-tage of being very flexible. The more wavelengths of spectral power density measured, the greater the accuracy to which the actual temperature distribution is measured.
12~05~7 12 50,972 The temperature range within which an optical fiber temperature sensor according to this invention can operate is a function of the heat resistant characteris-tics of the optical fiber's cladding and the impurities within the optical fiber itself which generate the thermal radiation spectral power output.
What has been described is a method and appar-atus for measuring a temperature profile with a sensor utilizing a single optical fiber by examining the spectral power output of the optical fiber.