EP3812712B1 - Fluid flow analysis - Google Patents

Fluid flow analysis Download PDF

Info

Publication number
EP3812712B1
EP3812712B1 EP19204453.5A EP19204453A EP3812712B1 EP 3812712 B1 EP3812712 B1 EP 3812712B1 EP 19204453 A EP19204453 A EP 19204453A EP 3812712 B1 EP3812712 B1 EP 3812712B1
Authority
EP
European Patent Office
Prior art keywords
measurement region
fluid
monitoring
fluid flow
variations
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
EP19204453.5A
Other languages
German (de)
French (fr)
Other versions
EP3812712C0 (en
EP3812712A1 (en
Inventor
Max Thele
Christian Mundt
Tobias Sander
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.)
Universitaet der Bundeswehr Muenchen
Original Assignee
Universitaet der Bundeswehr Muenchen
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 Universitaet der Bundeswehr Muenchen filed Critical Universitaet der Bundeswehr Muenchen
Priority to EP19204453.5A priority Critical patent/EP3812712B1/en
Publication of EP3812712A1 publication Critical patent/EP3812712A1/en
Application granted granted Critical
Publication of EP3812712B1 publication Critical patent/EP3812712B1/en
Publication of EP3812712C0 publication Critical patent/EP3812712C0/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/704Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow using marked regions or existing inhomogeneities within the fluid stream, e.g. statistically occurring variations in a fluid parameter
    • G01F1/708Measuring the time taken to traverse a fixed distance
    • G01F1/7086Measuring the time taken to traverse a fixed distance using optical detecting arrangements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/704Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow using marked regions or existing inhomogeneities within the fluid stream, e.g. statistically occurring variations in a fluid parameter
    • G01F1/708Measuring the time taken to traverse a fixed distance
    • G01F1/712Measuring the time taken to traverse a fixed distance using auto-correlation or cross-correlation detection means
    • 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
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/66Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by measuring frequency, phase shift or propagation time of electromagnetic or other waves, e.g. using ultrasonic flowmeters
    • G01F1/661Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by measuring frequency, phase shift or propagation time of electromagnetic or other waves, e.g. using ultrasonic flowmeters using light
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/41Refractivity; Phase-affecting properties, e.g. optical path length
    • G01N2021/416Visualising flow by index measurement
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/41Refractivity; Phase-affecting properties, e.g. optical path length
    • G01N21/45Refractivity; Phase-affecting properties, e.g. optical path length using interferometric methods; using Schlieren methods

Definitions

  • the present disclosure relates to a sensor arrangement for measuring properties of a fluid flow. More particularly, the present invention relates to a nonintrusive analysis of a fluid flow based on differential interferometry.
  • US 2004/113081 A1 relates to measuring the velocity of a flowing gas by detecting the pattern of infrared absorption of the gas stream at two locations spaced apart in the direction of flow and measuring the time lapse between the sensing of the same radiation pattern by the first and second sensors.
  • the velocity can be derived from the time lapse and the known distance between the sensors.
  • US 2019/0293543 A1 relates to measuring the properties of a liquid by exploiting the power modulation a laser light beam undergoes due to the retroreflection of the laser light beam itself towards the laser cavity from which the laser is generated when this laser light is directed towards a transparent conduit through which the liquid for which the properties are to be measured flows.
  • WO 01/90700 A2 relates to monitoring the flow rate of a fluid stream by heating the stream with a heating member to induce a change in the refractive index of the fluid, detecting the change in the refractive index of the fluid at a location remote from the heating member, and calculating the flow rate of the fluid from the change in the refractive index.
  • the present invention provides a method for analyzing a fluid flow according to claim 1 and a fluid flow analyzer according to claim 7.
  • the fluid flow may be analyzed non-intrusively.
  • the formulation "determining whether a first density variation, which is detected by monitoring the first measurement region, and a second density variation, which is detected by monitoring the second measurement region, match" particularly refers to determining whether (or when) a fluid volume that exhibits a (specific) non-uniform density distribution (or a specific density pattern) passes through the first measurement region and the second measurement region.
  • the method may further comprise determining a speed of the fluid flow between the first measurement region and the second measurement region based on a time period that passes between a first instant at which the first density variation is detected in the first measurement region and a second instant at which the second density variation is detected in the second measurement region and the distance between the centers of the measurement regions.
  • the time period may be determined based on calculating a temporal correlation between first density values detected (or measured) within the first measurement region and second density values detected (or measured) within the second measurement region.
  • the distance between the centers of the first measurement region and the second measurement region may be equal to or less than 100 times, 50 times, 20 times, 10 times or one time a size of the first or second measurement region (along an axis through the centers of the first measurement region and the second measurement region).
  • the first density variation and the second density variation may be caused, at least in part, by turbulence.
  • a turbulent flow which is to be analyzed may exhibit a non-uniform temperature distribution (without intrusive heating).
  • intra heating particularly refers to pulsed heating of a location upstream of the first measurement region for (no other purpose than) causing (artificial) temperature variations.
  • a density variation caused by a temperature variation travels at the speed of the fluid flow, as opposed to sound waves which travel at the speed of sound through the fluid.
  • monitoring the measurement regions for a density variation within the fluid which is caused by a temperature variation allows determining how long it took the fluid volume carrying the temperature variation to travel from the first measurement region to the second measurement region.
  • Monitoring the first measurement region for variations of the refractive index of the fluid is based on measuring an intensity of interfering electromagnetic waves which travel along different paths through the first measurement region
  • monitoring the second measurement region for variations of the refractive index of the fluid is based on measuring an intensity of interfering electromagnetic waves which travel along different paths through the second measurement region.
  • the distance between the centers of the first measurement region and the second measurement region may be equal to or less than 100 times, 50 times, 20 times, 10 times or one time a distance between focus points of the electromagnetic waves within the first measurement region, or a mean distance between the different paths through the first measurement region.
  • a shorter distance may facilitate matching, because the temperature gradient may decrease over time as a result of heat conduction from the (relatively) warmer fluid volume to the (relatively) cooler fluid volume.
  • Determining whether a first density variation, which is detected by monitoring the first measurement region, and a second density variation, which is detected by monitoring the second measurement region, match comprises calculating a temporal correlation between the intensity of the interfering electromagnetic waves which travel along the different paths through the first measurement region and the intensity of the interfering electromagnetic waves which travel along the different paths through the second measurement region.
  • the method may further comprise determining a speed of sound within the fluid flow from a width of a peak of a function of the calculated temporal correlation.
  • the method may further comprise monitoring a third measurement region for density variations within the fluid, as the fluid passes the third measurement region, monitoring a fourth measurement region for density variations within the fluid, as the fluid passes the fourth measurement region, and determining whether a third density variation, which is detected by monitoring the third measurement region, and a fourth density variation, which is detected by monitoring the fourth measurement region, match, wherein a first axis passes through the center of the first measurement region and the center of the second measurement region, a second axis passes through the center of the third measurement region and the center of the fourth measurement region, and the first axis and the second axis are not parallel and preferably perpendicular to each other to allow for a two-dimensional analysis of the fluid flow.
  • the method may further comprise monitoring a fifth measurement region for density variations within the fluid, as the fluid passes the fifth measurement region, monitoring a sixth measurement region for density variations within the fluid, as the fluid passes the sixth measurement region, and determining whether a fifth density variation, which is detected by monitoring the fifth measurement region, and a sixth density variation, which is detected by monitoring the sixth measurement region, match, wherein a third axis passes through the center of the fifth measurement region and the center of the sixth measurement region, the first axis and the third axis are not parallel and preferably perpendicular to each other, and the second axis and the third axis are not parallel and preferably perpendicular to each other to allow for a three-dimensional analysis of the fluid flow.
  • the computing device of the fluid flow analyzer may be configured to neglect intensity variations and/or local maxima caused by sound waves.
  • the emitter of the fluid flow analyzer may be a laser and the beams (of electromagnetic radiation) may be laser beams.
  • the fluid flow analyzer may comprise optical components.
  • the optical components may be configured to split a beam emitted from the emitter and/or focus the resulting beams to a spot in the first and second measurement region, respectively.
  • the spots may lie on a line of flow to increase the accuracy when determining the speed of the fluid flow.
  • the computing device may be configured to neglect intensity variations and/or local maxima caused by sound waves by filtering intensity variations that are below a threshold.
  • the computing device may be configured to determine a speed of the fluid flow between the first measurement region and the second measurement region based on the temporal shift and a distance between the centers of the first measurement region and the second measurement region.
  • the distance between the centers of the first measurement region and the second measurement region may be equal to or less than 100 times, 50 times, 20 times, 10 times or one time a size of the first measurement region or the second measurement region (along an axis through the centers of the first measurement region and the second measurement regions).
  • the first, second, third and fourth beam may be derived by splitting a beam emitted by the emitter.
  • the beam may be split by two Wollaston prisms.
  • the fluid flow analyzer may further comprise a third measurement region delimited by a fifth optical path of a fifth beam (of electromagnetic radiation) and a sixth optical path of a sixth beam (of electromagnetic radiation), a fourth measurement region delimited by a seventh optical path of a seventh beam (of electromagnetic radiation) and an eighth optical path of an eighth beam (of electromagnetic radiation), a third sensor for detecting a third intensity of the superimposed fifth and sixth beams, and a fourth sensor for detecting a fourth intensity of the superimposed seventh and eighth beams, wherein the computing device is configured to determine a temporal shift at which a correlation between the fourth intensity and the fifth intensity has a local maximum.
  • the first optical path and the fifth optical path may be not parallel and preferably perpendicular to each other to allow for a two-dimensional analysis of the fluid flow.
  • Fig. 1 shows a fluid flow analyzer which does not form part of the invention.
  • the fluid flow analyzer 10 comprises an emitter 12a (e.g., a laser) which emits a beam 14a of coherent electromagnetic radiation (e.g., a laser beam).
  • the beam 14a is split into a first beam 16a and a second beam 16b by optical components 18a arranged within the optical path of the beam 14a.
  • the first beam 16a and the second beam 16b traverse a fluid flow 20 and generate an interference pattern that is incident on a first sensor 22a.
  • the intensities detected by the first sensor 22a depend on the density of the fluid. This allows detecting a density variation within a fluid volume that crosses the beams 16a, 16b.
  • the bending angles of the first and second beams 16a, 16b depend on the refractive indices of the fluid and its surroundings and thus on the density of the fluid. Because the fluid flow 20 exhibits (relatively) cooler parts 20a and (relatively) warmer parts 20b, the bending angles of the first and second beams 16a, 16b change over time when the fluid passes through the measurement region 24a (schematically illustrated by the broken rectangle). As the first beam 16a and the second beam 16b are displaced from each other in the direction of the fluid flow 20 (indicated by the arrows), there is a delay between the bending pattern exhibited by the first beam 16a and the bending pattern exhibited by the second beam 16b.
  • a further measurement region 24b is established downstream by replicating the structure described above.
  • another emitter 12b e.g., another laser
  • a beam 14b of coherent electromagnetic radiation e.g., a laser beam
  • the third beam 16c and the fourth beam 16d traverse the fluid flow 20 and generate an interference pattern that is detected by a second sensor 22b.
  • the bending angles of the third and fourth beams 16c, 16d change over time when the fluid passes through the measurement region 24b (schematically illustrated by the broken rectangle).
  • the third beam 16c and the fourth beam 16d are displaced from each other in the direction of the fluid flow (indicated by the arrows), there is a delay between the bending pattern exhibited by the third beam 16c and the fourth beam 16d.
  • the fluid flow analyzer 10 further comprises a computing device 26 for determining a temporal shift ⁇ t at which a correlation C( ⁇ t) between the intensities detected by the sensors 22a, 22b has a local maximum.
  • the computing device 26 strives at finding matching density variations which allows calculating a speed of the fluid flow by dividing a distance between the measurement regions 24a, 24b by the temporal shift ⁇ t.
  • the computing device 26 may be configured to neglect intensity variations and/or local maxima caused by sound waves. As the density variations caused by sound waves are often smaller in amplitude than the density variations caused by temperature variations, the computing device 26 may be configured to filter intensity variations that are below a threshold. Moreover, sound waves caused by an isotropic turbulent flow may rather widen a correlation peak than adversely affect the calculation of the speed of the fluid.
  • Fig. 2 illustrates a modified fluid flow analyzer 10 which is based on two focused laser differential interferometers, FLDIs and differs from the fluid flow analyzer 10 of Fig. 1 in that the dimensions of the measurement regions 24a, 24b (perpendicular to the direction of flow) are reduced to facilitate matching and to improve spatial resolution (perpendicular to the direction of flow) with regard to determining the speed of the flow by dividing the distance D of the measurement regions 24a, 24b by the temporal shift ⁇ t.
  • the optical components 18a arranged within the optical path of the beam 14a are replaced by a lens 28a with a short focal length, a quarter wave plate 30a, a Wollaston prism 32a, and a lens 34a that focusses the first and second beams 16a, 16b to different spots within the fluid flow 20.
  • the first and second beam 16a, 16b are guided onto the sensor 22a by an inverse optical setting comprising a lens 36a, a Wollaston prism 38a, a quarter wave plate 40a, and a lens 42a with a short focal length.
  • a further measurement region 24b is established downstream by replicating said structure.
  • the optical components 18b arranged within the optical path of the beam 14b are replaced by a lens 28b with a short focal length, a quarter wave plate 30b, a Wollaston prism 32b, and a lens 34b that focusses the third and fourth beams 16c, 16d to different spots within the fluid flow 20.
  • the third and fourth beam 16c, 16d are guided onto the sensor 22b by an inverse setting comprising a lens 36b, a Wollaston prism 38b, a quarter wave plate 40b, and a lens 42b with a short focal length.
  • the intensity values detected by the sensors 22a, 22b are provided to the computing device 26 which is configured to determine a temporal shift ⁇ t at which a correlation C( ⁇ t) between the intensities detected by the sensors 22a, 22b has a local maximum.
  • the temporal shift indicates the time required for a fluid volume carrying a density variation to travel from the upstream measurement region 24a to the downstream measurement region 24b and allows calculating the speed of the fluid flow between the measurement regions 24a, 24b by dividing the distance D between the (centers of the) measurement regions 24a, 24b by the temporal shift ⁇ t.
  • the Distance (D) between the centers of the measurement regions (24a, 24b) is 0.01 m or less.
  • Fig. 3 and Fig. 3a illustrate a modified fluid flow analyzer 10 (that may be used for analyzing supersonic or hypersonic fluid flows 20).
  • the modified fluid flow analyzer 10 differs from the fluid flow analyzer 10 of Fig. 2 in that the distance D between the measurement regions 24a, 24b (in the direction of flow) is reduced to facilitate matching and increase spatial resolution (in the direction of flow) with regard to measuring the speed of the flow.
  • the beam 14 emitted by the emitter 12 e.g., a laser
  • the beam 14 is expanded by a lens 28 with a short focal length and the expanded beam 14 is circularly polarized by a first quarter wave plate 30 that has an optical axis which is rotated by 45° relative to the direction of polarization of the incident beam 14.
  • the circularly polarized expanded beam 14 is then split into two beams 14a, 14b by a first Wollaston prism 32 which may bend the beams 14a, 14b by less than 5 degrees (e.g., two degrees).
  • Beam 14a is circularly polarized by a second quarter wave plate 44 and split into a first beam 16a and a second beam 16b by a second Wollaston prism 46, wherein the direction of polarization of the first beam 16a may be perpendicular to the direction of polarization of the second beam 16b.
  • Beam 14b is circularly polarized by the second quarter wave plate 44 and split into a third beam 16c and a fourth beam 16d by the second Wollaston prism 46, wherein the direction of polarization of the third beam 16c may be perpendicular to the direction of polarization of the fourth beam 16d.
  • the second Wollaston prism 46 may bend the beams 16a, 16b, 16c, 16d by less than 2.5 degrees (e.g., one degree).
  • the beams 16a, 16b, 16c, 16d may be focused to spots within the fluid flow 20.
  • the spots may be placed on a line of flow. If a density variation travels along the line of flow, it will be detected subsequently by both sensors 22a, 22b.
  • the first and second beams 16a, 16b After traversing the first measurement region 24a, the first and second beams 16a, 16b incide on a third Wollaston prism 38a and a polarizer 48a and cause an interference pattern on the first sensor 22a (e.g., a photodiode which outputs a signal that is proportional to the intensity of the incident electromagnetic radiation).
  • the third and fourth beams 16c, 16d After traversing the second measurement region 24b, the third and fourth beams 16c, 16d incide on a fourth Wollaston prism 38b and a polarizer 48b and cause an interference pattern on the second sensor 22b (e.g., a photodiode which outputs a signal that is proportional to the intensity of the incident electromagnetic radiation).
  • a distance D between the (centers of the) measurement regions 24a, 24b may be equal to a distance between the focal points of the first beam 16a and the second beam 16b or to a distance of the focal points of the third beam 16c and the fourth beam 16d.
  • the distance D between the centers of the measurement regions (24a, 24b) is or less.
  • the speed of the fluid may be measured in two or three dimensions.
  • Fig. 3b shows graphs from two experiments that are based on the setting described with reference to Fig. 3 and Fig. 3a .
  • the top left graph of Fig. 3b shows a first signal S1 output by the first sensor 22a, a second signal S2 output by the second sensor 22b and a signal H indicating the onset of the fluid flow 20 (from the moment of signal rise).
  • the top right graph of Fig. 3b shows a first signal S1 output by the first sensor 22a and a second signal S2 output by the second sensor 22b.
  • the middle and bottom graphs show the local maxima (peaks) of a correlation function of the first signal S1 and the second signal S2.
  • the analyzer 10 is robust against mechanical shocks, as the frequency of such shocks is typically much smaller than the frequency of the signals S1, S2 which allows filtering mechanical shocks without degrading the accuracy of the analyzer 10.
  • Fig. 4 shows a flow chart of a method for analyzing a fluid flow that may be practiced based on any one of the analyzers 10 of Fig. 1 to Fig. 3 .
  • the method starts at step 50 with monitoring the first measurement region 24a for density variations within the fluid, as the fluid passes the first measurement region 24a. For instance, density variations of the fluid passing through the first measurement region 24a may be detected by interferometry.
  • the method continues at step 52 with monitoring the second measurement region 24b for density variations within the fluid, as the fluid passes the second measurement region 24b. For instance, density variations of the fluid passing through the second measurement region 24b may be detected based on interferometry.
  • the method ends with determining whether a first density variation, which is detected by monitoring the first measurement region 24a, and a second density variation, which is detected by monitoring the second measurement region 24b, match.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Measuring Volume Flow (AREA)

Description

    Technical Field
  • The present disclosure relates to a sensor arrangement for measuring properties of a fluid flow. More particularly, the present invention relates to a nonintrusive analysis of a fluid flow based on differential interferometry.
  • Background Art
  • Differential interferometry allows detecting density gradients in a fluid flow (cf. Smeets, Guenter and Avalos George. "Laser-Differential Interferometer Applications in Gas Dynamics." (1996)).
  • Moreover, US 2004/113081 A1 relates to measuring the velocity of a flowing gas by detecting the pattern of infrared absorption of the gas stream at two locations spaced apart in the direction of flow and measuring the time lapse between the sensing of the same radiation pattern by the first and second sensors. The velocity can be derived from the time lapse and the known distance between the sensors.
  • US 2019/0293543 A1 relates to measuring the properties of a liquid by exploiting the power modulation a laser light beam undergoes due to the retroreflection of the laser light beam itself towards the laser cavity from which the laser is generated when this laser light is directed towards a transparent conduit through which the liquid for which the properties are to be measured flows.
  • WO 01/90700 A2 relates to monitoring the flow rate of a fluid stream by heating the stream with a heating member to induce a change in the refractive index of the fluid, detecting the change in the refractive index of the fluid at a location remote from the heating member, and calculating the flow rate of the fluid from the change in the refractive index.
  • Further background can be found in WO 91/18280 A1 , US 2007/044572 A1 , US 2011/106397 A1 , EP1205736 A2 , EP0952431 A2 , DE3909621 A1 , US 2014/130606 A1 , and GB2352033 A .
  • Summary of Invention
  • The present invention provides a method for analyzing a fluid flow according to claim 1 and a fluid flow analyzer according to claim 7.
  • By monitoring the measurement regions for variations of the refractive index of the fluid, the fluid flow may be analyzed non-intrusively.
  • In this regard, the term "measurement region", as used throughout the description and the claims, particularly refers to a region around a fixed location through which (at least a part of) the fluid flows. Moreover, a measurement region may be a region (within the fluid flow) which is delimited by parallel or at least partially overlapping beams of electromagnetic radiation. Moreover, the formulation "monitoring a first measurement region for density variations within the fluid", as used throughout the description and the claims, particularly refers to determining whether a fluid volume that passes through the measurement region exhibits a non-uniform density distribution (or a specific density pattern). Hence, the formulation "determining whether a first density variation, which is detected by monitoring the first measurement region, and a second density variation, which is detected by monitoring the second measurement region, match", as used throughout the description and the claims, particularly refers to determining whether (or when) a fluid volume that exhibits a (specific) non-uniform density distribution (or a specific density pattern) passes through the first measurement region and the second measurement region.
  • The method may further comprise determining a speed of the fluid flow between the first measurement region and the second measurement region based on a time period that passes between a first instant at which the first density variation is detected in the first measurement region and a second instant at which the second density variation is detected in the second measurement region and the distance between the centers of the measurement regions.
  • For example, the time period may be determined based on calculating a temporal correlation between first density values detected (or measured) within the first measurement region and second density values detected (or measured) within the second measurement region.
  • The distance between the centers of the first measurement region and the second measurement region may be equal to or less than 100 times, 50 times, 20 times, 10 times or one time a size of the first or second measurement region (along an axis through the centers of the first measurement region and the second measurement region).
  • The first density variation and the second density variation may be caused, at least in part, by turbulence.
  • For example, a turbulent flow which is to be analyzed may exhibit a non-uniform temperature distribution (without intrusive heating).
  • In this regard, the term "intrusive heating" particularly refers to pulsed heating of a location upstream of the first measurement region for (no other purpose than) causing (artificial) temperature variations.
  • A density variation caused by a temperature variation travels at the speed of the fluid flow, as opposed to sound waves which travel at the speed of sound through the fluid. Hence, monitoring the measurement regions for a density variation within the fluid which is caused by a temperature variation allows determining how long it took the fluid volume carrying the temperature variation to travel from the first measurement region to the second measurement region.
  • Moreover, if the turbulence is isotropic, sound waves may travel in all directions and hence, the sound waves may not substantially affect (degrade) the monitoring.
  • Monitoring the first measurement region for variations of the refractive index of the fluid is based on measuring an intensity of interfering electromagnetic waves which travel along different paths through the first measurement region, and monitoring the second measurement region for variations of the refractive index of the fluid is based on measuring an intensity of interfering electromagnetic waves which travel along different paths through the second measurement region.
  • The distance between the centers of the first measurement region and the second measurement region may be equal to or less than 100 times, 50 times, 20 times, 10 times or one time a distance between focus points of the electromagnetic waves within the first measurement region, or a mean distance between the different paths through the first measurement region.
  • A shorter distance may facilitate matching, because the temperature gradient may decrease over time as a result of heat conduction from the (relatively) warmer fluid volume to the (relatively) cooler fluid volume.
  • Determining whether a first density variation, which is detected by monitoring the first measurement region, and a second density variation, which is detected by monitoring the second measurement region, match comprises calculating a temporal correlation between the intensity of the interfering electromagnetic waves which travel along the different paths through the first measurement region and the intensity of the interfering electromagnetic waves which travel along the different paths through the second measurement region.
  • The method may further comprise determining a speed of sound within the fluid flow from a width of a peak of a function of the calculated temporal correlation.
  • The method may further comprise monitoring a third measurement region for density variations within the fluid, as the fluid passes the third measurement region, monitoring a fourth measurement region for density variations within the fluid, as the fluid passes the fourth measurement region, and determining whether a third density variation, which is detected by monitoring the third measurement region, and a fourth density variation, which is detected by monitoring the fourth measurement region, match, wherein a first axis passes through the center of the first measurement region and the center of the second measurement region, a second axis passes through the center of the third measurement region and the center of the fourth measurement region, and the first axis and the second axis are not parallel and preferably perpendicular to each other to allow for a two-dimensional analysis of the fluid flow.
  • The method may further comprise monitoring a fifth measurement region for density variations within the fluid, as the fluid passes the fifth measurement region, monitoring a sixth measurement region for density variations within the fluid, as the fluid passes the sixth measurement region, and determining whether a fifth density variation, which is detected by monitoring the fifth measurement region, and a sixth density variation, which is detected by monitoring the sixth measurement region, match, wherein a third axis passes through the center of the fifth measurement region and the center of the sixth measurement region, the first axis and the third axis are not parallel and preferably perpendicular to each other, and the second axis and the third axis are not parallel and preferably perpendicular to each other to allow for a three-dimensional analysis of the fluid flow.
  • The computing device of the fluid flow analyzer may be configured to neglect intensity variations and/or local maxima caused by sound waves.
  • The emitter of the fluid flow analyzer may be a laser and the beams (of electromagnetic radiation) may be laser beams.
  • The fluid flow analyzer may comprise optical components. The optical components may be configured to split a beam emitted from the emitter and/or focus the resulting beams to a spot in the first and second measurement region, respectively.
  • For example, the spots may lie on a line of flow to increase the accuracy when determining the speed of the fluid flow.
  • The computing device may be configured to neglect intensity variations and/or local maxima caused by sound waves by filtering intensity variations that are below a threshold.
  • The computing device may be configured to determine a speed of the fluid flow between the first measurement region and the second measurement region based on the temporal shift and a distance between the centers of the first measurement region and the second measurement region.
  • The distance between the centers of the first measurement region and the second measurement region may be equal to or less than 100 times, 50 times, 20 times, 10 times or one time a size of the first measurement region or the second measurement region (along an axis through the centers of the first measurement region and the second measurement regions).
  • The first, second, third and fourth beam may be derived by splitting a beam emitted by the emitter.
  • For example, the beam may be split by two Wollaston prisms.
  • The fluid flow analyzer may further comprise a third measurement region delimited by a fifth optical path of a fifth beam (of electromagnetic radiation) and a sixth optical path of a sixth beam (of electromagnetic radiation), a fourth measurement region delimited by a seventh optical path of a seventh beam (of electromagnetic radiation) and an eighth optical path of an eighth beam (of electromagnetic radiation), a third sensor for detecting a third intensity of the superimposed fifth and sixth beams, and a fourth sensor for detecting a fourth intensity of the superimposed seventh and eighth beams, wherein the computing device is configured to determine a temporal shift at which a correlation between the fourth intensity and the fifth intensity has a local maximum.
  • The first optical path and the fifth optical path may be not parallel and preferably perpendicular to each other to allow for a two-dimensional analysis of the fluid flow.
  • It will be appreciated that the features and attendant advantages of the disclosed method may be realized by the disclosed apparatus and vice versa. Moreover, it is noted that throughout the description, features in brackets are to be regarded as optional.
  • Brief Description of Drawings
  • The foregoing aspects and many of the attendant advantages will become more readily appreciated as the same becomes better understood by reference to the following description of embodiments, when taken in conjunction with the accompanying drawings, wherein like reference numerals refer to like parts throughout the various views, unless otherwise specified.
    • Fig. 1 and Fig. 1a illustrate the operation of a fluid flow analyzer which may be useful for understanding the present invention, but does not form part of the invention.
    • Fig. 2 ad Fig. 2b illustrate the operation of a fluid flow analyzer according to a first embodiment;
    • Fig. 3, Fig. 3a, and Fig. 3b illustrate the operation of a fluid flow analyzer according to a second embodiment; and
    • Fig. 4 shows a flow chart of a method for analyzing a fluid flow.
  • Notably, the drawings are not drawn to scale and unless otherwise indicated, they are merely intended to conceptually illustrate the structures and procedures described herein.
  • Detailed Description
  • Fig. 1 shows a fluid flow analyzer which does not form part of the invention.
  • The fluid flow analyzer 10 comprises an emitter 12a (e.g., a laser) which emits a beam 14a of coherent electromagnetic radiation (e.g., a laser beam). The beam 14a is split into a first beam 16a and a second beam 16b by optical components 18a arranged within the optical path of the beam 14a. The first beam 16a and the second beam 16b traverse a fluid flow 20 and generate an interference pattern that is incident on a first sensor 22a.
  • The intensities detected by the first sensor 22a depend on the density of the fluid. This allows detecting a density variation within a fluid volume that crosses the beams 16a, 16b.
  • I.e., the bending angles of the first and second beams 16a, 16b depend on the refractive indices of the fluid and its surroundings and thus on the density of the fluid. Because the fluid flow 20 exhibits (relatively) cooler parts 20a and (relatively) warmer parts 20b, the bending angles of the first and second beams 16a, 16b change over time when the fluid passes through the measurement region 24a (schematically illustrated by the broken rectangle). As the first beam 16a and the second beam 16b are displaced from each other in the direction of the fluid flow 20 (indicated by the arrows), there is a delay between the bending pattern exhibited by the first beam 16a and the bending pattern exhibited by the second beam 16b.
  • This causes a characteristic change in intensity which can be detected by the first sensor 22a on which the interfering first and second beams 16a, 16b incide. Accordingly, the region through which the first and second beams 16a, 16b pass, forms a first measurement region 24a for detecting density variations within the fluid flow 20.
  • A further measurement region 24b is established downstream by replicating the structure described above. I. e., there is another emitter 12b (e.g., another laser) that emits a beam 14b of coherent electromagnetic radiation (e.g., a laser beam) which is split into a third beam 16c and a fourth beam 16d by optical components 18b arranged within the optical path of the beam 14b. The third beam 16c and the fourth beam 16d traverse the fluid flow 20 and generate an interference pattern that is detected by a second sensor 22b.
  • Because the fluid flow 20 exhibits (relatively) cooler parts 20a and (relatively) warmer parts 20b, the bending angles of the third and fourth beams 16c, 16d change over time when the fluid passes through the measurement region 24b (schematically illustrated by the broken rectangle). As the third beam 16c and the fourth beam 16d are displaced from each other in the direction of the fluid flow (indicated by the arrows), there is a delay between the bending pattern exhibited by the third beam 16c and the fourth beam 16d.
  • This causes a characteristic change in intensity which is detected by the second sensor 22b on which the interfering third and fourth beams 16c, 16d incide. Accordingly, the region through which the third and fourth beams 16c, 16d pass, forms a second measurement region 24b for detecting density variations within the fluid flow 20.
  • The fluid flow analyzer 10 further comprises a computing device 26 for determining a temporal shift Δt at which a correlation C(Δt) between the intensities detected by the sensors 22a, 22b has a local maximum. I. e., the computing device 26 strives at finding matching density variations which allows calculating a speed of the fluid flow by dividing a distance between the measurement regions 24a, 24b by the temporal shift Δt.
  • In order to avoid that the result is corrupted by the influence of sound waves that cause density variations that travel through the fluid at the speed of sound, the computing device 26 may be configured to neglect intensity variations and/or local maxima caused by sound waves. As the density variations caused by sound waves are often smaller in amplitude than the density variations caused by temperature variations, the computing device 26 may be configured to filter intensity variations that are below a threshold. Moreover, sound waves caused by an isotropic turbulent flow may rather widen a correlation peak than adversely affect the calculation of the speed of the fluid.
  • Fig. 2 illustrates a modified fluid flow analyzer 10 which is based on two focused laser differential interferometers, FLDIs and differs from the fluid flow analyzer 10 of Fig. 1 in that the dimensions of the measurement regions 24a, 24b (perpendicular to the direction of flow) are reduced to facilitate matching and to improve spatial resolution (perpendicular to the direction of flow) with regard to determining the speed of the flow by dividing the distance D of the measurement regions 24a, 24b by the temporal shift Δt.
  • Heretofore, the optical components 18a arranged within the optical path of the beam 14a are replaced by a lens 28a with a short focal length, a quarter wave plate 30a, a Wollaston prism 32a, and a lens 34a that focusses the first and second beams 16a, 16b to different spots within the fluid flow 20. After having traversed the measurement region 24a, the first and second beam 16a, 16b are guided onto the sensor 22a by an inverse optical setting comprising a lens 36a, a Wollaston prism 38a, a quarter wave plate 40a, and a lens 42a with a short focal length.
  • As in Fig. 1, a further measurement region 24b is established downstream by replicating said structure. I. e., the optical components 18b arranged within the optical path of the beam 14b are replaced by a lens 28b with a short focal length, a quarter wave plate 30b, a Wollaston prism 32b, and a lens 34b that focusses the third and fourth beams 16c, 16d to different spots within the fluid flow 20. After having traversed the measurement region 24b, the third and fourth beam 16c, 16d are guided onto the sensor 22b by an inverse setting comprising a lens 36b, a Wollaston prism 38b, a quarter wave plate 40b, and a lens 42b with a short focal length.
  • As a difference in optical path length which causes a change in regard to the interference pattern seen by the sensors 22a, 22b mainly results from different densities within the measurement regions 24a, 24b around the focal points, the spatial resolution (perpendicular to the direction of flow) is increased.
  • The intensity values detected by the sensors 22a, 22b are provided to the computing device 26 which is configured to determine a temporal shift Δt at which a correlation C(Δt) between the intensities detected by the sensors 22a, 22b has a local maximum. The temporal shift indicates the time required for a fluid volume carrying a density variation to travel from the upstream measurement region 24a to the downstream measurement region 24b and allows calculating the speed of the fluid flow between the measurement regions 24a, 24b by dividing the distance D between the (centers of the) measurement regions 24a, 24b by the temporal shift Δt. The Distance (D) between the centers of the measurement regions (24a, 24b) is 0.01 m or less.
  • Fig. 3 and Fig. 3a illustrate a modified fluid flow analyzer 10 (that may be used for analyzing supersonic or hypersonic fluid flows 20). The modified fluid flow analyzer 10 differs from the fluid flow analyzer 10 of Fig. 2 in that the distance D between the measurement regions 24a, 24b (in the direction of flow) is reduced to facilitate matching and increase spatial resolution (in the direction of flow) with regard to measuring the speed of the flow. To this end, the beam 14 emitted by the emitter 12 (e.g., a laser) is expanded by a lens 28 with a short focal length and the expanded beam 14 is circularly polarized by a first quarter wave plate 30 that has an optical axis which is rotated by 45° relative to the direction of polarization of the incident beam 14. The circularly polarized expanded beam 14 is then split into two beams 14a, 14b by a first Wollaston prism 32 which may bend the beams 14a, 14b by less than 5 degrees (e.g., two degrees).
  • Beam 14a is circularly polarized by a second quarter wave plate 44 and split into a first beam 16a and a second beam 16b by a second Wollaston prism 46, wherein the direction of polarization of the first beam 16a may be perpendicular to the direction of polarization of the second beam 16b. Beam 14b is circularly polarized by the second quarter wave plate 44 and split into a third beam 16c and a fourth beam 16d by the second Wollaston prism 46, wherein the direction of polarization of the third beam 16c may be perpendicular to the direction of polarization of the fourth beam 16d. The second Wollaston prism 46 may bend the beams 16a, 16b, 16c, 16d by less than 2.5 degrees (e.g., one degree).
  • The beams 16a, 16b, 16c, 16d may be focused to spots within the fluid flow 20. To increase the accuracy when determining the speed of the fluid flow 20, the spots may be placed on a line of flow. If a density variation travels along the line of flow, it will be detected subsequently by both sensors 22a, 22b.
  • After traversing the first measurement region 24a, the first and second beams 16a, 16b incide on a third Wollaston prism 38a and a polarizer 48a and cause an interference pattern on the first sensor 22a (e.g., a photodiode which outputs a signal that is proportional to the intensity of the incident electromagnetic radiation). After traversing the second measurement region 24b, the third and fourth beams 16c, 16d incide on a fourth Wollaston prism 38b and a polarizer 48b and cause an interference pattern on the second sensor 22b (e.g., a photodiode which outputs a signal that is proportional to the intensity of the incident electromagnetic radiation).
  • As shown in Fig. 3a, a distance D between the (centers of the) measurement regions 24a, 24b may be equal to a distance between the focal points of the first beam 16a and the second beam 16b or to a distance of the focal points of the third beam 16c and the fourth beam 16d. The distance D between the centers of the measurement regions (24a, 24b) is or less. Moreover, by using two or three analyzers that measure along non-parallel (e.g., perpendicular) axes, the speed of the fluid may be measured in two or three dimensions.
  • Fig. 3b shows graphs from two experiments that are based on the setting described with reference to Fig. 3 and Fig. 3a. The top left graph of Fig. 3b shows a first signal S1 output by the first sensor 22a, a second signal S2 output by the second sensor 22b and a signal H indicating the onset of the fluid flow 20 (from the moment of signal rise). The top right graph of Fig. 3b shows a first signal S1 output by the first sensor 22a and a second signal S2 output by the second sensor 22b. The middle and bottom graphs show the local maxima (peaks) of a correlation function of the first signal S1 and the second signal S2. As shown in Fig. 3b, the analyzer 10 is robust against mechanical shocks, as the frequency of such shocks is typically much smaller than the frequency of the signals S1, S2 which allows filtering mechanical shocks without degrading the accuracy of the analyzer 10.
  • Fig. 4 shows a flow chart of a method for analyzing a fluid flow that may be practiced based on any one of the analyzers 10 of Fig. 1 to Fig. 3. The method starts at step 50 with monitoring the first measurement region 24a for density variations within the fluid, as the fluid passes the first measurement region 24a. For instance, density variations of the fluid passing through the first measurement region 24a may be detected by interferometry. The method continues at step 52 with monitoring the second measurement region 24b for density variations within the fluid, as the fluid passes the second measurement region 24b. For instance, density variations of the fluid passing through the second measurement region 24b may be detected based on interferometry. At step 54, the method ends with determining whether a first density variation, which is detected by monitoring the first measurement region 24a, and a second density variation, which is detected by monitoring the second measurement region 24b, match.
  • Reference signs list
  • 10
    fluid flow analyzer
    12
    emitter
    12a
    emitter
    12b
    emitter
    14
    beam
    14a
    beam
    14b
    beam
    16a
    first beam
    16b
    second beam
    16c
    third beam
    16d
    fourth beam
    18a
    optical components
    18b
    optical components
    20
    fluid flow
    22a
    sensor
    22b
    sensor
    24a
    measurement region
    24b
    measurement region
    26
    computing device
    28
    lens
    28a
    lens
    28b
    lens
    30
    quarter-wave plate
    30a
    quarter-wave plate
    30b
    quarter-wave plate
    32
    Wollaston prism
    32a
    Wollaston prism
    32b
    Wollaston prism
    34
    lens
    34a
    lens
    34b
    lens
    36
    lens
    36a
    lens
    36b
    lens
    38a
    Wollaston prism
    38b
    Wollaston prism
    40a
    quarter-wave plate
    40b
    quarter-wave plate
    42a
    lens
    42b
    lens
    44
    quarter wave plate
    46
    Wollaston prism
    48a
    polarizer
    48b
    polarizer
    50
    process step
    52
    process step
    54
    process step
    D
    distance

Claims (10)

  1. A method for analyzing a fluid flow (20), the method comprising:
    monitoring (50) a first measurement region (24a) for density variations within the fluid, as the fluid passes the first measurement region (24a), wherein monitoring (50) the first measurement region (24a) for density variations within the fluid comprises monitoring the first measurement region (24a) for variations of the refractive index of the fluid, as the fluid passes the first measurement region (24a);
    monitoring (52) a second measurement region (24b) for density variations within the fluid, as the fluid passes the second measurement region (24b), wherein monitoring (52) the second measurement region (24b) for density variations within the fluid comprises monitoring the second measurement region (24b) for variations of the refractive index of the fluid, as the fluid passes the second measurement region (24b);
    wherein monitoring the first measurement region (24a) for variations of the refractive index of the fluid is based on interferometry and monitoring the second measurement region (24b) for variations of the refractive index of the fluid is based on interferometry;
    characterized in that
    monitoring the first measurement region (24a) for variations of the refractive index of the fluid is based on measuring an intensity of interfering electromagnetic waves which travel along different paths through the first measurement region (24a);
    monitoring the second measurement region (24b) for variations of the refractive index of the fluid is based on measuring an intensity of interfering electromagnetic waves which travel along different paths through the second measurement region (24b); and
    determining whether a first density variation, which is detected by monitoring (50) the first measurement region (24a), and a second density variation, which is detected by monitoring (52) the second measurement region (24b), match comprises calculating a temporal correlation between the intensity of the interfering electromagnetic waves which travel along the different paths through the first measurement region (24a) and the intensity of the interfering electromagnetic waves which travel along the different paths through the second measurement region (24b); wherein
    the first measurement region (24a) is delimited by a first optical path of a first beam (16a) and a second optical path of a second beam (16b), wherein the first optical path differs from the second optical path; wherein
    the second measurement region (24b) is delimited by a third optical path of a third beam (16c) and a fourth optical path of a fourth beam (16d) wherein the third optical path differs from the fourth optical path;
    wherein the first and second beams (16a, 16b) are focused to different spots within the fluid flow (20);
    wherein the third and fourth beams (16c, 16d) are focused to different spots within the fluid flow (20); and
    wherein a distance (D) between centers of the measurement regions (24a, 24b) is 0.01 m or less.
  2. The method of claim 1, further comprising:
    determining a speed of the fluid flow (20) between the first measurement region (24a) and the second measurement region (24b) based on a time period (Δt) that passes between a first instant at which the first density variation is detected in the first measurement region (24a) and a second instant at which the second density variation is detected in the second measurement region (24b) and the distance (D) between the centers of the measurement regions (24a, 24b).
  3. The method of claim 1 or 2,
    wherein the first density variation and the second density variation are caused, at least in part, by turbulence; and/or
    wherein the first density variation and the second density variation are caused by a temperature variation resulting from said turbulence.
  4. The method of claim 1, further comprising:
    determining a speed of sound within the fluid flow (20) from a width of a peak of a function of the calculated temporal correlation.
  5. The method of any one of claims 1 to 4, further comprising:
    monitoring a third measurement region for density variations within the fluid, as the fluid passes the third measurement region;
    monitoring a fourth measurement region for density variations within the fluid, as the fluid passes the fourth measurement region; and
    determining whether a third density variation, which is detected by monitoring the third measurement region, and a fourth density variation, which is detected by monitoring the fourth measurement region, match;
    wherein
    a first axis passes through the center of the first measurement region and the center of the second measurement region;
    a second axis passes through the center of the third measurement region and the center of the fourth measurement region; and
    the first axis and the second axis are not parallel and preferably perpendicular to each other to allow for a two-dimensional analysis of the fluid flow.
  6. The method of claim 5, further comprising:
    monitoring a fifth measurement region for density variations within the fluid, as the fluid passes the fifth measurement region;
    monitoring a sixth measurement region for density variations within the fluid, as the fluid passes the sixth measurement region; and
    determining whether a fifth density variation, which is detected by monitoring the fifth measurement region, and a sixth density variation, which is detected by monitoring the sixth measurement region, match;
    wherein
    a third axis passes through the center of the fifth measurement region and the center of the sixth measurement region;
    the first axis and the third axis are not parallel and preferably perpendicular to each other; and
    the second axis and the third axis are not parallel and preferably perpendicular to each other to allow for a three-dimensional analysis of the fluid flow.
  7. A fluid flow analyzer (10), adapted to perform the method for analyzing a fluid flow (20) of one of claims 1 to 6, comprising:
    an emitter (12, 12a, 12b) for emitting electromagnetic radiation;
    wherein the fluid flow analyzer (10) is characterized by
    a first measurement region (24a) delimited by a first optical path of a first beam (16a) and a second optical path of a second beam (16b), wherein the first optical path differs from the second optical path; and
    a second measurement region (24b) delimited by a third optical path of a third beam (16c) and a fourth optical path of a fourth beam (16d), wherein the third optical path differs from the fourth optical path;
    a first sensor (22a) for detecting a first intensity of first electromagnetic radiation of a first interference pattern caused by the superimposed first and second beams (16a, 16b) on the first sensor (22a);
    a second sensor (22b) for detecting a second intensity of second electromagnetic radiation of a second interference pattern caused by the superimposed third and fourth beams (16c, 16d) on the second sensor (22b); and
    a computing device (26) for determining a temporal shift (Δt) at which a correlation between the first intensity and the second intensity has a local maximum;
    wherein the first and second beams (16a, 16b) are focused to different spots within the fluid flow (20);
    wherein the third and fourth beams (16c, 16d) are focused to different spots within the fluid flow (20); and
    wherein a distance (D) between centers of the measurement regions (24a, 24b) is 0.01 m or less.
  8. The fluid flow analyzer (10) of claim 7,
    wherein the computing device (26) is configured to neglect intensity variations and/or local maxima caused by sound waves; or
    wherein the computing device (26) is configured to neglect intensity variations and/or local maxima caused by sound waves by filtering intensity variations that are below a threshold.
  9. The fluid flow analyzer (10) of claim 7 or 8,
    wherein the computing device (26) is configured to determine a speed of the fluid flow (20) between the first measurement region (24a) and the second measurement region (24b) based on the temporal shift (Δt) and the distance (D) between the centers of the first measurement region (24a) and the second measurement region (24b).
  10. The fluid flow analyzer (10) of claim 9,
    wherein the first, second, third and fourth beam (16a, 16b, 16c, 16d) are derived by splitting a beam (14) emitted by the emitter (12).
EP19204453.5A 2019-10-21 2019-10-21 Fluid flow analysis Active EP3812712B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
EP19204453.5A EP3812712B1 (en) 2019-10-21 2019-10-21 Fluid flow analysis

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
EP19204453.5A EP3812712B1 (en) 2019-10-21 2019-10-21 Fluid flow analysis

Publications (3)

Publication Number Publication Date
EP3812712A1 EP3812712A1 (en) 2021-04-28
EP3812712B1 true EP3812712B1 (en) 2024-04-24
EP3812712C0 EP3812712C0 (en) 2024-04-24

Family

ID=68296324

Family Applications (1)

Application Number Title Priority Date Filing Date
EP19204453.5A Active EP3812712B1 (en) 2019-10-21 2019-10-21 Fluid flow analysis

Country Status (1)

Country Link
EP (1) EP3812712B1 (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001090700A2 (en) * 2000-05-23 2001-11-29 Arizona Board Of Regents Novel method and apparatus for flow monitoring in micro-fluidic devices
US20190293543A1 (en) * 2018-03-26 2019-09-26 F-Lab S.R.L. Method and apparatus for measuring the properties of a liquid

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3909621A1 (en) * 1989-03-23 1990-09-27 Deutsche Forsch Luft Raumfahrt Speed measuring method
GB9011086D0 (en) * 1990-05-17 1990-07-04 Jackson Roger G Tomographic monitoring of fluid flows
US6128072A (en) * 1998-04-23 2000-10-03 Nova Gas Transmission Ltd. Optical flow meter integrally mounted to a rigid plate with direct optical access to the interior of a pipe
GB9916932D0 (en) * 1999-07-16 1999-09-22 Funes Gallanzi Marcelo Method and apparatus for investigating fluid flow variables element characteristics and near-surface temperature and forces
US6473705B1 (en) * 2000-10-10 2002-10-29 General Electric Company System and method for direct non-intrusive measurement of corrected airflow
GB0107184D0 (en) * 2001-03-22 2001-05-09 Yortech Instr Ltd Gas velocity measurement
US7503227B2 (en) * 2005-07-13 2009-03-17 Cidra Corporate Services, Inc Method and apparatus for measuring parameters of a fluid flow using an array of sensors
JP2009168688A (en) * 2008-01-17 2009-07-30 Mitsubishi Heavy Ind Ltd Fluid measuring device
US8960017B2 (en) * 2012-11-14 2015-02-24 Daniel Measurement And Control, Inc. System and method for ultrasonic metering using an orifice meter fitting

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001090700A2 (en) * 2000-05-23 2001-11-29 Arizona Board Of Regents Novel method and apparatus for flow monitoring in micro-fluidic devices
US20190293543A1 (en) * 2018-03-26 2019-09-26 F-Lab S.R.L. Method and apparatus for measuring the properties of a liquid

Also Published As

Publication number Publication date
EP3812712C0 (en) 2024-04-24
EP3812712A1 (en) 2021-04-28

Similar Documents

Publication Publication Date Title
Ceruzzi et al. Simultaneous velocity and density gradient measurements using two-point focused laser differential interferometry
CA2576978C (en) Optical sensor using low-coherence interferometry
US8040499B2 (en) Transmitted light refractometer
US8687204B2 (en) Method and apparatus for measuring refractive index based on a ratio between a number of second fringes divided by a difference of the number of second fringes minus a number of first fringes
Houpt et al. Focused and cylindrical-focused laser differential interferometer characterization of SBR-50 at Mach 2
Meier et al. Imaging laser Doppler velocimetry
Camillo et al. A low-effort and inexpensive methodology to determine beam separation distance of multi-foci FLDI
JPH1144641A (en) Method and apparatus for measuring refractive index distribution
EP3812712B1 (en) Fluid flow analysis
JP2005106835A (en) Method and apparatus for measuring refractive index distribution
Li et al. Numerical investigation of multi-beam laser heterodyne measurement with ultra-precision for linear expansion coefficient of metal based on oscillating mirror modulation
Fagan et al. Improvement in Rayleigh scattering measurement accuracy
JP2006189411A (en) Measuring instrument and measuring method for phase delay
EP3608653B1 (en) Apparatus and method for measuring particle size using backscattered light
KR102008253B1 (en) Multi channel optical profiler based on interferometer
Lucki et al. Fiber Optic and Free Space Michelson Interferometer—Principle and Practice
TWI843249B (en) Apparatus for fluid measurement
JPS63233305A (en) Displacement sensor
Mielke-Fagan et al. Rayleigh scattering measurements using a tunable liquid crystal Fabry-Perot interferometer
RU2497090C2 (en) Method for measurement of medium pulse pressure and device for its realisation (versions)
KR101617729B1 (en) Apparatus and methods for vibration using change of the amount of light, Analysis system for analyzing sample and Method thereof
Pfister et al. Heterodyne laser Doppler distance sensor with phase coding measuring stationary as well as laterally and axially moving objects
Kaebe Microsecond-Resolved Thermometry in Transient Shock Tube Flows using Tuneable Diode Laser Absorption Spectroscopy
Mielke et al. Multiple-point mass flux measurement system using Rayleigh scattering
KR20100032768A (en) Ultrasonic measuring apparatus using resonator length selective optical fiber interferometer

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE APPLICATION HAS BEEN PUBLISHED

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

AX Request for extension of the european patent

Extension state: BA ME

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE

17P Request for examination filed

Effective date: 20211025

RBV Designated contracting states (corrected)

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: EXAMINATION IS IN PROGRESS

17Q First examination report despatched

Effective date: 20220317

REG Reference to a national code

Ref country code: DE

Ref legal event code: R079

Ref document number: 602019050772

Country of ref document: DE

Free format text: PREVIOUS MAIN CLASS: G01F0001660000

Ipc: G01F0001708600

Ref country code: DE

Ref legal event code: R079

Free format text: PREVIOUS MAIN CLASS: G01F0001660000

Ipc: G01F0001708600

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: GRANT OF PATENT IS INTENDED

INTG Intention to grant announced

Effective date: 20231201

RIC1 Information provided on ipc code assigned before grant

Ipc: G01F 1/661 20220101ALI20231117BHEP

Ipc: G01N 21/41 20060101ALI20231117BHEP

Ipc: G01N 21/45 20060101ALI20231117BHEP

Ipc: G01P 5/26 20060101ALI20231117BHEP

Ipc: G01P 5/18 20060101ALI20231117BHEP

Ipc: G01F 1/712 20060101ALI20231117BHEP

Ipc: G01F 1/7086 20220101AFI20231117BHEP

RIN1 Information on inventor provided before grant (corrected)

Inventor name: SANDER, TOBIAS

Inventor name: MUNDT, CHRISTIAN

Inventor name: THELE, MAX

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE PATENT HAS BEEN GRANTED

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

REG Reference to a national code

Ref country code: DE

Ref legal event code: R096

Ref document number: 602019050772

Country of ref document: DE

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

U01 Request for unitary effect filed

Effective date: 20240508

U07 Unitary effect registered

Designated state(s): AT BE BG DE DK EE FI FR IT LT LU LV MT NL PT SE SI

Effective date: 20240517

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20240824

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: HR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20240424

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: GR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20240725

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: GB

Payment date: 20240912

Year of fee payment: 6