BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Luca Brandt

Luca Brandt contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.flu-dyncond-mat.softphysics.comp-phComputational Engineering, Finance, and Science

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
20works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2023arXiv

Intermittency in turbulent emulsions

We investigate the statistics of turbulence in emulsions of two-immiscible fluids of same density. We compute for the first time velocity increments between points conditioned to be located in the same phase or in different phases and examine their probability density functions (PDF) and the associated structure functions (SF). This enables us to demonstrate that the the presence of the interface reduces the skewness of the PDF at scales below the Kolmogorov-Hinze scale and therefore the magnitude of the energy flux towards the dissipative scales, which is quantified by the third-order SF. The analysis of the higher order SFs shows that multiphase turbulence is more intermittent than single-phase turbulence. In particular, the local scaling exponents of the SFs display a saturation about the Kolmogorov-Hinze scale and below, which indicates the presence of large velocity gradients across the interface. Interestingly, the statistics approach of classic homogeneous isotropic turbulence when significantly increasing the viscosity of the dispersed phase.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

FluTAS: A GPU-accelerated finite difference code for multiphase flows

We present the Fluid Transport Accelerated Solver, FluTAS, a scalable GPU code for multiphase flows with thermal effects. The code solves the incompressible Navier-Stokes equation for two-fluid systems, with a direct FFT-based Poisson solver for the pressure equation. The interface between the two fluids is represented with the Volume of Fluid (VoF) method, which is mass conserving and well suited for complex flows thanks to its capacity of handling topological changes. The energy equation is explicitly solved and coupled with the momentum equation through the Boussinesq approximation. The code is conceived in a modular fashion so that different numerical methods can be used independently, the existing routines can be modified, and new ones can be included in a straightforward and sustainable manner. FluTAS is written in modern Fortran and parallelized using hybrid MPI/OpenMP in the CPU-only version and accelerated with OpenACC directives in the GPU implementation. We present different benchmarks to validate the code, and two large-scale simulations of fundamental interest in turbulent multiphase flows: isothermal emulsions in HIT and two-layer Rayleigh-Bénard convection. FluTAS is distributed through a MIT license and arises from a collaborative effort of several scientists, aiming to become a flexible tool to study complex multiphase flows.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The interaction of droplet dynamics and turbulence cascade

The dynamics of droplet fragmentation in turbulence is described in the Kolmogorov-Hinze framework. Yet, a quantitative theory is lacking at higher concentrations when strong interactions between the phases and coalescence become relevant, which is common in most flows. Here, we address this issue through a fully-coupled numerical study of the droplet dynamics in a turbulent flow at high Reynolds number. By means of time-space spectral statistics, not currently accessible to experiments, we demonstrate that the characteristic scale of the process, the Hinze scale, can be precisely identified as the scale at which the net energy exchange due to capillarity is zero. Droplets larger than this scale preferentially break up absorbing energy from the flow; smaller droplets, instead, undergo rapid oscillations and tend to coalesce releasing energy to the flow. Further, we link the droplet-size-distribution with the probability distribution of the turbulent dissipation. This shows that key in the fragmentation process is the local flux of energy which dominates the process at large scales, vindicating its locality.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Modulation of homogeneous and isotropic turbulence in emulsions

We present a numerical study of emulsions in homogeneous and isotropic turbulence at $Re_λ=137$. The problem is addressed via Direct Numerical Simulations (DNS), where the Volume of Fluid (VOF) is used to represent the complex features of the liquid-liquid interface. We consider a mixture of two iso-density fluids, where fluid properties are varied with the goal of understanding their role in turbulence modulation, in particular the volume fraction ($0.03<α<0.5$), viscosity ratio ($0.01<μ_d/μ_c<100$) and large scale Weber number ($10.6<We_\mathcal{L}<106.5$). The analysis, performed by studying integral quantities and spectral scale-by-scale analysis, reveals that energy is consistently transported from large to small scales by the interface, and no inverse cascade is observed. Furthermore, the total surface is found to be directly proportional to the amount of energy transported, while viscosity and surface tension alter the dynamic that regulates energy transport. We also observe the $-10/3$ and $-3/2$ scaling on droplet size distributions, suggesting that the dimensional arguments which led to their derivation are verified in HIT conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A volume-of-fluid method for interface-resolved simulations of phase-changing two-fluid flows

We present a numerical method for interface-resolved simulations of evaporating two-fluid flows based on the volume-of-fluid (VoF) method. The method has been implemented in an efficient FFT-based two-fluid Navier-Stokes solver, using an algebraic VoF method for the interface representation, and extended with the transport equations of thermal energy and vaporized liquid mass for the single-component evaporating liquid in an inert gas. The conservation of vaporizing liquid and computation of the interfacial mass flux are performed with the aid of a reconstructed signed-distance field, which enables the use of well-established methods for phase change solvers based on level-set methods. The interface velocity is computed with a novel approach that ensures accurate mass conservation, by constructing a divergence-free extension of the liquid velocity field onto the entire domain. The resulting approach does not depend on the type of interface reconstruction (i.e.\ can be employed in both algebraic and geometrical VoF methods). We extensively verified and validated the overall method against several benchmark cases, and demonstrated its excellent mass conservation and good overall performance for simulating evaporating two-fluid flows in two and three dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Aerodynamically-driven rupture of a liquid film by turbulent shear flow

The rupture of a liquid film due to co-flowing turbulent shear flows in the gas phase is studied using a volume-of-fluid method. To simulate this multiphase problem, we use a simplified numerical setup where the liquid film is &#39;sandwiched&#39; between two fully developed boundary layers from a turbulent channel simulation. The film deforms and eventually ruptures within the shear zone created by the co-flows. This efficient setup allows systematic variation of physical parameters to gauge their role in the aerodynamically-driven deformation and rupture of a liquid film under fully developed sheared turbulence. The present work presents a detailed study of the developing pressure field over the deforming film and related aerodynamic effects, as previously suggested by other authors, in particular the role of the inviscid lift and drag forces. A cumulative lift force is introduced to capture the effect of the alternating pressure minima and maxima forming over the film which amplify and eventually rupture the film. A velocity scale derived from the lift-induced drag force reflects the state of the turbulent boundary layer over the film and collapses the temporal development of this cumulative lift force as well as the amplitude of film deformation with some success for the different film thicknesses and Reynolds numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An interface capturing method for liquid-gas flows at low-Mach number

Multiphase, compressible and viscous flows are of crucial importance in a wide range of scientific and engineering problems. Despite the large effort paid in the last decades to develop accurate and efficient numerical techniques to address this kind of problems, current models need to be further improved to address realistic applications. In this context, we propose a numerical approach to the simulation of multiphase, viscous flows where a compressible and an incompressible phase interact in the low-Mach number regime. In this frame, acoustics is neglected but large density variations of the compressible phase can be accounted for as well as heat transfer, convection and diffusion processes. The problem is addressed in a fully Eulerian framework exploiting a low-Mach number asymptotic expansion of the Navier-Stokes equations. A Volume of Fluid approach (VOF) is used to capture the liquid-gas interface, built on top of a massive parallel solver, second order accurate both in time and space. The second-order-pressure term is treated implicitly and the resulting pressure equation is solved with the eigenexpansion method employing a robust and novel formulation. We provide a detailed and complete description of the theoretical approach together with information about the numerical technique and implementation details. Results of benchmarking tests are provided for five different test cases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effects of surface nanostructure and wettability on pool boiling: A molecular dynamics study

We study the role of surface topology, surface chemistry, and wall superheat temperature on the onset of boiling, bubble nucleation and growth, and the possible formation of an insulating vapour film by means of large-scale MD simulations. In the numerical experiments, we control the system pressure by imposing a constant force on a moving piston. The simulations reveal that the presence of a nanostructure triggers the bubble formation, determines the nucleation site and facilitates the energy transfer from the hot substrate to the water. The surface chemistry, on the other hand, governs the shape of the formed bubble. A hydrophilic surface chemistry accelerates the bubble nucleation, however, decelerates the bubble expansion, thus postponing the formation of the film of vapour. Therefore, a hydrophilic surface provides better energy transfer from the hot wall to the water. By analysing the system energy, we show that irrespective of wall topology and chemistry, there is a wall temperature for which the amount of transferred energy is maximum.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

GPU acceleration of CaNS for massively-parallel direct numerical simulations of canonical fluid flows

This work presents the GPU acceleration of the open-source code CaNS for very fast massively-parallel simulations of canonical fluid flows. The distinct feature of the many-CPU Navier-Stokes solver in CaNS is its fast direct solver for the second-order finite-difference Poisson equation, based on the method of eigenfunction expansions. The solver implements all the boundary conditions valid for this type of problems in a unified framework. Here, we extend the solver for GPU-accelerated clusters using CUDA Fortran. The porting makes extensive use of CUF kernels and has been greatly simplified by the unified memory feature of CUDA Fortran, which handles the data migration between host (CPU) and device (GPU) without defining new arrays in the source code. The overall implementation has been validated against benchmark data for turbulent channel flow and its performance assessed on a NVIDIA DGX-2 system (16 Tesla V100 32Gb, connected with NVLink via NVSwitch). The wall-clock time per time step of the GPU-accelerated implementation is impressively small when compared to its CPU implementation on state-of-the-art many-CPU clusters, as long as the domain partitioning is sufficiently small that the data resides mostly on the GPUs. The implementation has been made freely available and open-source under the terms of an MIT license.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Implementation note on a minimal hybrid lubrication/granular dynamics model for dense suspensions

We describe and summarize a class of minimal numerical models emerged from recent development of simulation methods for dense particle suspensions in overdamped linear flows. The main ingredients include (i) a frame-invariant, short-range lubrication model for spherical particles, and (ii) a soft-core, stick/slide frictional contact model activated when particles overlap. We implement a version of the model using a modified velocity-Verlet algorithm that explicitly solves the $N$-body dynamical system in $\mathcal{O}(cN)$ operations, where $c$ is a kernel constant depending on the cutoff of particle interactions. The implementation is validated against literature results on jamming transition and shear thickening suspensions from 40% to 64% volume fractions. Potential strategies to extend the present methodology to non-spherical particles are also suggested for very concentrated suspensions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Low Reynolds number turbulent flows over elastic walls

We study the laminar and turbulent channel flow over a viscous hyper-elastic wall and show that it is possible to sustain an unsteady chaotic turbulent-like flow at any Reynolds number by properly choosing the wall elastic modulus. We propose a physical explanation for this effect by evaluating the shear stress and the turbulent kinetic energy budget in the fluid and elastic layer. We vary the bulk Reynolds number from 2800 to 10 and identify two distinct mechanisms for turbulence production. At moderate and high Reynolds numbers, turbulent fluctuations activate the wall oscillations, which, in turn, amplify the turbulent Reynolds stresses in the fluid. At very low Reynolds number, the only production term is due to the energy input from the elastic wall, which increases with the wall elasticity. This mechanism may be exploited to passively enhance mixing in microfluidic devices.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical simulations of a sphere settling in simple shear flows of yield stress fluids

We perform $3$D numerical simulations to investigate the sedimentation of a single sphere in the absence and presence of a simple cross shear flow in a yield stress fluid with weak inertia. In our simulations, the settling flow is considered to be the primary flow, whereas the linear cross shear flow is a secondary flow with amplitude $10\%$ of the primary flow. To study the effects of elasticity and plasticity of the carrying fluid on the sphere drag as well as the flow dynamics, the fluid is modeled using the elastovisco-plastic (EVP) constitutive laws proposed by \cite{saramito2009new}. The extra non-Newtonian stress tensor is fully coupled with the flow equation and the solid particle is represented by an immersed boundary (IB) method. Our results show that the fore-aft asymmetry in the velocity is less pronounced and the negative wake disappears when a linear cross shear flow is applied. We find that the drag on a sphere settling in a sheared yield stress fluid is reduced significantly as compared to an otherwise quiescent fluid. More importantly, the sphere drag in the presence of a secondary cross shear flow cannot be derived from the pure sedimentation drag law owing to the non-linear coupling between the simple shear flow and the uniform flow. Finally, we show that the drag on the sphere settling in a sheared yield-stress fluid is reduced at higher material elasticity mainly due to the form and viscous drag reduction.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical Simulations of Vorticity Banding of Emulsions in Shear Flows

Multiphase shear flows often show banded structures that affect the global behavior of complex fluids e.g. in microdevices. Here we investigate numerically the banding of emulsions, i.e. the formation of regions of high and low volume fraction, alternated in the vorticity direction and aligned with the flow (shear bands). These bands are associated with a decrease of the effective viscosity of the system. To understand the mechanism of banding experimentally observed we have performed interface resolved simulations of the two-fluid system. The experiments were perfomed starting with a random distribution of droplets which, under the applied shear, evolves in time resulting in a phase separation. To numerically reproduce this process, the banded structures are initialized in a narrow channel confined by two walls moving in opposite direction. We find that the initial banded distribution is stable when droplets are free to merge and unstable when coalescence is prevented. In this case, additionally, the effective viscosity of the system increases, resembling the rheological behavior of suspensions of deformable particles. Droplets coalescence, on the other hand, allows emulsions to reduce the total surface of the system and hence the energy dissipation associated to the deformation, which in turn reduces the effective viscosity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Regimes of heat transfer in particle suspensions

We present results of interface-resolved simulations of heat transfer in suspensions of finite-size neutrally-buoyant spherical particles for solid volume fractions up to 35% and bulk Reynolds numbers from 500 to 5600. An Immersed Boundary-Volume of Fluid method is used to solve the energy equation in the fluid and solid phase. We relate the heat transfer to the regimes of particle motion previously identified, i.e. a viscous regime at low volume fractions and low Reynolds number, particle-laden turbulence at high Reynolds and moderate volume fraction and particulate regime at high volume fractions. We show that in the viscous dominated regime, the heat transfer is mainly due to thermal diffusion with enhancement due to the particle-induced fluctuations. In the turbulent-like regime, we observe the largest enhancement of the global heat transfer, dominated by the turbulent heat flux. In the particulate shear-thickening regime, however, the heat transfer enhancement decreases as mixing is quenched by the particle migration towards the channel core. As a result, a compact loosely-packed core region forms and the contribution of thermal diffusion to the total heat transfer becomes significant once again. The global heat transfer becomes, in these flows at volume fractions larger than 25%, lower than in single-phase turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Weakly adhesive suspension shows rate-dependence in oscillatory but not steady shear flows

We report rheological measurements of a noncolloidal particle suspension in a Newtonian solvent at 40% solid volume fraction. An anomalous, frequency-dependent complex viscosity is found under oscillatory shear (OS) flow, whereas a constant dynamic viscosity is found under the same shear rates in steady shear (SS) flow. We show that this contradiction arises from the underlying microstructural difference between OS and SS, mediated by weak interparticle forces. Discrete element simulations of proxy particle suspensions confirm this hypothesis and reveal an adhesion-induced, shear thinning mechanism with a -1/5 slope, only in OS, in agreement with experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Finite-size spherical particles in a square duct flow of an elastoviscoplastic fluid: an experimental study

The present experimental study addresses the flow of a Yield Stress Fluid with some elasticity (Carbopol gel) in a square duct. The behaviour of two fluids with lower and higher yield stress is investigated at multiple Reynolds numbers $Re^* \in$ (1, 200) and Bingham numbers $Bi \in$ (0.01, 0.35). A secondary flow consisting of eight vortices is observed to recirculate the fluid from the corners to the core. Their extent and intensity grows with increasing $Re^*$. The second objective of this study is to explore the change in flow in the presence of particles. Almost neutrally-buoyant finite-size spherical particles with duct height, $2H$, to particle diameter, $d_p$, ratio of 12 are used at two volume fractions $ϕ$ = 5 and 10\%. Particle Tracking Velocimetry (PTV) is used to measure the velocity of these refractive-index-matched spheres, and PIV to extract the fluid velocity. Simple shadowgraphy is also used for qualitatively visualising the development of the particle distribution along the streamwise direction. The particle distribution pattern changes from being concentrated at the four corners, at low flow rates, to being focussed along a diffused ring between the center and the corners, at high flow rates. The presence of particles induces streamwise and wall-normal velocity fluctuations in the fluid phase; however, the primary Reynolds shear stress is still very small compared to turbulent flows. The size of the plug in the particle-laden cases appears to be smaller than the corresponding single phase cases. Similar to Newtonian fluids, the friction factor increases due to the presence of particles, almost independently of the suspending fluid matrix. Interestingly, predictions based on an increased effective suspension viscosity agrees quite well with the experimental friction factor for the concentrations used in this study.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Effects of the finite particle size in turbulent wall-bounded flows of dense suspensions

We use interface-resolved simulations to study finite-size effects in turbulent channel flow of neutrally-buoyant spheres. Two cases with particle sizes differing by a factor of 2, at the same solid volume fraction of 20% and bulk Reynolds number are considered. These are complemented with two reference single-phase flows: the unladen case, and the flow of a Newtonian fluid with the effective suspension viscosity of the same mixture in the laminar regime. As recently highlighted in Costa et al. (PRL 117, 134501), a particle-wall layer is responsible for deviations of the statistics from what is observed in the continuum limit where the suspension is modeled as a Newtonian fluid with an effective viscosity. Here we investigate the fluid and particle dynamics in this layer and in the bulk. In the particle-wall layer, the near wall inhomogeneity has an influence on the suspension micro-structure over a distance proportional to the particle size. In this layer, particles have a significant (apparent) slip velocity that is reflected in the distribution of wall shear stresses. This is characterized by extreme events (both much higher and much lower than the mean). Based on these observations we provide a scaling for the particle-to-fluid apparent slip velocity as a function of the flow parameters. We also extend the flow scaling laws in to second-order Eulerian statistics in the homogeneous suspension region away from the wall. Finite-size effects in the bulk of the channel become important for larger particles, while negligible for lower-order statistics and smaller particles. Finally, we study the particle dynamics along the wall-normal direction. Our results suggest that 1-point dispersion is dominated by particle-turbulence (and not particle-particle) interactions, while differences in 2-point dispersion and collisional dynamics are consistent with a picture of shear-driven interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Experimental investigation of turbulent suspensions of spherical particles in a square duct

We report experimental observations of turbulent flow with spherical particles in a square duct. Three particle sizes namely: $2H/d_{p}$ = 40, 16 and 9 ($2H$ being the duct full height and $d_{p}$ being the particle diameter) are investigated. The particles are nearly neutrally-buoyant with respect to the suspending fluid. Refractive Index Matched - Particle Image Velocimetry is used for fluid velocity measurement even at the highest particle volume fraction (20%) and Particle Tracking Velocimetry for the particle velocity statistics for the flows seeded with particles of the two largest sizes, whereas only pressure measurements are reported for the smallest particles. Settling effects are seen at the lowest bulk Reynolds number $Re_{2H}\approx$ 10000 whereas, at the highest $Re_{2H}\approx$ 27000, particles are in almost full suspension. The friction factor of the suspensions is found to be significantly larger than that of single-phase duct flow at the lower $Re_{2H}$ investigated; however, the difference decreases when increasing the flow rate and the total drag approaches the values of the single phase flow at the higher Reynolds number considered. The pressure drop is found to decrease with the particle diameter for volume fractions lower than $ϕ$ = 10% for nearly all $Re_{2H}$. However, at the highest volume fraction $ϕ$ = 20%, we report a peculiar non-monotonic behavior: the pressure drop first decreases and then increases with increasing particle size. The decrease of the turbulent drag with particle size at the lowest volume fractions is related to an attenuation of the turbulence. The drag increase for the two largest particle sizes at $ϕ$ = 20%, however, occurs despite this large reduction of the turbulent stresses, and it is therefore due to significant particle-induced stresses.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Turbulent channel flow of a dense binary mixture of rigid particles

We study turbulent channel flow of a binary mixture of finite-size neutrally-buoyant rigid particles by means of interface-resolved direct numerical simulations. We fix the bulk Reynolds number and total solid volume fraction, $Re_b = 5600$ and $Φ=20\%$, and vary the relative fraction of small and large particles. The binary mixture consists of particles of two different sizes, $2h/d_l=20$ and $2h/d_s=30$ where $h$ is the half channel height and $d_l$ and $d_s$ the diameter of the large and small particles. While the particulate flow statistics exhibit a significant alteration of the mean velocity profile and turbulent fluctuations with respect to the unladen flow, the differences between the mono-disperse and bi-disperse cases are small. However, we observe a clear segregation of small particles at the wall in binary mixtures, which affects the dynamics of the near wall region and thus the overall drag. This results in a higher drag in suspensions with a larger amount of large particles. As regards bi-disperse effects on the particle dynamics, a non-monotonic variation of the particle dispersion in the spanwise (homogeneous) direction is observed when increasing the percentage of small/large particles. Finally, we note that particles of the same size tend to cluster more at contact whereas the dynamics of the large particles gives highest collision kernels due to a higher approaching speed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical Study of the Sedimentation of Spheroidal Particles

The gravity-driven motion of rigid particles in a viscous fluid is relevant in many natural and industrial processes, yet this has mainly been investigated for spherical particles. We therefore consider the sedimentation of non-spherical (spheroidal) isolated and particle pairs in a viscous fluid via numerical simulations using the Immersed Boundary Method. The simulations performed here show that the critical Galileo number for the onset of secondary motions decreases as the spheroid aspect ratio departs from 1. Above this critical threshold, oblate particles perform a zigzagging motion whereas prolate particles rotate around the vertical axis while having their broad side facing the falling direction. Instabilities of the vortices in the wake follow when farther increasing the Galileo number. We also study the drafting-kissing-tumbling associated with the settling of particle pairs. We find that the interaction time increases significantly for non-spherical particles and, more interestingly, spheroidal particles are attracted from larger lateral displacements. This has important implications for the estimation of collision kernels and can result in increasing clustering in suspensions of sedimenting spheroids.

0 citations0 reviewsNo code linkedOpen access