Showing 1–13 of 13 results for author: de Tullio, M D

Cross-flow oscillations of a circular cylinder with mechanically coupled rotation

Authors: Alessandro Nitti, Giovanni De Cillis, Marco . D. de Tullio

Abstract: Flow-Induced Vibrations (FIV) of an elastically mounted circular cylinder are investigated by means of two-dimensional simulations. A mechanical coupling between cross-flow translation and rotation provides a single degree-of-freedom system in which the coupled rotational oscillations affect the fluid-structure dynamics. The structural response of this system is investigated exploring the design s… ▽ More Flow-Induced Vibrations (FIV) of an elastically mounted circular cylinder are investigated by means of two-dimensional simulations. A mechanical coupling between cross-flow translation and rotation provides a single degree-of-freedom system in which the coupled rotational oscillations affect the fluid-structure dynamics. The structural response of this system is investigated exploring the design space spanned by reduced velocity, coupling radius, and phase density ratio. The kinematic coupling introduces the rotation-induced shear layer modifications, as well as an equivalent inertia effect connected to the coupling force. Such a computational campaign is carried out by means of direct numerical simulations with immersed boundary forcing at a Reynolds number equal to 100. The investigated system exhibits the wake-body synchronisation features typical of lock-in for non-rotating cylinders. However, the kinematic coupling provides a novel FIV scenario, in which the oscillation amplitude is magnified in the locked configurations with respect to the forced rotation case. Furthermore, it is found a significant widening of the reduced velocity domain where the lock-in condition takes place. In view of the proposed analyses, it is determined that the coupled rotation guarantees the phase alignment between lift and displacement necessary to sustain the lock-in condition, making the oscillation amplitude grow indefinitely with the reduced velocity. This is inherently achieved due to the rotational shear layer and the added mass contribution, which prevent the exact match between oscillation frequency and system natural frequency in vacuum. The outcomes of this study might potentially lead to innovative water energy harvesters offering larger power outputs and extended optimal operating regions. △ Less

Submitted 10 July, 2024; originally announced July 2024.

An efficient active-stress electromechanical isogeometric shell model for muscular thin film simulations

Authors: Michele Torre, Simone Morganti, Alessandro Nitti, Marco Donato de Tullio, Josef Kiendl, Francesco Silvio Pasqualini, Alessandro Reali

Abstract: We propose an isogeometric approach to model the deformation of active thin films using layered, nonlinear, Kirchhoff Love shells. Isogeometric Collocation and Galerkin formulations are employed to discretize the electrophysiological and mechanical sub-problems, respectively, with the possibility to adopt different element and time-step sizes. Numerical tests illustrate the capabilities of the act… ▽ More We propose an isogeometric approach to model the deformation of active thin films using layered, nonlinear, Kirchhoff Love shells. Isogeometric Collocation and Galerkin formulations are employed to discretize the electrophysiological and mechanical sub-problems, respectively, with the possibility to adopt different element and time-step sizes. Numerical tests illustrate the capabilities of the active stress based approach to effectively simulate the contraction of thin films in both quasi-static and dynamic conditions. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Predicting the vascular adhesion of deformable drug carriers in narrow capillaries traversed by blood cell

Authors: A. Coclite, G. Pascazio, M. D. de Tullio, P. Decuzzi

Abstract: In vascular targeted therapies, blood-borne carriers should realize sustained drug release from the luminal side towards the diseased tissue. In this context, such carriers are required to firmly adhere to the vessel walls for a sufficient period of time while resisting force perturbations induced by the blood flow and circulating cells. Here, a hybrid computational model, combining a Lattice Bolt… ▽ More In vascular targeted therapies, blood-borne carriers should realize sustained drug release from the luminal side towards the diseased tissue. In this context, such carriers are required to firmly adhere to the vessel walls for a sufficient period of time while resisting force perturbations induced by the blood flow and circulating cells. Here, a hybrid computational model, combining a Lattice Boltzmann (LBM) and Immersed Boundary Methods (IBM), is proposed for predicting the strength of adhesion of particles in narrow capillaries (7.5 $μ\mathrm{m})$ traversed by blood cells. While flowing down the capillary, globular and biconcave deformable cells ( $7 μ\mathrm{m}$ ) encounter $2 μ\mathrm{m}$ discoidal particles, adhering to the vessel walls. Particles present aspect ratios ranging from $0.25$ to $1.0$ and a mechanical stiffness varying from rigid $(\mathrm{Ca}=0)$ to soft $\left(\mathrm{Ca}=10^{-3}\right)$. Cell-particle interactions are quantitatively predicted over time via three independent parameters: the cell membrane stretching $δp$; the cell-to-particle distance $r$, and the number of engaged ligand-receptor bonds $N_{\mathrm{L}}$. △ Less

Submitted 12 October, 2022; originally announced October 2022.

Comments: arXiv admin note: text overlap with arXiv:1907.00080

Journal ref: Journal of Fluid and Structures 82 (2018) 638-650

Characterization of micro-Capsules Deformation in Branching Channels

Authors: Alessandro Coclite, Marco D. de Tullio, Giuseppe Pascazio, Tiziano Politi

Abstract: In this paper, the dynamic of inertial capsules into microfluidic bifurcations is studied. The fluid evolution is based on the solution of the BGK -- lattice Boltzmann scheme including a forcing term accounting for immersed geometries. The dynamic-Immersed Boundary forcing strategy is adopted for imposing no-slip boundary conditions on moving deformable or rigid structures, while, on fixed immerse… ▽ More In this paper, the dynamic of inertial capsules into microfluidic bifurcations is studied. The fluid evolution is based on the solution of the BGK -- lattice Boltzmann scheme including a forcing term accounting for immersed geometries. The dynamic-Immersed Boundary forcing strategy is adopted for imposing no-slip boundary conditions on moving deformable or rigid structures, while, on fixed immersed geometries the Bouzidi-Firdaouss-Lallemand second-order bounce back technique is implemented. The proposed computational framework is employed to detail dynamics and deformation of rigid and deformable capsules traveling into a branching duct. This journey is characterized in terms of i) the capsule/bifurcation interaction depending on the sharpness of the branching channels junction; ii) daughter branches aperture angle; iii) occlusion ratio, the ratio between capsule size and main channel diameter; iv) flowing capsules stiffness; v) number of flowing particles. △ Less

Submitted 13 July, 2022; originally announced July 2022.

Journal ref: Applied Mathematics and Computation 434 (2022) 127445

A Lattice Boltzmann dynamic-Immersed Boundary scheme for the transport of deformable inertial capsules in low-Re flows

Authors: Alessandro Coclite, Sergio Ranaldo, Giuseppe Pascazio, Marco D. de Tullio

Abstract: In this work, a dynamic-Immersed--Boundary method combined with a BGK-Lattice--Boltzmann technique is developed and critically discussed. The fluid evolution is obtained on a three-dimensional lattice with 19 reticular velocities (D3Q19 computational molecule) while the immersed body surface is modeled as a collection of Lagrangian points responding to an elastic potential and a bending resistance… ▽ More In this work, a dynamic-Immersed--Boundary method combined with a BGK-Lattice--Boltzmann technique is developed and critically discussed. The fluid evolution is obtained on a three-dimensional lattice with 19 reticular velocities (D3Q19 computational molecule) while the immersed body surface is modeled as a collection of Lagrangian points responding to an elastic potential and a bending resistance. A moving least squares reconstruction is used to accurately interpolate flow quantities and the forcing field needed to enforce the boundary condition on immersed bodies. The proposed model is widely validated against well known benchmark data for rigid and deformable objects. Rigid transport is validated by computing the settling of a sphere under gravity for five different conditions. Then, the tumbling of inertial particles with different shape is considered, recovering the Jefferey orbit for a prolate spheroid. Moreover, the revolution period for an oblate spheroid and for a disk-like particle is obtained as a function of the Reynolds number. The existence of a critical Reynolds number is demonstrated for both cases above which revolution is inhibited. The transport of deformable objects is also considered. The steady deformation of a membrane under shear for three different mechanical stiffness is assessed. Then, the tumbling of a weakly-deformable spheroid under shear is systematically analyzed as a function strain stiffness, bending resistance and membrane mass. △ Less

Submitted 13 July, 2022; originally announced July 2022.

Journal ref: Published in Computer \& Mathematics with applications Volume 80, Issue 12, 15 December 2020, Pages 2860-2876

FSEI-GPU: GPU accelerated simulations of the fluid-structure-electrophysiology interaction in the left heart

Authors: Francesco Viola, Vamsi Spandan, Valentina Meschini, Joshua Romero, Massimiliano Fatica, Marco D. de Tullio, Roberto Verzicco

Abstract: The reliability of cardiovascular computational models depends on the accurate solution of the hemodynamics, the realistic characterization of the hyperelastic and electric properties of the tissues along with the correct description of their interaction. The resulting fluid-structure-electrophysiology interaction (FSEI) thus requires an immense computational power, usually available in large supe… ▽ More The reliability of cardiovascular computational models depends on the accurate solution of the hemodynamics, the realistic characterization of the hyperelastic and electric properties of the tissues along with the correct description of their interaction. The resulting fluid-structure-electrophysiology interaction (FSEI) thus requires an immense computational power, usually available in large supercomputing centers, and requires long time to obtain results even if multi-CPU processors are used (MPI acceleration). In recent years, graphics processing units (GPUs) have emerged as a convenient platform for high performance computing, as they allow for considerable reductions of the time-to-solution. This approach is particularly appealing if the tool has to support medical decisions that require solutions within reduced times and possibly obtained by local computational resources. Accordingly, our multi-physics solver has been ported to GPU architectures using CUDA Fortran to tackle fast and accurate hemodynamics simulations of the human heart without resorting to large-scale supercomputers. This work describes the use of CUDA to accelerate the FSEI on heterogeneous clusters, where both the CPUs and GPUs are used in synergistically with minor modifications of the original source code. The resulting GPU accelerated code solves a single heartbeat within a few hours (from three to ten depending on the grid resolution) running on premises computing facility made of few GPU cards, which can be easily installed in a medical laboratory or in a hospital, thus opening towards a systematic computational fluid dynamics (CFD) aided diagnostic. △ Less

Submitted 4 May, 2021; v1 submitted 28 March, 2021; originally announced March 2021.

Kinematic and dynamic forcing strategies for predicting the transport of inertial capsules via a combined lattice Boltzmann-Immersed Boundary method

Authors: A. Coclite, S. Ranaldo, M. D. De Tullio, P. Decuzzi, G. Pascazio

Abstract: Modeling the transport of deformable capsules under different flow regimens is crucial in a variety of fields, including oil rheology, blood flow and the dispersion of pollutants. The aim of this study is twofold. Firstly, a combined Lattice Boltzmann-Immersed Boundary (LBM-IB) approach is developed for predicting the transport of inertial deformable capsules. A Moving Least Squares (MLS) scheme h… ▽ More Modeling the transport of deformable capsules under different flow regimens is crucial in a variety of fields, including oil rheology, blood flow and the dispersion of pollutants. The aim of this study is twofold. Firstly, a combined Lattice Boltzmann-Immersed Boundary (LBM-IB) approach is developed for predicting the transport of inertial deformable capsules. A Moving Least Squares (MLS) scheme has been implemented to correlate the pressure, velocity and force fields of the fluid domain with the capsule dynamics. This computational strategy has been named LBM Dynamic IB. Secondly, this strategy is directly compared with a more conventional approach, named LBM Kinematic IB, where capsules move with the same velocity of the surrounding fluid. Multiple test cases have been considered for assessing the accuracy and efficiency of the Dynamic over Kinematic IB scheme, including the stretching of circular capsules in shear flow, the transport in a plane Poiseuille flow of circular and biconcave capsules, with and without inertia. By monitoring the capsule geometry over time, the two schemes have been documented to be in excellent agreement, especially for low Capillary numbers (Ca $\leq$ 0.01), in the case of non-inertial capsules. Despite a moderate increase in computational burden, the presented LBM Dynamic IB scheme is the sole capable of predicting the dynamics of both non-inertial and inertial deformable capsules. The proposed approach can be efficiently employed for studying the transport of blood cells, cancer cells and nano/micro capsules within a capillary flow. △ Less

Submitted 8 March, 2019; originally announced March 2019.

Journal ref: Computers and Fluids, 2019, 180, pp. 41-53

Predicting different adhesive regimens of circulating particles at blood capillary walls

Authors: A. Coclite, H. Mollica, S. Ranaldo, G. Pascazio, M. D. De Tullio, P. Decuzzi

Abstract: A fundamental step in the rational design of vascular targeted particles is the firm adhesion at the blood vessel walls. Here, a combined Lattice Boltzmann Immersed Boundary model is presented for predicting the near wall dynamics of circulating particles. A moving least squares algorithm is used to reconstruct the forcing term accounting for the immersed particle, whereas ligand-receptor binding… ▽ More A fundamental step in the rational design of vascular targeted particles is the firm adhesion at the blood vessel walls. Here, a combined Lattice Boltzmann Immersed Boundary model is presented for predicting the near wall dynamics of circulating particles. A moving least squares algorithm is used to reconstruct the forcing term accounting for the immersed particle, whereas ligand-receptor binding at the particle- wall interface is described via forward and reverse probability distributions. First, it is demonstrated that the model predicts with good accuracy the rolling velocity of tumor cells over an endothelial layer in a microfluidic channel. Then, particle-wall interactions are systematically analyzed in terms of particle geometries (circular, elliptical with aspect ratios 2 and 3), surface ligand densities (0.3, 0.5, 0.7 and 0.9), ligand-receptor bond strengths (1 and 2), and Reynolds numbers (Re = 0.01, 0.1 and 1.0). Depending on these conditions, four different particle-wall interaction regimens are identified, namely not adhering, rolling, sliding and firmly adhering particles. The proposed computational strategy can be efficiently used for predicting the near wall dynamics of particles with arbitrary geometries and surface properties and represents a fundamental tool in the rational design of particles for the specific delivery of therapeutic and imaging agents. △ Less

Submitted 4 August, 2018; originally announced August 2018.

Journal ref: Microfluidics and Nanofluidics 21.11 (2017), p. 168. issn: 1613-4990

A fast moving least squares approximation with adaptive Lagrangian mesh refinement for large scale immersed boundary simulations

Authors: Vamsi Spandan, Detlef Lohse, Marco D. de Tullio, Roberto Verzicco

Abstract: In this paper we propose and test the validity of simple and easy-to-implement algorithms within the immersed boundary framework geared towards large scale simulations involving thousands of deformable bodies in highly turbulent flows. First, we introduce a fast moving least squares (fast-MLS) approximation technique with which we speed up the process of building transfer functions during the simu… ▽ More In this paper we propose and test the validity of simple and easy-to-implement algorithms within the immersed boundary framework geared towards large scale simulations involving thousands of deformable bodies in highly turbulent flows. First, we introduce a fast moving least squares (fast-MLS) approximation technique with which we speed up the process of building transfer functions during the simulations which leads to considerable reductions in computational time. We compare the accuracy of the fast-MLS against the exact moving least squares (MLS) for the standard problem of uniform flow over a sphere. In order to overcome the restrictions set by the resolution coupling of the Lagrangian and Eulerian meshes in this particular immersed boundary method, we present an adaptive Lagrangian mesh refinement procedure that is capable of drastically reducing the number of required nodes of the basic Lagrangian mesh when the immersed boundaries can move and deform. Finally, a coarse-grained collision detection algorithm is presented which can detect collision events between several Lagrangian markers residing on separate complex geometries with minimal computational overhead. △ Less

Submitted 3 May, 2018; originally announced May 2018.

An efficient flamelet progress-variable method for modeling non-premixed flames in weak electric fields

Authors: M. Di Renzo, P. De Palma, M. D. de Tullio, G. Pascazio

Abstract: Combustion stabilization and enhancement of the flammability limits are mandatory objectives to improve nowadays combustion chambers. At this purpose, the use of an electric field in the flame region provides a solution which is, at the same time, easy to implement and effective to modify the flame structure. The present work describes an efficient flamelet progress-variable approach developed to… ▽ More Combustion stabilization and enhancement of the flammability limits are mandatory objectives to improve nowadays combustion chambers. At this purpose, the use of an electric field in the flame region provides a solution which is, at the same time, easy to implement and effective to modify the flame structure. The present work describes an efficient flamelet progress-variable approach developed to model the fluid dynamics of flames immersed in an electric field. The main feature of this model is that it can use complex ionization mechanisms without increasing the computational cost of the simulation. The model is based on the assumption that the combustion process is not directly influenced by the electric field and has been tested using two chemi-ionization mechanisms of different complexity in order to examine its behavior with and without the presence of heavy anions in the mixture. Using a one- and two-dimensional numerical test cases, the present approach has been able to reproduce all the major aspects encountered when a flame is subject to an imposed electric field and the main effects of the different chemical mechanisms. Moreover, the proposed model is shown to produce a large reduction in the computational cost, being able to shorten the time needed to perform a simulation up to 40 times. △ Less

Submitted 11 August, 2017; originally announced August 2017.

Comments: 26 pages, 13 figures, paper accepted for publication on Computers and Fluids

A parallel interaction potential approach coupled with the immersed boundary method for fully resolved simulations of deformable interfaces and membranes

Authors: Vamsi Spandan, Valentina Meschini, Rodolfo Ostilla-Monico, Detlef Lohse, Giorgio Querzoli, Marco D de Tullio, Roberto Verzicco

Abstract: In this paper we show and discuss the use of a versatile interaction potential approach coupled with an immersed boundary method to simulate a variety of flows involving deformable bodies. In particular, we focus on two kinds of problems, namely (i) deformation of liquid-liquid interfaces and (ii) flow in the left ventricle of the heart with either a mechanical or a natural valve. Both examples ha… ▽ More In this paper we show and discuss the use of a versatile interaction potential approach coupled with an immersed boundary method to simulate a variety of flows involving deformable bodies. In particular, we focus on two kinds of problems, namely (i) deformation of liquid-liquid interfaces and (ii) flow in the left ventricle of the heart with either a mechanical or a natural valve. Both examples have in common the two-way interaction of the flow with a deformable interface or a membrane. The interaction potential approach (de Tullio & Pascazio, Jou. Comp. Phys., 2016; Tanaka, Wada and Nakamura, Computational Biomechanics, 2016) with minor modifications can be used to capture the deformation dynamics in both classes of problems. We show that the approach can be used to replicate the deformation dynamics of liquid-liquid interfaces through the use of ad-hoc elastic constants. The results from our simulations agree very well with previous studies on the deformation of drops in standard flow configurations such as deforming drop in a shear flow or a cross flow. We show that the same potential approach can also be used to study the flow in the left ventricle of the heart. The flow imposed into the ventricle interacts dynamically with the mitral valve (mechanical or natural) and the ventricle which are simulated using the same model. Results from these simulations are compared with ad- hoc in-house experimental measurements. Finally, a parallelisation scheme is presented, as parallelisation is unavoidable when studying large scale problems involving several thousands of simultaneously deforming bodies on hundreds of distributed memory computing processors. △ Less

Submitted 11 December, 2016; originally announced December 2016.

A combined Lattice Boltzmann and Immersed Boundary approach for predicting the vascular transport of differently shaped particles

Authors: Alessandro Coclite, Marco Donato de Tullio, Giuseppe Pascazio, Paolo Decuzzi

Abstract: Modelling the vascular transport and adhesion of man-made particles is crucial for optimizing their efficacy in the detection and treatment of diseases. Here, a Lattice Boltzmann and Immersed Boundary methods are combined together for predicting the near wall dynamics of particles with different shapes in a laminar flow. For the lattice Boltzmann modelling, a Gauss-Hermite projection is used to de… ▽ More Modelling the vascular transport and adhesion of man-made particles is crucial for optimizing their efficacy in the detection and treatment of diseases. Here, a Lattice Boltzmann and Immersed Boundary methods are combined together for predicting the near wall dynamics of particles with different shapes in a laminar flow. For the lattice Boltzmann modelling, a Gauss-Hermite projection is used to derive the lattice equation, wall boundary conditions are imposed through the Zou-He framework, and a moving least squares algorithm accurately reconstructs the forcing term accounting for the immersed boundary. First, the computational code is validated against two well-known test cases: the sedimentation of circular and elliptical cylinders in a quiescent fluid. A very good agreement is observed between the present results and those available in the literature. Then, the transport of circular, elliptical, rectangular, square and triangular particles is analyzed in a Couette flow, at Re=20. All particles drifted laterally across the stream lines reaching an equilibrium position, independently of the initial conditions. For this large Reynolds number, the particle shape has no significant effect on the final equilibrium position but it does affect the absolute value and periodicity of the angular velocity. Specifically, elongated particles show longer oscillation periods and, most interestingly, larger variations in angular velocity. The longest particles exhibit a zero angular velocity for almost the whole rotational period. Collectively, this data demonstrates that the proposed approach can be efficiently used for predicting complex particle dynamics in biologically relevant flows. This computational strategy could have significant impact in the field of computational nanomedicine for optimizing the specific delivery of therapeutic and imaging agents. △ Less

Submitted 6 June, 2016; originally announced June 2016.

LES of the Sandia Flame D Using an FPV Combustion Model

Authors: Mario Di Renzo, Alessandro Coclite, Marco Donato de Tullio, Pietro De Palma, Giuseppe Pascazio

Abstract: The simulation of turbulent combustion phenomena is still an open problem in modern fluid dynamics. Considering the economical importance of hydrocarbon combustion in energy production processes, it is evident the need of an accurate tool with a relatively low computational cost for the prediction of this kind of reacting flows. In the present work, a comparative study is carried out among large e… ▽ More The simulation of turbulent combustion phenomena is still an open problem in modern fluid dynamics. Considering the economical importance of hydrocarbon combustion in energy production processes, it is evident the need of an accurate tool with a relatively low computational cost for the prediction of this kind of reacting flows. In the present work, a comparative study is carried out among large eddy simulations, performed with various grid resolutions, a Reynolds averaged Navier-Stokes simulation, and experimental data concerning the well-known Sandia D flame test case. In all the simulations, a flamelet progress variable model has been employed using various hypotheses for the joint probability density function closure. The filtered approach proved to be more accurate than the averaged one, even for the coarser grid used in this work. In fact both approaches have shown poorly accurate predictions in the first part of the combustion chamber, but only by the large eddy simulation one is capable to recover the inlet discrepancies with respect to the experimental data going along the streamwise direction. △ Less

Submitted 1 February, 2016; originally announced February 2016.

Journal ref: Energy Procedia, 82 (2015), Pag: 402-409