No Arabic abstract
The combination of microstructures and mixed wettability for enhancing nucleate boiling has attracted much attention in recent years. However, in the existing experimental and numerical studies, the tops of microstructures are entirely subjected to wettability modification, which makes the influences of mixed wettability dependant on the characteristic length of microstructures. In order to disclose the joint effects of surface structure and mixed wettability on nucleate boiling, in this work we propose an improved type of pillar-textured surface with mixed wettability, in which the tops of square pillars are partially subjected to wettability modification. Numerical investigation of the boiling heat transfer performance on the improved mixed-wettability surface is carried out using a three-dimensional thermal multiphase lattice Boltzmann model. The numerical results show that the width of the wettability-modified region plays an important role in the boiling performance of the improved mixed-wettability surface and the best boiling performance is achieved in the situation that the width of the wettability-modified region is sufficiently large but the bubble nucleated on the pillar top still does not interfere with the coalescence-departure mechanism of the bubbles nucleated around the pillar, which optimizes the joint effects of surface structure and mixed wettability for enhancing nucleate boiling. The influences of the shape of the wettability-modified region are also studied. Among the investigated shapes, the square is found to perform better than the other two shapes.
A hybrid thermal lattice Boltzmann (LB) model is presented to simulate thermal multiphase flows with phase change based on an improved pseudopotential LB approach [Q. Li, K. H. Luo, and X. J. Li, Phys. Rev. E 87, 053301 (2013)]. The present model does not suffer from the spurious term caused by the forcing-term effect, which was encountered in some previous thermal LB models for liquid-vapor phase change. Using the model, the liquid-vapor boiling process is simulated. The boiling curve together with the three boiling stages (nucleate boiling, transition boiling, and film boiling) is numerically reproduced in the LB community for the first time. The numerical results show that the basic features and the fundamental characteristics of boiling heat transfer are well captured, such as the severe fluctuation of transient heat flux in the transition boiling and the feature that the maximum heat transfer coefficient lies at a lower wall superheat than that of the maximum heat flux. Furthermore, the effects of the heating surface wettability on boiling heat transfer are investigated. It is found that an increase in contact angle promotes the onset of boiling but reduces the critical heat flux, and makes the boiling process enter into the film boiling regime at a lower wall superheat, which is consistent with the findings from experimental studies.
When boiling occurs in a liquid flow field, the phenomenon is known as forced-convection boiling. We numerically investigate such a boiling system on a cylinder in a flow at a saturated condition. To deal with the complicated liquid-vapor phase-change phenomenon, we develop a numerical scheme based on the pseudopotential lattice Boltzmann method (LBM). The collision stage is performed in the space of central moments (CMs) to enhance numerical stability for high Reynolds numbers. The adopted forcing scheme, consistent with the CMs-based LBM, leads to a concise yet robust algorithm. Furthermore, additional terms required to ensure thermodynamic consistency are derived in a CMs framework. The effectiveness of the present scheme is successfully tested against a series of boiling processes, including nucleation, growth, and departure of a vapor bubble for Reynolds numbers varying between 30 and 30000. Our CMs-based LBM can reproduce all the boiling regimes, i.e., nucleate boiling, transition boiling, and film boiling, without any artificial input such as initial vapor phase. We find that the typical boiling curve, also known as the Nukiyama curve, appears even though the focused system is not the pool boiling but the forced-convection system. Also, our simulations support experimental observations of intermittent direct solid-liquid contact even in the film-boiling regime. Finally, we provide quantitative comparison with the semi-empirical correlations for the forced-convection film boiling on a cylinder on the Nu-Ja diagram.
In recent years, the lattice Boltzmann (LB) method has been widely employed to simulate boiling phenomena [A. Markus and G. Hazi, Phys. Rev. E 83, 046705 (2011); Biferale et al., Phys. Rev. Lett. 108, 104502 (2012); Li et al., Phys. Rev. E 96, 063303 (2017); Wu et al., Int. J. Heat Mass Transfer 126, 773 (2018)]. However, a very important issue still remains open, i.e., how does boiling occur in the LB simulations? For instance, the existing LB studies showed that the boiling on a hydrophobic surface begins at a lower wall superheat than that on a hydrophilic surface, which qualitatively agrees well with experimental studies, but no one has yet explained how this phenomenon appears in the LB simulations and what happened in the simulations after changing the wettability of the heating surface. In this paper, the LB boiling mechanism is revealed by analyzing boiling on a flat surface with mixed wettability and boiling on a structured surface with homogeneous wettability. Through a theoretical analysis, we demonstrate that, when the same wall superheat is applied, in the LB boiling simulations the fluid density near the heating surface decreases faster on a hydrophobic surface than that on a hydrophilic surface. Accordingly, a lower wall superheat can induce the phase transition from liquid to vapor on a hydrophobic surface than that on a hydrophilic surface. Furthermore, a similar theoretical analysis shows that the fluid density decreases fastest at concave corners in the case of a structured surface with homogeneous wettability, which explains why vapor bubbles are nucleated at concave corners in the LB simulations of boiling on structured surfaces.
Based on mesoscale lattice Boltzmann (LB) numerical simulations, we investigate the effects of viscoelasticity on the break-up of liquid threads in microfluidic cross-junctions, where droplets are formed by focusing a liquid thread of a dispersed (d) phase into another co-flowing continuous (c) immiscible phase. Working at small Capillary numbers, we investigate the effects of non-Newtonian phases in the transition from droplet formation at the cross-junction (DCJ) to droplet formation downstream of the cross-junction (DC) (Liu $&$ Zhang, ${it Phys. ~Fluids.}$ ${bf 23}$, 082101 (2011)). We will analyze cases with ${it Droplet ~Viscoelasticity}$ (DV), where viscoelastic properties are confined in the dispersed phase, as well as cases with ${it Matrix ~Viscoelasticity}$ (MV), where viscoelastic properties are confined in the continuous phase. Moderate flow-rate ratios $Q approx {cal O}(1)$ of the two phases are considered in the present study. Overall, we find that the effects are more pronounced in the case with MV, where viscoelasticity is found to influence the break-up point of the threads, which moves closer to the cross-junction and stabilizes. This is attributed to an increase of the polymer feedback stress forming in the corner flows, where the side channels of the device meet the main channel. Quantitative predictions on the break-up point of the threads are provided as a function of the Deborah number, i.e. the dimensionless number measuring the importance of viscoelasticity with respect to Capillary forces.
The effects of viscoelasticity on the dynamics and break-up of fluid threads in microfluidic T-junctions are investigated using numerical simulations of dilute polymer solutions at changing the Capillary number ($mbox {Ca}$), i.e. at changing the balance between the viscous forces and the surface tension at the interface, up to $mbox{Ca} approx 3 times 10^{-2}$. A Navier-Stokes (NS) description of the solvent based on the lattice Boltzmann models (LBM) is here coupled to constitutive equations for finite extensible non-linear elastic dumbbells with the closure proposed by Peterlin (FENE-P model). We present the results of three-dimensional simulations in a range of $mbox{Ca}$ which is broad enough to characterize all the three characteristic mechanisms of breakup in the confined T-junction, i.e. ${it squeezing}$, ${it dripping}$ and ${it jetting}$ regimes. The various model parameters of the FENE-P constitutive equations, including the polymer relaxation time $tau_P$ and the finite extensibility parameter $L^2$, are changed to provide quantitative details on how the dynamics and break-up properties are affected by viscoelasticity. We will analyze cases with ${it Droplet ~Viscoelasticity}$ (DV), where viscoelastic properties are confined in the dispersed (d) phase, as well as cases with ${it Matrix ~Viscoelasticity}$ (MV), where viscoelastic properties are confined in the continuous (c) phase. Moderate flow-rate ratios $Q approx {cal O}(1)$ of the two phases are considered in the present study. Overall, we find that the effects are more pronounced in the case with MV, as the flow driving the break-up process upstream of the emerging thread can be sensibly perturbed by the polymer stresses.