Next Article in Journal
A High Gain Modified Quadratic Boost DC-DC Converter with Voltage Stress Half of Output Voltage
Previous Article in Journal
A Pilot Study of whether or Not Vegetable and Fruit Juice Containing Lactobacillus paracasei Lowers Blood Lipid Levels and Oxidative Stress Markers in Thai Patients with Dyslipidemia: A Randomized Controlled Clinical Trial
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Investigation of the Effect of Square and Sinusoidal Waves Vibration Parameters on Heat Sink Forced Convective Heat Transfer Enhancement

by
Ambagaha Hewage Dona Kalpani Rasangika
1,*,
Mohammad Shakir Nasif
1,*,
William Pao
1 and
Rafat Al-Waked
2
1
Department of Mechanical Engineering, Universiti Teknologi PETRONAS, Seri Iskandar 32610, Perak, Malaysia
2
Department of Mechanical and Maintenance Engineering, German Jordanian University, Amman 11180, Jordan
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2022, 12(10), 4911; https://doi.org/10.3390/app12104911
Submission received: 11 February 2022 / Revised: 20 April 2022 / Accepted: 9 May 2022 / Published: 12 May 2022
(This article belongs to the Section Mechanical Engineering)

Abstract

:
Among numerous electronic cooling methods, a vibrating heat sink using sinusoidal wave vibration effectively enhances the heat transfer by disturbing the thermal boundary layer. However, sinusoidal wave vibration has reached its limits in enhancing heat transfer. The present study utilizes a new square wave-shaped vibration and numerically investigates the thermal performance of a heat sink subjected to sinusoidal and square waves vibration. It is found that using the square wave vibration is more beneficial to the thermal performance of the heat sink than the sinusoidal wave. The sudden impulsive motion of square wave vibration induces a higher randomness of the airflow profile and recirculation zones than the sinusoidal wave, causing the air flow to impinge directly into the fin surfaces, and further enhances the heat transfer. Furthermore, increasing the frequency and amplitude leads to a higher heat transfer enhancement. Moreover, square wave vibration achieves a 25% increase in Nusselt values compared to the nonvibrating fins and it is 11% higher than the Nusselt number recorded by the sinusoidal vibration. Consequently, Reynolds number values can be reduced by 42.2% to achieve the Nusselt number values of nonvibrating fins, potentially reducing the cooling system or fin size. This reduction may contribute to solving the challenges of electronic systems compactness.

1. Introduction

Continued development in the electronic industry has led to a significant increase in energy consumption [1,2]. A higher percentage of this energy is converted into heat [3]. In addition, reducing the size of electronic components is on the rise, making the cooling of these components smaller whilst maintaining the compactness of the system a challenge. Among numerous cooling enhancements, one mechanism is the fan-based cooling system, a well-known method due to its high performance-to-cost ratio [4]. However, the thermal performance of fan-based electronic cooling systems is saturated by the formation of the thermal boundary layer, which is found to hinder the heat transfer between the cooling fluid and heated surface [5]. In addition, increasing cooling requires increasing fan size, which compromises the compactness of electronic components.
Consequently, researchers have implemented different mechanisms in electronic cooling systems to overcome the limited thermal performance. In the literature, the most popular passive thermal enhancement techniques in electronic cooling have been investigated, such as utilizing foam-based heat sinks [6,7], and optimizing the heat sink shape, size, arrangement, and artificial roughness [8,9,10]. Various approaches by enhancing the mixing process of working fluid with vortex generators [4,11,12], vortex tubes [13], mixing promoters [14], or flow-induced vibration [15,16] have also been attempted by many researchers. However, the benefits of electronic cooling by utilizing passive technologies have been found to be inefficient for high-power dissipation electronic systems.
Therefore, active cooling methods have been utilized to enhance the electronic cooling, including the projection of air at high speed toward a heated surface using air jets [17,18,19,20] or spray [14,21], investigated to enhance the thermal performance. Furthermore, using heat pipes [22] or adding oscillating heat pipes [23] and utilizing nanofluids in conjunction with a heat sink [24,25] as well as electronic cooling systems [26,27] have shown to be promising techniques for the thermal management of electronic systems.
However, another effective technique has been used; applying external vibrations to the heated object has also been studied as a heat transfer intensification technique. To prove the concept of vibrational heat transfer augmentation, experimental and numerical studies have been performed with different geometries, including gauge wire [28], spheres [29], and cylinders [30,31,32,33]. These studies noted that forced convective heat transfer is significantly influenced by external vibration. Due to the need to reduce the size of electronic components and cooling systems, applying vibration on fins and heat sinks could be a viable option to achieve a reduction in the electronic cooling system. Thus, several studies have considered the influence of vibrational frequency and amplitude variation on heat transfer augmentation. Gururatana and Li [34] investigated the heat transfer of a vibrating 2D pin fin numerically and observed an increase in heat transfer coefficient with frequency. An analysis detailing the influence of a vibrating 2D plate-fin on the flow field was reported by Rahman and Tafti [35]. The maximum heat transfer was obtained at the highest vibrational characteristics.
Dey and Chakrborty [5] numerically investigated the heat transfer augmentation of a vibrating actuated 3D fin. The authors obtained limited enhancement with lower vibrational characteristics. However, they noted an enhanced thermal performance with an increase in vibrational characteristics. Hussain, et al. [36] experimentally studied the thermal performance of piezo-actuated cooling fins and observed a higher heat transfer with the increase in vibrational frequency. Najim, et al. [37] conducted an experimental study indicating a potential increase in the thermal performance of the vibrating heat sink at higher frequencies.
Given that the application of external vibration is beneficial to thermal performance, it may also cause an adverse effect on the reliability of the system. There are studies that have investigated the reliability of the electronic systems under the vibration that occurs during shipping, handling, and the service life of the system. The reliability assessment of electronic assemblies under shock, harmonic, and random vibration has been performed experimentally [38] and analytically [39,40] to investigate the strains and stresses of printed circuit boards (PCBs) and the solder interconnect. Moreover, finite element modeling analysis (FEA) has been performed by several studies [41,42,43] to investigate the fatigue performance and failure modes of electronic products under vibrating conditions. However, no research has been performed to investigate the reliability of electronic systems due to the application of external vibration. Thus, a special consideration on the relatability of the system may be required in designing a vibrating system.
As mentioned, vibration has been a viable method to enhance heat transfer, and further investigation has been warranted to increase the enhancement. Although previous research has focused on the effect of sinusoidal vibration on forced convective heat transfer in heated geometries, new vibration methods by introducing different vibrational wave shapes on thermal systems are needed.
The vibrating heated surface enhances heat transfer mainly by setting the fluid in motion and consequently disrupting the thermal boundary layer [14]. However, the motion of the fluid particles depends not only on the vibrational characteristics but also on the type of vibrational wave shape. The different vibrational wave shapes may have a distinct effect on the flow profile, potentially affecting the thermal performance of the heat sink. This effect has not been considered in previous studies and is worth investigating.
In fact, Zhang, et al. [44] conducted an experimental study to investigate the effect of vibration modes on the detachment of low-rank coal particles from oscillating bubbles. The authors found that different vibration modes of oscillating bubbles have a distinct effect on path changes, and applying square wave excitations leads to the highest dramatic movement path changes due to its impulsive motion. In addition, the authors claimed that the square wave-shaped vibration induces the greatest acceleration within the particles compared to other vibration modes at similar vibrational frequencies and amplitudes. On the other hand, several experimental studies have been performed by utilizing a square wave-shaped vibration using a piezo-shaker [45], vibrators [44], and piezoelectrical actuators [46]. As applying square wave-shaped vibration has resulted in the highest movement path change as compared with the sinusoidal wave, utilizing square wave vibrations in heat sinks may result in a higher heat transfer enhancement, which has not been previously utilized in electronic systems and could be considered a new implementation to further enhance heat transfer.
As seen in the vibrational wave plot (Figure 1), square wave-shaped vibration causes a sudden impulsive motion of the heated geometry, potentially causing rapid velocity changes. Therefore, introducing the new square wave-shaped vibration in comparison with sinusoidal vibration may produce a higher heat transfer enhancement and consequently further reduce the size of the cooling system. To the best of the author’s knowledge, the application of square wave-shaped vibration on a heat sink to improve the thermal performance of the heat sink is new and has not yet been investigated by previous researchers.
The present study investigates the heat transfer enhancement resulting from a vibrating heat sink under the influence of the newly introduced square wave shape and sinusoidal wave shape using ANSYS/FLUENT CFD software. The examined heat sink structure is subjected to different vibration levels. Therefore, the dependence of heat transfer enhancement on the vibration frequency and amplitude of sinusoidal and square wave shapes is investigated. The possible reduction in both Reynolds number and size of the cooling system is discussed. Furthermore, the parameterization of the Nusselt number with frequency and amplitude under sinusoidal and square wave shapes is established to identify the operating envelope of vibration characteristics for enhanced thermal performance. Outcomes from the current study are expected to enable identification of the appropriate envelope of the vibrational characteristics, thereby enhancing the heat transfer of the heat sink and consequent reduction in the electronic cooling system and heat sink size.

2. Simulation Model and Methods

2.1. Model Development

The aluminum heat sink considered in the current study is presented in Figure 2 and the size of the rectangular fluid domain is 0.305 × 0.310 × 0.037 m3. A uniformly distributed power input of 9.5 W is supplied to the bottom of the heat sink, which corresponds to 6250 W/m2, which is usually generated in data centers [47,48]. The uniform velocity is employed at the air inlet with a temperature of 25 °C, and the airflow is laminar at a Reynolds number of 1000, which corresponds to an averaged air inlet velocity of 1.46 m/s.
Mass, momentum, and energy conservation laws form the governing equations for the current study. The flow is assumed to be a laminar incompressible ideal gas with constant thermo-physical properties. The body forces are neglected, and a no-slip condition is applied on the interface between the heat sink and fluid. Thus, the governing equations can be presented as:
Continuity equation:
u x + v y + w z = 0 ,
Momentum equations:
u t + ρ u u x + ρ v u y + ρ w w z = μ 2 u P x
ρ v t + ρ u v x + ρ v v y + ρ w w z = μ 2 v P y
ρ w t + ρ u w x + ρ v w y + ρ w w z = μ 2 w P z
Energy conservation equation:
ρ c p T t + ρ c p x u T + ρ c p y v T + ρ c p z w T = k f l u i d 2 T  
Energy equation for heat sink:
2 T = 1 α   T t  
where u, v, and w are the velocities along the x, y, and z directions, respectively, ρ is the fluid density, P is the hydrostatic pressure, and μ is the fluid viscosity, which is assumed to be constant. kfluid is the thermal conductivity of the fluid, cp is the constant specific heat, T is the temperature, and α is the thermal diffusivity of the fin material (kfin/ρcp). Past research has found that a higher heat transfer enhancement is obtained with horizontally vibrating heated objects [49,50]. Hence, the sinusoidal and square wave-shaped vibration is applied to the heat sink in the Y direction (Figure 2). To simulate vibration in FLUENT, User Define Function codes (UDFs) are developed, which include the displacement equations of sinusoidal and square wave-shaped vibration, given in Equations (7) and (8), respectively.
Y = A v s i n 2 π f t  
Y = 4 A v π   n = 1 , 3 , 5 . . N 1 n s i n 2 π n f t .  
where Y is the displacement of the heat sink at any time, f is the vibrational frequency, Av is the vibrational amplitude, and t is time.
The investigation is conducted over a range of vibrating frequencies 0–100 Hz with peak-to-peak amplitudes ranging from 0 m to 0.005 m. A Reynolds number of 1000 is used for the vibrating case. Theoretical Reynolds numbers under nonvibrating conditions have been calculated using Equations (9) and (10) [51,52], which are used for comparative analyses with the vibrating heat sink.
N u = tanh 2 N u i d e a l   k f l u i d K f i n   h b   h t f i n   t f i n L f i n + 1 2 N u i d e a l k f l u i d K f i n   h b   h t f i n   t f i n L f i n + 1 N u i d e a l  
N u i d e a l = R e . b L f i n P r 2 3 + 0.664 R e . b L f i n .   P r 1 / 3 1 + 3.65 R e . b L f i n 3 1 / 3  
where Nu is the Nusselt number under the normal condition with the fin effect, Nuideal is the ideal Nusselt number, Re is the Reynolds number, b is the fin spacing, and Lfin is the length of the fins. kfin is the thermal conductivity of the fin material and Pr is the Prandtl number. h is the height of the fins and tfin is the thickness of the fin.

2.2. Computational Details

The governing equations solved by ANSYS/FLUENT CFD software, discretized by the finite-volume method and the SIMPLE algorithm, is used to obtain the pressure–velocity coupling. The spatial discretization of momentum and energy equations is obtained by using a second-order upwind scheme. The transient formulation is obtained by using a second-order implicit scheme. The criteria for convergence of the iterative solution are specified as 10−6 for the residual of continuity, momentum, and energy.
Figure 3 shows the computational domain of the current study. While the heat sink is moving with sinusoidal and square wave-shaped vibration, the mesh of the fluid domain must be rearranged. A dynamic mesh setup is used to obtain the mesh motion, and a layering scheme is used to regenerate the mesh.
A mesh independency analysis is performed to ensure that the numerical results are independent of the number of cells. As shown in Figure 4, eight grids are considered in the current study with numbers of cells of 234,556, 324,568, 435,892, 498,460, 522,764, 580,356, 784,322, and 985,662. When the number of cells is increasing from 234,556 to 522,764, the Nusselt number is found to significantly decrease. However, a further increment in the number of cells leads to less than 1% of deviation in the Nusselt number. Thus, the grid with 522764 cells is used as the optimum mesh for all developed models.
The quality of a mesh used in the simulation is indicated by skewness and aspect ratio values. The maximum skewness and aspect ratio of the mesh for the present study are 0.62 and 6.12, respectively, meeting the requirements of acceptable mesh quality.
In addition, a time step independency investigation is carried out using three different time steps Δt = τ/8, τ/16, and τ/32 at a frequency of 100 Hz. The variations in the area-averaged Nusselt number of the heat sink with cycle time are shown in Figure 5. Reducing the time step (Δt) from τ/16 to τ/32 s gives less than 1% of variation in the peak-to-peak variation in the area-averaged Nusselt number. Here, τ is the cycle time of each vibration case. Thus, the τ/16 s time step is used in all simulated cases.
Validation of the CFD modeling is performed in terms of the sinusoidal vibrating heat sink, which is simulated under similar experimental conditions to Najim, et al. [37]. Numerical computations are performed under a range of vibrational frequencies, while the amplitude, supplied power, and air velocity are kept constant at 0.005 m, 17 W, and 0.6 m/s, respectively. As shown in Figure 6, a reasonable agreement is achieved in the validation against Najim, et al.’s [37] experimental measurements with a deviation of up to 8.3%. Therefore, CFD modeling in this study can be utilized with confidence to numerically analyze the heat transfer enhancement of the vibrating heat sink.
However, it can be seen that the average Nusselt number of the validation study is consistently lower than that of the experimental finding for the entire range of vibrational frequency. Hussain, et al. [36] performed numerical and experimental investigations and validated their numerical results against the experiments they performed. Their validation results were presented in terms of temperature of the fin, and they showed that their experimental temperature graph was consistently lower than the numerically obtained temperature. In other words, when the temperature graph corresponds to a higher Nusselt number, it can be concluded that the results obtained by the [36] is in agreement in terms of the trend with our validation results.

3. Results

CFD modelling is initially performed for the static heat sink (shown in Figure 2) at a Reynolds number of 1000. Then, Re = 1000 is used in the investigation of the vibration and amplitude variation effect on heat transfer. The airflow and vibrating directions are shown in the black arrows in the top view of the heat sink (Figure 7).
The vibration of the heat sink improves heat transfer by inducing randomness of the airflow profile and recirculation in the flow field. Thus, investigating the velocity profile shed from the vibrating heat sink is important. The effect of using the sinusoidal and square wave-shaped vibrations on the flow profile is expected to be different. Investigating the transient development of the velocity profile under each type of vibration compared with the velocity profile of the stationary heat sink is therefore important. The velocity profile of the nonvibrating heat sink and the transient development of air velocity vectors on a plane passing through the top section of the heat sink at different vibration cycle stages for sinusoidal and square wave shapes are presented. The air velocity vectors shown in Figure 8b–i are obtained with sinusoidal and square wave-shaped vibrations at frequencies of 100 Hz at a constant peak-to-peak amplitude of 0.003 m. In each figure, the two fins of the vibrating heat sink are shown as red surfaces.
After performing the validation of the model, the heat sink shown in Figure 2 is modeled under normal and vibration conditions under the range of vibrational frequencies and amplitudes for both sinusoidal and square wave shapes.
At the beginning of the sinusoidal motion, the heat sink moves to a positive Y direction. Figure 8b shows that, at the end of the first quarter of the cycle, the vibrating fins induce airflow recirculation near the trailing edges of the fins, whereas the flow between the two fins moves upward toward fin 1 following the vibration displacement direction and hitting the surface of fin 1.
At the end of the second and third quarter of the vibration cycle, Figure 8c,d show that the recirculation of air at the trailing edges of the fins continues. However, the airflow in the middle channel between the two fins changes its direction and hits the interior surfaces, especially the interior surface of fin 2. This induced flow change in direction becomes more noticeable at the end of the third quarter of the cycle. When the vibration continues at the end of the cycle, it can be seen in Figure 8e that the flow direction and recirculation become similar to Figure 8c.
At the beginning of the square wave-shaped motion, the heat sink reaches its maximum positive amplitude almost instantly. As a result, Figure 8f shows that the square wave excitation of the fins induces strong air recirculation near the trailing edges of fins compared to the sinusoidal motion at the end of the first quarter of the vibration cycle.
With increasing time, the heat sink moves to the end of the second quarter of the cycle (Figure 8g). The impulsive motion of the heat sink gives a higher air velocity magnitude to the air particles. Thus, the fluid adjacent to it gains noticeable movement path changes. This leads to the formation of strong clockwise and counter-clockwise recirculation zones near the leading and trailing edges of both fins.
Afterward, the heat sink moves in the negative y-direction continuously at the start of the third quarter of the vibration cycle. As indicated in Figure 8h, the flow circulation zones are found to be similar to those mentioned earlier at the end of the second quarter. However, the strength of the circulation zones near the leading and trailing edges of the fins decreases gradually at the end of the third quarter of the vibration cycle.
As vibration continuous, the heat sink moves to the positive y-direction at the end of the vibration cycle (Figure 8i), and the flow direction and recirculation zones are found to be opposite to those, as shown in Figure 8g.
Compared with the nonvibrating heat sink shown in Figure 8a, air velocity vectors under the vibration effect are not straight streamlines. When the heat sink vibrates with some frequency and amplitude, the airflow direction continuously changes toward the fins' surface, which is expected to improve the fins’ thermal performance. When the heat sink vibrates under a square wave shape, strong recirculation zones and secondary flows are formed around both fins at each quadrant of the vibration cycle owing to the sudden movement of the heat sink, causing a noticeable velocity increase in the flow field. By comparing air velocity vectors between the square shape and sinusoidal vibration, the velocity magnitude and strength of the recirculation zones are higher under square-shaped vibration. This phenomenon induces more disruption to the developed thermal boundary layer, thereby further enhancing the thermal performance of the heat sink.
Figure 9 shows the temperature contours of the nonvibrating heat sink compared with other cases of increasing frequency and amplitude at each type of vibrational wave shape. To avoid repetition, only one representing case for vibration frequency and amplitude under each wave shape is selected to show the effect of vibration on the thermal boundary layer. As seen in Figure 8a, the nonvibrating heat sink has a uniform airflow distribution in the direction of the flow field. Thus, the temperature gradient is high between the air stream and near the fin surfaces (Figure 9a), showing that the effect of the heat resistance of the thermal boundary layer causes and limits the cooling effect.
However, when vibration is applied, it causes a rapid change in the air velocity distribution and randomness and recirculation, causing the fluid to move directly toward the heated fin surfaces, as shown in Figure 8. This phenomenon consequently ruptures the thermal boundary layer. As seen in Figure 9b,c, applying vibration on the heat sink induces a disturbance to the continuous thermal wake. It leads to a widespread decrease in air temperature near the fin surfaces, which tends to have a higher cooling effect than the static fin. This randomness and recirculation are more noticeable with a higher air velocity magnitude when a square-shaped vibration is used. Figure 8f–i show that using a square wave-shaped vibration induces stronger recirculation zones, leading to a higher mixing effect within the flow field. Thus, the overall spread of thermal disturbance and lower air temperature becomes more noticeable when square wave excitation is applied compared with the sinusoidal wave shape.
The cooling effect of the square shape vibration can be seen in Figure 9c, where it causes a widespread low-temperature region at both fin surfaces, indicating a further cooling effect.
When the peak-to-peak amplitude is increased from 0.003 to 0.005 m at a constant frequency of 100 Hz, the recirculation zones and randomness of the airflow profile increase as the peak-to-peak amplitude increases. As a result, the thermal boundary layer is disturbed further, resulting in increased heat transfer. As seen in Figure 9b–e, increasing the amplitude leads to a widespread decrease in air temperature near the fin surfaces, tending to have a stronger cooling effect. However, as seen in Figure 9e, increasing the amplitude with square wave-shaped vibration leads to a higher decrease in air temperature near the fin surfaces due to stronger flow circulations, leading to a higher mixing effect within the flow field.
The discussion shows that vibrational frequency, amplitude, and wave shape are three important parameters that induce air mixing and rupture the thermal boundary layer. To evaluate the thermal effectiveness and consequent cooling enhancement of the system by varying those parameters, time-averaged Nusselt numbers are used as an indicator. As a result, the parameterization of time-averaged Nusselt number with vibration frequency and amplitude for each wave shape is performed and presented in Figure 10. The parametrization helps to more simply estimate the heat transfer enhancement and consequent reduction in the cooling system based on the vibrational parameters under each wave shape. The parametrization also aids in determining the optimum scenario of parameter combinations that provide the best thermal performance of the system.
Moreover, according to the Nusselt number value obtained from the modeling, Reynolds numbers are calculated from Equations (9) and (10), which gives the value of the Reynolds number under no vibrating conditions, and those values are presented at each data point in Figure 10.
As shown in Figure 10a below, when the vibrational frequency is below 30 Hz, the Nusselt number almost remains unchanged with each type of wave shape compared with the static heat sink. When the frequency is lower than the threshold value, the energy of vibration propagates through the solid interface, and the fluid may not be adequate to induce the strong recirculation zones or secondary flows within the flow field. Thus, applied vibration is insufficient to disturb the developed thermal resistance caused by the thermal boundary layer. Notably, this threshold is observed to be independent of the vibration wave shape. Thus, using sinusoidal or square wave-shaped vibration on the heat sink with a lower frequency (f < 30 Hz) does not influence the thermal performance of the system. However, after exceeding the threshold frequency, Nusselt values show a steady rise with increasing frequency at each vibrational wave shape.
This result can be attributed to the vibration effect, where a lower vibration frequency may not induce sufficient air perturbation, which cannot disrupt the thermal boundary layer. However, when the frequency increases beyond the threshold, it causes a higher agitation and air perturbation, causing disruption to the thermal boundary layer and consequently enhancing heat transfer.
As expected, the square wave excitation of the heat sink causes the Nusselt number to shift to an elevated value compared with the sinusoidal motion because square wave excitation of the heat sink attains the maximum amplitude instantly, which induces a strong impulsive motion, leading to a significant velocity change and stronger recirculation zones, thereby enhancing the cooling effect. Moreover, the relative Nusselt number enhancement difference in square and sinusoidal wave motion is improved with increased frequency. Compared with the static heat sink, the Nusselt number is improved up to 20% with square wave-shaped vibration, whereas a 12% increase in Nusselt number values is recorded when the sinusoidal frequency is used.
To investigate the effect of varying the amplitude at a constant frequency on heat transfer enhancement, a 100 Hz frequency is selected to show the effect on the Nusselt number. As seen in Figure 10b above, an increase in peak-to-peak amplitude at a frequency of 100 Hz leads to an increase in time-averaged Nusselt number in both types of excitations. The increase in the peak-to-peak amplitude at higher frequency imparts a higher fin displacement, resulting in a greater recirculation within the flow field and improved heat transfer.
However, the heat transfer enhancement is found to be higher with square wave excitations. Once the amplitude of the square-shaped vibration is increased, the impulse motion of the vibrating fins combined with the higher amplitude causes a higher air perturbation due to the impulsive nature of this vibration, resulting in a rapid change in airflow field direction near the fin surface, further enhancing the heat transfer. When the amplitude is increased from 0 to 0.005 m at 100 Hz, a 14.2% increase in Nusselt number values is recorded when the sinusoidal vibration is used. However, at 100 Hz frequency and 0.005 m amplitude, the Nusselt number magnitude increases by up to 25% with square wave-shaped vibration, providing the system with the highest thermal performance over the evaluated parameter range.
From Figure 10, the results show that increasing amplitude and frequency lead to a higher thermal performance of the vibrating heat sink. Therefore, it is also necessary to present the heat transfer enhancement in terms of the ratio of the Nusselt number of vibrating and nonvibrating conditions at each vibration level, which could be used as an indicator to perform the comparative analysis of the increase in heat transfer enhancement under vibrating conditions:
N u s s e l t   R a t i o = N u s s e l t   n u m b e r   w i t h   v i b r a t i n g   c o n d i t i o n s   N u s s e l t   n u m b e r   w i t h   n o n v i b r a t i n g   c o n d i t i o n s  
Figure 11a,b show the variation in Nusselt ratio with vibrational frequency and amplitude under sinusoidal and square wave shapes. Due to the insignificant heat transfer enhancement with a lower frequency (f < 30 Hz), the Nusselt ratio is found to be one at lower vibrational frequencies. However, after exceeding the critical frequency, the Nusselt ratio is found to steadily increase with increasing frequency. As expected, a higher enhancement of the Nusselt ratio is obtained with square wave-shaped excitation as compared to the sinusoidal vibration at each vibration frequency. As compared to the nonvibrating heat sink, 12% and 20% enhancements in Nusselt number ratios are obtainable by applying sinusoidal and square wave-shaped excitation, respectively, at f = 100 Hz and Amp = 0.003 m.
With the increase in peak-to-peak amplitude at constant frequency of 100 Hz, the Nusselt ratio also increases. As expected, an increase in Nusselt ratio is found to be higher with square wave-shaped vibration. As compared to the nonvibrating heat sink, 14.2 and 25% enhancements in Nusselt number ratios are obtained by utilizing the sinusoidal and square wave-shaped vibration, respectively, with the conventional heat sink at f = 100 Hz and Amp = 0.005 m.
Given that the parametric study shows the possibility of heat transfer intensification of the vibrating heat sink under sinusoidal and square wave shapes, this method can be used to reduce the size of the conventional cooling system. The Reynolds number values required under normal conditions (no vibration) to achieve the Nusselt number value under vibration conditions are presented in Figure 10 at each frequency and amplitude. Using the presented Reynolds number values, the percentage reduction in Reynolds number (PRRe) is calculated. It is used as an indicator to perform the comparative analysis of the reduction in Reynolds number under vibrating conditions. In this regard, investigating PRRe eventually estimates the possible reduction in the size of the cooling system based on the vibrational characteristics at each wave shape.
P R R e = 1 R e y n o l d s   n u m b e r   u n d e r   t h e   v i b r a t i n g   c o n d i t i o n   R e y n o l d s   n u m b e r   u n d e r   n o n v i b r a t i n g   c o n d i t i o n s   ×   100
Figure 12a,b show the variations in PRRe with vibrational frequency and amplitude at both sinusoidal and square wave shapes. Given that no heat transfer enhancement with a lower frequency range (f < 30 Hz) occurs, no opportunity for reductions in Reynolds number arises in this frequency range. However, after exceeding the critical frequency (f = 30 Hz), PRRe shows a steady rise with an increase in frequency for each type of wave shape. Given that the square wave-shaped vibration shows a higher heat transfer enhancement than the sinusoidal vibration at each frequency level, a higher percentage reduction in Reynolds number is obtained with increasing frequency under square wave-shaped excitation. As observed from Figure 12a, 25% and 36.5% reductions in Reynolds numbers are obtainable by applying sinusoidal and square wave-shaped excitation, respectively, at f = 100 Hz and Amp = 0.003 m compared to the conventional cooling system.
When the peak-to-peak amplitude increases at a constant frequency of 100 Hz, PRRe also increases. As expected, an increase in PRRe is found to be higher with square wave-shaped vibration. Figure 12b shows that PRRe rises to 28.8% with sinusoidal wave-shaped and 42.2% with square wave-shaped vibration at f = 100 Hz and Amp = 0.005 m. This result indicates that using vibration on heat sink structures can achieve a similar thermal performance at lower Reynolds numbers than can the conventional Reynolds numbers required to achieve the same thermal performance under nonvibrating conditions. Moreover, applying square wave-shaped vibration helps to achieve a further reduction in the size of the cooling system compared to sinusoidal vibration at every investigated frequency (f > 30 Hz) and amplitude level. Thus, the square-shaped vibration can achieve similar thermal performance with a smaller cooling system than a static heat sink with a much larger cooling system. Thus, this enhancement method might be well suited to the electronic industry because it is coherent with the challenges of compactness experienced in most modern electronic systems.

4. Discussion

Conclusively, the parameterization of Nusselt number with frequency and amplitude under sinusoidal and square wave shapes is established in the current study, helping to identify the appropriate operating envelope of vibration characteristics for enhanced thermal performance. Increasing the frequency or amplitude within the appropriate operating envelope region increases the heat sink thermal performance at each vibrational wave shape. However, square wave-shaped vibration is found to be more effective in the enhancement of heat transfer than sinusoidal vibration. Thus, applying square wave-shaped vibration on the forced convective cooling system is more beneficial to its thermal performance and further reduces the cooling system size compared with sinusoidal vibration at similar vibrational frequencies and amplitudes. By using square wave-shaped vibration at f = 100 Hz and 0.005 m amplitudes, the maximum enhancement of the Nusselt number is 25% compared with the nonvibrating heat sink, which is 11% higher than the Nusselt number recorded by the sinusoidal vibrating heat sink. This finding results in a 42.2% reduction in Reynolds number compared with the Reynolds number of nonvibrating cooling systems, which is 13.4% higher than the reduction in Reynolds number by the sinusoidal vibrating heat sink at similar vibrational frequencies and amplitudes.
Researchers have been continuously investigating new methods to enhance the cooling of electronic systems. Zhuang, et al. [53] used a fan-integrated microchannel heat sink with rhombus fractal-like units. They achieved an 18% increase in Nusselt number compared with a conventional parallel microchannel heat sink. Dey and Chakrborty [5] applied sinusoidal vibrations on a piezo-actuated fin. They achieved up to a 25.5% improvement in Nusselt number compared with a nonvibrating fin. Freegah, et al. [54] used plate-fin heat sinks with corrugated half-round pins in a vertical arrangement subject to parallel flow and recorded an approximately 34.5% increase in Nusselt number values compared with the traditional plate-fin heat sink. Masip, et al. [55] used an impinging air jet with a crossflow configuration. They recorded an up-to-40% increase in the average Nusselt number compared with the conventional cross-flow configuration. Moon, et al. [56] applied a pulsating flow through heated blocks. They recorded increases of up to approximately 23% in Nusselt number values compared with a nonpulsating flow.
As mentioned, using square-shaped wave vibration achieved an up-to-25% increase in Nusselt number and a 42.2% reduction in Reynolds number compared with nonvibrating fins. This result shows that square shape-waved vibration resulted in a higher cooling effect compared with other methods, making it achieve further reductions in the cooling system size of electronic components.
It is worth mentioning that the results presented in this research could be assumed applicable to heat sinks with more fins under the same heat flux and air velocity, because past research reported by Adhikari, et al. [57] found that the Nusselt number does not change with the increase in the number of static fins. They also reported that when they varied the number of fins, the flow and thermal contours indicated very small differences between the heat sink fin channels. Since the finding of [57] that the Nusselt number remains unchanged when the number of static fins increases, and in this research, the application of vibration on fins is equal, it is expected that the results presented in this research will be the same when number of fins increases under the same heat flux and air velocity. However, further research is required to confirm this assumption.
It is also important to mention that special consideration may be required in designing a vibrating system. Therefore, future research could be performed to investigate the reliability of the system demonstrating the mechanical stress factors and fatigue life due to the supplied vibration.

5. Conclusions

Numerical simulations are conducted to investigate the effect of vibration on the heat transfer performance of a heat sink. A parametric study is performed by using vibrational frequency, amplitude, and wave shapes to quantify the dependence of thermal performance on each parameter. When the vibrational frequency is lower than the critical frequency (f < 30 Hz), the thermal performance of the heat sink is not influenced by the vibration. However, after exceeding the critical frequency, an increase in amplitude and frequency leads to a higher enhancement in thermal performance at each vibrational wave shape. This result enables the identification of the appropriate operating envelope of the vibrational characteristics, thereby enhancing the heat transfer of the heat sink. Using square wave-shaped excitations provides more benefit to heat transfer than does sinusoidal vibration at every frequency and amplitude level. A maximum increase of 14.2% in Nusselt number values is recorded with sinusoidal vibration, while a maximum enhancement of 25% is obtained with square wave-shaped vibration. Thus, using a square wave-shaped vibration on the heat sink at the frequency of 100 Hz and amplitude of 0.005 m is found to achieve the highest thermal performance of the system over the range of the evaluated parameters. The process results in a 42.2% reduction in Reynolds number, making using the square wave-shaped vibration beneficial in reducing the conventional cooling system and achieve the compactness needed in modern electronic systems.
However, the current study is performed at a single Reynolds number and power input. Future research could be performed on a range of Reynolds numbers and power inputs. Moreover, this current study could be extended with experimental work and investigating the energy used to vibrate the fins as compared with heat transfer enhancement. This research could also be extended to investigate the reliability of the system due to the vibration of the heat sink.

Author Contributions

Conceptualization, A.H.D.K.R. and M.S.N.; methodology, A.H.D.K.R. and M.S.N.; software, M.S.N.; validation, A.H.D.K.R. and M.S.N.; formal analysis, A.H.D.K.R. and M.S.N.; investigation, A.H.D.K.R. and M.S.N.; resources, M.S.N.; data curation, A.H.D.K.R., M.S.N., W.P., and R.A.-W.; writing—original draft preparation, A.H.D.K.R.; writing—review and editing, M.S.N., W.P., and R.A.-W.; visualization, A.H.D.K.R. and M.S.N.; supervision, M.S.N., W.P., and R.A.-W.; project administration, M.S.N.; funding acquisition, M.S.N. All authors have read and agreed to the published version of the manuscript.

Funding

The Ministry of Higher Education, Malaysia (MoHE) through The Fundamental Research Grant (FRGS)—grant number FRGS/1/2018/TK03/UTP/02/7.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding authors. The data are not publicly available, due to privacy.

Acknowledgments

The authors would like to thank the Ministry of Higher Education (MOHE), Malaysia for providing financial assistance under the Fundamental Research grant Scheme (FRGS/1/2018/TK03/UTP/02/7) and Universiti Teknologi PETRONAS for providing the required facilities to conduct this research work.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature and Abbreviations

Nomenclature
ATotal heat transfer surface area of the heat sink (m2)
AmpPeak-to-peak amplitude (m)
AvVibrational amplitude (Amp/2) (m)
bFin spacing (m)
cpConstant specific heat (J/(kg K))
fVibrational frequency (Hz)
hHeight of the fins (m)
kfinThermal conductivity of the fin material (W/(m K))
kfluidThermal conductivity of the fluid (W/(m K))
LfinLength of the fins (m)
NuNusselt number of heat sink at the normal condition with fin effect
N u ¯ Time-averaged Nusselt number
NuAArea-averaged Nusselt number
NuidealIdeal Nusselt number
PHydrostatic pressure (Pa)
PrPrandtl number
ReReynolds number (ρUb/μ)
TTemperature (K)
tTime (s)
tfinThickness of the fin (m)
UInlet velocity (m/s)
u, v, wVelocities along the x, y, and z directions, respectively (m/s)
x, y, zCartesian coordinates (m)
YDisplacement at any time (m)
Greek Symbols
μFluid viscosity (Pa.s)
ρFluid density (Kg/m3)
τCycle time (s)
ΔtTime step (s)
αThermal diffusivity of the fin material (m2/s)
Abbreviations
PRRePercentage reduction in Reynolds number

References

  1. Gururatana, S. Heat transfer augmentation for electronic cooling. Am. J. Appl. Sci. 2012, 9, 436–439. [Google Scholar] [CrossRef] [Green Version]
  2. Fan, S.; Duan, F. A review of two-phase submerged boiling in thermal management of electronic cooling. Int. J. Heat Mass Transf. 2020, 150, 119324. [Google Scholar] [CrossRef]
  3. He, Z.; Xi, H.; Ding, T.; Wang, J.; Li, Z. Energy efficiency optimization of an integrated heat pipe cooling system in data center based on genetic algorithm. Appl. Therm. Eng. 2021, 182, 115800. [Google Scholar] [CrossRef]
  4. Li, H.Y.; Chen, C.L.; Chao, S.M.; Liang, G.F. Enhancing heat transfer in a plate-fin heat sink using delta winglet vortex generators. Int. J. Heat Mass Transf. 2013, 67, 666–677. [Google Scholar] [CrossRef]
  5. Dey, S.; Chakrborty, D. Enhancement of convective cooling using oscillating fins. Int. Commun. Heat Mass Transf. 2009, 36, 508–512. [Google Scholar] [CrossRef]
  6. Ur Rehman, T.; Ali, H.M.; Saieed, A.; Pao, W.; Ali, M. Copper foam/PCMs based heat sinks: An experimental study for electronic cooling systems. Int. J. Heat Mass Transf. 2018, 127, 381–393. [Google Scholar] [CrossRef]
  7. Ur Rehman, T.; Ali, H.M. Thermal performance analysis of metallic foam-based heat sinks embedded with RT-54HC paraffin: An experimental investigation for electronic cooling. J. Therm. Anal. Calorim. 2020, 140, 979–990. [Google Scholar] [CrossRef]
  8. Leon, O.; de Mey, G.; Dick, E.; Vierendeels, J. Comparison between the standard and staggered layout for cooling fins in forced convection cooling. J. Electron. Packag. Trans. ASME 2003, 125, 442–446. [Google Scholar] [CrossRef]
  9. Yang, K.S.; Li, S.L.; Chen, I.Y.; Chien, K.H.; Hu, R.; Wang, C.C. An experimental investigation of air cooling thermal module using various enhancements at low Reynolds number region. Int. J. Heat Mass Transf. 2010, 53, 5675–5681. [Google Scholar] [CrossRef]
  10. Leon, O.; de Mey, G.; Dick, E. Study of the optimal layout of cooling fins in forced convection cooling. Microelectron. Reliab. 2002, 42, 1101–1111. [Google Scholar] [CrossRef]
  11. Samadifar, M.; Toghraie, D. Numerical simulation of heat transfer enhancement in a plate-fin heat exchanger using a new type of vortex generators. Appl. Therm. Eng. 2018, 133, 671–681. [Google Scholar] [CrossRef]
  12. Cheng, J.C.; Tsay, Y.L.; Liu, C.T.; Chang, S. Heat transfer enhancement of microchannel heat sink with longitudinal vortex generators and bypass jet flow. Numer. Heat Transf. Part A Appl. 2020, 77, 807–819. [Google Scholar] [CrossRef]
  13. Moraveji, A.; Toghraie, D. Computational fluid dynamics simulation of heat transfer and fluid flow characteristics in a vortex tube by considering the various parameters. Int. J. Heat Mass Transf. 2017, 113, 432–443. [Google Scholar] [CrossRef]
  14. Léal, L.; Miscevic, M.; Lavieille, P.; Amokrane, M.; Pigache, F.; Topin, F.; Nogarède, B.; Tadrist, L. An overview of heat transfer enhancement methods and new perspectives: Focus on active methods using electroactive materials. Int. J. Heat Mass Transf. 2013, 61, 505–524. [Google Scholar] [CrossRef]
  15. Go, J.S.; Kim, S.J.; Lim, G.; Yun, H.; Lee, J.; Song, I.; Pak, Y.E. Heat transfer enhancement using flow-induced vibration of a microfin array. Sens. Actuators A Phys. 2001, 90, 232–239. [Google Scholar] [CrossRef]
  16. Sun, X.; Ye, Z.; Li, J.; Wen, K.; Tian, H. Forced convection heat transfer from a circular cylinder with a flexible fin. Int. J. Heat Mass Transf. 2019, 128, 319–334. [Google Scholar] [CrossRef]
  17. Mahalingam, R.; Glezer, A. Air cooled heat sinks integrated with synthetic jets. In Proceedings of the ITherm 2002. Eighth Intersociety Conference on Thermal and Thermomechanical Phenomena in Electronic Systems, San Diego, CA, USA, 30 May–1 June 2002. [Google Scholar] [CrossRef]
  18. Krishna, S.E.B.; Pradesh, A.; Prasad, K.L.; Dist, K.; Pradesh, A.; Kartheek, G.; Dist, K.; Pradesh, A.; Dist, K.; Pradesh, A. Experimental Analysis of Jet Impingement on Aluminium Heat Sink. Int. J. Heat Mass Transf. 2018, 9, 1129–1140. [Google Scholar]
  19. Kwon, B.; Foulkes, T.; Yang, T.; Miljkovic, N.; King, W.P. Air Jet Impingement Cooling of Electronic Devices Using Additively Manufactured Nozzles. IEEE Trans. Compon. Packag. Manuf. Technol. 2020, 10, 220–229. [Google Scholar] [CrossRef]
  20. Sansoucy, E.; Oosthuizen, P.H.; Refai-Ahmed, G. An experimental study of the enhancement of air-cooling limits for telecom/datacom heat sink applications using an impinging air jet. J. Electron. Packag. Trans. ASME 2006, 128, 166–171. [Google Scholar] [CrossRef]
  21. Wang, J.X.; Guo, W.; Xiong, K.; Wang, S.N. Review of aerospace-oriented spray cooling technology. Prog. Aerosp. Sci. 2020, 116, 100635. [Google Scholar] [CrossRef]
  22. Tang, H.; Tang, Y.; Wan, Z.; Li, J.; Yuan, W.; Lu, L.; Li, Y.; Tang, K. Review of applications and developments of ultra-thin micro heat pipes for electronic cooling. Appl. Energy 2018, 223, 383–400. [Google Scholar] [CrossRef]
  23. Ling, Y.Z.; Zhang, X.S.; Wang, F.; She, X.H. Performance study of phase change materials coupled with three-dimensional oscillating heat pipes with different structures for electronic cooling. Renew. Energy 2020, 154, 636–649. [Google Scholar] [CrossRef]
  24. Arasteh, H.; Mashayekhi, R.; Toghraie, D.; Karimipour, A.; Bahiraei, M.; Rahbari, A. Optimal arrangements of a heat sink partially filled with multilayered porous media employing hybrid nanofluid. J. Therm. Anal. Calorim. 2019, 137, 1045–1058. [Google Scholar] [CrossRef]
  25. Coşkun, T.; Çetkin, E. Heat Transfer Enhancement in a Microchannel Heat Sink: Nanofluids and/or Micro Pin Fins. Heat Transf. Eng. 2020, 41, 1818–1828. [Google Scholar] [CrossRef]
  26. Colangelo, G.; Favale, E.; Milanese, M.; de Risi, A.; Laforgia, D. Cooling of electronic devices: Nanofluids contribution. Appl. Therm. Eng. 2017, 127, 421–435. [Google Scholar] [CrossRef]
  27. Bahiraei, M.; Heshmatian, S. Electronics cooling with nanofluids: A critical review. Energy Convers. Manag. 2018, 172, 438–456. [Google Scholar] [CrossRef]
  28. Faircloth, J.M.; Schaetzle, W.J. Effect of vibration on heat transfer for flow normal to a cylinder. J. Heat Transfer. 1969, 91, 140–144. [Google Scholar] [CrossRef]
  29. Baxi, C.B.; Ramachandran, A. Effect of vibration on heat transfer from spheres. J. Heat Transfer. 1969, 91, 337–343. [Google Scholar] [CrossRef]
  30. Shalaby, M.A.; Elnegiry, E.A. Forced Convection Heat Transfer From Oscillating Horizontal Cylinder. Mansoura Eng. J. 2006, 28, 1–4. [Google Scholar]
  31. Pottebaum, T.S.; Gharib, M. Using oscillations to enhance heat transfer for a circular cylinder. Int. J. Heat Mass Transf. 2006, 49, 3190–3210. [Google Scholar] [CrossRef]
  32. Gau, C.; Wu, J.M.; Liang, C.Y. Heat transfer enhancement and vortex flow structure over a heated cylinder oscillating in the crossflow direction. J. Heat Transfer. 1999, 121, 789–795. [Google Scholar] [CrossRef]
  33. Karanth, D.; Rankin, G.W. A finite difference calculation of forced convective heat transfer from an oscillating cylinder. Int. J. Heat Mass Transf. 1994, 37, 1619–1630. [Google Scholar] [CrossRef]
  34. Gururatana, S.; Li, X. Heat transfer enhancement of small scale heat sinks using vibrating pin fin. Am. J. Sci. 2013, 10, 801–809. [Google Scholar] [CrossRef] [Green Version]
  35. Rahman, A.; Tafti, D. Characterization of heat transfer enhancement for an oscillating flat. Int. J. Heat Mass Transf. 2019, 147, 119001. [Google Scholar] [CrossRef]
  36. Hussain, I.Y.; Amori, K. Theoretical and Experimental Investigation of Heat Transfer Enhancement for Hot Base Using Oscillating Fins. Iraqi J. Mech. Mater. Eng. 2015, 15, 667–685. [Google Scholar]
  37. Najim, R.; Wahib, J.H.; Jalil, S.M.; Ibrahim, M. Experimental Study of the Effect of Vertical Oscillation on Forced Convection Coefficient from Vertical Channel. Anbar J. Sci. Eng. 2013, 260, 1–13. Available online: https://www.researchgate.net/publication/303565440_9_Experimental_Study_of_the_Effect_of_Vertical_Oscillation_on_Forced_Convection_Heat_Transfer_Coefficient_of_Vertical_Channel (accessed on 8 February 2022).
  38. Roucou, R.; Zaal, J.J.M.; Jalink, J.; de Heus, R.; Rongen, R. Effect of Environmental and Testing Conditions on Board Level Vibration. In Proceedings of the 2016 IEEE 66th Electronic Components and Technology Conference, Las Vegas, NA, USA, 31 May–3 June 2016. [Google Scholar] [CrossRef]
  39. Gharaibeh, M.A.; Pitarresi, J.M. Random vibration fatigue life analysis of electronic packages by analytical solutions and Taguchi method. Microelectron. Reliab. 2019, 102, 113475. [Google Scholar] [CrossRef]
  40. Gharaibeh, M. Reliability assessment of electronic assemblies under vibration by statistical factorial analysis approach. Solder. Surf. Mt. Technol. 2018, 30, 171–181. [Google Scholar] [CrossRef]
  41. Qin, S.; Li, Z.; Chen, X.; Shen, H. Comparing and modifying estimation methods of fatigue life for PCBA under random vibration loading by finite element analysis. In Proceedings of the 2015 Prognostics and System Health Management Conference, Beijing, China, 21–23 October 2015. [Google Scholar] [CrossRef]
  42. Xu, F.; Li, C.R.; Jiang, T.M.; Zhang, D.P. Fatigue Life Prediction for PBGA under Random Vibration Using Updated Finite Element Models. Exp. Technol. 2016, 40, 1421–1435. [Google Scholar] [CrossRef] [Green Version]
  43. Su, Q.T.; Gharaibeh, M.A.; Stewart, A.J.; Pitarresi, J.M.; Anselm, M.K. Accelerated Vibration Reliability Testing of Electronic Assemblies Using Sine Dwell with Resonance Tracking. J. Electron. Packag. Trans. ASME 2018, 140, 119–125. [Google Scholar] [CrossRef]
  44. Zhang, Y.; Xing, Y.; Ding, S.; Cao, Y.; Gui, X. Effect of vibration mode on detachment of low-rank coal particle from oscillating bubble. Powder Technol. 2019, 356, 880–883. [Google Scholar] [CrossRef]
  45. Thomas, B.P.; Pillai, S.A.; Narayanamurthy, C.S. Investigation on vibration excitation of debonded sandwich structures using time-average digital holography. Appl. Opt. 2017, 56, F7. [Google Scholar] [CrossRef] [PubMed]
  46. Razfar, M.R.; Behroozfar, A.; Ni, J. Study of the effects of tool longitudinal oscillation on the machining speed of electrochemical discharge drilling of glass. Precis. Eng. 2014, 38, 885–892. [Google Scholar] [CrossRef]
  47. Ebrahimi, K.; Jones, G.F.; Fleischer, A.S. A review of data center cooling technology, operating conditions and the corresponding low-grade waste heat recovery opportunities. Renew. Sustain. Energy Rev. 2014, 31, 622–638. [Google Scholar] [CrossRef]
  48. Siriwardana, J.; Jayasekara, S.; Halgamuge, S.K. Potential of air-side economizers for data center cooling: A case study for key Australian cities. Appl. Energy 2013, 104, 207–219. [Google Scholar] [CrossRef]
  49. Shokouhmand, H.; Abadi, S.M.A.N.R.; Jafari, A. The effect of the horizontal vibrations on natural heat transfer from an isothermal array of cylinders. Int. J. Mech. Mater. Des. 2011, 7, 313–326. [Google Scholar] [CrossRef]
  50. Shokouhmand, H.; Abadi, S.M.A.N.R. Finite element analysis of natural heat transfer from an isothermal array of cylinders in presence of vertical oscillations. Heat Mass Transf. Und Stoffuebertragung 2010, 46, 891–902. [Google Scholar] [CrossRef]
  51. Teertstra, P.; Yovanovich, M.M.; Culham, J.R. Analytical forced convection modeling of plate fin heat sinks. J. Electron. Manuf. 2000, 10, 253–261. [Google Scholar] [CrossRef]
  52. Kim, S.J.; Kim, D.K.; Oh, H.H. Comparison of fluid flow and thermal characteristics of plate-fin and pin-fin heat sinks subject to a parallel flow. Heat Transf. Eng. 2008, 29, 169–177. [Google Scholar] [CrossRef] [Green Version]
  53. Zhuang, D.; Yang, Y.; Ding, G.; Du, X.; Hu, Z. Optimization of Microchannel Heat Sink with Rhombus Fractal-like Units for Electronic Chip Cooling. Int. J. Refrig. 2020, 116, 108–118. [Google Scholar] [CrossRef]
  54. Freegah, B.; Hussain, A.A.; Falih, A.H.; Towsyfyan, H. CFD analysis of heat transfer enhancement in plate-fin heat sinks with fillet profile: Investigation of new designs. Therm. Sci. Eng. Prog. 2020, 17, 100458. [Google Scholar] [CrossRef]
  55. Masip, Y.; Campo, A.; Nuñez, S.M. Experimental analysis of the thermal performance on electronic cooling by a combination of cross-flow and an impinging air jet. Appl. Therm. Eng. 2020, 167, 114779. [Google Scholar] [CrossRef]
  56. Moon, J.W.; Kim, S.Y. Frequency-dependent heat transfer enhancement from rectangular heated block array in a pulsating channel flow. Int. J. Heat Mass Transf. 2005, 48, 4904–4913. [Google Scholar] [CrossRef]
  57. Adhikari, R.C.; Wood, D.H.; Pahlevani, M. An experimental and numerical study of forced convection heat transfer from rectangular fins at low Reynolds numbers. Int. J. Heat Mass Transf. 2020, 163, 120418. [Google Scholar] [CrossRef]
Figure 1. Vibration wave plot of square and sinusoidal wave-shaped vibration.
Figure 1. Vibration wave plot of square and sinusoidal wave-shaped vibration.
Applsci 12 04911 g001
Figure 2. Physical model. (a) Isometric view, (b) side view, (c) top view.
Figure 2. Physical model. (a) Isometric view, (b) side view, (c) top view.
Applsci 12 04911 g002
Figure 3. Computational domain. (a) Heat sink, (b) section view of heat sink and air domain, and (c) air domain.
Figure 3. Computational domain. (a) Heat sink, (b) section view of heat sink and air domain, and (c) air domain.
Applsci 12 04911 g003
Figure 4. Grid convergence.
Figure 4. Grid convergence.
Applsci 12 04911 g004
Figure 5. The variations in the area-averaged Nusselt numbers for different time steps.
Figure 5. The variations in the area-averaged Nusselt numbers for different time steps.
Applsci 12 04911 g005
Figure 6. Time-averaged Nusselt number validation against experimental measurements of Najim, et al. [37].
Figure 6. Time-averaged Nusselt number validation against experimental measurements of Najim, et al. [37].
Applsci 12 04911 g006
Figure 7. Top view of the heat sink.
Figure 7. Top view of the heat sink.
Applsci 12 04911 g007
Figure 8. Velocity vectors in the vicinity of the vibrating heat sink at different phases of a cycle of sinusoidal and square wave-shaped vibration at f = 100 Hz, Amp = 0.003 m, and Re = 1000 (a) without vibration, (b) sinusoidal at 1st quarter of the cycle, (c) sinusoidal at 2nd quarter of the cycle, (d) sinusoidal at 3rd quarter of the cycle, (e) sinusoidal at 4th quarter of the cycle, (f) square at 1st quarter of the cycle, (g) square at 2nd quarter of the cycle, (h) square at 3rd quarter of the cycle, and (i) square at 4th quarter of the cycle.
Figure 8. Velocity vectors in the vicinity of the vibrating heat sink at different phases of a cycle of sinusoidal and square wave-shaped vibration at f = 100 Hz, Amp = 0.003 m, and Re = 1000 (a) without vibration, (b) sinusoidal at 1st quarter of the cycle, (c) sinusoidal at 2nd quarter of the cycle, (d) sinusoidal at 3rd quarter of the cycle, (e) sinusoidal at 4th quarter of the cycle, (f) square at 1st quarter of the cycle, (g) square at 2nd quarter of the cycle, (h) square at 3rd quarter of the cycle, and (i) square at 4th quarter of the cycle.
Applsci 12 04911 g008
Figure 9. Temperature profiles on a plane passing across the tip of the heat sink at Re = 1000: (a) static heat sink; (b) sinusoidal wave shape, f = 100 Hz, Amp = 0.003 m; (c) square wave shape, f = 100 Hz, Amp = 0.003 m; (d) sinusoidal wave shape, f = 100 Hz, and Amp = 0.005 m; (e) square wave shape, f = 100 Hz, Amp = 0.005 m.
Figure 9. Temperature profiles on a plane passing across the tip of the heat sink at Re = 1000: (a) static heat sink; (b) sinusoidal wave shape, f = 100 Hz, Amp = 0.003 m; (c) square wave shape, f = 100 Hz, Amp = 0.003 m; (d) sinusoidal wave shape, f = 100 Hz, and Amp = 0.005 m; (e) square wave shape, f = 100 Hz, Amp = 0.005 m.
Applsci 12 04911 g009
Figure 10. Variation in time-averaged Nusselt number with vibrational parameters and corresponding Reynolds number under nonvibrating conditions: (a) frequency at constant Amp = 0.003 m and Re = 1000; (b) peak-to-peak amplitude at constant f = 100 Hz and Re = 1000.
Figure 10. Variation in time-averaged Nusselt number with vibrational parameters and corresponding Reynolds number under nonvibrating conditions: (a) frequency at constant Amp = 0.003 m and Re = 1000; (b) peak-to-peak amplitude at constant f = 100 Hz and Re = 1000.
Applsci 12 04911 g010
Figure 11. Variation in Nusselt ratio with vibrational parameters at Re = 1000: (a) frequency at constant Amp = 0.003 m; (b) peak-to-peak amplitude at constant f = 100 Hz.
Figure 11. Variation in Nusselt ratio with vibrational parameters at Re = 1000: (a) frequency at constant Amp = 0.003 m; (b) peak-to-peak amplitude at constant f = 100 Hz.
Applsci 12 04911 g011
Figure 12. Variation in percentage reduction in Reynolds number with vibrational parameters: (a) frequency at constant Amp = 0.003 m; (b) peak-to-peak amplitude at constant f = 100 Hz.
Figure 12. Variation in percentage reduction in Reynolds number with vibrational parameters: (a) frequency at constant Amp = 0.003 m; (b) peak-to-peak amplitude at constant f = 100 Hz.
Applsci 12 04911 g012
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rasangika, A.H.D.K.; Nasif, M.S.; Pao, W.; Al-Waked, R. Numerical Investigation of the Effect of Square and Sinusoidal Waves Vibration Parameters on Heat Sink Forced Convective Heat Transfer Enhancement. Appl. Sci. 2022, 12, 4911. https://doi.org/10.3390/app12104911

AMA Style

Rasangika AHDK, Nasif MS, Pao W, Al-Waked R. Numerical Investigation of the Effect of Square and Sinusoidal Waves Vibration Parameters on Heat Sink Forced Convective Heat Transfer Enhancement. Applied Sciences. 2022; 12(10):4911. https://doi.org/10.3390/app12104911

Chicago/Turabian Style

Rasangika, Ambagaha Hewage Dona Kalpani, Mohammad Shakir Nasif, William Pao, and Rafat Al-Waked. 2022. "Numerical Investigation of the Effect of Square and Sinusoidal Waves Vibration Parameters on Heat Sink Forced Convective Heat Transfer Enhancement" Applied Sciences 12, no. 10: 4911. https://doi.org/10.3390/app12104911

APA Style

Rasangika, A. H. D. K., Nasif, M. S., Pao, W., & Al-Waked, R. (2022). Numerical Investigation of the Effect of Square and Sinusoidal Waves Vibration Parameters on Heat Sink Forced Convective Heat Transfer Enhancement. Applied Sciences, 12(10), 4911. https://doi.org/10.3390/app12104911

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