Researcher profile

Marcello Sega

Marcello Sega contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
8topics
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

9 published item(s)

preprint2022arXiv

On the accuracy and performance of the lattice Boltzmann method with 64-bit, 32-bit and novel 16-bit number formats

Fluid dynamics simulations with the lattice Boltzmann method (LBM) are very memory-intensive. Alongside reduction in memory footprint, significant performance benefits can be achieved by using FP32 (single) precision compared to FP64 (double) precision, especially on GPUs. Here, we evaluate the possibility to use even FP16 and Posit16 (half) precision for storing fluid populations, while still carrying arithmetic operations in FP32. For this, we first show that the commonly occurring number range in the LBM is a lot smaller than the FP16 number range. Based on this observation, we develop novel 16-bit formats - based on a modified IEEE-754 and on a modified Posit standard - that are specifically tailored to the needs of the LBM. We then carry out an in-depth characterization of LBM accuracy for six different test systems with increasing complexity: Poiseuille flow, Taylor-Green vortices, Karman vortex streets, lid-driven cavity, a microcapsule in shear flow (utilizing the immersed-boundary method) and finally the impact of a raindrop (based on a Volume-of-Fluid approach). We find that the difference in accuracy between FP64 and FP32 is negligible in almost all cases, and that for a large number of cases even 16-bit is sufficient. Finally, we provide a detailed performance analysis of all precision levels on a large number of hardware microarchitectures and show that significant speedup is achieved with mixed FP32/16-bit.

preprint2022arXiv

Structure and scaling laws of liquid/vapor interfaces close to the critical point

To reach a deeper understanding of fluid interfaces it is necessary to identify a meaningful coarse-graining length that separates intrinsic fluctuations from capillary ones, given the lack of a proper statistical mechanical definition of the latter. Here, with the help of unsupervised learning techniques, we introduce a new length scale based on the local density of the fluid. This length scale follows a scaling law that diverges more mildly than the bulk correlation length upon approaching the critical point. This allows to distinguish regimes of correlated and uncorrelated capillary waves from that of intrinsic fluctuations.

preprint2020arXiv

Capillary interactions between soft capsules protruding through thin fluid films

When a suspension dries, the suspending fluid evaporates, leaving behind a dry film composed of the suspended particles. During the final stages of drying, the height of the fluid film on the substrate drops below the particle size, inducing local interface deformations that lead to strong capillary interactions among the particles. Although capillary interactions between rigid particles are well studied, much is still to be understood about the behaviour of soft particles and the role of their softness during the final stages of film drying. Here, we use our recently-introduced numerical method that couples a fluid described using the lattice Boltzmann approach to a finite element description of deformable objects to investigate the drying process of a film with suspended soft particles. Our measured menisci deformations and lateral capillary forces, which agree well with previous theoretical and experimental works in case of rigid particles, show the deformations become smaller with increasing particles softness, resulting in weaker lateral interaction forces. At large interparticle distances, the force approaches that of rigid particles. Finally, we investigate the time dependent formation of particle clusters at the late stages of the film drying.

preprint2014arXiv

Mesoscale Structures at Complex Fluid-Fluid Interfaces: a Novel Lattice Boltzmann / Molecular Dynamics Coupling

Complex fluid-fluid interfaces featuring mesoscale structures with adsorbed particles are key components of newly designed materials which are continuously enriching the field of soft matter. Simulation tools which are able to cope with the different scales characterizing these systems are fundamental requirements for efficient theoretical investigations. In this paper we present a novel simulation method, based on the approach of Ahlrichs and Dünweg [Ahlrichs and Dünweg, Int. J. Mod. Phys. C, 1998, 9, 1429], that couples the "Shan-Chen" multicomponent Lattice Boltzmann technique to off-lattice molecular dynamics to simulate efficiently complex fluid-fluid interfaces. We demonstrate how this approach can be used to study a wide class of challenging problems. Several examples are given, with an accent on bicontinuous phases formation in polyelectrolyte solutions and ferrofluid emulsions. We also show that the introduction of solvation free energies in the particle-fluid interaction unveils the hidden, multiscale nature of the particle-fluid coupling, allowing to treat symmetrically (and interchangeably) the on-lattice and off-lattice components of the system.

preprint2014arXiv

Regularization of the slip length divergence in water nanoflows by inhomogeneities at the Angstrom scale

We performed non-equilibrium Molecular Dynamics simulations of water flow in nano-channels with the aim of discriminating {\it static} from {\it dynamic} contributions of the solid surface to the slip length of the molecular flow. We show that the regularization of the slip length divergence at high shear rates, formerly attributed to the wall dynamics, is controlled instead by its static properties. Surprisingly, we find that atomic displacements at the Angstrom scale are sufficient to remove the divergence of the slip length and realize the no-slip condition. Since surface thermal fluctuations at room temperature are enough to generate these displacements, we argue that the no-slip condition for water can be achieved also for ideal surfaces, which do not present any surface roughness.

preprint2012arXiv

Identifying interfacial molecules in nonplanar interfaces: the generalized ITIM algorithm

We present a generalized version of the ITIM algorithm for the identification of interfacial molecules, which is able to treat arbitrarily shaped interfaces. The algorithm exploits the similarities between the concept of probe sphere used in ITIM and the circumsphere criterion used in the alpha-shapes approach, and can be regarded either as a reference-frame independent version of the former, or as an extended version of the latter that includes the atomic excluded volume. The new algorithm is applied to compute the intrinsic orientational order parameters of water around a DPC and a cholic acid micelle in aqueous environment, and to the identification of solvent-reachable sites in four model structures for soot. The additional algorithm introduced for the calculation of intrinsic density profiles in arbitrary geometries proved to be extremely useful also for planar interfaces, as it allows to solve the paradox of smeared intrinsic profiles far from the interface

preprint2010arXiv

Influence of pore dielectric boundaries on the translocation barrier of DNA

We investigate the impact of dielectric boundary forces on the translocation process of charged rigid DNA segments through solid neutral nanopores. We assess the electrostatic contribution to the translocation free energy barrier of a model DNA segment by evaluating the potential of mean force in absence and presence of polarization effects by means of coarse-grained molecular dynamics simulations. The effect of induced polarization charges has been taken into account by employing ICC*, a recently developed algorithm that can efficiently compute induced polarization charges induced on suitably discretized dielectric boundaries. Since water has a higher dielectric constant than the pore walls, polarization effects repel charged objects in the vicinity of the interface, with the effect of significantly increasing the free energy barrier. Another investigated side effect is the change of the counterion distribution around the charged polymer in presence of the induced pore charges. Furthermore we investigate the influence of adding salt to the solution.

preprint2010arXiv

The ICC* Algorithm: A fast way to include dielectric boundary effects into molecular dynamics simulations

We employ a fast and accurate algorithm to treat dielectric interfaces within molecular dynamics simulations and demonstrate the importance of dielectric boundary forces (DBFs) in two systems of interests in soft-condensed matter science. We investigate a salt solution confined to a slit pore, and a model of a DNA fragment translocating thorugh a narrow pore.

preprint2009arXiv

Mesoscopic simulations of the counterion-induced electroosmotic flow - a comparative study

We present mesoscopic simulations of the counterion-induced electroosmotic flow in different electrostatic coupling regimes. Two simulation methods are compared, Dissipative Particle Dynamics (DPD) and coupled Lattice-Boltzmann/Molecular Dynamics (LB/MD). A general mapping scheme to match DPD to LB/MD is developed. For the weak-coupling regime, analytic expressions for the flow profiles in the presence of partial-slip as well as no-slip boundary conditions are derived from the Poisson-Boltzmann and Stokes equations, which are in good agreement with the numerical results. The influence of electrofriction and partial slip on the flow profiles is discussed.