Molecular-dynamics studies of systems of confined dumbbell molecules

We present the results of molecular-dynamics simulations of systems of dumbbell molecules confined by parallel molecular walls. We have carried out systematic studies of three cases: freezing, steady flows, and stick-slip friction. We find that the molecular orientational degrees of freedom cause the surface layers to deviate from a planar configuration. Nevertheless, steady flows, in a channel as narrow as 15 molecular sizes, display continuum behavior. A range of mechanisms in the dynamics of the freezing of a confined fluid is found, as a function of the wall-fluid interactions and the bond length of the dumbbell molecules. The simple order-disorder transition associated with stick-slip motion in the presence of a layer of monoatomic lubricant molecules is supplanted by more complex behavior due to rotational degrees of freedom of the diatomic molecules.

Vibration and Stability of Fluid Conveying Pipes with Stochastic Parameters

Flexible cantilever pipes conveying fluids with high velocity are analysed for their dynamic response and stability behaviour. The Young's modulus and mass per unit length of the pipe material have a stochastic distribution. The stochastic fields, that model the fluctuations of Young's modulus and mass density are characterized through their respective means, variances and autocorrelation functions or their equivalent power spectral density functions. The stochastic non self-adjoint partial differential equation is solved for the moments of characteristic values, by treating the point fluctuations to be stochastic perturbations. The second-order statistics of vibration frequencies and mode shapes are obtained. The critical flow velocity is-first evaluated using the averaged eigenvalue equation. Through the eigenvalue equation, the statistics of vibration frequencies are transformed to yield critical flow velocity statistics. Expressions for the bounds of eigenvalues are obtained, which in turn yield the corresponding bounds for critical flow velocities.

Discontinuous shear thickening in confined dilute carbon nanotube suspensions

A monotonic decrease in viscosity with increasing shear stress is a known rheological response to shear flow in complex fluids in general and for flocculated suspensions in particular. Here we demonstrate a discontinuous shear-thickening transition on varying shear stress where the viscosity jumps sharply by four to six orders of magnitude in flocculated suspensions of multiwalled carbon nanotubes (MWNT) at very low weight fractions (approximately 0.5%). Rheooptical observations reveal the shear-thickened state as a percolated structure of MWNT flocs spanning the system size. We present a dynamic phase diagram of the non-Brownian MWNT dispersions revealing a starting jammed state followed by shear-thinning and shear-thickened states. The present study further suggests that the shear-thickened state obtained as a function of shear stress is likely to be a generic feature of fractal clusters under flow, albeit under confinement. An understanding of the shear-thickening phenomena in confined geometries is pertinent for flow-controlled fabrication techniques in enhancing the mechanical strength and transport properties of thin films and wires of nanostructured composites as well as in lubrication issues.

Instability, intermittency, and multiscaling in discrete growth models of kinetic roughening

We show by numerical simulations that discretized versions of commonly studied continuum nonlinear growth equations (such as the Kardar-Parisi-Zhangequation and the Lai-Das Sarma-Villain equation) and related atomistic models of epitaxial growth have a generic instability in which isolated pillars (or grooves) on an otherwise flat interface grow in time when their height (or depth) exceeds a critical value. Depending on the details of the model, the instability found in the discretized version may or may not be present in the truly continuum growth equation, indicating that the behavior of discretized nonlinear growth equations may be very different from that of their continuum counterparts. This instability can be controlled either by the introduction of higher-order nonlinear terms with appropriate coefficients or by restricting the growth of pillars (or grooves) by other means. A number of such ''controlled instability'' models are studied by simulation. For appropriate choice of the parameters used for controlling the instability, these models exhibit intermittent behavior, characterized by multiexponent scaling of height fluctuations, over the time interval during which the instability is active. The behavior found in this regime is very similar to the ''turbulent'' behavior observed in recent simulations of several one- and two-dimensional atomistic models of epitaxial growth.

Nonequilibrium phase transition in the kinetic Ising model: Critical slowing down and the specific-heat singularity

The nonequilibrium dynamic phase transition, in the kinetic Ising model in the presence of an oscillating magnetic field has been studied both by Monte Carlo simulation and by solving numerically the mean-field dynamic equation of motion for the average magnetization. In both cases, the Debye ''relaxation'' behavior of the dynamic order parameter has been observed and the ''relaxation time'' is found to diverge near the dynamic transition point. The Debye relaxation of the dynamic order parameter and the power law divergence of the relaxation time have been obtained from a very approximate solution of the mean-field dynamic equation. The temperature variation of appropriately defined ''specific heat'' is studied by the Monte Carlo simulation near the transition point. The specific heat has been observed to diverge near the dynamic transition point.

Nonequilibrium phase transition in the kinetic Ising model: Divergences of fluctuations and responses near the transition point

The nonequilibrium dynamic phase transition in the kinetic Ising model in the presence of an oscillating magnetic field is studied by Monte Carlo simulation. The fluctuation of the dynamic older parameter is studied as a function of temperature near the dynamic transition point. The temperature variation of appropriately defined ''susceptibility'' is also studied near the dynamic transition point. Similarly, the fluctuation of energy and appropriately defined ''specific heat'' is studied as a function of temperature near the dynamic transition point. In both cases, the fluctuations (of dynamic order parameter and energy) and the corresponding responses diverge (in power law fashion) near the dynamic transition point with similar critical behavior (with identical exponent values).

Polyethylene stibinite phosphate esters: Novel flame-retardant plasticizers for PVC

A new series of multielement flame-retardant plasticizers containing polyethylene stibinite phosphate esters have been prepared by bulk polymerization from ethylene glycol with various antimony (III) aryloxydichlorides and arylphosphorodichloridates possessing various combinations of substituent [Cl,Br,NO2]. All the polymers are pink-coloured viscous fluids. They were characterized by inherent viscosity, density, IR, H-1, C-13 and P-31 NMR spectroscopy. The thermal behaviour of the polymers was compared by thermogravimetric analysis and correlated with their structures. The flammability studies were carried out by the limiting oxygen index test. The polymers containing P, Sb, N and Pr elements in their backbone show superior thermal-and flame-retardant characteristics than the other polymers. A comparative study was carried out with one of the synthesized polymers as a polymeric flame-retardant additive to plasticized PVC. The results showed improved LOI and mechanical properties to that of the conventional flame-retardant additive composition. (C) 1997 Elsevier Science Ltd.

A finite element method for compressible viscous fluid flows

This work presents a mixed three-dimensional finite element formulation for analyzing compressible viscous flows. The formulation is based on the primitive variables velocity, density, temperature and pressure. The goal of this work is to present a `stable' numerical formulation, and, thus, the interpolation functions for the field variables are chosen so as to satisfy the inf-sup conditions. An exact tangent stiffness matrix is derived for the formulation, which ensures a quadratic rate of convergence. The good performance of the proposed strategy is shown in a number of steady-state and transient problems where compressibility effects are important such as high Mach number flows, natural convection, Riemann problems, etc., and also on problems where the fluid can be treated as almost incompressible. Copyright (C) 2010 John Wiley & Sons, Ltd.

Screw dynamo and the generation of nonaxisymmetric magnetic fields

A mechanism is presented here for the amplification of large-scale nonaxisymmetric magnetic fields as a manifestation of the dynamo effect. We generalize a result on restrictions of dynamo actions due to laminar flow originally derived by Zeldovich, Ruzmaikin, and Sokolov [Magnetic Fields in Astrophysics (Gordon and Breach, New York, 1983)]. We show how a screwlike motion having phi and z components of velocity can help to grow a magnetic field. This model postulates a large-scale flow having phi and z components with radial dependences (helical flow). Shear in the radial field, because of a near-flux-freezing condition, causes amplification of the phi component of the magnetic field. The radial and axial components grow due to the presence of turbulent diffusion. The shear in the large scale flow induces an indefinite growth of magnetic field without the a effect; nevertheless, turbulent diffusion forms an important part in the overall mechanism.

Micellar structures of dimeric surfactants with phosphate head groups and wettable spacers: A small-angle neutron scattering study

Dimeric or gemini surfactants consist of two hydrophobic chains and two hydrophilic head groups covalently connected by a hydrophobic or hydrophilic spacer. This paper reports the small-angle neutron scattering (SANS) measurements from aqueous micellar solutions of two different recently developed types of dimeric surfactants: (i) bis-anionic C16H33PO4--(CH2)(m)-PO4-C16H33,2Na(+) dimeric surfactants composed of phosphate head groups and a hydrophobic polymethylene spacer, referred to as 16-m-16,2Na(+), for spacer lengths m = 2, 4, 6, and 10, (ii) bis-cationic C16H33N+(CH3)(2)-CH2-(CH2-O-CH2)(p)-CH2-N+ (CH3)(2)C16H33,2Br(-) dimeric surfactants composed of dimethylammonium head groups and a wettable polyethylene oxide spacer, referred to as 16-CH2-p-CH2-16,2Br(-), for spacer lengths p = 1 - 3. The micellar structures of these surfactants are compared with the earlier studied bis-cationic C16H33N+ (CH3)(2)-(CH2)(m)-N+ (CH3)(2)C16H33,2Br(-) dimeric surfactants composed of dimethylammonium head groups and a hydrophobic polymethylene spacer, referred to as 16-m-16,2Br(-). It is found that 16-m-16,2Na(+), similar to 16-m-16,2Br(-), form various micellar structures depending on the spacer length. Micelles an disklike for rn = 2, rodlike for m = 4, and prolate ellipsoidal fur m = 6 and 10. The micelles of 16-CH2-p-CH2-16,2Br(-) are prolate ellipsoidal for all the values of p = 1 - 3. It is also found that micelles of 16-m-16,2Na(+) and 16-CH2-p-CH2-16,2Br(-) are large in comparison to those of 16-in-16,2Br(-) for similar spacer lengths. This is connected with the fact that both in 16-m-16,2Na(+) and 16-CH2-p-CH2-16,2Br(-), the head group or the spacer is more hydrated as compared to that in the 16-m-16,2Br(-). An increase in the hydration of the spacer or the head group increases the screening of the Coulomb repulsion between the charged head groups. This effect has been found to be more pronounced in the dimeric surfactants having wettable spacers. [S1063-651X(99)00303-7].

Phase separation in a simple model with dynamical asymmetry

We perform computer simulations of a Cahn-Hilliard model of phase separation that has dynamical asymmetry between the two coexisting phases. The dynamical asymmetry is incorporated by considering a mobility function that is order parameter dependent. Simulations of this model reveal morphological features similar to those observed in viscoelastic phase separation. In the early stages, the minority phase domains form a percolating structure that shrinks with time, eventually leading to the formation of disconnected regions that are characterized by the presence of random interfaces as well as isolated droplets. The domains grow as L(t)similar to t(1/3) in the very late stages. Although dynamical scaling is violated in the area shrinking regime, it is restored at late times. However, the form of the scaling function is found to depend on the extent of dynamical asymmetry. [S1063-651X(99)12101-9].

Pronounced hydrogel formation by the self-assembled aggregates of N-alkyl disaccharide amphiphiles

Six disaccharide amphiphiles were synthesized and their hydrogel-forming behavior was extensively studied. These amphiphiles were based on maltose and lactose. Since the gels formed from some of these systems showed the ability to "trap" water molecules upon gelation, these gels were described as "hydrogels". When these gels were heated to similar to 70 degrees C, the samples turned into clear, isotropic fluids, and upon gradual cooling, the hydrogels could be reproduced. Thus these systems were also "thermoreversible". The low molecular mass (MW 565) of the gelators compared to that of a typical polymeric gelator forming substance implies pronounced aggregation of the disaccharide amphiphiles into larger microstructures during gelation. To discern the aggregate textures and morphologies, the specimen hydrogel samples were examined by high-resolution scanning electron microscopy (SEM). A possible reason for the exceptionally high water gelating capacities (>6000 molecules of water per gelator molecule) exhibited by these N-alkyl disaccharide amphiphiles is the presence of large interlamellar spaces into which the water molecules get entrapped due to surface tension. In contrast to their single-chain counterparts, the double-chain lactosyl and maltosylamine amphiphiles upon solubilization in EtOH-H2O afforded hydrogels with reduced mechanical strengths. Interestingly, the corresponding microstructures were found to be quite different from the corresponding hydrogels of their single-chain counterparts. Rheological studies provided further insights into the behavior of these hydrogels. Varying the chain length of the alcohol cosolvent could modulate the gelation capacities, melting temperatures, and the mechanical properties of these hydrogels. To explain the possible reasons of gelation, the results of molecular modeling and energy minimization studies were also included.

Residence time distribution for a class of Gaussian Markov processes

We study the distribution of residence time or equivalently that of "mean magnetization" for a family of Gaussian Markov processes indexed by a positive parameter alpha. The persistence exponent for these processes is simply given by theta=alpha but the residence time distribution is nontrivial. The shape of this distribution undergoes a qualitative change as theta increases, indicating a sharp change in the ergodic properties of the process. We develop two alternate methods to calculate exactly but recursively the moments of the distribution for arbitrary alpha. For some special values of alpha, we obtain closed form expressions of the distribution function. [S1063-651X(99)03306-1].

Subquadratic quantum number dependence and other anomalies of vibrational dephasing in liquid nitrogen: Molecular dynamics simulation study from the triple point to the critical point and beyond

Vibrational phase relaxation near gas-liquid and liquid-solid phase coexistence has been studied by molecular dynamics simulations of N-N stretch in N-2. Experimentally observed pronounced insensitivity of phase relaxation from the triple point to beyond the boiling point is found to originate from a competition between density relaxation and resonant-energy transfer terms. The sharp rise in relaxation rate near the critical point (CP) can be attributed at least partly to the sharp, rise in vibration-rotation coupling contribution. Substantial subquadratic quantum number dependence of overtone dephasing rate is found near the CP and in supercritical fluids. [S0031-9007 (99)09318-7].

Temperature scaling in a dense vibrofluidized granular material

The leading order "temperature" of a dense two-dimensional granular material fluidized by external vibrations is determined. The grain interactions are characterized by inelastic collisions, but the coefficient of restitution is considered to be close to 1, so that the dissipation of energy during a collision is small compared to the average energy of a particle. An asymptotic solution is obtained where the particles are considered to be elastic in the leading approximation. The velocity distribution is a Maxwell-Boltzmann distribution in the leading approximation,. The density profile is determined by solving the momentum balance equation in the vertical direction, where the relation between the pressure and density is provided by the virial equation of state. The temperature is determined by relating the source of energy due to the vibrating surface and the energy dissipation due to inelastic collisions. The predictions of the present analysis show good agreement with simulation results at higher densities where theories for a dilute vibrated granular material, with the pressure-density relation provided by the ideal gas law, sire in error. [:S1063-651X(99)04408-6].