METHOD AND APPARATUS FOR MEASUREMENT OF AN ANALYTE
20230393064 · 2023-12-07
Inventors
Cpc classification
A61B5/097
HUMAN NECESSITIES
International classification
A61B5/097
HUMAN NECESSITIES
G01J3/44
PHYSICS
Abstract
The invention relates to a method for measurement and monitoring of analytes, for example ketones, particularly but not exclusively acetone in human or other mammalian breath. The invention also relates to apparatus for use in performance of the method. The present invention uses a spectroscopic technique known as CELIF which is a direct combination of the well-established and powerful laser-spectroscopic techniques cavity ring-down spectroscopy and laser-induced fluorescence. The method utilises a flow body to control a flow of sample gas through a laser beam.
Claims
1. An apparatus for quantifying an analyte in a gas or liquid phase sample comprising: a light source with an emission spectrum that overlaps with an absorption of the analyte, a pair of reflective mirrors located on an optical axis to form an optical cavity, the cavity having a sample inlet and a sample outlet; a fluorescence detector located at a location not on the cavity axis and arranged to provide a first signal in response to fluorescence within the cavity; a photon detector located axially external to the cavity and arranged to provide a second signal; wherein the apparatus including the light source, the cavity and the axial photon detector comprises a cavity-enhanced absorption spectrometer; wherein the apparatus including the fluorescence detector is configured to comprise a cavity-enhanced laser-induced fluorescence (CELIF) spectrophotometer; means for supplying an analyte-free gas sample or an analyte-containing gas sample to the cavity through the inlet; a processor adapted to receive a first signal from the fluorescence detector and a second signal from the axial photon detector, and further adapted to provide a measurement of the analyte concentration in the sample: wherein a flow body is located between the sample inlet and sample outlet, the flow body comprising a chamber extending along a direction of sample flow between the sample inlet and the sample outlet, the flow body further comprising a light source inlet and a light source outlet outlet arranged to provide a path for a light source through the chamber, the path extending transversely of the direction of sample flow; the flow body further comprising an aperture communicating with the fluorescence photon detector in a direction transverse to the path of the light source.
2. The apparatus of claim 1, wherein the light source is a laser and the path of the light source is a laser beam.
3. The apparatus of claim 2, wherein the laser is a sub-nanosecond-pulsed Nd:YAG laser.
4. The apparatus of claim 1, wherein the flow body further comprises one or more flow channels.
5. The apparatus of claim 1, wherein the chamber has a dimension perpendicular to the direction of sample flow, wherein the dimension increases from a minimum value at the sample inlet to a maximum value in the vicinity of the light source path and decreasing to a value smaller than the maximum value at the sample outlet.
6. The apparatus of claim 1, wherein the dimensions of the chamber increase or decrease smoothly.
7. The apparatus of claim 1, wherein the direction of sample flow, the laser beam path and the fluorescence detector aperture are perpendicular or orthogonal and intersect at a single point.
8. The apparatus of claim 1, wherein the surface of the chamber may have a smooth profile configured to minimize turbulent flow of the sample gas, or formation of eddies, during passage through the chamber.
9. The apparatus of claim 1, wherein the width of a cross section of the chamber taken parallel to the sample flow increases gradually from the inlet to a maximum value in the vicinity of the laser beam path (for example just before the laser beam path) and may decrease towards the outlet, so that the velocity of the sample gas decreases as it flows towards the vicinity of the laser beam and increases from the vicinity of the laser beam to the outlet
10. The apparatus of claim 1, wherein the chamber may comprise three sections: a first section adjacent to the sample inlet; an optional second, central section; and a third section adjacent to the sample outlet.
11. The apparatus of claim 10, wherein the first section has a circular cross section and the circular cross section of the first section increases in size in the direction of the sample flow.
12. The apparatus of claim 11, wherein the second section has a circular cross section, remaining at a constant size in the direction of sample flow.
13. The apparatus of claim 12, wherein the third section has a circular cross section and the circular cross section decreases in size in the direction of sample flow.
14. The apparatus of claim 13, wherein the first section is configured as an expanding cone starting from the narrowest point adjacent to the inlet and expanding to a widest point in the direction of the sample flow.
15. The apparatus of claim 12, wherein the second section is configured as a cylinder.
16. The apparatus of claim 13, wherein he third section is configured as a narrowing cone going from a widest point adjacent to the second section or the first section and decreasing in width to a narrowest point of the cone adjacent to the sample outlet.
17. A method of quantifying an analyte in a gas phase sample comprising the steps of: providing apparatus in accordance with the first aspect of the invention; introducing an analyte-free gas sample into the chamber obtaining a background signal; introducing a flowing reference sample containing a quantity of the analyte into the chamber and obtaining a reference signal; processing the background signal and reference signal to obtain an absolute absorption coefficient of the reference sample and to provide a calibration signal of the CELIF spectrophotometer; introducing a specimen gas sample containing the analyte into the chamber obtaining a specimen signal; processing the specimen signal and the calibration signal to obtain a measurement of the absolute concentration of the analyte in the specimen sample; wherein the measurement is made while the specimen sample is flowing through the chamber.
18. A method as claimed in claim 17, wherein the analyte concentration is measured repeatedly overtime intervals of 100 ns to 5 ms.
19. A method of claim 17, wherein an average value for the analyte concentration is measured over a period from 1.0 ms to 0.03 s.
20. A method of detecting or monitoring Type 1 diabetes in a patient comprising the step of measuring acetone in exhaled breath using the method as claimed in claim 17.
21. A method of detecting ketone acidosis in a patient comprising the step of: measuring acetone in exhaled breath using apparatus as claimed in claim 1.
22. A method as claimed in claim 21, wherein the apparatus further comprises a breath collector connected to an inlet of the apparatus.
23. A method as claimed in claim 22 wherein the breath collector comprises a manifold, face mask or breathing tube.
Description
[0068] The present invention is now discussed in more detail with reference to the Figures.
[0069]
[0070]
[0071]
[0072]
[0073]
[0074]
[0075]
[0076]
[0077]
[0078]
[0079]
[0080]
[0081]
[0082]
[0083]
[0084]
[0085]
[0086]
[0087]
[0088]
[0089]
[0090]
[0091]
[0092]
[0093]
[0094]
[0095]
[0096]
[0097]
[0098]
[0099]
[0100]
[0101]
[0102]
[0103]
[0104]
[0105] The width of the chamber or dimension of the chamber along the y axis and in the vicinity of the laser beam may have a maximum value so that a maximum amount of the gaseous sample is irradiated by the laser beam. In alternative embodiments, the width of the chamber or the internal dimension may reach a maximum value after or before the sample flow passes the laser beam path.
[0106] Expanding the chamber in the dimension of the laser/cavity axis (the y axis) and the fluorescence detection axis (z axis), compared to the chamber inlet and the chamber outlet, reduces the relative area of the two holes for the laser compared to the internal surface area of the chamber. If the internal surfaces of the chamber are designed well, then the flow will be parallel across the laser beam inlet and the laser beam outlet minimising loss of sample gas into the side arms that enclose the cavity.
[0107] The laser beam inlet and the laser beam outlet are preferably not sealed by windows. This is because windows could absorb and/or scatter light from the laser beam reducing the effectiveness of the ring-down aspect of the technique.
[0108] As shown in
[0109] In an embodiment, the apparatus may further comprise a second fluorescence detector located at a location not on the cavity axis and arranged to provide a first signal in response to fluorescence within the cavity. Such an embodiment is depicted in
[0110] Configuration of the Chamber
[0111] The chamber within the flow body controls the flow of sample across the laser beam path. Preferably, the chamber is configured such that the sample flow is laminar across the laser beam path with limited eddies or recirculation or capture of the sample within the beam path.
[0112] As discussed above, the width of a cross section of the chamber taken parallel to the sample flow (in the x-z plane and/or the x-y plane, referring to
[0113]
[0114] In embodiments, the chamber may follow the general configuration of the chamber discussed above, whereby the cross section of the chamber taken parallel to the sample flow (in the x-z plane and/or the x-y plane, referring to
[0115] In an embodiment, the width of the chamber (in the x-z plane and/or the x-y plane, referring to
[0116] The chamber may comprise three sections: a first section adjacent to the sample inlet; a second, central section; and a third section adjacent to the sample outlet. The second section may be absent, accordingly, in embodiments the chamber is formed of the first section and the third section.
[0117] In an embodiment the cross-sectional configuration of the chamber in the y-z plane (the plane parallel to the beam path and parallel to the axis of the fluorescence detector) is generally circular at the inlet. Accordingly, the first section may have a circular cross section. In an embodiment the circular cross section of the first section of the chamber increases in size in the direction of the sample flow. The second section may similarly have a circular cross section. However, the circular cross section of the second section remains a constant size. The third section also may have a circular cross section with the size of the cross section decreasing in the direction of sample flow until the third section meets the sample outlet.
[0118] In an embodiment the cross-sectional configuration of the chamber transitions from a circular configuration to a configuration where the cross section is circular with flattened sides. The point at which the circular cross section transitions into the flattened circular cross section marks the transition from the first section to the second section.
[0119] A flattened circular cross sectional configuration may take the form of a circle where a pair of opposing flat sides are present or where two pairs of opposing flat sides are present. Accordingly, the second section may have such a cross-sectional configuration. The cross-sectional configuration may transition from circular to quadrilateral, via the flattened circular shape previously mentioned. Such a transition may occur within the second section. Continuing in the direction of flow towards the outlet, the cross-sectional configuration may transition from a quadrilateral to a flattened circular shape to a circular shape or from a flattened circular shape to a circular shape, as appropriate, with the size of the cross section decreasing towards the outlet. The transition from a quadrilateral to a flattened circular shape to a circular shape or from a flattened circular shape to a circular shape marks the transition from the second section to the third section.
[0120] The discussion above revolved around a cross section of the chamber in the y-z plane (the plane parallel to the beam path and parallel to the axis of the fluorescence detector. We now consider the chamber cross section in the x-y plane (parallel to the laser beam path, parallel to the laser beam path, and perpendicular to the axis of the fluorescence detector).
[0121] The first section may be configured as an expanding cone starting from the narrowest point adjacent to the inlet and expanding to a widest point in the direction of the sample flow. The widest point of the cone is where the first section and second section meet or where the first section and the third section meet. The second section may be configures as a cylinder. The third section may be configured as a narrowing cone going from a widest point adjacent to the second section or the first section where the second section is absent and decreasing in width to a narrowest point of the cone adjacent to the sample outlet. The second section and the third section may meet at the widest point of the third section.
[0122] The first section may have a cross section in the x-y plane and/or the x-z plane shaped as an expanding cone starting from the narrowest point adjacent to the inlet and expanding to a widest point in the direction of the sample flow. The widest point of the cone is where the first section and second section meet or where the first section and the third section meet (in embodiments where there is not a second section). The second section may have a cross section in the x-y plane and/or the x-z plane shaped as quadrilateral shape or a quadrilateral shape with rounded corners. The third section may have a cross section in the x-y plane and/or the x-z plane shaped as a narrowing cone going from a widest point adjacent to the second section or the first section where the second section is absent and decreasing in width to a narrowest point of the cone adjacent to the sample outlet. The second section and the third section may meet at the widest point of the third section.
[0123] The cross-sectional shape of the chamber may have a plane of symmetry. The cross sectional shape of the chamber may have one or two planes of symmetry.
[0124] The second section may be the laser beam inlet and the laser beam outlet or may comprise the laser beam inlet and the laser beam outlet. The laser beam inlet and the laser beam outlet are on opposite sides of the second section and are positioned to provide a path for the laser beam through the chamber. The laser beam inlet and the laser beam outlet may alternatively be positioned in the first section or in the third section. Optionally, the laser beam outlet and the laser beam inlet are present at wider areas of the first or third section.
[0125] In an embodiment the chamber comprises a first section having a cross section parallel to the laser beam path and perpendicular to the axis of the fluorescence detector shaped as an expanding cone starting from the narrowest point adjacent to the inlet and expanding to a widest point and a third section having a cross section parallel to the laser beam path and perpendicular to the axis of the fluorescence detector shaped as a narrowing cone going from a widest point and decreasing in width to a narrowest point of the cone adjacent to the sample outlet, wherein the first section and the third section meet at their respect widest point.
[0126] It would be evident from the discussion above that the chamber may have an internal shape that is symmetrical. The plane of symmetry may be present at the laser beam path.
[0127] In embodiments the flow body further comprises one or more side flow channels. Preferably, the flow body comprises two side flow channels. The side flow channels may be placed at any position within the flow body yet outside of the chamber. In embodiments, two side flow channels flank the chamber running parallel to the sample flow direction. The side flow channels have an inlet to the chamber to allow gas to flow from the side flow channels into the chamber.
[0128] The side flow channels may be configured to provide air or inert gas inlet channels which may be provided at or upstream adjacent the laser beam inlet and outlet apertures in the flow body, so that air or inert gas flows across the laser beam inlet or outlet preventing ingress of sample, analyte-containing gas into the laser beam inlet or outlet apertures. This may reduce or avoid sample, analyte-containing gas becoming trapped in the inlet or outlet apertures and consequently being measured multiple times or concentrating in the cavity.
[0129] In embodiments the air or inert gas channels (which may be referred to as side flow channels) may join the chamber at the same point as the laser beam inlet and laser beam outlet. Alternatively, the air or inert gas channels may join the chamber proximate to the sample outlet.
[0130] In embodiments the side flow channels extend along the length of the flow body from the sample inlet end of the flow body to the sample outlet end of the flow body.
[0131] In embodiments the side flow channels may intersect with the optical cavity such that gas flows across the laser beam path to reduce or prevent ingress of sample, analyte-containing gas into the optical cavity. In addition, or in alternative embodiments, the side flow channels communicate with the chamber proximate the sample outlet via an inlet.
[0132] We have discovered that non-laminar flow of the analyte gas may give a lag in the response of the detector, particularly if repeated sampling of several breaths takes place with a resultant blurring of the observed output signal.
[0133] Further, transverse flow of the sample gas through the flow body may give a higher limit of detection of the CELIF device, of more than 2000 ppm, superior to the cavity ring down measurement and covering the dynamic range of the CELIF device.
[0134] In previously used methods for analysis of acetone in breath, a patient may exhale into a bag. This may provide a sort of average value but gives an artificial result.
[0135] Further, the procedure needs to be repeated whenever a measurement is required.
[0136] The present invention allows a patient to breathe naturally with a measurement being provided in real time for each exhaled breath. This facilitates measurements of acetone levels in breath of children or sleeping or unconscious patients. Also, measurements of animals' breath is facilitated.
[0137] The flow body of the present invention allows CELIF measurements of flowing samples. The flow body was designed so that it delivers the sample from the about 6 mm (about ¼ inch) tubing system before the laser beam path, to the measurement region with gradual increase in diameter, which may match the width of the field of view of the fluorescent detector. Then, gradually decrease the width of the flow body to match the 6 mm (about ¼ inch) tubing system at the sample outlet. As seen in
[0138] When viewed from the side (in the x-z plane), the sample inlet concavely expands from a narrow dimension to a wider dimension (however, this is less than the dimension that the flow body expands to when viewed from the top the x-y plane). The wider dimension in this case may be 8 mm). The flow body then returns to the initial diameter at the sample outlet.
[0139] At the top of the flow body there is an aperture (optionally in the size of 25 mm) for the LIF optics to be fitted.
[0140] The flow body is placed inside a six-way cross. The six-way cross consists of the laser beam path, the fluorescence detector and the sample flow.
[0141] Before testing the transverse flow measurements through the flow body it was essential to first compare CELIF measurements of static gas samples which fill the whole volume of the cavity with CELIF measurements of flowing gas samples through the cavity axis, which fill the same volume, to test whether the flow of gas affects the concentration measurements.
[0142] These tests showed that static CELIF measurements do not reflect the real acetone concentration in a sample. This may be due to optical bleaching of the stagnant acetone molecules in the cavity, which does not occur while flowing the sample as the molecules that interact with the laser get replaced quickly preventing bleaching. Thus, CRD and CELIF measurements must be made with flowing gas samples.
[0143] It was found that for the real-time CELIF measurements to be quick and to minimize the residence time of the samples in the CELIF device, gas samples should flow transversely through the flow body, from the sample inlet to the sample outlet at an angle to the beam path (optionally perpendicular) and not along the beam path, this can be seen in
[0144] Signals of different standard acetone-nitrogen gas mixtures were measured once with longitudinal flow and once with transverse flow, with a flow rate of 0.5 slpm in each flow direction. The measurements of each gas mixture were compared and are shown in
[0145]
[0146] CELIF calibration followed by CELIF measurements of different samples were made to explore the lower and higher limits of detection, dynamic range, of the CELIF device. Standard acetone-nitrogen gas mixtures of acetone concentration in the range between 2-100 ppm were used to do the CELIF calibration. The outlet of the standard gas bottles was connected to the mass flow controller that controls the pressure at the cavity inlet. The outlet of the CELIF device was connected to the second mass flow controller that controls the flow rate of the gas through the cavity.
[0147] Each CELIF calibration measurement was made by first measuring the laser independent background, γ, by blocking the laser then recording the LIF signal. Then measuring the background CELIF signal for an acetone-free sample, which was measured by longitudinally flowing 1 bar of N.sub.2 at 0.5 slpm. Then, the sample CELIF signal was measured by longitudinally flowing some of the acetone-nitrogen gas mixture at a pressure of 1 bar and a flow rate of 0.5 slpm. The sample CELIF measurement was recorded after 1 minute of the gas flow to ensure that the cavity pressure which was controlled by the mass flow controller was stable. Each background and sample CELIF measurement was an average of 1500 laser shots. The cavity was flushed with flowing N.sub.2 after each measurement. After each measurement the cavity is flushed with N.sub.2, for 2 minutes along the sample path, followed by 2 minutes along the laser path and then 1 minute along the sample path. A typical calibration measurement is based on five individual measurements of analyte concentrations ranging from 0 ppm to 100 ppm. Next, to produce a calibration curve the CELIF signals were plotted against the acetone concentration as measured by CRD. Then, a straight line with zero intercept was fitted to the data, yielding a calibration factor ′=1.49±0.03×10.sup.−5. The three standard deviation (3σ) limit of detection (LOD) of the cavity ring-down measurements was calculated to be 2 parts per million (ppm), when averaging over every 1500 laser shots in 100 ms, using
[0148] where, τ.sub.0 is the analyte-free ring-down time, δτ.sub.0 is the error in τ.sub.0 and σ is the absorption cross section of acetone at 266 nm. The higher limit of detection of the CRD measurements was estimated from the ability to fit a CRD transient to a single exponential decay to measure the ring-down time, which gave acetone concentration of 390 ppm as the higher limit.
[0149] Next, a high concentration acetone-nitrogen mixture was prepared in a 0.5 L gas cylinder by first pumping down the cylinder, then releasing some of the gas acetone into it, then topping up the cylinder with 9 bar of gas N.sub.2. The mixture then was left to mix thoroughly, for example for about 1 hour with the help of the heated side arms of the mixing cylinder. The gas line was connected to the mass flow controller and the flow body inlet, and the flow body exit was connected to the second mass flow controller. Then, γ was measured by blocking the laser and recording the LIF signal, and the background CELIF signal was measured by transversely flowing 1 bar of N.sub.2 through the flow-body at 0.5 slpm. After that, 1 bar of the acetone-nitrogen mixture was released into the chamber of the flow body with a flow rate of 0.5 slpm, and after 1 minute of the gas flow a CELIF measurement was recorded which was an average of 1500 laser shots. Next, the cavity was flushed with N.sub.2, and the rest of the gas mixture was diluted by topping it up with N.sub.2 and was left to mix for about 15 minutes. Then, another set of background and sample CELIF measurements were taken as above. The procedure was repeated many times until the limit of detection was reached.
[0150] ′, and was appended to the CELIF calibration graph. The acetone concentration CELIF limit of detection was calculated to be 1.6 ppm using
[0151]
[0152] For the apparatus of the present invention to be clinically accepted, the acetone CELIF measurement must be able to follow the real-time breathing pattern, by responding quickly to any change in the acetone concentration. The CELIF signal must rise and reach the maximum acetone concentration in less than 1 s, then fall and return to the baseline level before the next breath arrives. The rise time of the measurement is defined as the time taken for the signal to rise from 10% to 90% of the steady state signal, and the fall time of the measurement is the time taken for the signal to fall from 90% to 10% of the steady state of the signal.
EXAMPLES
[0153] Four flow bodies forming part of the invention were tested. Flow body 1 is shown in
Example 1
[0154] The performance of flow body 1, shown in
[0155] CELIF Response Time with Flow Body 1
[0156] A CELIF background measurement was made by measuring the laser independent background, γ, by blocking the laser then recording the LIF signal. Then measuring the background CELIF signal by flowing lab air through the flow body. Next, the CELIF signal was recorded every 100 ms by averaging 1500 laser shots, while the 3-way valve was set to flow lab air to record a baseline for the measurement, then after a few seconds, the 3-way valve was turned to the sample side and allowing the sample to flow for a few seconds before turning the 3-way valve back to the lab air side. The CELIF measurement was stopped when the baseline returned to the initial value.
[0157] The result of this measurement is shown in
[0158] The rise time of the CELIF signal was 6.2±0.2 s and the fall time was 6.6±0.2 s. Therefore, at least 12 s was needed for a single acetone measurement to be made using this flow body.
[0159] Sample Diffusion Through Side Arms with Flow Body 1
[0160] Time lapse acetone CELIF measurements were made to test whether the gas samples flowing thought the flow body diffuse into the side arms of the cavity, or if the samples being diluted by the stagnant air in the cavity side arms. A mixture of acetone in air was made using the method discussed above (creating a mixture of acetone vapour and air by putting a few drops of liquid acetone in a 3 L glass bulb which was open to the atmosphere). The mixture was flowed through the flow body at flow rate of 0.5 slpm for more than 1 hour, while the LabView program was recording the CRD, LIF and CELIF signals every 1 minute, averaging 1500 laser shots for each measurement.
[0161] The results showed that both the time-integrated CRD signal and the LIF signal decayed with time. As the gas extended into the side arms through the two laser holes in the flow body, the sample length increased and the light intensity at the middle section of the cavity decreased, which caused both the CRD and LIF signals to decay with time. Even though the sample diffused into the side arms of the cavity, the CELIF signal was constant for the entire period of the experiment as expected.
[0162] Computational Fluid Dynamics Simulations
[0163] To investigate the reason behind the long rise and fall time of the CELIF measurement made using flow body 1, and to find a suitable solution to this issue, Computational Fluid Dynamics (CFD) simulations were performed to visualise the theoretical flow through the flow body. Autodesk CFD was used for this purpose. These CFD simulations of the different designs of the flow body were made by J. Landes, “Optimisation of Analyte-Gas Flow Through a CELIF Sensor for Breath Acetone”, Masters thesis, Durham University, 2018, the contents of which is incorporated herein by reference.
[0164] Torun a CFD simulation for a flow body, a CAD drawing of the object must be created and imported into the CFD software. The materials of the various parts of the CAD object must be assigned including the flow volume. The boundary conditions of the flow system must be set including the volume flow rates of the sample gas at the inlets and the static pressure on the outlets. The samples used in the CELIF experiment contained only a small concentration of acetone, which would have a negligible impact on the flow of the gas. Thus, the fluid material “air” was used in all of the simulations. The following step was to create a mesh for the system for the 3D simulation. The mesh is where the simulation is actually performed, which is a series of tetrahedral elements which approximate the geometry of the CAD object. The partial differential equations used by the CFD software to describe the flow of fluid do not have solutions for complex systems, therefore the system is divided into smaller subsystems within which analytical solutions can be found. A fine mesh will match the geometry and result in an accurate simulation, however this costs a large amount of computer memory and long computational time. A coarse mesh will result in a fast simulation, but will produce inaccurate results.
[0165] The general method used to create a mesh for an object in the CFD software for the simulations performed in this project is as follows: first, using the auto-size function a simple coarse mesh was created. Then, the size adjustment was used to refine the mesh, which increases or decreases the size of the mesh elements by up to a factor of 5. Using this, the mesh was refined to give an element count between one and two million. Next, “surface refinement” and “gap refinement” were enabled to refine the mesh on the surfaces and between gaps in the CAD geometry. The final step was to adjust the wall layer settings. A wall layer is a mesh along the wall of the flow volume used to simulate the boundary layer. Therefore, it is important that the wall layer is thick enough and contains enough layers to resolve the whole boundary layer. The number of layers within the wall layer settings was set to 9, the layer factor, which determines the layer thickness, was set to 0.2 and the layer gradation, which controls the rate of growth of the wall layers, was kept at the automatic setting. This method generated the initial mesh for each CFD simulation.
[0166] Certain control settings had to be chosen to run the CFD simulations, such as if the flow is laminar or turbulent and if the fluid is incompressible. A flow is characterized to be laminer (smooth) or turbulent (rough) by a dimensionless parameter called Reynolds number: Re=(ρdυ)/η where ρ and η are the fluid's density and viscosity, dis the diameter of the flow channel, and υ is the flow velocity.
[0167] Values of Re<2000 predict a laminar flow, whereas values of Re>2500 usually indicate that flow will be turbulent. In a laminar flow system, some vortices, or localized swirling, can separate from the central streamlines within a flow channel, which unlike turbulent flow, fluid in vortices is composed of slowly moving currents and streamlines. The Mach number M is a dimensionless quantity representing the ratio of flow velocity to the speed of sound. At values of M<1, the compressibility of the fluid can be ignored.
[0168] The values of the Mach and Reynolds numbers of our flow system are significantly lower than the Mach and Reynolds number at which a fluid beings to show compressibility and turbulence, respectively. Therefore, it was determined that the simulations will use the approximations that the flow is laminar and incompressible to reduce the computational cost.
[0169] The simulation results must be mesh independent. Thus, mesh adaptation was enabled to ensure this was true for the simulations. This is a program which analyses the results and refines the mesh where needed. The results after running mesh adaptation were always mesh independent.
[0170] CFD Simulation of Flow Body 1
[0171] A CFD simulation was performed on flow body 1. The CAD geometry used for this simulation was as shown in
[0172] A set of traces spread evenly over the inlet were plotted to see the path the gas takes as it enters the flow body. The traces are a packet of data points which enter the system and map out the path the gas takes based on the velocity of the fluid.
[0173] The path of the traces of flow body 1 are shown in
Example 2
[0174] Flow body 2 was designed with a chamber similar to the shape of the laminar flow found in the simulation of flow body 1. The CAD drawing for this new design is shown in
[0175] Adding these two channels required the modification of the main flow tube body found in the simulation of flow body 1. This included increasing the width of the chamber after the laser beam path. The side channels were designed so that the air flow in the side channels would join the sample flow next to the beam inlet and beam outlet. This is so that at the point of measurement, minimal mixing would occur, reducing any effect the side channel flow would have on the sample concentration. This flow body prototype was 3D printed from a plastic material (PLA polylactic acid) after being studied with CFD simulations.
[0176] CFD Simulation of Flow Body 2
[0177] The design of flow body 2 shown in
[0178] Similar CFD simulation to that of flow body 1 was performed on the simplified flow body 2. First, only the main flow tube was being studied, so no volume flow rate boundary condition was placed on the side channel inlets. The trace profile mapping the flow through flow body 2 is shown in
[0179] Next, the effect of air flow through the side channels was studied using CFD simulations. The magnitude of the side channel flow rate was chosen so that the average velocity of the sample at the middle of the flow matches the velocity of the air in the side channels. From the previous simulations it was found that the averaged velocity of the sample at a sample flow rate of 0.5 L/min was 3 cm/s. Using this and the cross section area of the side channel, the flow rate at the inlet of each side channel was set to 0.013 L/min. The trace profile mapping the flow through the main flow chamber and the side channels of flow body 2 is shown in
Performance of Flow Body 2
[0180] CELIF Response Time with Flow Body 2
[0181] A CELIF response time test was made using flow body 2, first without flowing air through the side channels, and using the same experiment and sampling method as Example 1 as what was used with flow body 1. The inlets of the side channels were covered to prevent any flow through them. The CELIF signal from this measurement was plotted together with the CELIF signal measured with flow body 1 as shown in
[0182] The Effect of the Side Channel Flows
[0183] In order to test the effect of the side channel flows, the sampling method presented in Example 1 had to be slightly modified. As shown in
[0184] Time lapse acetone CELIF measurements were made to test whether the gas samples flowing through flow body 2 diffuse into the side arms of the cavity. A mixture of acetone in air was made and was flowed through the flow body at a flow rate of 0.5 slpm, first with no side channel flows, for more than 1 hour, while the CRD, LIF and CELIF signals were recorded every 1 minute, averaging 1500 laser shots for each measurement. This measurement was repeated with flowing air through the side channels at a flow rate of 0.025 slpm through each channel.
[0185] The results are shown in
[0186] Next, the effect of the side channel flow on the CELIF measurement was tested.
[0187] Using the sampling method shown in
[0188] Furthermore, the CRD signal of the measurement with side channel flows recovers faster than the CRD signal of the measurement without side channel flows when switching the 3-way valve to air after flowing the acetone mixture. This further reinforces the conclusion that the side channel flows helps in reducing the amount of the sample gas diffusing into the laser beam inlet and outlet of the chamber, and agrees with the results of the diffusion test. Also shown in the CRD signals in
[0189] The speed of the CELIF measurement with flow body 2 showed significant improvement compared to flow body 1. Even though the side channel flow proved to reduce the amount of the sample diffusing into the side arms, it caused reduction of the sample concentration. Taking these issues into account, a third flow body was designed as discussed in the following section.
Example 3
[0190] The vortex-free flow shape found in the CFD simulations of flow body 1 was modified in flow body 2 to fit in the side flow channels, which caused the formation of vortices in the sample flow in flow body 2. The design of the new flow body (flow body 3), is shown in
[0191] CFD Simulation of Flow Body-3
[0192] CFD simulation similar to the previous simulations was run for this flow body. The flow rate boundary condition for the sample inlet was set to 0.5 L/min, and the side channel flow rate was determined using the same method as flow body 2 and was set to 0.02 L/min through each channel. The boundary condition for the gas outlet was set to 0 Pa gauge pressure. The traces showed that no vortices are formed in the sample flow, and that no mixing occurs between the sample flow and the side channels flow at the middle of the flow body. The traces show that the sample would fill the total width of the main flow volume, thus, when this flow body was 3D printed from PLA plastic the LIF optics hole was set to be 19 mm.
[0193] Performance of Flow Body-3
[0194] CELIF Response Time with Flow Body-3
[0195] Similar to previous flow body designs, a CELIF response time test was made using the sampling method presented in
[0196] The Effect of the Side Channel Flows
[0197] Two CELIF response time measurements were made with flow body 3 using the same gas mixture. The first measurement was made with no side channel flows, and the second one was made with 0.05 slpm air flow through each channel. The two measurements are plotted together in
[0198] However, beneficially the slope in the CRD signal for the measurement without side channel flows, compared to the flat CRD signal in the measurement with side channel flows suggests that the side channel flows act as barriers between the flow regions and reduce the amount of the sample gas that diffuse into the side arms.
[0199] Time lapse acetone CELIF measurements were made with flow body 3 at different side channel flow rates. A mixture of acetone in air was made and was flowed through the flow body at a flow rate of 0.5 slpm for 15 minutes, while the CRD, LIF and CELIF signals were recorded every 1 minute, averaging 1500 laser shots for each measurement. This measurement was repeated flowing air through the side channels at flow rates of 0.025, 0.05 and 0.1 slpm through each channel.
[0200] The results are shown in
[0201] However, as shown in the acetone concentration plot, flowing air through the side channels with even a small flow rate caused dilution of the sample. Even with the improved design of the side flow channels, mixing between the sample and the side channel flows occurred, thus, the use of side channel flows should be avoided.
Example 4
[0202] CELIF Response Time
[0203] At this stage, no trials were made on patients. Thus, a gas sampling method which does not involve using a 3-way valve and interrupting the pressure inside the cavity was implemented.
[0204] The CELIF response time measurement was made by first allowing the lab air only to flow through the system, then valve b was open for a few seconds and then quickly closed to allow the sample to join the air flow into the system. This was repeated a second time, while the CRD, LIF and CELIF signals were recorded every 100 ms, averaging 1500 laser shots in each measurement.
[0205] The result is shown in
where, for the first acetone peak, A.sub.1 is the acetone concentration, tri is the centre of the rising edge, τ.sub.r1 is the time constant of the rising edge, t.sub.f1 is the centre of the falling edge and τ.sub.f1 is the time constant of the falling edge. The rest of the constants are the same for the second acetone peak. From that the 10% rise time
[0206] of the signal peak, t.sub.10, is
[0207] And the 90% rise time of the signal peak, t.sub.90, is
[0208] Thus, the 10-90% rise time of the signal is
2τ.sub.r ln(9)
[0209] and the same follows for the 90-10% fall time of the signal,
2τ.sub.f ln(9)
[0210] From the fit it was found that the 10-90% rise time of the CELIF measurement was 370±15 ms, and the 90-10% fall time was 850±21 ms, where the errors are from one standard deviation. The rise and fall time of the CELIF measurement are fast enough to follow a real breath pattern. While the two times should be equal, the longer fall time is likely to be caused by some acetone sticking into the tube fittings before the flow body which takes some time to clear up.
Example 5
Cavity Cleaning Between CELIF Measurements
[0211] To test how long it takes the cavity to return to CELIF background level after introducing a sample of high acetone concentration, a CELIF measurement was made by allowing lab air to flow into flow body 3 without side channel flows, flowing an acetone-air mixture for about 20 s, then turning the flow back to only lab air, as shown in
[0212] Next, the same measurement was repeated as shown in
[0213] While monitoring real patients breath, with a lower acetone concentration than the acetone used in this example, we could tolerate some acetone build up in the side arms as long as it does not affect the CELIF measurement. However, eventually the build up will become intolerable, and the cavity will need to be flushed by the method described above. The time after which the cavity will need to be cleaned depends on the amount of acetone concentration being monitored. Therefore, it is recommended to clean the cavity at least every 1 minute while monitoring a real patient's breath with acetone concentration of more than 30 ppm, and every 2 minutes if the acetone concentration is less than that.
Example 6
[0214] Flow Body 4
[0215] The design of the flow body is shown in
[0216]
in which [Ac](t) is the measured acetone concentration as a function of time, t, A.sub.i is the amplitude of one acetone pulse, t.sub.r,i and t.sub.f,i are the times of the rising and falling edges, respectively, and τ.sub.r,i and τ.sub.f,i are the rise and fall times, respectively. From the fit it was found that the 10-90% rise time of the CELIF measurement was 370±15 ms, and the 90-10% fall time was 850±21 ms, where the errors are one standard deviation.
[0217] In order to investigate the accuracy of the CELIF instrument and to fulfill the requirement for the CELIF instrument to be medically approved, the concentrations of acetone measured by the CELIF instrument must be compared with a proven analytical technique. A selected ion flow tube mass spectrometer (SIFT-MS) instrument (Voice200Ultra—Syft Technologies) was provided by Anatune to the Department of Chemistry at Durham University and was used to validate the CELIF instrument.
[0218] SIFT-MS is a new analytical technique for the simultaneous, real-time quantification of several trace compounds in air or breath samples as described in P. Španiěl and D. Smith, “Progress in SIFT-MS: Breath analysis and other applications,” Mass Spectrometry Reviews, vol. 30, pp. 236-267, July 2010; D. Smith and P. Španiěl, “SIFT-MS and FA-MS methods for ambient gas phase analysis: Developments and applications in the UK,” Analyst, vol. 140, no. 8, pp. 2573-2591, 2015; and D. Smith and P. Španiěl, “Selected ion flow tube mass spectrometry (SIFT-MS) for on-line trace gas analysis,” Mass Spectrometry Reviews, vol. 24, no. 5, pp. 661-700, 2005, the contents of which are hereby incorporated by reference.
[0219] CELIF Validation Measurements
[0220] A SIFT-MS instrument (Voice200Ultra—Syft Technologies) was provided by Anatune and was used for the CELIF validation measurements. It was placed such that concomitant readings could be taken on the same gas sample. The gas sampling setup is shown in
[0221] The CELIF device was first calibrated by using standard acetone-nitrogen gas mixtures in the range of 1-100 ppm (SIP Analytical, SIPCYL 110 LTR non-refillable can filled certified grade acetone in nitrogen). Each gas mixture was measured simultaneously by the two instruments. CRD and LIF measurements were recorded for CELIF calibration by flowing the gas longitudinally through the cavity axis and simultaneously the acetone concentration in the sample was measured by the SIFT-MS. Both the CELIF cavity and the SIFT-MS were flushed between the measurements. Each CELIF measurement was an average of 1500 laser shots and was recorded using a 0.5 slpm sample gas flow. Each SIFT-MS measurement was an average of acetone concentrations measured over 20 s of gas flow, with a measurement interval of 140 ms.
TABLE-US-00002 TABLE 1 A summary of the SIFT-MS instrument (Voice200Ultra - Syft Technologies) method used for the validation of the CELIF device. Carrier gast type Helium Reagent ion NO.sup.+ (30 amu) Production ion NO.sup.+CH.sub.3COCH.sub.3 (88 amu) Carrier Flow (tls) 5.0057 Sample Flow (tls) 0.3 Tube Temp (C.) 123.75 Tube Pressure (Torr) 0.6348 Reaction Time (ms) 8.203
[0222] Next, to validate the CELIF measurements with the SIFT-MS measurements, the procedure above was repeated using the same standard gas bottles but flowing the gas samples transversely through the flow body and using the previous CELIF calibration to calculate the acetone concentration. More validation measurements were also made using a homemade acetone-nitrogen mixture and a series of dilutions where each concentration was measured by CELIF and the SIFT-MS simultaneously. The SIFT-MS instrument was not optimised to measure acetone concentrations above about 100 ppm, thus the validation measurements were restricted to acetone concentrations between 1 ppm (CELIF limit of detection)-100 ppm.
[0223] The used SIFT-MS method is summarised in table 1. Surprisingly, a SIFT-MS instrument dependent correction factor was needed to correct for the acetone reaction rates which was important for accurate quantification of our validation measurements. In the (Anatune) lab, a gas sample that should give 6 ppm acetone was generated and measured with SIFT-MS using all the possible reagent ions, and from the data generated the reaction rates were amended. This yielded that a correction factor of 1.519 must be multiplied by the acquired SIFT-MS validation data.
[0224] The CELIF calibration measurements and the comparison between CRD and SIFT-MS measurements of acetone concentration are presented in
[0225] The acetone CELIF measurement procedure, the flow body design and the response time of the CELIF measurements were optimised. We achieved an acetone concentration dynamic range between 1.6-2000 ppm, covering the range of breath acetone concentration a DKA patient might have. The CELIF device was validated with the SIFT-MS.