Numerical characterization of the mechanical properties of metal nanowires

Nanowires (NWs) have attracted appealing and broad application owing to their remarkable mechanical, optical, electrical, thermal and other properties. To unlock the revolutionary characteristics of NWs, a considerable body of experimental and theoretical work has been conducted. However, due to the extremely small dimensions of NWs, the application and manipulation of the in situ experiments involve inherent complexities and huge challenges. For the same reason, the presence of defects appears as one of the most dominant factors in determining their properties. Hence, based on the experiments' deficiency and the necessity of investigating different defects' influence, the numerical simulation or modelling becomes increasingly important in the area of characterizing the properties of NWs. It has been noted that, despite the number of numerical studies of NWs, significant work still lies ahead in terms of problem formulation, interpretation of results, identification and delineation of deformation mechanisms, and constitutive characterization of behaviour. Therefore, the primary aim of this study was to characterize both perfect and defected metal NWs. Large-scale molecular dynamics (MD) simulations were utilized to assess the mechanical properties and deformation mechanisms of different NWs under diverse loading conditions including tension, compression, bending, vibration and torsion. The target samples include different FCC metal NWs (e.g., Cu, Ag, Au NWs), which were either in a perfect crystal structure or constructed with different defects (e.g. pre-existing surface/internal defects, grain/twin boundaries). It has been found from the tensile deformation that Young's modulus was insensitive to different styles of pre-existing defects, whereas the yield strength showed considerable reduction. The deformation mechanisms were found to be greatly influenced by the presence of defects, i.e., different defects acted in the role of dislocation sources, and many affluent deformation mechanisms had been triggered. Similar conclusions were also obtained from the compressive deformation, i.e., Young's modulus was insensitive to different defects, but the critical stress showed evident reduction. Results from the bending deformation revealed that the current modified beam models with the considerations of surface effect, or both surface effect and axial extension effect were still experiencing certain inaccuracy, especially for the NW with ultra small cross-sectional size. Additionally, the flexural rigidity of the NW was found to be insensitive to different pre-existing defects, while the yield strength showed an evident decrease. For the resonance study, the first-order natural frequency of the NW with pre-existing surface defects was almost the same as that from the perfect NW, whereas a lower first-order natural frequency and a significantly degraded quality factor was observed for NWs with grain boundaries. Most importantly, the <110> FCC NWs were found to exhibit a novel beat phenomenon driven by a single actuation, which was resulted from the asymmetry in the lattice spacing in the (110) plane of the NW cross-section, and expected to exert crucial impacts on the in situ nanomechanical measurements. In particular, <110> Ag NWs with rhombic, truncated rhombic, and triangular cross-sections were found to naturally possess two first-mode natural frequencies, which were envisioned with applications in NEMS that could operate in a non-planar regime. The torsion results revealed that the torsional rigidity of the NW was insensitive to the presence of pre-existing defects and twin boundaries, but received evident reduction due to grain boundaries. Meanwhile, the critical angle decreased considerably for defected NWs. This study has provided a comprehensive and deep investigation on the mechanical properties and deformation mechanisms of perfect and defected NWs, which will greatly extend and enhance the existing knowledge and understanding of the properties/performance of NWs, and eventually benefit the realization of their full potential applications. All delineated MD models and theoretical analysis techniques that were established for the target NWs in this research are also applicable to future studies on other kinds of NWs. It has been suggested that MD simulation is an effective and excellent tool, not only for the characterization of the properties of NWs, but also for the prediction of novel or unexpected properties.

Numerical simulations of particle deposition in metal foam heat exchangers

Australia is a high-potential country for geothermal power with reserves currently estimated in the tens of millions of petajoules, enough to power the nation for at least 1000 years at current usage. However, these resources are mainly located in isolated arid regions where water is scarce. Therefore, wet cooling systems for geothermal plants in Australia are the least attractive solution and thus air-cooled heat exchangers are preferred. In order to increase the efficiency of such heat exchangers, metal foams have been used. One issue raised by this solution is the fouling caused by dust deposition. In this case, the heat transfer characteristics of the metal foam heat exchanger can dramatically deteriorate. Exploring the particle deposition property in the metal foam exchanger becomes crucial. This paper is a numerical investigation aimed to address this issue. Two dimensional (2D) numerical simulations of a standard one-row tube bundle wrapped with metal foam in cross-flow are performed and highlight preferential particle deposition areas.

Response of a buried tunnel to surface blast using different numerical techniques

This paper presents a comparative study on the response of a buried tunnel to surface blast using the arbitrary Lagrangian-Eulerian (ALE) and smooth particle hydrodynamics (SPH) techniques. Since explosive tests with real physical models are extremely risky and expensive, the results of a centrifuge test were used to validate the numerical techniques. The numerical study shows that the ALE predictions were faster and closer to the experimental results than those from the SPH simulations which over predicted the strains. The findings of this research demonstrate the superiority of the ALE modelling techniques for the present study. They also provide a comprehensive understanding of the preferred ALE modelling techniques which can be used to investigate the surface blast response of underground tunnels.

Numerical simulations of particle deposition in metal foam heat exchangers

Australia is a high potential country for geothermal power with reserves currently estimated in the tens of millions of petajoules, enough to power the nation for at least 1000 years at current usage.However, these resources are mainly located in isolated arid regions where water is scarce. Therefore, wet cooling systems for geothermal plants in Australia are the least attractive solution and thus air-cooled heat exchangers are preferred. In order to increase the efficiency of such heat exchangers, metal foams have been used. One issue raised by this solution is the fouling caused by dust deposition. In this case, the heat transfer characteristics of the metal foam heat exchanger can dramatically deteriorate. Exploring the particle deposition property in the metal foam exchanger becomes crucial. This paper is a numerical investigation aimed to address this issue. Two-dimensional(2D numerical simulations of a standard one-row tube bundle wrapped with metal foam in cross-flow are performed and highlight preferential particle deposition areas.

Numerical modelling of load bearing light gauge steel frame wall systems exposed to realistic design fires

Fire safety has become an important part in structural design due to the ever increasing loss of properties and lives during fires. Conventionally the fire rating of load bearing wall systems made of Light gauge Steel Frames (LSF) is determined using fire tests based on the standard time-temperature curve in ISO834 [1]. However, modern commercial and residential buildings make use of thermoplastic materials, which mean considerably high fuel loads. Hence a detailed fire research study into the fire performance of LSF walls was undertaken using realistic design fire curves developed based on Eurocode parametric [2] and Barnett’s BFD [3] curves using both full scale fire tests and numerical studies. It included LSF walls without cavity insulation, and the recently developed externally insulated composite panel system. This paper presents the details of finite element models developed to simulate the full scale fire tests of LSF wall panels under realistic design fires. Finite element models of LSF walls exposed to realistic design fires were developed, and analysed under both transient and steady state fire conditions using the measured stud time-temperature curves. Transient state analyses were performed to simulate fire test conditions while steady state analyses were performed to obtain the load ratio versus time and failure temperature curves of LSF walls. Details of the developed finite element models and the results including the axial deformation and lateral deflection versus time curves, and the stud failure modes and times are presented in this paper. Comparison with fire test results demonstrate the ability of developed finite element models to predict the performance and fire resistance ratings of LSF walls under realistic design fires.

Numerical study of lid-driven square cavity with heat generation using LBM

In this study, the mixed convection heat transfer and fluid flow behaviors in a lid-driven square cavity filled with high Prandtl number fluid (Pr = 5400, ν = 1.2×10-4 m2/s) at low Reynolds number is studied using thermal Lattice Boltzmann method (TLBM) where ν is the viscosity of the fluid. The LBM has built up on the D2Q9 model and the single relaxation time method called the Lattice-BGK (Bhatnagar-Gross-Krook) model. The effects of the variations of non dimensional mixed convection parameter called Richardson number(Ri) with and without heat generating source on the thermal and flow behavior of the fluid inside the cavity are investigated. The results are presented as velocity and temperature profiles as well as stream function and temperature contours for Ri ranging from 0.1 to 5.0 with other controlling parameters that present in this study. It is found that LBM has good potential to simulate mixed convection heat transfer and fluid flow problem. Finally the simulation results have been compared with the previous numerical and experimental results and it is found to be in good agreement.

Numerical modelling of porous injection in a radical farming scramjet

A numerical investigation of the behaviour of fuel injection through a porous surface in an inlet-fuelled, radial-farming scramjet is presented. The performance of porous fuel injection is compared to discrete port hole injection at an equivalence ratio of φ ≈ 0.4 for both cases. The comparison is performed at a Mach 6.5 flow condition with a total specific enthalpy of 4.3 MJ/kg. The numerical results are compared to experiments performed in the T4 shock tunnel where available. The presented results demonstrate for the first time, that porous fuel injection has the potential to outperform port hole injectors in scramjet engines in terms of fuel-air mixing, ignition delays and achievable combustion efficiencies despite reduced fuel penetration heights.

Numerical investigation of the influence of topography on simulated downburst wind fields

A, dry, non-hydrostatic sub-cloud model is used to simulate an isolated stationary downburst wind event to study the influence topographic features have on the near-ground wind structure of these storms. It was generally found that storm maximum wind speeds could be increased by up to 30% because of the presence of a topographic feature at the location of maximum wind speeds. Comparing predicted velocity profile amplification with that of a steady flow impinging jet, similar results were found despite the simplifications made in the impinging jet model. Comparison of these amplification profiles with those found in the simulated boundary layer winds reveal reductions of up to 30% in the downburst cases. Downburst and boundary layer amplification profiles were shown to become more similar as the topographic feature height was reduced with respect to the outflow depth.

Numerical simulation of downburst winds

The wind field of an intense idealised downburst wind storm has been studied using an axisymmetric, dry, non-hydrostatic numerical sub-cloud model. The downburst driving processes of evaporation and melting have been paramaterized by an imposed cooling source that triggers and sustains a downdraft. The simulated downburst exhibits many characteristics of observed full-scale downburst events, in particular the presence of a primary and counter rotating secondary ring vortex at the leading edge of the diverging front. The counter-rotating vortex is shown to significantly influence the development and structure of the outflow. Numerical forcing and environmental characteristics have been systematically varied to determine the influence on the outflow wind field. Normalised wind structure at the time of peak outflow intensity was generally shown to remain constant for all simulations. Enveloped velocity profiles considering the velocity structure throughout the entire storm event show much more scatter. Assessing the available kinetic energy within each simulated storm event, it is shown that the simulated downburst wind events had significantly less energy available for loading isolated structures when compared with atmospheric boundary layer winds. The discrepancy is shown to be particularly prevalent when wind speeds were integrated over heights representative of tall buildings. A similar analysis for available full scale measurements led to similar findings.

Numerical simulation of idealised three-dimensional downburst wind fields

Convective downburst wind storms generate the peak annual gust wind speed for many parts of the non-cyclonic world at return periods of importance for ultimate limit state design. Despite this there is little clear understanding of how to appropriately design for these wind events given their significant dissimilarities to boundary layer winds upon which most design is based. To enhance the understanding of wind fields associated with these storms a three-dimensional numerical model was developed to simulate a multitude of idealised downburst scenarios and to investigate their near-ground wind characteristics. Stationary and translating downdraft wind events in still and sheared environments were simulated with baseline results showing good agreement with previous numerical work and full-scale observational data. Significant differences are shown in the normalised peak wind speed velocity profiles depending on the environmental wind conditions in the vicinity of the simulated event. When integrated over the height of mid- to high rise structures, all simulated profiles are shown to produce wind loads smaller than an equivalent 10 m height matched open terrain boundary layer profile. This suggests that for these structures the current design approach is conservative from an ultimate loading standpoint. Investigating the influence of topography on the structure of the simulated near-ground downburst wind fields, it is shown that these features amplify wind speeds in a manner similar to that expected for boundary layer winds, but the extent of amplification is reduced. The level of reduction is shown to be dependent on the depth of the simulated downburst outflow.

Numerical modelling and design of lipped channel beams subject to shear

Cold-formed steel members are increasingly used as primary structural elements in buildings due to the availability of thin and high strength steels and advanced cold-forming technologies. Cold-formed lipped channel beams (LCB) are commonly used as flexural members such as floor joists and bearers. Many research studies have been carried out to evaluate the behaviour and design of LCBs subject to pure bending actions. However, limited research has been undertaken on the shear behaviour and strength of LCBs. Hence a numerical study was undertaken to investigate the shear behaviour and strength of LCBs. Finite element models of simply supported LCBs with aspect ratios of 1.0 and 1.5 were considered under a mid-span load. They were then validated by comparing their results with test results and used in a detailed parametric study based on the validated finite element models. Numerical studies were conducted to investigate the shear buckling and post-buckling behaviour of LCBs. Experimental and numerical results showed that the current design rules in cold-formed steel structures design codes are very conservative for the shear design of LCBs. Improved design equations were therefore proposed for the shear strength of LCBs. This paper presents the details of this numerical study of LCBs and the results.

Numerical studies and design of hollow flange channel beams subject to combined bending and shear actions

LiteSteel beam (LSB) is a cold-formed steel hollow flange channel section produced using a patented manufacturing process involving simultaneous cold-forming and dual electric resistance welding. It is commonly used as floor joists and bearers in residential, industrial and commercial buildings. Design of the LSB is governed by the Australian cold-formed steel structures code, AS/NZS 4600. Due to the geometry of the LSB, as well as its unique residual stress characteristics and initial geometric imperfections resultant of manufacturing processes, currently available design equations for common cold-formed sections are not directly applicable to the LSB. Many research studies have been carried out to evaluate the behaviour and design of LSBs subject to pure bending actions and predominant shear actions. To date, however, no investigation has been conducted into the strength of LSB sections under combined bending and shear actions. Hence experimental and numerical studies were conducted to assess the combined bending and shear behaviour of LSBs. Finite element models of LSBs were developed to simulate their combined bending and shear behaviour and strength of LSBs. They were then validated by comparing the results with available experimental test results and used in a detailed parametric study. The results from experimental and finite element analyses were compared with current AS/NZS 4600 and AS 4100 design rules. Both experimental and numerical studies show that the AS/NZS 4600 design rule based on circular interaction equation is conservative in predicting the combined bending and shear capacities of LSBs. This paper presents the details of the numerical studies of LSBs and the results. In response to the inadequacies of current approaches to designing LSBs for combined bending and shear, two lower bound design equations are proposed in this paper.

Numerical estimation of 3D mechanical forces exerted by cells on non-linear materials

The exchange of physical forces in both cell-cell and cell-matrix interactions play a significant role in a variety of physiological and pathological processes, such as cell migration, cancer metastasis, inflammation and wound healing. Therefore, great interest exists in accurately quantifying the forces that cells exert on their substrate during migration. Traction Force Microscopy (TFM) is the most widely used method for measuring cell traction forces. Several mathematical techniques have been developed to estimate forces from TFM experiments. However, certain simplifications are commonly assumed, such as linear elasticity of the materials and/or free geometries, which in some cases may lead to inaccurate results. Here, cellular forces are numerically estimated by solving a minimization problem that combines multiple non-linear FEM solutions. Our simulations, free from constraints on the geometrical and the mechanical conditions, show that forces are predicted with higher accuracy than when using the standard approaches.

Numerical investigation of mechanical and thermal dynamic properties of the industrial transformer

Industrial transformer is one of the most critical assets in the power and heavy industry. Failures of transformers can cause enormous losses. The poor joints of the electrical circuit on transformers can cause overheating and results in stress concentration on the structure which is the major cause of catastrophic failure. Few researches have been focused on the mechanical properties of industrial transformers under overheating thermal conditions. In this paper, both mechanical and thermal properties of industrial transformers are jointly investigated using Finite Element Analysis (FEA). Dynamic response analysis is conducted on a modified transformer FEA model, and the computational results are compared with experimental results from literature to validate this simulation model. Based on the FEA model, thermal stress is calculated under different temperature conditions. These analysis results can provide insights to the understanding of the failure of transformers due to overheating, therefore are significant to assess winding fault, especially to the manufacturing and maintenance of large transformers.

Experimental and numerical analysis of the relationship between indoor and outdoor airborne particles in an operating room

This work investigated the impact of the HVAC filtration system and indoor particle sources on the relationship between indoor and outdoor airborne particle size and concentrations in an operating room. Filters with efficiency between 65% and 99.97% were used in the investigation and indoor and outdoor particle size and concentrations were measured. A balance mass model was used for the simulation of the impact of the surgical team, deposition rate, HVAC exhaust and air change rates on indoor particle concentration. The experimental results showed that high efficiency filters would not be expected to decrease the risk associated with indoor particles larger than approximately 1 µm in size because normal filters are relatively efficient for these large particles. A good fraction of outdoor particles were removed by deposition on the HVAC system surfaces and this deposition increased with particle size. For particles of 0.3-0.5 µm in diameter, particle reduction was about 23%, while for particles >10 µm the loss was about 78%. The modelling results showed that depending on the type of filter used, the surgical team generated between 93-99% of total particles, while the outdoor air contributed only 1-6%.