Experiments and Lagrangian simulations on the formation of droplets in continuous mode

Experimental and computational studies on the dynamics of millimeter-scale cylindrical liquid jets are presented. The influences of the modulation amplitude and the nozzle geometry on jet behavior have been considered. Laser Doppler anemometry (LDA) was used in order to extract the velocity field of a jet along its length, and to determine the velocity modulation amplitude. Jet shapes and breakup dynamics were observed via shadowgraph imaging. Aqueous solutions of glycerol were used for these experiments. Results were compared with Lagrangian finite-element simulations with good quantitative agreement. © 2011 The American Physical Society.

A Molecular Dynamics Simulation: the Effect of Finite Size on the Thermal Conductivity in a Single Crystal Silicon

Non-equilibrium molecular dynamics (NEMD) simulations are performed to calculate thermal conductivity. The environment-dependent interatomic potential (EDIP) potential on crystal silicon is adopted as a model system. The issues are related to nonlinear response, local thermal equilibrium and statistical averaging. The simulation results by non-equilibrium molecular dynamics show that the calculated thermal conductivity decreases almost linearly as the film thickness reduced at the nanometre scale. The effect of size on the thermal conductivity is also obtained by a theoretic analysis of the kinetic theory and formulas of the heat capacity. The analysis reveals that the contributions of phonon mean free path (MFP) and phonon number in a finite cell to thermal conductivity are very important.

A Constrained Particle Dynamics for Continuum-Particle Hybrid Method in Micro- and Nano-Fluidics

A hybrid method of continuum and particle dynamics is developed for micro- and nano-fluidics, where fluids are described by a molecular dynamics (MD) in one domain and by the Navier-Stokes (NS) equations in another domain. In order to ensure the continuity of momentum flux, the continuum and molecular dynamics in the overlap domain are coupled through a constrained particle dynamics. The constrained particle dynamics is constructed with a virtual damping force and a virtual added mass force. The sudden-start Couette flows with either non-slip or slip boundary condition are used to test the hybrid method. It is shown that the results obtained are quantitatively in agreement with the analytical solutions under the non-slip boundary conditions and the full MD simulations under the slip boundary conditions.

Size-Dependent Elastic Properties Of Ni Nanofilms By Molecular Dynamics Simulation

Size-dependent elastic properties of Ni nanofilms are investigated by molecular dynamics ( MD) simulations with embedded atom method (EAM). The surface effects are considered by calculating the surface relaxation, surface energy, and surface stress. The Young's modulus and yield stress are obtained as functions of thickness and crystallographic orientation. It is shown that the surface relaxation has important effects on the the elastic properties at nanoscale. When the surface relaxation is outward, the Young's modulus decreases with the film thickness decreasing, and vice versa. The results also show that the yield stresses of the films increase with the films becoming thinner. With the thickness of the nanofilms decreasing, the surface effects on the elastic properties become dominant.

Molecular Dynamics Simulation Of Adhesion Forces In A Dipalmitoylphosphatidylcholine Membrane

Adhesion forces of Dipalmitoylphosphatidylcholine ( DPPC) membrane in the gel phase are investigated by molecular dynamics ( MD) simulation. In the simulations, individual DPPC molecules are pulled out of DPPC membranes with different rates and we get the maximum adhesion forces of DPPC membrane. We find that the maximum adhesion forces increase with pull rate, from about 400 to 700 pN when pull rates are from 0.001 to 0.03 nm/ps. We analyze the relationship between pull rate and adhesion forces of different origins using Brownian dynamics and notice that viscosity of solvent plays an important role in adhesion forces. Then we simulate the motion of a single DPPC molecule in solvent and it elucidates that the maximum drag force is almost linear with respect to the pull rate. We use Stokes' relation to describe the motion of a single DPPC molecule and deduce the effective length of a DPPC molecule. Conformational analyses indicate that the free energy variation of a DPPC molecule inside and outside of the DPPC membrane is an essential part of adhesion energy.

Strengthening Mechanisms And Dislocation Dynamics In Twinned Metal Nanowires

Microtwins are frequently observed in face-centered-cubic (fcc) metal nanowires with low stacking fault energy. The authors have previously reported that the tensile Yield strength of copper nanowires can be increased by, the presence of twin boundaries. lit this work, simulations are carried out under both uniaxial tension and compression loading, to demonstrate that the strengthening effects are inherent to these nanowires, independent of the loading condition (tensile/compressive). It appears that the strengthening mechanism of the twinned nanowires can be attributed to stress redistribution due to the change of crystallographic orientations across twin boundaries, which requires larger external stress to make them Yield as compared to the twin-free wire.

Study Of Mechanical Properties Of Amorphous Copper With Molecular Dynamics Simulation

The formation and mechanical properties of amorphous copper are studied using molecular dynamics simulation. The simulations of tension and shearing show that more pronounced plasticity is found under shearing, compared to tension. Apparent strain hardening and strain rate effect are observed. Interestingly, the variations of number density of atoms during deformation indicate free volume creation, especially under higher strain rate. In particular, it is found that shear induced dilatation does appear in the amorphous metal.

A Note On The Numerical Simulations Of Flow Past A Wavy Square-Section Cylinder

The flow past a square-section cylinder with a geometric disturbance is investigated by numerical simulations. The extra terms, due to the introduction of mapping transformation simulating the effect of disturbance into the transformed Navier-Stokes equations, are correctly derived, and the incorrect ones in the previous literature are pointed out and analyzed. Furthermore, the relationship between the vorticity, especially on the cylinder surface, and the disturbance is derived and explained theoretically. The computations are performed at two Reynolds numbers of 100 and 180 and three amplitudes of waviness of 0.006, 0.025 and 0.167 with another aim to explore the effects of different Reynolds numbers and disturbance on the vortex dynamics in the wake and forces on the body. Numerical results have shown that, at the mild waviness of 0.025, the Karman vortex shedding is suppressed completely for Re = 100, while the forced vortex dislocation is appeared in the near wake at the Reynolds number of 180. The drag reduction is up to 21.6% at Re = 100 and 25.7% at Re = 180 for the high waviness of 0.167 compared with the non-wavy cylinder. The lift and the Strouhal number varied with different Reynolds numbers and the wave steepness are also obtained.

Adsorption of His-tagged peptide to Ni, Cu and Au (100) surfaces: Molecular dynamics simulation

Molecular dynamics (MD) simulations are performed to study the interaction of His-tagged peptide with three different metal surfaces in explicit water. The equilibrium properties are analyzed by using pair correlation functions (PCF) to give an insight into the behavior of the peptide adsorption to metal surfaces in water solvent. The intermolecular interactions between peptide residues and the metal surfaces are evaluated. By pulling the peptide away from the peptide in the presence of solvent water, peeling forces are obtained and reveal the binding strength of peptide adsorption on nickel, copper and gold. From the analysis of the dynamics properties of the peptide interaction with the metal surfaces, it is shown that the affinity of peptide to Ni surface is the strongest, while on Cu and An the affinity is a little weaker. In MD simulations including metals, the His-tagged region interacts with the substrate to an extent greater than the other regions. The work presented here reveals various interactions between His-tagged peptide and Ni/Cu/Au surfaces. The interesting affinities and dynamical properties of the peptide are also derived. The results give predictions for the structure of His-tagged peptide adsorbing on three different metal surfaces and show the different affinities between them, which assist the understanding of how peptides behave on metal surfaces and of how designers select amino sequences in molecule devices design. (c) 2007 Elsevier Ltd. All rights reserved.

Molecular dynamics studies on the local disordering of σ3 and σ11 grain boundaries in aluminium bicrystal

Molecular dynamics (MD) simulations using Morse interaction potential are performed in studies of [110] symmetrical tilt grain boundary (GB) structures with mis-orientation angles 50.5 degrees(Sigma 11), 129.5 degrees(Sigma 11), 70.5 degrees(Sigma 3) and 109.5 degrees(Sigma 3) at various tempratures. The GB structures are found to start local disordering at about 0.5T(m)(T-m is the melting point of aluminium) for 50.5 degrees(Sigma 11), 0.32T(m) for 129.5 degrees(Sigma 11) and 0.38T(m) for 70.5 degrees(Sigma 3), respectively. These results agree with conclusions deduced from the anelastic measurements. But, for twin-boundary structure 109.5 degrees(Sigma 3), this disordering has not been found even when temperature increases up to 0.9T(m).

Molecular Dynamics Simulation of Barnacle Cement

Barnacle cement is an underwater adhesive that is used for permanent settlement. Its main components are insoluble protein complexes that have not been fully studied. In present article, we chose two proteins of barnacle cement for study, 36-KD protein and Mrcp-100K protein. In order to investigate the characteristic of above two proteins, we introduced the method of molecular modeling. And the simulation package GROMACS was used to simulate the behavior of these proteins. In this article, before the simulations, we introduce some theories to predict the time scale for polymer relaxation. During the simulation, we mainly focus on two properties of these two proteins: structural stability and adhesive force to substrate. First, we simulate the structural stability of two proteins in water, and then the stability of 36-KD protein in seawater environment is investigated. We find that the stability varies in the different environments. Next, to study adhesive ability of two proteins, we simulate the process of peeling the two proteins from the substrate (graphite). Then, we analyze the main reasons of these results. We find that hydrogen bonds in proteins play an important role in the protein stability. In the process of the peeling, we use Lennard-Jones 12-6 potential to calculate the van der Waals interactions between proteins and substrate.

Size-dependent elastic properties of Ni nanofilms by molecular dynamics simulation

Size-dependent elastic properties of Ni nanofilms are investigated by molecular dynamics ( MD) simulations with embedded atom method (EAM). The surface effects are considered by calculating the surface relaxation, surface energy, and surface stress. The Young's modulus and yield stress are obtained as functions of thickness and crystallographic orientation. It is shown that the surface relaxation has important effects on the the elastic properties at nanoscale. When the surface relaxation is outward, the Young's modulus decreases with the film thickness decreasing, and vice versa. The results also show that the yield stresses of the films increase with the films becoming thinner. With the thickness of the nanofilms decreasing, the surface effects on the elastic properties become dominant.

Molecular dynamics simulation of adhesion forces in a dipalmitoylphosphatidylcholine membrane

Adhesion forces of Dipalmitoylphosphatidylcholine ( DPPC) membrane in the gel phase are investigated by molecular dynamics ( MD) simulation. In the simulations, individual DPPC molecules are pulled out of DPPC membranes with different rates and we get the maximum adhesion forces of DPPC membrane. We find that the maximum adhesion forces increase with pull rate, from about 400 to 700 pN when pull rates are from 0.001 to 0.03 nm/ps. We analyze the relationship between pull rate and adhesion forces of different origins using Brownian dynamics and notice that viscosity of solvent plays an important role in adhesion forces. Then we simulate the motion of a single DPPC molecule in solvent and it elucidates that the maximum drag force is almost linear with respect to the pull rate. We use Stokes' relation to describe the motion of a single DPPC molecule and deduce the effective length of a DPPC molecule. Conformational analyses indicate that the free energy variation of a DPPC molecule inside and outside of the DPPC membrane is an essential part of adhesion energy.

Quantum simulations of relativistic quantum physics in circuit QED

We present a scheme for simulating relativistic quantum physics in circuit quantum electrodynamics. By using three classical microwave drives, we show that a superconducting qubit strongly coupled to a resonator field mode can be used to simulate the dynamics of the Dirac equation and Klein paradox in all regimes. Using the same setup we also propose the implementation of the Foldy-Wouthuysen canonical transformation, after which the time derivative of the position operator becomes a constant of the motion.