Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Scientific Reports volume 14, Article number: 3953 (2024 ) Cite this article plate tube heat exchanger
A two-dimensional tube bundles fluid–structure coupling model was developed using the CFD approach, with a rigid body motion equation and the Newmark integral method. The numerical simulations were performed to determine the vibration coupling properties between various tube bundles of stiffness. Take the corner square tube bundles with a pitch ratio of 1.28 as the research object. The influence of adjacent tubes with different stiffness on the vibration of the central target tube was analyzed. The research results show that the vibration characteristic of tube bundles is affected by the flow field dominant frequency and the inherent frequency of tube bundles. The vibration of adjacent tube bundles significantly impacts the amplitude and frequency of the central target tube. The equal stiffness and large stiffness tubes upstream or downstream inhibit the vibration displacement of the target tube to some extent. The low-stiffness tubes upstream or downstream significantly enhanced the amplitude of the target tube. The findings can be used to provide a basis for reasonable design and vibration suppression of shell-and-tube heat exchangers.
Working as the most widely used heat exchange equipment, the Shell-and-tube heat exchanger has preferable characteristics such as good heat exchange performance and large heat transfer area1,2,3. It can work under high temperatures and pressure, with simple structure, low cost, convenient cleaning, and strong applicability. Heat exchangers occupy a large proportion of industrial investment, accounting for about 30% of the total investment in chemical equipment, 40% of the total investment in oil refining production equipment, and 20% of the total investment in nuclear power plant equipment. The shell-and-tube heat exchangers have the highest application rate in the industrial heat transfer field, accounting for about 70% of the current heat exchange equipment4,5.
For a shell-and-tube heat exchanger, the shell-side fluid scours the tube bundle laterally, leading to a flow-induced vibration, failing the heat exchange tube bundles6,7. Researchers have carried out much work in the past decades using theoretical, experimental, and numerical simulation methods. Their work has solved the problem of flow-induced vibration of heat exchanger tubes to a certain extent8,9,10,11. However, the failure problem caused by coupling vibration between bundles still occurs. Of the 128 heat exchanger tubes in a heat exchanger of Daqing Petrochemical Company, 15 failed and plugged (see Fig. 1). According to the detection and analysis of the failed tubes, the causes of tube failure are the touching abrasion and fatigue generated by the vibration of the tube bundle. The fluid-induced vibration in the heat exchanger is the main factor inducing the vibration, which causes severe economic and engineering losses. Therefore, it is necessary to investigate the coupling vibration between bundles to reduce the harm caused by such problems.
Profiles of failed tube bundles are shown in red rectangular (provided by Daqing Petrochemical Company). (a) Touching abrasion between tubes and baffle plate, (b) tube fatigue cracking.
To study the coupling vibration between tube bundles, Tang et al. used numerical simulation to learn the five primary vortex shedding modes of heat exchange tube bundles in two-dimensional unsteady flow under the action of saturated steam fluid, as well as the influence of shedding vortex on lift performance12,13,14, and coupled CFD method with computational structural dynamics (CSD) method. The fluid–structure coupling response can be predicted in the time domain. Tan15 et al. established a fluid–structure coupling vibration model to study the critical flow velocity of square-arranged bundles, and the calculated results agree well with the experimental data. Tang et al.13 proposed a high-precision, low-cost CFD/CSD high-order coupling algorithm. Fluent software was combined with the dynamic mesh method to systematically study the dynamic response of the structure of the tube bundles in the dense heat exchange tube bundles. Tube drag, lift, frequency, and amplitude were analyzed in the time and frequency domain.
In a typical shell-and-tube heat exchanger, the stiffness of tubes is different under different constraints. The tubes can be generally divided into A, B, and C types according to the different number of baffle plates. As shown in Fig. 2, the failure shell and tube heat exchanger have five baffle plates. A heat exchange tube passes through three baffle plates and two baffle gap areas while the stiffness is centered. B type of heat exchange tube through all five baffle plates, the large stiffness. The C type of heat exchange tubes passes through two baffle plates and three baffle notched areas with small stiffness.
Classification of heat exchange tubes.
The influences of coupling vibration of tube bundles on the central target tube are investigated when the stiffnesses of adjacent tubes are different, especially for the multiple flexibly-mounted tube bundles16,17. As shown in Fig. 1, the location of failure heat exchange tubes is the interface area of the two types of tubes, illustrating that the stiffnesses of tube bundles in the transverse flow are quite different. It should be considered when the coupling vibration characteristics between tubes are analyzed.
Fluent software is adopted to calculate the flow between tube bundles numerically. According to the actual cases of tube bundle damage, the vibration behavior of adjacent multiple bundles and the interaction mechanism between tubes were studied considering the non-uniform bundle stiffness combination. The coupling vibration of the central target tube and adjacent tubes with different stiffness was analyzed from the aspects of the flow field, structural response, and the movement path of the tube bundles to reveal the coupling vibration mechanism between non-uniform bundle stiffness.
Reynolds time-averaged method is used to calculate the transverse flow between bundles. The time-averaged continuity equation and momentum equation in the form of the tensor are as follows18:
It changes velocity and other variables into a time average, where the extra term is Reynolds stress, representing turbulence's influence. The turbulence model needs to be introduced to make the time-averaged equation closed.
Due to the small gaps between tubes, the fluid velocity in the main flow area significantly differs from the velocity between the gaps. Meanwhile, the flow around produces a large number of vortex structures. Therefore, it is essential to select a suitable turbulence model to simulate the flowfield between dense tube bundles to describe the interaction between the fluid and structure honestly.
Researchers have used different turbulence models to numerically calculate the flow field induced by tube beam vibration. For example, Bao19 used RNG k–ε to calculate the three-dimensional flow field between tube bundles with a pitch ratio 1.28. By comparing with the experimental data, Zhao20 believed that the Transition SST model is more suitable for numerical simulation of turbulent flow field between close-packed tube arrays. Darwish21 stated that the SST k–ω was more reliable in solving the flow between a rotated square tube array.
The turbulence model selected will be verified by comparing the computational results of different turbulence models with experimental data. Based on the velocity field experimental data of the staggered bundles of Balabani22, as shown in Fig. 3, a two-dimensional calculation model was established. The calculation mesh is shown in Fig. 4. There are 108 nodes on the tube bundle boundary. The number of encryption layers was 6, the initial height was 0.2 mm, and the mesh growth rate was 1.1. The RNG k–ε, SST k–ω, and Transition SST model were adapted for numerical simulation, and the velocity field distribution was obtained.
The geometric model for turbulence model verification.
Zoom on the numerical mesh with scales.
The comparison of x velocities calculated by different turbulence models is shown in Fig. 5. Total of 4 positions, x/d = 0, 2.5, 4.2, and 6.7, are compared, respectively. Maximum relative error of different turbulence models are shown in Table 1. At the initial position, all three models agree well with the measured results. The transition SST has the most significant calculation error, but the maximum relative error is only 6.3%. At x/d = 2.5, the overall trend of the calculated results of the three models is consistent with the experimental result. Still, in numerical terms, the SST k–ω model is the most compatible with the experiment results, with a maximum relative error of 13.8%. At the same time, the RNG k–ε and the transition SST have the maximum relative error of 27.6% and 38.7%, respectively. With flow development, the difference between different models is more significant. At x/d = 4.2, the maximum relative error of transition SST is 36.3%. The SST k–ω is also the closest to the actual measurement value, with a maximum error of 16.2%. At x/d = 6.7, the maximum error of SST k–ω is 9.3%. According to the above results, the error of the SST k–ω model is the smallest, and it is in the best agreement with the measured results. Therefore, the SST k–ω model is adopted in this work.
Comparison between model calculation and experiment.
The SST k–ω model was proposed by Menter23. Compared to most RANS models, the SST k–ω model provides better prediction of flow separation and performs well at lower pressure gradients, with high accuracy and economy. Different from standard k–ω, the SST k–ω model modifies the model constants and takes into account the transfer of turbulent shear stress, which uses k–ω in the inner region of the boundary layer and switches to k–ε in the free shear flow24.
where \(P = \tau_{ij} \frac{{\partial u_{j} }}{{\partial x_{j} }}\) , k is turbulence kinetic energy. ω is turbulence-specific dissipation rate. \(\mu_{t}\) is turbulence viscosity coefficient.
F1 and F2 are mixed functions.
According to Menter23, the coefficients for the model are shown in Table 2
The fluid-induced vibration of the tube bundles involves the interaction between the fluid and solid domain, which is challenging to calculate numerically25. Based on the computational fluid dynamics method, this paper establishes a simplified fluid-tube coupling vibration calculation model using the dynamic mesh method.
Because the calculation model that directly considers the coupling of tube bundle deformation and flow field change requires a lot of calculation and many iterations, the accuracy is difficult to guarantee. As shown in Fig. 6, the vibration model of the tube bundle is simplified into a two-dimensional spring-damp-mass system to improve the calculation efficiency26,27.
Simplified vibration flexible tube model (m—mass of the tube; k—stiffness of spring; c—damping of the spring-damp-mass system).
The stiffness, damping, and mass of the simplified tube bundles were calculated by referring to references15,28,29, respectively, and the values are shown in Table 3.
The fluid force on the tube bundle is obtained by solving the fluid motion equation. Taking the fluid force as known, the rigid body motion equation on the tube bundle can be solved, and the fluid–structure coupling of the tube bundle system is realized. The rigid body's motion equation at the bundle's boundary is solved by the Newmark-β integral method30,31.
The vibration differential equation of the tube is
where m, c, k, and F are mass, damping, stiffness of the spring-damp-mass system, and fluid force applied to the tube bundles. \(x\) , \(\dot{x}\) , \(\ddot{x}\) are the displacement, velocity and acceleration of tube32,33.
According to the Newmark-β method,
where \(x_{t}\) , \(\dot{x}_{t}\) and \(\ddot{x}_{t}\) are the displacement, velocity, and acceleration at t. \(x_{t + \Delta t}\) , \(\dot{x}_{t + \Delta t}\) and \(\ddot{x}_{t + \Delta t}\) are the displacement, velocity and acceleration at t + Δt, \(\gamma\) and \(\beta\) are coefficient of Newmark-β method. Δt is the time step.
Then, the vibration differential equation of the system at time t + Δt is
\(x_{t + \Delta t}\) can be expressed as \(K^{*} x_{t + \Delta t} = F_{t + \Delta t}^{*}\) .
where \(K^{*}\) is effective stiffness. \(F_{t + \Delta t}^{*}\) is effective load.
When the mass, stiffness, and damping of the vibrating tube bundle are obtained, the dynamic response of the next time can be obtained by substituting the previous time's dynamic response into the system's initial state.
Due to the large number of tubes in the heat exchanger, the cost of numerical calculation is significant, making it unrealistic to model and study the vibration of all tubes. It makes sense to reduce the number of tubes through reasonable simplification34,35. The three-dimensional model of the tube bundle is simplified into a two-dimensional model. At the same time, only a dozen tubes in typical positions are selected for targeted research, and the calculation amount and time are reduced under the premise of accuracy.
The geometric structure of the numerical calculation model is shown in Fig. 7. Take the corner square arrangement tube bundle as the research object. The heat exchange tube is simplified into a bundle with seven rows and columns. Fluid flows in from the left side and out from the right, and the top and bottom sides are symmetric. The distance from the inlet to the first tube bundle is ten times the tube bundle diameter, and from the last tube bundle to the outlet is 20 times the tube diameter. The tube bundle has a diameter of 25 mm and a pitch ratio of 1.28.
Physical model of vibration of different stiffness tube bundles.
Figure 7 also gives the scheme of the tube bundle stiffness combination. Tube 0 is the target tube, and tubes 1–8 are flexible tubes, which have been marked in the dotted box. The calculation method is the same as the target tube, and self-excited vibration can occur. The unlabeled tubes outside the dotted box are rigid, assuming no vibration occurs in the flow field and no particular definition is required. Three kinds of stiffness were considered: large, normal, and small. Large stiffness was taken as two times the normal stiffness, and the small stiffness was 1/2 of the normal stiffness. In the calculation, only the influence of stiffness change is considered. Thus, the mass of the tube bundle remains unchanged, and the inherent frequency changes due to the stiffness change. The stiffness and corresponding inherent frequency are shown in Table 4.
The stiffness of tube 0 is kept unchanged, and the tube bundles' stiffness is changed upstream or downstream of the target tube. The actual structure of the heat exchanger is referred to ensure that there are only two stiffnesses in the flexible tube bundle. The upstream and downstream are not repeated to determine the stiffness combination scheme.
Mesh independence study is very important for the numerical study authenticity. Based on the number of mesh cells, different mesh densities are used to find a compromise between computational cost and accuracy. Numerical calculations were performed on meshes of 140,000, 180,000, 250,000, 310,000 and 400,000 (corresponding to mesh spacings of 1 mm, 0.8 mm, 0.6 mm, 0.5 mm and 0.4 mm) with inlet velocity of 0.5 m/s. To simplify the calculations, all tube bundles are set to be rigid. The lift coefficients of the target tube (tube 0) are compared under different meshes to determine the appropriate mesh size. The results are shown in Fig. 8. After the number of meshes reaches 310,000, as shown in Fig. 8, the calculation results are almost unaffected by the sparsity of the grids. Relative error between number of 310,000 and 400,000 is just 2.7%. So this paper's research is conducted using 310,000 grid cells.
The computational mesh is shown in Fig. 9. Due to the calculation requirements of the dynamic mesh, the combining structure and unstructured mesh is used, in which an unstructured mesh was used in the tube bundles gaps. The near tube wall zones were divided into a structure mesh. The wall surface of the bundle was encrypted. The number of encryption layers was 8, the initial height was 0.2 mm, and the mesh growth rate was 1.1. The bundle's front and rear flow areas adopt a quadrilateral mesh structure.
Because of the need to simultaneously calculate the vibration characteristics of the tube bundle and realize coupling with the fluid, it is necessary to perform transient calculation36. The boundary condition and solution Settings are shown in Table 5.
A vibration experimental facility for a single flexible tube bundle is established to validate the fluid–structure interaction computational model based on dynamic mesh and linear solution methods. The experimental facility is shown in Fig. 10. It consists of three parts: a water circulation system, a data acquisition and analysis system, and the tube bundle vibration unit. The water circulation system comprises a pump, a supply adjustment system, flow and pressure measuring devices, a water tank, and pipe. Within the tube bundle vibration testing unit, the internal space measures 130 mm × 130 mm in cross-section. The experimental tube bundle is arranged in a 3 × 3 square configuration, with a flexible tube at the center, while the other tubes are firmly connected to the tube plate. Springs are employed to constrain the middle tube bundle to the fixed tube plate. Eight springs are evenly distributed at four positions—horizontally and vertically—at both ends of the flexible tube. Sensors can be fixed and connected through bolts at the ends of the central flexible tube. The data acquisition and analysis system is the DH5981 data acquisition system, using the complementary DHDAS software for signal analysis and processing. The sensors utilized are IEPE underwater piezoelectric acceleration sensors, with their relevant parameters detailed in Table 6.
In Figs. 11 and 12, the acceleration time course curves and spectral diagrams for both experimental and numerical simulations are presented under the condition of an inlet flow velocity of 0.1 m/s. It can be observed that, overall, the trends of the two are quite similar, with the simulated values slightly lower in magnitude compared to the experimental values. In the frequency spectrum, the main frequencies of the two are also relatively close.
Time course curve of acceleration.
Figure 13 shows the vorticity distribution between multiple flexibly-mounted tubes with different stiffness. SFM represents a single flexible mounted tube bundle with A target tube stiffness. The meanings of other labels are shown in Fig. 7. The vorticity distribution shows the process of vortex formation, shedding, and adhesion between the tube bundles37. The first tube bundles generate the vortices, which shed, attach, and move downstream along the flow direction. Due to the interaction between the tightly packed tube bundles, the vortices between the bundles could not fully develop until the last row of tube bundles. The vortices with different shedding characteristics fell off and formed a complete Karman vortex street38.
Vorticity of tube bundle with different stiffness combinations; (a) SFM; (b) AAA; (c) BAA; (d) AAB; (e) CAA; (f) AAC.
When the flexible tubes are added around the target tube, the vortex shape remains unchanged, the vortex separation point is more significant than 120°, and the vortex width is narrow. In general, the overall vortex strength of the multi-flexible tube bundles is lower than the single flexibly-mounted tube bundles, so the pressure fluctuation caused by the vortex shedding will also decrease. The vortex shedding patterns on the nine flexible mounted tubes are basically the same for AAA, BAA, and AAB. The influence of the upstream tube on vortex generation and shedding of the downstream tube is weakened, and the parallel vortex on the target tube is less than that on a single flexible mounted tube. When the upstream or downstream are large stiffness tubes, the vibration law of the target tube tends to a single flexible mounted tube because the vibration amplitude of a large stiffness tube decreases.
As a comparison, the upstream or dowstream are small stiffness for CAA and AAC. Under the same condition, the amplitude of the small stiffness tube is larger. The displacement of tube 1 in the lift direction is large, resulting in a certain distance between tube 1 and downstream tube 0 in the lift direction. When the red vortex at the lower part of tube 1 falls off and is attached to tube 0 again, it is already below the former freeze point of tube 0, and the pressure at the lower part of tube 0 increases. Meanwhile, the blue vortex at the upper part of tube 0 is in the shedding stage, and the upper pressure further decreases.
Combining the two causes a larger amplitude of the target tube 0. For AAC, when the low stiffness tube is added downstream, the large vibration of 4 or 6 tubes prevents the shedding of the vortex of 0 tubes to some extent, the pressure release is blocked, and the surface force of 0 tubes is increased. When the vortices above and below the central target tube finish shedding, it is a complete vortex-shedding cycle. It is obtained that the vortex shedding frequency between tubes is about 25 Hz, which is the dominant frequency of the flow field.
The amplitude of target tube 0 is analyzed by time–frequency. The result is shown in Fig. 14. From the amplitude time history curve, and it can be found that compared with the single flexible mounted tube, the amplitude of the target tube is smaller under the equal stiffness AAA, the upstream or downstream large stiffness BAA or AAB. The amplitude of the target tube is significantly suppressed. It should be noted that the upstream and downstream tube bundles are moderately stimulated under this working condition, and the vibration amplitude and phase of the upstream and downstream tubes are in a suitable range, which produces a good vibration suppression effect. For CAA and AAC, the amplitude of the target tube is equal to or larger than that of the single flexible tube, which can be considered to be caused by the flow block and vortex detachment deviation caused by the large motion of the low-stiffness tube. In conjunction with Table 7, it can be seen that the amplitude of AAA is the smallest, and its root mean square amplitude is reduced by more than 50% compared with a single flexible tube. The amplitude of BAA and AAB is about 10% higher than AAA's. The amplitude of the target tube is larger for CAA and AAC, with small stiffness.
Amplitude time–frequency analysis of target tube under different stiffness combinations; (a) time domain; (b) frequency domain.
From the spectrum diagram, we can find that the dominant frequency of the vibration of the target tube is about 25, which is the vortex-shedding frequency of the flow field. The addition of upstream and downstream flexible tubes complicates the target tube's amplitude frequency39. In addition to the primary frequency, there are multiple peaks, which should include the inherent frequency of the bundle and the exciting frequency of turbulence.
The motion trajectories of all flexible tubes were analyzed, and the motion trajectories of tube bundles within 0.3 s were selected, as shown in Figs. 15, 16, 17, 18, 19 and 20. For the equal stiffness AAA, each tube presents a typical eight or elliptic vibration trajectory40, and the amplitude of the tube bundle in the lift direction is slightly larger than that in the drag direction. The tube bundles move approximately periodically. Among them, the amplitudes of the upstream tube bundles 1, 2, 3, 7, and 8 are slightly larger than that of the downstream tube bundles 4, 5, and 6. The front of the upstream tube bundles are direct flow or fixed tubes, and the coupling characteristics between other tubes and a single flexible tube have similarities, so the amplitude is more prominent. The front flexible tubes affect the downstream tube bundles, and the vibration is suppressed.
Tube vibration trajectories of SFM.
Tube vibration trajectories of AAA.
Tube vibration trajectories of BAA.
Tube vibration trajectories of AAB.
Tube vibration trajectories of CAA.
Tube vibration trajectories of AAC.
For BAA, as shown in Fig. 17, the upstream tubes 1, 2, and 8 are large stiffness tubes. Under the premise of no resonance, the amplitude of large stiffness tubes is significantly lower than that of small stiffness tubes, as reflected in the trajectory. The amplitude of large stiffness tubes is one order of magnitude smaller than that of normal stiffness tubes. According to the trajectory, the amplitude of the target tube and other normal stiffness tube bundles downstream did not decrease after adding the upstream large stiffness tube, and the amplitude of the lift direction was more significant than the resistance direction.
The downstream tubes 4, 5, and 6 for AAB are large stiffness tubes. As shown in Fig. 18, the amplitude of the large stiffness tube is significantly lower than that of the normal stiffness tube, regardless of whether it is upstream or downstream. Adding the downstream tubes with large stiffness increases the amplitude of the target tube slightly. The influence of stiffness change of the downstream tube bundle on the target tube is more potent than that of the upstream tube bundle.
As shown in Fig. 19, when small stiffness tubes are added upstream, all flexible tube bundles amplify. Among them, the displacement of small stiffness tubes 1, 2, and 8 in the lift direction is about an order of magnitude larger than that of normal stiffness tubes, and the maximum amplitude on one side is greater than 2 mm, that is, more than 0.08 days. The three tubes with this stiffness are already in resonance from the trajectory. Because of the fluid excitation force generated by the larger amplitude of the upstream small stiffness tubes, the vibration of the downstream normal stiffness tubes is also enhanced, and the increase of the 3 and 7 tubes is the largest.
The AAC is a downstream small stiffness tube bundle. As shown in Fig. 20, under this condition, the small stiffness tubes 4, 5, and 6 have large amplitudes, and the maximum amplitude on one side is more significant than 2.3 mm, which is about 0.1 days. These tubes should be in a resonance state. The trajectories of the three small-stiffness tubes are also slightly different. Tube 5 shows noticeable vibration enhancement in the lift direction. However, the resistance direction is the same as that of the normal stiffness tubes, and the obstruction of the Tube 4 and the Tube 6 in the flow direction. Because of the influence of the small stiffness tube downstream, the vibration of the normal stiffness tube bundle upstream is enhanced.
The vibration response of the AAA combination is analyzed at different flow velocities, and the amplitude-time curves and frequency-domain curves are presented in Fig. 21. The inlet velocity range is 0.1–0.4 m/s.
Amplitude time–frequency analysis of target tube under different stiffness combinations; (a) time domain; (b) frequency domain.
As the flow velocity increases, the dominant frequency of the flow field also increases, showing a basic consistency with the single-flexible tube under the same flow conditions. At lower flow velocities, the vibration frequency of the target tube is relatively singular, exhibiting only one dominant frequency throughout the entire time domain, which is the dominant frequency of the flow field.
The amplitude and dominant frequency of the vibration of target tube 0 both increase with the increase in flow velocity. Additionally, non-dominant frequency components in the spectrum also increase, including both the inherent frequencies of the tube bundle and indistinct turbulent disturbance components. This suggests that the increase in flow velocity enhances the coupling vibration effects between the tube bundles and the degree of turbulence-induced vibration.
Figure 22 illustrates the maximum amplitude and root-mean-square amplitude of target tube 0 concerning various stiffness combinations at different incoming flow velocities. Under AAA, AAB or BAA, the target tube's amplitude notably decreases compared to configurations with lower stiffness in the upstream or downstream sections. Particularly, the upstream higher stiffness configuration consistently demonstrates a significant dampening effect on the downstream target tube across all flow velocities. Once the incoming flow velocity surpasses 0.2 m/s, there is a noticeable sharp increase in the amplitude of lower stiffness tube combinations, either upstream or downstream. This trend aligns with the characteristic behavior of fluid-elastic instability within tube bundles. The presence of lower stiffness tube combinations, whether upstream or downstream, to some extent, diminishes the critical velocity required for fluid-elastic instability, rendering the tube bundle more susceptible to instability.
Maximum amplitude and root-mean-square amplitude of target tube 0; (a) maximum amplitude; (b) root-mean-square amplitude.
The numerical calculation of coupling vibration of heat exchange tube bundles was completed using dynamic grid technology with rigid body motion equation. The vortex and vibration characteristics of tube bundles were analyzed, and the following conclusions were obtained:
Vortices shedding on the tube bundles are the primary cause of flow field fluctuation, which is the root cause of tube bundle vibration. The vortex shedding frequency on the bundles is the dominant frequency of the flow field. The tube bundles take the flow field dominant frequency as the primary vibration frequency and include the tube bundle's inherent frequency.
The vibration of multi-flexible tube bundles changes the vortex shedding distribution between tube bundles, and the vibration of the target tube differs from that of the single flexible tube bundle. On the premise that no resonance occurs, the upstream and downstream flexible tubes can suppress the vibration of the target tube to a certain extent, resulting in lower amplitude and more complex vibration frequency. However, the dominant frequency is the vortex shedding frequency.
The stiffness changes of the upstream and downstream tube bundles will affect the vibration of the central target tube. On the other hand, within the calculated range, no matter whether increasing the stiffness or decreasing the stiffness, the vibration of the central target tube cannot be better suppressed, and the amplitude of the central target tube is the smallest when the stiffness of the upstream and downstream and the target tube is the same. Reducing the stiffness of the tube bundles upstream or downstream will, to some extent, decrease the critical flow velocity of the target tube, making it more prone to experiencing fluid-elastic instability. This work can be used for better design and vibration suppression of shell-and-tube heat exchangers.
The data that support the findings of this study are available from the corresponding author upon reasonable request.
Kumar, A. & Bhattacharya, M. Numerical analysis of aseptic processing of a non-newtonian liquid food in a tubular heat exchanger. Chem. Eng. Commun. 103(1), 27–51 (2016).
Zamzari, F. et al. Numerical investigation of entropy generation and heat transfer of pulsating flow in a horizontal channel with an open cavity. J. Hydrodyn. 29(4), 632–646 (2017).
Ji, J.-D., Ge, P.-q & Bi, W.-b. Numerical analysis of shell-side flow-induced vibration of elastic tube bundle in heat exchanger. J. Hydrodyn. 30, 249–257 (2018).
Zhu, Y., Zeng, Q., Wan, L., Yang, Y. & Li, Z. Vibration response difference of caving mechanism under coal rock impact based on mechanical–hydraulic coupling. Sci. Rep. 13(1), 13794 (2023).
Article ADS CAS PubMed PubMed Central Google Scholar
Bhutta, M. M. A. et al. CFD applications in various heat exchangers design: A review. Appl. Therm. Eng. 32, 1–12 (2012).
Su, Y. et al. A study of the enhanced heat transfer of flow-induced vibration of a new type of heat transfer tube bundle—the planar bending elastic tube bundle. Nucl. Eng. Design 309(dec), 294–302 (2016).
Ji, J., Ge, P. & Bi, W. Numerical analysis on shell-side flow-induced vibration and heat transfer characteristics of elastic tube bundle in heat exchanger. Appl. Therm. Eng. Des. Processes Equip. Econ. 107, 544–551 (2016).
Sadek, O., Mohany, A. & Hassan, M. Numerical investigation of the cross flow fluidelastic forces of two-phase flow in tube bundle. J. Fluids Struct. 79, 171–186 (2018).
Ding, Z. et al. Influence of support gap on flow induced vibration of heat exchange tube. Ann. Nucl. Energy 180, 109443 (2023).
Ricciardi, G., Pettigrew, M. J. & Mureithi, N. W. Fluidelastic instability in a normal triangular tube bundle subjected to air-water cross-flow. J. Pressure Vessel Technol.ogy 133(6), 061301 (2011).
In-Cheol, C., Heung, J. C. & Seungtae, L. Flow-induced vibration of nuclear steam generator U-tubes in two-phase flow. Nucl. Eng. Design 241(5), 1508–1515 (2011).
Tang, D. et al. Investigation of shedding patterns and its influences on lift performances of a cylinder bundle in cross flow. J. Mech. Sci. Technol. 33(6), 2651–2663 (2019).
Tang, D. et al. A CFD/CSD coupled method with high order and its applications in flow induced vibrations of tube arrays in cross flow. Ann. Nucl. Energy 130, 347–356 (2019).
Tang, D. et al. On the number of tubes required to study oscillating vortices and frequency spectrums of tube arrays in cross flow. Ann. Nucl. Energy 124, 198–210 (2019).
Tan, W., Hao, Wu. & Zhu, G. Investigation of the vibration behavior of fluidelastic instability in closely packed square tube arrays. Trans. Tianjin Univ. 25, 124–142 (2019).
Ai, S. et al. Numerical simulation of flow-induced vibration of three-dimensional elastic heat exchanger tube bundle based on fluid-structure coupling. Shock Vib. 25, 1–17 (2022).
Ding, Z. et al. Numerical simulation research on the vibration of helical tube arrays under transverse flow. Energies 15(23), 9082 (2022).
Ansys Inc. Ansys Fluent Theory Guide (Ansys Inc, 2011).
Bao, M., Wang, L., Li, W. & Gao, T. The vibration analysis of tube bundles induced by fluid elastic excitation in shell side of heat exchanger. IOP Conf. Ser. Mater. Sci. Eng. 239, 012011 (2017).
Hai, Z., Puzhen, G., Ruifeng, T. & Xiaochang, L. A three-dimensional refined numerical simulation of cross-flow induced vibration mechanism in the tube bundle. Nucl. Eng. Design 405, 112223 (2023).
Darwish, S. et al. Flow-induced vibrations of a rotated square tube array subjected to single-phase cross-flow. J. Pressure Vessel Technol. 4, 144 (2022).
Balabani, S. & Yianneskis, M. An experimental study of the mean flow and turbulence structure of cross-flow over tube bundles. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2, 25 (2016).
Menter, F. R. Two-equation eddy-viscosity turbulence models for engineering applications. AIAA J. 32(8), 1598–1605 (1994).
Stergiannis, N., van Beeck, J. & Runacres, M. C. Full HAWT rotor CFD simulations using different RANS turbulence models compared with actuator disk and experimental measurements. Wind Energy Sci. Discuss. 2017, 1–20 (2017).
Upnere, S. Numerical study of flow-induced vibrations of multiple flexibly-mounted cylinders in triangular array. Latvian J. Phys. Tech. Sci. 55(5), 43–53 (2018).
Su, W., Tao, K. & Liu, F. Numerical analysis of vibration response of elastic tube bundle of heat exchanger based on fluid structure coupling analysis. Nonlinear Eng. 12(1), 235–239 (2023).
Shahzer, M. A. et al. A comprehensive investigation of vortex-induced vibrations and flow-induced rotation of an elliptic cylinder. Phys. Fluids 34, 3 (2022).
GB, T 151–2014. Heat Exchanger (Standards Press of China, 2014).
Meirovitch, L. Fundamentals of Vibrations (Mcgraw-Hill Higher Education, 2001).
Zheng, M. et al. Numerical investigation of bluff body for vortex induced vibration energy harvesting. Ocean Eng. 213, 107624 (2020).
Sattari, A. S. et al. Study of wave propagation in discontinuous and heterogeneous media with the dynamic lattice method. Sci. Rep. 12(1), 6343 (2022).
Article ADS CAS PubMed PubMed Central Google Scholar
Tewari, A. Atmospheric and Space Flight Dynamics 507–529 (Birkhũser, 2007).
Zhang, J. P., & Pan, L. Three-dimensional modeling and aeroelastic coupling analysis for the wind turbine blade. In 2009 World Non-Grid-Connected Wind Power and Energy Conference, 1–4 (IEEE, 2009).
Hafeez, A. et al. Analysis of flow-induced vibrations in a heat exchanger tube bundle subjected to variable tube flow velocity. Adv. Sci. Technol. Res. J. 17(2), 279–288 (2023).
De Pedro, B. et al. CFD modelling of the cross-flow through normal triangular tube arrays with one tube undergoing forced vibrations or fluidelastic instability. J. Fluids Struct. 64, 67–86 (2016).
Tan, W. et al. Fluidelastic instability research of tube bundles by a two-way fluid-structure interaction simulation. Int. J. Pressure Vessels Piping 199, 104705 (2022).
Joy, A., Joshi, V., Narendran, K. & Ghoshal, R. Piezoelectric energy extraction from a cylinder undergoing vortex-induced vibration using internal resonance. Sci. Rep. 13(1), 6924 (2023).
Article ADS CAS PubMed PubMed Central Google Scholar
Lee, Y. J. et al. Vortex-induced vibration wind energy harvesting by piezoelectric MEMS device in formation. Sci. Rep. 9(1), 20404 (2019).
Article ADS MathSciNet CAS PubMed PubMed Central Google Scholar
Qu, Y. et al. Numerical study on vortex-induced vibrations of a flexible cylinder subjected to multi-directional flows. Phys. Fluids 35, 3 (2023).
Ozgoren, M. & Rockwell, D. Interaction of a deep-water wave with a vertical cylinder: Effect of self-excited vibrations on quantitative flow patterns. J. Fluid Mech. 572, 189–217 (2007).
This research was funded by National Natural Science Foundation of China (Grant No. 11402051), Natural Science Foundation of Heilongjiang (Grant No. LH2022A004), Natural Science Foundation of Heilongjiang (Grant No. LH2022E016), 'Macao Young Scholars Program' (Project code: AM2022003), and Priority Postdoctoral Projects in Zhejiang Province (ZJ2023023).
College of Mechanical Science and Engineering, Northeast Petroleum University, Daqing, 163318, Heilongjiang, China
Yifang Yin, Zunce Wang, Jinglong Zhang, Yan Xu & Lidong Li
Heilongjiang Key Laboratory of Petroleum and Petrochemical Multiphase Treatment and Pollution Prevention, Daqing, 163318, Heilongjiang, China
Zunce Wang, Jinglong Zhang & Yan Xu
National Observation and Research Station of Coastal Ecological Environments in Macao, Macao Environmental Research Institute, Faculty of Innovation Engineering, Macau University of Science and Technology, Macao SAR, 999078, China
Mingyue Ma & Mingming Ge
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
Y.Y.: writing–original draft, investigation, formal analysis, validation. Z.W.: writing–review and editing, methodology, supervision, funding acquisition. M.M.: conceptualization, methodology. J.Z. and Y.X.: supervision, project administration. L.L.: investigation, validation. M.G.: data curation, visualization. All authors read and approved the final manuscript. The authors confirm that this work has not been published before, and its publication has been approved by all co-authors.
Correspondence to Jinglong Zhang, Yan Xu or Mingming Ge.
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/.
Yin, Y., Wang, Z., Ma, M. et al. Coupling vibration analysis of heat exchanger tube bundles under different stiffness conditions. Sci Rep 14, 3953 (2024). https://doi.org/10.1038/s41598-024-53855-x
DOI: https://doi.org/10.1038/s41598-024-53855-x
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
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)
stirred bioreactor Sign up for the Nature Briefing: AI and Robotics newsletter — what matters in AI and robotics research, free to your inbox weekly.