-
Langevin and Navier-Stokes Simulation of Three-Dimensional Protoplasmic Streaming
Authors:
Shuta Noro,
Satoshi Hongo,
Shinichiro Nagahiro,
Hisatoshi Ikai,
Hiroshi Koibuchi,
Madoka Nakayama,
Tetsuya Uchimoto,
Gildas Diguet
Abstract:
In this paper, we report the numerical results obtained using the Langevin Navier-Stokes (LNS) simulation of the velocity distribution of three-dimensional (3D) protoplasmic streaming in plant cells, such as those of {\it Nitella flexilis}. The LNS simulations are performed on 3D cylinders discretized by regular cubes in which fluid velocities are activated by boundary velocities parallel and nonp…
▽ More
In this paper, we report the numerical results obtained using the Langevin Navier-Stokes (LNS) simulation of the velocity distribution of three-dimensional (3D) protoplasmic streaming in plant cells, such as those of {\it Nitella flexilis}. The LNS simulations are performed on 3D cylinders discretized by regular cubes in which fluid velocities are activated by boundary velocities parallel and nonparallel to the longitudinal direction and a random Brownian force with strength $D$. We find that, for a finite $D$, the velocity distribution $h(V), V\!=\!|\vec{V}|$, has two different peaks at a small non-zero $V$ and a finite $V$, and the distribution $h(V_z)$ for $|V_z|$ along the longitudinal direction also has a peak at finite $V_z$. These results are in good agreement with the reported velocity distributions observed using laser Doppler velocimetry. Moreover, we study the effects of the Brownian force on biological material mixing and find that mixing along the $\vec{V}$ direction enhanced by the nonparallel circular motion is further improved by the Brownian force in the experimentally relevant region of $D$. In addition, the experimentally relevant $D$ is found to be consistent with the expectation from the fluctuation dissipation relation between the random stress and viscosity in the LNS equation of Landau and Lifschitz for incompressible fluids.
△ Less
Submitted 19 August, 2023; v1 submitted 26 December, 2021;
originally announced December 2021.
-
Langevin Navier-Stokes simulation of protoplasmic streaming by 2D MAC method
Authors:
Shuta Noro,
Satoshi Hongo,
Shinichiro Nagahiro,
Hisatoshi Ikai,
Madoka Nakayama,
Hiroshi Koibuchi,
Tetsuya Uchimoto,
Jean-Paul Rieu
Abstract:
We study protoplasmic streaming in plant cells such as chara brauni by simplifying the flow field to a two-dimensional Couette flow with Brownian random motion inside parallel plates. Protoplasmic streaming is receiving a lot of attention in many areas, such as agriculture-technology and biotechnology. The plant size depends on the velocity of streaming and the driving force originating in molecul…
▽ More
We study protoplasmic streaming in plant cells such as chara brauni by simplifying the flow field to a two-dimensional Couette flow with Brownian random motion inside parallel plates. Protoplasmic streaming is receiving a lot of attention in many areas, such as agriculture-technology and biotechnology. The plant size depends on the velocity of streaming and the driving force originating in molecular motors. Therefore, it is interesting to study detailed information on the velocity of streaming. Recently, experimentally observed peaks in the velocity distribution have been simulated by a 2D Langevin Navier-Stokes (LNS) equation for vortex and flow function. However, to simulate actual 3D flows, we have to use the NS equation for velocity, which, in the case of 2D flows, is not always equivalent to that for vorticity and stream function. In this paper, we report that a 2D LNS equation for velocity and pressure successfully simulates protoplasmic streaming by comparing the results with the experimental data and those obtained by 2D LNS simulations for vortex and flow function. Moreover, a dimensional analysis clarifies the dependence of numerical results on the strength $D$ of Brownian random force and physical parameters such as kinematic viscosity and cell size. We find from this analysis how the peak position in normalized velocity distribution moves depending on these parameters.
△ Less
Submitted 8 June, 2022; v1 submitted 20 December, 2021;
originally announced December 2021.
-
Stochastic Fluid Dynamics Simulations of the Velocity Distribution in Protoplasmic Streaming
Authors:
Vladislav Egorov,
Olga Maksimova,
Irina Andreeva,
Hiroshi Koibuchi,
Satoshi Hongo,
Shinichiro Nagahiro,
Toshiyuki Ikai,
Madoka Nakayama,
Shuta Noro,
Tetsuya Uchimoto,
Jean-Paul Rieu
Abstract:
Protoplasmic streaming in plant cells is directly visible in the cases of \textit{Chara corallina} and \textit{Nitella flexilis}, and this streaming is understood to play a role in the transport of biological materials. For this reason, related studies have focused on molecular transportation from a fluid mechanics viewpoint. However, the experimentally observed distribution of the velocity along…
▽ More
Protoplasmic streaming in plant cells is directly visible in the cases of \textit{Chara corallina} and \textit{Nitella flexilis}, and this streaming is understood to play a role in the transport of biological materials. For this reason, related studies have focused on molecular transportation from a fluid mechanics viewpoint. However, the experimentally observed distribution of the velocity along the flow direction $x$, which exhibits two peaks at $V_x\!=\!0$ and at a finite $V_x(\not=\!0)$, remains to be studied. In this paper, we numerically study whether this behavior of the flow field can be simulated by a 2D stochastic Navier-Stokes (NS) equation for Couette flow, in which random Brownian force is assumed. We present the first numerical evidence that these peaks are reproduced by the stochastic NS equation, which implies that the Brownian motion of the fluid particles plays an essential role in the emergence of these peaks in the velocity distribution. We also find that the position of the peak at $V_x(\not=\!0)$ moves with the variation in the strength $D$ of the random Brownian force, which also changes depending on physical parameters such as the kinematic viscosity, boundary velocity and diameter of the plant cells.
△ Less
Submitted 27 November, 2020; v1 submitted 22 June, 2020;
originally announced June 2020.
-
Negative pressure in shear thickening band of a dilatant fluid
Authors:
Shin-ichiro Nagahiro,
Hiizu Nakanishi
Abstract:
We perform experiments and numerical simulations to investigate spatial distribution of pressure in a sheared dilatant fluid of the Taylor-Couette flow under a constant external shear stress. In a certain range of shear stress, the flow undergoes the shear thickening oscillation around 20 Hz. We find that, during the oscillation, a localized thickened band rotates around the axis with the flow. {B…
▽ More
We perform experiments and numerical simulations to investigate spatial distribution of pressure in a sheared dilatant fluid of the Taylor-Couette flow under a constant external shear stress. In a certain range of shear stress, the flow undergoes the shear thickening oscillation around 20 Hz. We find that, during the oscillation, a localized thickened band rotates around the axis with the flow. {Based upon experiments and numerical simulations, we show that a major part of the thickened band is under negative pressure even in the case of discontinuous shear thickening}, which indicates that the thickening is caused by Reynolds dilatancy; the dilatancy causes the negative pressure in interstitial fluid, which generates contact structure in the granular medium., then frictional resistance hinders rearrangement of the structure and solidifies the medium.
△ Less
Submitted 13 December, 2016; v1 submitted 16 November, 2015;
originally announced November 2015.
-
Experimental observation of shear thickening oscillation
Authors:
Shin-ichiro Nagahiro,
Hiizu Nakanishi,
Namiko Mitarai
Abstract:
We report experimental observation of the shear thickening oscillation, i.e. the spontaneous macroscopic oscillation in the shear flow of severe shear thickening fluid. The shear thickening oscillation is caused by the interplay between the fluid dynamics and the shear thickening, and has been predicted theoretically by the present authors using a phenomenological fluid dynamics model for the dila…
▽ More
We report experimental observation of the shear thickening oscillation, i.e. the spontaneous macroscopic oscillation in the shear flow of severe shear thickening fluid. The shear thickening oscillation is caused by the interplay between the fluid dynamics and the shear thickening, and has been predicted theoretically by the present authors using a phenomenological fluid dynamics model for the dilatant fluid, but never been reported experimentally. Using a density-matched starch-water mixture, in the cylindrical shear flow of a few centimeters flow width, we observed strong vibrations of the frequency around 20 Hz, which is consistent with our theoretical prediction.
△ Less
Submitted 9 July, 2013; v1 submitted 12 November, 2012;
originally announced November 2012.
-
Fluid dynamics of dilatant fluid
Authors:
Hiizu Nakanishi,
Shin-ichiro Nagahiro,
Namiko Mitarai
Abstract:
Dense mixture of granules and liquid often shows a sever shear thickening and is called a dilatant fluid. We construct a fluid dynamics model for the dilatant fluid by introducing a phenomenological state variable for a local state of dispersed particles. With simple assumptions for an equation of the state variable, we demonstrate that the model can describe basic features of the dilatant fluid s…
▽ More
Dense mixture of granules and liquid often shows a sever shear thickening and is called a dilatant fluid. We construct a fluid dynamics model for the dilatant fluid by introducing a phenomenological state variable for a local state of dispersed particles. With simple assumptions for an equation of the state variable, we demonstrate that the model can describe basic features of the dilatant fluid such as the stress-shear rate curve that represents discontinuous severe shear thickening, hysteresis upon changing shear rate, instantaneous hardening upon external impact. Analysis of the model reveals that the shear thickening fluid shows an instability in a shear flow for some regime and exhibits {\it the shear thickening oscillation}, i.e. the oscillatory shear flow alternating between the thickened and the relaxed states. Results of numerical simulations are presented for one and two-dimensional systems.
△ Less
Submitted 20 December, 2011;
originally announced December 2011.
-
Model for coiling and meandering instability of viscous threads
Authors:
S. Nagahiro,
Y. Hayakawa
Abstract:
A numerical model is presented to describe both the transient and steady-state dynamics of viscous threads falling onto a plane. The steady-state coiling frequency w is calculated as a function of fall height H. In the case of weak gravity, w ~ H^{-1} and w ~ H are obtained for lower and higher fall heights respectively. When the effect of gravity is significant, the relation w ~ H^2 is observed…
▽ More
A numerical model is presented to describe both the transient and steady-state dynamics of viscous threads falling onto a plane. The steady-state coiling frequency w is calculated as a function of fall height H. In the case of weak gravity, w ~ H^{-1} and w ~ H are obtained for lower and higher fall heights respectively. When the effect of gravity is significant, the relation w ~ H^2 is observed. These results agree with the scaling laws previously predicted. The critical Reynolds number for coil-uncoil transition is discussed. When the gravity is weak, the transition occurs with hysteresis effects. If the plane moves horizontally at a constant speed, a variety of meandering oscillation modes can be observed experimentally. The present model also can describe this phenomenon. The numerically obtained state diagram for the meandering modes qualitatively agrees with the experiment.
△ Less
Submitted 4 March, 2009; v1 submitted 24 December, 2008;
originally announced December 2008.
-
Bending-Filament Model for the Buckling and Coiling Instability of Viscous Fluid Rope
Authors:
Shin-ichiro Nagahiro,
Yoshinori Hayakawa
Abstract:
A simple model is proposed for the buckling and coiling instability of a viscous "fluid rope" falling on a plane. By regarding a fluid rope as a one-dimensional flow, this model accounts for only the axial and shared viscous forces. Our model successfully reproduces several experiments with no adjustable parameters, such as the existence of three distinct coiling regimes reported in Phys. Rev. L…
▽ More
A simple model is proposed for the buckling and coiling instability of a viscous "fluid rope" falling on a plane. By regarding a fluid rope as a one-dimensional flow, this model accounts for only the axial and shared viscous forces. Our model successfully reproduces several experiments with no adjustable parameters, such as the existence of three distinct coiling regimes reported in Phys. Rev. Lett. 93, 214502 (2004). Our model allows for the discussion of unsteady motion. An expression for the critical fall height at which the coiling frequency changes from a decrease to increase was phenomenologically derived. It was found that the coil-uncoil transition shows remarkable hysteresis only for weak gravity condition.
△ Less
Submitted 23 May, 2008; v1 submitted 12 February, 2008;
originally announced February 2008.
-
Separation of long DNA chains using non-uniform electric field: a numerical study
Authors:
Shin-ichiro Nagahiro,
Satoyuki Kawano,
Hidetoshi Kotera
Abstract:
We study migration of DNA molecules through a microchannel with a series of electric traps controlled by an ac electric field. We describe the motion of DNA based on Brownian dynamics simulations of a beads-spring chain. Our simulation demonstrates that the chain captured by an electrode escapes from the binding electric field due to thermal fluctuation. We find that the mobility of chain would…
▽ More
We study migration of DNA molecules through a microchannel with a series of electric traps controlled by an ac electric field. We describe the motion of DNA based on Brownian dynamics simulations of a beads-spring chain. Our simulation demonstrates that the chain captured by an electrode escapes from the binding electric field due to thermal fluctuation. We find that the mobility of chain would depend on the chain length; the mobility sharply increases when the length of a chain exceeds a critical value, which is strongly affected by the amplitude of the applied ac field. Thus we can adjust the length regime, in which this microchannel well separates DNA molecules, without changing the structure of the channel. We also present a theoretical insight into the relation between the critical chain length and the field amplitude.
△ Less
Submitted 24 July, 2006; v1 submitted 17 July, 2006;
originally announced July 2006.
-
Collision of One-Dimensional Nonlinear Chains
Authors:
Shin-ichiro Nagahiro,
Yoshinori Hayakawa
Abstract:
We investigate one-dimensional collisions of unharmonic chains and a rigid wall. We find that the coefficient of restitution (COR) is strongly dependent on the velocity of colliding chains and has a minimum value at a certain velocity. The relationship between COR and collision velocity is derived for low-velocity collisions using perturbation methods. We found that the velocity dependence is ch…
▽ More
We investigate one-dimensional collisions of unharmonic chains and a rigid wall. We find that the coefficient of restitution (COR) is strongly dependent on the velocity of colliding chains and has a minimum value at a certain velocity. The relationship between COR and collision velocity is derived for low-velocity collisions using perturbation methods. We found that the velocity dependence is characterized by the exponent of the lowest unharmonic term of interparticle potential energy.
△ Less
Submitted 17 October, 2002;
originally announced October 2002.