CN110879301A - Method and system for simultaneously measuring two-dimensional distribution of liquid concentration and liquid movement - Google Patents

Method and system for simultaneously measuring two-dimensional distribution of liquid concentration and liquid movement Download PDF

Info

Publication number
CN110879301A
CN110879301A CN201910963960.0A CN201910963960A CN110879301A CN 110879301 A CN110879301 A CN 110879301A CN 201910963960 A CN201910963960 A CN 201910963960A CN 110879301 A CN110879301 A CN 110879301A
Authority
CN
China
Prior art keywords
liquid
concentration
light intensity
sub
dimensional distribution
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.)
Pending
Application number
CN201910963960.0A
Other languages
Chinese (zh)
Inventor
胡文成
张宝华
窦义涛
杨帅
杨桐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
AECC Shenyang Engine Research Institute
Original Assignee
AECC Shenyang Engine Research Institute
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by AECC Shenyang Engine Research Institute filed Critical AECC Shenyang Engine Research Institute
Priority to CN201910963960.0A priority Critical patent/CN110879301A/en
Publication of CN110879301A publication Critical patent/CN110879301A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/26Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting optical wave
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/30Measuring arrangements characterised by the use of optical techniques for measuring roughness or irregularity of surfaces
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/18Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the time taken to traverse a fixed distance

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Multimedia (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

The application relates to a method for simultaneously measuring two-dimensional distribution of liquid concentration and liquid movement, which comprises the following steps: generating at least two plane light beams which can be converged to form an effective measuring area parallel, wherein the effective measuring area comprises a plurality of sub-areas; respectively obtaining incident light intensity and emergent light intensity which are emitted by a light source and pass through the sub-regions, determining the liquid surface area in the sub-regions according to the incident light intensity and the emergent light intensity, and reconstructing the liquid surface area concentration in each voxel to form a first liquid surface concentration two-dimensional distribution and a second liquid surface concentration two-dimensional distribution; and mutually correlating the first liquid surface concentration two-dimensional distribution and the second liquid surface concentration two-dimensional distribution to form liquid motion two-dimensional distribution. The liquid-phase two-dimensional velocity measuring device can realize two-dimensional velocity measurement of a liquid phase on a certain section below the spray, is real transient measurement, can display the liquid concentration change trend in real time, expands the flow measurement function, and finally realizes simultaneous measurement of the concentration and the velocity of liquid.

Description

Method and system for simultaneously measuring two-dimensional distribution of liquid concentration and liquid movement
Technical Field
The application belongs to the technical field of imaging, and particularly relates to a method and a system for simultaneously measuring liquid concentration and liquid motion two-dimensional distribution.
Background
In an aircraft engine, a plurality of nozzles are provided, and the distribution of the concentration and the movement of liquid (fuel) sprayed by the nozzles needs to be measured in the design and test stages. By adopting the traditional SetScan technology, the method can only simultaneously measure the two-dimensional distribution of the surface diameter (SMD) of the plane and the surface area concentration of the particle diameter, but cannot measure the two-dimensional movement velocity distribution vertical to the axial plane of the spray, and has low time resolution and no function of recording images across frames.
Disclosure of Invention
It is an object of the present application to provide a method and system for simultaneously measuring a two-dimensional distribution of liquid concentration and liquid motion to solve or mitigate at least one of the problems of the background art.
In one aspect, the technical solution provided by the present application is: a method of simultaneously measuring a two-dimensional distribution of liquid concentration and liquid movement, the method comprising:
generating at least two parallel plane light beams, and converging the plane light beams to form an effective measuring area, wherein the effective measuring area is vertical to the flow direction of liquid, the effective measuring area comprises a plurality of sub-areas, and each plane light beam at least comprises two light sources which are sent at preset time intervals and converged;
acquiring incident light intensity and emergent light intensity which are emitted by at least one light source and pass through the sub-region, determining the liquid surface area in the sub-region according to the incident light intensity and the emergent light intensity, and reconstructing the liquid surface area concentration in each voxel to form first liquid surface concentration two-dimensional distribution;
acquiring incident light intensity and emergent light intensity which are emitted by at least one other light source and pass through the sub-region, determining the liquid surface area in the sub-region according to the incident light intensity and the emergent light intensity, and reconstructing the liquid surface area concentration in each voxel to form second liquid surface concentration two-dimensional distribution;
and mutually correlating the first liquid surface concentration two-dimensional distribution and the second liquid surface concentration two-dimensional distribution to form liquid motion two-dimensional distribution.
In an embodiment of the method of the present application, each planar light beam includes at least two light sources, and at least one light source is converged with another light source by refraction and then scattered to form a planar light beam.
In an embodiment of the method of the present application, the at least two light sources are emitted at predetermined time intervals by coordinated control.
In one embodiment of the method of the present application, the light source is a pulsed laser.
In an embodiment of the method of the present application, the incident light intensity and the emergent light intensity determine a liquid surface area in the sub-region, and are determined by the following relationship:
Iλe=Iλ0exp(-βλΔ)
in the formula IλeTo output a set light intensity, Iλ0β for incident light intensityλTo the attenuation coefficient, and in relation to the surface area of the liquid, Δ is the path taken by the light.
On the other hand, the technical scheme provided by the application is as follows: a system for simultaneously measuring a two-dimensional distribution of liquid concentration and liquid motion, the system comprising:
the device comprises at least two plane light beam forming devices, a light source, a liquid level sensor and a control device, wherein the at least two plane light beam forming devices are used for generating at least two parallel plane light beams and converging the plane light beams to form an effective measuring area, the effective measuring area is vertical to the flow direction of liquid, the effective measuring area comprises a plurality of sub-areas, and each plane light beam at least comprises two light sources which are sent according to a preset time interval and are converged;
the first liquid concentration distribution generating device is used for acquiring incident light intensity and emergent light intensity which are emitted by at least one light source and pass through the sub-area, determining the liquid surface area in the sub-area according to the incident light intensity and the emergent light intensity, and reconstructing the liquid surface area concentration in each voxel to form first liquid surface concentration two-dimensional distribution;
the second liquid concentration distribution generating device is used for acquiring incident light intensity and emergent light intensity which are emitted by at least one other light source and pass through the sub-area, determining the liquid surface area in the sub-area according to the incident light intensity and the emergent light intensity, and reconstructing the liquid surface area concentration in each voxel to form second liquid surface concentration two-dimensional distribution;
and the liquid motion distribution generating device is used for mutually correlating the first liquid surface concentration two-dimensional distribution and the second liquid surface concentration two-dimensional distribution to form liquid motion two-dimensional distribution.
In an embodiment of the system of the present application, the planar beam forming device includes:
at least two light source emitting modules for generating and emitting light sources;
at least one light source refraction module for reflecting or refracting the at least one light source at a predetermined angle to converge the light sources emitted by the at least two light source emission modules; and
the sheet light adjusting module is used for scattering the converged light source to form a plane light beam.
In an embodiment of the system of the present application, the planar beam forming apparatus further includes
And the control module is used for coordinating and controlling at least two light source emitting modules to emit light sources according to a preset time interval.
In an embodiment of the system of the present application, the light source generated by the light source emitting module is a pulsed laser.
In an embodiment of the system of the present application, the first liquid concentration generating device and/or the second liquid concentration generating device determines the liquid surface area in the sub-area according to the incident light intensity and the emergent light intensity by using the following relationship:
Iλe=Iλ0exp(-βλΔ)
in the formula IλeTo output a set light intensity, Iλ0β for incident light intensityλTo the attenuation coefficient, and in relation to the surface area of the liquid, Δ is the path taken by the light.
The method and the system for simultaneously measuring the liquid concentration and the liquid motion two-dimensional distribution can realize the two-dimensional speed measurement of the liquid phase of a certain section below the spray, can also display the liquid concentration change trend in real time for real transient measurement, expand the flow measurement function, and finally realize the simultaneous measurement of the concentration and the speed of the liquid.
Drawings
In order to more clearly illustrate the technical solutions provided by the present application, the following briefly introduces the accompanying drawings. It is to be expressly understood that the drawings described below are only illustrative of some embodiments of the invention.
FIG. 1 is a flow chart of a method for simultaneously measuring two-dimensional distributions of liquid concentration and liquid motion according to the present application.
FIG. 2 is a diagram of the system components for simultaneous measurement of two-dimensional distributions of liquid concentration and liquid motion according to the present application.
Fig. 3 is a schematic view of a scenario in which the present application is applied.
Fig. 4 is a sub-region diagram of an effective measurement region partition according to the present application.
Fig. 5 is a schematic diagram of effective measurement area division and imaging effect according to an embodiment of the present application.
Fig. 6 is a schematic diagram of an exposure timing sequence of a laser emitting and imaging device according to an embodiment of the present application.
FIG. 7 is a graph of two-plane liquid surface concentration distribution in an embodiment of the present application.
Fig. 8 is a vector diagram associated with the bipartite graph of fig. 7.
Detailed Description
In order to make the implementation objects, technical solutions and advantages of the present application clearer, the technical solutions in the embodiments of the present application will be described in more detail below with reference to the drawings in the embodiments of the present application.
In order to make the technical solution of the present application easier to understand, the following description will be made in conjunction with the method provided by the present application and the system provided by the present application.
As shown in fig. 1, the method 100 for simultaneously measuring two-dimensional distribution of liquid concentration and liquid movement provided by the present application includes the following steps:
step 110: generating at least two parallel plane light beams, converging the plane light beams to form an effective measuring area, wherein the effective measuring area is perpendicular to the flow direction of the liquid, the effective measuring area is divided into a plurality of sub-areas, and each plane light beam at least comprises two light sources which are sent at a preset time interval and are converged.
In order to generate at least two planar light beams according to the above method, at least two planar light beam forming devices 210 are provided in the system of the present application, the planar light beam forming devices 210 are used to generate planar light beams, and the planar light beams generated by each of the planar light beam forming devices 210 are parallel to each other and converge to form an effective measuring region 230.
In the embodiment shown in fig. 2, three sets of planar light beam forming devices 210 are given as an example, and in practical applications, only two planar light beam forming devices 210 are required to operate, so that when the number of planar light beam forming devices 210 is increased, the spatial resolution of the two-dimensional distribution measurement of the liquid concentration can be increased.
In the system of the present application, the planar light beam forming means 210 includes at least two light source emitting modules 211, at least one light source refracting module 212, and a sheet light adjusting module 213. Each light source emitting module 211 is for generating and emitting a light source. The light source refraction module 212 is configured to reflect or refract the light source generated by the at least one light source emission module 211 at a predetermined angle, so that the tube bundle generated by the light source emission module 211 is converged into one beam, and finally the light adjustment module 213 scatters the refracted light source to form a planar light beam.
In the embodiment shown in fig. 2, two light source emitting modules 211 and one light source refracting module 212 are provided. It is understood that in the system of the present application, the number of the light source emitting modules 211 may be three, four or more, and accordingly, the light source refraction module 212 may be added to make all the light emitted by the light source emitting modules 211 realize a common beam output.
In addition, since the light sources emitted by the light source emitting modules 211 in the present application need to be emitted at predetermined time intervals, the system in the present application further has a control module 214, and the control module 214 is connected to the light source emitting modules 211 to control each light source emitting module 211 to emit at predetermined time.
Further, the sheet light adjusting module 213 includes a first adjusting module 2131 and a second adjusting module 2132, the first adjusting module 2131 expands and expands the thinner common beam output light source to form a light source with a larger area, and the second adjusting module 2132 performs parallel constraint and outputs the expanded light source to form a planar light beam.
In some embodiments, the light source refraction module 212 is generally made of a material having a reflection or refraction function, such as reflective glass or a triangular prism. In some embodiments, the sheet light adjusting module 213 is an optical element made of a transparent material, such as a lens made of glass or crystal.
In order to make the light source easy to control and have better imaging effect, the light source in this application adopts a laser (pulse) light source, i.e. the light source emitting module 211 is a laser emitting device, such as a pulse laser diode. In order to ensure the stability of the emitted light intensity, microsecond-level pulse laser is adopted as a light source in the application. In order to ensure that the output wavelength has a narrow linewidth, the light source refraction module 212 employs a reflection grating to implement a wavelength adjustment function.
In this application, the effective measurement area is substantially perpendicular to the flow direction of the liquid to be measured.
As shown in fig. 3, which is a scenario applicable to the present application, after the liquid to be measured is ejected from the nozzle, an umbrella-shaped fog region is formed in a circumferential direction, which gradually spreads, and in the present application, the effective measurement area (the speed measurement section in the figure) is located below the nozzle and is perpendicular to the spraying direction.
When the planar light beams generated by the plurality of planar light beam forming devices 210 are shaped into parallel light beams by the lens group and then are converged in the effective measurement region 230, the effective measurement region 230 can measure the concentration of the sprayed liquid.
As shown in fig. 4, the sub-region 231 divides the effective measurement region 230 into several sub-regions 231, and each sub-region 231 (hereinafter referred to as voxel) has a certain thickness. When the effective measurement region 230 is divided, the division may be performed using a cartesian coordinate system or a polar coordinate system.
Step 120: and acquiring incident light intensity and emergent light intensity which are emitted by at least one light source and pass through the sub-region, determining the surface area of the liquid in the sub-region according to the incident light intensity and the emergent light intensity, and reconstructing the surface area concentration of the liquid in each voxel to form a first two-dimensional liquid surface concentration distribution.
Step 130: and acquiring incident light intensity and emergent light intensity which are emitted by at least one other light source and pass through the sub-region, determining the surface area of the liquid in the sub-region according to the incident light intensity and the emergent light intensity, and reconstructing the surface area concentration of the liquid in each voxel to form two-dimensional distribution of the surface concentration of the second liquid.
In order to realize the planar concentration distribution measurement of the liquid, the system of the present application is provided with a first liquid concentration distribution generating device and a second liquid concentration distribution generating device, wherein the first liquid concentration distribution generating device and the second liquid concentration distribution generating device are the same, and in actual use, the first liquid concentration distribution generating device and the second liquid concentration distribution generating device are all imaging devices 220 (such as a CCD camera, a COMS camera or an SCOMS camera, etc.), which generate the two-dimensional distribution of the liquid concentration by acquiring and recording the incident light intensity and the emergent light intensity of the planar light beam, and determining the attenuation coefficient according to the incident light intensity and the emergent light intensity.
During measurement, when the spray sprayed from the nozzle passes through an effective measurement area below the nozzle, the absorption of the light source in a voxel is completed to cause light intensity attenuation, and the light intensity attenuation meets the Beer's-Lambert law. I.e. the intensity of the light I emitted through the voxelλeAnd the incident light intensity Iλ0Satisfying an exponential decay. The attenuated laser light passes through the lens and is incident into the imaging device 220 for imaging.
Wherein the incident light intensity and the emergent light intensityAnd attenuation coefficient βλThe travel of the light in the voxel satisfies the following relation:
Iλe=Iλ0exp(-βλΔ) (1)
in the formula IλeTo output a set light intensity, Iλ0β for incident light intensityλΔ is the path taken by the light, which is the carbon dioxide absorption coefficient.
The calculation of the surface area concentration reconstruction is the same as the SetScan technology, and as in the left figure embodiment of FIG. 5, the region to be detected is divided into a plurality of voxels according to the number of pixels of the selected imaging device, and the parameters β in the voxelsλΔ is the unknown sought. According to the formula 1, the relation between the incident light intensity and the emergent light intensity of a certain voxel can be obtained to satisfy the formula 2 and the formula 3. The ratio of the emergent light intensity to the incident light intensity of a certain row/column of voxels satisfies formula 4, wherein n represents the nth voxel (n is related to the number of voxels divided by the detected region, i belongs to n).
An equation similar to that shown in equation 4 is satisfied for each image formed by the image forming apparatus 220. The results obtained for each imaging device 220
Figure BDA0002229873570000081
β corresponding to the voxel is obtained by iterative calculationλiΔiThe value is obtained. The iterative algorithm is similar to the ct (computed tomogry) algorithm. In the application, a MART (generalized adaptive reconstruction techniques) algorithm is adopted, and the method can be suitable for two-dimensional reconstruction calculation with large gradient change. For symmetric sprays, a group theory approach can be introduced to apply to iterative calculations in order to solve for a precise solution. For the above processes, details are not repeated herein, and related documents may be specifically referred to.
ΔiRelated to the refractive or relative refractive indices of the liquid and the ambient gas/liquid and the path taken by the light (i.e. the size of the voxel). The path traveled by light can be considered to be the same as a voxel, and the relative refractive index can be looked up from the properties of the fluid, thus Δiβ can be obtained according to the known quantityλiβ, as shown in equation 5λiRelated to the surface area of the liquid. From this, the liquid surface area (SetScan) within the voxel can be deducedTechnique, as shown on the right of fig. 5). And reconstructing the liquid surface concentration in each voxel to form a plane liquid surface concentration distribution.
Figure BDA0002229873570000082
Figure BDA0002229873570000083
Figure BDA0002229873570000084
βλΔ∝C (5)
Fig. 5 shows an example of a numerical simulation calculation result obtained in the case of using the effective measurement region division of the left diagram in the right diagram.
Finally, in step S140, the two-dimensional distribution of the surface concentration of the first liquid and the two-dimensional distribution of the surface concentration of the second liquid are correlated to form a two-dimensional distribution of the motion of the liquid.
The working principle of the transient two-dimensional velocity distribution measurement of the liquid motion in the application is based on the planar liquid surface concentration distribution obtained by the process. As shown in fig. 6, the light source emitting module 211 numbered 1 in each planar light beam forming device 210 simultaneously triggers the measured "planar liquid surface concentration distribution" result to form an image 1 (shown in fig. 7 a). After a time Δ T, the light source emitting module 211 numbered 2 in each planar light beam forming device 210 simultaneously triggers the measured "planar liquid surface concentration distribution" result to be image 2 (shown in fig. 7 b). The velocity measurement profile is calculated based on the cross-correlation of image 1 and image 2.
It should be noted that the exposure time of image 1 is substantially equal to the laser pulse width of number 1, and the exposure time of image 2 is greater than the laser pulse width of number 2, and is covered before and after the time.
The cross-correlation algorithm includes direct cross-correlation calculation and fourier transform. These two methods are now equivalent under less displacement conditions. The fourier Transform method generally uses a Fast Fourier Transform (FFT) method which is faster in computation speed.
The FFT method is briefly described as follows:
for image 1g1(x, y) and image 2g2(x, y) are respectively FFT transformed to obtain
Figure BDA0002229873570000091
Figure BDA0002229873570000092
In the formula
Figure BDA0002229873570000093
And
Figure BDA0002229873570000094
are respectively an image g1(x, y) and g2And (x, y) FFT transformation. x, y are image coordinates. And omega is a frequency-domain value corresponding to Fourier transform.
By using the translation characteristic of Fourier transform, the method can obtain
Figure BDA0002229873570000095
Wherein Δ x and Δ y are particle displacements.
Function(s)
Figure BDA0002229873570000096
By inverse Fourier transform
Figure BDA0002229873570000097
Substituting equation 8 into equation 9 results in
G(x,y)=g(x+△x,y+△y) (10)
Wherein g (x + △ x, y + y) is
Figure BDA0002229873570000101
G (x, y) has a maximum gray peak at (x, △ x, y + y), which is composed of 3 components (formula 11) due to the existence of background noise and other related quantities, the position of the maximum gray value (coherence peak) contains displacement information, so that the image displacement can be obtained by extracting the displacement of the peak center.
R(S)=RC(S)+RD(S)+RF(S) (11)
In the formula RDAnd (S) represents the maximum gray value and represents the displacement information. RC(S)+RFAnd (S) is random correlation and background noise correlation quantity.
By the above processing, the relative displacement of the density distributions of the image 1 and the image 2 can be obtained. From this, the moving speed of the concentration distribution, i.e., the moving speed of the liquid can be calculated. The calculation formula is shown as formula (12):
Figure BDA0002229873570000102
in the formula, Δ s is a displacement (vector) of the liquid, and Δ t is a time interval for recording images of "planar liquid surface concentration distribution" at different times. When delta t → 0, the instantaneous velocity vector v of the liquid is obtained.
As shown in fig. 8, the vector diagram obtained by the cross-correlation processing of the images of the "planar liquid surface concentration distribution" in fig. 7a and 7b is obtained, and Matlab is used as the processing software.
The method and the system for simultaneously measuring the liquid concentration and the liquid motion two-dimensional distribution can realize the two-dimensional speed measurement of the liquid phase of a certain section below the spray, can also display the liquid concentration change trend in real time for real transient measurement, expand the flow measurement function, and finally realize the simultaneous measurement of the concentration and the speed of the liquid.
The above description is only for the specific embodiments of the present application, but the scope of the present application is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present application should be covered within the scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.

Claims (10)

1. A method for simultaneously measuring two-dimensional distribution of liquid concentration and liquid motion, which is characterized by comprising
Generating at least two parallel plane light beams, and converging the plane light beams to form an effective measuring area, wherein the effective measuring area is vertical to the flow direction of liquid, the effective measuring area comprises a plurality of sub-areas, and each plane light beam at least comprises two light sources which are sent according to a preset time interval and are converged;
acquiring incident light intensity and emergent light intensity which are emitted by at least one light source and pass through the sub-region, determining the liquid surface area in the sub-region according to the incident light intensity and the emergent light intensity, and reconstructing the liquid surface area concentration in each voxel to form first liquid surface concentration two-dimensional distribution;
acquiring incident light intensity and emergent light intensity which are emitted by at least one other light source and pass through the sub-region, determining the liquid surface area in the sub-region according to the incident light intensity and the emergent light intensity, and reconstructing the liquid surface area concentration in each voxel to form second liquid surface concentration two-dimensional distribution;
and mutually correlating the first liquid surface concentration two-dimensional distribution and the second liquid surface concentration two-dimensional distribution to form liquid motion two-dimensional distribution.
2. The method of claim 1, wherein each planar beam comprises at least two light sources, at least one of which is refracted and converged with another light source before being scattered to form a planar beam.
3. The method of simultaneously measuring two-dimensional distributions of liquid concentration and liquid motion of claim 1, wherein said at least two light sources are emitted at predetermined time intervals by coordinated control.
4. A method for simultaneously measuring two-dimensional distributions of liquid concentration and liquid motion according to any one of claims 1 to 3, wherein said light source is a pulsed laser.
5. The method of claim 1, wherein the incident light intensity and the emergent light intensity determine the surface area of the liquid in the sub-region, and are determined by the following relationship:
Iλe=Iλ0exp(-βλΔ)
in the formula IλeTo output a set light intensity, Iλ0β for incident light intensityλTo the attenuation coefficient, and in relation to the surface area of the liquid, Δ is the path taken by the light.
6. A system for simultaneously measuring two-dimensional distributions of liquid concentration and liquid motion, the system comprising
The device comprises at least two plane light beam forming devices, a liquid level sensor and a control device, wherein the at least two plane light beam forming devices are used for generating at least two parallel plane light beams and converging the plane light beams to form an effective measuring area, the effective measuring area is vertical to the flow direction of liquid, the effective measuring area comprises a plurality of sub-areas, and each plane light beam at least comprises two light sources which are sent at preset time intervals and converged;
the first liquid concentration distribution generating device is used for acquiring incident light intensity and emergent light intensity which are emitted by at least one light source and pass through the sub-area, determining the liquid surface area in the sub-area according to the incident light intensity and the emergent light intensity, and reconstructing the liquid surface area concentration in each voxel to form first liquid surface concentration two-dimensional distribution;
the second liquid concentration distribution generating device is used for acquiring incident light intensity and emergent light intensity which are emitted by at least one other light source and pass through the sub-area, determining the liquid surface area in the sub-area according to the incident light intensity and the emergent light intensity, and reconstructing the liquid surface area concentration in each voxel to form second liquid surface concentration two-dimensional distribution;
and the liquid motion distribution generating device is used for mutually correlating the first liquid surface concentration two-dimensional distribution and the second liquid surface concentration two-dimensional distribution to form a liquid motion two-dimensional distribution.
7. The system for simultaneously measuring two-dimensional distributions of liquid concentration and liquid motion according to claim 6, wherein said planar beam forming means comprises
At least two light source emitting modules for generating and emitting light sources;
at least one light source refraction module for reflecting or refracting the at least one light source at a predetermined angle to converge the light sources emitted by the at least two light source emission modules; and
the sheet light adjusting module is used for scattering the converged light source to form a plane light beam.
8. The system for simultaneously measuring two-dimensional distributions of liquid concentration and liquid motion according to claim 6, wherein said planar beam forming means further comprises
And the control module is used for coordinating and controlling at least two light source emitting modules to emit light sources according to a preset time interval.
9. The system for simultaneously measuring the two-dimensional distribution of the concentration and the movement of the liquid as claimed in any one of claims 6 to 8, wherein the light source generated by the light source emitting module is a pulsed laser.
10. The system for simultaneously measuring two-dimensional distributions of liquid concentration and liquid motion according to claim 6, wherein said first liquid concentration generating means and/or said second liquid concentration generating means determines the liquid surface area in said sub-area based on said incident light intensity and said emergent light intensity by the following relationship:
Iλe=Iλ0exp(-βλΔ)
in the formula IλeTo output a set light intensity, Iλ0β for incident light intensityλTo the attenuation coefficient, and in relation to the surface area of the liquid, Δ is the path taken by the light.
CN201910963960.0A 2019-10-11 2019-10-11 Method and system for simultaneously measuring two-dimensional distribution of liquid concentration and liquid movement Pending CN110879301A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910963960.0A CN110879301A (en) 2019-10-11 2019-10-11 Method and system for simultaneously measuring two-dimensional distribution of liquid concentration and liquid movement

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910963960.0A CN110879301A (en) 2019-10-11 2019-10-11 Method and system for simultaneously measuring two-dimensional distribution of liquid concentration and liquid movement

Publications (1)

Publication Number Publication Date
CN110879301A true CN110879301A (en) 2020-03-13

Family

ID=69727981

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910963960.0A Pending CN110879301A (en) 2019-10-11 2019-10-11 Method and system for simultaneously measuring two-dimensional distribution of liquid concentration and liquid movement

Country Status (1)

Country Link
CN (1) CN110879301A (en)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1654962A (en) * 2005-01-18 2005-08-17 浙江大学 Two-phase flow digital particle image speed measurement method and device
CN101082561A (en) * 2007-07-03 2007-12-05 浙江大学 Method for measuring solid particle three-dimensional concentration field and velocity field in gas/solid two-phase stream
CN102156112A (en) * 2011-03-08 2011-08-17 中国科学院安徽光学精密机械研究所 Flue gas flow velocity measuring device and method
DE102010033027A1 (en) * 2010-07-29 2012-02-02 Iav Gmbh Ingenieurgesellschaft Auto Und Verkehr Method for measuring concentration of fuel in combustion chamber of internal combustion engine, involves vaporizing fuel so that phase boundary surface is formed between solid or liquid particles in fuel and fuel mixture
CN103472256A (en) * 2013-09-25 2013-12-25 东南大学 Flow two-dimensional velocity profile measuring method and device based on planar array CCD spatial filter
CN103889488A (en) * 2011-10-28 2014-06-25 皇家飞利浦有限公司 Analysis and control of aerosol flow
CN104019763A (en) * 2014-06-18 2014-09-03 重庆大学 Synchronous measurement device of fluid three-dimensional velocity field and deformable body three-dimensional shape
CN106018280A (en) * 2016-07-19 2016-10-12 江苏大学 Device and method capable of simultaneously measuring velocity field and concentration field
CN106290093A (en) * 2015-06-29 2017-01-04 天津同阳科技发展有限公司 Particle concentration based on mie scattering theory monitoring system and using method thereof
CN108020168A (en) * 2017-11-23 2018-05-11 哈尔滨工程大学 Nearly free surface gas-liquid two-phase flow field three-dimension measuring system and measuring method based on particle image velocimetry
CN208187912U (en) * 2018-01-25 2018-12-04 福州大学 The fibre-optical probe measuring instrument of particle speed and concentration in a kind of multiphase flow
CN109061220A (en) * 2018-09-04 2018-12-21 北京航空航天大学 A kind of air-flow two-dimension speed distribution measurement method based on laser absorption spectrum chromatography imaging technique

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1654962A (en) * 2005-01-18 2005-08-17 浙江大学 Two-phase flow digital particle image speed measurement method and device
CN101082561A (en) * 2007-07-03 2007-12-05 浙江大学 Method for measuring solid particle three-dimensional concentration field and velocity field in gas/solid two-phase stream
DE102010033027A1 (en) * 2010-07-29 2012-02-02 Iav Gmbh Ingenieurgesellschaft Auto Und Verkehr Method for measuring concentration of fuel in combustion chamber of internal combustion engine, involves vaporizing fuel so that phase boundary surface is formed between solid or liquid particles in fuel and fuel mixture
CN102156112A (en) * 2011-03-08 2011-08-17 中国科学院安徽光学精密机械研究所 Flue gas flow velocity measuring device and method
CN103889488A (en) * 2011-10-28 2014-06-25 皇家飞利浦有限公司 Analysis and control of aerosol flow
CN103472256A (en) * 2013-09-25 2013-12-25 东南大学 Flow two-dimensional velocity profile measuring method and device based on planar array CCD spatial filter
CN104019763A (en) * 2014-06-18 2014-09-03 重庆大学 Synchronous measurement device of fluid three-dimensional velocity field and deformable body three-dimensional shape
CN106290093A (en) * 2015-06-29 2017-01-04 天津同阳科技发展有限公司 Particle concentration based on mie scattering theory monitoring system and using method thereof
CN106018280A (en) * 2016-07-19 2016-10-12 江苏大学 Device and method capable of simultaneously measuring velocity field and concentration field
CN108020168A (en) * 2017-11-23 2018-05-11 哈尔滨工程大学 Nearly free surface gas-liquid two-phase flow field three-dimension measuring system and measuring method based on particle image velocimetry
CN208187912U (en) * 2018-01-25 2018-12-04 福州大学 The fibre-optical probe measuring instrument of particle speed and concentration in a kind of multiphase flow
CN109061220A (en) * 2018-09-04 2018-12-21 北京航空航天大学 A kind of air-flow two-dimension speed distribution measurement method based on laser absorption spectrum chromatography imaging technique

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
LIU C, ET AL: "Development of a fan-beam TDLAS-based tomographic sensor for rapid imaging of temperature and gas concentration", 《OPTICS EXPRESS》 *
刘存喜: "多级旋流空气雾化喷嘴雾化特性及光学测试方法研究", 《中国博士学位论文全文数据库》 *
张宝华 等: "运用相移多谱勒激光诊断技术研究燃烧室中的两相流动", 《航空发动机》 *

Similar Documents

Publication Publication Date Title
CN102323445B (en) Method for measuring flowing velocity of recycle gas in multiple laser cavity
Klinner et al. Tomographic shadowgraphy for three-dimensional reconstruction of instantaneous spray distributions
CN105043946A (en) Dual-wavelength based scattering angle self-calibration full-field rainbow measurement method and device
US20220043153A1 (en) Light Detection and Ranging
Manovski et al. 3D Lagrangian particle tracking of a subsonic jet using multi-pulse Shake-The-Box
EP3862787A1 (en) De-jitter of point cloud data for target recognition
US11624693B2 (en) System for estimating planar drop sizes
Eberhart et al. 3-D visualization of transparent fluid flows from snapshot light field data
CN110879301A (en) Method and system for simultaneously measuring two-dimensional distribution of liquid concentration and liquid movement
CN110879300A (en) Method and system for measuring velocity of flowing particles
Mayinger Advanced experimental methods
Dudorov et al. Retrieval of crosswind velocity based on the analysis of remote object images: Part 1—drift of a thin layer of turbulent inhomogeneities
Tislenko et al. Improving the Distance Measurement Accuracy of Active-Pulse Television Measuring Systems Using Polynomial Approximation
Mayinger The 1991 Max Jakob Memorial Award Lecture: Image-forming optical techniques in heat transfer: Revival by computer-aided data processing
Bencic et al. Advanced Optical Diagnostics for Ice Crystal Cloud Measurements in the NASA Glenn Propulsion Systems Laboratory
Savelyev et al. Stereo thermal marking velocimetry
Anders et al. Simultaneous in situ measurements of size and velocity of burning droplets
Ota et al. Laser interferometric CT measurement of the unsteady supersonic shock-vortex flow field discharging from two parallel and cylindrical nozzles
CN110726691A (en) Method and system for measuring two-dimensional distribution of carbon dioxide concentration
Wang et al. Sequential observation of rebound shock wave generated by collapse of vapor bubble in BOS system
Desantes et al. Characterization of local fuel concentrations and internal dynamics in DI Diesel sprays
JP2005308439A (en) Three dimensional geometry measuring arrangement by pattern projection method
Davis Tomographic plenoptic background oriented schlieren imaging
Komander et al. Data fusion of surface normals and point coordinates for deflectometric measurements
Atkinson et al. Evaluation and application of tomographic background oriented schlieren and PIV measurements to a heated jet

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20200313