Data processing: measuring – calibrating – or testing – Measurement system – Temperature measuring system
Reexamination Certificate
2000-10-18
2003-05-13
Bui, Bryan (Department: 2863)
Data processing: measuring, calibrating, or testing
Measurement system
Temperature measuring system
C702S094000, C702S196000
Reexamination Certificate
active
06564164
ABSTRACT:
BACKGROUND OF THE INVENTION
1. Field of the Invention
The present invention relates in general to a method for optimally determining sensor positions for an acoustic pyrometry, and more particularly to such a method which is capable of obtaining a minimum number of sensors satisfying a given error limit on the basis of an effective independence (EfI) method that is occasionally used for an acoustic holography or Indirect input measurement and determining positions of the obtained number of sensors suitable for error minimization.
2. Description of the Prior Art
As well known in the art, an acoustic pyrometry is a method useful to real-time measurement of temperature distributions on a two-dimensional sectional area. This pyrometry can be particularly useful in measuring a spatial temperature distribution and its time-based variations in a high-temperature place where a boiler is located or a harmful environment where a chemical reactor is located. For example, information regarding temperature distributions in a boiler used in a power plant is required in determining whether the boiler is abnormally operated and normally operating the boiler.
For measurement of temperature distributions in a boiler using the acoustic pyrometry, a plurality of acoustic emitters (also used as sensors) are mounted on the sidewall of the boiler and each of them senses sounds emitted from the others. Then, a sound propagation time based on the combination of the sensors is measured and subjected to an appropriate mathematical calculation process, thereby measuring a two-dimensional temperature distribution within the boiler, not a temperature at one point. At this time, an error of the measured temperature distribution is greatly influenced by the number and positions of the used sensors.
Now, a detailed description will be given of the basic principle of the above-stated acoustic pyrometry.
The acoustic pyrometry is a method for measuring a temperature in space using the fact that a propagation velocity of a sound is a function of the temperature. Namely, assuming that given conditions are a normal state, isentropic process, ideal gas, etc., it is well known from a continuity equation, momentum equation and state equation that the relation of the following equation 1 exists between the propagation velocity of the sound and the temperature.
C={square root over (&ggr;RT)}=K{square root over (T)}
[Equation 1]
where, C=acoustic velocity, T=temperature of medium [K],
K
=
γ
R
=
c
p
R
0
c
v
M
,
&ggr;=specific heat ratio, R
0
=gas constant and M=molecular mass.
Using the expression of the acoustic velocity in the above equation 1, it can be seen that the following equation 2 is established between a travel time t and travel distance dl of a sound pulse propagated within a sound field.
t
=
ⅆ
l
C
(
x
,
y
)
=
ⅆ
l
K
T
(
x
,
y
)
[
Equation
2
]
FIG. 1
is a graph illustrating the relation between sensor positions and a sound propagation path. Assuming that sensors are present respectively at two points A(x
0
, y
0
) and B(x
1
, y
1
) in a rectangular sectional area as shown in
FIG. 1
, calculation can be made with respect to time required for an acoustic wave emitted from the sensor at the point A to arrive at the sensor at the point B. Although the actual travel path of the acoustic wave is slightly curved from a straight line between the two points A and B because it is refracted due to a temperature gradient, it will be assumed to be a straight line. Hence, an arbitrary point on the path between the two points A and B can be expressed as in the following equation 3.
x
(
s
)=
x
0
+(
x
1
−x
0
)
s
y
(
s
)=
y
0
+(
y
1
−y
0
)
s
[Equation 3]
where,
s
=
l
L
,
[
Equation
3
]
L={square root over ((x
0
−x
1
)
2
+(y
0
−y
1
)
2
)} (0≦s≦1).
In the above equation 3, L is the distance between the two points A and B and l is the distance between the arbitrary point on the path and the point A. Dividing the x-coordinate system and y-coordinate system in the equation 3 respectively by lengths X and Y, dimensionless coordinate systems u and v can be obtained as in the below equation 4.
u
=
x
X
,
v
=
y
Y
[
Equation
4
]
where, u(s)=u
0
+(u
1
−u
0
)s, v(s)=v
0
+(v
1
−v
0
)s (0≦s≦1).
Because the temperature is a function of the position (x,y), or (u,v), the acoustic velocity is defined as a function of the position, too.
Accordingly, substituting (x, y) in the above equation 2 with (u, v) and integrating both sides on the basis of ds (or dl), the travel time of the sound can be expressed as in the following equation 5.
t
=
L
∫
0
1
ⅆ
s
C
(
u
,
v
)
=
L
∫
0
1
F
(
u
,
v
)
ⅆ
s
[
Equation
5
]
where,
F
(
u
,
v
)
=
1
C
(
u
,
v
)
.
Expressing F(u, v) in the above equation 5 as a Fourier series, the result is:
F
(
u
,
v
)
=
1
K
T
(
u
,
v
)
=
∑
m
∑
n
A
mn
G
(
mu
)
H
(
nv
)
[
Equation
6
]
In the above equation 6, G and H are admissible functions, which can be defined respectively as G=cos(mu) and H=cos(nv) in a rectangular sectional area. Substituting the equation 5 with the equation 6, the result is:
t
=
∑
m
∑
n
A
mn
f
mn
[
Equation
7
]
where, f
mn
=L ∫ G(mu)H(nv)ds, and m and n are dummy indexes.
Assuming that the number of sensors is N and the number of paths of each of the sensors is p, respective travel times t can be experimentally obtained as p simultaneous equations (p equations 7), which can be expressed in matrix as in the below equation 8.
[
f]
p×q
{A}
q×l
={t}
p×l
[Equation 8]
In the above equation 8, {A}={A
00
, A
01
, A
10
, A
02
, A
11
, A
20
, . . . }, and q is the number of coefficient terms taken for sufficient convergence of the Fourier coefficient A and must be smaller than or equal to p.
A Fourier coefficient vector {A} can be obtained from the measured time as in the following equation 9.
{[
f]
T
}
q×p
[f]
p×q
{A}
q×l
={[f]
T
}
q×p
{t}
p×l
{A}
q×l
={{[f]
T
}
q×p
[f]
p×q
}
q×q
−1
{{[f]
T
}
q×p
{t}
p×l
}
q×l
[Equation 9]
For example, for temperature measurement on a rectangular plane, f
mn
is expressed as in the following equation 10.
1
)
If
m
=
n
=
0
,
f
m
=
L
.
2
)
If
m
(
u
1
-
u
0
)
-
n
(
v
1
-
v
0
)
=
0
,
f
mn
=
L
2
{
cos
(
mu
0
-
nv
0
)
+
sin
(
mu
1
+
nv
1
)
-
sin
(
mu
0
+
nv
0
)
m
(
u
1
-
u
0
)
+
n
(
v
1
-
v
0
)
}
.
&IndentingNewLine;
3
)
If
m
(
u
1
-
u
0
)
+
n
(
v
1
-
v
0
)
=
0
,
f
mn
=
L
2
{
cos
(
mu
0
+
nv
0
)
+
sin
(
mu
1
-
nv
1
)
-
sin
(
mu
0
-
nv
0
)
m
(
u
1
-
u
0
)
-
n
(
v
1
-
v
0
)
}
.
4
)
In
other
cases
,
f
mn
=
L
2
{
sin
(
mu
1
+
nv
1
)
-
sin
(
mu
0
+
nv
0
)
m
(
u
1
-
u
0
)
+
n
(
v
1
-
v
0
)
}
+
sin
(
mu
1
-
nv
1
)
-
sin
(
mu
0
-
nv
0
)
m
(
u
1
-
u
0
)
-
n
(
v
1
-
v
0
)
.
[
Equation
&i
Ih Jeong Guon
Park Chul Min
LandOfFree
Method for optimally determining sensor positions for... does not yet have a rating. At this time, there are no reviews or comments for this patent.
If you have personal experience with Method for optimally determining sensor positions for..., we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Method for optimally determining sensor positions for... will most certainly appreciate the feedback.
Profile ID: LFUS-PAI-O-3075449