The collective dynamics of self-propelled particles

We propose a method for the dynamic simulation of a collection of self-propelled particles in a viscous Newtonian fluid. We restrict attention to particles whose size and velocity are small enough that the fluid motion is in the creeping flow regime. We propose a simple model for a self-propelled particle, and extended the Stokesian Dynamics method to conduct dynamic simulations of a collection of such particles. In our description, each particle is treated as a sphere with an orientation vector p, whose locomotion is driven by the action of a force dipole Sp of constant magnitude S0 at a point slightly displaced from its centre. To simplify the calculation, we place the dipole at the centre of the particle, and introduce a virtual propulsion force Fp to effect propulsion. The magnitude F0 of this force is proportional to S0. The directions of Sp and Fp are determined by p. In isolation, a self-propelled particle moves at a constant velocity u0 p, with the speed u0 determined by S0. When it coexists with many such particles, its hydrodynamic interaction with the other particles alters its velocity and, more importantly, its orientation. As a result, the motion of the particle is chaotic. Our simulations are not restricted to low particle concentration, as we implement the full hydrodynamic interactions between the particles, but we restrict the motion of particles to two dimensions to reduce computation. We have studied the statistical properties of a suspension of self-propelled particles for a range of the particle concentration, quantified by the area fraction φa. We find several interesting features in the microstructure and statistics. We find that particles tend to swim in clusters wherein they are in close proximity. Consequently, incorporating the finite size of the particles and the near-field hydrodynamic interactions is of the essence. There is a continuous process of breakage and formation of the clusters. We find that the distributions of particle velocity at low and high φa are qualitatively different; it is close to the normal distribution at high φa, in agreement with experimental measurements. The motion of the particles is diffusive at long time, and the self-diffusivity decreases with increasing φa. The pair correlation function shows a large anisotropic build-up near contact, which decays rapidly with separation. There is also an anisotropic orientation correlation near contact, which decays more slowly with separation. Movies are available with the online version of the paper.

Particle dynamics in a turbulent particle–gas suspension at high Stokes number. Part 1. Velocity and acceleration distributions

The effect of fluid velocity fluctuations on the dynamics of the particles in a turbulent gas–solid suspension is analysed in the low-Reynolds-number and high Stokes number limits, where the particle relaxation time is long compared with the correlation time for the fluid velocity fluctuations, and the drag force on the particles due to the fluid can be expressed by the modified Stokes law. The direct numerical simulation procedure is used for solving the Navier–Stokes equations for the fluid, the particles are modelled as hard spheres which undergo elastic collisions and a one-way coupling algorithm is used where the force exerted by the fluid on the particles is incorporated, but not the reverse force exerted by the particles on the fluid. The particle mean and root-mean-square (RMS) fluctuating velocities, as well as the probability distribution function for the particle velocity fluctuations and the distribution of acceleration of the particles in the central region of the Couette (where the velocity profile is linear and the RMS velocities are nearly constant), are examined. It is found that the distribution of particle velocities is very different from a Gaussian, especially in the spanwise and wall-normal directions. However, the distribution of the acceleration fluctuation on the particles is found to be close to a Gaussian, though the distribution is highly anisotropic and there is a correlation between the fluctuations in the flow and gradient directions. The non-Gaussian nature of the particle velocity fluctuations is found to be due to inter-particle collisions induced by the large particle velocity fluctuations in the flow direction. It is also found that the acceleration distribution on the particles is in very good agreement with the distribution that is calculated from the velocity fluctuations in the fluid, using the Stokes drag law, indicating that there is very little correlation between the fluid velocity fluctuations and the particle velocity fluctuations in the presence of one-way coupling. All of these results indicate that the effect of the turbulent fluid velocity fluctuations can be accurately represented by an anisotropic Gaussian white noise.

Dynamics of bound dyes in a nonpolymeric aqueous gel derived from a tripodal bile salt

Rotational dynamics of polarity sensitive fluorescent dyes (ANS and DPH) in a nonpolymertic aqueous gel derived from tripodal cholamide I was studied using ultrafast time-resolved fluorescence technique. Results were compared with that of naturally occurring di- and trihydroxy bile salts. ANS in the gel showed two rotational correlation time (phi) components, 13.2 ns (bound to the hydrophobic region of the gel) and 1.0 ns (free aqueous ANS), whereas DPH showed only one component (4.8 ns). In the sol state, faster rotational motion was observed, both for ANS and DPH. Our data revealed that dyes get encapsulated more tightly in the gel network when compared to the micellar aggregates. ANS has more restrained rotation compared to DPH. This was attributed to the interaction of the sulfonate group of ANS with water molecules and hydrophilic parts of the gelator molecule. No restricted rotation was observed for DPH in the gel state unlike when it is in the gel phase of lipid bilayer.

Dynamics of pulsed flow in an elastic tube

Internal haemorrhage, often leading to cardio-vascular arrest happens to be one of the prime sources of high fatality rates in mammals. We propose a simplistic model of fluid flow in our attempt to specify the location of the haemorrhagic spot, which, if located accurately, could possibly be operated leading to an instant cure. The model we employ for the purpose is basically fluid mechanical in origin and consists of a viscous fluid, pumped by a periodic force and flowing through an elastic tube. The analogy is with that of blood, pumped from the heart and flowing through an artery or vein. Our results, aided by graphical illustrations, match reasonably well with experimental observations.

Dynamics of water at the interface of a small protein, enterotoxin

Fully atomistic molecular dynamics simulations have been carried out to investigate the correlation of biological activity with dynamics of water molecules in an aqueous protein solution of the toxic domain of enterotoxin (PDB ID: 1ETN). This is a small protein of 13 amino acid residues. Our study of this water soluble protein clearly reveals that water dynamics slows down in the hydration layer. Despite this general slowing down, water molecules in the vicinity of the second beta turn of this protein exhibit faster dynamics than those near other regions of the protein. Since this beta turn is believed to play a critical role in the receptor binding of this protein, the faster dynamics of water near the beta turn m ay have biological significance. The collective orientational dynamics of the water molecules in the protein solution exhibits a characteristic long time component of 27 ps, which agrees well with dielectric relaxation experiments.

Identity, energetics, dynamics and environment of interfacial water molecules in a micellar solution

The structure and energetics of interfacial water molecules in the aqueous micelle of cesium perfluorooctanoate have been investigated, using large-scale atomistic molecular dynamics simulations, with the primary objective of classifying them. The simulations show that the water molecules at the interface fall into two broad classes: bound and free, present in a ratio of 9:1. The bound water molecules can be further categorized on the basis of the number of hydrogen bonds (one or two) that they form with the surfactant headgroups. The hydrogen bonds of the doubly hydrogen-bonded species are found to be, on the average, slightly weaker than those in the singly bonded species. The environment around interfacial water molecules is more ordered than that in the bulk. The surface water molecules have substantially lower potential energy, because of interaction with the micelle. In particular, both forms of bound water have energies that are lower by �2.5-4.0 kcal/ mol. Entropy is found to play an important role in determining the relative concentration of the species.

Study of the dynamics of protein folding through minimalistic models

Starting with the Levinthal paradox, a brief introduction to the protein folding problem is presented. The existing theories of protein folding, including the folding funnel scenario, are discussed. After briefly discussing different simulation studies of model proteins, we discuss our recent work on the dynamics of folding of the model HP-36 (the chicken villin headpiece) protein by using a simplified hydropathy scale. Special attention has been paid to the statics and dynamics of contact formation among the hydrophobic residues. The results obtained from this simple model appear to be surprisingly similar to several features observed in the folding of real proteins. The account concludes with a discussion of future problems.

Vortex dynamics and upper critical fields in ultrathin Bi films

Current-voltage (I-V) characteristics of quench condensed, superconducting, ultrathin Bi films in a magnetic field are reported. These I-V's show hysteresis for all films, grown both with and without thin Ge underlayers. Films on Ge underlayers, close to superconductor-insulator transition, show a peak in the critical current, indicating a structural transformation of the vortex solid. These underlayers, used to make the films more homogeneous, are found to be more effective in pinning the vortices. The upper critical fields (B-c2) of these films are determined from the resistive transitions in perpendicular magnetic field. The temperature dependence of the upper critical field is found to differ significantly from Ginzburg-Landau theory, after modifications for disorder.

Dynamics of pulled desorption with effects of excluded-volume interaction: The p-Laplacian diffusion equation and its exact solution

We analyze the dynamics of desorption of a polymer molecule which is pulled at one of its ends with force f, trying to desorb it. We assume a monomer to desorb when the pulling force on it exceeds a critical value f(c). We formulate an equation for the average position of the n-th monomer, which takes into account excluded-volume interaction through the blob-picture of a polymer under external constraints. The approach leads to a diffusion equation with a p-Laplacian for the propagation of the stretching along the chain. This has to be solved subject to a moving boundary condition. Interestingly, within this approach, the problem can be solved exactly in the trumpet, stem-flower and stem regimes. In the trumpet regime, we get tau = tau(0)n(d)(2), where n(d) is the number of monomers that have desorbed at the time tau. tau(0) is known only numerically, but for f close to f(c), it is found to be tau(0) similar to f(c)/(f(2/3) - f(c)(2/3)) If one used simple Rouse dynamics, this result would change to tau similar to f(c)n(d)(2)/(f - f(c)). In the other regimes too, one can find exact solution, and interestingly, in all regimes tau similar to n(d)(2). Copyright (C) EPLA, 2011

Procedural aspects of modeling the dynamics of a three wheeled vehicle using ADAMS-CAR

We describe in some detail the process of development of a dynamic model of a three wheeled vehicle using ADAMS-CAR. We first describe the rigid body model, and then the modeling of structural flexibilities. The aim of this report is to document procedural details of such modeling, with a view to presenting more research and development oriented investigations in the future. The contents of this report may also be of interest to practicing engineers engaged in multi-body dynamics modeling of wheeled vehicles.

Theory of Gas Hydrates: Effect of the Approximation of Rigid Water Lattice

One of the assumptions of the van der Waals and Platteeuw theory for gas hydrates is that the host water lattice is rigid and not distorted by the presence of guest molecules. In this work, we study the effect of this approximation on the triple-point lines of the gas hydrates. We calculate the triple-point lines of methane and ethane hydrates via Monte Carlo molecular simulations and compare the simulation results with the predictions of van der Waals and Platteeuw theory. Our study shows that even if the exact intermolecular potential between the guest molecules and water is known, the dissociation temperatures predicted by the theory are significantly higher. This has serious implications to the modeling of gas hydrate thermodynamics, and in spite of the several impressive efforts made toward obtaining an accurate description of intermolecular interactions in gas hydrates, the theory will suffer from the problem of robustness if the issue of movement of water molecules is not adequately addressed.

Modified lattice model for mode-I fracture analysis of notched plain concrete beam using probabilistic approach

A modified lattice model using finite element method has been developed to study the mode-I fracture analysis of heterogeneous materials like concrete. In this model, the truss members always join at points where aggregates are located which are modeled as plane stress triangular elements. The truss members are given the properties of cement mortar matrix randomly, so as to represent the randomness of strength in concrete. It is widely accepted that the fracture of concrete structures should not be based on strength criterion alone, but should be coupled with energy criterion. Here, by incorporating the strain softening through a parameter ‘α’, the energy concept is introduced. The softening branch of load-displacement curves was successfully obtained. From the sensitivity study, it was observed that the maximum load of a beam is most sensitive to the tensile strength of mortar. It is seen that by varying the values of properties of mortar according to a normal random distribution, better results can be obtained for load-displacement diagram.