Next Article in Journal
Well Testing Methodology for Multiple Vertical Wells with Well Interference and Radially Composite Structure during Underground Gas Storage
Previous Article in Journal
Fuzzy Logic–Based Decentralized Voltage–Frequency Control and Inertia Control of a VSG-Based Isolated Microgrid System
Previous Article in Special Issue
Review and Analysis of Electro-Magnetohydrodynamic Flow and Heat Transport in Microchannels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design Analysis of Heat Sink Using the Field Synergy Principle and Multitarget Response Surface Methodology

Department of Mechanical Engineering, National Kaohsiung University of Science and Technology, Kaohsiung 807618, Taiwan
*
Author to whom correspondence should be addressed.
Energies 2022, 15(22), 8399; https://doi.org/10.3390/en15228399
Submission received: 19 September 2022 / Revised: 31 October 2022 / Accepted: 7 November 2022 / Published: 10 November 2022
(This article belongs to the Special Issue Methods and Developments for Enhancement of Heat Transfer)

Abstract

:
This study describes a novel heat sink design approach employs the field synergy concept and multitarget response surface methodology (RSM). The multiobjective response surface methodology can be used to determine the simulation equations that will maximize the heat transfer of fins at various fin heights, fin angles, and fin circumferences, when considering the impact of jet flow heat exchange. The goal of the response value was to maintain the minimum possible average field coangle and fin temperature. The results show that the ideal heat sink size would be the following: fins with a height of 50 mm, an angle of 60 degrees, and the number of fins equal to five. We examined the impact of wall speed on the heat transfer caused by the field synergy angle. Our findings suggest that with the synchromesh of the display field, heat-dissipation efficiency rises.

1. Introduction

The Sustainable Development Goals (SDGs) noted that “economic growth”, “social advancement”, and “environmental protection” are universal problems for sustainable development. Taiwan’s industries are under intense pressure to meet worldwide net-zero emission standards and are urgently looking for a low-carbon business model. The heat problem has become an industry problem which academics urgently need to address [1,2,3,4,5,6]. The heat problem is due to the advancement of technology and chip processing to improve the miniaturization and multifunctionality of electronic components. To reduce the temperature of electronic products, heat sinks are a common and low-cost cooling method with effective heat dissipation. The size of the convective heat transfer coefficient is used as the enhancement criterion in classic convection-enhanced heat transfer technologies. The physical mechanism of convective heat transfer and control was reexamined by Bergles [7,8,9]. This new theory is mostly based on the energy equation. Guo and Wang [10,11], who have shown how synergistic convective heat transfer velocity and temperature field can also improve the heat transfer coefficient and raise heat exchanger efficiency, explain the basic concepts of the theory of convective heat transfer, also known as the synergy angle in the field. According to the field synergy theory, the efficiency of improved heat transfer technologies is influenced by the flow velocity vector, the temperature gradient between the angles, the speed field, and the heat flow field. The effect of energy savings is evident, the pump power increase is low, and the heat transfer is improved simultaneously. Numerous theoretical studies pertaining to optimization performance have been conducted as field synergy concept implementations on diverse thermal systems [12,13,14,15,16,17,18,19,20,21].
The link between several explanatory factors and one or more response variables is studied using the response surface methodology (RSM). A solution for one goal may affect other output parameters throughout the process of single objective optimization for individual solutions, leading to unintended results. As a result, the field synergy technique was employed to give effective to design parameters utilizing multiobjective optimization. For quality engineering, the RSM [22,23,24,25] used mathematical programming and the desired function to maximize control factor values. Axial fans were used to install a radiator for forced air cooling. The two major ways to boost forced convection heat dissipation are expanding the heat-dissipation area and raising the heat transfer coefficient, which results in a turbulent flow field. Additionally, increasing the radiator’s surface area to increase cooling effectiveness has frequently been used. In summary, jet flow heat exchange, a kind of heat exchange that makes use of the impact of airflow to enhance convective heat transfer, has the advantage of having a high heat exchange efficiency. In order to attain the anticipated maximum heat transfer efficiency and the lowest possible energy consumption, this study proposes a novel design method that applies axial flow and field synergy theory to the heat sink design without altering the size or shape of the heat sink. This research presents a numerical three-dimensional CFD (computational fluid dynamics) investigation of convective heat transfer for a novel design heat sink.
The main goal of this work is to propose a novel and practical approach for the construction of a heat source with improved heat transfer and to investigate the impact at different fin heights, fin angles, fin circumferences, and efficacy from the perspective of FSP. The heat transfer was improved, according to the results, and this work’s main focus was on numerical three-dimensional simulation of the actual cooling impact. This facilitates rapid heat removal from the system. The ability to heat the fins to alter their alignment, aliquot, height, and convection between them is crucial for this investigation. We provide a succinct overview of design analysis together with numerical analysis and test results for the standard heat sink.

2. Methodology

2.1. Field Synergy Principle (FSP)

Beginning with the energy equation of the three-dimensional convective heat transfer boundary layer, the energy equation is integrated inside the thermal boundary layer thickness. The FSP for incompressible flow and energy conservation is stated simply as follows in [10]:
ρ c p ( u T x + v T y + w T z ) = k 2 T
where cp and k are all constant, the dimensionless governing equations and boundary conditions were created using the following variables.
U ¯ = U U ;   y ¯ = y δ t ; x ¯ = x δ t   ;   z ¯ = z δ t   ; T ¯ = T T T w T ; δ = y = 0 ( 1 | U ¯ |   U ) d y ;   T = ( T w T ) δ t ¯ T ¯ ;   Q = ( k T ) y = 0 ;   and   ¯ = i ^   x ¯ + j ^   y ¯ + k ^   z ¯
The integral of Equation (1) over the thickness of the thermal boundary layer is:
ρ c p 0 δ t ( U · T ) d y = Q
In terms of these following nondimensional variables, the equations of FSP become:
R e = ρ U l μ , P r = v κ v = μ ρ , κ = k ρ C p , N u = h l k
Taking the dimensionless variables and rewriting Equation (3) as:
N u x = Re x Pr y ¯ = 0 1 ( U ¯ ¯ T ¯ ) d y ¯
The static temperature gradient and velocity dot product are given as:
U ¯ ¯ T ¯ = | U ¯ | | ¯ T ¯ | cos θ
The convection heat of the field synergy principle, which is based on the speed field and the heat flow field of synergy, determines whether the fluid flow is enhanced by heat transfer technology. These factors are the flow velocity vector and the temperature gradient between the angles. In the convective heat transfer of the radiator, there is a local angle between each node, and the local angle may signify the interaction between the local velocity vector and temperature gradient. By comparing the two vectors in their entirety, we can determine the average angle of the entire field. The angle between the velocity and temperature fields and the concept of field synergy both have an impact. From the preceding explanation, it is clear that when the angle is smaller and the heat transfer intensity is likewise good, the temperature difference between the flow field and the flow velocity is at its greatest. The formula for the synergistic angle can be found by using Equation (6), therefore it is represented as Equation (7) [26]:
θ = cos 1 ( U ¯ ¯ T ¯ | U ¯ | | ¯ T ¯ | )

2.2. Response Surface Methodology

Using the response surface approach, the relationship between a number of explanatory variables and one or more response variables is examined. The method was first introduced by George E. P. Box and K. B. Wilson in 1951. Getting out well prepared through second-order well-prepared examinations is the main tenet of RSM. The mathematical formulae for this inquiry were derived using second-order regression equations.
Y = β 0 + β 1 X 1 + β 2 X 2 + β 3 X 3 + β 11 X 1 2 + β 22 X 2 2 + β 33 X 3 2 + β 12 X 1 X 2 +   β 13 X 1 X 3 + β 23 X 2 X 3
where Y denotes the reaction variables of the heat sink, such as the field synergy angle and average fin temperature. β 0 is the constant term of this equation, β 1 ,   β 2 , β 3 , and β ij are the coefficients of the linear and the quadratic terms, respectively. Xi is the input factors namely fin height, heat sink angle, and fin aliquots. A solution for a single goal in single objective optimization for individual answers may impact other output parameters in the process, resulting in undesirable outcomes. As a result, multiobjective optimization was used to achieve effective design parameters using the Minitab software. Applying the desirability function, we converted the output value (yi) to independence desirability (di). We then calculated composite desirability (D) as follows:
D = d 1 × d 2 × d 3 × × d x x
where the x represents the number of output values (d). Composite desirability (D) is between 0 and 1 (the closer to 1, the better).

3. Design Analysis

3.1. Heat Sink Model

This model takes advantage of the geometric design features of COMSOL multiphysics software to create the model in Figure 1 and Figure 2, which has an aluminum heat sink to cool the square thermal section. The heat sink specifications were: length 68 mm, width 83 mm, thickness 1 mm, span 2 mm, forced fan ψ45 mm, and rpm 3000–4000.
A similar device was used to gauge the fins’ thermal resistance. The acrylic materials chamber circulated air, and the base of the heat fin was a square thermal bar with a fixed power source. The other external faces were covered with thermal insulation. The static pressure curve parameters were set for the fan, a fan was positioned near the backplane’s exit, and a fan was installed above the fins. The air entrance, by natural convection, was through the back side panel. A square insulating body was placed in the heat source to cover the area around the heating rod, minimizing heat loss.
In order to reduce the complexity of the problem and the amount of time required for computation, the finite element approach of the software itself was used to calculate heat transfer in this study. The governing equation was described, the heat transfer phenomena in the solution region was calculated, and the basic physical model assumptions were made. The study makes the following assumptions:
  • There will be a steady-state heat transfer.
  • The fluid in the chamber is three-dimensional forced convection.
  • The working fluid has an incompressible density.
  • The model space’s walls are all nonslip surfaces.
  • There is little radiation impact.
In this numerical analysis, these solvers are step-steady-state and steady-state. The convergence criterion for every iteration of the procedure was specified so that the residual value is below each range. Relative convergence was the criterion in this study. The convergence condition was satisfied when the value falls below 10−4.

3.2. Boundary Conditions

The turbulent flow field is typically used by the numerical simulation flow field. The analysis of the fin heat-dissipation system most frequently uses the k-mode and the k-mode, according to literature. There is not currently a single turbulent flow mode that can handle all issues. Through numerical simulations and experiments, the two groups of models were compared, and in the same construction model, the model that best matched the outcomes of the experiments was chosen. Following an experimental comparison, it was decided to conduct a simulation analysis of the k- turbulent flow model because the results were found to be more accurate. The boundary condition is the setting that matters the most in the numerical simulation. A more complete set of conditions will produce results that are close to the experimental value, but too many conditions will cause the calculation to become overly complex and put the cart before the horse. For this reason, it is crucial to be able to set the proper boundary conditions. Table 1 lists the boundary and physical conditions discussed in this study.
Physical control and user control were separated in COMSOL software’s mesh design. The latter was employed in this work, and finer mesh configurations were created for various locations, as indicated in Table 2 and Figure 3.

3.3. Experiment Measuring

The data from the thermal imager and anemometer were compared in order to validate the accuracy of the simulation analysis. The basic specifications of the thermal imager (Thermo GEAR G120) were the following:
(A)
Measuring range: −40 °C to 500 °C
(B)
Resolution: 0.06 °C
(C)
Accuracy: ±2 °C or ±2% of reading, whichever was greater
Additionally, the basic specifications of the anemometer (UT-362-USB) were:
(A)
Measuring range: 2~30 m/s
(B)
Resolution: 0.1 m/s
(C)
Accuracy: 2–10 m/s (±3% ± 0.5); 10–30 m/s (±3% ± 0.8)
Figure 4a shows the actual temperature field of the heat sink using an infrared thermal image apparatus and measures the airflow velocity with an anemometer shown in Figure 4b.

3.4. Results and Discussion

In order to create models and ideal designs, the response surface method (RSM), an expanded experiment method, combines statistics and optimization techniques. The multiobjective response surface method is the best improvement model. After establishing the design parameters for the heat sink, we used the response surface method’s central composite design to create a factor distribution. The finite element method was then used to analyze the location of the factors. Twenty experiments for the heat sink design were used to provide the results of the simulation study summarized in Table 3 and Figure 5a–c. Based on the above data, we calculated the actual synergistic angle using Equation (9).
Finding the impact of a particular factor’s magnitude requires regression analysis. As indicated in Table 4 for the field synergy angle and Table 5 for the average fin temperature, if (p < 0.05), the factor could be eliminated and the RSM could be reconstructed.
After significance analysis, the expression for the field synergy angle of response surface equation Y1 was as follows:
Y1 = 84.4161 + 0.555736A − 0.113884B − 2.21681C − 0.004475A2 + 0.00100056B2 + 0.2955C2 − 0.00171042A ×
B − 0.0212188A × C + 0.00760208B × C
Y2 represents the expression for Fin average temperature:
Y 2   = 4.78346 0.0171883 A 8.31712   ×   10 4 B 0.0870512 C   +   4.80948   ×   10 5 A 2     1.18980   ×   10 5 B 2   +   0.00686988 C 2   +   4.51505   10 5 A   ×   B   +   0.000818858 A   ×   C     1.79596
The distribution of the field synergy angle from the bottom to the top of the fin is shown by the XY section in Figure 6a–c. The field synergy angle is approximately 20–50 degrees when viewed from above since the first encounter has a higher flow velocity, but the field synergy angle rises as one approaches the bottom of the fin because the wind speed naturally decreases and the heat transfer enhancement becomes relatively minimal. The opposite has changed with regard to the field synergy angle in the second half of the fin because of the temperature field. The field synergy angle of the fins in the second half is higher because the gradient change and the velocity field change are not coordinated, causing the velocity field change to be unable to keep up with the temperature difference between the fins.
Through a variety of response surface approach experimental designs, the best model was discovered. Table 6 displays the field coangle, velocity, and temperature at the greatest temperature for a single chip. The degree of heat transmission may significantly increase when the fluid enters a jet, and the common interpretation is that as the temperature boundary layer thins, the wall temperature gradient increases. However now, from the standpoint of field synergy, we could look into how wall velocity affects heat transfer. The efficiency of heat dissipation increased with the synchromesh of the display field.

4. Conclusions

The field synergy idea and multitarget response surface methodology (RSM) were used in a revolutionary heat sink design strategy to maximize the heat transfer. There are some conclusions that can be made:
  • We investigated the effects of wall velocity on heat transport due to the field synergy angle. The effectiveness of heat dissipation increased with the synchromesh of the display field.
  • The gradient and velocity field changes were not coordinated, which prevented the velocity field change from keeping up with the temperature difference between the fins. This resulted in a higher field synergy angle of the fins in the second half.
  • When the synergy angle was employed to show the fins’ effects in heat transfer, the research findings supported the design combination. It is important to show the effectiveness of the heat-dissipation effect.

Author Contributions

Conceptualization, M.-C.L.; methodology, M.-C.L.; software, R.-F.L.; formal analysis, M.-C.L.; investigation, M.-C.L.; data curation, R.-F.L.; writing—original draft preparation, M.-C.L.; supervision, M.-C.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Kaohsiung University of Science and Technology under Grant No. 111E9010BA02 and the APC was funded by National Kaohsiung University of Science and Technology.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

cpHeat capacity, J/(kg∙K)
hHeat transfer coefficient, W/(m2∙K)
kThermal conductivity, W/(m∙K)
lLength, m
NuNusselt number, dimensionless
PrPrandtl number, dimensionless
QConduction heat source, W/m2
ReReynold’s number, dimensionless
tTime, s
TTemperature, K
UFluid velocity field, m/s.
U ¯ Mean velocity vector, m/s.
uVelocity in the x direction, m/s.
vVelocity in the y direction, m/s.
wVelocity in the z direction, m/s.
xSpatial coordinate
XiInput factor
ySpatial coordinate
yiOutput value
YReaction variables
zSpatial coordinate
Greek symbols
δ Velocity boundary layer thickness, m
δ t Temperature boundary layer thickness, m
ρ Fluid density, kg/m3
μDynamic viscosity, Pa/s
νKinematic viscosity, m2/s
θ Field synergistic angle
▽TTemperature gradient vector, K
β ijCoefficient of the linear and quadratic terms

References

  1. Wang, P.; Cui, X.; Weng, J.; Cai, Z.; Cai, R. Experimental investigation of the heat transfer performance of an oscillating heat pipe with LiCl salt solution. Int. J. Heat Mass Transf. 2020, 158, 120033. [Google Scholar] [CrossRef]
  2. Ji, X.; Yang, X.; Zhang, Y.; Zhang, Y.; Wei, J. Experimental study of ultralow flow resistance fractal microchannel heat sinks for electronics cooling. Int. J. Therm. Sci. 2022, 179, 107723. [Google Scholar] [CrossRef]
  3. See, Y.; Ho, J.; Leong, K.; Wong, T. Experimental investigation of a topology-optimized phase change heat sink optimized for natural convection. Appl. Energy 2022, 314, 118984. [Google Scholar] [CrossRef]
  4. Shaeri, M.R.; Sarabi, S.; Randriambololona, A.M.; Shadlo, A. Machine learning-based optimization of air-cooled heat sinks. Therm. Sci. Eng. Prog. 2022, 34, 101398. [Google Scholar] [CrossRef]
  5. Muneeshwaran, M.; Lee, Y.-J.; Wang, C.-C. Performance improvement of heat sink with vapor chamber base and heat pipe. Appl. Therm. Eng. 2022, 215, 118932. [Google Scholar] [CrossRef]
  6. Leonardo, M.; Tapas, K.M.; Reddy, K.S. General correlations among geometry, orientation and thermal performance of natural convective micro-finned heat sinks. Int. J. Heat Mass Trans. Sci. Direct. 2015, 91, 711–724. [Google Scholar]
  7. Bergles, A.E. Heat Transfer Enhancement—The Maturing of Second-Generation Heat Transfer Technology. Heat Transf. Eng. 1997, 18, 47–55. [Google Scholar] [CrossRef]
  8. Bergles, A.E. Techniques to Enhance Heat Transfer. In Handbook of Heat Transfer, 3rd ed.; Rohsenow, W.M., Hartentt, J.P., Cho, Y.I., Eds.; McGraw-Hill: New York, NY, USA, 1998. [Google Scholar]
  9. Bergles, A.E. Enhanced Heat Transfer: Endless Frontier, or Maarmature Routine. Enhanced Heat Transfer. 1999, 6, 79–88. [Google Scholar] [CrossRef]
  10. Guo, Z.Y.; Li, D.Y.; Wang, B.X. A novel concept for convective heat transfer enhancement. Int. J. Heat Mass Transf. 1998, 41, 2221–2225. [Google Scholar] [CrossRef]
  11. Wang, S.; Li, Z.X.; Guo, Z.Y. Novel concept and device of heat transfer augmentation. In Proceedings of the Eleventh International Conference on Heat Transfer, Kyongju, Korea, 1 November 1998; pp. 405–408. [Google Scholar]
  12. Li, X.; He, Y.-L.; Tao, W.-Q. Analysis and extension of field synergy principle (FSP) for compressible boundary-layer heat transfer. Int. J. Heat Mass Transf. 2015, 84, 1061–1069. [Google Scholar] [CrossRef]
  13. Zhao, J.; Huang, S.; Gong, L.; Huang, Z. Numerical study and optimizing on micro square pin-fin heat sink for electronic cooling. Appl. Therm. Eng. 2016, 93, 1347–1359. [Google Scholar] [CrossRef]
  14. Zhanga, Y.; Liu, X. Application of Field Synergy Principle for Fin Reshaping of a Natural Convection Radiator. In Proceedings of the 9th International Symposium on Heating, Ventilation and Air Conditioning (ISHVAC) and the 3rd International Conference on Building Energy and Environment (COBEE), Tianjin, China, 12–15 July 2015; Volume 121, pp. 1726–1733. [Google Scholar]
  15. Yao, Y.; Weiwei, W.; Yang, K. Mechanism study on the enhancement of silica gel regeneration by power ultrasound with field synergy principle and mass diffusion theory. J. Heat Mass Transf. 2015, 90, 769–780. [Google Scholar] [CrossRef]
  16. Hamid, O.A.M.; Zhang, B.; Yang, L. Application of field synergy principle for optimization fluid flow and convective heat transfer in a tube bundle of a pre-heater. Energy 2014, 76, 241–253. [Google Scholar] [CrossRef]
  17. Ou, J.-J.; Li, L.-F.; Cui, T.; Chen, Z.-M. Application of field synergy principle to analysis of flow field in underhood of LPG bus. Comput. Fluids 2014, 103, 186–192. [Google Scholar] [CrossRef]
  18. Kundu, B.; Yook, S.J. An accurate approach for thermal analysis of porous longitudinal, spine and radial fins with all nonlinearity effects–analytical and unified assessment. Appl. Math. Comput. 2021, 402, 126124. [Google Scholar] [CrossRef]
  19. Wankhade, P.A.; Kundu, B.; Das, R. Establishment of non-Fourier heat conduction model for an accurate transient thermal response in wet fins. Int. J. Heat Mass Transf. 2018, 126, 911–923. [Google Scholar] [CrossRef]
  20. Das, R.; Kundu, B. Forward and Inverse Analyses of Two-Dimensional Eccentric Annular Fins for Space-Restriction Circumstances. J. Thermophys. Heat Transf. 2021, 35, 80–91. [Google Scholar] [CrossRef]
  21. Kundu, B. Exact Method for Annular Disc Fins with Heat Generation and Nonlinear Heating. J. Thermophys. Heat Transf. 2017, 31, 337–345. [Google Scholar] [CrossRef]
  22. Liu, C.; Yang, H. Multi-objective optimization of a concrete thermal energy storage system based on response surface methodology. Appl. Therm. Eng. 2021, 202, 117847. [Google Scholar] [CrossRef]
  23. Mohammed, B.S.; Khed, V.C.; Nuruddin, M.F. Rubbercrete mixture optimization using response surface methodology. J. Clean. Prod. 2018, 171, 1605–1621. [Google Scholar] [CrossRef]
  24. Mohammed, B.S.; Fang, O.C.; Anwar, H.; Lachemi, M. Mix proportioning of concrete containing paper mill residuals using response surface methodology. Construct. Build. Mater. 2012, 35, 63–68. [Google Scholar] [CrossRef]
  25. Abdulkadir, I.; Mohammed, B.S.; Liew, M.; Wahab, M. Modelling and multi-objective optimization of the fresh and mechanical properties of self-compacting high volume fly ash ECC (HVFA-ECC) using response surface methodology (RSM). Case Stud. Constr. Mater. 2021, 14, e00525. [Google Scholar] [CrossRef]
  26. Lin, M.-C.; Lin, R.-F. Innovative design for an active uniform heat dissipation system. Results Phys. 2020, 19, 103355. [Google Scholar] [CrossRef]
Figure 1. Thermal resistance measurement apparatus for CPU coolers.
Figure 1. Thermal resistance measurement apparatus for CPU coolers.
Energies 15 08399 g001
Figure 2. Model for measuring the thermal resistance of heat dissipation.
Figure 2. Model for measuring the thermal resistance of heat dissipation.
Energies 15 08399 g002
Figure 3. Mesh of numerical model.
Figure 3. Mesh of numerical model.
Energies 15 08399 g003
Figure 4. (a) actual temperature field of the heat sink using an infrared thermal image apparatus; (b) actual airflow velocity using an anemometer.
Figure 4. (a) actual temperature field of the heat sink using an infrared thermal image apparatus; (b) actual airflow velocity using an anemometer.
Energies 15 08399 g004
Figure 5. (a) Fin height, (b) heat sink angle, (c) fin aliquots.
Figure 5. (a) Fin height, (b) heat sink angle, (c) fin aliquots.
Energies 15 08399 g005
Figure 6. (a) distribution of the fin’s top field synergy angles; (b) distribution of the fin’s middle field synergy angles; and (c) distribution of the fin’s bottom field synergy angles.
Figure 6. (a) distribution of the fin’s top field synergy angles; (b) distribution of the fin’s middle field synergy angles; and (c) distribution of the fin’s bottom field synergy angles.
Energies 15 08399 g006
Table 1. Boundary circumstances Setting.
Table 1. Boundary circumstances Setting.
1. Material
NameMaterialThermal conductivity [W/(m·K)]
Heat SinkAluminum 6063-T83201
ChamberAcrylic plastic0.18
Heat sourceAl2O3-Aluminum oxide30
Thermal insulation coating materialOwens–Corning fiberglass0.018
Heat sticksRed copper385
FluidAir0.026
2. Turbulence mode: k-ω
Condition attributeParameterPosition
k-ωIncompressible flowInside the entire chamber
WallWall function (no slippage)Except for entrance, exit, fan
Entrance1 m/s Wall rear
Export1 atmWall
Fan 1PQ curveAbove the fin
Fan 2PQ curveChamber fan
3. Heat transfer mode
Condition attributeParameterPosition
Heat transferAmbient temperature: 25 °CThe whole chamber
Solid heat transferDepending on the material typeFins, fans, thermal insulation, shells
Fluid heat transferAirInside the chamber
Initial value25 °COverall
Heat source80 WHeat source
Thin layerHeat paste:
Layer thickness: 50 μm
With ST-350
Between the fins and the
thermal bars
TemperatureEntry temperature: 25 °CChamber rear
OutflowExportFan and drain on both sides
Table 2. Mesh construction parameters.
Table 2. Mesh construction parameters.
Node number20743
Tetrahedron elements number 122097
Triangle elements number34126
Minimum element quality0.005516
Average element quality0.5344
Element volume ratio9.048 × 10−8
Mesh volume1.944 × 108 mm 3
Table 3. Response surface experiment design.
Table 3. Response surface experiment design.
ItemA
Fin Height
(mm)
B
Heat Sink Angle (Degree)
C
Fin
Aliquots
Field Synergy Angle
(Y1)
Fin Average
Temperature (°C)
(Y2)
1300390.60567.715
24030390.08959.309
34030587.95359.793
43030587.95359.793
55030587.95359.793
6300791.30368.211
7400592.77958.428
85060387.00754.403
9500393.40951.529
104030587.95359.793
11500792.56857.225
124030789.07663.077
134030587.95359.793
143060788.77867.924
155060787.83258.260
164030587.95359.793
174030587.95359.793
184030587.95359.793
194060585.82359.320
203060386.09765.586
Table 4. Experimental regression analysis for the field synergy angle.
Table 4. Experimental regression analysis for the field synergy angle.
Response Surface Regression: Y1 versus A, B, C
Estimated Regression Coefficients for Y1
TermCoefficientp value
A88.13200
B0.40330.095
C−2.51270
A × A0.23500.308
B × B−0.44750.309
C × C0.90050.056
A × B1.18200.018
A × C−0.51310.062
B × C−0.42440.113
R-Sq(adj) = 89.69%
Table 5. Experimental regression analysis for the average fin temperature.
Table 5. Experimental regression analysis for the average fin temperature.
Response Surface Regression: Y2 versus A, B, C
Estimated Regression Coefficients for Y2
TermCoefficientp value
A4.088960
B−0.078920
C0.005120.696
A × A0.027730.054
B × B0.004810.847
C × C−0.010710.668
A × B0.027480.284
A × C0.013550.364
B × C0.016380.277
R-Sq(adj) = 66.99%
Table 6. Best optimization combinations.
Table 6. Best optimization combinations.
ItemFin HeightHeat Sink AngleFin
Aliquots
Fin’s
Bottom
Field
Synergy
Fin’s Middle
Field Synergy
Fin’s
Top
Field
Synergy
Best
combination
50 mm60°553.47673.58142.039
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lin, M.-C.; Lin, R.-F. Design Analysis of Heat Sink Using the Field Synergy Principle and Multitarget Response Surface Methodology. Energies 2022, 15, 8399. https://doi.org/10.3390/en15228399

AMA Style

Lin M-C, Lin R-F. Design Analysis of Heat Sink Using the Field Synergy Principle and Multitarget Response Surface Methodology. Energies. 2022; 15(22):8399. https://doi.org/10.3390/en15228399

Chicago/Turabian Style

Lin, Ming-Che, and Ruei-Fong Lin. 2022. "Design Analysis of Heat Sink Using the Field Synergy Principle and Multitarget Response Surface Methodology" Energies 15, no. 22: 8399. https://doi.org/10.3390/en15228399

APA Style

Lin, M. -C., & Lin, R. -F. (2022). Design Analysis of Heat Sink Using the Field Synergy Principle and Multitarget Response Surface Methodology. Energies, 15(22), 8399. https://doi.org/10.3390/en15228399

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