scholarly journals Comparison of the Volume of Fluid and CLSVOF Methods for the Assessment of Flow Boiling in Silicon Microgaps

2017 ◽  
Vol 139 (11) ◽  
Author(s):  
Daniel Lorenzini ◽  
Yogendra Joshi

The three-dimensional (3D) stacking of integrated circuits (ICs), and emergent microelectronic technologies require low-profile cooling solutions for the removal of relatively high heat fluxes. The flow boiling of dielectric refrigerants represents a feasible alternative to such applications by providing compatibility with the electrical interconnections, relatively uniform temperature profiles, and higher heat transfer coefficients than those obtained with single phase-cooling. Despite important experimental evidence in this area has been recently reported in the literature, the modeling of such has remained in basic and limited forms due to the associated complexities with the physics of two-phase flow with phase-change. In an effort to expand the studied possibilities for the modeling of flow boiling, the present investigation compares two different phase-tracking methods for the analysis of such phenomena: the volume of fluid (VOF) and the coupled level set—volume of fluid (CLSVOF) techniques. These interface tracking and reconstruction techniques are coupled with a phase change model that accounts for the mass and energy transfer source terms to the governing equations. The geometric domain is constituted by a silicon microgap 175 μm high with a substrate thickness of 50 μm, and populated with circular pin fins of 150 μm diameter, where the heat conduction is simultaneously solved with temperature dependent properties. The flow boiling regimes and their spatial and temporal evolution are compared between both methods by maintaining the operating conditions. Results indicate that both methods provide a good capability to predict major two-phase flow regimes observed in experimental studies with these types of arrangements. However, the CLSVOF offers a sharper interface reconstruction than the standard VOF method by predicting bubble nucleation and departure mechanisms more closely to experimental observations.

Author(s):  
Nishant Tiwari ◽  
Manoj Kumar Moharana

Flow boiling in microchannel heat sink offers an effective cooling solution for high power density micro devices. A three-dimensional numerical study based on volume of fraction model (VOF) coupled with evaporation condensation model accounting for the liquid-vapor phase change is undertaken to recreate vapor bubble formation in saturated flow boiling in wavy microchannel. Constant wall heat flux imposed at the bottom surface of the substrate while other faces are insulated. To understand the conjugate effects, simulations has been carried out for substrate thickness to channel depth ratio (δsf ∼ 1–5), substrate wall to fluid thermal conductivity ratio (ksf ∼ 22–300) and waviness (γ ∼ 0.008–0.04). Bubble nucleation, growth, and departure of bubble plays a significant role in heat transfer and pressure drop characteristics in two-phase flow in wavy microchannel. The coolant (water) temperature at the inlet is taken to be 373 K while flow was at atmospheric pressure. This makes shorter waiting period of bubble nucleation, and the number density of bubbles on the solid surface increases. This results in enhancement of the boiling effect, and thus with the presence of bubbles, the mixing of laminar boundary layers improves and enhances the overall heat transfer coefficient. Channel amplitude play an important factor that can suitably reduce the friction factor and enhances the heat transfer coefficient.


Author(s):  
Jacqueline Barber ◽  
Khellil Sefiane ◽  
David Brutin ◽  
Lounes Tadrist

Boiling in microchannels remains elusive due to the lack of full understanding of the mechanisms involved. A powerful tool in achieving better comprehension of the mechanisms is detailed imaging and analysis of the two phase flow at a fundamental level. We induced boiling in a single microchannel geometry (hydraulic diameter 727 μm), using a refrigerant FC-72, to investigate several flow patterns. A transparent, metallic, conductive deposit has been developed on the exterior of rectangular microchannels, allowing simultaneous uniform heating and visualisation to be conducted. The data presented in this paper is for a particular case with a uniform heat flux of 4.26 kW/m2 applied to the microchannel and inlet liquid mass flowrate, held constant at 1.33×10−5 kg/s. In conjunction with obtaining high-speed images and videos, sensitive pressure sensors are used to record the pressure drop profiles across the microchannel over time. Bubble nucleation, growth and coalescence, as well as periodic slug flow, are observed in the test section. Phenomena are noted, such as the aspect ratio and Reynolds number of a vapour bubble, which are in turn correlated to the associated pressure drops over time. From analysis of our results, images and video sequences with the corresponding physical data obtained, it is possible to follow visually the nucleation and subsequent both ‘free’ and ‘confined’ growth of a vapour bubble over time.


2017 ◽  
Vol 140 (1) ◽  
Author(s):  
Daniel Lorenzini ◽  
Yogendra K. Joshi

The computational fluid dynamics (CFD) modeling of boiling phenomena has remained a challenge due to numerical limitations for accurately simulating the two-phase flow and phase-change processes. In the present investigation, a CFD approach for such analysis is described using a three-dimensional (3D) volume of fluid (VOF) model coupled with a phase-change model accounting for the interfacial mass and energy transfer. This type of modeling allows the transient analysis of flow boiling mechanisms, while providing the ability to visualize in detail temperature, phase, and pressure distributions for microscale applications with affordable computational resources. Results for a plain microchannel are validated against benchmark correlations for heat transfer (HT) coefficients and pressure drop as a function of the heat flux and mass flux. Furthermore, the model is used for the assessment of two-phase cooling in microelectronics under a realistic scenario with nonuniform heat fluxes at localized regions of a silicon microchannel, relevant to the cooling layer of 3D integrated circuit (IC) architectures. Results indicate the strong effect of two-phase flow regime evolution and vapor accumulation on HT. The effects of reduced saturation pressure, subcooling, and flow arrangement are explored in order to provide insight about the underlying physics and cooling performance.


2021 ◽  
Vol 895 (1) ◽  
pp. 012002
Author(s):  
V S Alekseev ◽  
R S Seryi

Abstract Currently sluice washing devices are the most common in alluvial gold mining. Their use provides a sufficiently high performance, relatively low power consumption, and acceptable recovery of valuable components. The theoretical provisions of traditional hydraulics make it possible to determine all the main parameters of the movement of particles of rocks and gold in the pulp, however, in real operating conditions of the sluice box, their actual values will differ greatly from the calculated ones, especially if there are solid fractions in the pulp with a particle size of more than 20 mm. This is explained by significant fluctuations in the values of the surface, average and bottom velocities of the two-phase flow, vertical pulsation velocity in conditions of constrained movement of the different fractional composition of rocks. The article presents the results of experimental studies to identify the dependence of the distance traveled by an individual gold particle and host rocks in a two-phase flow through a sluice, the bottom of which is lined with trapping coatings, on the design and technological parameters of the flushing device. The mathematical model for determining this distance formed the basis of the Gold Enriching program. The program allows, in a wide range of initial data, to determine the zones of concentration of gold of a certain size at the sluice boxes.


Author(s):  
Gae¨l Gue´don ◽  
Emanuela Colombo ◽  
Fabio Inzoli

Several engineering plants and components, such as nuclear reactors, boilers, bubble columns, chemical reactors or oil/gas wells may be characterized by two-phase flows. Appropriate and safe operation of these systems may be supported by the accurate prediction of the multiphase flow pattern with consistent estimation of the void fraction and flooding phenomenon. A preliminary approach for modeling flow patterns in countercurrent two-phase flows in vertical systems, based on the volume-of-fluid (VOF) method, is here presented. The general objective of the study is to investigate the feasibility of large scale two-phase systems simulations using the VOF method. The specific objective is to provide a first set of qualitative information about the fluid dynamics structures in this flow at specific operating conditions. Three-dimensional simulations are performed using a finite volume commercial CFD code. Validation of the numerical approach is achieved with experimental data taken from literature for typical air-water flows in bubble columns. The interaction of the downward water recirculation with the bubbles of air is indeed representative of a bubbly countercurrent two-phase flow and therefore may be a valuable test case. The use of the VOF method is here privileged, since no further closure relations are needed for exchange coefficients between the continuous and dispersed phases, and for breakup and coalescence of bubbles/drops. The validated model is therefore most likely applicable in complex situations where the flow behavior is unknown. Finally a specific analysis of a countercurrent two-phase flow is presented with the objective of simulating the flow within an industrial pipe where two immiscible fluids, with different density are injected. A fluid “A” is injected at the bottom of the pipe and it is supposed to exit the pipe in the upper part. A fluid “B”, immiscible in “A”, is injected at the top of the pipe in countercurrent and it is supposed not to influence the path of fluid “A” from inlet to the exit. Fluid “B” may be used for reacting somehow with fluid “A”, but in this paper only the fluid dynamic condition is considered.


Energies ◽  
2019 ◽  
Vol 12 (23) ◽  
pp. 4483 ◽  
Author(s):  
Rita Mastrullo ◽  
Alfonso William Mauro

This paper presents a critical analysis of possible data reduction procedures for the evaluation of local heat transfer coefficient during flow boiling experiments. The benchmark method using one-dimensional (1-D) heat transfer in a heated tube was compared to a new data reduction method in which both radial and circumferential contributions to the conductive heat transfer inside a metal tube are considered. Using published experimental flow boiling data, the circumferential profiles of the wall superheat, inner wall heat flux, and heat transfer coefficients were independently calculated with the two data reduction procedures. The differences between the two methods were then examined according to the different heat transfer behavior observed (symmetric or asymmetric), which in turn was related to the two-phase flow regimes occurring in a channel during evaporation. A statistical analysis using the mean absolute percentage error (MAPE) index was then performed for a database of 417 collected flow boiling data taken under different operating conditions in terms of working fluid, saturation temperature, mass velocity, vapor quality, and imposed heat flux. Results showed that the maximum deviations between the two methods could reach up to 130% in the case of asymmetric heat transfer. Finally, the possible uses of the two data reduction methods are discussed, pointing out that the two-dimensional (2-D) model is the most reliable method to be employed in the case of high-level modeling of two-phase flow or advanced design of heat exchangers and heat spreader systems.


Author(s):  
Novitrian ◽  
Minoru Takahashi

Experimental studies of Pb-Bi-Water direct contact boiling two-phase flow have been conducted (in the chimney 29.9mm in diameter and 2,367mm in length). The rated operating conditions of system pressure was 7MPa, the liquid metal temperature was 460°C and injected water temperature was 220°C. Pb-Bi flow rates were measured and analytically estimated using balance condition between buoyancy force and pressure losses. The drift-flux model was used to determine the void fraction with assumption that bubbles sizes were dependent on the fluid surface tension and the density ratio of Pb-Bi to steam and water mixture. The result of analysis shows that the deviation from the experimental data was less than 20%.


Sign in / Sign up

Export Citation Format

Share Document