Abstract

Geometrically, the single-blade centrifugal impeller, commonly used today as a sewage pump, is not axially symmetric. For this reason, the static pressure around the impeller fluctuates greatly when the impeller is rotating, and not only the radial thrust but also the axial thrust shows large fluctuations. Therefore, it is extremely important for the improvement of pump reliability to quantitatively grasp these fluctuating hydrodynamic forces. In this study, we investigated the unsteady hydrodynamic forces in a closed-type centrifugal pump with a single blade for different blade outlet angles using a numerical analysis that takes into account both experiment and the leakage flow. The results clearly showed the effect of the blade outlet angle on that act on the impeller. The root-mean-square value of the fluctuating component of the total radial thrust was roughly the same for whichever impeller at low flow rate, but at high flow rates, the value increased for impellers with larger blade outlet angles. Moreover, when the leakage flow rate increased with increasing static pressure around the impeller, such that the rear and front shroud parts were subject to high pressure, the absolute value of the axial thrust on both these parts increased.

1. Introduction

Because a sewage pump is used for the transfer of liquids containing foreign bodies, such as solid and fibrous material, it must not only be highly efficient but also allow the easy passage of foreign bodies. The passed particle size relative to the pump bore (the minimum particle size in the flow path inside the pump) exists as an index of the ease with which these foreign bodies can pass through [1]. It therefore follows that the single-blade centrifugal pump, consisting of just a single blade of broad blade width, has been widely used for this type of pump so that a large passed particle size is formed. Nevertheless, even for centrifugal pumps with many blades, geometrical manufacturing errors cause disproportional fluid force [2]. In contrast, with the single-blade centrifugal pump, by its very nature, the impeller itself is not axially symmetric from a geometrical viewpoint. Additionally, because a volute casing is installed on the outside of the impeller that is not axially symmetric, the static pressure distribution surrounding the impeller is extremely nonuniform, and an excessively fluctuating radial thrust acts on the impeller [3, 4]. Because such excessively fluctuating radial thrust can lead to a drop in the axial strength, a shortening of the lifetime of the bearings and shaft sealing equipment, and the generation of axial vibrations, elucidating the generation mechanism behind this has become an important challenge.

Regarding the radial thrust in single-blade centrifugal pumps, such subjects have been investigated as the effect of the number of blades [3, 5], winding angle of the blade [3], distribution of the blade angle [6] and its relationship with the static pressure between the blades [4], effect of the casing shape [3, 6], effectiveness of adding extra mass [4, 7], and relationship to vibration of the casing [7]. In addition, the relationship between the hydrodynamic forces acting on the impeller and the whirling of the impeller has been investigated [5, 810]. We have also proposed a design method for single-blade centrifugal impellers and have discovered, for a closed-type centrifugal impeller with a single blade that was designed using this design method as a basis, how the performance characteristics relate to the radial thrust [11]. While modeling the components of the radial thrust using an unsteady conservation-of-momentum law to illustrate their behavior [12], we also explained how they were affected by the blade outlet angle [13]. However, hardly any research has been made on the axial thrust. In our previous studies, the hydrodynamic forces acting on the rear and front shroud parts were neglected in the numerical analyses. Accordingly, by performing a numerical analysis that takes into account the leakage flow on the rear surface of the rear and front shroud parts, we were able to demonstrate that not only did the radial thrust fluctuate with fluctuations in the static pressure distribution around the impeller but also the axial thrust fluctuated greatly [14]. Furthermore, modeling the components of the axial thrust, we showed how they behaved [15].

In this study, we implemented a numerical analysis that takes into account the leakage flow on the rear surface of the rear and front shroud parts of a closed-type single-blade centrifugal pump for different blade outlet angles. And with a comparative investigation, we could clarify how the blade outlet angle affects both the radial and axial thrust in the closed-type single-blade centrifugal pump.

2. Experimental Apparatus and Method

In this study, we used three types of single-blade centrifugal impellers with respective blade outlet angles of 8°, 16°, and 24°. However, in the experiments, we used only impellers with = 8 and = 16. An overview of the test impeller ( = 8) is shown in Figure 1, and its specifications are given in Table 1. An overview of the test pump is also given in Figure 2. The impeller with = 8 had a passed particle size of 56 mm and had been designed on the basis of the design method reported previously [11]. This design method stipulates the passed particle size beforehand, and the impeller design is implemented using, as design parameters, first, the relative velocity ratio of the blade inlet and blade outlet and, second, the absolute velocity ratio of the impeller suction inlet and the blade inlet. The impellers with = 16 and = 24 have meridian planes that are equivalent to that for the impeller with = 8, and, as shown in Table 1, only the blade outlet angle is varied. This pump has three narrow gaps between the impeller and the casing, as seen in Figure 2. These are the liner-ring gap (the gap with the casing wearing ring), the gap with the rear shroud part, and the gap with the front shroud part. The gap widths are, respectively, 0.75, 0.5, and 0.5 mm and the corresponding gap lengths are 17, 10.1, and 10.5 mm. It should be noted that, because a mechanical tool has been adopted for the shaft sealing equipment in this pump, almost no leakage to the exterior occurs on the rear shroud part side.

Figure 3 gives an overview of the experimental apparatus. The experiment was conducted at a rotational speed of  min−1, and the flow rate was adjusted with the gate valve mounted on the discharge side and measured using an electromagnetic flow rate meter. The total head was found from the static pressure measured before and after the pump with a strain-gauge-type pressure transducer and flow rate. The shaft power was found by measuring the rotational speed and torque with a magnetoelectric-type rotation detector and torque detector, respectively.

The coordinate system is defined by the -, -, and -axes, as shown in Figure 2. The blade phase angle is taken in the anticlockwise direction from the -axis, and the point where the end of the blade outlet passes through the -axis is assumed as = 0° (the position shown in Figure 2). Eight monitor points were provided on the wall surface on the side of the front shroud part at positions 10 mm radially outwards from the impeller outlet in an anticlockwise direction from the -axis as = 0°. The pressure measurement was synchronized with the rotation of the impeller using the torque signal for the = 0° position of the end of the blade outlet as a reference.

For the radial thrust measurement, an - load cell was incorporated in the periphery of the bearings, as shown in Figures 4 and 5, and it was synchronized with the rotation of the impeller in the same way as that for the pressure measurement [1115]. The “four-active-gauge” method was used for the - load cell by attaching the connected strain gauges together. The dynamic calibration method [16] was adopted for calibration. The relationship between the output voltage and the force was tested in air for a known centrifugal force acting on a calibrating circular plate, which has a mass found by subtracting from the impeller mass; that mass corresponds to the buoyancy force that normally acts on it. The calibration curves for the - and -directions when the force acts in the -direction are shown in Figure 6. The output voltage and the force are proportionally related. Moreover, when the force acts in the -direction, its influence on the -direction is small. It should be noted that the dynamic balance was taken between the impeller under test and calibrating circular plate. The effect of the heat from the bearings was eliminated using a method where the pump was stopped over a short interval of around 15 s under thermally stable conditions and was operated immediately afterwards, adjusting the zero point to take measurements under the same conditions.

3. Numerical Analysis Method

For the numerical analysis, we performed a three-dimensional unsteady-flow analysis using the general-purpose thermofluid analysis code ANSYS CFX13.0. Reynolds-averaged Navier-Stokes equations were used as a basis, and the SST (shear stress transport) model was used to model the turbulent flow.

The computational domain and computational grid for the pump part are shown, respectively, in Figures 7(a) and 7(b). This study differs from past research [1113] in that the gaps between the outer surfaces of the rear and front shroud parts of the impeller and the casing are actually modeled. The computation domain consists of the suction pipe, impeller, volute casing, regions in the gaps with the rear and front shroud parts, and discharge pipe. The computational grids for the impellers with = 8, = 16, and = 24 have 690,000, 705,000, and 707,000 elements, respectively, where the total number of elements in the grids is 2,328,000, 2,343,000, and 2,345,000, respectively. Here, the computational grids for the suction pipe, volute casing, regions in the gaps with the rear and front shroud parts, and the discharge pipe are identical in all cases, with, respectively, 149,000, 462,000, 340,000, 572,000, and 115,000 elements. The number of elements in the gaps counted radially is five for the liner-ring gap and four for both the rear shroud part and front shroud part gaps. A nonstructured tetrahedral grid is used as the computational grid for the impeller and volute casing, whereas a nonstructured hexahedral grid is used for the others.

To study the dependence of the computational model on the grid used, we performed analyses for the = 8 case, where there were, respectively, 8, 6, and 6 elements in the radial direction in the liner-ring gap, rear shroud part gap, and front shroud part gap, and computational grids for the rear shroud part gap and front shroud part gap regions had 561,000 and 1,017,000 elements, respectively, with a total grid number of approximately 2,994,000 elements. The results showed that the errors in the flow rate, at the maximum efficiency point in the analysis results for the = 8 case, were 0.2% for the head, 0.3% for the shaft power, 0.7% for the root-mean-square value of the total radial thrust, and 4.1% for the time-averaged value of the total axial thrust. For an impeller with = 8, we performed analyses with volute casing computational grid numbers of 1,035,000 and 693,000 elements, with a total grid number of approximately 2,904,000 elements. The results showed that the errors in the flow rate, at the maximum efficiency point in the analysis results for the = 8 case, were 1.8% for the head, 1.6% for the shaft power, 0.6% for the root-mean-square value of the total radial thrust, and 3.7% for the time-averaged value of the total axial thrust. As mentioned above, we verified that the grid dependence for this computational model was comparatively small [15].

The boundary conditions for the suction pipe inlet were set as those on the inlet boundary and for the discharge pipe as those on the outlet boundary and were given as the mass flow rate on the inlet boundary, a static pressure (gauge pressure) on the outlet boundary of 0 Pa, and the rotational speed in the impeller domain. The wall surface was assumed to be a nonslip condition, and the boundary between the rotating and static domains was conjugated using the transient rotor-stator technique [17]. The time step was set such that 360 steps were a single rotation of the impeller, and the calculations were continued until the fluctuations in the flow became practically constant.

4. Experimental and Analysis Results and Considerations

4.1. Comparison of Pump Performances

The pump performance for each different impeller based on both experiment and numerical analysis is shown in Figures 8(a)8(c). Qualitatively, the experimental and calculated values for the cases = 8 and = 16 showed good agreement for the head coefficient , shaft power coefficient , and pump efficiency . The calculated values of are seen to increase in impellers for which the increases with increasing flow rate. With the flow rate at the point of maximum efficiency, = 0.030, for = 16 and = 24 increased by 13.7% and 33.8%, respectively, relative to = 8. This could be attributed to the increase in the circumferential component of the absolute velocity of the blade outlet as the blade outlet angle increases, whereas, for the calculated values of , the flow rate at the point of maximum efficiency with = 8 is = 0.030 and that for = 16 and = 24 migrated to the larger side is = 0.039. The for = 16 and = 24 at = 0.039 was 0.735 and 0.767, respectively, and we find that increases for impellers with larger . The flow rate at the point of maximum efficiency is thought to be moving toward a larger flow rate because the theoretical pump head for a finite number of blades increases as the blade outlet angle increases, and the matching point with Worster’s characteristic formula for the volute [18] then moves toward a larger flow rate.

4.2. Comparison of Radial Thrust

In this study, the impeller surface is divided into three regions, as shown in Figure 9. The blade surface is assumed to be the hydraulic part, the outer surface of the rear shroud part including its outer periphery as the rear shroud part, and the outer surface of the front shroud part including its outer periphery as the front shroud part. The total radial thrust acting on the impeller is expressed by the following formula [14, 15]: Here, is the radial thrust on the hydraulic part [N], is the radial thrust on the rear shroud part [N], and is the radial thrust on the front shroud part [N].

The radial thrust can be defined by the following formula as the vector sum of the time-averaged value (an average component) and a fluctuating component [3, 4, 16]: Here, is the time-averaged value and is the fluctuating component.

The time-averaged value and the root-mean-square value of fluctuating component of the total radial thrust coefficient on the impeller from both experiment and numerical analysis are shown in Figure 10. The calculated value for radial thrust was found by integrating, over each of the impeller parts, the pressure acting on the wall surface and the shear stress on the wall surface. From Figure 10, the experimental and calculated values for and the root-mean-square value of are seen to agree qualitatively. Although all the impellers show roughly the same values for the root-mean-square value of with a small flow rate of = 0.008, at large flow rates, it increases for impellers with larger . This is thought to be because, although it also depends on the load distribution on the blade, for large flow rates, there is a relative increase in the load on the blade for impellers with larger . On the other hand, with , it increases for either a low flow rate or a high flow rate for whichever of the impellers, and there are comparatively small differences due to changes in . This is believed to be due to the circumferential distribution of the time-averaged pressure coefficient at the monitor points in the impeller outlet shown in Figure 11, which are the results of numerical analysis; gets higher for impellers with larger for any of the flow rates, and that value increases uniformly. Thus, we have confirmed that there is good agreement between the experimental and calculated values of at the impeller outlet [15].

Subsequently, we investigate the behavior of the radial thrust acting on each part of the different impellers. The Lissajous figures for the radial thrust coefficient for the hydraulic part, rear shroud part, and front shroud part of each impeller are shown in Figures 12(a)12(c) for = 0.030. For each radial thrust coefficient, the values at = 0, 90, 180, and 270 are shown, respectively, by ●, ▲, ■, and ◆, and the time-averaged values are given using the same symbols as those for the radial thrust coefficient. These are represented similarly in Figures 13(a) and 13(b). From Figures 12(a)12(c), the time-averaged value of the radial thrust acting on each part is located near the origin for any of the impellers. The fluctuating component of the radial thrust for both the hydraulic and rear shroud parts increases for impellers with larger . However, in the case of radial thrust for the front shroud part, it is smaller for the impeller with = 24 than for the one with = 16. Because a leakage flow exists on the side of the front shroud part, the effect of an increase in the leakage flow rate for the impeller with = 24 (discussed later) can reasonably be considered as the reason for this finding. Moreover, we find that the radial thrust for the rear shroud and front shroud parts is small in relation to the radial thrust for the hydraulic part for all the impellers. It therefore follows that the radial thrust for the hydraulic part dominates the total radial thrust.

Thus, concentrating on the radial thrust coefficient for the hydraulic part, its Lissajous figure is shown in Figures 13(a) and 13(b) for the low flow rate = 0.008 and high flow rate = 0.047. With the low flow rate = 0.008, the radial thrust for the hydraulic part is large in the vicinity of = 180°, and it is largest for the impeller with = 24. On the contrary, the radial thrust for the hydraulic part is small in the vicinity of = 0°, and it is smallest for the impeller with = 24.

On the other hand, with the high flow rate = 0.047, the radial thrust for the hydraulic part is large in the vicinity of = 0°, and it is largest for the impeller with = 24. Conversely, the radial thrust for the hydraulic part is small in the vicinity of = 180° and is smallest for the impeller with = 8. As shown with the numerical analyses in Figures 14(a)14(d), 15(a)15(d), 16(a)16(d), and 17(a)17(d), from the static pressure distribution on the impellers with = 8 and = 24, at both = 0.008 and = 0.047, the following observations can be made. The cross-section is in the center of the blade width. At the low flow rate = 0.008, the time-averaged increases as the end of the volute winding is approached. Moreover, if = 180°, because for a pressure surface with 180° < < 360° increases more, the radial thrust for the hydraulic part acts in the -axis direction. For a large blade outlet angle of = 24, both the for the pressure surface and time-averaged at the end of the volute winding are the largest, and the radial thrust for the hydraulic part in that case is the largest. However, if = 0°, the for a pressure surface with 0° < < 180° is increased, and the for = 24 is the largest. However, owing to the effect of the time-averaged , which increases as the end of the volute winding is approached, the time-averaged value of the radial thrust for the hydraulic part acts largely in the +-axis direction. Therefore, the radial thrust for the hydraulic part can be considered to act in the –-axis direction in terms of its time-averaged value. For = 24, the time-averaged is the largest, and the radial thrust for the hydraulic part is the smallest. On the other hand, with a flow rate = 0.047, the time-averaged increases as the beginning of the volute winding is approached. Moreover, if = 0°, because the for a pressure surface with 0° < < 180° increases more, the radial thrust for the hydraulic part acts in the direction. For = 24, both the for the pressure surface and time-averaged at the beginning of the volute winding are the largest, and the radial thrust for the hydraulic part is the largest. On the contrary, if = 180°, for a pressure surface with 180° < < 360° is increased. However, due to the effect of the time-averaged , which decreases as the end of the volute winding is approached, the time-averaged value of the radial thrust for the hydraulic part acts largely in the −-axis direction. Therefore, the radial thrust for the hydraulic part can be considered to act in the +-axis direction in terms of its time-averaged value. For = 8, both the for the pressure surface and time-averaged at the end of the volute winding are the smallest, and the radial thrust for the hydraulic part can then be considered to have its smallest value.

4.3. Comparison of Axial Thrust

Similar to the radial thrust in the previous section, the total axial thrust acting on the impeller is expressed by the following formula [14, 15]: Here, is the axial thrust for the hydraulic part [N], is the axial thrust for the rear shroud part [N], and is the axial thrust for the front shroud part [N].

The time-averaged values of the axial thrust coefficient for the hydraulic part, rear shroud part, and front shroud part are shown for each of the impellers in Figures 18(a) and 18(b). The calculated value for the axial thrust was found in a similar manner to that for the radial thrust, where the pressure acting on the wall surfaces and the shear stress on the wall surfaces is integrated overall into each part of the impeller. Comparing the absolute values of the axial thrust coefficient for each part and for all the impellers, the axial thrust coefficient for the hydraulic part is the smallest. The axial thrust coefficients for the rear shroud part and front shroud part are dominant, unlike the cases for the radial thrust. and are in opposing directions relative to each other, and their absolute values increase for smaller flow rates. Moreover, is larger than , from which it follows that the total axial thrust coefficient acts in the −-axis direction and its absolute value increases for smaller flow rates.

Comparing the individual impellers, the absolute value of the total axial thrust coefficient is the largest for = 24. Little difference exists between = 8 and = 16; however, for small flow rates, the difference is slightly larger for = 8. The larger the for the impeller is, the more the absolute values of and increase at high flow rates. However, when the blade outlet angle is varied from 8° to 16°, the increase in the absolute value of is relatively small compared with the increase in the absolute value of , and for = 0.008 in particular, the absolute values of for = 8 and = 16 approximately agree.

Concentrating on the total axial thrust and the axial thrust on the rear and front shroud parts for = 0.008 and = 0.047, the fluctuation of , , and during a single rotation of the impeller is shown in Figures 19(a)19(c). From Figure 19(a), the axial thrust coefficient for all the impellers shows a similar behavior and fluctuates greatly during a single rotation of the impeller. The absolute value of is largest in the vicinity of = 90°–180° for a low flow rate = 0.008, but conversely, for a high flow rate = 0.047, the largest value is at = 0°(360°). From Figures 19(b) and 19(c), the absolute values of and also behave in a manner similar to the absolute value of . Therefore, for all the impellers, the total axial thrust is partially offset by the axial thrust of the front shroud part, but we see a behavior roughly similar to the axial thrust for the rear shroud part.

For impellers with = 8 and = 24, we investigated the axial thrust for the rear and front shroud parts in relation to the pressure field. For a low flow rate = 0.008, the static pressure distribution in the clearance domain for the rear and front shrouds for the = 8 case is shown in Figures 20(a)20(d) and 21(a)21(d), respectively, whereas the same pressure distributions for the case = 24 are shown in Figures 22(a)22(d) and 23(a)23(d), respectively. Furthermore, the same pressure distributions are shown for a high flow rate = 0.047 in Figures 24(a)24(d) and 25(a)25(d) for = 8 and in Figures 26(a)26(d) and 27(a)27(d) for = 24. Here, the static pressure distribution in the clearance domain for the front and rear shrouds is for the cross-sections C-C and D-D, respectively, of Figure 3. For any of the impellers, the in the clearance domain of either the front or rear shroud fluctuates greatly during rotation of the impeller in the same way as the static pressure distribution in the interior of the pump, as shown in Figures 1417. For a low flow rate = 0.008, the for the clearance domain of the front shroud takes a high value for all the impellers in the region = 90°–180°. Owing to the effect of leakage flow, the for the clearance domain of the front shroud drops overall below that for the clearance domain of the rear shroud but similarly takes a high value in the region = 90°–180°, which corresponds to the behavior of the axial thrust for the rear and front shroud parts, shown, respectively, in Figures 19(b) and 19(c). At this point, if we examine the static pressure distribution in the interior of the pump, as shown in Figures 14 and 15, the static pressure around the impeller is greater overall in the region = 90°–180° than in the region = 180°–270°. That is to say, we can reason that when the static pressure around the impeller increases overall, the static pressure in the clearance domain for the front and rear shrouds increases in accordance with the increase in the leakage flow rate. This is evident from the fact that if we consider the fluctuation in the leakage flow rate coefficient from the liner ring as shown in Figure 28(a) with = 0.008, for any of the impellers attains its maximum value in the vicinity of = 90°–180°, which is larger; the larger value is the for the impeller.

On the other hand, for a large flow rate = 0.047, the for the clearance domain of the front and rear shrouds takes a high value at = 0°(360°), which corresponds to the behavior of the axial thrust for the rear and front shroud parts, as shown in Figures 19(b) and 19(c). At this point, looking at the static pressure distribution in the interior of the pump as shown in Figures 16 and 17, the static pressure around the impeller at = 0° is higher overall than in the region = 180°–270°. It therefore follows that, as shown in Figure 28(b), for = 0.047, attains its maximum value in the vicinity of = 0°, which is larger; the larger value is for the impeller.

From the above considerations, because the time-averaged at the end of the volute winding for low flow rates is large, in the region = 90°–180°, and because the time-averaged at the beginning of the volute winding for large flow rates is large, at = 0°, the high of the pressure side propagates in the circumferential direction and the around the impeller increases. We can assume that accompanying the increase in this , the leakage flow rate also increases, and owing to the high pressure acting on the outer surface of the rear and front shrouds, the axial thrust for the rear and front shroud parts therefore increases. Because the larger the blade outlet angle for the impeller, the higher the of the pressure side and time-averaged at either the beginning or end of the volute winding. We can argue that this tendency is very pronounced.

5. Conclusions

From the results of investigating the effect of the blade outlet angle on the fluctuating hydrodynamic forces in a closed-type centrifugal pump with a single blade, using both experiment and numerical analysis, we have found the following.(1)At low flow rates, all the impellers give a root-mean-square value for the fluctuating component of the total radial thrust that is roughly the same, but for high flow rates, it is larger for impellers with larger blade outlet angles. On the other hand, although the time-averaged value of the total radial thrust also increases with the flow rate for all the impellers, for either low or high flow rates, the difference due to variations in the blade outlet angle is comparatively small.(2)At low flow rates, the radial thrust for the hydraulic part increases for a blade phase angle in the vicinity of = 180° and decreases in the vicinity of = 0°. These tendencies are more significant in impellers with larger blade outlet angles. On the contrary, for high flow rates, the radial thrust for the hydraulic part increases in the vicinity of = 0° and is largest in impellers with large blade outlet angles. Conversely, the radial thrust for the hydraulic part decreases in the vicinity of = 180° and is smallest in impellers with small blade outlet angles.(3)For the time-averaged value of the axial thrust for the rear and front shroud parts, the absolute value gets larger for impellers with large blade outlet angles at high flow rates. Therefore, for the time-averaged value of the total axial thrust, the absolute value is largest for the impeller with a blade outlet angle of 24°, but hardly any difference is seen between the impellers of 8° and 16°.(4)The absolute value of the axial thrust for the rear and front shroud parts is largest, for low flow rates, in the vicinity of = 90°–180° and for high flow rates, in the vicinity of = 0°. This is because the leakage flow rate increases with an increase in the static pressure around the impeller, so that a high pressure acts on the outer surface of the rear and front shrouds and there is a greater tendency for this to occur in impellers with large blade outlet angles.

Nomenclature

:Blade width (m)
:Pressure coefficient
:Leakage flow rate coefficient =
:Impeller diameter (m)
:Axial thrust coefficient
:Radial thrust coefficient
:Axial thrust (N)
:Radial thrust (N)
:Gravitational acceleration (m/s2)
:Total head (m)
:Shaft power (W)
:Rotational speed (min−1)
:Static pressure (gauge pressure) (Pa)
:Leakage flow rate (m3/s)
:Flow rate (m3/s)
:Impeller radius (m)
:Circumferential velocity (m/s)
:Blade angle (°).
Greek Letters
:Pump efficiency
:Circumferential angle (°)
:Blade phase angle (°)
:Shaft power coefficient
:Fluid density (kg/m3)
:Flow rate coefficient
:Head coefficient .
Subscripts
0:Impeller suction inlet
1:Blade inlet
2:Blade outlet
:Axial component
:Radial component
:Suction side reference position
: component
: component
: component.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This work was supported by JSPS KAKENHI Grant no. 24760130.