Dynamics of viscous vesicles in shear flow

The dynamics of giant lipid vesicles under shear flow is experimentally investigated. Consistent with previous theoretical and numerical studies, two flow regimes are identified depending on the viscosity ratio between the interior and the exterior of the vesicle, and its reduced volume or excess surface. At low viscosity ratios, a tank-treading motion of the membrane takes place, the vesicle assuming a constant orientation with respect to the flow direction. At higher viscosity ratios, a tumbling motion is observed in which the whole vesicle rotates with a periodically modulated velocity. When the shear rate increases, this tumbling motion becomes increasingly sensitive to vesicle deformation due to the elongational component of the flow and significant deviations from simpler models are observed. A good characterization of these various flow regimes is essential for the validation of analytical and numerical models, and to relate microscopic dynamics to macroscopic rheology of suspensions of deformable particles, such as blood.     

DYNAMIC ANALYSIS OF A ROTATING CANTILEVER BEAM BY USING THE FINITE ELEMENT METHOD

A finite element analysis for a rotating cantilever beam is presented in this study. Based on a dynamic modelling method using the stretch deformation instead of the conventional axial deformation, three linear partial differential equations are derived from Hamilton's principle. Two of the linear differential equations are coupled through the stretch and chordwise deformations. The other equation is an uncoupled one for the flapwise deformation. From these partial differential equations and the associated boundary conditions, are derived two weak forms: one is for the chordwise motion and the other is for the flapwise motion. The weak forms are spatially discretized with newly defined two-node beam elements. With the discretized equations, the behaviours of the natural frequencies are investigated for the variation of the rotating speed. In addition, the time responses and distributions of the deformations and stresses are computed when the rotating speed is prescribed. The effects of the rotating speed profile on the vibrations of the beam are also investigated.     

Transport theory for quantum crystals

A correlation function approach is developed to treat non-equilibrium phenomena of quantum crystals at low frequency and long wavelength within the renormalized harmonic approximation (RHA). The derivation of the transport equations is carried out by studying the hierarchy of equations of motion for the retarded Green's functions of a pure, nonprimitive, nonionic, anharmonic lattice. Using a factorization technique to take into account the most important terms due to the particle fluctuations and the leading contributions to the hydrodynamic singularities of the phonon self-energy, we find a differential equation for the displacement field and a generalized transport equation for the phonon gas. The microscopic RHA expressions for the local temperature, the local heat density and the energy current are derived; the quasiparticle parameters (elastic constants, generalized Grüneisen parameters, quasiparticle interaction) entering the equations of motion are shown to be consistent with the RHA. In the hydrodynamic regime the general equations are reduced to two coupled differential equations for the lattice deformations and for the local temperature. Then only the displacement-displacement, the displacement-energy density and the energy density-energy density correlation functions show macroscopic fluctuations; for these functions thermodynamical sum-rules are derived.     

New Results for Directed Vesicles and Chains near an Attractive Wall

In this paper we present new exact results for single fully directed walks and fully directed vesicles near an attractive wall. This involves a novel method of solution for these types of problems. The major advantage of this method is that it, unlike many other single-walker methods, generalizes to an arbitrary number of walkers. The method of solution involves solving a set of partial difference equations with a Bethe Ansatz. The solution is expressed as a “constant-term” formula which evaluates to sums of products of binomial coefficients. The vesicle critical temperature is found at which a binding transition takes place, and the asymptotic forms of the associated partition functions are found to have three different entropic exponents depending on whether the temperature is above, below, or at its critical value. The expected number of monomers adsorbed onto the surface is found to become proportional to the vesicle length at temperatures below critical. Scaling functions near the critical point are determined.     

An elastic,plastic, viscous model for slow shear of a liquid foam

We suggest a scalar model for deformation and flow of an amorphous material such as a foam or an emulsion. To describe elastic, plastic and viscous behaviours, we use three scalar variables: elastic deformation, plastic deformation rate and total deformation rate; and three material-specific parameters: shear modulus, yield deformation and viscosity. We obtain equations valid for different types of deformations and flows slower than the relaxation rate towards mechanical equilibrium. In particular, they are valid both in transient or steady flow regimes, even at large elastic deformation. We discuss why viscosity can be relevant even in this slow shear (often called “quasi-static”) limit. Predictions of the storage and loss moduli agree with the experimental literature, and explain with simple arguments the non-linear large amplitude trends.     

基于离散变分原理的耗散动力学模拟方法：模拟三维囊泡形状

将基于离散变分原理的耗散动力学模拟方法应用到三维囊泡体系,通过优化囊泡的弯曲能求解其平衡态形状. 该方法的优点之一是不需要预先假定对称性. 针对特定约化自发曲率的囊泡体系,该方法模拟获得了一系列轴对称形状,模拟结果与文献中预先假定轴对称条件的计算方法所报道的结果一致,这验证了该模拟方法的可靠性及精确性. 此外,使用该方法研究了两个差别巨大的平衡态形状之间的转变动力学,在转变过程中观察到了多个非轴对称的中间形状. 研究结果表明该方法不仅可以模拟囊泡的非轴对称结构,而且具备模拟囊泡在剧烈形变下演化过程的能力. 为研究更复杂的囊泡体系,特别是生物膜的形变提供了一个重要的理论模拟方法. 关键词： 生物膜 离散空间变分法 耗散动力学 三角网格划分     

A boundary integral method for simulating the dynamics of inextensible vesicles suspended in a viscous fluid in 2D

We present a new method for the evolution of inextensible vesicles immersed in a Stokesian fluid. We use a boundary integral formulation for the fluid that results in a set of nonlinear integro-differential equations for the vesicle dynamics. The motion of the vesicles is determined by balancing the non-local hydrodynamic forces with the elastic forces due to bending and tension. Numerical simulations of such vesicle motions are quite challenging. On one hand, explicit time-stepping schemes suffer from a severe stability constraint due to the stiffness related to high-order spatial derivatives and a milder constraint due to a transport-like stability condition. On the other hand, an implicit scheme can be expensive because it requires the solution of a set of nonlinear equations at each time step. We present two semi-implicit schemes that circumvent the severe stability constraints on the time step and whose computational cost per time step is comparable to that of an explicit scheme. We discretize the equations by using a spectral method in space, and a multistep third-order accurate scheme in time. We use the fast multipole method (FMM) to efficiently compute vesicle–vesicle interaction forces in a suspension with a large number of vesicles. We report results from numerical experiments that demonstrate the convergence and algorithmic complexity properties of our scheme.     

Buckling of single layer graphene sheet based on nonlocal elasticity and higher order shear deformation theory

Higher order shear deformation theory (HSDT) is reformulated using the nonlocal differential constitutive relations of Eringen. The equations of motion of the nonlocal theories are derived. The developed equations of motion have been applied to study buckling characteristics of nanoplates such as graphene sheets. Navier's approach has been used to solve the governing equations for all edges simply supported boundary conditions. Analytical solutions for critical buckling loads of the graphene sheets are presented. Nonlocal elasticity theories are employed to bring out the small scale effect on the critical buckling load of graphene sheets. Effects of (i) nonlocal parameter, (ii) length, (iii) thickness of the graphene sheets and (iv) higher order shear deformation theory on the critical buckling load have been investigated. The theoretical development as well as numerical solutions presented herein should serve as reference for nonlocal theories as applied to the stability analysis of nanoplates and nanoshells.     

Thermodynamical Bending of FGM Sandwich Plates Resting on Pasternak's Elastic Foundations

The analysis of thermoelastic deformations of a simply supported functionally graded material (FGM) sandwich plates subjected to a time harmonic sinusoidal temperature field on the top surface and varying through-the-thickness is illustrated in this paper. The FGM sandwich plates are assumed to be made of three layers and resting on Pasternak's elastic foundations. The volume fractions of the constituents of the upper and lower layers and, hence, the effective material properties of them are assumed to vary in the thickness direction only whereas the core layer is still homogeneous. When in-plane sinusoidal variations of the displacements and the temperature that identically satisfy the boundary conditions at the edges, the governing equations of motion are solved analytically by using various shear deformation theories as well as the classical one. The influences of the time parameter, power law index, temperature exponent, top-to-bottom surface temperature ratio, side-to-thickness ratio and the foundation parameters on the dynamic bending are investigated.     

Vacillating breathing and tumbling of vesicles under shear flow

The dynamics of vesicles under a shear flow are analyzed analytically in the small deformation regime. We derive two coupled nonlinear equations which describe the vesicle orientation in the flow and its shape evolution. A new type of motion is found, namely, a "vacillating-breathing" mode: the vesicle orientation undergoes an oscillation around the flow direction, while the shape executes breathing dynamics. This solution coexists with tumbling. Moreover, we provide an explicit expression for the tumbling threshold. A rheological law for a dilute vesicle suspension is outlined.     

Deformation parameter changes in fission mass yields within the systematic statistical scission-point model

The fission fragment mass-yields are evaluated for pre-actinide and actinide isotopes using a systematic statistical scission point model. The total potential energy of the fissioning systems at the scission point is presented in approximate relations as functions of mass numbers,deformation parameters and the temperature of complementary fission fragments. The collective temperature, Tcoll, and the temperature of fission fragments, Ti, are separated and the effect of collective temperature on mass yields results is investigated. The fragment temperature has been calculated with the generalized superfluid model. The sum of deformation parameters of complementary fission fragments has been obtained by fitting the calculated results with the experimental data. To investigate the transitions between symmetric and asymmetric modes mass yields for pre-actinide and heavy actinides are calculated with this model. The transition from asymmetric to symmetric fission is well reproduced using this systematic statistical scission point model. The calculated results are in good agreement with the experimental data with Tcoll= 2 Me V at intermediate excitation energy and with T_(coll)= 1MeV for spontaneous fission.Despite the Langevin model, in the scission point model, a constraint on the deformation parameters of fission fragments has little effect on the results of the mass yield.     

Effect of temperature and geometrical nonlinearity on wave propagation characteristics of carbon nanotubes

Considering the effect of temperature and geometrical nonlinearity in the constitutive relation, the equation of motion for a carbon nanotube is obtained based on the Euler–Bernouli beam model. Also, the effect of van der Waals forces is taken into account in the formulation. The carbon nanotube is assumed to be under the application of a constant distributed external load. At any temperature, the equilibrium solutions of the governing equations for a single-walled carbon nanotube (SWCNT) and a double-walled carbon nanotube (DWCNT) are obtained. A small perturbation is assumed around the equilibrium solution. Using this perturbation, the nonlinear equations of motion are linearized. Using the linearized form of the equations of motion, the characteristic equations and dispersion relations are obtained. It is shown that in the linear case and for the case of high temperature there exists a temperature beyond which the phase velocity does not exist. It is shown that in the case of room or low temperature there is no critical value for temperature. Based on the dispersion equation, a relation for the critical value of temperature is obtained. It is found that when the large deformation effect is taken into account, the critical value for temperature does not exist. Also, the effect of large deformations on phase velocities and lateral deformations of single-walled and double-walled carbon nanotube beams are studied. It is found that unlike the linear theory, the nonlinear theory predicts a non-zero phase velocity at the temperature corresponding to linear critical temperature.     

Macroscopic equations of motion and dynamic polarizability in the study of nuclear collective excitations

Time dependent Hartree-Fock equations are derived using a variational principle over a restricted part of the space of the Slater determinants, in the limit of small deformations (RPA). When an external oscillating field interacts with the nucleus, this method leads to an explicit expression of the nuclear response function (dynamic polarizability) as a function of the external frequency and of the deformation field, defining the nuclear deformation induced by the interaction. A linear differential equation for the deformation field is also obtained: in the limit ω → ∞ it has analytical solutions which satisfy the energy-weighted sum rule, evaluated in a HF ground state, in both isoscalar and isovector modes.     

Lifting the singular nature of a model for peeling of an adhesive tape

We investigate the dynamics of peeling of an adhesive tape subjected to a constant pull speed. Due to the constraint between the pull force, peel angle and the peel force, the equations of motion derived earlier fall into the category of differential-algebraic equations (DAE) requiring an appropriate algorithm for its numerical solution. By including the kinetic energy arising from the stretched part of the tape in the Lagrangian, we derive equations of motion that support stick-slip jumps as a natural consequence of the inherent dynamics itself, thus circumventing the need to use any special algorithm. In the low mass limit, these equations reproduce solutions obtained using a differential-algebraic algorithm introduced for the earlier singular equations. We find that mass has a strong influence on the dynamics of the model rendering periodic solutions to chaotic and vice versa. Apart from the rich dynamics, the model reproduces several qualitative features of the different waveforms of the peel force function as also the decreasing nature of force drop magnitudes.     

Molecular dynamic simulations of elastomer structure and its influence on anisotropic stress under time-varying strain

We investigate the evolution of polymer structure and its influence on uniaxial anisotropic stress under time-varying uniaxial strain, and the role of external control variables such as temperature, strain rate, chain length, and density, using molecular dynamics simulation. At temperatures higher than glass transition, stress anisotropy in the system is reduced even though the bond stretch is greater at higher temperatures. There is a significant increase in the stress level with increasing density. At higher densities, the uncoiling of the chains is suppressed and the major contribution to the deformation is by internal deformation of the chains. At faster rates of loading stress anisotropy increases. The deformation mechanism is mostly due to bond stretch and bond bending rather than overall shape and size. Stress levels increase with longer chain length. There is a critical value of the functionality of the cross-linkers beyond which the uniaxial stress developed increases caused primarily by bond stretching due to increased constraint on the motion of the monomers. Stacking of the chains in the system also plays a dominant role in the behaviour in terms of excluded volume interactions. Low density, high temperature, low values of functionality of cross-linkers, and short chain length facilitate chain uncoiling and chain slipping in cross-linked polymers.     

Application of a fusion model using collective variables and microscopically related mass parameters

A calculation of nucleus-nucleus collisions is presented, using a model which starts from a TDHF equation and leads to classical equations of motion for a set of four collective variables. Restricting to axial symmetry and assuming the liquid drop mass formula to hold, a differential equation is derived, which describes nuclear deformations and energies and is used to construct a potential energy surface for the collective variables. The nuclear deformations are obtained without the need of shape parameters. The equations of motion for the collective variables are solved numerically.     

Motor proteins transporting cargos

Processive motor proteins such as kinesin and myosin-V are enzymes that use the energy of ATP hydrolysis to travel along polar cytoskeletal filaments. One of the functions of these proteins is the transport of vesicles and protein complexes that are linked to the light chains of the motors. Modeling the light chain by a linear elastic spring, and using the two-state model for one- and two-headed molecular motors, we study the influence of thermal fluctuations of the cargo on the motion of the motor-cargo complex. We solve numerically the Fokker-Planck equations of motor motion, and find that the mean velocity of the motor-cargo complex decreases monotonously as the spring becomes softer. This effect is due to the random force of thermal fluctuations of the cargo disrupting the operation of the motor. Increasing the size (thus, the friction coefficient) of the cargo also decreases the velocity. Surprisingly, we find that for a given size of the cargo, the velocity has a maximum for a certain friction of the motor. We explain this effect by the interplay between the characteristic length of thermal fluctuations of the cargo on a spring, the motor diffusion length, and the filament period. Our results may be relevant for the interpretation of single-molecule experiments with molecular motors (bead assays), where the motor motion is observed by tracking of a bead attached to the motor.     

The equations of motion of initially stressed Timoshenko tubular beams conveying fluid

In this paper the equations of motion of an initially stressed Timoshenko tubular beam subjected to a tensile follower load and conveying fluid are derived by using the appropriate statement of Hamilton's principle. This latter is obtained first for “open” systems, the instantaneous total mass of which does not necessarily remain constant in the course of deformation—“open” denoting that there is momentum transport in and out of the system. The equations of motion are derived separately for a cantilevered system and for one with both extremities of the tube clamped. Yet another derivation for the cantilevered tube is presented with the system considered to be quasi-closed, where all flow-induced effects are incorporated through the virtual work, as if they were “external” forces. All three sets of equations are found to be identical. These equations are then compared with those obtained, more simply, by the Newtonian force-balance approach. Some differences between them are found to exist, the principal of which are associated with the follower or other tensile forces; these are discussed at some length, and the equations of motion obtained here are compared to those obtained by other researchers for Timoshenko beams subjected to follower or tensile forces.     

A modified digital phase-shift moiré technique for impact deformation measurements

Due to the short loading times and high deformation rates inherent to impact experiments, measurement of the occurring deformations is not straightforward. Presented in this paper is a technique to obtain the displacements and deformations of a specimen subjected to a uni-axial impact load. During the experiment the deformation of a line grating attached to the specimen is captured using a streak camera. From the recorded deforming grating the specimen displacements are automatically derived using an advanced numerical algorithm, based on the interference between the specimen grating and a virtual reference grating. Numerical interference is considered because it allows that the pitch of the reference grating is adapted to the changing amplitude of the deformation. Indeed, at each moment of the deformation process the pitch of the reference grating is chosen such that the highest possible accuracy and sensitivity is guaranteed. Because of this, large changes in deformation amplitude are allowed, and the technique is applicable to a wide range of materials. Eventual imperfections of the specimen grating and temperature effects are taken into account. Specimen displacements are extracted automatically by means of a phase-shifting technique.

The non-contact measurement technique yields high resolution, quantitative information on the specimen deformation, along the entire length of the specimen and during the full duration of the experiment. Interaction by the operator is excluded. Results are presented of a high strain rate tensile test on a steel sheet specimen showing local deformations up to about 170%.