Next Article in Journal
A Fouling Comparison Study of Algal, Bacterial and Humic Organic Matters in Seawater Desalination Pretreatment Using Ceramic UF Membranes
Next Article in Special Issue
Enhanced Antifouling in Flat-Sheet Polyphenylsulfone Membranes Incorporating Graphene Oxide–Tungsten Oxide for Ultrafiltration Applications
Previous Article in Journal
Effect of Methacrylic Acid Monomer on UV-Grafted Polyethersulfone Forward Osmosis Membrane
Previous Article in Special Issue
Membranes Coated with Graphene-Based Materials: A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Vibration Characteristic Analysis of Hollow Fiber Membrane for Air Dehumidification Using Fluid-Structure Interaction

1
Key Laboratory of Microelectronic Packaging and Assembly Technology of Guangxi Department of Education, School of Mechanical and Electrical Engineering, Guilin University of Electronic Technology, Guilin 541004, China
2
Guangdong Provincial Key Laboratory of Distributed Energy Systems, Dongguan University of Technology, Dongguan 523808, China
3
School of Energy and Built Environment, Guilin University of Aerospace Technology, Guilin 541004, China
4
Department of Science and Technology, Nanning College for Vocational Technology, Nanning 530008, China
5
Guangxi Beitou Environmental Protection & Water Group Co., Ltd., Nanning 530025, China
6
Key Laboratory of Enhanced Heat Transfer and Energy Conservation of Education Ministry, School of Chemistry and Chemical Engineering, South China University of Technology, Guangzhou 510641, China
*
Authors to whom correspondence should be addressed.
Membranes 2023, 13(2), 233; https://doi.org/10.3390/membranes13020233
Submission received: 11 December 2022 / Revised: 19 January 2023 / Accepted: 20 January 2023 / Published: 15 February 2023
(This article belongs to the Special Issue Separation Principles and Applications of Membrane Technology)

Abstract

:
Hollow fiber membrane dehumidification is an effective and economical method of air dehumidification. The hollow fiber membrane module is the critical component of the dehumidification system, which is formed by an arrangement of several hollow fiber membranes. The air stream crosses over the fiber bundles when air dehumidification is performed. The fibers vibrate with the airflow. To investigate the characteristics of the fluid-induced vibration of the hollow fiber membrane, the two-way fluid-structure interaction model under the air-induced condition was established and verified by experiments. The effect of length and air velocity on the vibration and modal of a single hollow fiber membrane was studied, as well as the flow characteristics using the numerical simulation method. The results indicated that the hollow fiber membrane was mainly vibrated by fluid impact in the direction of the airflow. When the air velocity was 1.5 m/s~6 m/s and the membrane length was 100~400 mm, the natural frequency of the membrane was negatively correlated with length and positively correlated with air velocity. Natural frequencies were more sensitive to changes in length than changes in air velocity. The maximum equivalent stress and total deformation increased with air velocity and length. The maximum equivalent stress was concentrated at both ends, and the maximum deformation occurred in the middle. The research results provided a basis for the structural design of hollow fiber membranes under flow-induced vibration conditions.

1. Introduction

People in modern society spend most of their time indoors. Therefore, there is an increasing concern for indoor air quality (IAQ) in houses and offices. Especially with the outbreak of the novel coronavirus in the spring of 2020, people are more aware of IAQ. Indoor air humidity is one of the most important factors affecting IAQ. High or low indoor air humidity can harm human health and comfort. The recommended relative humidity range is 40–60% in dwellings and air-conditioned buildings in summer [1]. During summer, fresh air relative humidity is very high, and sometimes up to 90% in a coastal area. Therefore, it is necessary to dehumidify the fresh air and then supply it to the room. Fresh air dehumidification consumes a large amount of energy, which makes up 20% to 40% of the energy consumed by air conditioning systems [2,3]. The independent temperature and humidity control system can significantly decrease the energy consumption of all fresh air conditioning systems [4,5,6]. The independent temperature and humidity control system can separately handle air sensible heat and latent heat, which can significantly reduce the energy consumption of the air conditioning system. It is one of the key techniques to achieve the goal of “carbon peaking and carbon neutrality” in the field of refrigeration and air conditioning.
There are various air dehumidification approaches, including condensation dehumidification, solid adsorption dehumidification, solution dehumidification, and membrane-based liquid dehumidification. Membrane-based liquid dehumidification has been gaining much attention from both academic and industrial fields because of its high efficiency, environmental friendliness, and energy savings [7,8,9]. Membrane-based liquid dehumidification uses a selective semi-permeable membrane to selectively penetrate vapor, but not liquid and other gases, to realize the dehumidification objective [10,11,12]. Flat membrane dehumidification and hollow fiber membrane dehumidification systems are very effective methods to solve the droplet-carrying problem in traditional liquid dehumidification systems. Compared with the flat membrane, the hollow fiber membrane has the advantages of high-packing density, compactness, and high efficiency, and has been widely concerned. [13,14]. Ghidossi et al. [15] developed a mass transfer model for membranes under laminar and turbulent flow conditions and analyzed the mass transfer behavior of membranes under different flow regimes. Zhang [16] studied the hollow fiber membrane humidifier component of seawater desalination systems. A simulation model was built considering the flow, and heat and mass transfer characteristics of both the solution and air sides, and related experiments were designed to validate the model. Huang et al. [17] studied the performance of elliptical hollow fiber membrane tubes in the coupled flow and heat transfer by modeling elliptical hollow fiber membrane tubes and discovered that as the elliptical semi-axis ratio (b/a) decreased, the mean friction coefficient and the Nusselt number also decreased. Fluids flowing in the short semi-axial direction have improved combined heat transfer and fluid flow properties. Ma et al. [18] investigated the effect of membrane arrangement on heat and mass transfer. The simulation results predicted the effect of fiber distance on membrane module hydrodynamics and filtration performance. Most of these studies are from the point of view of the flow field and ignore the vibration and deformation of the membrane.
Hollow fiber membranes are flexible tubes, so during operation, flow-induced vibrations are generated when airflow crosses the membrane tubes. Flow-induced vibrations transform the shape of hollow fiber membrane tubes, thus affecting fluid flow, heat and mass transfer performance, mechanical performance, and so on. Li et al. [19] studied the effectiveness of vibration in anti-fouling control of a submerged hollow fiber membrane bioreactor. Their results showed that low-frequency continuous transverse vibration is very effective for fouling control, and it is more effective than longitudinal vibration. Low et al. [20] experimentally investigated transverse vibration on hollow fiber membrane tubes to effectively reduce the external concentration polarization and increase forward osmosis. Chai et al. [21] analyzed the performance of rotational vibration in a submerged hollow fiber membrane system for the bioseparation of high-concentration yeast suspensions. Their investigation also showed that the system performance on critical flux for the rotational vibration system is insensitive to packing density. During membrane operation, the fluid-structure interaction cannot be ignored. Fluid-structure interaction has been widely used in heat exchangers [22,23,24]. The effect of fluid-structure interaction on hollow fiber membrane performance has been studied. Zamani et al. [25] studied the effect of longitudinal vibration on the antifouling properties of hollow fiber membranes by the fluid-structure interaction (FSI) method. The effect of membrane surface wall shear rate was analyzed, and it was found that dynamic shear-enhanced filtration by vibration can reduce concentration polarization and membrane contamination. Liu et al. [26] found that the FSI was effectively evaluated in terms of the effects of fiber tube material, size, and relaxation on the shear effect and flux on the membrane tube surface, and the FSI simulated fiber tube displacements show good agreement with experimentally measured fiber displacement data. Sangeetha. [27] investigated the performance of hollow fiber membranes by coupling flow characteristics with structural mechanics to study the interaction between blood flow and membrane structure, and the results showed that the convoluted corrugated fiber model would have better mechanical stability. The majority of these studies have concentrated on the effect of solid domain changes on the entire membrane module. In considering the effects of two-way fluid-structure interaction, Jeong et al. [28] analyzed the hollow fiber membrane modules used for dehumidification in pneumatic systems. The fluid flow and dehumidification characteristics were determined by the FSI method for different shapes of deformed baffles. The difference between the results of the experiment and low analysis was about 5%. Li et al. [29] developed a two-way fluid-structure interaction model for fluid-induced fiber membrane vibration, and the fluid-induced vibration could result in heat and mass transfer enhancement factors of up to 68.9% and 96.2% compared to the non-fluid-induced vibration state. Jang et al. [30] simulated the behavior of fluid and shear stresses in a membrane reactor. The results of the fluid-structure interaction analysis showed that reducing slack can increase the high membrane tension of a membrane. The two-way fluid-structure interaction better reflects the impact of fluid interactions on overall component performance during actual operation.
As mentioned above, the FSI method improves the antifouling ability as well as the heat and mass transfer ability of the hollow fiber membrane to some extent. The vibration characteristics are very important for membrane modules, affecting the mechanical properties of the membrane and process preparation of the module. Few studies have investigated the effect of fluid-structure interaction on the vibrational properties of the membrane itself. The two-way fluid-structure interaction method can provide more accurate results to guide the module design, so the vibration characteristics of hollow fiber membranes for air dehumidification using two-way fluid-structure interaction are investigated in this paper.
To explore the vibration characteristics of the hollow fiber membrane under flow-induced vibrations, this paper establishes the two-way fluid-structure interaction vibration model of the hollow fiber membrane and verifies the mathematical model by the air-induced membrane tube vibration experiment. The effects of membrane length and air velocity on vibration deformation, equivalent stress, and natural frequency of membranes are then investigated using the numerical simulation method. According to the research results of hollow fiber membrane vibration, the goal of this paper is to provide theoretical guidance for membrane module structural design.

2. Numerical Model

2.1. Fluid-Structure Interaction Model Equations

2.1.1. Fluid Domain

The air is an incompressible Newtonian fluid, and the continuity equation is:
u x + v y + w z = 0
The momentum equations are:
u t + u u x + v u y + w u z = 1 ρ a p x + υ 2 u x 2 + 2 u y 2 + 2 u z 2 + G x
u t + u v x + v v y + w v z = 1 ρ a p y + υ 2 v x 2 + 2 v y 2 + 2 v z 2 + G y
w t + u w x + v w y + w w z = 1 ρ a p z + υ 2 w x 2 + 2 w y 2 + 2 w z 2 + G z
where u, v, w represent the velocity components in the x, y, and z directions (m2/s), ρa is the density of the air (kg/m3), p is the air pressure (N), t is time (s), ν is the kinematic viscosity (m2/s), G is the gravity (N), subscripts “x”, “y”, and “z” are the different directions, respectively.

2.1.2. Solid Domain

The solid control equation is given by Newton’s second law; its equation is:
ρ s a s = · σ s + f s
where ρs is the density of the hollow fiber membrane (kg/m3), as is the acceleration vector, σs is the stress tensor, fs is the volume force vector.
The dynamic equation of the membrane under the air-induced condition can be written in the equation: [31]
M D + C V + K A = F
where [M], [C], [K] represent the mass matrix, damping matrix, and stiffness matrix, respectively; {D}, {V}, {A} represent the displacement vector, velocity vector, and acceleration vector, respectively; and {F} represents the force vector.

2.1.3. Coupling Boundary Equations of the Air and Membrane Tube

In the fluid-structure interaction calculation, the displacement (d) and stress (τ) on the fluid-structure interaction interface should meet the equilibrium condition on both sides. In other words, the following equation is satisfied: [32]
d f = d s
n f τ f = n s τ s
where d is the displacement (m), n is the unit normal vector of stress, which is perpendicular to the interface of the fluid and solid. Subscripts “f”, “s” are the fluid domain and solid domain, respectively.

2.2. Geometry and Meshing

Figure 1 depicts the geometric model of a single hollow fiber membrane, and Table 1 depicts the physical parameters of the fiber membrane. After specifying the material properties, the model is meshed. The fluid domain is meshed by the 3D fluid element, with a mesh cell count of 114,000. The solid domain is made of unstructured tetrahedral mesh with a solid187 unit type. The mesh number is 87,057. The meshing is shown in Figure 2.
The mesh of the computational model directly determines the efficiency and accuracy of the computation. This paper uses the maximum deformation and maximum equivalent stress of the hollow fiber membrane as indicators to verify mesh independence. As shown in Table 2, the changes in maximum deformation and maximum equivalent stress of hollow fiber membranes with different mesh numbers (55,110, 108,300, 433,000) of 100 mm hollow fiber membranes at an air velocity of 1.5 m/s are compared. The mesh number rises from 108,300 to 433,000, and the maximum deformation and maximum equivalent stress increase is less than 0.2%. To balance the efficiency and accuracy of the calculation results, a mesh number of 108,300 is selected for solving.

2.3. Boundary Conditions

Considering the influence of gravity, the fluid domain adopts the velocity inlet, the air velocities are 1.5 m/s, 3 m/s, 4.5 m/s, and 6 m/s, respectively, and the corresponding Reynolds numbers are 166, 332, 498, and 664, respectively. When the air flows through the membrane tube, the fluid flow is in a transition state, and the K-epsilon turbulence model is selected [33]; the pressure outlet is set to 0 Pa; the non-slip wall is set; The SIMPLE algorithm is chosen by the velocity-pressure coupling; the interface is set to FSI, and the dynamic mesh model of the diffusion smoothing method is selected; The ends of the membrane are set as fixed constraints in the solid domain, and the outer walls of the membrane are set as coupling surfaces.

2.4. Fluid-Structure Interaction Solution Process

Firstly, the fluid domain and solid domain of a single hollow fiber membrane are established and specified, respectively. Both the solid domain and fluid domain are separated into meshes. The governing equations in the fluid domain are solved using the finite volume approach, and the pressure distribution on the fluid-structure interface is obtained. The system coupling module then transfers this to the solid domain. The solid domain is applied to the solid surface through the set fluid-solid interface based on the pressure obtained from the transfer solution, and the transient dynamic analysis is performed using the finite element method to obtain the displacement of the membrane tube. Finally, the module feeds the analysis results back into the fluid domain. The fluid domain’s governing equations are solved using the new boundary conditions, and the system coupling exchange calculation is repeated until the pre-specified time limit is reached. The data transfer interface in the system coupling module exchanges fluid and solid data with a time step of 0.005 s during the fluid-structure coupling calculation, which takes one second. The effective value determines the convergence of this time step result, which is judged to be convergent when it is less than 1e−4.

2.5. Modal Analysis and Solution

In a two-way fluid-structure interaction modal analysis, the computational fluid dynamics (CFD) flow field calculation and finite element method (FEM) transient coupling analysis are performed first. The velocity and pressure of the fluid are obtained from the fluid domain solution, displacement and stress-strain are obtained from the solid domain solution, and coupling analysis is performed after convergence to a certain accuracy for the pre-stress modal analysis. Fixed constraints are set at both ends of the fiber membrane; the large deformation is turned on to increase the calculation accuracy, and the first six natural frequencies and modes of the hollow fiber membrane are solved.

3. Model Verification

The experimental validation bench was built by members of this group. Figure 3 depicts the hollow fiber membrane fixing module. Figure 4 depicts the whole experimental setup of the hollow fiber membrane flow-induced vibrations. The air velocity was controlled by a three-phase variable frequency fan (SIEMENS® 1LE001, Munich, Germany), and the experimental flow channel was a 5 m long circular air duct with a built-in flow equalizing plate. The fan’s air is subjected to two rectifier flows, and the flow is fully developed when it reaches the end. At this time, it is considered that the air velocity has reached a relatively stable state. The hollow fiber membranes were then measured and observed. During the experiment, to improve the numerical accuracy of the experiment, the air intake velocity was tested in the range of 1.5~10.5 m/s. To ensure the accuracy of the air velocity measurement and reduce the experimental error, the data was collected after the fan frequency was adjusted for 60 s to stabilize it, and a hot-wire anemometer (Testo-425, Testo SE, Lake Titi, Germany) was used for collection five times at each air velocity. The uncertainty of the air velocity measurement is ±0.5%. DENSHIGIKENQ Co., Ltd.’s laser vibrometer (melectro® V-100lm type-Da, DENSHIGIKENQ Co., Ltd., Shoshima City, Tokyo) was used to measure flow-induced fiber deformation (Osaka, Japan), the laser vibrometer (melectro® V-100lm type-DA, DENSHIGIKENQ Co., Ltd., Shoshima City, Tokyo) can only measure vibration displacement in one plane. Since the axial direction of the hollow fiber membrane has been fixedly constrained, we only need to measure the in-plane XZ and plane YZ. After measuring the displacement vibration on a plane, we only need to adjust the placement of the fiber membrane to keep it relative to the laser vibrometer to complete the next experiment.
Table 3 shows the comparison of the results between experimental and simulated values at different air velocities. From the table, it can be seen that the maximum errors between the simulated and experimental data of membrane tube amplitude in the x and y directions are 3% and 5%, respectively. The results also show that the vibration mainly occurs in the x direction. An accurate two-way fluid-structure interaction model can be demonstrated by this. It provides a guarantee for the further simulation analysis below.

4. Results and Discussion

4.1. Effect of Length and Air Velocity on the Hollow Fiber Membrane’s Inherent Frequency and Mode Shape

At air velocities of 1.5 m/s, 3 m/s, 4.5 m/s, and 6 m/s, four hollow fiber membranes with varying lengths of 100 mm, 200 mm, 300 mm, and 400 mm, respectively, are chosen for a two-way fluid-structure interaction modal analysis. Table 4, Table 5, Table 6 and Table 7 show the first 6-order natural frequencies with different fiber lengths at different air velocities.
From Table 4, Table 5, Table 6 and Table 7, under different air velocities, the first six natural frequencies characterize the same features, decreasing with the increase in length and the speed of decrease, gradually slowing down with the decrease itself. At a 1.5 m/s air velocity, the variation in the natural frequency of each order at membrane length increasing from 100 to 200 mm is 9–12 times the variation at the membrane length ascending from 300 to 400 mm. While at an air velocity of 6 m/s, this number is 14–15 times comparing the variation in the natural frequency of each order at different increases of membrane length from 100 to 200 mm and from 300 to 400 mm, respectively. At the same time, the variation in higher order frequency with length is larger than those of lower order frequency. The increase in the length will affect the mass and stiffness of the membrane and eventually vary in natural frequency. The analyses indicate that the natural frequency is negatively correlated with the length.
Under different lengths, the first six natural frequencies characterize the same features, increasing with the increase in air velocity. The increase in air velocity will increase the shear force between the walls and lead to the rise in natural frequency. This rise is small compared to the effect of length on the natural frequency. The natural frequencies of each order at a membrane length of 100 mm at an air velocity of 6 m/s only rise by 1–3 Hz in comparison to that at a 1.5 m/s air velocity. The law increases with the increase in the length of the natural frequencies of each order. While the membrane length is 400 mm, the natural frequencies of each order at the air velocity of 6 m/s increase by 4–8 Hz compared with those at 1.5 m/s air velocity. Through analysis, we get a positive correlation between natural frequency and air velocity. Compared with air velocity variation, the natural frequency is more sensitive to length variation.
With the change in the length and air velocity, the mode shape region corresponding to each order natural frequency of the hollow fiber membrane does not change, but the corresponding value is changed. Therefore, the first six modes are given at different lengths and air velocities for an air velocity of 3 m/s and membrane length of 200 mm under the two-way fluid-structure interaction. As shown in Figure 5, the first two modes are closer to the resonant frequency due to their lower frequencies. Therefore, the shape is similar, and the distribution of deformation positions is consistent; due to the low frequency, which is closer to the resonant frequency, and the fact that the air flow is mainly concentrated in the middle, the vibration in the middle part is more intense. With the increase in the order, the frequency increases, distortion can easily occur in a certain range, the area of severe deformation increases, and the maximum deformation creeps from the middle to both ends. Compared to the first order, the maximum deformation of the sixth order is reduced by 5%. This demonstrates that there is a connection between deformation and order.

4.2. Effect of Length and Air Velocity on the Hollow Fiber Membrane Vibration

Four hollow fiber membranes with varying lengths of 100 mm, 200 mm, 300 mm, and 400 mm were chosen for a two-way fluid-structure interaction vibration study at the air velocities of 1.5 m/s, 3 m/s, 4.5 m/s, and 6 m/s. Figure 6 shows the maximum deformation and maximum equivalent stress changes of the two-way fluid-structure interaction membranes with different fiber lengths at different air velocities.
From Figure 6, under different air velocities, with increasing length, the maximum deformation and maximum equivalent stress of the fiber increase. At 1.5 m/s air velocity, the maximum deformation and maximum equivalent stress of 100 mm fiber length are 0.057 mm and 0.044 MPa, respectively. The maximum deformation and maximum equivalent stress of 400 mm fiber length are 5.304 mm and 0.605 MPa, respectively. The maximum deformation and maximum equivalent stress of 100 mm fiber length are 1% and 7% of the maximum deformation and maximum equivalent stress of 400 mm fiber length, respectively. While at 6 m/s air velocity, the maximum deformation and maximum equivalent stress of 100 mm fiber length are 0.406 mm and 0.317 MPa, respectively. At a 400 mm fiber length, the maximum deformation and maximum equivalent stress are 9.2 mm and 1.694 MPa, respectively. The maximum deformation and maximum equivalent stress of 100 mm fiber length are 4% and 18% of those of 400 mm fiber length, respectively. This indicates that the maximum deformation and maximum equivalent stress of the membrane increase with the increase in the length at different air velocities.
Under different lengths, the maximum deformation and maximum equivalent stress of the fiber increase with the increase in air velocity. At the length 100 mm, the maximum deformation and maximum equivalent stress of the fiber at 1.5 m/s air velocity are 14% and 13% of the maximum deformation and maximum equivalent stress of the fiber at a 6 m/s speed, respectively. When the length is 400 mm, the membrane’s maximum deformation and maximum equivalent stress at 1.5 m/s air velocity are, respectively, 57% and 36% of those values at 6 m/s air velocity. This indicates that as air velocity increases, the maximum deformation and maximum equivalent stress of the membrane do as well. The same law is observed for different lengths.
With the change in the length and air velocity, the deformation and stress regions of the membrane do not change, but the corresponding values are changed. Therefore, the total deformation and equivalent stress distribution cloud diagrams obtained by a two-way fluid-structure interaction at 3 m/s air velocity and 200 mm length of hollow fiber membrane are given under different lengths and different air velocities, as shown in Figure 7. The maximum deformation is mainly in the middle position. The maximum stress is concentrated at both ends. The fixed constraint at both ends improves the stiffness of both ends of the membrane and effectively inhibits the vibration of both ends, so the deformation amount is small but the stress is maximum. Long periods fixed at both ends are likely to lead to stress concentration as air velocity increases, which may affect the life of the membrane. Attention should be paid to the design.

4.3. Flow Characteristics

With the change of the length in air velocity, the velocity vector and pressure distribution area of the fluid flow do not change; only the corresponding values are changed. Therefore, the pressure distribution cloud diagram and flow field velocity vector distribution cloud diagram obtained by a two-way fluid-structure interaction at the membrane length of 200 mm at the air velocity of (t = 1 s) 3 m/s are given under different lengths and different air velocities. As shown in Figure 8, the pressure value in the direction of the windward inlet is significantly greater than the pressure value in the direction of the leeward outlet. This is also the main reason for the fluid-induced vibration of the membrane tube. The velocity density near the wall of the windward inlet is increasing. The velocity gradient changes are increasing. The velocity of the fluid domain on the leeward side of the membrane tube is relatively lower than that of other regions. The data is processed for the pressure on the fluid-structure interaction contact surface. z = 0 cross-section and z = 100 cross-section of the pressure variation with the position in the x direction, as shown in Figure 9. It can be seen that the pressure value at the leftmost end of the contact surface is the largest, and the pressure value decreases along the x direction (airflow direction) and then increases in the opposite direction. The cross-sectional z = 100 position is shifted to the right to the z = 0 position. The flow field characterization shows a trend in vibration change consistent with the experiment.

5. Conclusions

This paper takes a single hollow fiber membrane as the research object, establishes the hollow fiber membrane of a two-way fluid-structure interaction model under the condition of airflow, and verifies the accuracy of the hollow fiber membrane in a two-way fluid-structure interaction model through experiments. The vibration and modal characteristics of single-root hollow fiber membranes with different lengths and air velocities are investigated by numerical simulation, and the flow characteristics are analyzed. The analysis leads to the following conclusions.
In the range of air velocity 1.5 m/s~6 m/s and membrane length 100~400 mm, the natural frequency of the hollow fiber membrane is negatively correlated with its length, and positively correlated with air velocity; compared with the change in air velocity, the natural frequency of each order is more sensitive than the change in length. The trend in higher order frequency variation with length is greater than that of lower order frequency variation. With the increase in the order, the dramatic area of the hollow fiber membrane mode increases, and the maximum deformation creeps from the middle to both ends.
The maximum deformation of the hollow fiber membrane occurs in the middle position, and the maximum equivalent stress is concentrated at both ends. The long time fixed at both ends can easily cause stress concentration, which may affect the life of the membrane. The maximum deformation and equivalent stress of the membrane increase with the increase in length and air velocity. The longer the length and the higher the air velocity, the more obvious the effect. The deformation and stress distribution caused by membrane tube vibration can be predicted based on fluid-structure interaction. The results can guide the design of membrane components.
Through the analysis of the pressure field and velocity field flow, the causes of hollow fiber membrane vibration are shown. The hollow fiber membrane vibration is mainly in the direction of airflow, which is consistent with the experimental law of the laser vibrometer.

Author Contributions

Conceptualization, S.Z.; methodology, T.Z.; software, J.C.; validation, N.L. and Y.D.; investigation, J.C.; resources, C.D.; data curation, J.C.; writing—original draft preparation, J.C.; review and editing, C.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Science Foundation of China, Project No. 51906071. Guangxi Province of Natural Science Foundation, Project No. AD22088044. The Guangdong Provincial Key Laboratory of Distributed Energy Systems, Project No. 2020B1212060075.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Qi, Y.; Liu, J.; Lia, D.; Zhang, H.; Cao, X.; Wei, S.; Yoshino, H. Large-scale and long-term monitoring of the thermal environments and adaptive behaviors in Chinese urban residential buildings. Build. Sci. 2020, 168, 106524. [Google Scholar] [CrossRef]
  2. Zhang, L. Heat and mass transfer in a cross-flow membrane-based enthalpy exchanger under naturally formed boundary conditions. Int. J. Heat Mass Transf. 2007, 50, 151–162. [Google Scholar] [CrossRef]
  3. Zhang, L.; Liang, C.; Pei, L. Conjugate heat and mass transfer in membrane-formed channels in all entry regions. Int. J. Heat Mass Transf. 2010, 53, 815–824. [Google Scholar] [CrossRef]
  4. Huang, S.; Zhang, L. Researches and trends in membrane-based liquid desiccant air dehumidification. Renew. Sustain. Energy Rev. 2013, 28, 425–440. [Google Scholar] [CrossRef]
  5. Zhang, L.; Zhang, N. A heat pump driven and hollow fiber membrane-based liquid desiccant air dehumidification system: Modeling and experimental validation. Energy 2014, 65, 441–451. [Google Scholar] [CrossRef]
  6. Lin, J.; Huang, S.; Wang, R.; Chua, K.J. On the in-depth scaling and dimensional analysis of a cross-flow membrane liquid desiccant dehumidifier. Appl. Energy 2019, 250, 786–800. [Google Scholar] [CrossRef]
  7. Zhang, W.K.; Yang, M.; Chen, J.C.; Tao, S.; Huang, X.; Hu, B.; Huang, S.M. Quasi-counter flow parallel-plate membrane con-tactors (QCPMC) for liquid desiccant air dehumidification: Conjugate heat and mass transfer. Int. J. Therm. Sci. 2018, 134, 665–672. [Google Scholar] [CrossRef]
  8. Qi, R.; Dong, C.; Zhang, L.Z. A review of liquid desiccant air dehumidification: From system to material manipulations. Energy 2020, 215, 109897. [Google Scholar] [CrossRef]
  9. Iulianelli, A.; Drioli, E. Membrane engineering: Latest advancements in gas separation and pre-treatment processes, petrochemical industry and refinery, and future perspectives in emerging applications. Fuel Process. Technol. 2020, 206, 106464. [Google Scholar] [CrossRef]
  10. Liang, C. Research on a refrigeration dehumidification system with membrane-based total heat recovery. Heat Transf. Eng. 2014, 35, 1043–1049. [Google Scholar] [CrossRef]
  11. Liu, X.; Qu, M.; Liu, X.; Wang, L. Membrane-based liquid desiccant air dehumidification: A comprehensive review on materials, components, systems and performances. Renew. Sustain. Energy Rev. 2019, 110, 446–466. [Google Scholar] [CrossRef]
  12. Liu, Y.; Zhang, S.; Cui, X.; Su, J.; Yan, W.; Jin, L. Flow behavior and mass transfer of humid air across fiber membrane bundles. Appl. Therm. Eng. 2022, 207, 118200. [Google Scholar] [CrossRef]
  13. Zhang, L. Heat and mass transfer in a randomly packed hollow fiber membrane module: A fractal model approach. Int. J. Heat Mass Transf. 2011, 54, 2921–2931. [Google Scholar] [CrossRef]
  14. Nazif, A.; Karkhanechi, H.; Saljoughi, E.; Mousavi, S.M.; Matsuyama, H. Effective Parameters on Fabrication and Modification of Braid Hollow Fiber Membranes: A Review. Membranes 2021, 11, 884. [Google Scholar] [CrossRef]
  15. Ghidossi, R.; Veyret, D.; Moulin, P. Computational fluid dynamics applied to membranes: State of the art and opportunities. Chem. Eng. Process. 2006, 45, 437–454. [Google Scholar] [CrossRef]
  16. Zhang, L. An analytical solution to heat and mass transfer in hollow fiber membrane contactors for liquid desiccant air dehumidification. J. Heat Transf. 2011, 133, 092001. [Google Scholar] [CrossRef]
  17. Huang, S.; Yang, M.; Yang, M.; Yang, X. Fluid flow and heat transfer across an elliptical hollow fiber membrane tube bank for air humidification. Int. J. Therm. Sci. 2013, 73, 28–37. [Google Scholar] [CrossRef]
  18. Ma, C.; Liu, Y.; Li, F.; Shen, C.; Huang, M.; Wang, Z.; Cao, C.; Zhou, Q.; Sheng, Y.; Sand, W. CFD simulations of fiber-fiber interaction in a hollow fiber membrane bundle: Fiber distance and position matters. Sep. Purif. Technol. 2020, 209, 707–713. [Google Scholar] [CrossRef]
  19. Li, T.; Law, A.W.; Jiang, Y.; Harijanto, A.K.; Fane, A.G. Fouling control of submerged hollow fibre membrane bioreactor with transverse vibration. J. Membr. Sci. 2016, 505, 216–224. [Google Scholar] [CrossRef]
  20. Low, K.S.; Wang, Y.; Ng, D.Y.F.; Goh, K.; Li, Y.; Wang, R. Understanding the effect of transverse vibration on hollow fiber membranes for submerged forward osmosis processes. J. Membr. Sci. 2020, 610, 118211. [Google Scholar] [CrossRef]
  21. Chai, M.; Ye, Y.; Chen, V. Application of rotational vibration in a submerged hollow fibre membrane system for bioseparation of high concentration yeast suspensions. J. Membr. Sci. 2019, 573, 145–156. [Google Scholar] [CrossRef]
  22. Izhar, A.; Qureshi, A.H.; Khushnood, S. Simulation of fluid structure interaction of heat exchanger tube using ANSYS CFX. Life Sci. J. 2013, 10, 2317–2328. [Google Scholar]
  23. 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. 2016, 107, 544–551. [Google Scholar] [CrossRef]
  24. Khushnood, S.; Nizam, L.A. Experimental study on cross-flow induced vibrations in heat exchanger tube bundle. China Ocean Eng. 2017, 31, 91–97. [Google Scholar] [CrossRef]
  25. Zamani, F.; Law, A.W.K.; Fane, A.G. Hydrodynamic analysis of vibrating hollow fibre membranes. J. Membr. Sci. 2013, 429, 304–312. [Google Scholar] [CrossRef]
  26. Liu, X.; Wang, Y.; Waite, T.D.; Leslie, G. Fluid Structure Interaction analysis of lateral fibre movement in submerged membrane reactors. J. Membr. Sci. 2016, 504, 240–250. [Google Scholar] [CrossRef]
  27. Sangeetha, M.S.; Kandaswamy, A. Certain Investigations on the Performance of Hemodialyser Membrane Based on Morphological Characteristics and Modeling and Simulation Techniques. Ph.D. Thesis, Anna University, Tamil Nadu, India, 2018. [Google Scholar]
  28. Jeong, E.; Khan, H.A.; Yun, S.; Lee, K. Fluid-Solid Interaction Analysis for Improvement in the Dehumidification Characteristics of a Hollow Fiber Membrane Module for Use in a Pneumatic Power Unit. J. Korean Phys. Soc. 2019, 75, 791–800. [Google Scholar] [CrossRef]
  29. Li, Z.; Chen, B.; Liang, C.; Li, N.; Zhao, Y.; Dong, C. The Effect of Flow-Induced Vibration on Heat and Mass Transfer Performance of Hollow Fiber Membranes in the Humidification/Dehumidification Process. Membranes 2021, 11, 918. [Google Scholar] [CrossRef]
  30. Jang, Y.; Ryoo, H.; Shin, Y.; Park, H. Simulation of fluid and shear stress behaviors in reciprocating membrane bioreactors for fouling control using fluid–structure interaction analysis. Sep. Purif. Technol. 2022, 294, 121220. [Google Scholar] [CrossRef]
  31. Hou, H.P.; Rui, T.T.; Deng, R.; Xu, Z.F.; Liu, S.H. Modal analysis of central impression cylinder based on fluid–solid coupling method. J. Low Freq. Noise. Vib. Act. Control 2021, 40, 772–783. [Google Scholar] [CrossRef] [Green Version]
  32. Khayyer, A.; Shimizu, Y.; Gotoh, H.; Hattori, S. Multi-resolution ISPH-SPH for accurate and efficient simulation of hydroelastic fluid-structure interactions in ocean engineering. Ocean Eng. 2021, 226, 108652. [Google Scholar] [CrossRef]
  33. Huang, S.M.; Qiu, D.; Wu, Z.; Yang, M.; Zeng, S.; Chen, Y. Fluid flow and heat transfer across a curved hollow fiber membrane tube bank (CHFMTB): Effects of the tube deformations. Int. J. Heat Mass Transf. 2018, 116, 471–483. [Google Scholar] [CrossRef]
Figure 1. Geometric model.
Figure 1. Geometric model.
Membranes 13 00233 g001
Figure 2. Meshing Diagram: (a) Fluid domain local mesh; (b) Solid domain local grid.
Figure 2. Meshing Diagram: (a) Fluid domain local mesh; (b) Solid domain local grid.
Membranes 13 00233 g002
Figure 3. Hollow fiber membrane fixation module.
Figure 3. Hollow fiber membrane fixation module.
Membranes 13 00233 g003
Figure 4. Diagram of hollow fiber membrane flow-induced vibration experimental apparatus.
Figure 4. Diagram of hollow fiber membrane flow-induced vibration experimental apparatus.
Membranes 13 00233 g004
Figure 5. The first 6 modes of the fluid-structure interaction mode of hollow fiber membrane (length 200 mm, air velocity 3 m/s). (a) the first order; (b) the second order; (c); the third order; (d) the fourth order; (e) the fifth order; (f) the sixth order.
Figure 5. The first 6 modes of the fluid-structure interaction mode of hollow fiber membrane (length 200 mm, air velocity 3 m/s). (a) the first order; (b) the second order; (c); the third order; (d) the fourth order; (e) the fifth order; (f) the sixth order.
Membranes 13 00233 g005aMembranes 13 00233 g005b
Figure 6. The maximum deformation and equivalent stress of fluid-structure interaction at different air velocities for different fiber lengths. (a) 100 mm; (b) 200 mm; (c) 300 mm; (d) 400 mm.
Figure 6. The maximum deformation and equivalent stress of fluid-structure interaction at different air velocities for different fiber lengths. (a) 100 mm; (b) 200 mm; (c) 300 mm; (d) 400 mm.
Membranes 13 00233 g006aMembranes 13 00233 g006b
Figure 7. Cloud diagram of total deformation and equivalence force distribution (length 200 mm, air velocity 3 m/s). (a) Equivalent stress (MPa); (b) Total deformation (mm).
Figure 7. Cloud diagram of total deformation and equivalence force distribution (length 200 mm, air velocity 3 m/s). (a) Equivalent stress (MPa); (b) Total deformation (mm).
Membranes 13 00233 g007
Figure 8. Flow field cloud diagram at z = 100 section position (length 200 mm, air velocity 3 m/s): (a) pressure distribution; (b) velocity vector distribution.
Figure 8. Flow field cloud diagram at z = 100 section position (length 200 mm, air velocity 3 m/s): (a) pressure distribution; (b) velocity vector distribution.
Membranes 13 00233 g008
Figure 9. Scatter plot of pressure variation with the position in x-direction at z = 0 section and z = 100 section position.
Figure 9. Scatter plot of pressure variation with the position in x-direction at z = 0 section and z = 100 section position.
Membranes 13 00233 g009
Table 1. Physical properties parameters.
Table 1. Physical properties parameters.
ParameterValueUnit
Air density1.23kg/m3
Aerodynamic viscosity1.89 × 10−5kg/m·s−1
Membrane material density1610kg/m−3
Membrane Young’s modulus280MPa
Membrane Poisson’s ratio0.38-
Fiber outer diameter1.7mm
Fiber inner diameter1.5mm
Fiber length100mm
Table 2. Mesh independence verification.
Table 2. Mesh independence verification.
Mesh NumberMaximum Deformation
(mm)
Maximum Equivalent Stress
(MPa)
55,1100.04280.0526
108,3000.04330.0536
433,0000.04340.0538
Table 3. Experimental measurements and simulation results of hollow fiber membrane amplitudes in x and y directions at different air velocities.
Table 3. Experimental measurements and simulation results of hollow fiber membrane amplitudes in x and y directions at different air velocities.
Air Velocity (m2/s)Amplitude in the x Direction (mm)Amplitude in the y Direction (mm)
Experimental
(mm)
Numerical
(mm)
Deviation
(%)
Experimental
(mm)
Numerical
(mm)
Deviation
(%)
1.50.02660.0271.510.04460.0475.38
60.3120.322.560.02260.0231.76
10.50.7580.760.260.0510.05−1.96
Table 4. The first 6-order natural frequencies with different fiber lengths at air velocity 1.5 m/s.
Table 4. The first 6-order natural frequencies with different fiber lengths at air velocity 1.5 m/s.
Frequency (Hz)Air Velocity (1.5 m/s)
100 (mm)200 (mm)300 (mm)400 (mm)
f183.9422.0911.955.23
f283.9828.5315.827.26
f3230.1560.0129.5414.48
f4230.1660.0129.5414.50
f5448.1117.6354.7428.37
f6448.1117.7754.7428.41
Table 5. The first 6-order natural frequencies with different fiber lengths at air velocity 3 m/s.
Table 5. The first 6-order natural frequencies with different fiber lengths at air velocity 3 m/s.
Frequency (Hz)Air Velocity (3 m/s)
100 (mm)200 (mm)300 (mm)400 (mm)
f184.0124.7812.477.37
f284.1730.8216.939.93
f3230.2463.2330.3516.30
f4230.2463.2330.3516.30
f5448.19119.2755.6930.31
f6448.21119.856.0531.20
Table 6. The first 6-order natural frequencies with different fiber lengths at air velocity 4.5 m/s.
Table 6. The first 6-order natural frequencies with different fiber lengths at air velocity 4.5 m/s.
Frequency (Hz)Air Velocity (4.5 m/s)
100 (mm)200 (mm)300 (mm)400 (mm)
f184.2926.1012.778.25
f285.0133.6017.4411.95
f3230.6265.0630.8320.78
f4230.6265.0730.8320.30
f5448.61121.3756.2635.89
f6448.66121.9256.2936.17
Table 7. The first 6-order natural frequencies with different fiber lengths at air velocity 6 m/s.
Table 7. The first 6-order natural frequencies with different fiber lengths at air velocity 6 m/s.
Frequency (Hz)Air Velocity (6 m/s)
100 (mm)200 (mm)300 (mm)400 (mm)
f184.9726.5713.419.48
f287.0434.8218.9014.39
f3231.5465.0631.8621.17
f4231.5665.9431.8621.17
f5449.63121.9457.4936.38
f6449.8122.7757.7236.61
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liang, C.; Chen, J.; Li, N.; Dong, Y.; Zhong, T.; Zeng, S.; Dong, C. Vibration Characteristic Analysis of Hollow Fiber Membrane for Air Dehumidification Using Fluid-Structure Interaction. Membranes 2023, 13, 233. https://doi.org/10.3390/membranes13020233

AMA Style

Liang C, Chen J, Li N, Dong Y, Zhong T, Zeng S, Dong C. Vibration Characteristic Analysis of Hollow Fiber Membrane for Air Dehumidification Using Fluid-Structure Interaction. Membranes. 2023; 13(2):233. https://doi.org/10.3390/membranes13020233

Chicago/Turabian Style

Liang, Caihang, Jiaxing Chen, Nanfeng Li, Yanfang Dong, Tao Zhong, Si Zeng, and Chuanshuai Dong. 2023. "Vibration Characteristic Analysis of Hollow Fiber Membrane for Air Dehumidification Using Fluid-Structure Interaction" Membranes 13, no. 2: 233. https://doi.org/10.3390/membranes13020233

APA Style

Liang, C., Chen, J., Li, N., Dong, Y., Zhong, T., Zeng, S., & Dong, C. (2023). Vibration Characteristic Analysis of Hollow Fiber Membrane for Air Dehumidification Using Fluid-Structure Interaction. Membranes, 13(2), 233. https://doi.org/10.3390/membranes13020233

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

Article Metrics

Back to TopTop