Research article Topical Sections

Solar energy storage in evacuated tubes solar collector using nanofluid embedded in a saturated porous media in the fully developed region: Al2O3 nanofluid embedded in graphite as a saturated porous media

  • Received: 17 April 2021 Accepted: 21 June 2021 Published: 02 July 2021
  • The evacuated tube solar collector is considered an efficient, convenient, and economical option used to convert solar energy into heat. In this work, enhancement of evacuated tubes solar collector performance and the potential for energy storage by using Al2O3 water-based nanofluid embedded in Graphite as a saturated porous media was presented and studied theoretically. The Governing equations derived from the principles of conservation of mass, momentum, and energy were written in a dimensionless form, and solved these equations analytically for the fully developed region, and numerically throughout the entrance region, and the heat flux data for Amman, Jordan was used in this study. The analysis of the effect of different parameters as porosity, pore diameter, nanoparticles solid volume fraction, pressure, and radius of conduit on temperature variation through the pipe and Nusselt number was done. The results show that the flow of nanofluid in porous media enhances the evacuated tube's performance compared to just using water in the same medium under the same conditions. In August, the temperature variation in the evacuated tube at the same condition for Al2O3 water-based nanofluid for 1.5% volume fraction reach 42 ℃ and 73 ℃ for 6% volume fraction. In contrast, In January for the 1.5% volume fraction the temperature variation reached 25 ℃ and 43 ℃ for 6% volume fraction. On the other hand, the stored energy approximately reached 30 kJ/kg and 45 kJ/kg for solid volume fraction 1.5% and 6% respectively at the same conditions.

    Citation: Mohannad B. Khair, Hamzeh M. Duwairi. Solar energy storage in evacuated tubes solar collector using nanofluid embedded in a saturated porous media in the fully developed region: Al2O3 nanofluid embedded in graphite as a saturated porous media[J]. AIMS Energy, 2021, 9(4): 854-881. doi: 10.3934/energy.2021040

    Related Papers:

    [1] Danyang Wang, Ning Chen . Research on Hopf bifurcation of vehicle air compressor system with single-pendulum TMD. Electronic Research Archive, 2025, 33(5): 3285-3304. doi: 10.3934/era.2025145
    [2] Yang Jin, Wanxiang Li, Hongbing Zhang . Transition characteristics of the dynamic behavior of a vehicle wheel-rail vibro-impact system. Electronic Research Archive, 2023, 31(11): 7040-7060. doi: 10.3934/era.2023357
    [3] Yu Chen, Qingyang Meng, Zhibo Liu, Zhuanzhe Zhao, Yongming Liu, Zhijian Tu, Haoran Zhu . Research on filtering method of rolling bearing vibration signal based on improved Morlet wavelet. Electronic Research Archive, 2024, 32(1): 241-262. doi: 10.3934/era.2024012
    [4] Juping Yang, Junguo Wang, Yongxiang Zhao . Simulation of nonlinear characteristics of vertical vibration of railway freight wagon varying with train speed. Electronic Research Archive, 2022, 30(12): 4382-4400. doi: 10.3934/era.2022222
    [5] Jian-Tao Liu, Mei-Tong Fu . Theoretical modeling of thin binder interphases in highly filled PBX composites together with the closed form expression of the effective isotropic moduli of a simplified PBX model. Electronic Research Archive, 2025, 33(2): 1045-1069. doi: 10.3934/era.2025047
    [6] Wen Li, Jian Wang, Zhanpeng Du, Hongfeng Ma, Lijie Zhang, Libin Duan . Lightweight design method and application of MEWP bracket based on multi-level optimization. Electronic Research Archive, 2022, 30(12): 4416-4435. doi: 10.3934/era.2022224
    [7] Yanlong Zhang, Rui Zhang, Li Wang . Oblique impact dynamic analysis of wedge friction damper with Dankowicz dynamic friction. Electronic Research Archive, 2024, 32(2): 962-978. doi: 10.3934/era.2024047
    [8] Xinfeng Ge, Keqin Su . Stability of thermoelastic Timoshenko system with variable delay in the internal feedback. Electronic Research Archive, 2024, 32(5): 3457-3476. doi: 10.3934/era.2024160
    [9] Shuang Cai, Jiwang Zhang, Cong Li, Zequn He, Zhimin Wang . A RUL prediction method of rolling bearings based on degradation detection and deep BiLSTM. Electronic Research Archive, 2024, 32(5): 3145-3161. doi: 10.3934/era.2024144
    [10] Qian He, Wenxin Du, Feng Shi, Jiaping Yu . A fast method for solving time-dependent nonlinear convection diffusion problems. Electronic Research Archive, 2022, 30(6): 2165-2182. doi: 10.3934/era.2022109
  • The evacuated tube solar collector is considered an efficient, convenient, and economical option used to convert solar energy into heat. In this work, enhancement of evacuated tubes solar collector performance and the potential for energy storage by using Al2O3 water-based nanofluid embedded in Graphite as a saturated porous media was presented and studied theoretically. The Governing equations derived from the principles of conservation of mass, momentum, and energy were written in a dimensionless form, and solved these equations analytically for the fully developed region, and numerically throughout the entrance region, and the heat flux data for Amman, Jordan was used in this study. The analysis of the effect of different parameters as porosity, pore diameter, nanoparticles solid volume fraction, pressure, and radius of conduit on temperature variation through the pipe and Nusselt number was done. The results show that the flow of nanofluid in porous media enhances the evacuated tube's performance compared to just using water in the same medium under the same conditions. In August, the temperature variation in the evacuated tube at the same condition for Al2O3 water-based nanofluid for 1.5% volume fraction reach 42 ℃ and 73 ℃ for 6% volume fraction. In contrast, In January for the 1.5% volume fraction the temperature variation reached 25 ℃ and 43 ℃ for 6% volume fraction. On the other hand, the stored energy approximately reached 30 kJ/kg and 45 kJ/kg for solid volume fraction 1.5% and 6% respectively at the same conditions.



    With the development of the engineering equipment, the detrimental effects of low-frequency vibrations have become more pronounced, leading to increasing interest in low-frequency vibration isolation technologies with a broader frequency ranges [1]. In passive vibration isolation technology, structures or materials with nonlinear stiffness are employed to lower the natural frequency while maintaining a high static load-bearing capacity, thereby achieving effective low-frequency isolation [2]. Among these methods, quasi-zero stiffness (QZS) isolation is a representative nonlinear low-frequency isolation technique, characterized by relatively low dynamic stiffness, which enables effective vibration isolation of low frequencies [3]. In practical applications, Li et al. [4] applied QZS isolators to a propulsion shaft system of underwater vehicles, providing a feasible method for effectively isolating longitudinal vibrations transferred from the propulsion shaft to the housing. Wu et al. [5] introduced a quasi-zero stiffness structure into a triboelectric vibration sensor (TVS), achieving a new level of ultra-wide frequency response and high sensitivity and offering insights for the development of high-resolution TVSs. Abuabiah et al. [6] used a QZS structure consisting of two oblique springs and a vertical spring in vehicle seats, resulting in a wider isolation frequency range for the seats.

    The traditional QZS isolators are constructed by a parallel configuration of positive and negative stiffness structures [7,8,9]. Negative stiffness refers to the structural characteristic where the required force decreases as the displacement induced by a load increases. With advancements of quasi-zero stiffness theory, geometric nonlinear structures and novel nonlinear materials have been incorporated into quasi-zero stiffness isolators [10]. Geometric nonlinear structures leverage the unstable states and nonlinear dynamic responses of the structure to achieve negative stiffness [11,12,13,14,15]. Liu et al. [16] studied the nonlinear dynamic response of an origami systems and achieved quasi-zero stiffness vibration isolation. It was demonstrated that superior isolation performance of the designed isolator could be achieved, and the role of geometric nonlinearity effects has been illustrated by simulation studies. Yang et al. [17] established an analytical model for a 3D re-entrant honeycomb auxetic cellular structure, and the negative Poisson's ratio characteristics of the 3D honeycomb model was investigated. Fully recoverable negative stiffness honeycombs were fabricated using thermoplastic and metallic parent materials. It offers nearly ideal shock isolation performance and provides a new direction for the development of quasi-zero stiffness vibration isolators. Figure 1 illustrates three types of QZS structure. The emergence of magnetic negative stiffness technology also provides a new direction for achieving QZS vibration isolation [18]. By utilizing magnetic negative stiffness materials or devices, it is possible to achieve very low stiffness, effectively isolating low-frequency vibrations, especially in cases where traditional elastic isolation systems cannot meet the requirements [19]. Yang et al. [20] improved the vibration isolation performance of a semi-active suspension by integrating the negative stiffness mechanism with a magnetorheological (MR) damper. Zhu et al. [21] designed a novel six-degree-of-freedom vibration isolator using magnetic levitation as the payload support mechanism, which provides QZS in multiple directions and effectively enhances the vibration isolation performance.

    Figure 1.  (a) The traditional QZS isolator [8], (b) origami QZS isolator [16], and (c) honeycomb structure used in vibration isolator [17].

    Smart materials such as the shape memory alloy (SMA) exhibit inherent nonlinear stiffness, making them suitable for direct application in the construction of isolators. Moghadam et al. [22] verified the feasibility of using SMA for vibration isolation in their research. Rustighi et al. [23] implemented a tuned vibration absorber (TVA) using SMA. Du et al. [24] applied SMA in a vibration isolation model and conducted theoretical analysis. Song et al. [25] applied SMA to a civil structure, while Santos and Nunes [26] used real-time temperature modulation of SMA to control the elastic modulus of the elements, enabling continuous adjustment of stiffness over a wide frequency range.

    SMA exhibits unique thermo-mechanical properties, enabling the conversion of mechanical energy into thermal energy and vice versa due to phase transformation. This distinctive, reversible change in physical properties could be induced under appropriate mechanical and thermal loads [27,28].

    As illustrated in Figure 2, the stress-strain relationship of SMA reveals the phenomenon of a "stress plateau" [29]. Within this region, the dynamic stiffness of the SMA is exceptionally low, making it suitable for application in vibration isolation systems as a replacement for traditional linear elastic components. This characteristic, known as "super-elasticity" [30,31], is primarily caused by the phase transformation from austenite to martensite, which will induce significant strain to occur with minimal stress variation. However, this phenomenon is confined to a strain range. Outside of this range, SMA behaves similarly to conventional metals, undergoing elastic deformation in either the austenitic or martensitic phases [26,32,33]. For SMA-based vibration isolation systems, the high static stiffness ensures load-bearing capacity, while the low dynamic stiffness broadens the effective vibration-isolating range.

    Figure 2.  Stress plateau in SMA.

    Here, we propose a quasi-zero stiffness vibration isolation system that primarily utilizes an external SMA springs accumulator as the main nonlinear component. The system will be referred to as SMAS-QZS. The system leverages the "stress plateau" caused by the phase transformation within SMA springs. The constitutive equations for the SMA spring are introduced, and a dynamic model for the SMAS-QZS is developed. Harmonic balance method is employed to derive analytical solutions for the amplitude-frequency response of the SMAS-QZS, and the analytical results are compared with those of a linear vibration isolation system to demonstrate the effectiveness of the quasi-zero stiffness system in isolating low-frequency vibrations. Finally, a prototype of SMAS-QZS for experimental purposes is setup and tested, and the amplitude-frequency response of the SMAS-QZS system under low-frequency sinusoidal displacement excitation is measured and compared with its theoretic counterparts. Both the accuracy of the simulation results and the effectiveness of the SMAS-QZS are verified.

    In this paper, the Landau theory of phase transitions is utilized to describe the phase transformation of SMA spring [34,35]. According to the Landau theory, a fundamental characteristic of phase transitions is the change in symmetry. Typically, high-temperature phases exhibit greater symmetry of lattice structure, while specific symmetrical lattice structures are absent in low-temperature phases. Conversely, low-temperature phases display reduced symmetry, a phenomenon referred to as "symmetry loss", which arises due to the absence of certain symmetry structure present in high-temperature phase.

    Changes in lattice structure symmetry reflect variations in the degree of order of the system, where a reduction in symmetry corresponds to an increase in order. The degree of order is quantified by a parameter known as order parameter, which is a fundamental parameter in characterizing the phase transition process. In high-symmetry phases, the order parameter is equal to zero, while in low-symmetry phases, it takes on a non-zero value.

    Landau proposed that the free energy function of a system can be expressed as a power series expansion in terms of the order parameter and temperature, retaining only the first few terms for practical applicability. The coefficients of the power series are determined based on the specific properties of the material. In this way, the phase transition process can be phenomenologically described through the free energy function. The essential part of modeling phase transition is to construct a free energy function where its minima represent the stable states of the material. Consequently, the phase transition process is characterized by the changes in the shape of the free energy function at different temperatures.

    In this study, the Landau theory is applied to the SMA springs model, in which the Landau free energy function is constructed to model the transformation between the martensitic and austenitic phases. In the models proposed by Wang and Melnik [36], the strain ε is used as the order parameter to describe the phase transformation. The Landau free energy function is expressed as a function of the selected order parameter ε and temperature θ:

    Ψ(ε,θ)=k1(θθ0)2ε2+k24ε4+k36ε6, (1)

    where, Ψ(ε,θ) is the Landau free energy, k1, k2, and k3 are basic material parameters, and θ0 is the reference transformation temperature.

    For SMA springs, the helical structure undergoes deformation under tensile and compressive loading. The axial force generates torque and shear stress on the cross-section of the spring. As a result, it is essential to consider the phase transformation of the SMA under torsional stress. In this paper, the conventional strain is replaced by torsional strain as the order parameter in the Landau free energy function. Consequently, Eq (1) is modified as shown in the following equation:

    Ψ(γ,θ)=k1(θθ0)2γ2+k24γ4+k26γ6, (2)

    where γ is torsional strain.

    By plotting the free energy function with different temperatures, it is illustrated in Figure 3 that the free energy is strongly influenced by the material temperature [37]. It is evidenced that changes in SMA temperature induce phase transformations within the material. This is because, within a specific temperature range, certain phases lose stability and is transformed into other phases.

    Figure 3.  Landau free energy under different temperatures.

    According to the constitutive relationship:

    T=Ψγ. (3)

    Since the interest of the current investigation is about the hysteretic characteristics of the SMA spring as a whole, the distribution of torsional strain across the cross-section of the SMA spring is disregarded. In this context, the torsional strain is simplified to the torsional angle. The constitutive function can be expressed as follows:

    T=k1(θθ0)α+k2α3+k3α5, (4)

    where T represents the torque, and α is the torsional angle. For simplicity, k1, k2, andk3 are redefined parameters that incorporate basic material and geometric properties. To facilitate subsequent calculations of vibration isolation performance, it is necessary to convert the torsional angle into displacement, thereby establishing the force-displacement relationship for the SMA spring. When an external forceF is applied, the spring undergoes a deformation x. The relationship between x and the torsional angleα, as well as the external force F and the torque T of the SMA wire, is expressed as:

    T=FR,x=αR, (5)

    where R is spring mean radius; thus, the constitutive function of SMA spring could be expressed as:

    F=k1(θθ0)R2x+k2R4x3+k3R6x5. (6)

    In Section 2.1, the static constitutive equation for the SMA spring was derived. However, for the dynamic analysis of the vibration isolation system, it is essential to establish the dynamic model of the spring to accurately capture the hysteretic characteristics of the SMA spring. To achieve this, the Euler-Lagrange equation is employed to formulate the equation [38]. The general Euler-Lagrange equation is expressed as follows:

    ddt(L˙q)=Lq, (7)

    where L is the Lagrangian, q denotes the generalized coordinates, and ˙q is generalized velocity. For a general physical system, the Lagrangian is typically defined as the difference between the kinetic energy and the potential energy, expressed as:

    L=VΨ, (8)

    where V is the kinetic energy, and Ψ is the potential energy.

    When the SMA spring is compressed, the system accumulates potential energy, which can be expressed as:

    Ψ=0xF(x)dx. (9)

    Considering the phase transformation characteristics of the SMA spring, it is crucial to account for dissipative potential energy of the system. Therefore, a Rayleigh dissipative term is introduced, which is expressed as:

    R=12υ(dxdt)2, (10)

    where υ is the dissipation parameter of the material.

    Therefore, L can be expressed as:

    L=l0(ρ2(˙x)2Ψ)dx, (11)

    where ρ represents the mass density of the material, x is the deformation of SMA spring, and ρ2(˙x)2 is the kinetic energy density of the system.

    Substituting the potential energy Ψ and the dissipative term R into the Lagrange equation, the dynamic equation of the SMA spring is expressed as:

    d2xdt2+νdxdt=A1(θθ0)x+A2x3+A3x5+f(t). (12)

    To facilitate subsequent expressions, A1, A2, andA3 are newly defined parameters that encompass material and geometric properties, and f(t) is the external force. In this study, the inertial term has a negligible effect on the system compared to the dissipative term [39], so that Eq (12) can be simplified as:

    ν˙x(A1(θθ0)x+A2x3+A3x5)=f(t). (13)

    For further calculation, the parameters A1, A2, andA3 need to be obtained by experimental measurement, and two types of SMA wires with different diameters were selected for the purpose. Due to experimental limitations, only SMA tension springs were fabricated. The SMA wire was wound around a screw rod and secured with clamps after winding. The assembly was then placed in a furnace and heated to 773.15 K for one hour. When it was cooled in air, the formed SMA spring was extracted from the screw rod. The SMA spring prepared for this study is shown in Figure 4. The parameters of SMA spring are shown in Table 1.

    Figure 4.  Two types of SMA springs.
    Table 1.  Parameters of two types of SMA springs.
    The Parameter of SMA spring (NiTi alloy) 1.0 mm wire SMA spring 1.5 mm wire SMA spring
    Wire diameter (mm) 1 1.5
    Mean diameter (mm) 7.5 8.2
    Effective number of turns 30 30
    Austenite phase transformation completely temperature (K) 258 258

     | Show Table
    DownLoad: CSV

    For newly formed SMA springs, it is crucial to perform multiple cycles of tensile training until the stimulation-response behavior of the spring reaches a stable state before proceeding with further analysis [29]. During each tensile cycle, it must be ensured that the SMA spring is transformed from the austenitic phase to the martensitic phase during the loading process and that it reverts back to the austenitic phase during the unloading process. As illustrated in Figure 5, one end of the spring is fixed to the frame, while the other end is attached to a force gauge. The force gauge is driven by a stepper motor to enforce cyclic tensile loading. A total of 25 loading cycles of training were performed.

    Figure 5.  Cyclic tension test.

    It is shown in Figure 6 that both types of SMA springs exhibit significant hysteresis behavior, and the temperature is 285 K. As the number of cycles increases, the residual strain of the SMA springs continues to grow. However, after a certain number of cyclic trainings, the residual strain stabilizes, and the hysteresis loop for each cycle becomes nearly identical. This indicates that the training process of the SMA spring is completed, and the measured stimulation-response curve is appropriate for parameter identification and further analysis. For clarification purposes, the experimental tension-deformation curve with hysteresis loop of the final cycle is presented in Figure 6 by thick solid lines.

    Figure 6.  Cyclic loading of SMA springs. (a) 1.0 mm SMA wire and (b) 1.5 mm SMA wire.

    In order to estimate the model parameters, the final cycle of the tensile training process for both sets of springs are taken for parameter identification purposes. Concerning the numerical stability of the hysteretic dynamics, the mentioned constitutive model is approximated by a lower piecewise polynomial instead of a high order polynomial. Moreover, the nonlinear least square approximation method is employed to estimate the parameters in the model, aiming to minimize the discrepancy between the simulated values and the experimental data. Consequently, the parameter identification problem is reformulated as a nonlinear optimization problem, expressed by the following equation:

    minp,q,n,vG=Mi=1( CiCi)2, (14)

    where p,q,andn represent the values on the y-axis at the interpolation points, v is the time constant, M is the number of selected experimental data points,  Ci denotes the displacement value corresponding to the i-th data point from the experimental results, and Ci is the simulated displacement. To achieve precise identification of the four parameters simultaneously, the Nelder-Mead nonlinear optimization algorithm is employed and implemented using the MATLAB's built-in fminsearch function [40].

    For the 1.0 mm wire SMA spring, the coordinates of the interpolation points are: (0, 0), (12.03, p), (13.91, q), and (19.60, n). For the 1.5 mm wire SMA spring, the coordinates of interpolation points are: (0, 0), (7.92, p), (11.74, q), and (18.90, n). p,q,andn are the y-coordinates of these interpolation points. Those points are marked in Figure 7, where p indicates the beginning of the transformation of austenite to martensite within the SMA spring, n represents the end of the transformation, and q indicates the reverse process.

    Figure 7.  The result of parameter identification: (a) 1.0 mm wire SMA spring and (b) 1.5 mm wire SMA spring.

    In Figure 7, the experimental values and simulated values are compared. It can be observed that the constitutive model, given as the differential equation, effectively captures the super-elasticity effect of the SMA spring and accurately predicts the forces at various states during the tensile cycles. The results of the parameter identification based on the experimental data are presented in Table 2 and 3. From Figure 7, it is evident that the hysteresis phenomenon in the 1.5 mm wire SMA spring is particularly pronounced. The upper plateau demonstrates a strong load-bearing capacity, while the lower plateau could support a smaller load, with a less distinct stress plateau effect. In contrast, the 1.0 mm wire SMA spring exhibits a more distinct stress plateau phenomenon, and the distance between the upper and lower plateau of the hysteresis loop is smaller; this characteristic makes it more suitable for use in the SMAS-QZS vibration isolator system. Therefore, the 1.0 mm wire SMA spring is selected for vibration isolation.

    Table 2.  Parameters identified of 1.0 mm wire SMA spring.
    v p q n
    Initial value 0.01 23.60 19.72 30.65
    Optimized value 0.22 21.50 21.68 31.41

     | Show Table
    DownLoad: CSV
    Table 3.  Parameters identified for the 1.5 mm wire SMA spring.
    v p q n
    Initial value 0.01 29.1 16.16 59.8
    Optimized value 0.157 22.615 19.354 58.650

     | Show Table
    DownLoad: CSV

    Due to fabrication constraints, only tension springs were fabricated, which differ from linear compression springs commonly used in traditional passive vibration isolation systems. To integrate the tension springs into a quasi-zero stiffness system, the SMA springs are installed within an external accumulator, as shown in Figure 8. The function of the accumulator is to provide displacement compensation for the plunger and ensure adequate load-bearing capacity. It consists of a plunger cylinder and several SMA springs. The accumulator is connected via a tube to a plunger, and a mass block m is mounted on top of the plunger. One end of each spring is fixed to the frame and the other end is attached to the plunger. These springs generate both static supporting force and dynamic force, and the force is transmitted to the mass block m through hydraulic oil. A cam exciter is installed beneath the plunger to provide sinusoidal displacement excitation. To enhance the quasi-zero stiffness characteristics of the system, the springs are pre-stretched into the stress plateau range, as illustrated in the range marked in Figure 7(a).

    Figure 8.  The model of the proposed SMAS-QZS vibration isolator.

    The dynamic equation for the mass block can be expressed as follows:

    m¨x+FC=FQZS, (15)

    where x is the displacement of the mass, FC denotes the system damping force, and FQZS is the force provided by the external SMA springs accumulator through the hydraulic oil. Due to the difference in cross-sectional areas of the plunger and the accumulator, the displacements of the mass block and the SMA springs are not identical, and the equation can be written as follows:

    (x0x)S1=x1S2FQZSS1=FSMAS2FSMA=n(ν˙x1A1(θθ0)x1A2x13A3x15)FC=c˙(xx0), (16)

    where x0 is the sinusoidal displacement excitation, ω is the excitation frequency, x1 denotes the displacement of the plunger in accumulator and also the displacement of the SMA springs, FSMA is the force of SMA springs, S1 and S2 are the cross-sectional areas of the plunger and the accumulator separately, n is the number of SMA springs, and c is the system damping coefficient, which is taken from reference [36]. By substituting Eq (16) into Eq (15), the dynamic equation can be rewritten as:

    m¨x+c˙(xx0)=nS[υS˙(x0x)A1(θθ0)S(x0x)A2S3(x0x)3A3S5(x0x)5], (17)

    where S=S1/S2. To simplify the process of calculation, the relative displacement y=xx0 is substituted into Eq (17):

    ¨y+C˙y+a1(θθ0)y+a2y3+a3y5=¨x0, (18)

    where C=nυS2+cm, a1=nS2A1m, a2=nS4A2m, and a3=nS6A3m.

    The excitation part of the model is an eccentric cam structure. The cam is designed as a cylinder with a radius of r2=28mm, featuring an eccentric hole with a distance of d=10mm. The top surface of the cam is in contact with a roller, which has a radius of r1=6.5mm. The motion trajectory is illustrated in Figure 9.

    Figure 9.  Geometric analysis of the cam: (a) Cam and driven part, (b) low position, (c) random position, and (d) high position.

    According to the sine theorem, the relationship between those geometrical parameters of the cam can be expressed as follows:

    r2sin(πη1)=dsinη2=hsin(η1η2),sinη2=dr2sinη1. (19)

    The displacement of the excitation can be expressed as follows (set the initial height as 0):

    h=r22(dsinη1)2dcosη1(r2d). (20)

    As illustrated in Figure 10, the excitation curve fits the theoretical sinusoidal curve.

    Figure 10.  Excitation curve.

    Therefore, the sinusoidal displacement excitation can be represented as: x0=xmcos(ωt), and xm is the amplitude of excitation.

    The HBM [41] is used to obtain analytical solution of Eq (18). The solution of the function is expanded into Fourier series as follows:

    y=y0+Nn=1[αncos(nω0t)+βnsin(nω0t)]. (21)

    The system starts from the static state, so y0=0. The HBM considers only the first-order harmonic approximation and neglects the higher-order harmonic. It is assumed that the relative displacement y shares the same frequency as the external excitation:

    y=αcos(ωt)+βsin(ωt). (22)

    Substitute Eq (22) into Eq (18) and make the coefficients of the sin(ωt) and cos(ωt) terms on either side of the equation equal:

    ω2α+Cωβ+a1(θθ0)α+34a2(α2+β2)α+58a3(α2+β2)2α=ω2xmω2βCωα+a1(θθ0)β+34a2(α2+β2)β+58a3(α2+β2)2β=0. (23)

    The value of response amplitude Y=α2+β2 is calculated by Eq (23):

    Y=ω4xm2(ω2+a1(θθ0)+34a2Y2+58a3Y4)2+C2ω2. (24)

    The response amplitude of the model is represented in a non-dimensional displacement form:

    X=xxm. (25)

    To better evaluate the performance of the SMAS-QZS in isolating low-frequency vibrations, a damped linear vibration isolation system is introduced for comparison. The linear system utilizes the same damping value and an appropriate stiffness. The calculation process of the linear system is not elaborated here, and the result is presented as follows:

    Xlinear=xlinearxm=1(1r2ω)2+(2ξrω)2, (26)

    where rω is the frequency ratio, and ξ denotes the damping ratio.

    As illustrated in Figure 11, compared to the linear system, the response amplitude curve of SMAS-QZS exhibits strong nonlinear characteristics. Due to the inherent stress plateau characteristics of the SMA spring, the natural frequency of the system is significantly reduced. A lower natural frequency of system could effectively broaden the vibration isolation frequency range, especially the low-frequency. Therefore, the SMAS-QZS could effectively isolate the low-frequency vibration. Additionally, the phase transition between austenite and martensite within SMA during the stretching process dissipates vibration energy and contributes to amplitude suppression. The reference temperature is θ0=208K, and the values of the parameters are: a1=229.8,a2=1.8×103,a3=1.5×106,andC=70.

    Figure 11.  The response amplitude: (a) Comparison with the linear system, and (b) comparison with different damping ratios.

    In Figure 11(a), due to the nonlinear characteristics of the system, the response amplitude exhibits multiple solutions within a certain region; thus, the system is unstable within this region. The impact of different damping ratios on the response amplitude of the system is illustrated in Figure 11(b). As the damping ratio increases, the response amplitude decreases.

    In order to validate the influence of temperature on the isolation performance of SMAS-QZS, the response amplitude at different temperatures is compared. As illustrated in Figure 12, it can be observed that as the temperature increases, the degree of nonlinearity in the system decreases. When the temperature exceeds a certain threshold, the characteristics of response amplitude become nearly identical to those of a linear system. This is because the austenite in the material at high temperatures cannot transform into martensite under the stress, and the stress plateau characteristics disappear as a result [29]. Consequently, the quasi-zero stiffness characteristics of the SMAS-QZS nearly vanishes at high temperature.

    Figure 12.  The response amplitude at different temperatures.

    To validate the accuracy of the simulation results, an experimental system is established. The theoretical structure of the experimental system is illustrated in Figure 13.

    Figure 13.  Theoretical structure of the experimental system: (a) Excitation component, and (b) hydraulic component.

    The excitation component is illustrated in Figure 14 and consists of a stepper motor and a cam mechanism. The motor is controlled by the Arduino microcontroller and provides the rotational motion in different frequencies for the cam mechanism. The cam mechanism converts the rotational motion of the motor into a sinusoidal displacement excitation.

    Figure 14.  Excitation component.

    The SMA springs accumulator proposed in this study is illustrated in Figure 15(a). It leverages the stress plateau characteristics of SMA to achieve quasi-zero stiffness characteristics. Due to the limitation of using only tension SMA springs, an experimental design is developed. The SMA springs accumulator mostly consists of a plunger and several SMA springs. The plunger within the accumulator is connected to the load-bearing plunger cylinder via a pipe to enable force transmission.

    Figure 15.  (a) SMA springs accumulator, and (b) SMAS-QZS experimental model.

    The force comes from a certain number of SMA springs arranged circumferentially around the plunger cylinder. The number of SMA springs can be adjusted to 2, 4, 6, or 8 depending on the load weight. To ensure the piston moves only vertically, an external guiding mechanism is incorporated consisting of sliders and shafts. As shown in Figure 15(b), sliders and shafts are utilized to ensure the load-bearing platform moves only in the vertical direction.

    A retractable displacement sensor is mounted with one end fixed to the frame and the other end attached to the load-bearing platform, measuring the response displacement of the platform during the experimental process. All sensor signals are processed through an amplifier before being transmitted to a data acquisition card, ensuring the accuracy and reliability of the data. The experiment setup and prototype of SMAS-QZS is shown in Figure 16.

    Figure 16.  SMAS-QZS experiment system.

    The response displacement and simulation results of SMAS-QZS at 0.5, 1.0, 1.5, and 2.0 Hz are illustrated in Figure 17; the data is processed using filters to remove noise signals. Due to the output error of the motor, there is a certain phase difference between the experimental data and the simulation results. From Figure 17, it can be observed that experimental response displacement aligns well with the simulation results across excitation frequencies.

    Figure 17.  Response displacement at different frequencies.

    As illustrated in Figure 17, the response displacement exceeds the amplitude of the excitation at an input frequency of 0.5 Hz. This is because 0.5 Hz is within the resonance region, and the system is unable to isolate vibrations at this frequency. At frequencies above 1.0 Hz, the vibration isolation capability of the system is significantly enhanced, and the vibration suppression effect for the excitation displacement reaches over 50%. This phenomenon verifies the accuracy of the displacement response simulation results. When the excitation frequency exceeds 0.5 Hz, it surpasses the natural frequency of the system. Therefore, the resonance phenomenon disappears, and the vibration isolation capability of the system becomes apparent. Consequently, the system effectively isolates vibrations with the frequency above 0.5 Hz, and the effectiveness of the system in isolating low-frequency vibrations also demonstrates the quasi-zero stiffness characteristic.

    In this paper, a novel quasi-zero stiffness vibration isolation model using SMA springs is proposed. First, the Landau theory is utilized to construct the static constitutive equation of the SMA spring. The static constitutive equation can describe the characteristics of the phase transformation from austenite to martensite. Subsequently, the Euler-Lagrange equation is utilized to develop the dynamic model.

    To obtain the parameters of the dynamic model, tensile tests are conducted on SMA springs with different wire diameters to obtain the force-displacement characteristics, and the parameters are identified using the nonlinear least square optimization method. To effectively leverage the stress plateau characteristics of the SMA spring, an external SMA springs accumulator is employed to develop a theoretical quasi-zero stiffness vibration isolation model (SMAS-QZS), and the SMA springs are pre-stretched to a specific length. The HBM is used to derive an analytical solution for the amplitude-frequency response. The analytical results are compared with those of a linear vibration isolation system to demonstrate the effectiveness of the quasi-zero stiffness system in isolating low-frequency vibrations. Furthermore, we examine the effects of different parameters on the vibration isolation performance of SMAS-QZS, such as damping ratio and temperature.

    Finally, a prototype of SMAS-QZS for the experimental purpose is set up and tested, and the amplitude frequency response of the experimental system under low-frequency sinusoidal displacement excitation is measured and compared with its theoretic counterparts. By applying external excitations within the frequency range of 0.5–2.0 Hz, the experimental data aligns well with the simulation results. The system is proven to effectively isolate vibrations at frequencies above 0.5 Hz.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    This work has been supported by the National Natural Science Foundation of China (Grant Nos. 51575478 and 61571007).

    The authors declare there is no conflict of interest.



    [1] Sabiha MA, Saidur R, Mekhilef S, et al. (2015) Progress and latest developments of evacuated tube solar collectors. Renewable Sustainable Energy Rev 51: 1038-1054. doi: 10.1016/j.rser.2015.07.016
    [2] Chopra K, Tyagi VV, Pandey AK, et al. (2018) Global advancement on experimental and thermal analysis of evacuated tube collector with and without heat pipe systems and possible applications. Appl Energy 228: 351-389. doi: 10.1016/j.apenergy.2018.06.067
    [3] Choi SUS (1995) Enhancing thermal conductivity of fluids with nanoparticles. American Society of Mechanical Engineers, Fluids Eng Div 231: 99-105.
    [4] Baqaie Saryazdi A, Talebi F, Armaghani T, et al. (2016) Numerical study of forced convection flow and heat transfer of a nanofluid flowing inside a straight circular pipe filled with a saturated porous medium. Eur Phys J Plus 131.
    [5] Kasaeian A, Daneshazarian R, Mahian O, et al. (2017) Nanofluid flow and heat transfer in porous media: a review of the latest developments. Int J Heat Mass Transfer 107: 778-791. doi: 10.1016/j.ijheatmasstransfer.2016.11.074
    [6] Gautam A, Saini RP (2020) A review on technical, applications and economic aspect of packed bed solar thermal energy storage system. J Energy Storage 27: 101046. doi: 10.1016/j.est.2019.101046
    [7] Nield, DA, Bejan A (2017) Convection in Porous Media. Available from: https://doi.org/10.1007/978-3-319-49562-0.
    [8] Mahmoudi Y, Hooman K, Vafai K (2019) Convective Heat Transfer in Porous Media.
    [9] Muhammad MJ, Muhammad IA, Sidik NAC, et al. (2016) Thermal performance enhancement of flat-plate and evacuated tube solar collectors using nanofluid: a review. Int Commun Heat Mass Transfer 76: 6-15. doi: 10.1016/j.icheatmasstransfer.2016.05.009
    [10] Rashidi S, Esfahani JA, Rashidi A (2017) A review on the applications of porous materials in solar energy systems. Renewable Sustainable Energy Rev 73: 1198-1210. doi: 10.1016/j.rser.2017.02.028
    [11] Ahmed OK, Mohammed ZA (2017) Influence of porous media on the performance of hybrid PV/Thermal collector. Renewable Energy 112: 378-387. doi: 10.1016/j.renene.2017.05.061
    [12] Tajik Jamal-Abad M (2017) Experimental investigation on the effect of partially metal foam inside the absorber of parabolic trough solar collector. Int J Eng 30: 281-287.
    [13] AlMasa'deh HA, Duwairi HM. Modeling of Fluid Flow and Heat Transfer inside a Saturated Porous Conduit at Constant Surface Heat Flux. Available from: https://doi.org/10.18280/mmc_b.860309.
    [14] Saedodin SAHZS, Zamzamian SAH, Nimvari ME, et al. (2017) Performance evaluation of a flat-plate solar collector filled with porous metal foam: Experimental and numerical analysis. Energy Convers Manage 153: 278-287. doi: 10.1016/j.enconman.2017.09.072
    [15] Esfe MH, Bahiraei M, Hajbarati H, et al. (2020) A comprehensive review on convective heat transfer of nanofluids in porous media: Energy-related and thermohydraulic characteristics. Appl Therm Eng, 115487. doi: 10.1016/j.applthermaleng.2020.115487
    [16] Alihosseini S, Jafari A (2020) The effect of porous medium configuration on nanofluid heat transfer. Appl Nanosci 10: 895-906. doi: 10.1007/s13204-019-01192-1
    [17] Esmaeili M, Karami M, Delfani S (2020) Performance enhancement of a direct absorption solar collector using copper oxide porous foam and nanofluid. Int J Energy Res 44: 5527-5544. doi: 10.1002/er.5305
    [18] Tay NHS, Liu M, Belusko M, et al, (2017) Review on transportable phase change material in thermal energy storage systems. Renewable Sustainable Energy Rev 75: 264-277. doi: 10.1016/j.rser.2016.10.069
    [19] Liu L, Su D, Tang Y, et al. (2016) Thermal conductivity enhancement of phase change materials for thermal energy storage: A review. Renewable Sustainable Energy Rev 62: 305-317. doi: 10.1016/j.rser.2016.04.057
    [20] Lohrasbi S, Miry SZ, Gorji-Bandpy M, et al. (2017) Performance enhancement of finned heat pipe assisted latent heat thermal energy storage system in the presence of nano-enhanced H2O as phase change material. Int J Hydrogen Energy 42: 6526-6546. doi: 10.1016/j.ijhydene.2017.01.045
    [21] Bazri S, Badruddin IA, Naghavi MS, et al. (2019) An analytical and comparative study of the charging and discharging processes in a latent heat thermal storage tank for solar water heater system. Sol Energy 185: 424-438. doi: 10.1016/j.solener.2019.04.046
    [22] Naghavi MS, Ang BC, Rahmanian B, et al. (2020) On-demand dynamic performance of a thermal battery in tankless domestic solar water heating in the tropical region. Appl Therm Eng 167: 114790. doi: 10.1016/j.applthermaleng.2019.114790
    [23] Mahdi RA, Mohammed HA, Munisamy KM, et al. (2015) Review of convection heat transfer and fluid flow in porous media with nanofluid. Renewable Sustainable Energy Rev 41: 715-734. doi: 10.1016/j.rser.2014.08.040
    [24] Khanafer K, Tavakkoli F, Vafai K, et al. (2015) A critical investigation of the anomalous behavior of molten salt-based nanofluids. Int Commun Heat Mass Transfer 69: 51-58. doi: 10.1016/j.icheatmasstransfer.2015.10.002
    [25] Corcione M (2011) Empirical correlating equations for predicting the effective thermal conductivity and dynamic viscosity of nanofluids. Energy Convers Manage 52: 789-793. doi: 10.1016/j.enconman.2010.06.072
    [26] Kandelousi MS (2014) KKL correlation for simulation of nanofluid flow and heat transfer in a permeable channel. Phys Lett A 378: 3331-3339. doi: 10.1016/j.physleta.2014.09.046
    [27] Doerr M, Frommherz M (2002) Graphite (C) Classifications, Properties & Applications.
    [28] Sanchez-Coronado J, Chung DDL (2003) Thermomechanical behavior of a graphite foam. Carbon 41: 1175-1180. doi: 10.1016/S0008-6223(03)00025-3
    [29] Vafai K (2015) Handbook of porous media. Crc Press.
    [30] Karaipekli A, Biçer A, Sarı A, et al. (2017) Thermal characteristics of expanded perlite/paraffin composite phase change material with enhanced thermal conductivity using carbon nanotubes. Energy Convers Manage 134: 373-381. doi: 10.1016/j.enconman.2016.12.053
    [31] Kim JH, Jeong E, Lee YS (2015) Preparation and characterization of graphite foams. J Ind Eng Chem 32: 21-33. doi: 10.1016/j.jiec.2015.09.003
    [32] Bianco V, Manca O, Nardini S, et al. (2015) Heat Transfer Enhancement with Nanofluids.
    [33] Oosthuizen P, Naylor D (1999) Introduction to Convective Heat Transfer Analysis. WCB/McGrawHill.
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2663) PDF downloads(130) Cited by(5)

Figures and Tables

Figures(22)  /  Tables(3)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog