Optimization design of hump phenomenon of low specific speed centrifugal pump based on CFD and orthogonal test
HomeHome > Blog > Optimization design of hump phenomenon of low specific speed centrifugal pump based on CFD and orthogonal test

Optimization design of hump phenomenon of low specific speed centrifugal pump based on CFD and orthogonal test

Dec 02, 2023

Scientific Reports volume 12, Article number: 12121 (2022) Cite this article

1266 Accesses

Metrics details

Aiming to eliminate the hump phenomenon in a low specific speed centrifugal pump, its structural parameters were optimized using the computational fluid dynamics method. Based on the \(k - \varepsilon\) turbulence model, a 3D steady analysis of the internal flow field was carried out. The \(L_{9} \left( {3^{4} } \right)\) orthogonal table was established, and four structural parameters, including the impeller outlet diameter, impeller outlet width, number of blades, and blade outlet angle, were selected as influencing factors. Nine orthogonal test schemes were developed and the results were analyzed through the weight matrix analysis method, obtaining the weight of the selected factors on the test results. The optimal scheme was selected according to the weight and the weight matrix analysis results have shown that the impeller outlet width had the dominant influence on the head, shaft power, and efficiency. Further, the blade number was the main influencing factor for shaft power and efficiency. The centrifugal pump flow control test bench was built to carry out the numerical simulation and test all the prototype and optimization pump indexes. Through the external characteristic test, it can be seen that the \(\beta_{2} Z^{0.773}\) of the optimized pump is 87.889, which is 24.89% lower than that of the prototype pump, which effectively optimizes the hump phenomenon of the centrifugal pump. The experimental results have shown that in underrated working conditions, the working performance of the optimized pump was improved significantly. The head size was reduced by 1.424%, and the efficiency was increased by 7.896%. By optimizing the structural pump parameters, its jet-wake hydraulic loss was reduced, and the head curve hump phenomenon was effectively eliminated. All the performance indexes of the optimized pump were higher than those of the prototype, verifying both the accuracy and reliability of the orthogonal test and weight matrix analysis method. Finally, obtained results provide a reference for the structural design of high-performance centrifugal pumps.

A low specific speed centrifugal pump is a type of centrifugal pump with a specific speed ranging between 20 and 80. It is characterized by low flow, high head, and low volume, and is widely used in production and life1. Using a low specific speed centrifugal pump, it is easy to produce an unstable hump phenomenon under the low flow rate condition. This, in turn, increases vibration and noise, thus shortening the pump service life and reducing its working reliability. Currently, the hump phenomenon mechanism of low specific speed centrifugal pumps is not clear, and the hump phenomenon of the head curve cannot be eliminated from the design. Therefore, besides studying the centrifugal pump working mechanism, it is also necessary to optimize the critical pump structural parameters to improve its working performance.

For a long time, reducing and eliminating the hump phenomenon of centrifugal pump head curve with low specific speed becomes an important segment of centrifugal pump research. Zhang Desheng et al.2 established ten design schemes and performed numerical simulation and performance prediction on low specific speed centrifugal pumps. The authors obtained distributions of static pressure, streamlines, velocity, and turbulent kinetic energy in the pump, and improved the internal flow characteristics. Zhang et al.3 used the SAS turbulence model to carry out the full-channel numerical simulation of the pump-turbine, determining the effect of the pump flow structure mechanism on the hump characteristics. Zhang Peifang et al.4 analyzed the hydraulic performance of low specific speed centrifugal pumps, observing the cause of the head curve hump phenomenon aiming to propose a solution. Li et al.5 used the three-dimensional constant value simulation equation to design the pump blades. Further, they analyzed the impeller flow characteristics and found the pump energy discharge characteristics in the hump zone mode. Chen et al.6 experimented by adding two partitions in the pump suction section. Experimental results have shown that the proposed method could effectively improve the IS centrifugal pump performance curve and eliminate the head curve hump.

With the continuous development of computational fluid dynamics (CFD), numerical simulations have become one of the most important methods to study the internal flow field of centrifugal pumps7,8,9,10,11,12,13,14. In this paper, a subtype of the centrifugal pump was used as the research object, and the CFD method was used to calculate the internal flow field of the centrifugal pump. Accounting for the influence of various structural parameters on the pump performance, structural parameters were optimized via orthogonal tests. The centrifugal pump head, shaft power, and efficiency were calculated using range analysis. The optimal parameter combination was obtained through weight matrix analysis. A centrifugal pump flow control test bench was built to verify the design method. Finally, the experimental results have confirmed that the centrifugal pump performance indicators improved after carrying out the surface optimization; the head curve humps were eliminated, and the desired effect was achieved.

The structural parameters of a centrifugal pump are as follows: flow \(Q = 102{\text{m}}^{3} /h\), head \(H = 100{\text{m}}\), rotational speed \(n = 2900{\text{r}} /min\), and specific speed \(n_{s} = 56.343\). The semi-open impeller structure was adopted and the main structural parameters included: blade number \(Z = 7\), blade outlet angle \(\beta_{2} { = }26^\circ\), blade wrap angle \(\varphi { = }110^\circ\), impeller inlet \(D_{j} = 100{\text{mm}}\) , and outlet \(D_{2} = 290{\text{mm}}\) diameters, impeller outlet width \(b_{2} = 13{\text{mm}}\), and volute base circle diameter \(D_{3} = 302{\text{mm}}\).

The centrifugal pump calculation domain is shown in Fig. 1. The fluid flows into the centrifugal pump along the \(z\) axis, while the volute circumference is located in the \(xoy\) plane. The overall calculation domain includes both the impeller and volute water body.

Calculation domain.

During the operation, the centrifugal pump is affected by the impact loss. There is a difference between the theoretical \(H_{t} - Q_{t}\) and the actual \(H - Q\) centrifugal pump performance curves. According to theoretical analysis, theoretical centrifugal pump head \(H_{t}\) and flow \(Q_{t}\) satisfy the following equation:

where \(u_{2}\) is the circumferential blade outlet velocity calculated as:

Further, in Eq. (1), \(\sigma\) represents the Stoddard slip coefficient obtained as:

By sorting equations given above, theoretical centrifugal pump head \(H_{t}\) and flow \(Q_{t}\) can be obtained via:

According to Eq. (4), the theoretical centrifugal pump head \(H_{t}\) and flow \(Q_{t}\) have a linear relationship. The theoretical \(H_{t} - Q_{t}\) curve shows a gentle trend, while the actual \(H - Q\) curve is prone to the hump phenomenon15. When aiming to find the curve without a hump phenomenon, its slope should be increased. Using Eq. (4), impeller outlet diameter, impeller outlet width, the number of blades, and blade outlet angle were selected as critical factors affecting the hump phenomenon in the actual \(H - Q\) curve.

Under the small flow conditions, a low specific speed centrifugal pump easily generates vortex inside the flow passage which blocks the flow, causing large hydraulic loss and head reduction. This makes it easy to produce the hump phenomenon in the \(H - Q\) curve.

An orthogonal test is an analysis method for studying the influence of multiple factors and levels on test results by using an orthogonal table16. Using the orthogonality, representative combinations can be selected to experiment, and the influence of each parameter is analyzed to obtain the optimal combination. Thus, the orthogonal test is an efficient, fast, and economical method for designing experiments.

Based on the \(H_{t} - Q_{t}\) curve equation analysis and by considering the influence of various structural parameters, the impeller outlet diameter, its width, blade number, and blade outlet angle were selected as the experimental factors for the orthogonal test. The orthogonal test table was used to combine all the factors and carry out the numerical calculation, among which three horizontal values were set for each factor, as shown in Table 1.

According to \(L_{9} (3^{4} )\) the orthogonal test table, nine test scheme groups were set up (see Table 2).

The GAMBIT professional grid generation software was used to divide the computational domain and obtain high-quality hexahedral grids. Additionally, it was used to refine local grids for complex structures. Four groups of grids with rather various grid numbers were selected based on head to carry out the grid independence analysis and to ensure the reliability and accuracy of the calculation results. The results are shown in Table 3. When the number of grids exceeded 2.14 million, the calculated head changed slightly, proving that grids are not relevant. Thus, the final number of grids was 2.14 million.

Figures 2 and 3 show schematic diagrams of calculation domain grids.

Impeller grids.

Volute grids.

The internal flow field was calculated using Fluent software. The numerical solution was calculated based on the finite volume method of the completely unstructured grid. The standard \(k - \varepsilon\) model was selected as the turbulence model to solve the internal flow field. Moreover, the SIMPLE algorithm was used as the coupling method for internal pressure and velocity. The speed inlet was taken as the inlet condition, while the outlet boundary condition was free outflow. The non-slip wall was taken as the fixed-wall boundary condition, and its rotating counterpart was taken as the impeller wall. The convergence accuracy was 0.0001.

The CFD software was used to simulate the internal flow field of nine test combination groups. A mathematical model using head, shaft power, and efficiency as objective functions was established. Their respective values were calculated respectively using Eq. (5)17. The calculation results are shown in Table 4.

Range analysis was carried out using the orthogonal test results; the analysis results are shown in Table 5.

As shown in Table 5, when considering head, the main sequence of influencing indexes was \(BACD\), with the optimal combination as follows: \(D_{2} { = }300{\text{mm}}\), \(b_{2} { = }13{\text{mm}}\), \(z{ = }7\), and \(\beta_{2} = 26^\circ\). For shaft power, the main sequence of influencing indexes was \(BCAD\) and the optimal combination was \(D_{2} { = }290{\text{mm}}\), \(b_{2} { = }9{\text{mm}}\), \(z{ = }5\), and \(\beta_{2} = 24^\circ\). Finally, for the efficiency, influencing indexes of various factors were arranged as \(BCDA\), with combination \(D_{2} { = }300{\text{mm}}\), \(b_{2} { = }9{\text{mm}}\), \(z{ = }6\), and \(\beta_{2} = 24^\circ\) being optimal. Aiming to rapidly find the optimal combination and the main sequence of influencing indexes for each factor, weight matrix analysis was carried out. The weight matrix was obtained for three objective functions and structural parameters were selected according to the weight.

The multi-objective optimization of the weight matrix was carried out to establish a three-layer structure model based on the orthogonal test scheme18, as shown in Table 6.

The first layer represents the test objective function layer and is defined; assuming there are \(l\) influencing factors in the orthogonal test, each influencing factor has \(m\) levels. Moreover, the objective function value of the factor \(A_{i}\) at level \(j\) is \(k_{ij}\). The larger value of the orthogonal experiment objective function represents a better specimen; hence \(K_{ij} = k_{ij}\). If the objective function value is smaller, the better, then \(K_{ij} = 1/k_{ij}\). The matrix was established, as shown in Eq. (6).

The second layer represents all the factor layers and is defined as \(T_{i} = 1/\sum\nolimits_{j = 1}^{m} {K_{ij} }\). The matrix shown in Eq. (7) was established next:

The third layer is horizontal; the extreme differences in values of the orthogonal test factors were denoted as \(s_{i}\), which is defined as \(S_{ij} = s_{i} /\sum\nolimits_{i = 1}^{l} {s_{i} }\). Finally, the matrix shown in Eq. (8) can be written:

The weight matrix affects the objective function and is defined as: \(\omega = MTS\) and was established as:

In Eq. (9), \(\omega_{1} = K_{11} T_{1} S_{1}\), where \(K_{11} T_{1}\) is the ratio between the factor \(A_{1}\) target value at the first level to the target value of all the levels. Further, \(S_{1}\) is the ratio between the polar difference of factor \(A_{1}\) and the total polar difference. The product result reflects the influence of factors \(A_{1}\) on the objective function at the first level, along with the magnitude of extreme factor \(A_{1}\) differences. Through weight matrix analysis, the influence of each factor level on the objective function can be obtained. The primary and secondary order and the optimal combination of factor influences on the objective function can be quickly obtained through weight.

Using the above-presented expressions, weight matrices of three objective functions were calculated. For the head and efficiency, the higher objective function values represent the better solution, meaning that the corresponding values are \(K_{ij} = k_{ij}\), \(T_{i} = 1/\sum\nolimits_{j = 1}^{m} {K_{ij} }\), \(S_{ij} = s_{i} /\sum\nolimits_{i = 1}^{l} {s_{i} }\). For the shaft power, a lower objective function value is better. The corresponding values are \(K_{ij} = 1/k_{ij}\), \(T_{i} = 1/\sum\nolimits_{j = 1}^{m} {K_{ij} }\), \(S_{ij} = s_{i} /\sum\nolimits_{i = 1}^{l} {s_{i} }\).

The first objective function weight matrix (head) is as follows:

The second objective function weight matrix (shaft power) is as follows:

The third objective function weight matrix (efficiency) is:

Finally, the total weight matrix of the orthogonal test objective function is the average of the weight matrices of three objective functions:

According to the weight matrix calculation results, factor influences on the orthogonal test objective function are ordered as \(BCAD\); weights of the horizontal values for each factor are \(A_{3}\), \(B_{3}\), \(C_{2}\), and \(D_{3}\). The optimal orthogonal test combination scheme is \(A_{3} B_{3} C_{2} D_{3}\), namely: \(D_{2} { = }300\,{\text{mm}}\), \(b_{2} { = }9\,{\text{mm}}\), \(Z{ = }6\), and \(\beta_{2} = 22^\circ\).

To verify the optimization model feasibility, the full-flow CFD numerical simulation of the optimization model was carried out. The outputs were compared to results for the prototype pump.

Counter-clockwise numbering was used near the tongue flow passage (No. 1 flow passage). The internal pressure surface speed of each flow passage was greater than the suction surface speed. Additionally, the outlet speed reached the maximum value in the last flow passage. Due to the hydraulic impact between the flow passage and the tongue, the liquid flow direction changed. The relative velocity distribution at the intermediate centrifugal pump interface is shown in Fig. 4.

The relative velocity distribution at the intermediate interface of centrifugal pump.

When the small flow rate condition is active, the low-speed zone is mainly found on the blade pressure surface. Simultaneously, the liquid flow direction is disorderly between the blade inlet and outlet. Once the flow rate drops under \(0.8Q_{t}\), the degree of turbulence increases continuously. The vortices appear in multiple flow channels, and the number of vortices increases with the decrease in flow rate. Vortices within flow passages are continuously expanding in the clockwise impeller direction, blocking by low-speed liquid. This results in larger hydraulic losses, causing the prototype pump to produce a hump of centrifugal pump head curve under the condition of a small flow rate.

For the optimization pump, the degree of liquid flow direction disorder can be optimized above the \(0.4Q_{t}\) flow condition. In that case, the liquid flow direction in each impeller passage is more stable. Furthermore, under the \(0.4Q_{t}\) flow condition, a small number of vortices appear in the flow passage, similar to the prototype pump under the small flow condition. Thus, the head located at the closed dead center of the optimization pump is decreased.

A flow control test bench was built based on the NGL002 centrifugal pump test device to verify the accuracy of optimization results, as shown in Fig. 5. The prototype pump and the optimization pump were tested and verified for multi-condition external characteristics.

Flow control test bench.

The prototype and the optimization pump impellers are shown in Fig. 6.

Impeller.

The experimental results containing the prototype and the optimization pump external characteristics are shown in Figs. 7 and 8.

Head curve.

Efficiency curve.

As can be seen from the test results, the variation law of the external characteristic curve practically remains the same despite varying working conditions. The numerical simulation results are in good agreement with the experimental values, reflecting the changing trend of all indexes under different flow rates. The head and efficiency obtained via numerical simulation are higher than the experimental values, mostly because the energy loss of each part and the impeller manufacturing errors were not considered in the numerical simulation. Using the verification results obtained for external characteristics, it can be seen that the optimization pump performance index is higher than that of the prototype. The \(\beta_{2} Z^{0.773}\) prototype pump and optimization pump are 117.013 and 87.889, respectively. Since the optimized blade outlet angle and number of blades are reduced, the hump phenomenon was effectively optimized19.

The external characteristic test results were extracted for various indexes, as shown in Table 7.

Based on Table 7, it is evident that under rated conditions, the head of the optimized centrifugal pump was reduced by 1.424%, while the efficiency increased by 7.896%. As the prototype pump produced a hump under \(0.8Q_{t}\) working conditions, the head under \(0.6Q_{t} - 1.0Q_{t}\) working conditions was greater than the optimization pump head. The actual optimization pump \(H - Q\) curve does not show a hump phenomenon, meaning that the optimization goal is achieved: the index values are effectively improved. To sum up, the orthogonal test and weight matrix analysis method are feasible, with the optimized design scheme accuracy being verified.

Hydraulic design of low specific speed centrifugal pump was carried out. Three-dimensional impeller and volute models were created.

Nine groups of test schemes were designed using the orthogonal test, and the influence order of each factor on each of the indexes was obtained through range analysis. Further, by using the weight matrix analysis and the weight relationship between the factor levels, a set of optimization models were obtained: \(D_{2} { = }300{\text{mm}}\), \(b_{2} { = }9{\text{mm}}\), \(Z{ = }6\), and \(\beta_{2} = 22^\circ\). The internal flow field of the prototype and the optimization pump was numerically simulated. The simulation results have shown that the hump phenomenon was significantly improved in the optimization pump. Hence, the feasibility of the optimization scheme was verified.

The centrifugal pump flow control test bench was built and the simulation and test values of each prototype and optimization pump index were obtained under different working conditions. Test results have shown that the performance indexes were higher in the optimization pump, eliminating the hump phenomenon and improving its hydraulic performance. Accuracies of the design process and optimization method were further verified.

In this paper, a type of centrifugal pump was taken as the research object and hydraulic design was carried out based on the CFD. The main structural parameters were optimized via an orthogonal test. A test platform was built to enable carrying out verification tests, which verified the orthogonal test scheme accuracy and improving the working performance of the centrifugal pump. Compared to the prototype pump, the outlet impeller diameter of the optimized pump was increased, while the impeller outlet width, the number of blades, and the blade outlet angle were smaller. With the increase in the impeller outlet diameter and the decrease of its width, the flow passage area increased, reducing the hydraulic losses. In that case, the overload would be avoided, and the "jet-wake" phenomenon and hump would be eliminated. At the same time, with the decrease in the number of blades and the blade outlet angle, the area of each flow channel increased. Thus, the static pressure at the blade outlet decreased, mostly to avoid local backflow and flow separation. It can be seen that, by optimizing the structural centrifugal pump parameters, the pump head was reduced, and the internal flow stability and working efficiency were improved. Thus, a series of negative problems such as vibration, noise, and high energy consumption of pump and pipeline systems caused by the hump phenomenon could be mitigated effectively. In the next step, multi-objective optimization design will be carried out around the external characteristics of centrifugal pump head, shaft power and efficiency, focusing on in-depth research on centrifugal pump flow-induced vibration noise, in order to comprehensively improve the working performance of centrifugal pump.

Ma, H., Ding, R. & Yang, D. CFD research on hump phenomenon of low specific speed centrifugal pump. Fluid Mach. 41(12), 43–47 (2013) ((in chinese)).

Google Scholar

Zhang, D. et al. Turbulence analysis and experiments of low-specific-speed centrifugal pump. Trans. CSAE 26(11), 108–113 (2010) ((in chinese)).

Google Scholar

Zhang, C., Xia, L. & Diao, W. Influence of flow structures evolution on hump characteristics of a model pump-turbine in pump mode. J. Zhejiang Univ. (Eng. Sci.) 51(11), 162–171 (2017) ((in chinese)).

Google Scholar

Zhang, P., Yun, S. & Yu, J. Present status and development trends of low-specific-speed centrifugal pumps. J. Drain. Irrig. Mach. Eng. 21(6), 42–45 (2003) ((in chinese)).

CAS Google Scholar

Li, D. Y. et al. Fluid flow analysis of drooping phenomena in pump mode for a given guide vane setting of a pump-turbine model. J. Zhejiang Univ. Sci. A 16(11), 851–863 (2015).

Article CAS Google Scholar

Yan, C. et al. Experiment investigation into measures for eliminating hump in head curve of IS type centrifugal pumps. J. Drain. Irrig. Mach. Eng. 38(1), 34–37 (2020).

Google Scholar

Pei, J., Wang, W. J. & Yuan, S. Q. Statistical analysis of pressure fluctuations during unsteady flow for low-specific-speed centrifugal pumps. J. of Central South Univ. 21(3), 1017–1024 (2014).

Article CAS Google Scholar

Suqing, W. Three-dimensional design method of induction wheel and its performance verification. General Mach. 1, 61–64 (2018) ((in chinese)).

Google Scholar

Ni, X. D. et al. Improved similarity criterion for seepage erosion using mesoscopic coupled PFC–CFD model. J. Central South Univ. 22(8), 3069–3078 (2015).

Article Google Scholar

Chen, Z. et al. A new correlation for predicting the thermal conductivity of liquid refrigerants. J. Therm. Anal. Calorim. 143, 795–800 (2021).

Article CAS Google Scholar

Murlidhar, B. R. et al. The effects of particle swarm optimisation and genetic algorithm on ANN results in predicting pile bearing capacity. Int. J. Hydromech. 3(1), 69–87 (2020).

Article Google Scholar

Jiang, D. et al. A fully convolutional neural network-based regression approach for effective chemical composition analysis using near-infrared spectroscopy in cloud. J. Artif. Intell. Technol. 1, 74–82 (2021).

Article Google Scholar

Canan, K. & Burak, M. Optimisation design and operation parameters of a photovoltaic thermal system integrated with natural zeolite. Int. J. Hydromech. 3(2), 128–139 (2020).

Article Google Scholar

Al-Rasheda, A. A. A. A. et al. Finite volume simulation of mixed convection in an inclined lid-driven cavity filled with nanofluids: effects of a hot elliptical centric cylinder, cavity angle and volume fraction of nanoparticles. Phys. A 121122, 1–35 (2019).

MathSciNet Google Scholar

Xingfan, G. Theory and Design of Pump (Aerospace Press, 2011) ((in chinese)).

Google Scholar

Yuqin, W. et al. Optimization design of centrifugal pump cavitation performance based on orthogonal test. J. Shaoyang Univ. (Nat. Sci. Edit.) 15(3), 26–33 (2018) ((in chinese)).

Google Scholar

Houlin, L. & Minggao, T. Modern Design Methods for Centrifugal Pumps (China Machine Press, 2013) ((in chinese)).

Google Scholar

Xiaoling, W., Bingjun, X. & Qiang, Z. Optimization design of the stability for the plunger assembly of oil pumps based on multi-target orthogonal test design. J. Hebei Univ. Eng. (Nat. Sci. Edit.) 27(3), 95–99 (2010) ((in chinese)).

Google Scholar

Fu, L. & Hong, Z. Optimization design of lower speed pump based on genetic algorithm. J. Chin. Agric. Mech. 37(2), 233–236 (2016) ((in chinese)).

Google Scholar

Download references

This project is supported by Anhui Province University Discipline (Professional) Top Talent Academic Funding Project (No. gxbjZD2021076). This project is supported by Key Project of Natural Science Research in Colleges and Universities of Anhui Province (No. KJ2021A1026). This project is supported by Key Project of Natural Science Foundation of Chaohu University (No. XLZ-201902). This project was funded by National College Student Innovation and Entrepreneurship Training Program Project (No. 202210380038). This project is supported by Anhui Province Excellent Talent Cultivation Innovation Project (No. 2018zygc031).

School of Mechanical Engineering of Chaohu University, Anhui Chaohu, 238000, China

Wang Yu-qin & Ding Ze-wen

Technological University of the Philippines, 1106, Manila, Philippines

Wang Yu-qin

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

Z.D. is responsible for the numerical simulation of the centrifugal pump and the calculation of the weight matrix. Y.W. puts forward the overall design ideas and orthogonal test schemes, and is responsible for test verification and other work.

Correspondence to Wang Yu-qin.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Yu-qin, W., Ze-wen, D. Optimization design of hump phenomenon of low specific speed centrifugal pump based on CFD and orthogonal test. Sci Rep 12, 12121 (2022). https://doi.org/10.1038/s41598-022-16430-w

Download citation

Received: 27 July 2021

Accepted: 11 July 2022

Published: 15 July 2022

DOI: https://doi.org/10.1038/s41598-022-16430-w

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.