Source author record

Leonid Berlyand

Leonid Berlyand appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

13 published item(s)

preprint2022arXiv

Dynamics and steady state of squirmer motion in liquid crystal

We analyze a nonlinear PDE system describing the motion of a microswimmer in a nematic liquid crystal environment. For the microswimmer's motility, the squirmer model is used in which self-propulsion enters the model through the slip velocity on the microswimmer's surface. The liquid crystal is described using the well-established Beris-Edwards formulation. In previous computational studies, it was shown that the squirmer, regardless of its initial configuration, eventually orients itself either parallel or perpendicular to the preferred orientation dictated by the liquid crystal. Furthermore, the corresponding solution of the coupled nonlinear system converges to a steady state. In this work, we rigorously establish the existence of steady state and also the finite-time existence for the time-dependent problem. Finally, using a two-scale asymptotic expansion we derive a reduced model for the collective swimming of squirmers as they reach their steady state orientation and speed.

preprint2021arXiv

Stability of Contraction-Driven Cell Motion

We consider motility of keratocyte cells driven by myosin contraction and introduce a 2D free boundary model for such motion. This model generalizes a 1D model from [12] by combining a 2D Keller-Segel model and a Hele-Shaw type boundary condition with the Young-Laplace law resulting in a boundary curvature term which provides a regularizing effect. We show that this model has a family of traveling solutions with constant shape and velocity which bifurcates from a family of radially symmetric stationary states. Our goal is to establish observable steady motion of the cell with constant velocity. Mathematically, this amounts to establishing stability of the traveling solutions. Our key result is an explicit asymptotic formula for the stability-determining eigenvalue of the linearized problem. This formula greatly simplifies the task of numerically computing the sign of this eigenvalue and reveals the physical mechanisms of stability. The derivation of this formula is based on a special ansatz for the corresponding eigenvector which exhibits an interesting singular behavior such that it asymptotically (in the small-velocity limit) becomes parallel to another eigenvector. This reflects the non-self-adjoint nature of the linearized problem, a signature of living systems. Finally, our results describe the onset of motion via a transition from unstable radial stationary solutions to stable asymmetric traveling solutions.

preprint2016arXiv

Flagella bending affects macroscopic properties of bacterial suspensions

To survive in harsh conditions, motile bacteria swim in complex environment and respond to the surrounding flow. Here we develop a PDE model describing how the flagella bending affects macroscopic properties of bacterial suspensions. First, we show how the flagella bending contributes to the decrease of the effective viscosity observed in dilute suspension. Our results do not impose tumbling (random re-orientation) as it was done previously to explain the viscosity reduction. Second, we demonstrate a possibility of bacterium escape from the wall entrapment due to the self-induced buckling of flagella. Our results shed light on the role of flexible bacterial flagella in interactions of bacteria with shear flow and walls or obstacles.

preprint2016arXiv

Phase field model of cell motility: sharp interface limit in sub-critical case

We consider a system of two PDEs arising in modeling of motility of eukariotic cells on substrates. This system consists of the Allen-Cahn equation for the scalar phase field function coupled with another vectorial parabolic equation for the orientation of the actin filament network. The two key properties of this system are (i) presence of gradients in the coupling terms (gradient coupling) and (ii) mass (volume) preservation constraints. We first prove that the sharp interface property of initial conditions is preserved in time. Next we formally derive the equation of the motion of the interface, which is the mean curvature motion perturbed by a nonlinear term that appears due to the properties (i)-(ii). This novel term leads to surprising features of the the motion of the interface. Because of these properties maximum principle and classical comparison techniques do not apply to this system. Furthermore, the system can not be written in a form of gradient flow, which is why recently developed Gamma-convergence techniques also can not be used for the justification of the formal derivation. Such justification is presented in a one-dimensional model problem and it leads to a stability result in a class of 'sharp interface' initial data.

preprint2016arXiv

Phase-Field Model of Cell Motility: Traveling Waves and Sharp Interface Limit

This letter is concerned with asymptotic analysis of a PDE model for motility of a eukaryotic cell on a substrate. This model was introduced in [1], where it was shown numerically that it successfully reproduces experimentally observed phenomena of cell-motility such as a discontinuous onset of motion and shape oscillations. The model consists of a parabolic PDE for a scalar phase-field function coupled with a vectorial parabolic PDE for the actin filament network (cytoskeleton). We formally derive the sharp interface limit (SIL), which describes the motion of the cell membrane and show that it is a volume preserving curvature driven motion with an additional nonlinear term due to adhesion to the substrate and protrusion by the cytoskeleton. In a 1D model problem we rigorously justify the SIL, and, using numerical simulations, observe some surprising features such as discontinuity of interface velocities and hysteresis. We show that nontrivial traveling wave solutions appear when the key physical parameter exceeds a certain critical value and the potential in the equation for phase field function possesses certain asymmetry.

preprint2016arXiv

Sharp interface limit in a phase field model of cell motility

We consider a system of two coupled parabolic PDEs introduced in [1] to model motility of eukaryotic cells. We study the asymptotic behavior of solutions in the limit of a small parameter related to the width of the interface in phase field function (sharp interface limit). We formally derive an equation of motion of the interface, which is mean curvature motion with an additional nonlinear term. In a 1D model parabolic problem we rigorously justify the sharp interface limit. To this end, a special form of asymptotic expansion is introduced to reduce analysis to a single nonlinear PDE. Further stability analysis reveals a qualitative change in the behavior of the system for small and large values of the coupling parameter. Using numerical simulations we also show discontinuities of the interface velocity and hysteresis. Also, in the 1D case we establish nontrivial traveling waves when the coupling parameter is large enough.

preprint2015arXiv

Effective Rheological Properties in Semidilute Bacterial Suspensions

Interactions between swimming bacteria have led to remarkable experimentally observable macroscopic properties such as the reduction of the effective viscosity, enhanced mixing, and diffusion. In this work, we study an individual based model for a suspension of interacting point dipoles representing bacteria in order to gain greater insight into the physical mechanisms responsible for the drastic reduction in the effective viscosity. In particular, asymptotic analysis is carried out on the corresponding kinetic equation governing the distribution of bacteria orientations. This allows one to derive an explicit asymptotic formula for the effective viscosity of the bacterial suspension in the limit of bacterium non-sphericity. The results show good qualitative agreement with numerical simulations and previous experimental observations. Finally, we justify our approach by proving existence, uniqueness, and regularity properties for this kinetic PDE model.

preprint2015arXiv

On an evolution equation in a cell motility model

This paper deals with the evolution equation of a curve obtained as the sharp interface limit of a non-linear system of two reaction-diffusion PDEs. This system was introduced as a phase-field model of (crawling) motion of eukaryotic cells on a substrate. The key issue is the evolution of the cell membrane (interface curve) which involves shape change and net motion. This issue can be addressed both qualitatively and quantitatively by studying the evolution equation of the sharp interface limit for this system. However, this equation is non-linear and non-local and existence of solutions presents a significant analytical challenge. We establish existence of solutions for a wide class of initial data in the so-called subcritical regime. Existence is proved in a two step procedure. First, for smooth ($H^2$) initial data we use a regularization technique. Second, we consider non-smooth initial data that are more relevant from the application point of view. Here, uniform estimates on the time when solution exists rely on a maximum principle type argument.

preprint2014arXiv

Complexity reduction in many particles systems with random initial data

We consider the motion of interacting particles governed by a coupled system of ODEs with random initial conditions. Direct computations for such systems are prohibitively expensive due to a very large number of particles and randomness requiring many realizations in their locations in the presence of strong interactions. While there are several approaches that address the above difficulties, none addresses all three simultaneously. Our goal is to develop such a computational approach in order to capture the experimentally observed emergence of correlations in the collective state (patterns due to strong interactions). Our approach is based on the truncation of the BBGKY hierarchy that allows one to go beyond the classical Mean Field limit and capture correlations while drastically reducing the computational complexity. Finally, we provide an example showing a numerical solution of this nonlinear and non-local system.

preprint2012arXiv

Boundary conditions in PDE model of collisions of swimmers

The goal of the paper is to determine boundary conditions in PDE models of collisions of microswimmers in a viscous fluid. We consider two self-propelled spheres (microswimmers) moving towards each other in viscous fluid. We first show that under commonly used no-slip boundary conditions on the fluid-solid interface the microswimmers do not collide which is a generalization of the well-known no-collision paradox for solid bodies (with no self-propulsion) in a viscous fluid. Secondly, we show that the microswimmers do collide when the no-slip boundary conditions are replaced by the Navier boundary conditions which therefore provides an adequate model of microswimmers such as swimming bacteria. The self-propulsion mechanism generates a drag force pulling a bacterium backwards and the collision problem is reduced to the analysis of competition between the drag and self-propulsion. For no-slip this is done by utilizing the Lorentz Reciprocal Theorem and the analytical solution for two solid spheres in the fluid. The analysis for the Navier boundary conditions is based on the variational formulation of Stokes problem. A Poincare type inequality for symmetrized gradient is introduced in this work.

preprint2010arXiv

Flux norm approach to finite dimensional homogenization approximations with non-separated scales and high contrast

We consider divergence-form scalar elliptic equations and vectorial equations for elasticity with rough ($L^\infty(Ω)$, $Ω\subset \R^d$) coefficients $a(x)$ that, in particular, model media with non-separated scales and high contrast in material properties. We define the flux norm as the $L^2$ norm of the potential part of the fluxes of solutions, which is equivalent to the usual $H^1$-norm. We show that in the flux norm, the error associated with approximating, in a properly defined finite-dimensional space, the set of solutions of the aforementioned PDEs with rough coefficients is equal to the error associated with approximating the set of solutions of the same type of PDEs with smooth coefficients in a standard space (e.g., piecewise polynomial). We refer to this property as the {\it transfer property}. A simple application of this property is the construction of finite dimensional approximation spaces with errors independent of the regularity and contrast of the coefficients and with optimal and explicit convergence rates. This transfer property also provides an alternative to the global harmonic change of coordinates for the homogenization of elliptic operators that can be extended to elasticity equations. The proofs of these homogenization results are based on a new class of elliptic inequalities which play the same role in our approach as the div-curl lemma in classical homogenization.

preprint2010arXiv

Radially symmetric minimizers for a $p$-Ginzburg Landau type energy in $\R^2$

We consider the minimization of a p-Ginzburg-Landau energy functional over the class of radially symmetric functions of degree one. We prove the existence of a unique minimizer in this class, and show that its modulus is monotone increasing and concave. We also study the asymptotic limit of the minimizers as p \rightarrow \infty. Finally, we prove that the radially symmetric solution is locally stable for $p$ in the interval $(2,4]$.