Modelling the nano-scale phenomena in condensed matter physics via computer-based numerical simulations

A review of the atomistic modelling of the behaviour of nano-scale structures and processes via molecular dynamics (MD) simulation method of a canonical ensemble is presented. Three areas of application in condensed matter physics are considered. We focus on the adhesive and indentation properties of the solid surfaces in nano-contacts, the nucleation and growth of nano-phase metallic and semi-conducting atomic and molecular films on supporting substrates, and the nano- and multi-scale crack propagation properties of metallic lattices. A set of simulations selected from these fields are discussed, together with a brief introduction to the methodology of the MD simulation. The pertinent inter-atomic potentials that model the energetics of the metallic and semi-conducting systems are also given.

Nanoscopic modelling of the adhesion, indentation and fracture characteristics of metallic systems via molecular dynamics simulations

Large-scale molecular dynamics simulations have been performed on canonical ensembles to model the adhesion and indentation characteristics of 3-D metallic nano-scale junctions in tip-substrate geometries, and the crack propagation in 2-D metallic lattices. It is shown that irreversible flows in nano-volumes of materials control the behaviour of the 3-D nano-contacts, and that local diffusional flow constitutes the atomistic mechanism underlying these plastic flows. These simulations show that the force of adhesion in metallic nano-contacts is reduced when adsorbate monolayers are present at the metal—metal junctions. Our results are in agreement with the conclusions of very accurate point-contact experiments carried out in this field. Our fracture simulations reveal that at low temperatures cleavage fractures can occur in both an elemental metal and an alloy. At elevated temperatures, the nucleation of dislocations is shown to cause a brittle-to-ductile transition. Limiting crack propagation velocities are computed for different strain rates and a dynamic instability is shown to control the crack movement beyond this limiting velocity, in line with the recent experimental results.

The effects of thermo-electrically induced convection in alloy solidification

The effects of a constant uniform magnetic field on thermoelectric currents during dendritic solidification were investigated using a 2-dimensional enthalpy based numerical model. Using an approximation of the dendrite growing in free space it was found that the resulting Lorentz force generates a circulating flow influencing the solidification pattern. As the magnetic field strength increases it was found that secondary growth on the clockwise side of the primary arm of the dendrite was encouraged, while the anticlockwise side is suppressed due to a reduction in local free energy. The preferred direction of growth rotated in the clockwise sense under an anti-clockwise flow for both the binary alloy and pure material. The tip velociy is significantly increased compared to growth in stagnant flow. This is due to a small recirculation that follows the tip of the dendrite; bringing in colder liquid and lower concentrations of solute. The recirculation being not normally incident on the tip is most likely the cause for the rotation. Grain growth consisting of multiple seeds with the same anisotropy growing in the same plane, gives a competition to release latent heat resulting in stunted growth. The initial growth for each dendrite is very similar to the single seed cases indicating that dendrites must become before the thermoelectric interactions are significant.

Modelling the tilt-casting process for the tranquil filling of titanium alloy turbine blades

The tilt-casting method is used to achieve tranquil filling of gamma-TiAl turbine blades. The reactive alloy is melted in a cold crucible using an induction coil and then the complete crucible-mould- running system assembly is rotated through 180degrees to transfer the metal into the mould. The induction current is ramped down gradually as the rotation starts and the mould is preheated to maintain superheat. The liquid metal then enters the mould and the gas within it (argon) escapes through the inlet aperture and through auxiliary vents. Solidification starts as soon the metal enters the mould and it is important to account for this effect to predict and prevent misruns. The rotation rate has to be controlled carefully to allow sufficient time for gas evacuation, but at the same time preserve superheat. This 3-phase system is modelled using the FV method, with a fast implicit numerical scheme used to capture the transient liquid free surface. The enthalpy method is used to model solidification and predict defects such as trapped bubbles, macro-porosity or surface connected porosity. Modeling is used to support an experimental program for the development of a production method for gamma-TiAl blades, with a target length of 40cm. The experiments provide validation for the model and the model in turn optimizes the tilt-casting process. The work is part of the EU project IMPRESS.

Effects of magnetic fields on crystal growth

The effects of a constant uniform magnetic field on thermoelectric currents during dendritic solidification were investigated using a two-dimensional enthalpy based numerical model. Using an approximation for three-dimensional unconstricted growth, the resulting Lorentz forces generate a circulating flow influencing the solidification pattern. Under the presence of a strong magnetic field secondary growth on the clockwise side of the primary arm of the dendrite was encouraged, whereas the anticlockwise side is suppressed due to a reduction in local free energy. The preferred direction of growth rotated in the clockwise sense under an anticlockwise flow. The tip velocity is significantly increased compared with growth in stagnant flow. This is due to a small recirculation at the tip of the dendrite; bringing in colder liquid and lowering the concentration of solute.

Investigation of wall-slip effect on lead-free solder paste and isotropic conductive adhesives

Wall-slip plays an important role in characterising the flow behaviour of solder paste materials. The wall slip arises due to the various attractive and repulsive forces acting between the solder particles and the walls of the measuring geometry.These interactions could lead to the presence of a thin solvent layer adjacent to the wall, which gives rise to slippage. The wall slip effect can play an important role in ensure successfulpaste release after the printing process. Wall-slip plays animportant role in characterising the flow behaviour of solderpastes and isotropic conductive adhesives. The study investigates the wall-slip formation in solder paste andisotropic conductive adhesives using flow visualisation technique. The slip distance was measured for parallel plate with different surface roughness in order to quantify the wallslip formations in these paste materials. An ink marker line was drawn between the parallel plate and the free surface of the sample. The parallel was rotated slowly at a constant shear rate of 0.05 sec-1 and the displacement of the ink marker was observed using a video microscope and image capturing software was utilised to capture the displacement of ink marker. From this study, it was found that the wall-slip effect was evident in all the paste materials. In addition, the different surface roughness of the parallel plates did not prevent the formation of wall-slip. This study has revealed that the wallslip effect could used to understand the flow behaviour of the paste in the stencil printing process.