Ultrasonic and Impact-Echo Testing for the Detection of Scaling in Geothermal Pipelines

The efficiency and longevity of components in geothermal plants are significantly reduced by mineral deposition in pipes, also known as scaling [1]. A number of methods are used to monitor fouling (a broader term for all types of organic and inorganic deposition) in other industries, such as the oil, gas and food industries. These include direct and indirect, offline, in-situ and online, as well as inline methods, each of which have pros and cons [2, 3]. An effective monitoring solution specific to the high temperature conditions in geothermal plants, which optimises maintenance and cleaning measures in dependency of scaling, has yet to be developed. In this paper, two non-destructive testing techniques: contact ultrasonic testing (UT) and impact-echo (IE) testing are investigated in their viability as methods to detect scaling growth in geothermal plants. A descaling measurement was conducted on a heavily scaled segment of an obsolete pipeline from the production well of the geothermal power plant in Sauerlach, Germany. The pipeline segment was inserted in a test rig and the scaling was gradually etched away using an acidic solution. At regular time intervals, the scaling thickness was measured mechanically and contact UT as well as IE measurements were carried out. The testing apparatus was designed to withstand high temperatures (140∘C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${140}\,^{\circ }{\text {C}}$$\end{document} at the inlet pipe surface) and to be easy to install whilst being cost-efficient. Both techniques yielded usable results with submillimeter resolution. The advantages and limitations of the two methods are discussed. Impact-echo testing, in particular, can be automated and presents a simple and cost-efficient scaling monitoring option.


Introduction
With the increasing tendency toward renewable energy, geothermal plants are a growing industry [1][2][3]. From 1995 to 2005 a worldwide increase in geothermal energy production from 112441 to 587786.4 TJ/year was documented by the international geothermal association (IGA) [4]. The German geothermal association (BVG) lists a total of 38 running Scaling presents a major cause for loss of efficiency and part failure (such as in heat exchangers or pumps). Whilst a number of monitoring concepts are used to detect scaling in other industries, an effective scaling monitoring solution for the high temperature conditions in geothermal plants does not yet exist. Most often, scaling in geothermal plants is addressed with regular cleaning measures using high doses of aggressive chemicals with limited knowledge of the present degree of scaling [5].
In the oil and gas industries, offline methods commonly used to detect scaling include chemical analysis of the produced water or brine, and monitoring residual scaling inhibitor or the total suspended solids present. Whilst costeffective and readily available, these indirect methods are limited in their accuracy as samples taken are not necessarily representative of the overall state of the pipeline and can deteriorate by the time they reach the lab [3]. In-situ and online techniques such as the Thickness Shear Mode Res-onator (TSMR) or the fiber optic detection principle require sensors to be inserted into the pipe and then measure scaling growth locally on the instrumentation itself [6,7]. TSMR measures the change in resonance frequency of a piezoelectric mass sensor due to the mass of the adherent scaling. The fiber optic detection system measures radiation transmitted through an optical fiber with a partially exposed core, whereby a reduction thereof would imply scaling growth on the fiber surface. For geothermal power plants in particular, a noninvasive approach would be favourable as the high temperature and large volume flow rate would require very robust instrumentation and careful calibration leading to down time. Pipeline Inspection Gauges (PIGS) can be used for inline inspection of pipelines and can map their inside geometry to determine scaling thickness [8]. A PIG which is adapted to the high temperature of geothermal water would, however, present quite a challenge. Scaling can also be detected indirectly by calculating the fouling resistance via temperature and flow rate measurements [9] or monitoring the pressure drop at the inlet and outlet of heat exchangers [10]. These methods are, however, insensitive to submillimeter fouling layers, can have a considerable error (± 20%) and the measured phenomena are not exclusively caused by scaling, so results could be misinterpreted [3,9].
The aim of this research is to develop a practicable monitoring solution which optimises the administration of cleaning measures in geothermal plants in dependency of the amount of scaling present. To this end, two potential non-destructive techniques are investigated and discussed: UT and IE.
In standard contact UT, piezoelectric sensors are used to transmit ultrasound pulses into the test piece and record any reflected echoes. Two challenges for contact UT in this application are the high temperature and curved surface of geothermal pipes. Whilst there are a number of piezoelectric materials around [11], with sufficiently high Curie temperatures (beyond which they are depolarised, losing their piezoelectric properties) which could be tailor-made to fit a given curvature radius, the simpler and much less costly option is to place a temperature-resistant wedge between a standard contact transducer and the pipe surface. In a prototype version of IE testing, the response of the pipeline system to a manually applied impact using an impact hammer is measured using a lower frequency force sensor coupled directly to the pipe. The recorded signal is then analysed in the frequency domain [12].
The predominant type of scaling in the Bavarian Molasse basin, where the Sauerlach facilities are located, is calcite [2]. The research described in the following therefore focuses on this locally most common pipe-scaling combination of steelcalcite.

Experimental Setup
Scaling growth generally takes a considerable amount of time. For the purpose of our experiments, the time span in question, even using measures to accelerate the process, was unfeasible. A new 'top-down' rather than 'bottomup' approach was devised in which existing heavily scaled pipeline segments were gradually descaled using an acidic solution. During the descaling process, ultrasound and impact-echo measurements were carried out. A descaling test rig was designed with the project partners Z & H Wassertechnik, measX, the Chair of Hydrogeology and the Chair of Non-Destructive Testing, Technical University of Munich, assembled by Z & H Wassertechnik and modified by TUM ZfP. The test rig and the two measurement techniques applied in the descaling experiment are described in this section.

Descaling Test Rig
The descaling test rig was designed to gradually descale a piece of heavily scaled pipe by pumping an acidic solution through it. A scaled pipe segment is clamped between two 4 cm thick plastic plates, with an inlet in the lower plate and an outlet in the upper plate, standing in a chemical-resistant glass fibre reinforced polymer (GFRP) tub with an outlet leading to a pump. The acidic solution used to etch away the scaling was the cleaning agent RS 990 by Z & H Wassertechnik GmbH, St. Wendel, with active components 50% nitric acid and 85% phosphoric acid. 20 L of RS 990 were diluted with desalinated water to a 50% solution, filled into the GFRP tub and circulated through the pipe segment via the pump. The course of the acidic solution through the test rig is shown by the arrows in Fig. 1a.
The test piece was a 30 cm piece of DN200 black steel pipeline with a wall thickness of 11.5 mm from the production well at the geothermal power plant in Sauerlach, Germany. The pipe had grown a ∼ 6.6 mm-thick scaling layer over the course of 6 months. The rate of scaling growth is highly dependent on factors such as how often the submersible pump in the production well is switched off and on. A cross section of the pipe with the initial scaling layer is shown in Fig. 1b at the top. A microscope image of a fine-cut scaling sample is shown in Fig. 1b bottom. The applicability of the methods described to other types of scaling is currently under investigation.

Contact Ultrasonic Testing Setup
Two factors which need to be considered in coupling contact ultrasound transducers to geothermal pipelines are the high temperature (∼ 140 • C at the production well in Sauerlach) and the curvature of the surface. In the following measurement, a 2.7 cm-thick temperature-resistant polyphenylen- Fig. 1 a Descaling test rig. An acidic solution was pumped in the direction of the arrows shown, through the pipe segment, gradually etching away the scaling layer. b Top: Cross section of scaling layer in the pipe segment. Bottom: Microscope image of a fine-cut scaling sample sulphon (PPSU) wedge, milled to fit the curvature of the pipe in the test rig, was used. Heat transfer through the PPSU wedge was tested by placing it on a sand bath (for even temperature distribution) on a hotplate at 150 • C for 15 h.
A high temperature couplant (Echotrace HT, Karl Deutsch) was used between probe and wedge and between wedge and steel pipe. A constant contact pressure was ensured by a combination of a suitably designed and 3D printed holding bracket and an elastic band (Fig. 2). Preliminary experiments comparing the signals recorded on pipe segments with and without scaling for a range of frequencies (0.5 MHz, 1 MHz, 2.25 MHz, 5 MHz, 10 MHz and 15 MHz) showed that 5 MHz was a good compromise between ultrasound range and resolution.
The ultrasound instrumentation used was an Olympus Omniscan MX2 in impulse-echo mode with a 5 MHz, 6 mm diameter Olympus V201 RM contact transducer. The parameters were set to 100 V pulse excitation voltage, 49 dB gain, zero noise reduction, no averaging and no filter. Over the course of the descaling experiment a set of 100 A-scans was recorded every 3 min. Aligning such a set of A-scans to resemble a B-scan showed no relevant phase change in the signals over the recording period which was on the order of 2 s. A much lower number of A-scans, for example 10, would suffice for the same SNR in future measurements. Data analysis was carried out in Matlab: The 100 A-scans recorded at each time step were averaged, for improved signal to noise ratio (SNR), and the time value corresponding to the amplitude peak of the pipe backwall echo was subtracted from that of the scaling backwall echo giving the two-way traveltime (TWT) in the scaling itself. The results are detailed in Sect. 3.1.

Impact-Echo Setup
In addition to the US contact measurements, an impact-echo measurement was carried out at each 3 min time interval.
The pipe was impacted manually with a steel impact hammer, the resulting signal was recorded using a force sensor and the frequency response of the system was analysed. On the receiving end we used a TiePie USB oscilloscope, connected via a signal conditioner to a broadband force sensor (both PCB Piezotronics) with an upper frequency limit of 36 kHz. The TiePie multichannel oscilloscope software was used for data acquisition. The measurement parameters were set to a range of 200 mV, 2 MHz acquisition rate and 50 kS per recording, giving a frequency resolution of 20 Hz. In the first descaling measurement a ∅ 5 mm steel impact hammer was used for excitation. Because the peak resonant frequency was not consistent over the course of the descaling experiment we also measured with a ∅ 10 mm impact hammer during the second experiment. The diameter and therefore, mass of the tip influences the stimulated frequency spectrum, whereby a smaller diameter leads to a higher and broader frequency spectrum [13]. The ∅ 10 mm impact hammer diameter yielded promising results: A clear and continuous correlation between scaling thickness and the excited resonant frequencies is observed, whereby resonance frequencies decrease with scaling thickness in the pipe. The different frequency peaks can be attributed to the excited wave modes [14,15].
The scaling thickness was also measured mechanically, for reference, using a threaded rod with a pointer at one end, a rubber tip at the other end and a 360 • protractor. One full rotation of the threaded rod corresponds to a depth change of 1.5 mm. The concept is shown in Fig. 3. The scaling thickness at the beginning of the experiment was measured with a caliper.

Solution Analysis
For mass balance and calculating the scaling thickness via the calcium ion concentration, samples of the solution were taken after 3, 9 and 15 min, then every 9 min until 105 min, then every 30 min. The concentration of sodium, potassium, calcium and magnesium ions was determined by ion chromatography (IC; 881 Compact IC pro, Metrohm, Germany; column: C4-150, Metrohm, Germany). Before the measurement, the samples were diluted with ultra pure water (dilution factor of 500 and 1000 respectively) and filtrated (0.45 µm). The pH and EC of the solution were measured continuously with probes (Atlas Scientific, NY, USA) during the experiment.

Results of Contact Ultrasound Measurement
The first contact US descaling measurement yielded a linear plot of scaling thickness against TWT. Once the relevant peaks were identified, the shift in their position could be tracked over the course of the measurement, revealing the scaling thickness inside the pipe. Figure 4 shows the averaged A-scans for two points in the measurement time series: One was recorded near the beginning and one at the end of the first descaling measurement, corresponding to scaling thicknesses of 6.6 mm and 3.35 mm, respectively. The pipe and scaling backwall echoes are indicated.
The evolution of the TWT with scan time is shown in Fig. 5a and b for the first and second descaling measurement, respectively. The plateaus in the scatter plots shown are due to the limited time resolution of the instrumentation with the specified measurement parameters (0.05 µs). A higher resolution could be achieved using the described instrumentation by reducing the A-scan time frame. Whilst the pipe and scaling can easily be distinguished during the first measurement, the progressively thinner scaling layer thickness in the second descaling measurement makes superposition of the pipe and scaling backwall echoes a problem. A-scans in which the relevant backwall echoes could no longer be distinguished were removed in Fig. 5b. This included all Ascans taken after 69 min into the experiment, corresponding to a scaling thickness of ∼ 1.5 mm.
The travel time plotted against scaling thickness for the first and second measurements is shown in Fig. 6. The mechanically measured scaling thickness had a considerable error of ± 0.4 mm (± 0.16 mm due to acid corrosion of the threaded hole in the steel pipe over the course of the experiment, and a further ± 0.24 mm due to the uneven surface of the scaling, derived from the maximum measured difference between scaling crest and vale of 0.48 mm). An approximate longitudinal wave velocity in the scaling of 7330 m s −1 , with 95 % confidence bounds, can be derived from the gradient in Fig. 6. This is coherent with literature values for the p-wave velocity in calcite single crystals which range from 5700 m s −1 to ∼ 7700 m s −1 depending on the crystal orientation [16].
Given the linear correlation of TWT to scaling thickness, a thickness resolution of 0.18 mm could be attained with the current measurement parameters.
The preliminary experiment investigating the heat transfer through the wedge yielded a temperature of 60 • C at the upper

Results of Impact-Echo Measurement
Impact-echo signals for manual excitation with a ∅ 10 mm impact hammer, recorded at the beginning and end of the second descaling experiment are shown in Fig. 7a and c. 20 such signals were recorded at each time interval. These were detrended, autocorrelated, fast Fourier transformed and averaged to give frequency spectra such as those shown in Fig. 7b and d. The corresponding frequency spectra generated using the ∅ 5 mm impact hammer are shown in Fig. 8. The simplest scenario for automated evaluation of the scaling thickness from the IE data would be to track the peak resonant frequency over time. This requires an informed choice of hammer. Whilst excitation with a ∅ 5 mm impact hammer, for example, showed a clear decrease in resonant frequencies (Fig. 9a, b) with scan time, the peak resonant frequency was not consistent over the course of the descaling experiment.
Excitation with a 10 mm impact hammer, on the other hand, yielded an almost continuous curve for the peak resonant frequency over the course of the descaling process  Fig. 10). By implementing a moving window, which successively shifts with a particular peak value in the preceding frequency spectrum, the desired frequency can be tracked, provided it is continuous. This was done for Fig. 10. Here, the peak resonant frequency decreases with scaling thickness from an initial value of 16.6 kHz, for a (3.350 ± 0.004) mm thick scaling layer, to 15.28 kHz when no scaling is present. A linear trend is shown in the corresponding spectrogram plotted against scaling thickness (Fig. 11). With the given frequency resolution of 20 Hz, a peak frequency range of 15.28-16.6 kHz and a linear correlation between the two, a thickness resolution of (0.1 ± 0.1) mm could theoretically be reached for the scaling thickness, whereby a statistical error of 30 Hz is accounted for.
The frequency domain measured comprises modal oscillations of the pipe. As the scaling layer is etched away, the stiffness of the oscillating system decreases. A lower stiffness leads to lower resonant frequencies [18]. The mass also decreases with decreasing layer thickness, counteracting this effect, but the effect of the change in stiffness dominates [19].  Fig. 11 a Plotted against the mechanically measured scaling thickness, the frequency spectra show a linear trend between frequency and layer thickness. Higher peak frequencies show a larger shift in frequency with layer thickness. b shows an outtake around the peak frequency

Calculation of the Dissolved Scaling Layer Thickness
When the RS990 was added to the scaled test rig, the protons of the acid dissolved the CaCO 3 and the pH increased slowly (see Fig. 12a). Since the solution was not changed during one experiment, the capacity of the acid (= free available protons) for dissolving solid CaCO 3 decreased and the solution became saturated. This was accompanied by a continuous increase in pH value after approx. 60 min. The electrical conductivity decreased over time (see Fig. 12b) as the protons of the acids were consumed and the dissolved CaCO 3 has less influence on the conductivity than H + . For mass balance and the back calculation of the scaling thickness, IC measurements were carried out. The mass of dissolved CaCO 3 during the experiment was calculated using the concentration of Ca 2+ ions. Tables 1 and 2 show the evolution of the calcium and magnesium ion concentrations over time. Since the concentration of magnesium ions was negligible, the assumption of pure calcite scaled on the test rig is valid and in accordance with the findings of Köhl et al. [2]. Before the first experiment, a scaling layer of (6.62 ± 0.24) mm was measured by a sliding caliper. This means the amount of precipitated CaCO 3 on the test rig was (3.640 ± 0.018) kg (or: (1.460 ± 0.007) kg Ca 2+ ; simplification: CaCO 3 layer = hollow cylinder). After the first experiment, the scaling layer measured with a sliding caliper showed that (3.27 ± 0.24) mm of the CaCO 3 scaling were dissolved by the acidic solution. This equals to (1.770 ± 0.018) kg calcite (or (695.200 ± 7.327) g Ca 2+ ). The solution had a total volume of (40 ± 5) L, i.e. (17.380 ± 2.180) g L −1 Ca 2+ should be dissolved after the first experiment. In the second experiment, (3.35 ± 0.24) mm of the scaling layer were dissolved according to the sliding caliper, which is equivalent to (1.870 ± 0.018) g Ca 2+ ). calcite (or (734.440 ± 7.327) g Ca 2+ ) corresponding to a Ca 2+ concentration of (18.361 ± 2.302) g L −1 . Although the experiments were run in succession, they should be considered as independent experiments: The set-up was openend between the experiments and the solution was exchanged.
At the end of the first experiment the measured Ca 2+ concentration of the solution was (17.753 ± 0.065) g L −1 . It confirmed that half of the scaling layer was dissolved during the first experiment. At the end of the second experiment the measured Ca 2+ concentration was (15.784 ± 0.252) g L −1 . The back calculation of the scaling thickness showed that during the first experiment a scaling layer of (3.294 ± 0.412)mm was dissolved which is in really good agreement with the mechanically measured dissolved scaling (3.25 ± 0.40) mm, see Fig. 13a. At the end of the second experiment a scaling thickness of (2.858 ± 0.361) mm was dissolved according to the ionic measurements (mechanical measurement: (3.33 ± 0.40) mm). Both, mechanical measurement and chemical measurement show the same development over time for the first and for the second experiment (see Fig. 13a and b respectively). With regard of the uncertainties of the experiment and measurements there is no significant difference between the measurements. Nevertheless, the mechanical measurements give a higher value for the thickness compared to the chemical measurements (see Fig. 13b), because the probe is tipping on a rough surface (see Fig. 1b), thus capturing only the peaks of the precipitates. With increasing roughness, this deviation increases. The calculated scaling layer dissolved during the experiment is shown in Tables 1 and 2

Conclusion and Discussion
Whilst both measurement techniques show potential for insitu application, there are a number of factors which need to be considered. For contact US testing, the prerequisites are knowledge of the longitudinal wave velocity in the pipe material as well as in the scaling; a suitably milled wedge to optimise coupling to the pipe and protect the ultrasound probe from high temperatures and access to a suitable location for the measurement apparatus. Given this information, submillimeter resolution can be reached for scaling layers thick enough for pipe and scaling echoes to be resolved. A large gain (48 dB) was necessary to make the scaling signals visible. A more inhomogeneous type of scaling (with many scattering centres) would have a lower SNR and may be more difficult to detect. The nature of the bonding between pipe and scaling is also important: Scaling signals could not be made visible in the presence of pores or cavities between pipe and scaling. Reproducibility is a further point of query which needs to be addressed: Slight shifts in coupling can effect significant changes in signal amplitude and local disbonding between pipe wall and scaling may also falsify results. Although the scaling signal is clearly visible in the experiments described, it is not necessarily detectable for other measurement constellations [20]. As the critical components (pump, heat exchangers) are typically not directly accessible, one way to determine their degree of scaling would be to model the scaling development in relevant regions of the plant and extrapolate the values of interest from local US measurements made in accessible locations nearby. Testing should be carried out in several positions to reduce susceptibility to local deviation. For the measurement parameters described, superposition of the scaling and pipe signals for thin scaling layers limits detectability to layers ≥ 1.5 mm. Reducing the pulse width, for example by increasing the transducer frequency, using shear wave transducers or decreasing the number of cycles, could enable the detection of thinner layers. It should be noted, however, that a higher frequency would lead to increased scattering and a reduced penetration depth [17], the attenuation of shear waves is significantly higher than that of longitudinal waves [21], whilst a shorter pulse would reduce the energy introduced into the test piece, all leading to a lower SNR. The minimum detectable scaling layer would also vary for different types of scaling: A lower US velocity would lead to an increased resolution. Geothermal power plants in the Bavarian Molasse Basin mostly implement preventive measures such as reducing the temperature of the thermal water from the production well by some degrees or injecting carbon dioxide to counteract limescale [22,23]. Layer thicknesses of ≥ 1.5 mm are therefore rare on the production side of the plant. The injection well, on the other hand, can have scaling layers on the order of 1 cm, and could be an interesting monitoring location. Improved measureability might be achieved by comparison to simulation results or using machine learning algorithms but this requires further research.
Theoretically, an air-coupled ultrasonic testing setup would also withstand the high temperature of geothermal pipe surfaces. Due to the high impedance mismatch between air and steel, however, the signal amplitude introduced into the specimen would be negligible (∼ 0.001%) [24]. This would result in a very low SNR. Another interesting approach would be phased array ultrasonic testing with full matrix capture (FMC) acquisition and post-processing, for example using the total focusing method (TFM) for maximal resolution. This was also tested and has comparable limitations to single element contact UT, but is more expensive (financially as well as computationally) and the probe introduces less energy into the specimen: without TFM post-processing the scaling signal could not be observed at all.
The correlation of the obtained signals to scaling layer thickness could be validated by the back calculation of the dissolved scaling layer via the measured Ca 2+ ion concentration of the solution: the calculation only deviates by 1.35% for the first and by 14.26% for the second experiment from the mechanical measurement and shows the same development over time as the contact US and IE measurements. It makes sense that the deviation between the mechanical and chemical measurement of the second experiment exceeds that of the first experiment as the play in the threaded hole caused by acid corrosion increased over the course of the experiments. Furthermore, due to the probe tipping on the rough surface of the scaling, only the peaks of the precipitates are detected, which is why the mechanically measured scaling thickness should generally slightly exceed the one calculated via the Ca 2+ ion concentration of the solution (see also Figs. 1b, 13b). With increasing roughness, the deviation also increases.
A scaling monitoring using IE could resolve much thinner scaling layers. It does, however, need a calibration of the system, correlating a given pipe geometry and scaling type to the corresponding frequency spectrum. For in-situ measurements this will involve recording signals and creating a reference spectrogram for the geothermal power plant in question during one revision cycle. The frequency spectrum recorded at the beginning of such a cycle, so directly after cleaning, will correspond to a bare pipe with no scaling. The frequency spectrum at the end of the cycle can be attributed to the corresponding scaling thickness determined via mechanical measurement provided the pipe can be opened at this time. For future determination of scaling thickness the frequency spectrum at a given time can be correlated via linear regression to this reference spectrogram.
The length of the pipe segment studied in this work was chosen for practiacl reasons (weight and manoeuvrability). The resonance frequencies would be expected to decrease with increasing ratio of pipe length to outer radius, whereby the amount by which a resonance frequency decreases varies depending on vibration mode [19,25]. It can be seen in the spectrograms presented, that the majority of peaks exhibit a consistently decreasing trend with reduced scaling thickness. The occurence of the maximum peak, however, may be affected by a number of factors, such as the excited frequency range (dependent on contact time, diameter of impact hammer) [26], in-situ boundary conditions, impedance ratio of pipe to scaling and differences in attenuation depending on material properties and vibration mode [27]. Furthermore, lower frequencies shift less than higher frequencies with scaling layer thickness. The dependency of frequency shift on frequency is yet to be determined (a simulation to this end is currently in progress). Nevertheless, once a specific peak, or peak group is chosen as the indicator for scaling thickness, based on parametric studies, the described method can be effectively applied.
For optimal excitation (where a peak resonant frequency can simply be tracked over time), a suitable hammer size/ material and impact energy have to be established. Automation of the excitation using a solenoid-based impactor is currently being tested. Increasing the frequency resolution could further improve accuracy.
Laser excitation presents an alternative broadband, contactfree means of excitation [28,29], but comes with higher costs and extensive safety regulations.
Provided a calibration of the measuring apparatus can be carried out, UT and IE have a comparable resolution and are relatively inexpensive and simple to implement. The limiting factor for contact UT is the minimum scaling layer thickness which can be detected, the ultrasound properties of the scaling and optimal coupling for a usable SNR. That of IE testing is the calibration needed to evaluate the results. Combining data from contact UT and IE testing could improve resolution and the reliability of the diagnosis. The next steps towards a reliable scaling monitoring will involve a more detailed analysis of the effect of temperature and comparison to simulation results.
Data Availability Not applicable.

Declarations
Competing interests The authors declare that they have no conflict of interest.
Ethical Approval and Consent to Participate Not applicable.

Authors' Information Not applicable.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.