Search

We investigate the dynamics of plasma-based acceleration processes with collisionless particle dynamics and non negligible thermal effects. We aim at assessing the applicability of fluid-like models, obtained by suitable closure assumptions applied to the relativistic kinetic equations, thus not suffering of statistical noise, even in presence of a finite temperature. The work here presented focuses on the characterization of pressure anisotropies, which crucially depend on the adopted closure scheme, and hence are useful to discern the appropriate thermal fluid model. To this aim, simulation results of spatially resolved fluid models with different thermal closure assumptions are compared with the results of particle-in-cell (PIC) simulations at changing temperature and amplitude of plasma oscillations., Comment: The following article has been submitted to Physics of Plasma

The steady-state behavior of a single drop under shear flow has been extensively investigated in the limit of small deformation and negligible inertia effects. In this work, we combine the calculations proposed by Flumerfelt [R. W. Flumerfelt, J. Colloid Interface Sci. 76, 330 (1980)] for unconfined drops with interface viscosity, with those by Shapira & Haber [M. Shapira and S. Haber, Int. J. Multiph. Flow. 16, 305 (1990)] for confined drops without interface viscosity. By merging these two approaches, we provide comprehensive analytical predictions for steady-state drop deformation and inclination angle across a wide range of physical conditions, from confined to unconfined droplets, including or excluding the effect of interface viscosity. The proposed analytical predictions are also robust concerning variations in the viscosity ratio, making our model general enough to include any of the above conditions.

We develop a numerical method for simulating the dynamics of a droplet immersed in a generic time-dependent velocity gradient field. This approach is grounded on the hybrid coupling between the lattice Boltzmann (LB) method, employed for the flow simulation, and the immersed boundary (IB) method, utilized to couple the droplet with the surrounding fluid. We show how to enrich the numerical scheme with a mesh regularization technique, allowing droplets to sustain large deformations. The resulting methodology is adapted to simulate the dynamics of droplets in homogeneous and isotropic turbulence, with the characteristic size of the droplet being smaller than the characteristic Kolmogorov scale of the outer turbulent flow. We report on statistical results for droplet deformation and orientation, collected from an ensemble of turbulent trajectories, as well as comparisons with theoretical models in the limit of small deformation.

A comprehensive characterization of lattice Boltzmann (LB) schemes to perform warm fluid numerical simulations of particle wakefield acceleration (PWFA) processes is discussed in this paper. The LB schemes we develop hinge on the moment matching procedure, allowing the fluid description of a warm relativistic plasma wake generated by a driver pulse propagating in a neutral plasma. We focus on fluid models equations resulting from two popular closure assumptions of the relativistic kinetic equations, i.e., the local equilibrium and the warm plasma closure assumptions. The developed LB schemes can thus be used to disclose insights on the quantitative differences between the two closure approaches in the dynamics of PWFA processes. Comparisons between the proposed schemes and available analytical results are extensively addressed., Comment: Published on Physics of Plasmas

The wetting dynamics of liquid particles, from coated droplets to soft capsules, holds significant technological interest. Motivated by the need to simulate liquid metal droplet with an oxidize surface layer, in this work we introduce a computational scheme that allows to simulate droplet dynamics with general surface properties and model different levels of interface stiffness, describing also cases that are intermediate between pure droplets and capsules. Our approach is based on a combination of the immersed boundary (IB) and the lattice Boltzmann (LB) methods. Here, we validate our approach against the theoretical predictions in the context of shear flow and static wetting properties and we show its effectiveness in accessing the wetting dynamics, exploring the ability of the scheme to address a broad phenomenology., Comment: 14 pages, 6 figures

Membrane viscosity is known to play a central role in the transient dynamics of isolated viscoelastic capsules by decreasing their deformation, inducing shape oscillations and reducing the loading time, that is, the time required to reach the steady-state deformation. However, for dense suspensions of capsules, our understanding of the influence of the membrane viscosity is minimal. In this work, we perform a systematic numerical investigation based on coupled immersed boundary -- lattice Boltzmann (IB-LB) simulations of viscoelastic spherical capsule suspensions in the non-inertial regime. We show the effect of the membrane viscosity on the transient dynamics as a function of volume fraction and capillary number. Our results indicate that the influence of membrane viscosity on both deformation and loading time strongly depends on the volume fraction in a non-trivial manner: dense suspensions with large surface viscosity are more resistant to deformation but attain loading times that are characteristic of capsules with no surface viscosity, thus opening the possibility to obtain richer combinations of mechanical features.

We propose an extension of the Maffettone-Minale (MM) model to predict the dynamics of an ellipsoidal droplet in a shear flow. The parameters of the MM model are traditionally retrieved in the framework of the perturbation theory for small deformations, i.e., small capillary numbers ($\mbox{Ca} \ll 1$) applied to Stokes equations. In this work, we take a novel route, in that we determine the model parameters at finite capillary numbers ($\mbox{Ca}\sim {\cal O}(1)$) without relying on perturbation theory results, while retaining a realistic representation in creeping time and steady deformation attained by the droplet for different realizations of the viscosity ratio $\lambda$ between the inner and the outer fluids. This extended MM (EMM) model hinges on an independent characterization of the process of droplet deformation via numerical simulations of Stokes equations employing the Immersed Boundary - Lattice Boltzmann (IB-LB) numerical techniques. Issues on droplet breakup are also addressed and discussed within the EMM model., Comment: 26 pages, 10 figures

We explore a novel simulation route for Plasma Wakefield Acceleration (PWFA) by using the computational method known as the Lattice Boltzmann Method (LBM). LBM is based on a discretization of the continuum kinetic theory while assuring the convergence towards hydrodynamics for coarse-grained fields (i.e., density, velocity, etc.). LBM is an established numerical analysis tool in computational fluid dynamics, able to efficiently bridge between kinetic theory and hydrodynamics, but its application in the context of PWFA has never been investigated so far. This paper aims at filling this gap. Results of LBM simulations for PWFA are discussed and compared with those of a code (Architect) implementing a Cold Fluid (CF) model for the plasma. In the hydrodynamic framework, we discuss the importance of regularization effects related to diffusion properties intrinsic of the LBM, allowing to go beyond the traditional CF approximations. Issues on computational efficiency are also addressed.

We use mesoscale numerical simulations to investigate the unsteady dynamics of a single red blood cell (RBC) subjected to an external mechanical load. We carry out a detailed comparison between the {\it loading} (L) dynamics, following the imposition of the mechanical load on the RBC at rest, and the {\it relaxation} (R) dynamics, allowing the RBC to relax to its original shape after the sudden arrest of the mechanical load. Such a comparison is carried out by analyzing the characteristic times of the two corresponding dynamics, i.e., $t_L$ and $t_R$. When the intensity of the mechanical load is small enough, the two kinds of dynamics are {\it symmetrical} ($t_L \approx t_R$) and independent of the typology of mechanical load (intrinsic dynamics); otherwise, in marked contrast, an {\it asymmetry} is found, wherein the loading dynamics is typically faster than the relaxation one. This asymmetry manifests itself with non-universal characteristics, e.g., dependency on the applied load and/or on the viscoelastic properties of the RBC membrane. To deepen such a non-universal behaviour, we consider the viscosity of the erythrocyte membrane as a variable parameter and focus on three different typologies of mechanical load (mechanical stretching, shear flow, elongational flow): this allows to clarify how non-universality builds up in terms of the deformation and rotational contributions induced by the mechanical load on the membrane. Finally, we also investigate the effect of the elastic shear modulus on the characteristic times $t_L$ and $t_R$. Our results provide crucial and quantitative information on the unsteady dynamics of RBC and its membrane response to the imposition/cessation of external mechanical loads.

The tumbling to tank-treading (TB-TT) transition for red blood cells (RBCs) has been widely investigated, with a main focus on the effects of the viscosity ratio $\lambda$ (i.e., the ratio between the viscosities of the fluids inside and outside the membrane) and the shear rate $\dot{\gamma}$ applied to the RBC. However, the membrane viscosity $\mu_m$ plays a major role in a realistic description of RBC's dynamics, and only a few works have systematically focused on its effects on the TB-TT transition. In this work, we provide a parametric investigation on the effect of membrane viscosity $\mu_m$ on the TB-TT transition, for a single RBC. It is found that, at fixed viscosity ratios $\lambda$, larger values of $\mu_m$ lead to an increased range of values of capillary number at which the TB-TT transition occurs. We systematically quantify such an increase by means of mesoscale numerical simulations based on the lattice Boltzmann models.