Pub Date : 2024-11-18DOI: 10.1140/epje/s10189-024-00458-z
Francesco Michele Ventrella, Guido Boffetta, Massimo Cencini, Filippo De Lillo
We propose a simple numerical model for the motion of microswimmers based on the immersed boundary method. The swimmer, either pusher or puller, is represented by a distribution of point forces corresponding to the body and the flagellum. We study in particular the minimal model consisting of only three beads (two for the body and one for the flagellum) connected by rigid, inextensible links. When the beads are collinear, standard straight swimming is realized and, in the absence of propulsion, we demonstrate that the model recovers Jeffery’s equation for a thin rod. Conversely, by imposing an angle between body and flagellum the swimmer moves on circular orbits. We discuss how two swimmers, in collinear or non-collinear geometry, scatter upon encounter. Finally, we explore the dynamics of a large number of swimmers reacting to one another only via hydrodynamic interactions, and exemplify their complex collective dynamics in both straight and circular swimmers.
{"title":"Modeling straight and circle swimmers: from single swimmer to collective motion","authors":"Francesco Michele Ventrella, Guido Boffetta, Massimo Cencini, Filippo De Lillo","doi":"10.1140/epje/s10189-024-00458-z","DOIUrl":"10.1140/epje/s10189-024-00458-z","url":null,"abstract":"<p>We propose a simple numerical model for the motion of microswimmers based on the immersed boundary method. The swimmer, either pusher or puller, is represented by a distribution of point forces corresponding to the body and the flagellum. We study in particular the minimal model consisting of only three beads (two for the body and one for the flagellum) connected by rigid, inextensible links. When the beads are collinear, standard straight swimming is realized and, in the absence of propulsion, we demonstrate that the model recovers Jeffery’s equation for a thin rod. Conversely, by imposing an angle between body and flagellum the swimmer moves on circular orbits. We discuss how two swimmers, in collinear or non-collinear geometry, scatter upon encounter. Finally, we explore the dynamics of a large number of swimmers reacting to one another only via hydrodynamic interactions, and exemplify their complex collective dynamics in both straight and circular swimmers.</p>","PeriodicalId":790,"journal":{"name":"The European Physical Journal E","volume":"47 11","pages":""},"PeriodicalIF":1.8,"publicationDate":"2024-11-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142646615","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-10-16DOI: 10.1140/epje/s10189-024-00456-1
Jean-Baptiste Fournier
The pointlike curvature constraint (PCC) model and the disk detachment angle (DDA) model for the deformation-mediated interaction of conical integral protein inclusions in biomembranes are compared in the small deformation regime. Given the radius of membrane proteins, which is comparable to the membrane thickness, it is not obvious which of the two models should be considered the most adequate. For two proteins in a tensionless membranes, the PCC and DDA models coincide at the leading-order (sim r^{-4}) in their separation but differ at the next order. Yet, for distances larger than twice the proteins diameter, the difference is less than (10%). Like the DDA model, the PCC model includes all multibody interactions in a non-approximate way. The asymptotic (sim r^{-4}) many-body energy of triangular and square protein clusters is exactly the same in both models. Pentagonal clusters, however, behave differently; they have a vanishing energy in the PCC model, while they have a non-vanishing weaker (sim r^{-6}) asymptotic power law in the DDA model. We quantify the importance of multibody interactions in small polygonal clusters of three, four and five inclusions with identical or opposite curvatures in tensionless or tense membranes. We find that the pairwise approximation is almost always very poor. At short separation, the three-body interaction is not sufficient to account for the full many-body interaction. This is confirmed by equilibrium Monte Carlo simulations of up to ten inclusions.
{"title":"Multibody interactions between protein inclusions in the pointlike curvature model for tense and tensionless membranes","authors":"Jean-Baptiste Fournier","doi":"10.1140/epje/s10189-024-00456-1","DOIUrl":"10.1140/epje/s10189-024-00456-1","url":null,"abstract":"<p>The pointlike curvature constraint (PCC) model and the disk detachment angle (DDA) model for the deformation-mediated interaction of conical integral protein inclusions in biomembranes are compared in the small deformation regime. Given the radius of membrane proteins, which is comparable to the membrane thickness, it is not obvious which of the two models should be considered the most adequate. For two proteins in a tensionless membranes, the PCC and DDA models coincide at the leading-order <span>(sim r^{-4})</span> in their separation but differ at the next order. Yet, for distances larger than twice the proteins diameter, the difference is less than <span>(10%)</span>. Like the DDA model, the PCC model includes all multibody interactions in a non-approximate way. The asymptotic <span>(sim r^{-4})</span> many-body energy of triangular and square protein clusters is exactly the same in both models. Pentagonal clusters, however, behave differently; they have a vanishing energy in the PCC model, while they have a non-vanishing weaker <span>(sim r^{-6})</span> asymptotic power law in the DDA model. We quantify the importance of multibody interactions in small polygonal clusters of three, four and five inclusions with identical or opposite curvatures in tensionless or tense membranes. We find that the pairwise approximation is almost always very poor. At short separation, the three-body interaction is not sufficient to account for the full many-body interaction. This is confirmed by equilibrium Monte Carlo simulations of up to ten inclusions.</p>","PeriodicalId":790,"journal":{"name":"The European Physical Journal E","volume":"47 10","pages":""},"PeriodicalIF":1.8,"publicationDate":"2024-10-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142443254","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-10-15DOI: 10.1140/epje/s10189-024-00457-0
Amin Ullah, Jianzhong Lin, Yuxiang Yin
Sedimentation characteristics of a squirmer in a power-law fluid within a vertical channel are studied numerically using the two-dimensional lattice Boltzmann method. The effects of swimming type (− 5 ≤ β ≤ 5), self-propelling strength (0.5 ≤ α ≤ 1.1), power-law indexes (0.5 ≤ n ≤ 1.5), and the density ratio of the squirmer to the fluid (γ = 1.01, 1.5 and 2.3) on the sedimentation of the squirmer are analyzed. Four settlement patterns are identified: steady falling in the center, downward along the wall, oscillating with large amplitude and oscillating around the centerline. The squirmer in the channel exhibits more fluctuations in shear-thinning (n < 1) and Newtonian (n = 1) fluids compared to shear-thickening fluids (n > 1). Additionally, a puller (β > 0) settles faster than a pusher (β < 0) in shear-thinning and Newtonian fluids. Puller generates flow towards their head and away from their tail, exhibiting small amplitude oscillations. Pushers exhibit higher amplitude oscillations throughout the channel, creating flow towards their tail and away from their head. At lower γ, the fluctuation of the squirmer is less pronounced compared to higher γ.
{"title":"Research on sedimentation characteristics of squirmer in a power-law fluid","authors":"Amin Ullah, Jianzhong Lin, Yuxiang Yin","doi":"10.1140/epje/s10189-024-00457-0","DOIUrl":"10.1140/epje/s10189-024-00457-0","url":null,"abstract":"<div><p>Sedimentation characteristics of a squirmer in a power-law fluid within a vertical channel are studied numerically using the two-dimensional lattice Boltzmann method. The effects of swimming type (− 5 ≤ <i>β</i> ≤ 5), self-propelling strength (0.5 ≤ <i>α</i> ≤ 1.1), power-law indexes (0.5 ≤ <i>n</i> ≤ 1.5), and the density ratio of the squirmer to the fluid (<i>γ</i> = 1.01, 1.5 and 2.3) on the sedimentation of the squirmer are analyzed. Four settlement patterns are identified: steady falling in the center, downward along the wall, oscillating with large amplitude and oscillating around the centerline. The squirmer in the channel exhibits more fluctuations in shear-thinning (<i>n</i> < 1) and Newtonian (<i>n</i> = 1) fluids compared to shear-thickening fluids (<i>n</i> > 1). Additionally, a puller (<i>β</i> > 0) settles faster than a pusher (<i>β</i> < 0) in shear-thinning and Newtonian fluids. Puller generates flow towards their head and away from their tail, exhibiting small amplitude oscillations. Pushers exhibit higher amplitude oscillations throughout the channel, creating flow towards their tail and away from their head. At lower <i>γ</i>, the fluctuation of the squirmer is less pronounced compared to higher <i>γ</i>.</p><h3>Graphical abstract</h3><div><figure><div><div><picture><source><img></source></picture></div></div></figure></div></div>","PeriodicalId":790,"journal":{"name":"The European Physical Journal E","volume":"47 10","pages":""},"PeriodicalIF":1.8,"publicationDate":"2024-10-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142434831","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-10-14DOI: 10.1140/epje/s10189-024-00455-2
I. Tarras, A. Eddakoun, A. Hader, S. Moushi, I. Bakassi, R. Et Touizi, I. Achik, M. Eddahby, A. El Bachiri, Y. Boughaleb
Self-propelled collective motion is a highly complex phenomenon, necessitating advanced practical and theoretical tools for comprehension. The significance of studying collective motion becomes apparent in its diverse applications. For instance, addressing evacuation challenges in scenarios with multiple agents can be achieved through an examination of collective motion. Research indicates that the transition of individuals (such as birds, fish, etc.) from a state of rest to equilibrium constitutes a phase transition. Our interest of the issue is to delve into the nature of this transitional phase and elucidate the parameters that shape it. Hence, the primary aim of this paper is to grasp the kinetic phase transition by examining how initial velocity and repulsive interactions impact the dynamics of the system. To gain insight into the complex behavior of multi-agent systems, we apply an extended version of the classical Vicsek model. This extension includes an additional interaction zone, the repulsive zone, where particles repel each other at close range to avoid collisions. Our study uses numerical simulations to explore the system's behavior under various conditions. The focus of this study is the impact of initial velocity on the collective movement of particles. The importance of this research lies in comprehending how velocity affects the overall movement. The conclusion we can draw from these results is that the initial velocity affects both the noise and the density. The novelty of the work is the transition phase, yet it lacks universal characteristics because the critical noise depends on the initial velocity system and the repulsion radius zone. Notably, the repulsion radius and particle density play pivotal roles in achieving a phase transition from one equilibrium state to another aligned equilibrium state.