CN112268622A - Simultaneous reconstruction algorithm for flame three-dimensional temperature and soot volume fraction distribution - Google Patents
Simultaneous reconstruction algorithm for flame three-dimensional temperature and soot volume fraction distribution Download PDFInfo
- Publication number
- CN112268622A CN112268622A CN202011099349.7A CN202011099349A CN112268622A CN 112268622 A CN112268622 A CN 112268622A CN 202011099349 A CN202011099349 A CN 202011099349A CN 112268622 A CN112268622 A CN 112268622A
- Authority
- CN
- China
- Prior art keywords
- flame
- algorithm
- temperature
- volume fraction
- new
- 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.)
- Granted
Links
- 239000004071 soot Substances 0.000 title claims abstract description 38
- 230000005855 radiation Effects 0.000 claims abstract description 74
- 238000000034 method Methods 0.000 claims abstract description 39
- 239000000779 smoke Substances 0.000 claims abstract description 31
- 230000003595 spectral effect Effects 0.000 claims abstract description 25
- 230000008569 process Effects 0.000 claims abstract description 20
- 230000005540 biological transmission Effects 0.000 claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims abstract description 11
- 239000013598 vector Substances 0.000 claims description 46
- 238000010521 absorption reaction Methods 0.000 claims description 44
- 238000000137 annealing Methods 0.000 claims description 23
- 239000011159 matrix material Substances 0.000 claims description 23
- 238000001228 spectrum Methods 0.000 claims description 9
- 230000005457 Black-body radiation Effects 0.000 claims description 6
- 238000012804 iterative process Methods 0.000 claims description 4
- 239000007787 solid Substances 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 238000002922 simulated annealing Methods 0.000 abstract description 11
- 230000008901 benefit Effects 0.000 abstract description 8
- 238000005457 optimization Methods 0.000 description 12
- 238000009529 body temperature measurement Methods 0.000 description 7
- 239000002245 particle Substances 0.000 description 7
- 238000002485 combustion reaction Methods 0.000 description 6
- 238000005259 measurement Methods 0.000 description 5
- 238000000691 measurement method Methods 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000008033 biological extinction Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000002939 conjugate gradient method Methods 0.000 description 1
- 238000002247 constant time method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000000295 emission spectrum Methods 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000019637 foraging behavior Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/0014—Radiation pyrometry, e.g. infrared or optical thermometry for sensing the radiation from gases, flames
- G01J5/0018—Flames, plasma or welding
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/52—Radiation pyrometry, e.g. infrared or optical thermometry using comparison with reference sources, e.g. disappearing-filament pyrometer
- G01J5/53—Reference sources, e.g. standard lamps; Black bodies
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/80—Calibration
- G01J5/802—Calibration by correcting for emissivity
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Engineering & Computer Science (AREA)
- Plasma & Fusion (AREA)
- Radiation Pyrometers (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
The invention discloses a flame three-dimensional temperature and soot volume fraction distribution simultaneous reconstruction algorithm, which comprises the following steps: inputting emergent spectral radiation intensity information of flame radiation rays in different directions, dispersing flame micro-elements, and establishing a flame radiation transmission equation; and setting a target function, performing iterative calculation by using a constructed simultaneous reconstruction algorithm, and reconstructing to obtain flame three-dimensional temperature and smoke black volume fraction distribution. In the algorithm, in the iteration process, the non-negative least square algorithm is nested in the simulated annealing algorithm, so that the search efficiency of the SA algorithm in the process of searching the target function is obviously improved. Meanwhile, the SA algorithm has the advantages of high quality, strong initial value robustness, universality, easiness in implementation and the like, so that the NNLS-SA algorithm has good global search characteristics and is not easy to fall into a local optimal value, and has high precision and good convergence in the process of simultaneously solving the three-dimensional flame temperature and smoke black volume fraction distribution.
Description
Technical Field
The invention relates to a flame three-dimensional temperature and soot volume fraction distribution simultaneous reconstruction algorithm, and belongs to the technical field of flame temperature measurement.
Background
The combustion process is an extremely complex process in which physical and chemical changes couple and interact with each other. The temperature, the volume fraction of soot and the spatial distribution of other related products in the combustion flame are directly related to the safety, effectiveness and cleanness of the energy conversion and utilization process, and the measurement work of the key parameters has important significance for the safe operation, energy conservation and emission reduction of equipment.
In recent years, the development of laser, photoelectron, information and other technologies greatly promotes the deep research and application of non-contact temperature measurement technology. The non-contact temperature measurement technology has the advantage of not interfering the flow field to be measured, and can be roughly divided into two main categories, namely an active non-contact temperature measurement method and a passive non-contact temperature measurement method. The active temperature measurement method is to apply external excitation signals such as light and sound to high-temperature combustion flame, and to measure the temperature by detecting the interaction result of the combustion process and the applied external signals. The passive temperature measurement method does not adopt any external signal, and only detects information generated in the combustion process to measure the flame temperature, mainly refers to a radiation imaging method. Non-invasive soot particle measurement methods in combustion flames can also be divided into two types according to the source of an optical signal in the measurement process, wherein one type is a measurement method based on laser or other externally added light sources, such as an extinction method and a laser-induced incandescent light method; another class is measurement methods based on the emission spectrum of the flame or soot particles themselves, such as bicolor methods and emission CT methods. Laser-based optical measurement methods are widely applied in the field of measurement of flame temperature and smoke volume fraction, however, because the methods all need an additional light source to give an initial signal, the methods bring inconvenience when applied to industrial field measurement, and complicated field conditions bring difficulty to the arrangement of the light source and the detector and the calibration of the light path. These problems can be effectively avoided by measuring the flame's own emission spectral signal to directly calculate the temperature and soot particle volume fraction distribution in the target flame.
The radiation image method is used for measuring the three-dimensional flame temperature and the volume fraction distribution of the smoke black, and the temperature and the volume fraction of the smoke black in a radiation transmission system are solved according to the boundary radiation intensity distribution, which is a typical inverse problem. The existing solving algorithm of the radiation transmission inverse problem mainly comprises the following four types: the first type is a traditional regularization method, which comprises a Tikhonov regularization algorithm, a truncated singular value decomposition algorithm, a least square QR decomposition algorithm and the like; the second category is gradient-based optimization algorithms, including conjugate gradient methods, Levenberg-Marquardt (L-M) algorithms, and the like; the third category is an intelligent optimization algorithm based on random search, which comprises a genetic algorithm, simulated annealing, a particle swarm algorithm, an ant colony algorithm and the like; the fourth category is hybrid optimization algorithms, i.e., combining two or more different types of algorithms into a hybrid algorithm to improve the overall performance of the algorithm.
The regularization algorithm is a type of algorithm which is widely applied in the radiation transmission inverse problem solution, but the method can only process simple inverse problems about the solution of a linear equation system. For the problem of simultaneous reconstruction of flame temperature and smoke volume fraction, the problem is a typical nonlinear problem, and therefore, the simultaneous reconstruction of flame temperature and smoke volume fraction distribution cannot be realized by adopting a regularization method only. Gradient-based optimization algorithms are a class of nonlinear optimization algorithms. The algorithm has the advantages of high convergence rate and good result stability, but the algorithm generally has the defects of dependence on an initial value, difficulty in processing a local optimal value and the like. Compared with the traditional method, the intelligent optimization algorithm is developed later, but has the advantages of not depending on an initial value, not needing to differentiate the optimized function and the like.
Therefore, by taking advantage of different algorithms, the different algorithms are organically combined to form a hybrid optimization algorithm, respective advantages are exerted, and a new idea is provided for simultaneously solving the volume fraction distribution of the temperature and the smoke black. Although the intelligent optimization algorithm provides more possibilities for the organic combination of different algorithms, due to the defects of premature convergence of a genetic algorithm, stagnation of an ant colony algorithm searching process, low precision of a particle swarm algorithm, easiness in dispersion and the like, the application of the existing hybrid optimization algorithm in the research of the radiation transmission inverse problem is less, and particularly, the problem of simultaneous reconstruction of flame three-dimensional temperature and smoke black volume fraction distribution is rarely researched, so that the proposal of a new hybrid optimization algorithm and the improvement of the existing hybrid algorithm are still to be further researched.
Disclosure of Invention
The invention aims to solve the technical problem of providing an algorithm for simultaneously reconstructing three-dimensional flame temperature and smoke black volume fraction distribution aiming at the defects of the existing solving algorithm.
In order to solve the technical problems, the invention adopts the technical scheme that:
an algorithm for simultaneously reconstructing three-dimensional flame temperature and smoke black volume fraction distribution is characterized by comprising the following steps:
firstly, inputting emergent spectral radiation intensity information of flame radiation rays in different directions, dispersing flame micro-elements, and establishing a flame radiation transmission equation;
the first step comprises the following steps:
(11) the emergent spectral radiation intensity information of the input flame radiation rays in different directions comprises Iλ(R)(x, y, z, θ, ψ) and Iλ(G)(x,y,z,θ,ψ);
Wherein, Iλ(R)(x, y, z, theta, psi) as a starting point coordinate (x, y, z), a zenith angle of a direction coordinate theta, and a circumferential angle psi; i isλ(G)(x, y, z, theta, psi) as initial point coordinates (x, y, z), zenith angle of direction coordinates theta, and circumferential angle psi;
(12) dispersing the flame micro-elements, and performing three-dimensional grid division on the flame micro-elements;
(13) recording the length and the number of each light ray passing through different flame micro-elements according to the radiation intensity information input in the step (11), and establishing the following equation:
in the formula, nKRepresenting the total number of flame micro-elements penetrated by the light with the serial number K; k is lightThe serial number of (2); i isλK(s) is the lambda spectral radiation intensity of the radiation ray with sequence number K in the s direction;is serial number nKThe absorption coefficient of the infinitesimal body of (2); i.e. ip,jpRepresents the flame micro-element body through which the light passes along the transmission direction;the black body radiation intensity of the lambda spectrum of the flame radiation light with the serial number K passing through the vn-th flame micro element body;the light ray has a serial number nKThe geometric length of the flame micro-elements of (a); λ is the wavelength of the flame radiation;
the radiation transmission equation of all directions of light rays is combined to obtain the following equation:
Iλ=AIbλ (2)
wherein, IλVector [ I ] formed by spectral emergent radiation intensity of flame in all directionsλ1,Iλ2,…,Iλm,…,IλK](ii) a A is a coefficient matrix formed by calculating coefficients A, Ibλvector formed by black body radiation intensity of flame micro-elementV is the total number of flame control bodies (number of each infinitesimal body is 1,2, V, …, V), the total number of flame control bodies passing through is n for the flame radiation ray with the serial number K, and the number of each infinitesimal body passing through along the ray direction is V1, V2, vip,…,vn。
The detailed calculation formula of A is shown as formula (3),
setting a target function, performing iterative calculation by using a constructed simultaneous reconstruction algorithm (formula 4-14), and reconstructing to obtain flame three-dimensional temperature and smoke black volume fraction distribution;
(21) setting an objective function F, and initializing an absorption coefficient vector k ═ k0Inner loop maximum iteration number loopmaxThe threshold epsilon for whether the iteration process is terminated, the number of initialization inner loops Loop is 1, the number of outer loops Loop is 1, and the initial annealing temperature T of the annealing plans,0;
F=||I′λ(G)-Iλ(G)|| (4)
Wherein, Iλ(G)、I′λ(G)Respectively the measured and calculated spectral radiation intensity of G wavelength and the initial value vector kappa of the absorption coefficient0Is a vector which is randomly and continuously distributed;
(22) initializing a coefficient matrix a ═ a according to equation (3)0Calculating an initial objective function Fg,
Fg=||A0Ibλ(G)-Iλ(G)|| (5)
Randomly selecting the p-th component kappa in the absorption coefficient vector kappapDisturbing by using the formulas (6) and (7) to obtain a new absorption coefficient vector kappanewAnd updating the coefficient matrix A according to the formula (3) to obtain Anew:
κnew,p=κp+delta (6)
Wherein delta follows a standard normal distribution of equation 7:
delta~N(0,1) (7)
by the updated coefficient matrix AnewAnd the radiation intensity I of the blackbody spectrum under the current iteration numberbλ(G)Calculating an estimate I 'of the radiation intensity vector at the flame G wavelength for the current iteration number'λ(G)(κnew):
I′λ(G)(κnew)=AnewIbλ(G) (8)
Calculating a current state objective function Fl:
Fl=||I′λ(G)(κnew)-Iλ(G)|| (9)
Wherein I in the iterative processbλ(G)Obtained by the following method:
firstly, the spectral radiation intensity vector I of the measured R wavelength is utilizedλ(R)The vector kappa of the absorption coefficient updated in the iterative processnewCoefficient matrix AnewSolving the formula (2) by using a non-negative least square algorithm to obtain the blackbody spectral radiation intensity I of each infinitesimal body of the flame under the R wavelengthbλ(R)Solving the temperature value T of each infinitesimal body of the flame by the formula (10)f;
Tf=c2/λ(R)ln{c1/[λ(R)5πIbλ(R)]+1} (10)
In the formula, c1Is a first radiation constant, c2Is a second radiation constant; t isfIs the temperature of each flame infinitesimal body; λ (R) is the R wavelength of the light;
then, the temperature of each flame infinitesimal body obtained by the formula (10) is used for calculating the blackbody spectral radiation intensity I of each flame infinitesimal body under the G wavelength according to the formula (11)bλ(G);
Wherein λ (G) is the G wavelength of light;
the solving steps of the non-negative least square algorithm are as follows:
the method comprises the following steps: setting the sequence gamma as an empty set, and setting the sequence Z as [1,2,3 …, M]And I isbλ=0。
Step two: calculating M-dimensional vector, eta ═ AT(Iλ-AIbλ)。
Step three: if the sequence Z is an empty set or for all sequences i e Z, ω is satisfiediLess than or equal to 0, turn to
And step eleven, otherwise, carrying out the next step.
Step four: find the exponent q, so that ηq=max{ωi:i∈Z}。
Step five: the exponent q is shifted from the sequence Z into the sequence Γ.
Step six: a. theΓDenoted as an S M matrix. If i ∈ Γ, the ith column of matrix A is the new matrix AΓColumn i. If i ∈ Z, then the new matrix AΓIs a zero vector and calculates a least squares problemDefines only the solution Z satisfying the correspondence of the element i in the sequence ΓiZ corresponding to element i in the sequence ZiIs a zero solution.
Step seven: if all Z corresponding to the element i in the sequence gamma are satisfiediAre both greater than 0, then IbλZ and go to step two, otherwise proceed to the next step.
Step eight: the exponent u is found in the sequence Γ elements so that it satisfies the following relationship.
Step eleven: and finishing the calculation.
(23) Calculating Δ F ═ Fl-FgUpdating the objective function value FgAnd an absorption coefficient vector κ;
if Δ F<0, the annealing process is in the stateAs an important state, the absorption coefficient vector k ═ knewAn objective function Fg=Fl;
If Δ F>If the state is 0 and the state is important, the determination is made based on the probability w that the solid is in the state, and the calculation of the probability function is shown in equation (12). Generating a [0,1 ] by a random data generator]Random number xi of interval, if w>Xi, the new state is the important state, and the updated absorption coefficient vector k is equal to knewAn objective function Fg=Fl(ii) a Otherwise, the absorption coefficient vector and the target function are not updated;
w=exp[-ΔF/(kbTs)] (12)
wherein k isbIs the Boltzmann constant, TsAn annealing temperature for an annealing plan;
(24) judging whether the inner circulation is stopped:
judging whether the loop of the iteration times of the internal loop is less than the loop of the maximum iteration times of the internal loopmaxAnd if so, looping to step (22); otherwise, stopping the internal circulation and carrying out the next step;
(25) judging whether the external circulation is stopped:
judging whether the iteration reaches the convergence condition according to the formula (13), if not, updating the outer Loop frequency Loop to Loop +1, and updating the annealing temperature T of the annealing plan according to the annealing mode of the formula (14)sLooping to step (22); otherwise, the circulation is stopped, and the next step is carried out;
Ts=Ts,0·αLoop (14)
in the formula, α is a temperature decay rate, and is usually selected to be 0.7. ltoreq. α.ltoreq.1.0.
(27) Calculating the temperature and the volume fraction of the smoke black, and outputting the result:
calculating the temperature T of each micro element of the flame when the iteration is stopped according to the formula (13)endAnd absorption coefficient kappaendOutputting the temperature and the absorption coefficient, and according to the following formula,calculating the volume fraction of the smoke black;
a=1.811+0.1263lnλ+0.0217ln2λ+0.0417ln3λ (17)
b=0.5821+0.1213lnλ+0.2309ln2λ+0.0011ln3λ
in the formula (f)vIs the soot volume fraction of the flame infinitesimal body; kappaendIs the absorption coefficient of the flame infinitesimal body; and E (m) is the complex refractive index m of the soot varying with the wavelength, which is a function of a-bi, wherein a and b are the real part and the imaginary part of the complex refractive index m of the soot, respectively.
Has the advantages that:
1. the problem of simultaneous reconstruction of flame three-dimensional temperature and soot volume fraction distribution is a typical nonlinear solving problem, and the temperature and soot volume fraction need to be decoupled in the solving process, and meanwhile, the accuracy of the solving result needs to be ensured. The invention utilizes a Simulated Annealing (SA) algorithm to carry out global search on the objective function value, and simultaneously embeds a Non-Negative Least square (NNLS) algorithm in the SA algorithm to solve the flame temperature field in the searching process so as to obtain the reconstruction result of high-precision temperature and smoke volume fraction.
In the global search process of the SA algorithm, the absorption coefficient values of one or more micro-elements are updated every time, so that the global search time is longer, but the NNLS algorithm embedded in the SA can ensure the non-negativity of the temperature reconstruction of each micro-element of the flame, and the search efficiency of the SA algorithm in the process of searching the minimum value of the target function is remarkably improved. Meanwhile, the SA algorithm has the advantages of high quality, strong initial value robustness, universality, easiness in implementation and the like, so that the NNLS-SA algorithm has good global search characteristics and is not easy to fall into a local optimal value, and has high precision and good convergence in the process of simultaneously solving the three-dimensional flame temperature and smoke black volume fraction distribution.
Drawings
FIG. 1 is a flow chart of an algorithm for simultaneous reconstruction of flame three-dimensional temperature and soot volume fraction distribution;
FIG. 2 is a comparison of reconstructed flame temperatures using the NNLS-SA and NNLS-PSO algorithms, respectively;
FIG. 3 is a comparison of the reconstructed flame soot volume fractions using the NNLS-SA and NNLS-PSO algorithms, respectively.
Detailed Description
The invention is further illustrated with reference to the following figures and specific examples. It is to be understood that these examples are for illustrative purposes only and are not intended to limit the scope of the present invention, which is to be given the full breadth of the appended claims and any and all equivalent modifications thereof which may occur to those skilled in the art upon reading the present specification.
Step one, setting flame parameters, inputting emergent spectral radiation intensity information of flame radiation rays in different directions, dispersing flame micro elements, and establishing a flame radiation transmission equation. Assuming a cylindrical flame model, the height H is 0.03m, the bottom radius R is 0.008m, and the flame internal temperature and soot volume fraction distribution obey equations (1) and (2). Divide the flame into(V64) infinitesimals. Selecting 10000 emergent spectrum radiation intensity signals I of light rays in different directionsλ(R)(x, y, z, θ, ψ) and Iλ(G)(x, y, z, θ, ψ); wherein, (x, y, z) is the initial point coordinate of the light ray, x is-0.0160, the variation range of y is-0.0074, and the variation range of z is 0-0.03. (theta, psi) is the direction coordinate of the light, theta is the zenith angle, psi is the circumferential angle, theta is in the range of-5-1.2 deg., psi is in the range of 80-92 deg., lambda (R) and lambda (G) are respectively the wavelengths of R and G channels 610nm and 530 nm;
according to the input radiation intensity information, recording the length and the number of each light ray passing through different flame micro-elements, and establishing the following equation:
in the formula, nKRepresenting the total number of flame micro-elements penetrated by the light with the serial number K; k is the serial number of the ray; i isλK(s) is the radiation intensity of the lambda spectrum of the radiation ray with the serial number K in the s direction;is serial number nKThe absorption coefficient of the infinitesimal body of (2); i.e. ip,jpRepresents the flame micro-element body through which the light passes along the transmission direction;is the black body radiation intensity of lambda spectrum of flame radiation light with serial number K passing through the vn-th flame micro element,the light ray has a serial number nKThe geometric length of the flame micro-elements of (a); λ is the wavelength of the flame radiation; wherein the absorption coefficient kappaλAnd volume fraction f of sootvSee formulas (4) - (6);
a=1.811+0.1263lnλ+0.0217ln2λ+0.0417ln3λ (6)
b=0.5821+0.1213lnλ+0.2309ln2λ+0.0011ln3λ
where e (m) is the complex refractive index m of soot as a function of wavelength, a-bi, where a and b are the real and imaginary parts of the complex refractive index m of soot, respectively.
The radiation transmission equation of all directions of light rays is combined to obtain the following equation:
Iλ=A·Ibλ (7)
wherein, IλVector [ I ] formed by spectral emergent radiation intensity of flame in all directionsλ1,Iλ2,…,Iλm,…,IλK](ii) a A is a coefficient matrix formed by calculating coefficients A, Ibλvector formed by black body radiation intensity of flame micro-elementV is the total number of flame control bodies (number of each infinitesimal body is 1,2, V, …, V), the total number of flame control bodies passing through is n for the flame radiation ray with the serial number K, and the number of each infinitesimal body passing through along the ray direction is V1, V2, vip,…,vn。
The detailed calculation formula of A is shown as formula (8),
step two: and setting an objective function, and carrying out iterative calculation by using the constructed simultaneous reconstruction algorithm, and reconstructing three-dimensional flame temperature and smoke black volume fraction distribution at the same time.
(21) Setting an objective function F, as shown in equation (9), and initializing an absorption coefficient vector k ═ k0The number of inner loop loops is 1, and the maximum number of inner loop iterations loopmax50, the Loop number Loop is 1, and the threshold value epsilon for whether the iteration process is terminated is 10-5Annealing planned initial temperature Ts,0=2000K;
F=||I′λ(G)-Iλ(G)|| (9)
Wherein, Iλ(G)、I′λ(G)The initial value vector kappa of the spectral radiation intensity and the absorption coefficient obtained by measurement and forward calculation respectively0Are vectors that are randomly and continuously distributed.
(22) Initializing a coefficient matrix a ═ a according to equation (8)0Calculating an initial objective function Fg,
Fg=||A0Ibλ(G)-Iλ(G)|| (10)
Randomly selecting the p-th component kappa in the absorption coefficient vector kappapPerturbing by equations (11) and (12) to obtain a new absorption coefficient vector kappanewAnd updating the coefficient matrix A according to the formula (8) to obtain Anew:
κnew,p=κp+delta (11)
Wherein delta follows a standard normal distribution of equation 12:
delta~N(0,1) (12)
by the updated coefficient matrix AnewAnd the radiation intensity I of the blackbody spectrum under the current iteration numberbλ(G)Calculating an estimate I 'of the radiation intensity vector at the flame G wavelength for the current iteration number'λ(G)(κnew):
I′λ(G)(κnew)=AnewIbλ(G) (13)
Calculating a current state objective function Fl:
Fl=||I′λ(G)(κnew)-Iλ(G)|| (14)
The NNLS algorithm is realized by utilizing the lsqnanneg function of Matlab to solve the formula (7), and the blackbody spectral radiation intensity I of each infinitesimal body of the flame under the R wavelength is obtainedbλ(R)Solving the temperature value T of each infinitesimal body of the flame by the formula (15)f,
Tf=c2/λ(R)ln{c1/[λ(R)5πIbλ(R)]+1} (15)
In the formula, c1Is a first radiation constant, c2Is a second radiation constant of 3.7418 × 10-16W·m2And 1.4388 × 10-2m·K;TfIs the temperature of each flame infinitesimal body, K; λ (R) is the R wavelength of the light;
the blackbody spectral radiant intensity I of each flame element at the G wavelength was calculated from the temperature of each flame element obtained by the formula (15) according to the Planck's law (formula (16))bλ(G):
Wherein λ (G) is the G wavelength of light; according to the formula (7), the radiation intensity I is determined by the current coefficient matrix A and the blackbody spectrumbλ(G)Calculating an estimate I 'of the radiation intensity vector at the current flame G wavelength'λ(G)(κnew) Calculating the current state objective function F according to equation (14)l。
(24) Calculating Δ F ═ Fl-FgUpdating the objective function value FgAnd an absorption coefficient vector κ;
if Δ F<0, the state received by the annealing process is an important state, and the absorption coefficient vector k is equal to knewAn objective function Fg=Fl;
If Δ F>If the state is 0 and the state is important, the determination is made based on the probability w that the solid is in the state, and the calculation of the probability function is shown in equation (17). Generating a [0,1 ] by a random data generator]Random number xi of interval, if w>Xi, the new state is the important state, and the updated absorption coefficient vector k is equal to knewAn objective function Fg=Fl(ii) a Otherwise, the absorption coefficient vector and the target function are not updated;
w=exp[-ΔF/(kbTs)] (17)
wherein k isbIs the Boltzmann constant, kb=1.3806488×10-23,TsAn annealing temperature for an annealing plan;
(25) judging whether the inner circulation is stopped:
judging whether the iteration number k of the inner loop is less than the maximum iteration number k of the inner loopmaxAnd if so, looping to step (22); otherwise, stopping the inner loop and proceeding to the next step;
(26) judging whether the external circulation is stopped:
judging whether the iteration reaches the convergence condition according to the equation (18), if not, updating the outer Loop frequency Loop to Loop +1, and updating the annealing temperature T of the annealing plan according to the annealing mode of the equation (19)sLooping to step (22); otherwise, the circulation is stopped, and the next step is carried out;
Ts=Ts,0·αLoop (19)
in the formula, α is a temperature decay rate, and is usually selected to be 0.7. ltoreq. α.ltoreq.1.0.
(27) Calculating the temperature and the volume fraction of the smoke black, and outputting the result:
calculating the temperature T of each micro element of the flame at the time of iteration termination according to the formula (18)endAnd absorption coefficient kappaendOutputting the temperature and the absorption coefficient, and calculating the volume fraction of the smoke black according to the following formula;
in the formula (f)vIs the volume fraction of soot of the flame micro-elements, ppm; kappaendIs the absorption coefficient of the flame infinitesimal body, m-1(ii) a E (m) is a function of the complex refractive index m of soot.
In order to evaluate the performance of the NNLS-SA algorithm, the NNLS-PSO algorithm is introduced for comparison, the temperature of flames and the volume fraction of soot are respectively reconstructed by using the NNLS-SA algorithm and the NNLS-PSO algorithm, wherein a Particle Swarm Optimization algorithm (PSO) is proposed according to the Swarm foraging behavior of birds, and is also an intelligent Optimization algorithm based on random search.
The reconstruction relative error of the flame soot volume fraction is calculated according to equation (21), the reconstruction relative error of the flame temperature is calculated according to equation (22), and the reconstruction results are shown in fig. 2 and 3.
In the formula (f)v,rstFor the reconstructed volume fraction of soot, fv,setFor a set volume fraction of soot, TrstFor reconstructed flame temperature values, TsetIs a set flame temperature value.
The NNLS-PSO algorithm mainly utilizes a Monte Carlo method, updates the whole value of the absorption coefficient vector kappa every time of iteration, searches a solution space through a particle group, and obtains the optimal values of the volume fraction of the smoke black and the temperature by comparing the historical optimal values of the objective function. In each iteration process of the NNLS-SA algorithm, the absorption coefficient values of one or more micro-elements are only updated, and then the change of the objective function before and after updating is compared to obtain the optimal value of the volume fraction of the smoke black and the temperature.
The NNLS-SA algorithm embeds the NNLS algorithm in the SA algorithm, guarantees the temperature property of temperature field solving, and meanwhile solves the volume fraction of the smoke black by utilizing the good global search characteristic of the SA algorithm to the nonlinear problem. Compared with the NNLS-PSO algorithm, the algorithm does not blindly update the absorption coefficients of all the infinitesimal bodies at random, but selects the absorption coefficients of one or more infinitesimal bodies to update, so that the algorithm has higher precision, is not easy to fall into a local extreme value, has better stability, and has higher solving precision of the NNLS-SA algorithm under the same iteration times. From the reconstruction results of the flame temperatures in fig. 2, it can be seen that in grids 46 to 48, the true values of the temperatures are 901K, the temperatures reconstructed using the NNLS-SA algorithm are 894K, 894K and 901K, respectively, and the temperatures reconstructed using the NNLS-PSO algorithm are 863K, 862K and 863K, respectively, and it can be seen that the flame temperatures reconstructed using the NNLS-SA algorithm are more accurate. It can also be seen from the reconstructed results of the volume fractions of the soot in fig. 3 that in grids No. 10 to 12, the true values of the volume fractions of the soot are 0.8031ppm, the volume fractions of the flame soot reconstructed by the NNLS-SA algorithm are 0.8197ppm, 0.8042ppm and 0.8082ppm, respectively, and the volume fractions of the flame soot reconstructed by the NNLS-PSO algorithm are 0.7283ppm, 0.9224ppm and 1.3466ppm, respectively, and the reconstructed results of the NNLS-SA algorithm on the volume fractions of the soot are more accurate than those reconstructed by the NNLS-PSO algorithm. Meanwhile, the average relative errors of the reconstructed temperature and the volume fraction of the smoke black of the two algorithms are calculated, the average relative errors of the reconstructed temperature and the volume fraction of the smoke black of the NNLS-SA algorithm are respectively 0.7% and 6.05%, and the average relative errors of the reconstructed temperature and the volume fraction of the smoke black of the NNLS-PSO algorithm are respectively 10.87% and 43.85%. In conclusion, the NNLS-SA algorithm has higher reconstruction precision no matter the reconstruction of the temperature or the volume fraction of the smoke black.
Claims (7)
1. An algorithm for simultaneously reconstructing three-dimensional flame temperature and smoke black volume fraction distribution is characterized by comprising the following steps:
firstly, inputting emergent spectral radiation intensity information of flame radiation rays in different directions, dispersing flame micro-elements, and establishing a flame radiation transmission equation;
and step two, setting a target function, performing iterative calculation by using a constructed simultaneous reconstruction algorithm, and reconstructing to obtain flame three-dimensional temperature and smoke black volume fraction distribution.
2. The algorithm for simultaneously reconstructing flame three-dimensional temperature and soot volume fraction distribution according to claim 1, wherein the first step comprises:
(11) the emergent spectral radiation intensity information of the input flame radiation rays in different directions comprises emergent spectral radiation intensity Iλ(R)(x, y, z, theta, psi) and the intensity of the emergent spectral radiation Iλ(G)(x,y,z,θ,ψ);
Wherein, Iλ(R)(x, y, z, theta, psi) as a starting point coordinate (x, y, z), a zenith angle of a direction coordinate theta, and a circumferential angle psi; i isλ(G)(x, y, z, theta, psi) as initial point coordinates (x, y, z), zenith angle of direction coordinates theta, and circumferential angle psi;
(12) dispersing the flame micro-elements, and performing three-dimensional grid division on the flame micro-elements;
(13) recording the length and the number of each light ray passing through different flame micro-elements according to the radiation intensity information input in the step (11), and establishing the following equation:
in the formula, nKRepresenting the total number of flame micro-elements penetrated by the light with the serial number K; k is the serial number of the ray; i isλK(s) is the lambda spectral radiation intensity of the radiation ray with sequence number K in the s direction;is serial number nKThe absorption coefficient of the infinitesimal body of (2); i.e. ip,jpRepresents the flame micro-element body through which the light passes along the transmission direction;the black body radiation intensity of the lambda spectrum of the flame radiation light with the serial number K passing through the vn-th flame micro element body;the light ray has a serial number nKThe geometric length of the flame micro-elements of (a); λ is the wavelength of the flame radiation;
the radiation transmission equation of all directions of light rays is combined to obtain the following equation:
Iλ=AIbλ (2)
wherein, IλVector [ I ] formed by spectral emergent radiation intensity of flame in all directionsλ1,Iλ2,…,Iλm,…,IλK](ii) a A is a coefficient matrix formed by calculating coefficients A, Ibλvector formed by black body radiation intensity of flame micro-elementV is the total number of flame control bodies, each infinitesimal body is numbered 1,2, V, …, V, for the flame radiation ray with the serial number K, the total number of the penetrated flame control bodies is n, and the number of the penetrated infinitesimal bodies along the ray direction is V1, V2, vip,…,vn。
4. the algorithm for simultaneously reconstructing flame three-dimensional temperature and soot volume fraction distribution according to claim 3, wherein the second step comprises:
(21) setting an objective function F, and initializing an absorption coefficient vector k ═ k0Inner loop maximum iteration number loopmaxThe threshold epsilon for whether the iteration process is terminated, the number of initialization inner loops Loop is 1, the number of outer loops Loop is 1, and the initial annealing temperature T of the annealing plans,0;
F=||I′λ(G)-Iλ(G)|| (4)
Wherein, Iλ(G)、I′λ(G)Respectively the measured and calculated spectral radiation intensity of G wavelength and the initial value vector kappa of the absorption coefficient0Is a vector which is randomly and continuously distributed;
(22) initializing a coefficient matrix a ═ a according to equation (3)0Calculating an initial objective function Fg,
Fg=||A0Ibλ(G)-Iλ(G)|| (5)
Randomly selecting the p-th component kappa in the absorption coefficient vector kappapDisturbing by using the formulas (6) and (7) to obtain a new absorption coefficient vector kappanewAnd updating the coefficient matrix A according to the formula (3) to obtain Anew:
κnew,p=κp+delta (6)
Wherein delta follows a standard normal distribution of equation 7:
delta~N(0,1) (7)
by the updated coefficient matrix AnewAnd the radiation intensity I of the blackbody spectrum under the current iteration numberbλ(G)Calculating an estimate I 'of the radiation intensity vector at the flame G wavelength for the current iteration number'λ(G)(κnew):
I′λ(G)(κnew)=AnewIbλ(G) (8)
Calculating a current state objective function Fl:
Fl=||I′λ(G)(κnew)-Iλ(G)|| (9)
Wherein I in the iterative processbλ(G)Obtained by the following method:
firstly, the spectral radiation intensity vector I of the measured R wavelength is utilizedλ(R)The vector kappa of the absorption coefficient updated in the iterative processnewCoefficient matrix AnewSolving the formula (2) by using a non-negative least square algorithm to obtain the blackbody spectral radiation intensity I of each infinitesimal body of the flame under the R wavelengthbλ(R)By the formula (10),solving the temperature value T of each infinitesimal body of the flamef;
Tf=c2/λ(R)ln{c1/[λ(R)5πIbλ(R)]+1} (10)
In the formula, c1Is a first radiation constant, c2Is a second radiation constant; t isfIs the temperature of each flame infinitesimal body; λ (R) is the R wavelength of the light;
then, the temperature of each flame infinitesimal body obtained by the formula (10) is used for calculating the blackbody spectral radiation intensity I of each flame infinitesimal body under the G wavelength according to the formula (11)bλ(G);
Wherein λ (G) is the G wavelength of light;
(23) calculating Δ F ═ Fl-FgUpdating the objective function value FgAnd an absorption coefficient vector κ;
if Δ F<0, the state received by the annealing process is an important state, and the absorption coefficient vector k is equal to knewAn objective function Fg=Fl;
If Δ F>0, judging whether the solid is accepted as an important state according to the probability w of the solid in the state; generating a [0,1 ] by a random data generator]Random number xi of interval, if w>Xi, the new state is the important state, and the updated absorption coefficient vector k is equal to knewAn objective function Fg=Fl(ii) a Otherwise, the absorption coefficient vector and the target function are not updated;
(24) judging whether the inner circulation is stopped:
judging whether the loop of the iteration times of the internal loop is less than the loop of the maximum iteration times of the internal loopmaxAnd if so, looping to step (22); otherwise, stopping the internal circulation and carrying out the next step;
(25) judging whether the external circulation is stopped:
judging whether the iteration reaches the convergence condition according to the formula (12), and if not, updating the outer loop timesThe Loop +1 is set to the Loop, and the annealing temperature T of the annealing plan is updated in accordance with the annealing method of the formula (13)sLooping to step (22); otherwise, the circulation is stopped, and the next step is carried out;
Ts=Ts,0·αLoop (13)
wherein α is a temperature decay rate;
(26) calculating the temperature and the volume fraction of the smoke black, and outputting the result:
calculating the temperature T of each micro element of the flame when the iteration is stopped according to the formula (12)endAnd absorption coefficient kappaendOutputting the temperature and the absorption coefficient, and calculating the volume fraction of the smoke black according to the following formula;
in the formula (f)vIs the soot volume fraction of the flame infinitesimal body; kappaendIs the absorption coefficient of the flame infinitesimal body; and E (m) is the complex refractive index m of the soot varying with the wavelength, which is a function of a-bi, wherein a and b are the real part and the imaginary part of the complex refractive index m of the soot, respectively.
5. The algorithm for simultaneous reconstruction of flame three-dimensional temperature and soot volume fraction distribution according to claim 4, wherein the probability function w is:
w=exp[-ΔF/(kbTs)] (17)
wherein k isbIs the Boltzmann constant, TsThe annealing temperature for the annealing plan.
6. The algorithm for simultaneously reconstructing three-dimensional flame temperature and soot volume fraction distribution according to claim 4, wherein the temperature decay rate α is a value: alpha is more than or equal to 0.7 and less than or equal to 1.0.
7. The algorithm for simultaneous reconstruction of flame three-dimensional temperature and soot volume fraction distribution according to claim 4, wherein the solving of the non-negative least squares algorithm comprises the following steps:
the method comprises the following steps: setting the sequence gamma as an empty set, and setting the sequence Z as [1,2,3 …, M]And I isbλ=0;
Step two: calculating M-dimensional vector, eta ═ AT(Iλ-AIbλ);
Step three: if the sequence Z is an empty set or for all sequences i e Z, ω is satisfiediTurning to the step eleven if the temperature is less than or equal to 0, otherwise, carrying out the next step;
step four: find the exponent q, so that ηq=max{ωi:i∈Z};
Step five: moving the exponent q from the sequence Z into the sequence Γ;
step six: a. theΓDenoted as an S M matrix. If i ∈ Γ, the ith column of matrix A is the new matrix AΓColumn i. If i ∈ Z, then the new matrix AΓIs a zero vector and calculates a least squares problemDefines only the solution Z satisfying the correspondence of the element i in the sequence ΓiZ corresponding to element i in the sequence ZiIs a zero solution;
step seven: if all Z corresponding to the element i in the sequence gamma are satisfiediAre both greater than 0, then IbλZ, and go to step two, otherwise go to the next step;
step eight: finding the exponent u in the sequence Γ elements such that it satisfies the following relationship:
step eleven: and finishing the calculation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011099349.7A CN112268622B (en) | 2020-10-14 | 2020-10-14 | Simultaneous reconstruction algorithm for flame three-dimensional temperature and soot volume fraction distribution |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011099349.7A CN112268622B (en) | 2020-10-14 | 2020-10-14 | Simultaneous reconstruction algorithm for flame three-dimensional temperature and soot volume fraction distribution |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112268622A true CN112268622A (en) | 2021-01-26 |
CN112268622B CN112268622B (en) | 2021-12-24 |
Family
ID=74338091
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011099349.7A Active CN112268622B (en) | 2020-10-14 | 2020-10-14 | Simultaneous reconstruction algorithm for flame three-dimensional temperature and soot volume fraction distribution |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112268622B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114046883A (en) * | 2021-10-27 | 2022-02-15 | 东南大学 | Method and device for simultaneously reconstructing flame temperature and soot concentration distribution and storage medium |
CN115014775A (en) * | 2022-04-27 | 2022-09-06 | 中国人民解放军国防科技大学 | Method for measuring temperature and emissivity of combustion chamber of scramjet engine |
CN115014776A (en) * | 2022-04-27 | 2022-09-06 | 中国人民解放军国防科技大学 | Method for reconstructing two-dimensional distribution of flame temperature and emissivity of combustion chamber of scramjet engine |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6850862B1 (en) * | 1999-10-22 | 2005-02-01 | Meta Instruments Srl | Methods and equipment for the measurement of the three-dimensional distribution of the temperatures within dielectric means |
CN102881041A (en) * | 2012-08-21 | 2013-01-16 | 中国科学院计算技术研究所 | Multi-source measured data-based flame modeling method and system |
CN107421644A (en) * | 2017-08-28 | 2017-12-01 | 南京大学 | The air remote sensing evaluation method of the complete surface temperature in city |
CN108225577A (en) * | 2017-12-19 | 2018-06-29 | 东南大学 | Flame three dimensional temperature distribution method for reconstructing based on double Spectral Radiation Informations |
-
2020
- 2020-10-14 CN CN202011099349.7A patent/CN112268622B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6850862B1 (en) * | 1999-10-22 | 2005-02-01 | Meta Instruments Srl | Methods and equipment for the measurement of the three-dimensional distribution of the temperatures within dielectric means |
CN102881041A (en) * | 2012-08-21 | 2013-01-16 | 中国科学院计算技术研究所 | Multi-source measured data-based flame modeling method and system |
CN107421644A (en) * | 2017-08-28 | 2017-12-01 | 南京大学 | The air remote sensing evaluation method of the complete surface temperature in city |
CN108225577A (en) * | 2017-12-19 | 2018-06-29 | 东南大学 | Flame three dimensional temperature distribution method for reconstructing based on double Spectral Radiation Informations |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114046883A (en) * | 2021-10-27 | 2022-02-15 | 东南大学 | Method and device for simultaneously reconstructing flame temperature and soot concentration distribution and storage medium |
CN114046883B (en) * | 2021-10-27 | 2023-12-22 | 东南大学 | Simultaneous reconstruction method and device for flame temperature and soot concentration distribution and storage medium |
CN115014775A (en) * | 2022-04-27 | 2022-09-06 | 中国人民解放军国防科技大学 | Method for measuring temperature and emissivity of combustion chamber of scramjet engine |
CN115014776A (en) * | 2022-04-27 | 2022-09-06 | 中国人民解放军国防科技大学 | Method for reconstructing two-dimensional distribution of flame temperature and emissivity of combustion chamber of scramjet engine |
Also Published As
Publication number | Publication date |
---|---|
CN112268622B (en) | 2021-12-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112268622B (en) | Simultaneous reconstruction algorithm for flame three-dimensional temperature and soot volume fraction distribution | |
Fu et al. | Temperature measurements of diesel fuel combustion with multicolor pyrometry | |
Sun et al. | Three-dimensional temperature field measurement of flame using a single light field camera | |
Huang et al. | Simultaneous reconstruction of 3D temperature distribution and radiative properties of participating media based on the multi-spectral light-field imaging technique | |
CN108225577B (en) | Flame three-dimensional temperature distribution reconstruction method based on double-spectrum radiation information | |
Xia et al. | Two-step tomographic reconstructions of temperature and species concentration in a flame based on laser absorption measurements with a rotation platform | |
Niu et al. | A novel parametric level set method coupled with Tikhonov regularization for tomographic laser absorption reconstruction | |
Ren et al. | Machine learning applied to retrieval of temperature and concentration distributions from infrared emission measurements | |
Liu et al. | Reconstruction of axisymmetric temperature and gas concentration distributions by combining fan-beam TDLAS with onion-peeling deconvolution | |
CN111238663B (en) | Method for measuring biomass flame temperature and emissivity based on Rayleigh approximation | |
CN101907492A (en) | Object emission spectrum-based temperature measurement method | |
Shi et al. | Simultaneous measurement of flame temperature and species concentration distribution from nonlinear tomographic absorption spectroscopy | |
Zhang et al. | A multispectral thermometry based on multi-objective constraint optimization | |
Zhang et al. | Nonlinear acoustic tomography for measuring the temperature and velocity fields by using the covariance matrix adaptation evolution strategy algorithm | |
Niu et al. | Three-dimensional rapid visualization of flame temperature field via compression and noise reduction of light field imaging | |
Liu et al. | Measurement of temperature and wavelength-dependent emissivity distributions using multi-wavelength radiation thermometry | |
Li et al. | Meshed axisymmetric flame simulation and temperature reconstruction using light field camera | |
Ertürk et al. | Inverse methods in thermal radiation analysis and experiment | |
Shui et al. | Linear multispectral absorption tomography based on regularized iterative methods | |
Gao et al. | Radial basis function coupled SART method for dynamic LAS tomography | |
Zhang et al. | Reconstruction of combustion temperature and gas concentration distributions using line-of-sight tunable diode laser absorption spectroscopy | |
CN114046883B (en) | Simultaneous reconstruction method and device for flame temperature and soot concentration distribution and storage medium | |
CN116893152A (en) | Online detection method and system for carbon dioxide concentration in flue gas | |
Zheng et al. | Measurement of temperature and emissivity of biomass candle flame using spectral thermometry | |
CN109100045B (en) | Gas temperature probability density distribution reconstruction method based on single light path multispectral |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |