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

Vaporization dynamics of functional droplets in a hot laminar air jet

The vaporization characteristics of pendant droplets of various chemical compositions (like conventional fuels, alternative fuels and nanosuspensions) subjected to convective heating in a laminar air jet have been analyzed. Different heating conditions were achieved by controlling the air temperature and velocity fields around the droplet. A hybrid timescale has been proposed which incorporates the effects of latent heat of vaporization, saturation vapor pressure and thermal diffusivity. This timescale in essence encapsulates the different parameters that influence the droplet vaporization rate. The analysis further permits the evaluation of the effect of various parameters such as surrounding temperature, Reynolds number, far-field vapor presence, impurity content and agglomeration dynamics (nanosuspensions) in the droplet. Flow visualization has been carried out to understand the role of internal recirculation on the vaporization rate. The visualization indicates the presence of a single vortex cell within the droplet on account of the rotation and oscillation of the droplet due to aerodynamic load. External heating induced agglomeration in nanofluids leads to morphological changes during the vaporization process. These morphological changes and alteration in vaporization behavior have been assessed using high speed imaging of the diameter regression and Scanning Electron Microscopy images of the resultant precipitate. (C) 2012 Elsevier Ltd. All rights reserved.

Diversity-Rate Trade-off for Fading Channels with Power Control

We consider a complex, additive, white Gaussian noise channel with flat fading. We study its diversity order vs transmission rate for some known power allocation schemes. The capacity region is divided into three regions. For one power allocation scheme, the diversity order is exponential throughout the capacity region. For selective channel inversion (SCI) scheme, the diversity order is exponential in low and high rate region but polynomial in mid rate region. For fast fading case we also provide a new upper bound on block error probability and a power allocation scheme that minimizes it. The diversity order behaviour of this scheme is same as for SCI but provides lower BER than the other policies.

Force Biased Molecular Dynamics Simulation Study of Effect of Dendrimer Generation on Interaction with DNA

We have studied the effect of dendrimer generation on the interaction between dsDNA and the PAMAM dendrimer using force biased simulation of dsDNA with three generations of dendrimer: G3, G4, and G5. Our results for the potential of mean force (PMF) and the dendrimer asphericity along the binding pathway, combined with visualization of the simulations, demonstrate that dendrimer generation has a pronounced impact on the interaction. The PMF increases linearly with increasing generation of the dendrimer. While, in agreement with previous results, we see an increase in the extent to which the dendrimer bends the dsDNA with increasing dendrimer generation, we also see that the deformation of the dendrimer is greater with smaller generation of the dendrimer. The larger dendrimer forces the dsDNA to conform to its structure, while the smaller dendrimer is forced to conform to the structure of the dsDNA. Monitoring the number of bound cations at different values of force bias distance shows the expected effect of ions being expelled when the dendrimer binds dsDNA.

A computational study of flame dynamics on soot production and transport for an acetylene-air laminar diffusion flame in the field of a convecting line vortex

The present work involves a computational study of soot formation and transport in case of a laminar acetylene diffusion flame perturbed by a convecting line vortex. The topology of the soot contours (as in an earlier experimental work [4]) have been investigated. More soot was produced when vortex was introduced from the air side in comparison to a fuel side vortex. Also the soot topography was more diffused in case of the air side vortex. The computational model was found to be in good agreement with the experimental work [4]. The computational simulation enabled a study of the various parameters affecting soot transport. Temperatures were found to be higher in case of air side vortex as compared to a fuel side vortex. In case of the fuel side vortex, abundance of fuel in the vort ex core resulted in stoichiometrically rich combustion in the vortex core, and more discrete soot topography. Overall soot production too was low. In case of the air side vortex abundance of air in the core resulted in higher temperatures and more soot yield. Statistical techniques like probability density function, correlation coefficient and conditional probability function were introduced to explain the transient dependence of soot yield and transport on various parameters like temperature, a cetylene concentration.

Interface capturing simulations of acoustically driven bubble dynamics in and near tissue

Tissue injury during therapeutic ultrasound or lithotripsy is thought, in cases, to be due to the action of cavitation bubbles. Assessing this and mitigating it is challenging since bubble dynamics in the complex confinement of tissues or in small blood vessels are challenging to predict. Simulations tools require specialized algorithms to simultaneously represent strong acoustic waves and shocks, topologically complex liquid‐vapor phase boundaries, and the complex viscoelastic material dynamics of tissue. We discuss advances in a simulation tool for such situations. A single‐mesh Eulerian solver is used to solve the governing equations. Special sharpening terms maintain the liquid‐vapor interface in face of the finite numerical dissipation included in the scheme to accurately capture shocks. A recent enhancement to this formulation has significantly improved this interface capturing procedure, which is demonstrated for simulation of the Rayleigh collapse of a bubble. The solver also transports elastic stresses and can thus be used to assess the effects of elastic properties on bubble dynamics. A shock‐induced bubble collapse adjacent to a model elastic tissue is used to demonstrate this and draw some conclusions regarding the injury suppressing role that tissue elasticity might play.

Water treatment potential in relation to phytoplankton dynamics in tropical sewage-fed urban water bodies

Urban water bodies frequently receive untreated sewage and water levels in such water bodies are maintained by daily inputs of sewage. They function as “variable-zone” anaerobic-aerobic lagoons suffering several macrophyte, biotic and abiotic stresses. We have studied two such lakes in Bangalore (Bellandur-360 ha and Varthur-220 ha) to understand whether such an occurrence could be made beneficial (maintaining water levels as well as treatment). Such hypertrophic water body receives sewage at 180-250mg/L and is discharged at 25-80mg/L COD/BOD in different seasons. In an earlier study we reported macrophyte altering the purification function of the water body. In this paper we studied the impact of phytoplankton dynamics and macrophyte cover on the functions such as organic load removal. Algal community analysis, algal biomass, macrophyte cover, water quality, nutrient status was studied seasonally during 2009-2010. Oxygen deficiency and sometimes anoxia, recorded from surface samples resulted in high quantities of NH4+-N (30-40mg/L) and phosphate (0.5-4mg/L)-characteristics of anoxic hypertrophic urban lakes. The productiveness favoured high phytoplanktonic community characterized by small cells (<10μm; Chlorella sp. - highest numbers). The lake could be clearly demarcated into an initial anaerobic zone (40% area), a facultative zone (20%) and an aerobic zone (40%) based on redox values and GIS/bathymetry. During summer the lake is covered by floating macrophytes converting the lake into an anoxic/anaerobic water pool subduing the water purification function as well as aesthetics. When macrophytes are controlled such sewage fed water bodies can be used for treating urban wastewater while also maintaining water sustainability in these semi-arid ecosystems. This paper reports the community dynamics of phytoplankton, their function and competition with macrophytes.

Diffusion of pentane isomers in faujasite-type zeolites : NMR and molecular dynamics study

Diffusion of pentane isomers in zeolites NaX has been investigated using pulsed field gradient nuclear magnetic resonance (PFG-NMR) and molecular dynamics (MD) techniques respectively. Temperature and concentration dependence of diffusivities have been studied. The diffusivities obtained from NMR are roughly an order of magnitude smaller than those obtained from MD. The dependence of diffusivity on loading at high temperatures exhibits a type I behavior according to the classification of Karger and Pfeifer 1]. NMR diffusivities of the isomers exhibit the order D(n-pentane) > D(isopentane) > D(neopentane). The results from MD suggest that the diffusivities of the isomers follow the order D(n-pentane) < D(isopentane) < D(neopentane). The activation energies from NMR show E-a(n-pentane) < E-a(isopentane) < E-a(neopentane) whereas those from MD suggest the order E-a(n-pentane) > (isopentane) > E-a(neopentane). The latter follows the predictions of levitation effect whereas those of NMR appears to be due to the presence of defects in the zeolite crystals. The differences between diffusivities estimated by NMR and MD are attributed to the longer time and length scales sampled by the NMR technique, as compared to MD. (C) 2012 Elsevier Inc. All rights reserved.

Microscopic dynamics of nanoparticle monolayers at air-water interface

We present results of surface mechanical and particle tracking measurements of nanoparticles trapped at the air-water interface as a function of their areal density. We monitor both the surface pressure (II) and isothermal compression modulus (epsilon) as well as the dynamics of nanoparticle clusters, using fluorescence confocal microscopy while they are compressed to very high density near the two dimensional close packing density Phi similar to 0.82. We observe non-monotonic variation in both epsilon and the dynamic heterogeneity, characterized by the dynamical susceptibility chi(4) with Phi, in such high density monolayers. We provide insight into the underlying nature of such transitions in close packed high density nanoparticle monolayers in terms of the morphology and flexibility of these soft colloidal particles.. We discuss the significance our results in the context of related studies on two dimensional granular or colloidal systems. (C) 2013 Elsevier Inc. All rights reserved.

Gap-Dependent Quasiparticle Dynamics and Coherent Acoustic Phonons in CaFe2As2 across Spin Density Wave Phase Transition

We report ultrafast quasiparticle (QP) dynamics and coherent acoustic phonons in undoped CaFe2As2 iron pnictide single crystals exhibiting spin-density wave (SDW) and concurrent structural phase transition at temperature T-SDW similar to 165K using femtosecond time-resolved pump-probe spectroscopy. The contributions in transient differential reflectivity arising from exponentially decaying QP relaxation and oscillatory coherent acoustic phonon mode show large variations in the vicinity of T-SDW. From the temperature-dependence of the QP recombination dynamics in the SDW phase, we evaluate a BCS-like temperature dependent charge gap with its zero-temperature value of similar to(1.6 perpendicular to 0.2)k(B)T(SDW), whereas, much above T-SDW, an electron-phonon coupling constant of similar to 0.13 has been estimated from the linear temperature-dependence of the QP relaxation time. The long-wavelength coherent acoustic phonons with typical time-period of similar to 100 ps have been analyzed in the light of propagating strain pulse model providing important results for the optical constants, sounds velocity and the elastic modulus of the crystal in the whole temperature range of 3 to 300 K.

A unified model for the dynamics of driven ribbon with strain and magnetic order parameters

We develop a unified model to explain the dynamics of driven one dimensional ribbon for materials with strain and magnetic order parameters. We show that the model equations in their most general form explain several results on driven magnetostrictive metallic glass ribbons such as the period doubling route to chaos as a function of a dc magnetic field in the presence of a sinusoidal field, the quasiperiodic route to chaos as a function of the sinusoidal field for a fixed dc field, and induced and suppressed chaos in the presence of an additional low amplitude near resonant sinusoidal field. We also investigate the influence of a low amplitude near resonant field on the period doubling route. The model equations also exhibit symmetry restoring crisis with an exponent close to unity. The model can be adopted to explain certain results on magnetoelastic beam and martensitic ribbon under sinusoidal driving conditions. In the latter case, we find interesting dynamics of a periodic one orbit switching between two equivalent wells as a function of an ac magnetic field that eventually makes a direct transition to chaos under resonant driving condition. The model is also applicable to magnetomartensites and materials with two order parameters. (C) 2013 American Institute of Physics. http://dx.doi.org/10.1063/1.4790845]

Nucleation dynamics and pool boiling characteristics of high pressure refrigerant using thermochromic liquid crystals

This paper reports an experimental investigation of oscillating temperature field beneath a single isolated nucleation site using a non-invasive TLC (thermochromic liquid crystal) based thermography technique. Empirical correlations are presented to demonstrate the influence of system pressure and wall heat flux on different ebullition characteristics in the nucleate pool boiling regime of refrigerant R-134a. TLC transient response and two-phase flow structure are captured using synchronized, high resolution imaging. It is observed that the area of influence of nucleation site exhibits a two-part distinct transient behavior during the bubble growth period and broadens to a maximum of 1.57 times the bubble diameter at the instant of bubble departure. This is accompanied by a sharp fall of 2.5 degrees C in the local excess temperature at the nucleation site, which results in momentary augmentation (similar to 40%) in the local heat transfer coefficient at the nucleation origin. The enhanced heat transfer rate observed during the bubble peel-off event is primarily due to transient micro-convection in the wake of the retreating bubble. Further, the results indicate that a slight increase in system pressure from 813.6 to 882.5 kPa has no considerable effect on either the wall superheat or the overall heat transfer coefficient and ebullition frequency. In addition, correlations have been obtained for bubble Reynolds number, Jackob number and the dimensionless bubble generation frequency in terms of modified boiling number.

Predicting the lithospheric stress field and plate motions by joint modeling of lithosphere and mantle dynamics

The way in which basal tractions, associated with mantle convection, couples with the lithosphere is a fundamental problem in geodynamics. A successful lithosphere-mantle coupling model for the Earth will satisfy observations of plate motions, intraplate stresses, and the plate boundary zone deformation. We solve the depth integrated three-dimensional force balance equations in a global finite element model that takes into account effects of both topography and shallow lithosphere structure as well as tractions originating from deeper mantle convection. The contribution from topography and lithosphere structure is estimated by calculating gravitational potential energy differences. The basal tractions are derived from a fully dynamic flow model with both radial and lateral viscosity variations. We simultaneously fit stresses and plate motions in order to delineate a best-fit lithosphere-mantle coupling model. We use both the World Stress Map and the Global Strain Rate Model to constrain the models. We find that a strongly coupled model with a stiff lithosphere and 3-4 orders of lateral viscosity variations in the lithosphere are best able to match the observational constraints. Our predicted deviatoric stresses, which are dominated by contribution from mantle tractions, range between 20-70 MPa. The best-fitting coupled models predict strain rates that are consistent with observations. That is, the intraplate areas are nearly rigid whereas plate boundaries and some other continental deformation zones display high strain rates. Comparison of mantle tractions and surface velocities indicate that in most areas tractions are driving, although in a few regions, including western North America, tractions are resistive. Citation: Ghosh, A., W. E. Holt, and L. M. Wen (2013), Predicting the lithospheric stress field and plate motions by joint modeling of lithosphere and mantle dynamics.

Computational fluid dynamics simulation of open ended thermoacoustic prime mover with straight and tapered resonator

A thermoacoustic refrigerator driven by a thermoacoustic primemover is an effective way to produce durable and long lasting refrigeration due to high reliability, no exotic materials, and no moving parts. Resonator geometry is also one of the important factors that influence the performance of a thermoacoustic prime mover, namely, frequency. Computational fluid dynamics simulation of performance comparison of thermoacoustic prime mover with a straight and tapered resonator is chosen for the present study under an identical stack condition with the air as a working fluid. The frequency and pressure amplitude of oscillations obtained from simulation results were found to be more in the tapered resonator than the straight resonator. Apart from computational fluid dynamics simulation, the simulation studies have also been conducted using design environment for low-amplitude thermoacoustic energy conversion, which predicts the performance of thermoacoustic primemover comparatively well. Simulation results from computational fluid dynamics and design environment for low-amplitude thermoacoustic energy conversion were compared and found to be matching well, representing the good validity of computational fluid dynamics modeling.

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.