Simulations of impinging droplets with surfactant-dependent dynamic contact angle

An arbitrary Lagrangian-Eulerian (ALE) finite element scheme for computations of soluble surfactant droplet impingement on a horizontal surface is presented. The numerical scheme solves the time-dependent Navier-Stokes equations for the fluid flow, scalar convection-diffusion equation for the surfactant transport in the bulk phase, and simultaneously, surface evolution equations for the surfactants on the free surface and on the liquid-solid interface. The effects of surfactants on the flow dynamics are included into the model through the surface tension and surfactant-dependent dynamic contact angle. In particular, the dynamic contact angle (theta(d)) of the droplet is defined as a function of the surfactant concentration at the contact line and the equilibrium contact angle (theta(0)(e)) of the clean surface using the nonlinear equation of state for surface tension. Further, the surface forces are included into the model as surface divergence of the surface stress tensor that allows to incorporate the Marangoni effects without calculating the surface gradient of the surfactant concentration on the free surface. In addition to a mesh convergence study and validation of the numerical results with experiments, the effects of adsorption and desorption surfactant coefficients on the flow dynamics in wetting, partially wetting and non-wetting droplets are studied in detail. It is observed that the effects of surfactants are more in wetting droplets than in the non-wetting droplets. Further, the presence of surfactants at the contact line reduces the equilibrium contact angle further when theta(0)(e) is less than 90 degrees, and increases it further when theta(0)(e) is greater than 90 degrees. Nevertheless, the presence of surfactants has no effect on the contact angle when theta(0)(e) = 90 degrees. The numerical study clearly demonstrates that the surfactant-dependent contact angle has to be considered, in addition to the Marangoni effect, in order to study the flow dynamics and the equilibrium states of surfactant droplet impingement accurately. The proposed numerical scheme guarantees the conservation of fluid mass and of the surfactant mass accurately. (C) 2015 Elsevier Inc. All rights reserved.

Estimation of activation energy for electroporation and pore growth rate in liquid crystalline and gel phases of lipid bilayers using molecular dynamics simulations

Molecular dynamics simulations of electroporation in POPC and DPPC lipid bilayers have been carried out at different temperatures ranging from 230 K to 350 K for varying electric fields. The dynamics of pore formation, including threshold field, pore initiation time, pore growth rate, and pore closure rate after the field is switched off, was studied in both the gel and liquid crystalline (L-alpha) phases of the bilayers. Using an Arrhenius model of pore initiation kinetics, the activation energy for pore opening was estimated to be 25.6 kJ mol(-1) and 32.6 kJ mol(-1) in the L-alpha phase of POPC and DPPC lipids respectively at a field strength of 0.32 V nm(-1). The activation energy decreases to 24.2 kJ mol(-1) and 23.7 kJ mol(-1) respectively at a higher field strength of 1.1 V nm(-1). At temperatures below the melting point, the activation energy in the gel phase of POPC and DPPC increases to 28.8 kJ mol(-1) and 34.4 kJ mol(-1) respectively at the same field of 1.1 V nm(-1). The pore closing time was found to be higher in the gel than in the L-alpha phase. The pore growth rate increases linearly with temperature and quadratically with field, consistent with viscosity limited growth models.

Notch sensitivity in nanoscale metallic glass specimens: Insights from continuum simulations

Recent experiments have shown that nano-sized metallic glass (MG) specimens subjected to tensile loading exhibit increased ductility and work hardening. Failure occurs by necking as opposed to shear banding which is seen in bulk samples. Also, the necking is generally observed at shallow notches present on the specimen surface. In this work, continuum finite element analysis of tensile loading of nano-sized notched MG specimens is conducted using a thermodynamically consistent non-local plasticity model to clearly understand the deformation behavior from a mechanics perspective. It is found that plastic zone size in front of the notch attains a saturation level at the stage when a dominant shear band forms extending across the specimen. This size scales with an intrinsic material length associated with the interaction stress between flow defects. A transition in deformation behavior from quasi-brittle to ductile becomes possible when this critical plastic zone size is larger than the uncracked ligament length. These observations corroborate with atomistic simulations and experimental results. (C) 2015 Elsevier Ltd. All rights reserved.

Ab Initio MD Simulations of the Bronsted Acidity of Glutathione in Aqueous Solutions: Predicting pK(a) Shifts of the Cysteine Residue

The tripeptide glutathione (GSH) is one of the most abundant peptides and the major repository for nonprotein sulfur in both animal and plant cells. It plays a critical role in intracellular oxidative stress management by the reversible formation of glutathione disulfide with the thiol-disulfide pair acting as a redox buffer. The state of charge of the ionizable groups of GSH can influence the redox couple, and hence the pK(a) value of the cysteine residue of GSH is critical to its functioning. Here we report ab initio Car-Parrinello molecular dynamics simulations of glutathione solvated by 200 water molecules, all of which are considered in the simulation. We show that the free-energy landscape for the protonation-deprotonation reaction of the cysteine residue of GSH computed using metadynamics sampling provides shift in the dissociation constant values as compared with the isolated accurate estimates of the pK(a) and correctly predicts the cysteine amino acid.

Meridional gradients in aerosol vertical distribution over Indian Mainland: Observations and model simulations

Multi-year observations from the network of ground-based observatories (ARFINET), established under the project `Aerosol Radiative Forcing over India' (ARFI) of Indian Space Research Organization and space-borne lidar `Cloud Aerosol Lidar with Orthogonal Polarization' (CALIOP) along with simulations from the chemical transport model `Goddard Chemistry Aerosol Radiation and Transport' (GOCART), are used to characterize the vertical distribution of atmospheric aerosols over the Indian landmass and its spatial structure. While the vertical distribution of aerosol extinction showed higher values close to the surface followed by a gradual decrease at increasing altitudes, a strong meridional increase is observed in the vertical spread of aerosols across the Indian region in all seasons. It emerges that the strong thermal convections cause deepening of the atmospheric boundary layer, which although reduces the aerosol concentration at lower altitudes, enhances the concentration at higher elevations by pumping up more aerosols from below and also helping the lofted particles to reach higher levels in the atmosphere. Aerosol depolarization ratios derived from CALIPSO as well as the GOCART simulations indicate the dominance of mineral dust aerosols during spring and summer and anthropogenic aerosols in winter. During summer monsoon, though heavy rainfall associated with the Indian monsoon removes large amounts of aerosols, the prevailing southwesterly winds advect more marine aerosols over to landmass (from the adjoining oceans) leading to increase in aerosol loading at lower altitudes than in spring. During spring and summer months, aerosol loading is found to be significant, even at altitudes as high as 4 km, and this is proposed to have significant impacts on the regional climate systems such as Indian monsoon. (C) 2015 Elsevier Ltd. All rights reserved.

Pattern formation in thin films of polymer solutions: Theory and simulations

Self-assembly has been recognized as an efficient tool for generating a wide range of functional, chemically, or physically textured surfaces for applications in small scale devices. In this work, we investigate the stability of thin films of polymer solutions. For low concentrations of polymer in the solution, long length scale dewetting patterns are obtained with wavelength approximately few microns. Whereas, for concentrations above a critical value, bimodal dispersion curves are obtained with the dominant wavelength being up to two orders smaller than the usual dewetting length scale. We further show that the short wavelength corresponds to the phase separation in the film resulting in uniformly distributed high and low concentration regions. Interestingly, due to the solvent entropy, at very high concentration values of polymer, a re-entrant behaviour is observed with the dominant length scale now again corresponding to the dewetting wavelength. Thus, we show that the binary films of polymer solutions provide additional control parameters that can be utilized for generating functional textured surfaces for various applications. (C) 2016 AIP Publishing LLC.

Anthropogenic aerosol fraction over the Indian region: model simulations versus multi-satellite data analysis

Anthropogenic aerosols play a crucial role in our environment, climate, and health. Assessment of spatial and temporal variation in anthropogenic aerosols is essential to determine their impact. Aerosols are of natural and anthropogenic origin and together constitute a composite aerosol system. Information about either component needs elimination of the other from the composite aerosol system. In the present work we estimated the anthropogenic aerosol fraction (AF) over the Indian region following two different approaches and inter-compared the estimates. We espouse multi-satellite data analysis and model simulations (using the CHIMERE Chemical transport model) to derive natural aerosol distribution, which was subsequently used to estimate AF over the Indian subcontinent. These two approaches are significantly different from each other. Natural aerosol satellite-derived information was extracted in terms of optical depth while model simulations yielded mass concentration. Anthropogenic aerosol fraction distribution was studied over two periods in 2008: premonsoon (March-May) and winter (November-February) in regard to the known distinct seasonality in aerosol loading and type over the Indian region. Although both techniques have derived the same property, considerable differences were noted in temporal and spatial distribution. Satellite retrieval of AF showed maximum values during the pre-monsoon and summer months while lowest values were observed in winter. On the other hand, model simulations showed the highest concentration of AF in winter and the lowest during pre-monsoon and summer months. Both techniques provided an annual average AF of comparable magnitude (similar to 0.43 +/- 0.06 from the satellite and similar to 0.48 +/- 0.19 from the model). For winter months the model-estimated AF was similar to 0.62 +/- 0.09, significantly higher than that (0.39 +/- 0.05) estimated from the satellite, while during pre-monsoon months satellite-estimated AF was similar to 0.46 +/- 0.06 and the model simulation estimation similar to 0.53 +/- 0.14. Preliminary results from this work indicate that model-simulated results are nearer to the actual variation as compared to satellite estimation in view of general seasonal variation in aerosol concentrations.

Molecular dynamics and simulations study on the vibrational and electronic solvatochromism of benzophenone

Solvent plays a key role in diverse physico-chemical and biological processes. Therefore, understanding solute-solvent interactions at the molecular level of detail is of utmost importance. A comprehensive solvatochromic analysis of benzophenone (Bzp) was carried out in various solvents using Raman and electronic spectroscopy, in conjunction with Density Functional Theory (DFT) calculations of supramolecular solute-solvent clusters generated using classical Molecular Dynamics Simulations (c-MDSs). The >C=O stretching frequency undergoes a bathochromic shift with solvent polarity. Interestingly, in protic solvents this peak appears as a doublet: c-MDS and ad hoc explicit solvent ab initio calculations suggest that the lower and higher frequency peaks are associated with the hydrogen bonded and dangling carbonyl group of Bzp, respectively. Additionally, the dangling carbonyl in methanol (MeOH) solvent is 4 cm(-1) blue-shifted relative to acetonitrile solvent, despite their similar dipolarity/polarizability. This suggests that the cybotactic region of the dangling carbonyl group in MeOH is very different from its bulk solvent structure. Therefore, we propose that this blue-shift of the dangling carbonyl originates in the hydrophobic solvation shell around it resulting from extended hydrogen bonding network of the protic solvents. Furthermore, the 1(1)n pi* (band I) and 1(1)pi pi* (band II) electronic transitions show a hypsochromic and bathochromic shift, respectively. In particular, these shifts in protic solvents are due to differences in their excited state-hydrogen bonding mechanisms. Additionally, a linear relationship is obtained for band I and the >C=O stretching frequency (cm(-1)), which suggests that the different excitation wavelengths in band I correspond to different solvation states. Therefore, we hypothesize that the variation in excitation wavelengths in band I could arise from different solvation states leading to varying solvation dynamics. This will have implications for ultrafast processes associated with electron-transfer, charge transfer, and also the photophysical aspects of excited states. (C) 2016 AIP Publishing LLC.

Depletion of nonlinearity in magnetohydrodynamic turbulence: Insights from analysis and simulations

It is shown how suitably scaled, order-m moments, D-m(+/-), of the Elsasser vorticity fields in three-dimensional magnetohydrodynamics (MHD) can be used to identify three possible regimes for solutions of the MHD equations with magnetic Prandtl number P-M = 1. These vorticity fields are defined by omega(+/-) = curl z(+/-) = omega +/- j, where z(+/-) are Elsasser variables, and where omega and j are, respectively, the fluid vorticity and current density. This study follows recent developments in the study of three-dimensional Navier-Stokes fluid turbulence Gibbon et al., Nonlinearity 27, 2605 (2014)]. Our mathematical results are then compared with those from a variety of direct numerical simulations, which demonstrate that all solutions that have been investigated remain in only one of these regimes which has depleted nonlinearity. The exponents q(+/-) that characterize the inertial range power-law dependencies of the z(+/-) energy spectra, epsilon(+/-)(k), are then examined, and bounds are obtained. Comments are also made on (a) the generalization of our results to the case P-M not equal 1 and (b) the relation between D-m(+/-) and the order-m moments of gradients of magnetohydrodynamic fields, which are used to characterize intermittency in turbulent flows.

Cold gas in cluster cores: global stability analysis and non-linear simulations of thermal instability

We perform global linear stability analysis and idealized numerical simulations in global thermal balance to understand the condensation of cold gas from hot/virial atmospheres (coronae), in particular the intracluster medium (ICM). We pay particular attention to geometry (e.g. spherical versus plane-parallel) and the nature of the gravitational potential. Global linear analysis gives a similar value for the fastest growing thermal instability modes in spherical and Cartesian geometries. Simulations and observations suggest that cooling in haloes critically depends on the ratio of the cooling time to the free-fall time (t(cool)/t(ff)). Extended cold gas condenses out of the ICM only if this ratio is smaller than a threshold value close to 10. Previous works highlighted the difference between the nature of cold gas condensation in spherical and plane-parallel atmospheres; namely, cold gas condensation appeared easier in spherical atmospheres. This apparent difference due to geometry arises because the previous plane-parallel simulations focused on in situ condensation of multiphase gas but spherical simulations studied condensation anywhere in the box. Unlike previous claims, our non-linear simulations show that there are only minor differences in cold gas condensation, either in situ or anywhere, for different geometries. The amount of cold gas depends on the shape of tcool/tff; gas has more time to condense if gravitational acceleration decreases towards the centre. In our idealized plane-parallel simulations with heating balancing cooling in each layer, there can be significant mass/energy/momentum transfer across layers that can trigger condensation and drive tcool/tff far beyond the critical value close to 10.

FEM Simulations and Optimization about Residual Stresses in Coating Structures with Functionally Graded Materials Layer

An elasto-plastic finite element method is developed to predict the residual stresses of thermal spraying coatings with functionally graded material layer. In numerical simulations, temperature sensitivity of various material constants is included and mix

Formation of Fivefold Deformation Twins in Nanocrystalline Face-Centered-Cubic Copper Based on Molecular Dynamics Simulations

Fivefold deformation twins were reported recently to be observed in the experiment of the nanocrystalline face-centered-cubic metals and alloys. However, they were not predicted previously based on the molecular dynamics (MD) simulations and the reason was thought to be a uniaxial tension considered in the simulations. In the present investigation, through introducing pretwins in grain regions, using the MD simulations, the authors predict out the fivefold deformation twins in the grain regions of the nanocrystal grain cell, which undergoes a uniaxial tension. It is shown in their simulation results that series of Shockley partial dislocations emitted from grain boundaries provide sequential twining mechanism, which results in fivefold deformation twins. (c) 2006 American Institute of Physics.

The Uniaxial Tensile Deformation of Ni Nanowire: Atomic-Scale Computer Simulations

The mechanical deformations of nickel nanowire subjected to uniaxial tensile strain at 300 K are simulated by using molecular dynamics with the quantum corrected Sutten-Chen many-body force field. We have used common neighbor analysis method to investigate the structural evolution of Ni nanowire during the elongation process. For the strain rate of 0.1%/ps, the elastic limit is up to about 11% strain with the yield stress of 8.6 GPa. At the elastic stage, the deformation is carried mainly through the uniform elongation of the distances between the layers (perpendicular to the Z-axis) while the atomic structure remains basically unchanged. With further strain, the slips in the {111} planes start to take place in order to accommodate the applied strain to carry the deformation partially, and subsequently the neck forms. The atomic rearrangements in the neck region result in a zigzag change in the stress-strain curve; the atomic structures beyond the region, however, have no significant changes. With the strain close to the point of the breaking, we observe the formation of a one-atom thick necklace in Ni nanowire. The strain rates have no significant effect on the deformation mechanism, but have some influence on the yield stress, the elastic limit, and the fracture strain of the nanowire.