Pore-scale analysis of the minimum liquid film thickness around elongated bubbles in confined gas-liquid flows

© 2017 Elsevier Ltd The fluid mechanics of elongated bubbles in confined gas-liquid flows in micro-geometries is important in pore-scale flow processes for enhanced oil recovery and mobilization of colloids in unsaturated soil. The efficiency of such processes is traditionally related to the thickness of the liquid film trapped between the elongated bubble and the pore's wall, which is assumed constant. However, the surface of long bubbles presents undulations in the vicinity of the rear meniscus, which may significantly decrease the local thickness of the liquid film, thus impacting the process of interest. This study presents a systematic analysis of these undulations and the minimum film thickness induced in the range Ca=0.001-0.5 and Re=0.1-2000. Pore-scale Computational Fluid Dynamics (CFD) simulations are performed with a self-improved version of the opensource solver ESI OpenFOAM which is based on a Volume of Fluid method to track the gas-liquid interface. A lubrication model based on the extension of the classical axisymmetric Bretherton theory is utilized to better understand the CFD results. The profiles of the rear meniscus of the bubble obtained with the lubrication model agree fairly well with those extracted from the CFD simulations. This study shows that the Weber number of the flow, We=CaRe, is the parameter that best describes the dynamics of the interfacial waves. When We 0.1, a larger number of wave crests becomes evident on the surface of the rear meniscus of the bubble. The liquid film thickness at the crests of the undulations thins considerably as the Reynolds number is increased, down to less than 60% of the value measured in the flat film region. This may significantly influence important environmental processes, such as the detachment and mobilization of micron-sized pollutants and pathogenic micro-organisms adhering at the pore's wall in unsaturated soil

Morphology of long gas bubbles propagating in square capillaries

We present the results of a systematic analysis of the morphology of the thin lubrication film surround- ing a long gas bubble transported by a liquid flow in a square capillary. Direct numerical simulations of the flow are performed using the Volume-Of-Fluid method implemented in OpenFOAM, for a range of capillary and Reynolds numbers Ca = 0 . 002 −0 . 5 and Re = 1 −20 0 0 , and very long bubbles, up to 20 times the hydraulic diameter of the channel. The lubrication film surrounding the bubbles is always re- solved by the computational mesh, and therefore the results are representative of a fully-wetting liquid. This study shows that when Ca ≥0.05, the long gas bubble exhibits an axisymmetric shape on the chan- nel cross-section, whereas for lower capillary numbers the bubble flattens at the centre of the channel wall and thick liquid lobes are left at the corners. When Ca ≤0.01, the thin film at the centre of the wall assumes a saddle-like shape, which leads to the formation of two constrictions at the sides of the liquid film profile, where minimum cross-sectional values of the film thickness are observed. The result- ing cross-stream capillary pressure gradients drain liquid out of the thin-film, whose thickness decreases indefinitely as a power-law of the distance from the bubble nose. Therefore, the film thickness depends on the length of the bubble, unlike flow in circular channels. We report detailed values of the centre- line, diagonal and minimum film thickness along the bubble, bubble speed, and cross-sectional gas area fraction, at varying Ca and Re. Inertial effects retard the formation of the saddle-shaped thin-film at the channel centre, which may never form if the bubble is not sufficiently long. However, the film thins at a faster rate towards the bubble rear as the Reynolds number of the flow is increased

Numerical study of the impact of the channel shape on microchannel boiling heat transfer

Flow boiling in multi-microchannel evaporators is recognised as one of the most efficient cooling solutions for high-performance electronics, and has therefore received increasing attention during the recent years. Despite the extensive literature, there is no general agreement yet about the effect of the channel cross-sectional shape on the boiling heat transfer performance, which results on a limited availability of thermal design guidelines and tools. This article presents the results of a systematic analysis of the impact of the channel shape on the bubble dynamics and heat transfer, under flow boiling conditions. Simulations are carried out using a customised version of OpenFOAM, and the Volume-Of-Fluid method is chosen to capture the liquid-vapour interface dynamics. A benchmark flow model is utilised, where a single isolated bubble is seeded at the channel upstream and transported by a liquid flow across the diabatic section, which is heated by a constant and uniform heat flux. Flow conditions that apply well to the flow boiling of water or refrigerant fluids in sub-millimetre channels at low heat flux ( ~ 10 kW/m2) are investigated, with cross-section width-to-height aspect-ratios ranging from 1 to 8, while the hydraulic diameter of the channel is fixed. This study emphasises that the heat transfer performances for different channel shapes are closely related to the perimetral distribution of the liquid film surrounding the very long bubbles. Square channels exhibit the highest heat transfer coefficients at low flow rates, due to a very thin liquid film that forms at the centre of the wall, but are more at risk of film dryout. High aspect-ratio rectangular channels may be beneficial at larger flow rates, as they promote the formation of an extended liquid film that covers up to 80 % of the cross-section perimeter. At larger aspect-ratios, the average heat transfer coefficient along the shorter wall becomes orders of magnitude smaller than the value detected along the longer wall, owing to a strong asymmetry in the liquid film perimetral distribution

Undulations on the surface of elongated bubbles in confined gas-liquid flows

© 2017 American Physical Society. A systematic analysis is presented of the undulations appearing on the surface of long bubbles in confined gas-liquid flows. CFD simulations of the flow are performed with a self-improved version of the open-source solver ESI OpenFOAM (release 2.3.1), for Ca=0.002-0.1 and Re=0.1-1000, where Ca=μU/σ and Re=2ρUR/μ, with μ and ρ being, respectively, the viscosity and density of the liquid, σ the surface tension, U the bubble velocity, and R the tube radius. A model, based on an extension of the classical axisymmetric Bretherton theory, accounting for inertia and for the curvature of the tube's wall, is adopted to better understand the CFD results. The thickness of the liquid film, and the wavelength and decay rate of the undulations extracted from the CFD simulations, agree well with those obtained with the theoretical model. Inertial effects appear when the Weber number of the flow We=CaRe=O(10-1) and are manifest by a larger number of undulation crests that become evident on the surface of the rear meniscus of the bubble. This study demonstrates that the necessary bubble length for a flat liquid film region to exist between the rear and front menisci rapidly increases above 10R when Ca>0.01 and the value of the Reynolds number approaches 1000

Il-6 Serum Levels and Production Is Related to an Altered Immune Response in Polycystic Ovary Syndrome Girls with Insulin Resistance

Polycystic ovarian syndrome (PCOS) is frequently characterized by obesity and metabolic diseases including hypertension, insulin resistance, and diabetes in adulthood, all leading to an increased risk of atherosclerosis. The present study aimed to evaluate serum and production of inflammatory markers in adolescent Sardinian PCOS. On the basis of HOMA findings, patients were divided into noninsulin resistant (NIR) and insulin resistant (IR), and were weight- and age-matched with healthy girls. Inflammatory cytokines (TNF-α, IL-6, Il-10, TGF-β) and lipokines (leptin, adiponectin), the reactant hs-CRP, and in vitro inflammatory lympho-monocyte response to microbial stimulus were evaluated. In healthy and PCOS subjects, leptin and hs-CRP were correlated with BMI, whereas adiponectin was significantly reduced in all PCOS groups. Although cytokines were similar in all groups, Interleukin-6 (IL-6) was significantly higher in IR PCOS. Moreover, in the latter group lipopolysaccharide-activated monocytes secreted significantly higher levels of IL-6 compared to NIR and control subjects. To conclude, IR PCOS displayed increased IL-6 serum levels and higher secretion in LPS-activated monocytes, whilst revealing no differences for other inflammatory cytokines. These results suggest that in PCOS patients an altered immune response to inflammatory stimuli is present in IR, likely contributing towards determining onset of a low grade inflammation

Dynamics of long gas bubbles rising in a vertical tube in a cocurrent liquid flow

© 2019 American Physical Society. When a confined long gas bubble rises in a vertical tube in a cocurrent liquid flow, its translational velocity is the result of both buoyancy and mean motion of the liquid. A thin film of liquid is formed on the tube wall and its thickness is determined by the interplay of viscous, inertial, capillary and buoyancy effects, as defined by the values of the Bond number (Bo≡ρgR2/σ with ρ being the liquid density, g the gravitational acceleration, R the tube radius, and σ the surface tension), capillary number (Cab≡μUb/σ with Ub being the bubble velocity and μ the liquid dynamic viscosity), and Reynolds number (Reb≡2ρUbR/μ). We perform experiments and numerical simulations to investigate systematically the effect of buoyancy (Bo=0-5) on the shape and velocity of the bubble and on the thickness of the liquid film for Cab=10-3-10-1 and Reb=10-2-103. A theoretical model, based on an extension of Bretherton's lubrication theory, is developed and utilized for parametric analyses; its predictions compare well with the experimental and numerical data. This study shows that buoyancy effects on bubbles rising in a cocurrent liquid flow make the liquid film thicker and the bubble rise faster, when compared to the negligible gravity case. In particular, gravitational forces impact considerably the bubble dynamics already when B

Machine-learning blends of geomorphic descriptors: value and limitations for flood hazard assessment across large floodplains

Recent literature shows several examples of simplified approaches that perform flood hazard (FH) assessment and mapping across large geographical areas on the basis of fast-computing geomorphic descriptors. These approaches may consider a single index (univariate) or use a set of indices simultaneously (multivariate). What is the potential and accuracy of multivariate approaches relative to univariate ones? Can we effectively use these methods for extrapolation purposes, i.e., FH assessment outside the region used for setting up the model? Our study addresses these open problems by considering two separate issues: (1) mapping flood-prone areas and (2) predicting the expected water depth for a given inundation scenario. We blend seven geomorphic descriptors through decision tree models trained on target FH maps, referring to a large study area (∼ 105 km2). We discuss the potential of multivariate approaches relative to the performance of a selected univariate model and on the basis of multiple extrapolation experiments, where models are tested outside their training region. Our results show that multivariate approaches may (a) significantly enhance flood-prone area delineation (accuracy: 92%) relative to univariate ones (accuracy: 84%), (b) provide accurate predictions of expected inundation depths (determination coefficient ∼0.7), and (c) produce encouraging results in extrapolation

Dynamics of long bubbles propagating through cylindrical micro-pin fin arrays

The dynamics of two-phase flows confined within complex and non-straight geometries is of interest for a variety of applications such as micro-pin fin evaporators and flow in unsaturated porous media. Despite the propagation of bubbles in straight channels of circular and noncircular cross-sections has been studied extensively, very little is known about the fluid dynamics features of bubbles and liquid films deposited upon the inner walls of complex geometries. In this work, we investigate the dynamics of long gas bubbles and thin films as bubbles propagate through arrays of in-line cylindrical pins of circular shape in cross-flow, for a range of capillary and Reynolds numbers relevant to heat transfer applications and flow in porous media, different pitch of the cylinders and bubble lengths. Three-dimensional numerical simulations of the two-phase flow are performed using the open-source finite-volume library OpenFOAM v.1812, using a geometric Volume of Fluid (VOF) method to capture the interface dynamics. Systematic analyses are conducted for a range of capillary numbers Ca = 0. 04 --2R, Reynolds numbers Re = 1 -- 1000, streamwise pitch of the cylinders sx = 0. 125R, with R being the radius of the pin fins, and initial bubble length Lb = 2. 5R --12R. The simulations reveal that when bubbles propagate through pin fin arrays, they tend to partially coat the cylinders with a thin liquid film and to expand in the cross-stream direction within the gap left between adjacent cylinders. The liquid film deposited on the cylinders is significantly thinner than that reported for straight channels and similar geometrical constraints. As the streamwise distance between the cylinders is decreased, the flow configuration tends towards that for a straight channel, whereas larger distances cause the bubble to expand excessively in the cross-stream direction, and to eventually arrest when sx > 2R. Inertial effects have a strong impact on the bubble shape and dynamics when Re > 500, triggering time-dependent patterns that lead to bubble fragmentation and much thicker liquid films

Coupled atomistic–continuum simulations of nucleate boiling

Data availability: The software used to generate the data is publicly available on github and linked to the submission.Boiling is a striking example of a multiscale process, where the dynamics of bubbles is governed by the interplay between the molecular interactions responsible for nucleation, and the macroscale hydrodynamic and thermal boundary layers. A complete description of this phenomenon requires coupling molecular- and continuum-scale fluid mechanics into a single modelling framework. This article presents a hybrid atomistic–continuum computational model for coupled simulations of nucleate boiling. A domain decomposition coupling method is utilised, where the near-wall region is solved by a Molecular Dynamics description, which handles nucleation and the moving contact lines, while the bulk flow region is solved by a continuum-scale description based on the Navier–Stokes equations. The latter employs a Volume Of Fluid method to track the evolution of the liquid–vapour interface and the interphase mass transfer is computed via the Hertz–Knudsen–Schrage relationship. Boiling of a Lennard-Jones fluid over a heated wall is simulated and the hybrid solution is validated against a fully molecular solution. The results obtained with the coupled framework in terms of time-dependent bubble volume, phase-change rates, bubble dynamics and evolution of the temperature field agree quantitatively with those achieved by a MD-only simulation. The coupled framework reproduces the bubble growth rate over time from nucleation until a bubble diameter of about 70 nm, demonstrating the accuracy and robustness of the coupling architecture. This also demonstrates that the fluid dynamics description based on the Navier–Stokes equations is capable of correctly capturing the main heat and mass transfer mechanisms responsible for bubble growth at the nanoscale. The proposed modelling framework paves the way towards multiscale simulations of boiling, where the necessary molecular-level physics is retained in a computational fluid dynamics solver.This work was funded under the embedded CSE programme of the ARCHER2 UK National Supercomputing Service (https://www.archer2.ac.uk), project ARCHER2-eCSE06-1 “Hybrid Atomistic–Continuum Simulations of Boiling Across Scales”

Computational study of bubble, thin-film dynamics and heat transfer during flow boiling in non-circular microchannels

Flow boiling in multi-microchannel evaporators is one of the most efficient thermal management solutions for high-power-density applications. However, there is still a lack of understanding of the governing two-phase heat and mass transfer processes that occur in these devices, which has resulted in a limited availability of applicable boiling heat transfer prediction methods based on first principles, and of reliable thermal design tools. This article presents a systematic analysis of the dynamics of bubbles and the surrounding liquid film during flow boiling in three-side-heated non-circular microchannels. The study is performed using a custom version of ESI OpenFOAM v2106 with a geometric volume-of-fluid method to capture the interface dynamics, also incorporating conjugate heat transfer through the evaporator walls. The hydraulic diameter of the channel is fixed to ℎ = 0.229 mm and the range of width-to-height aspect ratios = 0.25−4 is examined. We investigate different fluids, namely water, HFE7100, R1233zd(E), R1234ze(E), and evaporator materials, namely copper, aluminium, silicon, stainless steel, with base heat fluxes in the range = 50 − 200 kW∕m2. The results show that conjugate heat transfer acts to make the temperature distributions around the perimeter of the channel cross-section more uniform, and that the topography of the lubricating film and the extension of the dry vapour patches that develop while the film is depleted both depend on the cross-sectional channel shape and influence the heat transfer performance significantly. For highly wetting conditions, channels with = 0.25 tend to allow enhanced heat transfer rates, with a spatially-averaged Nusselt number that is 50% higher than that obtained for = 1 (square channels) and 10% higher than that for = 4. This arises thanks to an extended evaporating film that covers the vertical walls which, owing to the three-side-heated configuration, contribute twice to the spatially-averaged heat transfer performance. For more hydrophobic conditions, large dry patches develop over the vertical walls for = 0.25 due to the lower evaporator temperatures, leading to reduced heat transfer, with thermal performance weakly dependent on in the range = 0.5 − 2