Fluid Mechanics of Aquatic Locomotion at Large Reynolds Numbers

Abstract | There exist a huge range of fish species besides other aquatic organisms like squids and salps that locomote in water at large Reynolds numbers, a regime of flow where inertial forces dominate viscous forces. In the present review, we discuss the fluid mechanics governing the locomotion of such organisms. Most fishes propel themselves by periodic undulatory motions of the body and tail, and the typical classification of their swimming modes is based on the fraction of their body that undergoes such undulatory motions. In the angulliform mode, or the eel type, the entire body undergoes undulatory motions in the form of a travelling wave that goes from head to tail, while in the other extreme case, the thunniform mode, only the rear tail (caudal fin) undergoes lateral oscillations. The thunniform mode of swimming is essentially based on the lift force generated by the airfoil like crosssection of the fish tail as it moves laterally through the water, while the anguilliform mode may be understood using the “reactive theory” of Lighthill. In pulsed jet propulsion, adopted by squids and salps, there are two components to the thrust; the first due to the familiar ejection of momentum and the other due to an over-pressure at the exit plane caused by the unsteadiness of the jet. The flow immediately downstream of the body in all three modes consists of vortex rings; the differentiating point being the vastly different orientations of the vortex rings. However, since all the bodies are self-propelling, the thrust force must be equal to the drag force (at steady speed), implying no net force on the body, and hence the wake or flow downstream must be momentumless. For such bodies, where there is no net force, it is difficult to directly define a propulsion efficiency, although it is possible to use some other very different measures like “cost of transportation” to broadly judge performance.

Free Energy Barriers for Escape of Water Molecules from Protein Hydration Layer

Free energy barriers separating interfacial water molecules from the hydration layer at the surface of a protein to the bulk are obtained by using the umbrella sampling method of free energy calculation. We consider hydration layer of chicken villin head piece (HP-36) which has been studied extensively by molecular dynamics simulations. The free energy calculations reveal a strong sensitivity to the secondary structure. In particular, we find a region near the junction of first and second helix that contains a cluster of water molecules which are slow in motion, characterized by long residence times (of the order of 100 ps or more) and separated by a large free energy barrier from the bulk water. However, these ``slow'' water molecules constitute only about 5-10% of the total number of hydration layer water molecules. Nevertheless, they play an important role in stabilizing the protein conformation. Water molecules near the third helix (which is the important helix for biological function) are enthalpically least stable and exhibit the fastest dynamics. Interestingly, barrier height distributions of interfacial water are quite broad for water surrounding all the three helices (and the three coils), with the smallest barriers found for those near the helix-3. For the quasi-bound water molecules near the first and second helices, we use well-known Kramers' theory to estimate the residence time from the free energy surface, by estimating the friction along the reaction coordinate from the diffusion coefficient by using Einstein relation. The agreement found is satisfactory. We discuss the possible biological function of these slow, quasi-bound (but transient) water molecules on the surface.

Structural Transformations, Composition Anomalies and a Dramatic Collapse of Linear Polymer Chains in Dilute Ethanol-Water Mixtures

Water-ethanol mixtures exhibit many interesting anomalies, such as negative excess partial molar volume of ethanol, excess sound absorption coefficient at low concentrations, and positive deviation from Raoult's law for vapor pressure, to mention a few. These anomalies have been attributed to different, often contradictory origins, but a quantitative understanding is still lacking. We show by computer simulation and theoretical analyses that these anomalies arise from the sudden emergence of a bicontinuous phase that occurs at a relatively low ethanol concentration of x(eth) approximate to 0.06-0.10 (that amounts to a volume fraction of 0.17-0.26, which is a significant range!). The bicontinuous phase is formed by aggregation of ethanol molecules, resulting in a weak phase transition whose nature is elucidated. We find that the microheterogeneous structure of the mixture gives rise to a pronounced nonmonotonic composition dependence of local compressibility and nonmonotonic dependence in the peak value of the radial distribution function of ethyl groups. A multidimensional free energy surface of pair association is shown to provide a molecular explanation of the known negative excess partial volume of ethanol in terms of parallel orientation and hence better packing of the ethyl groups in the mixture due to hydrophobic interactions. The energy distribution of the ethanol molecules indicates additional energy decay channels that explain the excess sound attenuation coefficient in aqueous alcohol mixtures. We studied the dependence of the solvation of a linear polymer chain on the composition of the water-ethanol solvent. We find that there is a sudden collapse of the polymer at x(eth) approximate to 0.05-a phenomenon which we attribute to the formation of the microheterogeneous structures in the binary mixture at low ethanol concentrations. Together with recent single molecule pulling experiments, these results provide new insight into the behavior of polymer chain and foreign solutes, such as enzymes, in aqueous binary mixtures.

Relaxation and jump dynamics of water at the mica interface

The orientational relaxation dynamics of water confined between mica surfaces is investigated using molecular dynamics simulations. The study illustrates the wide heterogeneity that exists in the dynamics of water adjacent to a strongly hydrophilic surface such as mica. Analysis of the survival probabilities in different layers is carried out by normalizing the corresponding relaxation times with bulk water layers of similar thickness. A 10-fold increase in the survival times is observed for water directly in contact with the mica surface and a non-monotonic variation in the survival times is observed moving away from the mica surface to the bulk-like interior. The orientational relaxation time is highest for water in the contact layer, decreasing monotonically away from the surface. In all cases the ratio of the relaxation times of the 1st and 2nd rank Legendre polynomials of the HH bond vector is found to lie between 1.5 and 1.9 indicating that the reorientational relaxation in the different water layers is governed by jump dynamics. The orientational dynamics of water in the contact layer is particularly novel and is found to undergo distinct two-dimensional hydrogen bond jump reorientational dynamics with an average waiting time of 4.97 ps. The waiting time distribution is found to possess a long tail extending beyond 15 ps. Unlike previously observed jump dynamics in bulk water and other surfaces, jump events in the mica contact layer occur between hydrogen bonds formed by the water molecule and acceptor oxygens on the mica surface. Despite slowing down of the water orientational relaxation near the surface, life-times of water in the hydration shell of the K ion are comparable to that observed in bulk salt solutions. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4717710]

Analysis of Parameter Values in the van der Waals and Platteeuw Theory for Methane Hydrates Using Monte Carlo Molecular Simulations

The van der Waals and Platteuw (vdVVP) theory has been successfully used to model the thermodynamics of gas hydrates. However, earlier studies have shown that this could be due to the presence of a large number of adjustable parameters whose values are obtained through regression with experimental data. To test this assertion, we carry out a systematic and rigorous study of the performance of various models of vdWP theory that have been proposed over the years. The hydrate phase equilibrium data used for this study is obtained from Monte Carlo molecular simulations of methane hydrates. The parameters of the vdWP theory are regressed from this equilibrium data and compared with their true values obtained directly from simulations. This comparison reveals that (i) methane-water interactions beyond the first cage and methane-methane interactions make a significant contribution to the partition function and thus cannot be neglected, (ii) the rigorous Monte Carlo integration should be used to evaluate the Langmuir constant instead of the spherical smoothed cell approximation, (iii) the parameter values describing the methane-water interactions cannot be correctly regressed from the equilibrium data using the vdVVP theory in its present form, (iv) the regressed empty hydrate property values closely match their true values irrespective of the level of rigor in the theory, and (v) the flexibility of the water lattice forming the hydrate phase needs to be incorporated in the vdWP theory. Since methane is among the simplest of hydrate forming molecules, the conclusions from this study should also hold true for more complicated hydrate guest molecules.

A Series of Trifacial Pd6 Molecular Barrels with Porphyrin Walls

Three new nanoscopic trigonal prisms, (tmen)6Pd6(H2L)3](NO3)12 (1), (Meen)6Pd6(H2L)3](NO3)12 (2), and (2,2'-bipy)6Pd6(H2L)3](NO3)12 (3), have been synthesized in excellent yields through single-step metalligand-coordination-driven self-assembly using 5,10,15,20-tetrakis(3-pyridyl)porphyrin (H2L) as a donor and cis-blocked PdII 90 degrees acceptors. These complexes were fully characterized by spectroscopic studies and single-crystal X-ray diffraction. All of these barrels quantitatively bind ZnII ions in the N4 pockets of the porphyrin walls at room temperature. Their corresponding zinc-embedded complexes, (tmen)6Pd6(ZnL)3](NO3)12 (1?a), (Meen)6Pd6(ZnL)3](NO3)12 (2?a), and (2,2'-bipy)6Pd6(ZnL)3](NO3)12 (3?a), were synthesized under ambient conditions by the post-synthetic binding of ZnII ions into the H2N4 pockets of the porphyrin walls of these complexes. These zinc-embedded complexes were characterized by electronic absorption, fluorescence emission, 1H NMR spectroscopy, as well as elemental analysis. Complexes 13 exhibited considerable microporosity in their solid state. Complex 1 was an efficient adsorbent for nitrogen gas and EtOH, MeOH, and water vapors.

Thermal vibration analysis of orthotropic nanoplates based on nonlocal continuum mechanics

This paper presents the thermal vibration analysis of orthotropic nanoplates such as graphene, using the two variable refined plate theory and nonlocal continuum mechanics for small scale effects. The nanoplate is modeled based on two variable refined plate theory and the axial stress caused by the thermal effects is also considered. The two variable refined plate theory takes account of transverse shear effects and parabolic distribution of the transverse shear strains through the thickness of the plate, hence it is unnecessary to use shear correction factors. Nonlocal governing equations of motion for the nanoplate are derived from the principle of virtual displacements. The closed form solution for thermal-vibration frequencies of a simply supported rectangular nanoplate has been obtained by using Navier's method of solution. Numerical results obtained by the present theory are compared with available solutions in the literature and the molecular dynamics results. The influences of the small scale coefficient, the room or low temperature, the high temparature, the half wave number and the aspect ratio of nanoplate on the natural frequencies are considered and discussed in detail. It can be concluded that the present theory, which does not require shear correction factor, is not only simple but also comparable to the first-order and higher order shear deformation theory. The present analysis results can be used for the design of the next generation of nanodevices that make use of the thermal vibration properties of the nanoplates. (C) 2012 Elsevier B.V. All rights reserved.

Role of Specific Cations and Water Entropy on the Stability of Branched DNA Motif Structures

DNA three-way junctions (TWJs) are important intermediates in various cellular processes and are the simplest of a family of branched nucleic acids being considered as scaffolds for biomolecular nanotechnology. Branched nucleic acids are stabilized by divalent cations such as Mg2+, presumably due to condensation and neutralization of the negatively charged DNA backbone. However, electrostatic screening effects point to more complex solvation dynamics and a large role of interfacial waters in thermodynamic stability. Here, we report extensive computer simulations in explicit water and salt on a model TWJ and use free energy calculations to quantify the role of ionic character and strength on stability. We find that enthalpic stabilization of the first and second hydration shells by Mg2+ accounts for 1/3 and all of the free energy gain in 50% and pure MgCl2 solutions, respectively. The more distorted DNA molecule is actually destabilized in pure MgCl2 compared to pure NaCl. Notably, the first shell, interfacial waters have very low translational and rotational entropy (i.e., mobility) compared to the bulk, an entropic loss that is overcompensated by increased enthalpy from additional electrostatic interactions with Mg2+. In contrast, the second hydration shell has anomalously high entropy as it is trapped between an immobile and bulklike layer. The nonmonotonic entropic signature and long-range perturbations of the hydration shells to Mg2+ may have implications in the molecular recognition of these motifs. For example, we find that low salt stabilizes the parallel configuration of the three-way junction, whereas at normal salt we find antiparallel configurations deduced from the NMR. We use the 2PT analysis to follow the thermodynamics of this transition and find that the free energy barrier is dominated by entropic effects that result from the decreased surface area of the antiparallel form which has a smaller number of low entropy waters in the first monolayer.

Nonlocal continuum mechanics based ultrasonic flexural wave dispersion characteristics of a monolayer graphene embedded in polymer matrix

Wave propagation in graphene sheet embedded in elastic medium (polymer matrix) has been a topic of great interest in nanomechanics of graphene sheets, where the equivalent continuum models are widely used. In this manuscript, we examined this issue by incorporating the nonlocal theory into the classical plate model. The influence of the nonlocal scale effects has been investigated in detail. The results are qualitatively different from those obtained based on the local/classical plate theory and thus, are important for the development of monolayer graphene-based nanodevices. In the present work, the graphene sheet is modeled as an isotropic plate of one-atom thick. The chemical bonds are assumed to be formed between the graphene sheet and the elastic medium. The polymer matrix is described by a Pasternak foundation model, which accounts for both normal pressure and the transverse shear deformation of the surrounding elastic medium. When the shear effects are neglected, the model reduces to Winkler foundation model. The normal pressure or Winkler elastic foundation parameter is approximated as a series of closely spaced, mutually independent, vertical linear elastic springs where the foundation modulus is assumed equivalent to stiffness of the springs. For this model, the nonlocal governing differential equations of motion are derived from the minimization of the total potential energy of the entire system. An ultrasonic type of flexural wave propagation model is also derived and the results of the wave dispersion analysis are shown for both local and nonlocal elasticity calculations. From this analysis we show that the elastic matrix highly affects the flexural wave mode and it rapidly increases the frequency band gap of flexural mode. The flexural wavenumbers obtained from nonlocal elasticity calculations are higher than the local elasticity calculations. The corresponding wave group speeds are smaller in nonlocal calculation as compared to local elasticity calculation. The effect of y-directional wavenumber (eta(q)) on the spectrum and dispersion relations of the graphene embedded in polymer matrix is also observed. We also show that the cut-off frequencies of flexural wave mode depends not only on the y-direction wavenumber but also on nonlocal scaling parameter (e(0)a). The effect of eta(q) and e(0)a on the cut-off frequency variation is also captured for the cases of with and without elastic matrix effect. For a given nanostructure, nonlocal small scale coefficient can be obtained by matching the results from molecular dynamics (MD) simulations and the nonlocal elasticity calculations. At that value of the nonlocal scale coefficient, the waves will propagate in the nanostructure at that cut-off frequency. In the present paper, different values of e(0)a are used. One can get the exact e(0)a for a given graphene sheet by matching the MD simulation results of graphene with the results presented in this article. (c) 2012 Elsevier Ltd. All rights reserved.

Effect of pressure on the ionic conductivity of Li+ and Cl- ions in water

A molecular dynamics simulation study of aqueous solution of LiCl is reported as a function of pressure. Experimental measurements of conductivity of Li+ ion as a function of pressure shows an increase in conductivity with pressure. Our simulations are able to reproduce the observed trend in conductivity. A number of relevant properties have been computed in order to understand the reasons for the increase in conductivity with pressure. These include radial distribution function, void and neck distributions, hydration or coordination numbers, diffusivity, velocity autocorrelation functions, angles between ion-oxygen and dipole of water as well as OH vector, mean residence time for water in the hydration shell, etc. These show that the increase in pressure acts as a structure breaker. The decay of the self part of the intermediate scattering function at small wave number k shows a bi-exponential decay at 1 bar which changes to single exponential decay at higher pressures. The k dependence of the ratio of the self part of the full width at half maximum of the dynamic structure factor to 2Dk(2) exhibits trends which suggest that the void structure of water is playing a role. These support the view that the changes in void and neck distributions in water can account for changes in conductivity or diffusivity of Li+ with pressure. These results can be understood in terms of the levitation effect. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4756909]

Simulation of Influence of Bilayer Melting on Dynamics and Thermodynamics of Interfacial Water

We investigate the effect of bilayer melting transition on thermodynamics and dynamics of interfacial water using molecular dynamics simulation with the two-phase thermodynamic model. We show that the diffusivity of interface water depicts a dynamic crossover at the chain melting transition following an Arrhenius behavior until the transition temperature. The corresponding change in the diffusion coefficient from the bulk to the interface water is comparable with experimental observations found recently for water near 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) vesicles Phys. Chem. Chem. Phys. 13, 7732 (2011)]. The entropy and potential energy of interfacial water show distinct changes at the bilayer melting transition, indicating a strong correlation in the thermodynamic state of water and the accompanying first-order phase transition of the bilayer membrane. DOI: 10.1103/PhysRevLett.110.018303

Thermal vibration analysis of monolayer graphene embedded in elastic medium based on nonlocal continuum mechanics

This paper presents the thermal vibration analysis of single-layer graphene sheet embedded in polymer elastic medium, using the plate theory and nonlocal continuum mechanics for small scale effects. The graphene is modeled based on continuum plate theory and axial stress caused by the thermal effects is also considered. Nonlocal governing equation of motion for this graphene sheet system is derived from the principle of virtual displacements. The closed form solution for thermal-vibration frequencies of a simply supported rectangular nanoplate has been obtained by using the Navier's method of solution. Numerical results obtained by the present theory are compared with available solutions in the literature and the molecular dynamics results. The influences of the small scale coefficient, the room or low temperature, the high temperature, the half wave number and the aspect ratio of nanoplate on the natural frequencies are considered and discussed in detail. The thermal vibration analysis of single- and double-layer graphene sheets are considered for the analysis. The mode shapes of the respective graphene system are also captured in this work. The present analysis results can be used for the design of the next generation of nanodevices that make use of the thermal vibration properties of the graphene.

Structure, molecular dynamics, and stress in a linear polymer

We present a study correlating uniaxial stress in a polymer with its underlying structure when it is strained. The uniaxial stress is significantly influenced by the mean-square bond length and mean bond angle. In contrast, the size and shape of the polymer, typically represented by the end-to-end length, mass ratio, and radius of gyration, contribute negligibly. Among externally set control variables, density and polymer chain length play a critical role in influencing the anisotropic uniaxial stress. Short chain polymers more or less behave like rigid molecules. Temperature and rate of loading, in the range considered, have a very mild effect on the uniaxial stress.

Ionomer based blend as water vapor barrier material for organic device encapsulation

Blends of poly (ethylene-co-methacrylic acid) (PEMA) and poly (vinyl alcohol-co-ethylene) (EVOH) were studied for encapsulating Schottky structured organic devices. A calcium degradation test was used to determine water vapor transmission rates and to determine the moisture barrier performance of neat and blend films. Moisture barrier analysis for the neat and blend compositions was discussed concerning the interactions in the blend, diffusivity of water molecules through the unit cell systems, and the occupiable free volumes available in the unit cells using molecular dynamics simulations. The experimental results of water vapor permeation were correlated with diffusion behavior predicted from molecular dynamics simulations results. The effectiveness of the blend as a suitable barrier material in increasing the lifetime of an encapsulated Schottky structured organic device was determined.

EFFECT OF CORE-SHELL STRUCTURE OF HYDROGEL BEADS ON THE THRESHOLD CONCENTRATION OF WATER FOR SWELLING AND ITS PH SENSITIVITY

In this paper we investigate the effect of core-shell structure of Sodium Alginate based hydrogel beads and their size on certain activation threshold concentration of water for applications in swelling and pH sensing. This type of hydrogel experiences diffusive pressure due to transport of certain free charges across its interface with a solvent or electrolyte. This process is essentially a dynamic equilibrium of the electric force field, stress in the polymeric network with cage like structure and molecular diffusion including phase transformation due to pressure imbalance between the hydrogel and its surroundings. The effect of pH of the solvant on the swelling rate of these beads has been studied experimentally. A mathematical model of the swelling process has been developed by considering Nernst-Planck equation representing the migration of mobile ions and Er ions, Poisson equation representing the equilibrium of the electric field and mechanical field equation representing swelling of the gel. An attempt has been made to predict the experimentally observed phenomena using these numerical simulations. It is observed experimentally that certain minimum concentration called activation threshold concentration of the water molecules must be present in the hydrogel in order to activate the swelling process. For the required activation threshold concentration of water in the beads, the pH induced change in the rate of swelling is also investigated. This effect is analyzed for various different core-shell structures of the beads.