Parallelization of dissipative particle dynamics simulation

The dissipative particle dynamics simulation is usually used to study polymer in mesoscopic space. The traditional methods are resource intensive, especially when the scale of research is large. Therefore, improving computing efficiency is a key point in this research area. Two major issues are addressed in this paper. First, the DPD methods are analysed and the most time-consuming parts are identified: conservative force, dissipative force and random force. Second, we describe how to parallelize the existing serial application in the Windows Compute Cluster Server (WCCS) platform. The results show that the parallel algorithm not only effectively shortens the computing time, but also improves the resource utilization rate.

Equilibration strategies in molecular dynamics simulation of polyethers

A comparison of the NVT and NVE ensemble simulations of tetraglyme showed that, in terms of energy, temperature and most of the structural features the results were very similar. However, major differences were observed in dynamic properties, ie in the mean square displacement and in the O--O distances. A fast equilibration method suitable for amorphous polymer systems is also detailed. This was accomplished by the reassignment of the velocity distribution to the particles, after the interruption of an NVE simulation.

Many-body dissipative particle dynamics simulation of wetting phenomena

With the development of the simulation of particle dynamics, the traditional dissipative particle dynamics (DPD) method can not satisfy the needs of research in static or dynamic wetting phenomena. However, the Many-body DPD approach extends the ability of the traditional method to simulate the interface between solid and liquid or some other situation. In this paper, we propose a Many-body DPD program to simulate the solid-liquid interface and get satisfactory results.

Study of growth mechanism of TiC cluster in ferrite via molecular dynamics simulation

Abstract Atomistic simulations were used to investigate the evolution process of titanium carbide clusters to mature precipitates in ferrite. The typical kinetic of carbide cluster growth was studied in detail through analyzing the atomic interactions of a carbide cluster with scattered carbon atoms. The driving force required for cluster growth was calculated along with the atomic diffusivity in the iron matrix, exploring the change in response as two main growth steps. The growth kinetic improved the understanding of precipitate evolution at the atomic level.

A molecular dynamics simulation of alloy carbide clusters formation

 The formation of alloy carbide cluster in ferrite was investigated via molecular dynamics simulation, which disclosed the cluster property and formation mechanism. These together provided a better fundamental understanding of the cluster formation and firm information for the evolution of cluster and precipitate in high-strength low-alloy steel.

CISR-ODE, A C plus plus framework with ODE solver for code based system dynamics simulation

Ordinary differential equations are used for modelling a wide range of dynamic systems. Even though there are many graphical software applications for this purpose, a fully customised solution for all problems is code-level programming of the model and solver. In this project, a free and open source C++ framework is designed to facilitate modelling in native code environment and fulfill the common simulation needs of control and many other engineering and science applications. The solvers of this project are obtained from ODEINT and specialised for Armadillo matrix library to provide an easy syntax and a fast execution. The solver code is minimised and its modification for users have become easier. There are several features added to the solvers such as controlling maximum step size, informing the solver about sudden input change and forcing custom times into the results and calling a custom method at these points. The comfort of the model designer, code readability, extendibility and model isolation have been considered in the structure of this framework. The application manages the output results, exporting and plotting them. Modifying the model has become more practical and a portion of corresponding codes are updated automatically. A set of libraries is provided for generation of output figures, matrix hashing, control system functions, profiling, etc. In this paper, an example of using this framework for a classical washout filter model is explained.

A full body musculoskeletal model based on flexible multibody simulation approach utilised in bone strain analysis during human locomotion

Load-induced strains applied to bone can stimulate its development and adaptation. In order to quantify the incident strains within the skeleton, in vivo implementation of strain gauges on the surfaces of bone is typically used. However, in vivo strain measurements require invasive methodology that is challenging and limited to certain regions of superficial bones only such as the anterior surface of the tibia. Based on our previous study [Al Nazer et al. (2008) J Biomech. 41:1036–1043], an alternative numerical approach to analyse in vivo strains based on the flexible multibody simulation approach was proposed. The purpose of this study was to extend the idea of using the flexible multibody approach in the analysis of bone strains during physical activity through integrating the magnetic resonance imaging (MRI) technique within the framework. In order to investigate the reliability and validity of the proposed approach, a three-dimensional full body musculoskeletal model with a flexible tibia was used as a demonstration example. The model was used in a forward dynamics simulation in order to predict the tibial strains during walking on a level exercise. The flexible tibial model was developed using the actual geometry of human tibia, which was obtained from three-dimensional reconstruction of MRI. Motion capture data obtained from walking at constant velocity were used to drive the model during the inverse dynamics simulation in order to teach the muscles to reproduce the motion in the forward dynamics simulation. Based on the agreement between the literature-based in vivo strain measurements and the simulated strain results, it can be concluded that the flexible multibody approach enables reasonable predictions of bone strain in response to dynamic loading. The information obtained from the present approach can be useful in clinical applications including devising exercises to prevent bone fragility or to accelerate fracture healing.

Molecular dynamics study of response of liquid N,N-dimethylformamide to externally applied electric field using a polarizable force field

The behavior of Liquid N,N-dimethylformamide subjected to a wide range of externally applied electric fields (from 0.001 V/nm to 1 V/nm) has been investigated through molecular dynamics simulation. To approach the objective the AMOEBA polarizable force field was extended to include the interaction of the external electric field with atomic partial charges and the contribution to the atomic polarization. The simulation results were evaluated with quantum mechanical calculations. The results from the present force field for the liquid at normal conditions were compared with the experimental and molecular dynamics results with non-polarizable and other polarizable force fields. The uniform external electric fields of higher than 0.01 V/nm have a significant effect on the structure of the liquid, which exhibits a variation in numerous properties, including molecular polarization, local cluster structure, rotation, alignment, energetics, and bulk thermodynamic and structural properties.

Molecular dynamics modelling of diffusional formation of titanium carbide clusters in iron matrix

Molecular dynamics simulation was employed to study the atomic interactions in titanium carbides and iron matrix containing carbon and titanium, which are significant for understanding the formation of titanium carbide cluster during precipitate process. The atoms trajectory and diffusion coefficients of carbon in titanium carbide were analyzed to provide a vacancy-exchanging mechanism and clarify the carbon concentration dependence of carbon diffusion in titanium carbide. The dependence of the formation of titanium carbide cluster in iron matrix on carbon was determined from the study of atoms diffusivity, cluster formation and formation energy of titanium carbide cluster. The simulation results provided insight into the carbon diffusion process and improved the understanding of the formation of titanium carbide cluster.

Numerical simulation of a Savonius turbine above an infinite-width forward facing step

The Savonius turbine, although simple in construction, typically has a maximum power coefficient (cP) of about 0.2. This is significantly lower than the cP of the axial flow propeller-type turbine which typically can be as high as 0.5. However, a simple means to improve the cP of a Savonius turbine is to install it above a forward facing step, for example, a cliff or a building. In this work, prior experimental results of the tow testing of a Savonius turbine installed above a finite-width bluff body were used to validate computational fluid dynamics simulation of the same experimental conditions. The validated simulation settings were then used to obtain the maximum cP of a similar turbine of finite width but installed above an infinite-width forward facing step over a range of installation positions above and behind the step.

Effects of quenching rate on amorphous structures of Cu46Zr54 metallic glass

Metallic glass shows some superior properties different from crystalline, but the nature of amorphous structure and structural change during glass transition have not been completely understood yet. Molecular dynamics simulation provides intuitive insight into the microstructure and properties at atomistic level. Before probing into the microstructures of metallic glass with molecular dynamics (MD) simulation, it is important to obtain amorphous state first. In the current work, we reproduce the process of manufacturing metallic glass in laboratory including the melting, equilibrating and quenching procedure with molecular dynamics simulations. The structure changing at melting point and glass transition temperature are investigated with the different cooling processing. The partial radial distribution function (PRDF) is applied as a criterion to judge the final amorphous state obtained considering the quenching at different cooling rates and the effects of cooling rate on the formation of amorphous structures are further discussed.

Effects of quenching rate on amorphous structures of Cu46Zr54 metallic glass

Metallic glass shows some superior properties different from crystalline, but the nature of amorphous structure and structural change during glass transition have not been completely understood yet. Molecular dynamics simulation provides intuitive insight into the microstructure and properties at atomistic level. Before probing into the microstructures of metallic glass with molecular dynamics (MD) simulation, it is important to obtain amorphous state first. In the current work, we reproduce the process of manufacturing metallic glass in laboratory including the melting, equilibrating and quenching procedure with molecular dynamics simulations. The structure changing at melting point and glass transition temperature are investigated with the different cooling processing. The partial radial distribution function (PRDF) is applied as a criterion to judge the final amorphous state obtained considering the quenching at different cooling rates and the effects of cooling rate on the formation of amorphous structures are further discussed.

C-BN single-walled nanotubes from hybrid connection of BN/C nanoribbons : prediction by ab initio density functional calculations

We demonstrated for the first time by ab initio density functional calculation and molecular dynamics simulation that C0.5(BN)0.5 armchair single-walled nanotubes (NT) are gapless semiconductors and can be spontaneously formed via the hybrid connection of graphene/BN Nanoribbons (GNR/BNNR) at room temperature. The direct synthesis of armchair C0.5(BN)0.5 via the hybrid connection of GNR/BNNR is predicted to be both thermodynamically and dynamically stable. Such novel armchair C0.5(BN)0.5 NTs possess enhanced conductance as that observed in GNRs. Additionally, the zigzag C0.5(BN)0.5 SWNTs are narrow band gap semiconductors, which may have potential application for light emission. In light of recent experimental progress and the enhanced degree of control in the synthesis of GNRs and BNNR, our results highlight an interesting avenue for synthesizing a novel specific type of C0.5(BN)0.5 nanotube (gapless or narrow direct gap semiconductor), with potentially important applications in BNC-based nanodevices.

A VMD plugin for NAMD simulations on Amazon EC2

VMD and NAMD are two major molecular dynamics simulation software packages, which can work together for mining structural information of bio-molecules. Carrying out such molecular dynamics simulations can help researchers to understand the roles and functions of various bio-molecules in life science research. Recently, clouds have provided HPC clusters on demand that allow users to benefit from their flexibility, elasticity, and lower costs. Although cloud computing promises to provide seamless access to HPC clusters through the abstraction of services, which hide the details of the underlying software and hardware infrastructure, users without in depth computing knowledge are still forced to cope with many low level system and programming details. Therefore, we have designed and developed a software plugin of VMD, which can provide an integrated framework for NAMD to be executed on Amazon EC2. The proposed Amazon EC2 Plugin for VMD frees users from performing many tedious computing tasks such as launching, connecting and terminating Amazon EC2 compute instances; configuring a HPC cluster; and installing middleware and software applications before the system is readily available for any scientific investigation. This allows VMD/NAMD users to spend less time getting applications to work on HPC clusters but more time for bio-research.