Next Article in Journal
Water Modification by Cold Plasma Jet with Respect to Physical and Chemical Properties
Previous Article in Journal
A Spectrophotometry Evaluation of Clear Aligners Transparency: Comparison of 3D-Printers and Thermoforming Disks in Different Combinations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Conjugate Heat Transfer Analysis and Heat Dissipation Design of Nucleic Acid Detector Instrument

1
School of Mechanical and Automotive Engineering, Xiamen University of Technology, Xiamen 361024, China
2
Xiamen Key Laboratory of Intelligent Manufacturing Equipment, Xiamen 361005, China
3
School of Public Health, Xiamen University, Xiamen 361005, China
4
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, Xiamen 361005, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(23), 11966; https://doi.org/10.3390/app122311966
Submission received: 25 October 2022 / Revised: 15 November 2022 / Accepted: 21 November 2022 / Published: 23 November 2022
(This article belongs to the Section Applied Thermal Engineering)

Abstract

:
Temperature affects both the stability of nucleic acid detectors and efficiency of DNA amplification. In this study, temperature and flow inside a nucleic acid detector were simulated and the results were used to design vents for the instrument casing. A test platform was constructed to collect experimental temperature data that were used for simulation validation. The experimental and simulation results showed that the temperature error was less than ±3 K. A total of two heat-dissipation schemes were designed based on the simulation and a new instrument casing was fabricated based on the scheme with the best results. Nucleic acid amplification was performed continuously for 120 min using a prototype with the new casing. The temperatures of the monitoring points were stable and the maximum temperature measured only 307.76 K (34.61 °C). Therefore, waste heat was effectively eliminated, which ensured safety of the electronic components and stability of the nucleic acid detection process.

1. Introduction

Quantitative real-time polymerase chain reaction (RT-PCR) is a vital technique in the diagnostic industry as it facilitates rapid detection, the importance of which was emphasized by the COVID-19 outbreak in December 2019. In particular, the use of PCR technology in combination with nucleic acid detectors has attracted considerable attention from researchers [1,2,3]. These instruments provide the specific temperatures required for each stage of the PCR process: (1) DNA separation at temperatures greater than 363.15 K, (2) primer annealing at 323.15–348.15 K, and (3) optimal extension at 345.15–351.15 K [4,5]. The precise temperature depends on the reagent system. The stages are repeated until the DNA is sufficiently amplified for fluorescence detection.
However, the internal temperature of the instrument increases as the operation time increases, which increases the failure rate of the components; therefore, heat dissipation is an important consideration in the design of these devices [6,7,8]. Furthermore, efficient heat dissipation increases the cooling rate during temperature cycling, which reduces the time required for DNA amplification and improves the sample detection rate during extensive batch testing. Therefore, it is necessary to research heat dissipation to improve the reliability and stability of PCR devices.
Many studies on PCR instruments have focused on nucleic acid amplification and fluorescence detection. Research on the amplification stage has typically focused on temperature homogeneity [9], temperature control [10], and the rate of temperature changes [11]. By contrast, research on the fluorescence detection stage has typically focused on the fluorescent light pathway, structure of the detection system, and signal-processing algorithms [12,13]. Furthermore, researchers have investigated the heat dissipation and thermal design of PCR instruments. For example, Aziz et al. [14] simulated the 3D heat transfer of a continuous-flow PCR model and used ANSYS CFX to investigate the effect of the fluids on the temperature distribution in the canals. Naghdloo et al. [15] used a numerical simulation of a fully automated PCR system to improve the structure of the operation unit and reduce the operation time. Subsequently, they optimized the time consumption by applying a mixed convective heat transfer mechanism. Parng et al. [16] developed an innovative design in which multilayer microchannels were built in a single tube; thus, droplet generation and PCR could be accomplished in the same tube. Their experiments showed that this design had high heating and cooling rates. Zhou et al. [17] analyzed the pedestal heat transfer of a PCR device and proposed a mathematical model to describe the pedestal temperature field, which was simplified into a simple 2D heat transfer model. Thus, they analyzed the effects of various parameters on the temperature field of the pedestal during uneven heating. Finally, Raza et al. [18] designed a well-shaped sample block with a volume of 50 µL and demonstrated that the fluid sample exhibited a higher temperature response rate and greater uniformity. This pedestal can be manufactured without any special process, paving the way for low-cost PCR.
However, there has been little research on the design of nucleic acid detectors concerning heat dissipation. Therefore, this study aimed to design an instrument ventilation scheme that could provide stable and safe ambient temperatures. This was achieved by using the finite-volume method to simulate the temperature distribution and flow, validating the simulation using experimental data, and using the simulation to design a new instrument casing. Finally, the experimental results prove that the designed heat-dissipation scheme demonstrates high heat-dissipation capability.

2. Geometry and Experimental Method

2.1. Geometry

The structure of the portable automatic nucleic acid detector used in this study is shown in Figure 1. The overall length L, width W, and height H were 0.360, 0.326, and 0.44 m, respectively, and the casing was 0.005 m thick. Numerous techniques can be used to design the vents on the device; however, a closed model was selected for this study, and the vents were designed according to the simulation results.

2.2. Experimental Method

As shown in Figure 1c,d, the DNA amplification module used a thermoelectric cooler (TEC) to increase and decrease the temperature of the biochip; therefore, the TEC was considered to be a heat source. A test platform, shown in Figure 2, was constructed to identify other heat sources. An infrared thermal imager was used to determine the positions of the heat sources (FLUKE Ti401 PRO, resolution: 640 × 480, range: 253.15–923.15 K, accuracy: ±2 K) and a thermocouple was used to measure their temperatures (thermocouple: KAIPUSEN K, range: 253.15–473.15 K, accuracy: ±0.1 K). Imaging was conducted at room temperature (approximately 299.15 K) and the results are shown in Figure 3. The other heat sources (called primary heat sources) in the instrument were located on the printed circuit boards (PCBs). Thus, there were three components acting as heat sources in the instrument: TEC and two PCBs.
The temperature measured at the surface of a heat source was assumed to be the temperature of the heat source. The temperatures on each side of the TEC and the outer surface temperatures of the primary heat sources on the PCBs were recorded with a data acquisition interval of 2 s. The temperatures of the heat sources throughout the DNA amplification process are shown in Figure 4. As shown in Figure 4a, a maximum temperature of 368.15 K was recorded on TEC side1 during the amplification cycling stage. Figure 4b shows the temperatures of the primary heat sources in the PCBs, which increased initially before stabilizing.

3. Simulation Process

The heat-source temperature data were used as source terms for the simulation. To improve the efficiency of the simulation, the duration of the simulation was half that of the entire DNA amplification process, that is, 600 s. The simulation and experimental data were compared at 30 s intervals.
All heat sources transfer heat to solids they are in contact with via conduction and with the air via convection [19,20]. However, Fedorov and Viskanta [21] suggested that radiation may be ignored when the typical temperature is less than 373.15 K. The maximum temperature in this study was 368.15 K (TEC side1); therefore, conduction and convection had a dominant effect, and radiation was ignored.

3.1. Material Parameters

The materials used in the modules and electronic components of the instrument included aluminum, steel, acrylic, and copper. PCBs generally have a “mezzanine” design consisting of copper and FR4. Considering the equivalent modeling method for PCBs reported by Mark [22], the plane thermal conductivity k p and normal thermal conductivity k n were used to represent the thermal conductivity of the PCBs. These terms are given by:
k p = δ c u j = 1 N k j + ( N 1 ) δ F k F δ t
k n = δ t δ c u j = 1 N 1 k j + ( N 1 ) δ F k F
where k j , k F , and N are the thermal conductivities of the copper-clad layer, FR4, and the number of layers of the copper cladding, respectively; δ t is the total thickness of the PCB; δ c u and δ F are the thicknesses of the copper classing and the unclad layer, respectively; PCBs typically consist of FR4 and copper; therefore, k j is between k F and k c u . That is:
k j A j A t k c u
where A j is the area of the j-th layer of the copper cladding, A t is the area of the PCB, k c u is the thermal conductivity of the copper. In an equivalent model, δ c u and δ F are considered to be equal; that is:
δ t = N δ c u + ( N 1 ) δ F
The parameters of the PCBs are given in Table 1 and Table 2. Using these parameters, the thermal conductivities of PCB1 and PCB2 were determined to be k p , p c b 1 = 31.83 W·m−1·K−1; k n , p c b 1 = 0.38 W·m−1·K−1; k p , p c b 2 = 30.83 W·m−1·K−1; k n , p c b 2 = 0.38 W·m−1·K−1. The physical parameters of other materials are listed in Table 3. The contact pin between the component and base plate was not considered in the modeling. The component materials used corresponded to the packaging materials listed in Table 3.

3.2. Computational Domain and Boundary Conditions

ANSYS determines the fluid temperature and velocity fields by solving the governing equations, which include the mass, momentum, and energy conservation equations [23,24]. For the simulation, the corresponding boundary conditions must be defined to resolve the corresponding equations [25,26]. The fluid and solid domains are illustrated in Figure 5. To reduce the complexity of the analysis, it was assumed that:
  • The flow remains incompressible, ideal, and single phase.
  • The initial (ambient) temperature remained constant.
  • The heat lost from the enclosure to the surroundings was negligible.
  • The connections between parts were considered seamless.
Conjugate heat transfer control equations based on the above assumptions are solved by using the Reynold-average [27,28,29].
Continuity equation:
ρ f t + d i v ( ρ f u ) = 0
where ρ f is the density of the fluid, u is the velocity vector. When the fluid is incompressible, ρ f t = 0 .
Momentum in x-direction:
( ρ f u ) t + d i v ( ρ f u u ) = d i v ( μ f g r a d u ) p x + S u + [ ( ρ f u 2 ¯ ) x ( ρ f u v ¯ ) y ( ρ f u w ¯ ) z ]
Similarly, the momentum equation in y- and z-directions can be written. Where u , v and w are velocity components in the x-, y- and z-directions, respectively. μ f is the dynamic viscosity of the fluid, p is the fluid pressure, S u is source term, u i u j ¯ is the specific Reynolds stress term, and i and j are Cartesian coordinates.
Energy equation for fluid domain is given by:
( ρ f T ) t + d i v ( ρ f u T ) = d i v ( k f C p , f g r a d T ) + S f + [ ( ρ f u T ¯ ) x ( ρ f v T ¯ ) y ( ρ f w T ¯ ) z ]
Energy equation for solid domain is given by:
( ρ s C p , s T ) t = d i v ( k s T ) + Q
where k f and k s are the thermal conductivities of the fluid and solid, respectively, T is temperature, ρ s is the density of the solid, C p , f and C p , s are the specific heats of the fluid and solid, respectively, S f and Q are source terms.
The k−ω shear stress transport (SST) turbulence model was adopted. The k−ω SST model includes all the refinements of the k−ω BSL model, and accounts for the transportation of the turbulence shear stress in the definition of turbulent viscosity. The transport equations for the k−ω BSL model are [30,31,32]:
( ρ k ) t + ( ρ k u i ) x i = x j ( Γ k k x j ) + G k Y k + S k + G b
( ρ ω ) t + ( ρ ω u i ) x i = x j ( Γ ω ω x j ) + G ω Y ω + D ω + S ω + G ω b
where G k represents the production of turbulence kinetic energy, G ω represents the generation of w , Γ k and Γ ω represent the effective diffusivity of k and ω , respectively. Y k and Y ω represent the dissipation of k and ω due to turbulence, D ω represents the cross-diffusion term, S k and S ω are user-defined source terms, G b and G ω b account for buoyancy terms. A detailed description of k−ω equations can be found in [30,31,32].
In the simulation, a fan provided the sole source of the flow. The fan was considered to be infinitely thin and the discontinuous increase in pressure across it was specified as a function of the velocity of the fluid through the fan [33,34]. As shown in Figure 6, a polynomial was used to define the pressure gradient caused by the fan. That is:
Δ P = 71.053 8.6303 v 0.4238 v 2
where Δ P is the pressure jump, and v is the magnitude of the local fluid velocity normal to the fan.
Figure 5 shows a model of the PCB plate and the corresponding heat sources, heat sources were used as the temperature source term, the values are shown in Figure 4b. The experimental data obtained from TEC side1 and side2 were included the final temperature under the effect of the entire TEC, side1 and side2 were used as the temperature wall, The values are shown in Figure 4a.
At the fluid–solid interface [35],
u = 0 ; T f = T s ; k f T f n = k s T s n
At the outer wall of the housing,
T s n = 0
where T f and T s are fluid temperature and solid temperature, respectively.
In this study, the solar-type pressure-based method was used. The simulation settings are listed in Table 4. A total of twenty iterations were performed for each time step. The simulations were considered to converge when the energy residual curve was less than 1 × 10−6, and other residual curves were less than 1 × 10−3.

3.3. Grid Independence

Mesh independence is an important aspect of the computational fluid dynamics simulation. Although the model used in this study was simplified, it still exhibited a high degree of complexity. The grid independence of the model was verified to ensure the accuracy of the simulation [36]. The average fluid temperature T a v e in the last moment of the simulation was used as the verification index for grid independence. The relative change rate η was defined as:
η i = | T a v e , i T a v e , i 1 T a v e , i 1 | × 100 %
where T a v e , i is the average value of the simulation of T a v e at the current grid scale and T a v e , i 1 is the average value at the previous grid scale. The calculated value of T a v e is shown in Figure 7. When the mesh size was less than 0.008 m, the T a v e was approximately unchanged and η was less than 0.5%. This agrees with the grid independence [37]. The mesh size of the last refinement was used to determine the simulation accuracy. As shown in Figure 8, the maximum element size was 0.005 m and the number of grids was 1,271,934.

4. Simulation Results and Discussion

In this section, the temperatures of all the monitoring points (K1 and K2, B1–B4, p1–p9) were measured using the platform shown in Figure 2 (using thermocouples). The measurements were repeated five times and the system was allowed to cool naturally to room temperature (approximately 299.15 K) after each measurement.

4.1. Amplification Module

The simulation results for the amplification module are shown in Figure 9. Figure 9a,c show that the temperature decreased as the distance outward from the heat source increased. In addition, there was an apparent temperature gradient across the fins. The frame and fins conducted and dissipated heat. Figure 9b shows that the temperature was higher in the central region, owing to the position of the amplification module.
The temperatures at points K1 and K2 were measured to verify the accuracy of the amplification model. A comparison of the experimental and simulation data is shown in Figure 10. During the reverse transcription phase, the temperature of TEC side2 was more consistent. However, as the instrument entered the pre-denaturation and cycle phases, the overall simulated temperatures of TEC side2 were slightly higher than the experimental values. In the model, there are gaps between the parts and the air in these gaps increases the thermal resistance. The simplified simulation model ignores the thermal resistance of the connections between the parts; hence, the thermal conductivity of the simulation is higher than that of the actual model. Therefore, when the temperature changes, the simulation value is always greater than the experimental value. Overall, the trends of the simulated and experimental results were consistent.

4.2. PCB

Figure 11 shows the simulated and actual temperature distributions for the two PCBs, which were similar in both cases. There were a large number of relatively centralized heat sources on PCB1 which were clearly shown in the simulation. However, there were some differences between the simulated and experimental results for PCB2. In particular, a minor heat source on the right side of PCB2 was captured by the thermal imager but not considered in the simulation.
To verify the accuracy of the PCB model simulation, points B1–B4 were studied without affecting the operation of the PCBs. A comparison of the experimental and simulation results is shown in Figure 12. During the recording process, the initial value was always higher than room temperature (approximately 299.15 K). This is because there were many components and copper wires in the PCBs which caused Joule heating during operation and affected the PCBs. In addition, the secondary heat source was ignored in the model. Therefore, the simulated temperature at each point was lower than the experimental value. Moreover, it was difficult to maintain a stable room temperature during the tests, which presented another source of error. We predict that the simulated temperature will still be lower than the experimental temperature after 600 s; however, the error will gradually decrease. Overall, the simulation and experimental results show similar increasing trends, which indicates that the model and simulation can be used for reference.

4.3. Overall Device

Figure 13 shows the simulated temperature distribution for the entire solid between 60 and 240 s. To show greater distinction, the upper limit for the displayed temperature was adjusted to 333.15 K (60 °C). The figure clearly shows that heat was conducted through the solid over time and a temperature gradient formed over the entire solid domain. Under the action of the fan, the temperature also increased in nearby objects that were not in contact with the heat source.
As shown in Figure 14, two sections were created to show more detail regarding the temperature inside the casing. The temperatures at z = 0.135 m and x = 0.036 m are shown in Figure 15 and Figure 16, respectively. The temperature field shows that the temperature was highest at the heat source and that the heat spread through the entire space with the aid of the fan. The flow velocities of the two sections and the temperature of the air domain at the last moment of the simulation are shown in Figure 17. The air passed through the fan and over the cooling fins. The flow rate was highest at the fan outlet and gradually decreased as the heat-dissipation structure blocked the air and caused it to diffuse. Figure 17a,b show the cross-sectional flow velocity and the path of heat diffusion. The temperature spread along the wall under the action of the velocity field and finally gathered at the top near the rear-left corner.
The temperatures at points p1–p9 (indicated by red crosses in Figure 18) were used to verify the utility of the simulation. The electronic components and screen control panel in the upper part of the instrument are vulnerable to heat; therefore, the points were selected from this region. The simulation and experimental results are presented in Figure 19. The errors in the simulation were mainly caused by the model simplification. In addition, it was difficult to model the wiring in the instrument for the simulation process. The wires exchange heat with the air, which affects the local airflow and temperature. Hence, there were some differences between the experiments and simulations. The temperature error (obtained by subtracting the experimental value from the simulated value) of each point is shown in Figure 20. The errors were within ±3 K and the maximum relative error for p4 was 0.9%. Therefore, the results of the finite-volume simulation can be used to guide the optimal design of the temperature field for nucleic acid detection instruments.

4.4. Heat-Dissipation Scheme

During operation, the fan circulated the air inside the instrument, and over time the ambient internal temperature increased. High ambient temperatures threaten the stability of the electronic components. The analysis in Section 4.3 showed that the heat was concentrated at the top-back section of the instrument; based on this, two heat-dissipation schemes were designed, as shown in Figure 21. In Plan 1, an inlet was positioned at the bottom of the front of the instrument and an outlet was positioned at the top of the back (Figure 21a). In Plan 2, inlets were positioned around the bottom sides of the instrument and an outlet was positioned at the top (Figure 21b). In both cases, the area of the inlet was approximately equal to the area of the outlet.
There are fans for the amplification modules inside the instrument. According to the pressure difference obtained by simulation, combined with the pressure polynomial in Section 3.2, it is estimated that the velocity is between 5.5–6 m·s−1 (flow rate is between 0.014–0.016 m3·s−1). Considering that the volume of the fluid domain inside the entire instrument was 0.038 m3, a silent fan with a maximum flow rate of 0.028 m3·s−1 was preliminarily selected (Delta Electronics, AFB0924HH). Thus, the pressure polynomial was:
Δ P = 53.935 12.953 v 0.3001 v 2
The temperatures at p1–p9 at the last moment of the simulation were used as the evaluation indexes. Figure 22 shows the temperatures at each point under each scheme. At each point, the temperature under Plan 2 was lower than that under Plan 1; therefore, Plan 2 had a superior heat-dissipation effect. Figure 23 shows the temperature and flow velocity of the section and the overall streamline of the Plan-2 model at the last moment of the simulation. The fluid entered the instrument at a high velocity and then slowed as it became blocked by the internal structure of the device. The flow velocity was increased to remove heat under the action of the fan on the amplification module.
The trends in the simulated and experimental temperatures at points p1–p9 under Plan 2 are shown in Figure 24. The relative error of each point did not exceed ±2 K and the maximum relative error of p9 was 0.54%. The acrylic casing of the prototype based on Plan 2 was used to evaluate the heat-dissipation capability. Figure 25a shows the prototype of the casing and instrument body. Nucleic acid amplification was conducted continuously for 120 min and the temperatures of p1–p9 were measured simultaneously, as shown in Figure 25b. The maximum temperature was 307.76 K. During the first amplification process, the temperature increased rapidly and then stabilized in subsequent amplification cycles. After each amplification, the temperature decreased significantly, which indicates that the residual heat was removed from the instrument immediately after each amplification. This shows that the proposed heat-dissipation scheme for nucleic acid detectors is feasible.

5. Conclusions

In this study, a vent scheme was designed to provide a safe working environment for the reliable operation of a nucleic acid detector. A simplified 3D model of the instrument was produced and the finite-volume method was used to simulate the temperature and velocity fields within it. The simulation showed that heat accumulated in certain areas of the instrument. A temperature-measurement platform was constructed to measure the temperature experimentally and the results were used to validate the simulation and its accuracy was confirmed. A comparison between the experimental and simulation results showed that the temperature error did not exceed ± 3 K and that the maximum relative error was 0.9%. Therefore, the finite-volume simulation results can be used to guide the design of the temperature field for nucleic acid detectors. A total of two ventilation schemes were designed based on the simulation results and the best scheme was employed to fabricate a new casing for the instrument. In a long-term nucleic acid amplification experiment, the monitoring points reached a maximum temperature of 307.76 K. Therefore, the vents in the designed casing showed excellent heat-dissipation performance; hence, they can remove waste heat and create a safe operating environment for the instrument.

Author Contributions

Methodology, X.L.; Validation, S.H.; Formal analysis, W.C.; Investigation, W.C.; Writing—original draft, S.Z.; Writing—review & editing, X.L.; Project administration, D.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China, grant number 62003284; Natural Science Foundation of Fujian Province, China, grant number 2021J011200; The presidential research foundation of Xiamen University grant number 20720210089.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

We would like to thank the editors and anonymous reviewers for their reviews and comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wu, Y.; Peng, W.; Zhao, Q.; Piao, J.; Zhang, B.; Wu, X.; Wang, H.; Shi, Z.; Gong, X.; Chang, J. Immune fluorescence test strips based on quantum dots for rapid and quantitative detection of carcino-embryonic antigen. Chin. Chem. Lett. 2017, 28, 1881–1884. [Google Scholar] [CrossRef]
  2. Li, H.T.; Pan, J.Z.; Fang, Q. Development and Application of Digital PCR Technology. Prog. Chem. 2020, 32, 581–593. [Google Scholar]
  3. Mackay, I.M.; Arden, K.E.; Nitsche, A. Real-time PCR in virology. Nucleic Acids Res. 2002, 30, 1292–1305. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Valasek, M.A.; Repa, J.J. The power of real-time PCR. Adv. Physiol. Educ. 2005, 29, 151–159. [Google Scholar] [CrossRef]
  5. Klein, D. Quantification using real-time PCR technology: Applications and limitations. Trends Mol. Med. 2002, 8, 257–260. [Google Scholar] [CrossRef]
  6. Lall, P.; Pecht, M.; Hakim, E.B. Influence of Temperature on Microelectronics and System Reliability: A Physics of Failure Approach; Taylor & Francis Group: Florida, FL, USA, 1997. [Google Scholar]
  7. Cooper, M. Investigation of Arrhenius acceleration factor for integrated circuit early life failure region with several failure mechanisms. IEEE Trans. Compon. Packag. Technol. 2005, 28, 561–563. [Google Scholar] [CrossRef]
  8. Lall, P.; Pecht, M.; Hakim, E.B. Characterization of functional relationship between temperature and microelectronic reliability. Microelectron. Reliab. 1995, 35, 377–402. [Google Scholar] [CrossRef]
  9. Son, J.H.; Hong, S.G.; Haack, A.J.; Gustafson, L.; Song, M.S.; Hoxha, O.; Lee, L.P. Rapid Optical Cavity PCR. Adv. Healthc. Mater. 2016, 5, 167–174. [Google Scholar] [CrossRef]
  10. Li, J.; Liu, S.; Li, S.; Zhang, P.; Liu, H.; Shao, L. PCR Instrument Temperature Control System Based on Multimodal control. In Proceedings of the 2020 Chinese Control And Decision Conference (CCDC), Hefei, China, 22–24 August 2020; pp. 149–154. [Google Scholar] [CrossRef]
  11. Jalili, A.; Bagheri, M.; Shamloo, A.; Ashkezari, A.H.K. A plasmonic gold nanofilm-based microfluidic chip for rapid and inexpensive droplet-based photonic PCR. Sci. Rep. 2021, 11, 1–13. [Google Scholar] [CrossRef]
  12. Ramezani, A. CtNorm: Real time PCR cycle of threshold (Ct) normalization algorithm. J. Microbiol. Methods 2021, 187, 106267. [Google Scholar] [CrossRef]
  13. Nomada, H.; Morita, K.; Higuchi, H.; Yoshioka, H.; Oki, Y. Carbon–polydimethylsiloxane-based integratable optical technology for spectroscopic analysis. Talanta 2017, 166, 428–432. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Aziz, I.; Jamshaid, R.; Zaidiand, T.; Akhtar, I. Numerical simulation of heat transfer to optimize DNA amplification in Polymerase Chain Reaction. In Proceedings of the 2016 International Conference on Applied Mathematics, Simulation and Modelling, Beijing, China, 28–29 May 2016; pp. 456–462. [Google Scholar] [CrossRef]
  15. Naghdloo, A.; Ghazimirsaeed, E.; Shamloo, A. Numerical simulation of mixing and heat transfer in an integrated centrifugal microfluidic system for nested-PCR amplification and gene detection. Sens. Actuators B Chem. 2018, 283, 831–841. [Google Scholar] [CrossRef]
  16. Parng, S.-H.; Wu, P.-J.; Tsai, Y.-Y.; Hong, R.-S.; Lee, S.-J. A 3D tubular structure with droplet generation and temperature control for DNA amplification. Microfluid. Nanofluidics 2021, 25, 1–8. [Google Scholar] [CrossRef]
  17. Zhou, A.; Liu, B.; Xu, X.; Li, W. An Analysis of the Two-Dimensional Heat Transfer Model of the PCR Instrument Pedestal. In Proceedings of the An Analysis of the Two-Dimensional Heat Transfer Model of the PCR Instrument Pedestal, Shanghai, China, 28–29 May 2016. [Google Scholar] [CrossRef] [Green Version]
  18. Raza, M.; Malik, I.R.; Akhtar, I. Design and Thermal Analysis of Sample Block PCR for DNA Amplification. In Proceedings of the 2020 17th International Bhurban Conference on Applied Sciences and Technology, Islamabad, Pakistan, 14–18 January 2020; pp. 210–217. [Google Scholar]
  19. Hao, H.; Guan, X.; Chen, J.; Pan, Z. Internal cooling design of PCB based on thermal model. Electronic Compon. Mater. 2016, 35, 76–79. [Google Scholar]
  20. Mesforush, S.; Jahanshahi, A.; Zadeh, M.K. Finite element simulation of isothermal regions in serpentine shaped PCB electrodes of a micro-PCR device. In Proceedings of the 2019 27th Iranian Conference on Electrical Engineering, Yazd, Iran, 30 April 2019–2 May 2019; pp. 280–284. [Google Scholar] [CrossRef]
  21. Fedorov, A.G.; Viskanta, R. Three-dimensional conjugate heat transfer in the microchannel heat sink for electronic packaging. Int. J. Heat Mass Transf. 2000, 43, 399–415. [Google Scholar] [CrossRef]
  22. Shankaran, G.V.; Singh, R.K. Selection of Appopriate Thermal Model For Printed Circut Boards in CFD Analysis. In Proceedings of the thermal and Thermomechanical Proceedings 10th Intersociety Conference on Phenomena in Electronics Systems, San Diego, CA, USA, 30 May–2 June 2006; pp. 234–242. [Google Scholar]
  23. Jun, Y.; Dakui, F.; Hang, Z. Numerical simulation of pump jet propulsor and analysis of scale effect. Shipbuild. China 2020, 61, 91–99. [Google Scholar]
  24. Versteeg, H.K.; Malalasekera, W. An Introduction to Computational Fluid Dynamics: The Finite Volume Method; Pearson Education: New York, NY, USA, 2007. [Google Scholar]
  25. Behera, R.C.; Dutta, P.; Srinivasan, K. Numerical Study of Interrupted Impinging Jets for Cooling of Electronics. IEEE Trans. Compon. Packag. Technol. 2007, 30, 275–284. [Google Scholar] [CrossRef] [Green Version]
  26. Boulet, M.; Marcos, B.; Dostie, M.; Moresoli, C. CFD modeling of heat transfer and flow field in a bakery pilot oven. J. Food Eng. 2010, 97, 393–402. [Google Scholar] [CrossRef]
  27. Jo, J.C.; Kang, D.G. CFD analysis of thermally stratified flow and conjugate heat transfer in a PWR pressurizer surgeline. J. Press. Vessel. Technol. 2010, 132, 021301. [Google Scholar] [CrossRef]
  28. Widiawaty, C.D.; Siswantara, A.I.; Budiarso; Gunadi, G.G.R.; Pujowidodo, H.; Syafei, M.H.G. A CFD simulation and experimental study: Predicting heat transfer performance using SST k-ω turbulence model. IOP Conf. Ser. Mater. Sci. Eng. 2020, 909, 012004. [Google Scholar] [CrossRef]
  29. Han, J.W.; Zhao, C.J.; Yang, X.T.; Qian, J.P.; Fan, B.L. Computational modeling of airflow and heat transfer in a vented box during cooling: Optimal package design. Appl. Therm. Eng. 2015, 91, 883–893. [Google Scholar] [CrossRef]
  30. Menter, F.R. Two-equation eddy-viscosity turbulence models for engineering applications. AIAA J. 1994, 32, 1598–1605. [Google Scholar] [CrossRef] [Green Version]
  31. Menter, F.R.; Kuntz, M.; Langtry, R. Ten years of industrial experience with the SST turbulence model. Turbul. Heat Mass Transf. 2003, 4, 625–632. [Google Scholar]
  32. Wilcox, D.C. Turbulence Modeling for CFD; DCW industries: La Canada, CA, USA, 1998. [Google Scholar]
  33. Król, A.; Król, M. Study on numerical modeling of jet fans. Tunn. Undergr. Space Technol. 2018, 73, 222–235. [Google Scholar] [CrossRef]
  34. Fan, J.; Dong, H.; Xu, X.; Teng, D.; Yan, B.; Zhao, Y. Numerical Investigation on the Influence of Mechanical Draft Wet-Cooling Towers on the Cooling Performance of Air-Cooled Condenser with Complex Building Environment. Energies 2019, 12, 4560. [Google Scholar] [CrossRef] [Green Version]
  35. Paramanandam, K.; Venkatachalapathy, S.; Srinivasan, B. Heat transfer enhancement in microchannels using ribs and secondary flows. Int. J. Numer. Methods Heat Fluid Flow 2021, 32, 2299–2319. [Google Scholar] [CrossRef]
  36. Luo, H.; Lu, B.; Zhang, J.; Wu, H.; Wang, W. A grid-independent EMMS/bubbling drag model for bubbling and turbulent fluidization. Chem. Eng. J. 2017, 326, 47–57. [Google Scholar] [CrossRef] [Green Version]
  37. Zaidan, M.H.; Alkumait, A.A.; Ibrahim, T.K. Assessment of heat transfer and fluid flow characteristics within finned flat tube. Case Stud. Therm. Eng. 2018, 12, 557–562. [Google Scholar] [CrossRef]
Figure 1. Structure of the portable automatic nucleic acid detector: (a) original 3D model; (b) simulation model; (c) amplification module; (d) heat transfer form of amplification module.
Figure 1. Structure of the portable automatic nucleic acid detector: (a) original 3D model; (b) simulation model; (c) amplification module; (d) heat transfer form of amplification module.
Applsci 12 11966 g001aApplsci 12 11966 g001b
Figure 2. Overview of test platform.
Figure 2. Overview of test platform.
Applsci 12 11966 g002
Figure 3. Imaging showing locations of the heat sources: (a) Location of heat source inside the instrument; (b) Location of heat source on PCB.
Figure 3. Imaging showing locations of the heat sources: (a) Location of heat source inside the instrument; (b) Location of heat source on PCB.
Applsci 12 11966 g003
Figure 4. Measurement results of the heat sources: (a) TEC; (b) PCB1 and PCB2.
Figure 4. Measurement results of the heat sources: (a) TEC; (b) PCB1 and PCB2.
Applsci 12 11966 g004
Figure 5. Fluid domain, solid domain, and heat source.
Figure 5. Fluid domain, solid domain, and heat source.
Applsci 12 11966 g005
Figure 6. Fitting curve of the fan.
Figure 6. Fitting curve of the fan.
Applsci 12 11966 g006
Figure 7. Grid independence verification.
Figure 7. Grid independence verification.
Applsci 12 11966 g007
Figure 8. Grid with maximum element size of 0.005 m.
Figure 8. Grid with maximum element size of 0.005 m.
Applsci 12 11966 g008
Figure 9. Simulation results of amplification module: (a) Overall temperature distribution; (b) Temperature distribution on the back of heat-dissipation structure; (c) Temperature distribution on the cross section.
Figure 9. Simulation results of amplification module: (a) Overall temperature distribution; (b) Temperature distribution on the back of heat-dissipation structure; (c) Temperature distribution on the cross section.
Applsci 12 11966 g009
Figure 10. Comparison of the simulation and experimental results for the amplification module.
Figure 10. Comparison of the simulation and experimental results for the amplification module.
Applsci 12 11966 g010
Figure 11. Experimental and simulated temperature distributions: (a) PCB2; (b) PCB1.
Figure 11. Experimental and simulated temperature distributions: (a) PCB2; (b) PCB1.
Applsci 12 11966 g011
Figure 12. Comparison of the simulation and experimental results for PCBs: (a) PCB2; (b) PCB1.
Figure 12. Comparison of the simulation and experimental results for PCBs: (a) PCB2; (b) PCB1.
Applsci 12 11966 g012
Figure 13. Temperature distribution for the entire solid.
Figure 13. Temperature distribution for the entire solid.
Applsci 12 11966 g013aApplsci 12 11966 g013b
Figure 14. Cross sections at x = 0.038 m and z = 0.135 m.
Figure 14. Cross sections at x = 0.038 m and z = 0.135 m.
Applsci 12 11966 g014
Figure 15. Section temperature at z = 0.135 m.
Figure 15. Section temperature at z = 0.135 m.
Applsci 12 11966 g015
Figure 16. Section temperature at x = 0.038 m.
Figure 16. Section temperature at x = 0.038 m.
Applsci 12 11966 g016
Figure 17. (a) Section velocity at z = 0.135 m; (b) section velocity at x = 0.038 m; (c,d) temperature of the air domain.
Figure 17. (a) Section velocity at z = 0.135 m; (b) section velocity at x = 0.038 m; (c,d) temperature of the air domain.
Applsci 12 11966 g017
Figure 18. Schematic diagram showing the locations of the verification points.
Figure 18. Schematic diagram showing the locations of the verification points.
Applsci 12 11966 g018
Figure 19. Simulation and experimental temperatures at each verification point over time.
Figure 19. Simulation and experimental temperatures at each verification point over time.
Applsci 12 11966 g019
Figure 20. Relative deviation in temperature at verification points p1–p9.
Figure 20. Relative deviation in temperature at verification points p1–p9.
Applsci 12 11966 g020
Figure 21. Heat-dissipation schemes: (a) plan 1; (b) plan 2.
Figure 21. Heat-dissipation schemes: (a) plan 1; (b) plan 2.
Applsci 12 11966 g021
Figure 22. Temperature at points p1–p9 under the different schemes.
Figure 22. Temperature at points p1–p9 under the different schemes.
Applsci 12 11966 g022
Figure 23. (a) Section temperature at x = 0.038 m; (b) section velocity at x = 0.038 m; (c) overall streamline.
Figure 23. (a) Section temperature at x = 0.038 m; (b) section velocity at x = 0.038 m; (c) overall streamline.
Applsci 12 11966 g023
Figure 24. Variation trends in simulation value and experimental value of points p1–p9.
Figure 24. Variation trends in simulation value and experimental value of points p1–p9.
Applsci 12 11966 g024
Figure 25. (a) Prototype; (b) temperature curve of points p1–p9.
Figure 25. (a) Prototype; (b) temperature curve of points p1–p9.
Applsci 12 11966 g025
Table 1. Parameters of the PCBs.
Table 1. Parameters of the PCBs.
δ p (m) A 1 (m2) A 2 (m2) A 3 (m2) A 4 (m2)
PCB10.00160.0120.0120.0120.01
PCB20.00160.0230.020.0220.019
Table 2. Parameters of the PCBs.
Table 2. Parameters of the PCBs.
δ c u (m) δ F (m) A P (m2) k c u (W·m−1·K−1) k F (W·m−1·K−1) N
PCB10.0000350.00004870.0133880.354
PCB20.0000350.00004870.0243880.354
Table 3. Material parameters.
Table 3. Material parameters.
MaterialDensity
(kg·m−3)
Specific Heat
(J·kg−1·K−1)
Thermal Conductivity
(W·m−1·K−1)
6061 aluminum alloy2750896155
Stainless steel805548013.8
Packaging plastic20909000.67
FR4190013690.3
Pom14301.4650.23
Si2329700130
Ceramics1345710150
Acrylic119014640.18
Red copper8960385401
Wood70023100.173
Table 4. Simulation settings.
Table 4. Simulation settings.
ObjectType
Fluid and solid contact surfacescouple wall
Fan boundarypolynomial
TEC side1 and side2temperature wall
Heat source on the PCBvalues
Fluidair, incompressible
Inlet and outletpressure inlet/outlet
Pressure–velocity couplesimple
The discretization of convection–diffusionsecond-order upwind
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, X.; Cai, W.; Huang, S.; Zhu, S.; Zhang, D. Conjugate Heat Transfer Analysis and Heat Dissipation Design of Nucleic Acid Detector Instrument. Appl. Sci. 2022, 12, 11966. https://doi.org/10.3390/app122311966

AMA Style

Lin X, Cai W, Huang S, Zhu S, Zhang D. Conjugate Heat Transfer Analysis and Heat Dissipation Design of Nucleic Acid Detector Instrument. Applied Sciences. 2022; 12(23):11966. https://doi.org/10.3390/app122311966

Chicago/Turabian Style

Lin, Xiaohui, Weihuang Cai, Shaolei Huang, Sijie Zhu, and Dongxu Zhang. 2022. "Conjugate Heat Transfer Analysis and Heat Dissipation Design of Nucleic Acid Detector Instrument" Applied Sciences 12, no. 23: 11966. https://doi.org/10.3390/app122311966

APA Style

Lin, X., Cai, W., Huang, S., Zhu, S., & Zhang, D. (2022). Conjugate Heat Transfer Analysis and Heat Dissipation Design of Nucleic Acid Detector Instrument. Applied Sciences, 12(23), 11966. https://doi.org/10.3390/app122311966

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