Next Article in Journal
Electroreforming of Biomass for Value-Added Products
Previous Article in Journal
Main Problems Using DEM Modeling to Evaluate the Loose Soil Collection by Conceptual Machine as a Background for Future Extraterrestrial Regolith Harvesting DEM Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Fast Bubble Detection Method in Microtubes Based on Pulsed Ultrasound

1
School of Mechatronic Engineering, Xi’an Technological University, Xi’an 710021, China
2
School of Mechanical Engineering, Xi’an Jiaotong University, Xi’an 710049, China
*
Author to whom correspondence should be addressed.
Micromachines 2021, 12(11), 1402; https://doi.org/10.3390/mi12111402
Submission received: 3 September 2021 / Revised: 9 November 2021 / Accepted: 12 November 2021 / Published: 15 November 2021
(This article belongs to the Section E:Engineering and Technology)

Abstract

:
In the process of biological microfluidic manipulation, the bubbles generated in the tube will seriously reduce the gauging accuracy. This paper introduces an improving method that can estimate the size of microbubbles in real time. Hence, the measurement data of the liquid volume can be modified according to this method. A microbubble detector based on the pulsed-ultrasound method was studied, including the device structure and the working principle. The assessment formula of the microbubbles in the tube was derived from the simulation results, which adopted the two-phase theory. The digital image processing method was applied to fulfill the microbubble calibration. This detection method was applied to measure the microbubbles in the tube and to modify the flow volume in a timely manner. The results of the experiments showed that this method is effective at improving the microflow gauging accuracy.

1. Introduction

In clinical medicine and biological research, we often encounter situations requiring accurate fluid measurement, such as quantitative blood transfusion, biological microfluidic sample testing, etc. [1]. The bubbles in the flow undoubtedly disturb the fluid action, especially in the case of microfluidic processing. The bubbles can obstruct the semipermeable membrane in blood dialysis, disorder the indexes in the blood-gas analysis, and even put a patient under transfusion in danger [2,3]. In biomedical microfluidic experiments, the bubbles existing in the tube will seriously reduce the measurement accuracy. However, whether bubbles are generated depends on many factors, such as the pumping actuation, pressure differences along their path, different pressures between the inside and outside of the liquid surface, etc., so they are difficult to avoid and eliminate [4,5].
Although it is difficult to eliminate bubbles in the flow, it is possible to measure and control them [6]. Especially in microflow gauging, the bubble size and velocity are monitored in real time and compensated for in the fluid volume, which is a relatively effective method. Many technologies are currently adopted for the detection of microbubbles in tubes, such as capacitive detection [7], photoelectric detection [8], ultrasonic detection [9], image processing [10,11], etc. For capacitive detection, a high detection accuracy is difficult to obtain, and the photoelectric detection method is affected by the liquid’s color.
In recent years, due to the development of machine-learning techniques, bubble-detection methods using image processing have been widely studied, and such methods are able to obtain good results [12,13]. However, ultrasound is still a good choice to detect the bubbles in microtubes because of its good directional capability, penetrating ability, easy acquirement, greater accuracy, etc. Ultrasound techniques for bubble detection in flows can be classified into three groups: transmission techniques, Doppler techniques, and pulse techniques [14]. The transmission technique can measure bubble sizes. However, this technique has a low accuracy and position resolution. The Doppler technique is suitable for velocity field measurement, but is less used to measure the bubble size. The pulse technique measures the axial (along an ultrasonic beam) distance between the ultrasonic transmitter and a receiver or a reflecting object [15]. Applying the pulsed-ultrasound technique, the sizes and velocities of microbubbles in the flow in the tube can be obtained in real time. The detection process is not affected by the liquid color, nor the material of the microtubes.
In this paper, a fast bubble-detection method based on ultrasound technology is proposed. To verify the effectiveness of the method, a bubble detector was fabricated. This detector was applied to measure microbubbles in the tube and to modify the flow volume in a timely manner. The results of the experiment showed that the proposed method is effective at improving the microflow gauging accuracy.

2. Materials and Methods

2.1. The Principles of Bubble Movement in Microtubes and the Ultrasound Probe

The movement of a bubble with the flow in a microtube is different from that in a macrotube. If we observe microbubbles in the flow, there is a clear interface between the bubbles and the liquid. We considered the liquid as a continuous phase and the bubbles as a disperse phase. Adopting the FLUNT software, we built a traceable volume of fluid (VOF) model of the two-phase (gas–liquid) interface and studied its movement properties. From the simulation, we knew that during the process, the bubble shape had a slight change, similar to an ellipsoid, because of the effect of surface tension and viscous drag. The flow velocities of the bubble and liquid were almost equal if we kept the injection rate constant. The conclusion of the simulation corresponded to the observation experiment. The flow movement characteristics showed that the ultrasound was adapted to detect the microbubbles in the flow. The designed detection system was based on the ultrasonic transmission-medium-attenuation principle. Ultrasonic attenuation can be classified into three types: diffusion attenuation, absorptive attenuation, and scattering attenuation. The bubbles in the microtube would lead to the latter two attenuations. If the ultrasonic plane wave normal traverses the interface of the liquid–gas phase, its energy would be reflected and transmitted. We assumed the characteristic sound impedance of two phases as Z 1 and Z 2 . The sound intensity reflection ratio R and sound intensity transmission ratio T are defined as:
R = ( Z 2 Z 1 Z 2 + Z 1 ) 2
T = 4 Z 1 Z 2 ( Z 1 + Z 2 ) 2
The characteristic acoustic impedance of water and air are 1.5 × 10 6 kg/m 2 s and 430 kg/m 2 s, respectively. According to Equations (1) and (2), if the ultrasonic wave is incident from water to air, their interface sound intensity reflection ratio R is 0.999, and the sound intensity transmission ratio T is 0.001. This means that when the ultrasonic wave is incident from water to air, almost all of its energy is reflected on the boundary surface. In the detection system, we launched a high-frequency ultrasound pulse wave on one side of a pipeline with flow, and the wave through the pipeline was received on the other side. Once the bubbles appeared in the flow, the received wave energy decreased significantly because of absorptive attenuation and scattering attenuation. If the bubbles are large enough, ultrasound can only diffract at the bubbles’ edges and form an acoustic shadow. Most ultrasonic energy is absorbed by the bubble, so the received signal amplitude would be reduced to a very low value, and the received signal frequency would change. In the signal-processing circuit, adopting D/A conversion and signal interruption, we can obtain the variation in time of the ultrasonic signal frequency and voltage amplitude. Thus, we can obtain the interval time Δ t , which is the bubble going through the ultrasonic probe [16]. Because the speed of the bubble in the pipeline is approximately equal to the flow speed, we took the flow injection speed v as the bubble speed. If the bubble in the micropipeline is in the form of an air embolism, we can take the microtube diameter D as the bubble diameter d. The bubble volume V can be estimated simply as:
V = π 4 d 2 v Δ t
If the bubble is a sphere ( v Δ t < d ), the volume V can be calculated with the parameter of the high-accuracy calibration.

2.2. The Ultrasonic Probe Structure

The ultrasonic probe is a device that can transmit and receive ultrasound. Therefore, it is the key part of the microbubble-detection system. We selected the T40-160 and R40-160 piezoelectric ultrasonic devices (Audio Well Electronics Ltd., Guangzhou, China) as the transmitter and receiver. A higher ultrasonic frequency will detect smaller bubbles because of its smaller diffusion angle and better sonic wave directional characteristic, but the higher frequency signal has a lower anti-interference ability. In this paper, the ultrasonic frequency was 40 kHz, and the distance between the receiver and pipeline was no more than 10 mm. The structure of the probe is shown in Figure 1. We designed a sleeve joint between the transmitter and receiver. The micropipeline went through the side of the sleeve perpendicularly. The sleeve’s inner diameter was fit to the outer circumference of the micropipeline. This sleeve structure can lessen the sonic energy dissipation during the transmitting process. The prototype of the probe is shown in Figure 2, which is a top view.
In order to avoid contamination and facilitate the fit to the pipeline, we selected Delrin as the probe mounting material. Delrin possesses many excellent characteristics, such as acid and alkali resistance, wear resistance, good electrical properties, easy fabrication, low sensitivity to variations of the environmental temperature, etc. It is unfavorable if air exists in the installation gap of the parts, which would generate reflected wave clutter. We used PDMS to encapsulate the joints of the parts and utilized its flow-curing characteristic to eliminate air.

2.3. Design of the Circuit of the Microbubble Detection System

The circuit of the microbubble-detection system was composed of a pulsed-ultrasonic-signal-transmitting unit, an ultrasonic-signal-receiving unit, and a signal-processing unit. Its information flowchart is shown in Figure 3 In the pulsed-signal-emission unit, we adopted an SCM crystal system to launch square waves at 40 kHz, which had good integration and reliability, while we took advantage of the 74LS04 reverser to improve the driving capability. The AD620 (Analog Devices, Inc., Wilmington, MA, USA) operational amplifier, which can achieve an enlarged effect of 1000-times, was used as the signal-processing device of the ultrasonic receiver. When the signal was amplified, the noise component in the communication process was removed by the second-order filter circuit.
The key devices of the signal-processing unit were two AT89S51 SCMs (Atmel Corporation). The ADC chip was the AD7810 (Analog Devices, Inc.), whose upper and lower biases of the amplitude of the voltage were −0.4 V and +0.4 V. The storage unit was the CY7C132 (Cypress Semiconductor Corporation, San Jose, CA, USA). SCM 1 was responsible for the AD conversion and deposited the conversion results into the storage unit. SCM 2 was used to read the data of the storage and obtain the amplitude changes of the ultrasonic signal. Then, it could calculate the time that the bubbles took to go through the pipeline and the volume of the bubbles by comparing the voltage peak size. It could send messages to the PC, such as the magnitude of the voltage, the interrupt time, and the bubble volume. In the detection process, the energy attenuation of the ultrasonic signal within the bubbles was influenced by the temperature. In the experiment, we mounted a DS18B20 (Dallas Semiconductor, San Jose, CA, USA) intelligent temperature sensor to measure the temperature and fulfill the compensation calculation.

3. Experiments

In our developed microbubble-detection system, the signal processing was so complex that many factors could affect the measurement precision, such as the accuracy of the A/D conversion, the modulation enlargement, the shaping filter circuit, etc. In the experiment, we took advantage of the digital-image-processing technology to calibrate this system accurately. First, distilled water was injected and aspirated into a 1.5 mm-diameter transparent pipeline repeatedly via a syringe pump in order to produce bubbles randomly. The syringe pump was the CARVO XL-3000 (Tecan Group Ltd., Mannedorf, Switzerland). Then, we allowed the bubble–fluid to flow through the pipeline smoothly. We used a high-speed camera to acquire a bubble image before the bubbles went through the ultrasonic probe. The ultrasonic signal was calibrated using the data of the bubble characteristics from the images. Two formulas were used to calculate the volume of the bubbles. When v Δ t > 1.5 mm (diameter of the pipeline), the bubble volume was calculated with a cylinder; when v Δ t < 1.5 mm (diameter of the pipeline), the bubble volume was calculated with a sphere volume. Figure 4 shows one of the serial images and its calculating process. There are two bubbles in the image. Combining the mathematical morphology algorithm and edge-extraction algorithm, the image was enhanced by removing the background, and the bubble edges were extracted. An ellipses was drawn to fit the bubble outline in order to estimate the bubble volume. Figure 5 shows the ellipse-fitting result of the larger bubble projection in Figure 4.
Then, the image pixels were converted to the real size according to the camera parameters. The bubble projection area was equal to the ellipse area. The bubble volume was calculated with a sphere volume, the projection area of which was equal to the ellipse area. The bubble projection outline in Figure 4 was converted as follows: The long axis of the fitting ellipse was 0.649 mm, and the short axis was 0.527 mm. The bubble volume was calculated as 1.158 mm 3 . When the bubble was captured by the camera, we recorded the interrupt time of the ultrasonic detection system and the voltage signal from the ultrasonic probe. According to the bubble volume, which was calculated by the imaging system, we could calibrate the ultrasonic detection system. In the experiment, many bubbles of different sizes were detected and measured. Five bubbles of similar sizes were chosen to form a combination of sizes, and finally, ten sets with different sizes of bubbles were used for the calculation and analysis. The bubble size from the image was compared to the ultrasonic voltage signal. The time for bubbles to pass through the probe ( Δ t ) was recorded.

4. Results and Discussions

Table 1 shows one bubble for each size. The experimental results of the image-processing method and the ultrasonic method are included. The serial numbers represent 10 different sizes of bubbles. When the volumes of the bubbles were 1.158 μ L and 2.202 μ L, the interrupt time of the ultrasonic detection system was zero. When the volumes of the bubbles were 2.832 μ L, the voltage that came from the receiver of the ultrasonic detection system changed. However, the interrupt time of the system was still zero due to the error of the AD conversion and timing control of the AT89S51.
To avoid these errors, we determined the threshold of the ultrasonic detection system to be 3 μ L. When the volumes of the bubbles were less than 3 μ L, the voltage that came from the receiver of the ultrasonic detection system changed, and the interrupt time was zero. According to the experimental data, as long as the bubble volumes exceeded 2.832 μ L, there was a significant change in the signal of the ultrasonic detection system. Therefore, the detection threshold of the bubbles was 3 μ L. The results shown in Figure 6 also showed that the detection results had a small deviation from the image-processing results. The maximum deviation was 14.88%, and the minimum deviation was 5.94%.
In order to verify the effectiveness of the method under different parameter conditions, a sample-adding experiment was designed. That is, a given volume of liquid sample was added to a target vessel using a sample-adding system, which is usually a syringe pump. The volume of the liquid sample was compensated by measuring the volume of bubbles during the whole process using the proposed bubble-detection method. In the experiment, the bubbles were generated in the tube using the same method as above. Then, a 4 mL volume of the liquid was added using the sample-adding system via the tube. According to the volume of bubbles detected by the bubble-detection system, an additional corresponding volume was added to compensate the sample-adding system. Two sample-adding rates were used, i.e., different speeds of the liquid and bubbles when performing the bubble-detection method. After adding the sample, the liquid was weighed before and after the bubble detection using an electronic balance. Ten experiments were performed for each speed. The results are shown in Table 2 and Table 3.
From the results, it can be seen that the average volume of detected bubbles was 84.53 µL at a flow rate of 0.667 mL/s. The accuracy of the sample-adding system was improved by 1.876% after compensation by the bubble detection method. At flow rate of 0.067 mL/s, the average volume of detected bubbles was 73.48 µL, and the accuracy of spiking was improved by 1.363% after compensation. When the flow speed became slower, the bubble generation decreased. Regardless of the amount of bubbles generated, the proposed method can effectively detect the volume of bubbles and compensate the sample-adding system.

5. Conclusions

In this paper, we applied the ultrasonic-detection technique to the detection of bubbles in biological microsamples and developed an ultrasonic-bubble-detection sensor with the features of good applicability, easy manufacturing and installation, and high precision. We calibrated the bubble-detection system by the machine-microvision-measurement technology and determined the detection threshold of the bubble volume. The D/A conversion precision of the microcontroller and the interrupt timing error were the key parameters that affected the detection limit of the system.

Author Contributions

Conceptualization, Y.L. and L.F.; methodology, Y.L.; validation, J.W. (Junwu Wu); writing—original draft preparation, Y.L.; writing—review and editing, Y.L.; visualization, L.F.; experiment J.W. (Jinju Wang). All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Natural Science Foundation of Shaanxi Province (No. 2020JQ-813) and the Key Research and Development Program of Shaanxi Province (No. 2020GY-147).

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lorenz, M.G. Liquid-handling robotic workstations for functional genomics. JALA J. Assoc. Lab. Autom. 2004, 9, 262–267. [Google Scholar] [CrossRef]
  2. Bottenus, N.; Straube, T.; Farling, S.; Vesel, T.; Klitzman, B.; Deshusses, M.A.; Cheifetz, I.M. Ultrasonic bubble detection and tracking using spatial coherence and motion modeling. In Medical Imaging 2021: Ultrasonic Imaging and Tomography; International Society for Optics and Photonics: Bellingham, WA, USA, 2021; Volume 11602, p. 116020F. [Google Scholar]
  3. De Somer, F.M.; Vetrano, M.R.; Van Beeck, J.P.; Van Nooten, G.J. Extracorporeal bubbles: A word of caution. Interact. Cardiovasc. Thorac. Surg. 2010, 10, 995–1001. [Google Scholar] [CrossRef] [Green Version]
  4. Pereiro, I.; Khartchenko, A.F.; Petrini, L.; Kaigala, G.V. Nip the bubble in the bud: A guide to avoid gas nucleation in microfluidics. Lab Chip 2019, 19, 2296–2314. [Google Scholar] [CrossRef]
  5. Li, Y.; Liu, X.; Huang, Q.; Ohta, A.T.; Arai, T. Bubbles in microfluidics: An all-purpose tool for micromanipulation. Lab Chip 2021, 21, 1016–1035. [Google Scholar] [CrossRef] [PubMed]
  6. Mizushima, Y.; Sakamoto, A.; Saito, T. Measurement technique of bubble velocity and diameter in a bubble column via single-tip optical-fiber probing with judgment of the pierced position and angle. Chem. Eng. Sci. 2013, 100, 98–104. [Google Scholar] [CrossRef] [Green Version]
  7. Quoc, T.V.; Dac, H.N.; Quoc, T.P.; Dinh, D.N.; Duc, T.C. A printed circuit board capacitive sensor for air bubble inside fluidic flow detection. Microsyst. Technol. 2015, 21, 911–918. [Google Scholar] [CrossRef]
  8. Jin, T.; Xu, B.; Tang, K.; Hong, J.P. Bubble counter based on photoelectric technique for leakage detection of cryogenic valves. J. Zhejiang Univ.-Sci. A 2008, 9, 88–92. [Google Scholar] [CrossRef]
  9. Fletcher, S.M.P.; Ogrodnik, N.; O’Reilly, M.A. Enhanced detection of bubble emissions through the intact spine for monitoring ultrasound-mediated blood-spinal cord barrier opening. IEEE Trans. Biomed. Eng. 2019, 67, 1387–1396. [Google Scholar] [CrossRef] [PubMed]
  10. Haas, T.; Schubert, C.; Eickhoff, M.; Pfeifer, H. BubCNN: Bubble detection using Faster RCNN and shape regression network. Chem. Eng. Sci. 2020, 216, 115467. [Google Scholar] [CrossRef]
  11. Wang, Q.; Li, X.; Xu, C.; Yan, T.; Li, Y. Bubble recognizing and tracking in a plate heat exchanger by using image processing and convolutional neural network. Int. J. Multiph. Flow 2021, 138, 103593. [Google Scholar] [CrossRef]
  12. Kim, Y.; Park, H. Deep learning-based automated and universal bubble detection and mask extraction in complex two-phase flows. Sci. Rep. 2021, 11, 8940. [Google Scholar] [CrossRef]
  13. Cerqueira, R.F.; Paladino, E.E. Development of a deep learning-based image processing technique for bubble pattern recognition and shape reconstruction in dense bubbly flows. Chem. Eng. Sci. 2021, 230, 116163. [Google Scholar] [CrossRef]
  14. Povolny, A.; Batsaikhan, M.; Ihara, T.; Takahashi, H.; Kikura, H. Bubble Size Measurement by the Ultrasonic Pulse Echo with Tracking Technique. J. Flow Control. Meas. Vis. 2019, 7, 11–27. [Google Scholar] [CrossRef] [Green Version]
  15. Povolny, A.; Kikura, H.; Ihara, T. Ultrasound Pulse-Echo Coupled with a Tracking Technique for Simultaneous Measurement of Multiple Bubbles. Sensors 2018, 18, 1327. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Wang, T.; Wang, J.; Ren, F.; Jin, Y. Application of Doppler ultrasound velocimetry in multiphase flow. Chem. Eng. J. 2003, 92, 111–122. [Google Scholar] [CrossRef]
Figure 1. The structure of probe (section plane).
Figure 1. The structure of probe (section plane).
Micromachines 12 01402 g001
Figure 2. The prototype of the ultrasonic probe (top view).
Figure 2. The prototype of the ultrasonic probe (top view).
Micromachines 12 01402 g002
Figure 3. Flowchart of the bubble-detection system.
Figure 3. Flowchart of the bubble-detection system.
Micromachines 12 01402 g003
Figure 4. The processing of the bubble image. (a) Original image. (b) Background removed and brightness enhanced. (c) Edge detection. (d) Mathematical morphology expansion. (e) Mathematical morphology corrosion. (f) Bubble-edge-detection results.
Figure 4. The processing of the bubble image. (a) Original image. (b) Background removed and brightness enhanced. (c) Edge detection. (d) Mathematical morphology expansion. (e) Mathematical morphology corrosion. (f) Bubble-edge-detection results.
Micromachines 12 01402 g004
Figure 5. Ellipse fitting of the bubble.
Figure 5. Ellipse fitting of the bubble.
Micromachines 12 01402 g005
Figure 6. Comparative experiments of the image-processing and ultrasonic-bubble-detection system.
Figure 6. Comparative experiments of the image-processing and ultrasonic-bubble-detection system.
Micromachines 12 01402 g006
Table 1. The detected bubble volume of the bubble-detection system.
Table 1. The detected bubble volume of the bubble-detection system.
No.Bubble Volume from
Image Processing ( μ L)
Probe Amplitude (V)Bubble through Time (s)Bubble Volume from
Ultrasonic Detector ( μ L)
11.1580.41200
22.2020.39200
32.8320.38603.000
43.4150.3541.563.212
54.6360.3452.495.326
66.4490.3222.487.328
77.8210.2874.078.414
89.4830.2456.1110.162
910.2120.2346.3911.340
1012.110.2127.2013.071
Note: temperature: 20 °C; flow rate: 0.667 mL/s; the sample liquid is distilled water with a density of 0.9982 mg/mL.
Table 2. Results of the sample-adding experiments using the bubble-detection method for compensation (flow rate: 0.667 mL/s).
Table 2. Results of the sample-adding experiments using the bubble-detection method for compensation (flow rate: 0.667 mL/s).
No.Liquid Mass
without Compensation
(g)
Bubble Volume ( μ L)Liquid Mass
with Compensation
(g)
13.8946105.23.9857
23.8953109.63.9921
33.881384.23.9611
43.896098.13.9810
53.8903101.23.9803
63.892892.33.9728
73.887363.23.9413
83.887663.13.9425
93.882968.23.9429
103.901260.23.9536
Note: temperature: 20 °C; the sample liquid is distilled water with a density of 0.9982 mg/mL.
Table 3. Results of the sample-adding experiments using the bubble-detection method for compensation (flow rate: 0.067 mL/s).
Table 3. Results of the sample-adding experiments using the bubble-detection method for compensation (flow rate: 0.067 mL/s).
No.Liquid Mass
without Compensation
(g)
Bubble Volume ( μ L)Liquid Mass
with Compensation
(g)
13.962787.24.0416
23.952154.23.9914
33.946294.54.0293
43.975565.34.0127
53.9205108.34.0165
63.973758.64.0104
73.943658.93.9842
83.923699.64.0065
93.942156.33.9649
103.943651.93.9711
Note: temperature: 20 °C; the sample liquid is distilled water with a density of 0.9982 mg/mL.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, Y.; Wu, J.; Fu, L.; Wang, J. A Fast Bubble Detection Method in Microtubes Based on Pulsed Ultrasound. Micromachines 2021, 12, 1402. https://doi.org/10.3390/mi12111402

AMA Style

Li Y, Wu J, Fu L, Wang J. A Fast Bubble Detection Method in Microtubes Based on Pulsed Ultrasound. Micromachines. 2021; 12(11):1402. https://doi.org/10.3390/mi12111402

Chicago/Turabian Style

Li, Yiqing, Junwu Wu, Leijie Fu, and Jinju Wang. 2021. "A Fast Bubble Detection Method in Microtubes Based on Pulsed Ultrasound" Micromachines 12, no. 11: 1402. https://doi.org/10.3390/mi12111402

APA Style

Li, Y., Wu, J., Fu, L., & Wang, J. (2021). A Fast Bubble Detection Method in Microtubes Based on Pulsed Ultrasound. Micromachines, 12(11), 1402. https://doi.org/10.3390/mi12111402

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop