US10240988B2 - Fully automated calibration for acoustic pyrometry - Google Patents
Fully automated calibration for acoustic pyrometry Download PDFInfo
- Publication number
- US10240988B2 US10240988B2 US14/894,010 US201414894010A US10240988B2 US 10240988 B2 US10240988 B2 US 10240988B2 US 201414894010 A US201414894010 A US 201414894010A US 10240988 B2 US10240988 B2 US 10240988B2
- Authority
- US
- United States
- Prior art keywords
- measuring vessel
- acoustic pyrometry
- transceivers
- acoustic
- pyrometry measuring
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active, expires
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K15/00—Testing or calibrating of thermometers
- G01K15/005—Calibration
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K11/00—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
- G01K11/22—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using measurement of acoustic effects
- G01K11/24—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using measurement of acoustic effects of the velocity of propagation of sound
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K13/00—Thermometers specially adapted for specific purposes
- G01K13/02—Thermometers specially adapted for specific purposes for measuring temperature of moving fluids or granular materials capable of flow
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K13/00—Thermometers specially adapted for specific purposes
- G01K13/02—Thermometers specially adapted for specific purposes for measuring temperature of moving fluids or granular materials capable of flow
- G01K13/024—Thermometers specially adapted for specific purposes for measuring temperature of moving fluids or granular materials capable of flow of moving gases
-
- G01K2013/024—
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K2205/00—Application of thermometers in motors, e.g. of a vehicle
- G01K2205/04—Application of thermometers in motors, e.g. of a vehicle for measuring exhaust gas temperature
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/02881—Temperature
Definitions
- This disclosure is directed to methods and systems for acoustic pyrometry.
- Acoustic pyrometry is the process of using sound waves to measure temperature inside an object.
- a pyrometer measures the time-of-flight (TOF) of acoustic waves traveling between a fixed sound source (transmitter) and microphone (receiver) pair.
- TOF time-of-flight
- the speed of sound C traveling in the medium can be computed, and the average temperature T in Kelvin of the path traversed by the sound wave can be computed using
- tomographic reconstruction can be used to generate a temperature distribution map of the acoustic traveling plane.
- Exemplary embodiments of the invention as described herein are directed to a fully automatic calibration system that can accurately measure the locations or coordinates of end points of each sound path. Furthermore, a calibration algorithm and system according to an embodiment of the invention can be directly applied to other pyrometry applications when system calibration is necessary.
- a method of calibrating transceiver positions inside an acoustic pyrometry measuring vessel that contains a plurality of transceivers including acquiring time-of-flight (TOFs) ⁇ ti ,j measurements from a plurality of pairs i,j of transceivers inside the acoustic pyrometry measuring vessel, and using an estimated radius of the acoustic pyrometry measuring vessel to estimate errors ⁇ j of displacement angles of the transceivers.
- TOFs time-of-flight
- the method includes estimating a radius of the acoustic pyrometry measuring vessel from an average of the acquired TOF measurements.
- the method includes repeating the steps of estimating a radius of the acoustic pyrometry measuring vessel and using the estimated radius to estimate errors ⁇ j of displacement angles until convergence.
- the radius r i of each transceiver the acoustic pyrometry measuring vessel is estimated from
- the method includes determining a speed of sound in the acoustic pyrometry measuring vessel from a temperature and gas composition of a gas inside the acoustic pyrometry measuring vessel.
- the speed of sound C in the acoustic pyrometry measuring vessel is determined from
- T the temperature T in degrees Kelvin
- ⁇ a ratio between specific heats of the gas at constant pressure and constant volume
- R the gas constant
- m a molecular weight of the gas.
- the errors ⁇ j of displacement angles of the transceivers are estimated using a least squares fit.
- a system for calibrating transceiver positions for acoustic including an acoustic pyrometry measuring vessel, a plurality of transceivers disposed about a perimeter of the measuring vessel, the transceivers configured to transmit an acoustic wave signal upon receipt of an electric signal and to receive acoustic wave signals and convert a received acoustic wave signal into an electric signal, and a computer processor configured to transmit and receive the electric signals sent to and received from the transceivers, to determine time-of-flight (TOFs) ⁇ ti ,j measurements from each pair i,j of transceivers inside the acoustic pyrometry measuring vessel, to estimate a radius of the acoustic pyrometry measuring vessel from an average of the acquired TOF measurements, and to use the estimated radius of the acoustic pyrometry measuring vessel to estimate errors ⁇ j of displacement angles of the transceivers.
- TOFs time-of-flight
- the computer processor is configured to determine a speed of sound in the acoustic pyrometry measuring vessel from a temperature and gas composition of a gas inside the acoustic pyrometry measuring vessel.
- a non-transitory program storage device readable by a computer, tangibly embodying a program of instructions executed by the computer to perform the method steps for calibrating transceiver positions inside an acoustic pyrometry measuring vessel that contains a plurality of transceivers.
- FIG. 1 illustrates the distribution of eight transceivers inside a vessel of radius R, where ⁇ i depicts the location of transceiver i, according to an embodiment of the invention.
- FIG. 2 illustrates the temperature estimation convergence without autocalibration, according to an embodiment of the invention.
- FIG. 3 illustrates temperature estimation convergence with auto-calibration, according to an embodiment of the invention.
- FIG. 4 is a flow chart of a method for auto-calibration of the locations or coordinates of end points of each sound path, according to an embodiment of the invention.
- FIG. 5 is a block diagram of an exemplary computer system for implementing a method for auto-calibration of the locations or coordinates of end points of each sound path, according to an embodiment of the invention.
- Exemplary embodiments of the invention as described herein generally include systems and methods for measuring the locations or coordinates of end points of each sound path in an acoustic pyrometry system. Accordingly, while the invention is susceptible to various modifications and alternative forms, specific embodiments thereof are shown by way of example in the drawings and will herein be described in detail. It should be understood, however, that there is no intent to limit the invention to the particular forms disclosed, but on the contrary, the invention is to cover all modifications, equivalents, and alternatives falling within the spirit and scope of the invention.
- FIG. 1 is a schematic illustration of a plurality of transceivers that are distributed along the inside of an outer wall of a vessel 10 according to an embodiment of the invention.
- the transceivers are in signal communication with an external controller, such as a computer.
- An exemplary computer system for controlling the transceivers is depicted in FIG. 5 , described below.
- Let ⁇ t i,j denote the TOF measurement from transmitter i to receiver j.
- the TOF measurement and constant temperature information can be used to automatically calibrate the transceiver coordinates and exhaust geometry to improve the accuracy of the temperature estimation.
- all transceivers may be assumed to be the same distance r from the center of the exhaust, i.e., there is no displacement distance for the transceivers.
- An experiment to test and validate an automated calibration system can be designed as follows.
- a gas turbine is placed in a room with a temperature of 8.76° C., equivalent to 281.76° K.
- the TOF measurements among all different sound traversed paths are estimated using an approach disclosed in U.S. patent application Ser. No. 13/961,292, “Noise Robust Time of Flight Estimation for Acoustic Pyrometry”, filed on Aug. 7, 2013, assigned to the assignee of the present application, the contents of which are herein incorporated by reference in their entirety.
- FIG. 4 is a flow chart of a method for auto-calibration of the locations or coordinates of end points of each sound path, according to an embodiment of the invention.
- a container or vessel containing a gas of known composition is placed inside a room of known temperature and allowed to equilibrate.
- a method begins at step 40 by determining the speed of sound C 0 in the vessel from the temperature.
- the container includes a plurality of transceivers on the perimeter in positions that are approximately known.
- a electronic controller sends electronic signals to the transceivers to initiate transmission of acoustic signals, and receives an electronic signal from each transceiver when the transceiver receives an acoustic signal.
- An exemplary, non-limiting electronic controller is a digital computer, and the electronic signals may be digital electronic signals.
- the digital computer may be configured to perform the following computations.
- TOFs ⁇ ti ,j are acquired from all pairs i,j of transceivers.
- an acoustic wave is sent from one transceiver at a time, and signals are recorded from each of the remaining transceivers upon receipt of the acoustic wave, before sending an acoustic wave from another transceiver.
- the TOFs ⁇ ti ,j can be determined from the time a transmit signal was sent to the transmitting transceiver, and the time a signal was received from each receiving transceiver indicating receipt of an acoustic signal.
- the radius r i of each transceiver can be estimated by an average over all TOF measurements:
- the estimated values of all r i can be used at step 43 to estimate errors ⁇ j of the displacement angles of the transceivers from EQ. (3), above.
- An exemplary, non-limiting method of estimating the errors ⁇ j from EQ. (3) is a linear least squares.
- Step 42 can then be repeated, adjusting the displacement angles ⁇ j with the corrections calculated in step 43 , and step 43 can be repeated with the new radius value r i of transceiver i
- steps can be repeated from step 44 until convergence.
- FIG. 2 shows estimated temperatures 21 over time without auto-calibration, i.e. using the 9 input parameters describing location and geometry, without further verifying with the TOF measurements.
- embodiments of the present invention can be implemented in various forms of hardware, software, firmware, special purpose processes, or a combination thereof.
- the present invention can be implemented in software as an application program tangible embodied on a computer readable program storage device.
- the application program can be uploaded to, and executed by, a machine comprising any suitable architecture.
- FIG. 5 is a block diagram of an exemplary computer system for controlling the transceivers depicted in FIG. 1 and implementing a method for auto-calibration of the locations or coordinates of end points of each sound path, according to an embodiment of the invention.
- a computer system 51 for implementing the present invention can comprise, inter alia, a central processing unit (CPU) 52 , a memory 53 and an input/output (I/O) interface 54 .
- the computer system 51 is generally coupled through the I/O interface 54 to a display 55 and various input devices 56 such as a mouse and a keyboard.
- the support circuits can include circuits such as cache, power supplies, clock circuits, and a communication bus.
- the memory 53 can include random access memory (RAM), read only memory (ROM), disk drive, tape drive, etc., or a combinations thereof.
- RAM random access memory
- ROM read only memory
- the present invention can be implemented as a routine 57 that is stored in memory 53 and executed by the CPU 52 to process the signal from the signal source 58 .
- the computer system 51 is a general purpose computer system that becomes a specific purpose computer system when executing the routine 57 of the present invention.
- the computer system 51 also includes an operating system and micro instruction code.
- the various processes and functions described herein can either be part of the micro instruction code or part of the application program (or combination thereof) which is executed via the operating system.
- various other peripheral devices can be connected to the computer platform such as an additional data storage device and a printing device.
Abstract
Description
where γ is the ratio between the specific heats of the gas at constant pressure and constant volume, R is the gas constant and m is the molecular weight of the gas. In principle, if a sufficient number of TOF measurements are available in conjunction with knowledge of the paths traversed by sound waves, tomographic reconstruction can be used to generate a temperature distribution map of the acoustic traveling plane.
where C0 is a speed of sound, N is a number of transceivers, and θi is a displacement angle of transceiver i.
where T is the temperature T in degrees Kelvin, γ is a ratio between specific heats of the gas at constant pressure and constant volume, R is the gas constant, and m is a molecular weight of the gas.
where C0 is a speed of sound, r is the estimated radius, and {circumflex over (θ)}i is an intended angle of transceiver i.
where i must be greater than j for θi>θj and value on the left side of EQ. (2) to be non-negative, thus EQ. (2) is valid for all i and j. To take into account the uncertainty or allowed error margin during installation, each transceiver location can be modeled as
θi={circumflex over (θ)}i+Δθi ,i=1, . . . ,8,
where {circumflex over (θ)}i denotes an expected or intended location for θi, for example, {circumflex over (θ)}i=(i−1)×45°, i=1, . . . 8, i.e., all transceivers are uniformly distributed along a circle. After keeping linear terms of a Taylor expansion of EQ. (2) at ({circumflex over (θ)}i−{circumflex over (θ)}j)/2 and re-arranging the terms, the first order approximation becomes:
which is a linear function of Δθi for i=1, . . . , 8. Hence with EQ. (3), one can easily estimate Δθi using least squares to fit, therefore θi, i=1, . . . , 8, for a given radius r.
where the brackets < > indicate an arithmetic mean over transceiver i. The estimated values of all ri can be used at
Claims (13)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US14/894,010 US10240988B2 (en) | 2013-05-30 | 2014-05-29 | Fully automated calibration for acoustic pyrometry |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201361828936P | 2013-05-30 | 2013-05-30 | |
US14/894,010 US10240988B2 (en) | 2013-05-30 | 2014-05-29 | Fully automated calibration for acoustic pyrometry |
PCT/US2014/039971 WO2014194056A1 (en) | 2013-05-30 | 2014-05-29 | Fully automated calibration for acoustic pyrometry |
Publications (2)
Publication Number | Publication Date |
---|---|
US20160109304A1 US20160109304A1 (en) | 2016-04-21 |
US10240988B2 true US10240988B2 (en) | 2019-03-26 |
Family
ID=51033529
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US14/894,010 Active 2035-11-16 US10240988B2 (en) | 2013-05-30 | 2014-05-29 | Fully automated calibration for acoustic pyrometry |
Country Status (3)
Country | Link |
---|---|
US (1) | US10240988B2 (en) |
DE (1) | DE112014002598T5 (en) |
WO (1) | WO2014194056A1 (en) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2738372B1 (en) * | 2012-11-29 | 2018-02-28 | Ansaldo Energia Switzerland AG | Gas turbine temperature measurement |
US9581087B1 (en) | 2015-09-17 | 2017-02-28 | Siemens Energy, Inc. | Method and system for thermal expansion compensation in heated flow characterization |
CN105953947B (en) * | 2016-07-07 | 2018-06-29 | 青岛华迪科技有限公司 | The calibration equipment of boiler sound wave temperature measuring equipment |
US11788904B2 (en) * | 2020-12-04 | 2023-10-17 | Perceptive Sensor Technologies, Inc. | Acoustic temperature measurement in layered environments |
US11946905B2 (en) | 2020-12-30 | 2024-04-02 | Perceptive Sensor Technologies, Inc. | Evaluation of fluid quality with signals |
US11860014B2 (en) | 2022-02-11 | 2024-01-02 | Perceptive Sensor Technologies, Inc. | Acoustic signal detection of material composition in static and dynamic conditions |
US11940420B2 (en) | 2022-07-19 | 2024-03-26 | Perceptive Sensor Technologies, Inc. | Acoustic signal material identification with nanotube couplant |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6834992B2 (en) * | 1997-07-05 | 2004-12-28 | Combustion Specialists, Inc. | Acoustic pyrometer |
EP1785701A1 (en) | 2005-11-11 | 2007-05-16 | Siemens Milltronics Process Instruments Inc. | Apparatus and method for determining a temperature of a volume of gas |
US20120150413A1 (en) * | 2010-12-14 | 2012-06-14 | Bunce Richard H | Gas turbine engine control using acoustic pyrometry |
US8256953B2 (en) * | 2008-10-31 | 2012-09-04 | Yuhas Donald E | Methods and apparatus for measuring temperature and heat flux in a material using ultrasound |
US20130306375A1 (en) * | 2008-06-27 | 2013-11-21 | Wajid Rasheed | Drilling tool, apparatus and method for underreaming and simultaneously monitoring and controlling wellbore diameter |
US9453784B2 (en) * | 2013-09-04 | 2016-09-27 | Siemens Energy, Inc. | Non-intrusive measurement of hot gas temperature in a gas turbine engine |
US9702768B2 (en) * | 2012-09-05 | 2017-07-11 | Siemens Corporation | Noise robust time of flight estimation for acoustic pyrometry |
US9945737B2 (en) * | 2014-03-13 | 2018-04-17 | Siemens Energy, Inc. | Method for determining waveguide temperature for acoustic transceiver used in a gas turbine engine |
-
2014
- 2014-05-29 DE DE112014002598.4T patent/DE112014002598T5/en active Pending
- 2014-05-29 WO PCT/US2014/039971 patent/WO2014194056A1/en active Application Filing
- 2014-05-29 US US14/894,010 patent/US10240988B2/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6834992B2 (en) * | 1997-07-05 | 2004-12-28 | Combustion Specialists, Inc. | Acoustic pyrometer |
EP1785701A1 (en) | 2005-11-11 | 2007-05-16 | Siemens Milltronics Process Instruments Inc. | Apparatus and method for determining a temperature of a volume of gas |
US20130306375A1 (en) * | 2008-06-27 | 2013-11-21 | Wajid Rasheed | Drilling tool, apparatus and method for underreaming and simultaneously monitoring and controlling wellbore diameter |
US8256953B2 (en) * | 2008-10-31 | 2012-09-04 | Yuhas Donald E | Methods and apparatus for measuring temperature and heat flux in a material using ultrasound |
US20120150413A1 (en) * | 2010-12-14 | 2012-06-14 | Bunce Richard H | Gas turbine engine control using acoustic pyrometry |
US8565999B2 (en) * | 2010-12-14 | 2013-10-22 | Siemens Energy, Inc. | Gas turbine engine control using acoustic pyrometry |
US9702768B2 (en) * | 2012-09-05 | 2017-07-11 | Siemens Corporation | Noise robust time of flight estimation for acoustic pyrometry |
US9453784B2 (en) * | 2013-09-04 | 2016-09-27 | Siemens Energy, Inc. | Non-intrusive measurement of hot gas temperature in a gas turbine engine |
US9945737B2 (en) * | 2014-03-13 | 2018-04-17 | Siemens Energy, Inc. | Method for determining waveguide temperature for acoustic transceiver used in a gas turbine engine |
Non-Patent Citations (1)
Title |
---|
PCT International Search Report dated Nov. 6, 2014 corresponding to PCT International Application No. PCT/US2014/039971 filed May 29, 2014 (13 pages). |
Also Published As
Publication number | Publication date |
---|---|
DE112014002598T5 (en) | 2016-04-14 |
US20160109304A1 (en) | 2016-04-21 |
WO2014194056A1 (en) | 2014-12-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10240988B2 (en) | Fully automated calibration for acoustic pyrometry | |
CN111435159B (en) | Method and equipment for positioning | |
JP6525325B2 (en) | Method and device for determining device location | |
JP2009198454A (en) | Position detection system, position detection server, and terminal | |
BR102015001678B1 (en) | System and test method for a radar unit | |
TWI544822B (en) | Signal strength distribution establishing method and wireless positioning system | |
CN104076348A (en) | Radar beyond visual range base line passive cooperative localization method | |
EP3290944B1 (en) | Hybrid tdoa closed form hyperbolic and spherical iteration geo-location technique | |
US20210003436A1 (en) | Time-of-flight generating circuit and chip, flow meter and method of the same | |
JP2021508367A (en) | Phase comparison of multi-frequency transmission to assist in determining position or time | |
US20140327580A1 (en) | Using Measured Angular Coordinates of an Object Relative to a Directional Transceiver | |
KR101930137B1 (en) | Estimation Apparatus for Sound source and Estimation Method for Sound source Using the same | |
JP6659054B2 (en) | Radio wave monitoring device, radio wave monitoring system, radio wave monitoring method and program | |
US20170026930A1 (en) | Localization apparatus and localization method | |
Yu et al. | Indoor localization algorithm for TDOA measurement in NLOS environments | |
KR20090040083A (en) | Method and server for estimating position of terminal | |
JP2007263614A (en) | Device and method for calibrating echo sounder receiver position | |
US10557915B2 (en) | Providing an indication of a heading of a mobile device | |
CN109798917B (en) | Positioning precision inspection method | |
KR101241926B1 (en) | Conformity evaluation method of radome for 3 dimensional interferometric array by phase difference error estimation | |
JP6920596B2 (en) | IFB correction value estimation method, device and server | |
JP2008267973A (en) | Positioning system, mobile communication terminal, positioning method, positioning server, positioning ic chip and positioning program | |
US11275171B2 (en) | Location using round trip time vectors using an airborne platform | |
US20230358851A1 (en) | System and method for determining the orientation in space of a linear array of radio receviers | |
EP4006569B1 (en) | Two-way signal positioning method and two-way signal positioning system thereof |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: SIEMENS ENERGY, INC., FLORIDA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:DESILVA, UPUL P.;ULERICH, NANCY H.;REEL/FRAME:037245/0555 Effective date: 20151124 |
|
AS | Assignment |
Owner name: SIEMENS CORPORATION, NEW JERSEY Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:CLAUSSEN, HEIKO;FANG, TONG;ROSCA, JUSTINIAN;AND OTHERS;SIGNING DATES FROM 20151124 TO 20151207;REEL/FRAME:037255/0290 |
|
AS | Assignment |
Owner name: SIEMENS CORPORATION, NEW JERSEY Free format text: CORRECTIVE ASSIGNMENT TO CORRECT THE FORTH ASSIGNEE'S NAME PREVIOUSLY RECORDED AT REEL: 037255 FRAME: 0290. ASSIGNOR(S) HEREBY CONFIRMS THE ASSIGNMENT;ASSIGNORS:CLAUSSEN, HEIKO;FANG, TONG;ROSCA, JUSTINIAN;AND OTHERS;SIGNING DATES FROM 20151124 TO 20151207;REEL/FRAME:048089/0041 |
|
AS | Assignment |
Owner name: UNITED STATES DEPARTMENT OF ENERGY, DISTRICT OF CO Free format text: CONFIRMATORY LICENSE;ASSIGNOR:SIEMENS ENERGY, INC.;REEL/FRAME:048188/0156 Effective date: 20181106 |
|
AS | Assignment |
Owner name: SIEMENS ENERGY, INC., FLORIDA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:SIEMENS CORPORATION;REEL/FRAME:048264/0143 Effective date: 20190121 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 4TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1551); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY Year of fee payment: 4 |