Stabilized SPH-based simulations of impact dynamics using acceleration-corrected artificial viscosity

Artificial viscosity in SPH-based computations of impact dynamics is a numerical artifice that helps stabilize spurious oscillations near the shock fronts and requires certain user-defined parameters. Improper choice of these parameters may lead to spurious entropy generation within the discretized system and make it over-dissipative. This is of particular concern in impact mechanics problems wherein the transient structural response may depend sensitively on the transfer of momentum and kinetic energy due to impact. In order to address this difficulty, an acceleration correction algorithm was proposed in Shaw and Reid (''Heuristic acceleration correction algorithm for use in SPH computations in impact mechanics'', Comput. Methods Appl. Mech. Engrg., 198, 3962-3974) and further rationalized in Shaw et al. (An Optimally Corrected Form of Acceleration Correction Algorithm within SPH-based Simulations of Solid Mechanics, submitted to Comput. Methods Appl. Mech. Engrg). It was shown that the acceleration correction algorithm removes spurious high frequency oscillations in the computed response whilst retaining the stabilizing characteristics of the artificial viscosity in the presence of shocks and layers with sharp gradients. In this paper, we aim at gathering further insights into the acceleration correction algorithm by further exploring its application to problems related to impact dynamics. The numerical evidence in this work thus establishes that, together with the acceleration correction algorithm, SPH can be used as an accurate and efficient tool in dynamic, inelastic structural mechanics. (C) 2011 Elsevier Ltd. All rights reserved.

Impact dynamics of high Weber number drops on chemically modified metallic surfaces

Hydrophobic/superhydrophobic metallic surfaces prepared via chemical treatment are encountered in many industrial scenarios involving the impingement of spray droplets. The effectiveness of such surfaces is understood through the analysis of droplet impact experiments. In the present study, three target surfaces with aluminum (Al-6061) as base material-acid-etched, Octadecyl Trichloro Silane (OTS) coated, and acid-etched plus OTS-coated-were prepared. Experiments on the impact of inertia dominated water drops on these chemically modified aluminum surfaces were carried out with the objective to highlight the effect of chemical treatment on the target surfaces on key sub-processes occurring in drop impact phenomenon. High speed videos of the entire drop impact dynamics were captured at three Weber number (We) conditions representative of high We (We > 200) regime. During the early stages of drop spreading, the drop impact resulted in ejection of secondary droplets from spreading drop front on the etched surfaces resembling prompt splash on rough surfaces whereas no such splashing was observable on untreated aluminum surface. Prominent development of undulations (fingers) were observed at the rim of drop spreading on the etched surfaces; between the etched surfaces the OTS-coated surface showed a subdued development of fingers than the uncoated surface. The impacted drops showed intense receding on OTS-coated surfaces whereas on the etched surface a highly irregular receding, with drop liquid sticking to the surface, was observed. Quantitative analyses were performed to reveal the effect of target surface characteristics on drop impact parameters such as temporal variation of spread factor of drop lamella, temporal variation of average finger length during spreading phase, maximum drop spreading, time taken to attain maximum spreading, sensitivity of maximum spreading to We, number of fingers at maximum spreading, and average receding velocity of drop lamella. Existing models for maximum drop spreading showed reasonably good agreement with the experimental measurements on the target surfaces except the acid-etched surface. (C) 2014 Elsevier B.V. All rights reserved.

Feedback control of impact dynamics: The bouncing ball revisited

We study the the design of a tracking controller for the popular bouncing ball model: the continuous-time actuation of a table is used to control the impacts of the table with a bouncing ball. The proposed control law uses the impact times as the sole feedback information. We show that the acceleration of the table at impact plays no role in the stability analysis but is an important parameter for the robustness of the feedback system to model uncertainty, in particular to the uncertainty on the coefficient of restitution. © 2006 IEEE.

Physicochemical characteristics and droplet impact dynamics of superhydrophobic carbon nanotube arrays.

The physicochemical and droplet impact dynamics of superhydrophobic carbon nanotube arrays are investigated. These superhydrophobic arrays are fabricated simply by exposing the as-grown carbon nanotube arrays to a vacuum annealing treatment at a moderate temperature. This treatment, which allows a significant removal of oxygen adsorbates, leads to a dramatic change in wettability of the arrays, from mildly hydrophobic to superhydrophobic. Such change in wettability is also accompanied by a substantial change in surface charge and electrochemical properties. Here, the droplet impact dynamics are characterized in terms of critical Weber number, coefficient of restitution, spreading factor, and contact time. Based on these characteristics, it is found that superhydrophobic carbon nanotube arrays are among the best water-repellent surfaces ever reported. The results presented herein may pave a way for the utilization of superhydrophobic carbon nanotube arrays in numerous industrial and practical applications, including inkjet printing, direct injection engines, steam turbines, and microelectronic fabrication.

Molecular dynamics studies of nanoparticle impacts

This thesis concerns the dynamics of nanoparticle impacts on solid surfaces. These impacts occur, for instance, in space, where micro- and nanometeoroids hit surfaces of planets, moons, and spacecraft. On Earth, materials are bombarded with nanoparticles in cluster ion beam devices, in order to clean or smooth their surfaces, or to analyse their elemental composition. In both cases, the result depends on the combined effects of countless single impacts. However, the dynamics of single impacts must be understood before the overall effects of nanoparticle radiation can be modelled. In addition to applications, nanoparticle impacts are also important to basic research in the nanoscience field, because the impacts provide an excellent case to test the applicability of atomic-level interaction models to very dynamic conditions. In this thesis, the stopping of nanoparticles in matter is explored using classical molecular dynamics computer simulations. The materials investigated are gold, silicon, and silica. Impacts on silicon through a native oxide layer and formation of complex craters are also simulated. Nanoparticles up to a diameter of 20 nm (315000 atoms) were used as projectiles. The molecular dynamics method and interatomic potentials for silicon and gold are examined in this thesis. It is shown that the displacement cascade expansionmechanism and crater crown formation are very sensitive to the choice of atomic interaction model. However, the best of the current interatomic models can be utilized in nanoparticle impact simulation, if caution is exercised. The stopping of monatomic ions in matter is understood very well nowadays. However, interactions become very complex when several atoms impact on a surface simultaneously and within a short distance, as happens in a nanoparticle impact. A high energy density is deposited in a relatively small volume, which induces ejection of material and formation of a crater. Very high yields of excavated material are observed experimentally. In addition, the yields scale nonlinearly with the cluster size and impact energy at small cluster sizes, whereas in macroscopic hypervelocity impacts, the scaling 2 is linear. The aim of this thesis is to explore the atomistic mechanisms behind the nonlinear scaling at small cluster sizes. It is shown here that the nonlinear scaling of ejected material yield disappears at large impactor sizes because the stopping mechanism of nanoparticles gradually changes to the same mechanism as in macroscopic hypervelocity impacts. The high yields at small impactor size are due to the early escape of energetic atoms from the hot region. In addition, the sputtering yield is shown to depend very much on the spatial initial energy and momentum distributions that the nanoparticle induces in the material in the first phase of the impact. At the later phases, the ejection of material occurs by several mechanisms. The most important mechanism at high energies or at large cluster sizes is atomic cluster ejection from the transient liquid crown that surrounds the crater. The cluster impact dynamics detected in the simulations are in agreement with several recent experimental results. In addition, it is shown that relatively weak impacts can induce modifications on the surface of an amorphous target over a larger area than was previously expected. This is a probable explanation for the formation of the complex crater shapes observed on these surfaces with atomic force microscopy. Clusters that consist of hundreds of thousands of atoms induce long-range modifications in crystalline gold.

Impact of water drops onto the junction of a hydrophobic texture and a hydrophilic smooth surface

Wettability gradient surfaces play a significant role in control and manipulation of liquid drops. The present work deals with the analysis of water drops impacting onto the junction line between hydrophobic texture and hydrophilic smooth portions of a dual-textured substrate made using stainless steel material. The hydrophobic textured portion of the substrate comprised of unidirectional parallel groove-like and pillar-like structures of uniform dimensions. A high-speed video camera recorded the spreading and receding dynamics of impacting drops. The drop impact dynamics during the early inertia driven impact regime remains unaffected by the dual-texture feature of the substrate. A larger retraction speed of drop liquid observed on the hydrophobic portion of the substrate during the impact of low velocity drops makes the drop liquid on the higher wettability portion to advance further (secondary drop spreading). The net horizontal drop velocity towards the hydrophilic portion of the dual-textured substrate decreases with increasing drop impact velocity. The available experimental results suggest that the movement of bulk drop liquid away from the impact point during drop impact on the dual-textured substrate is larger for the impact of low inertia drops. (C) 2010 Elsevier B.V. All rights reserved.

Optimised form of acceleration correction algorithm within SPH-based simulations of impact mechanics

In the context of SPH-based simulations of impact dynamics, an optimised and automated form of the acceleration correction algorithm (Shaw and Reid, 2009a) is developed so as to remove spurious high frequency oscillations in computed responses whilst retaining the stabilizing characteristics of the artificial viscosity in the presence of shocks and layers with sharp gradients. A rational framework for an insightful characterisation of the erstwhile acceleration correction method is first set up. This is followed by the proposal of an optimised version of the method, wherein the strength of the correction term in the momentum balance and energy equations is optimised. For the first time, this leads to an automated procedure to arrive at the artificial viscosity term. In particular, this is achieved by taking a spatially varying response-dependent support size for the kernel function through which the correction term is computed. The optimum value of the support size is deduced by minimising the (spatially localised) total variation of the high oscillation in the acceleration term with respect to its (local) mean. The derivation of the method, its advantages over the heuristic method and issues related to its numerical implementation are discussed in detail. (C) 2011 Elsevier Ltd. All rights reserved.

Experimental investigation of biofuel drop impact on stainless steel surface

Blends of conventional fuels such as Jet-A1 (aviation kerosene) and diesel with bio-derived components, referred to as biofttels, are gradually replacing the conventional fuels in aircraft and automobile engines. There is a lack of understanding on the interaction of spray drops of such biofuels with solid surfaces. The present study is an experimental investigation on the impact of biofuel drops onto a smooth stainless steel surface. The biofuel is a mixture of 90% commercially available camelina-derived biofuel and 10% aromatics. Biofuel drops were generated using a syringe-hypodermic needle arrangement. On demand, the needle delivers an almost spherical drop with drop diameter in the range 2.05-2.15 mm. Static wetting experiments show that the biofuel drop completely wets the stainless steel surface and exhibits an equilibrium contact angle of 5.6. High speed video camera was used to capture the impact dynamics of biofuel drops with Weber number ranging from 20 to 570. The spreading dynamics and maximum spreading diameter of impacting biofuel drops on the target surface were analyzed. For the impact of high Weber number biofuel drops, the spreading law suggests beta similar to tau(0.5) where beta is the spread factor and tau, the nondimensionalized time. The experimentally observed trend of maximum spread factor, beta(max) of camelina biofuel drop on the target surface with We compares well with the theoretically predicted trend from Ukiwe-Kwok model. After reaching beta(max), the impacting biofuel drop undergoes a prolonged sluggish spreading due to the high wetting nature of the camelina biofuel-stainless steel system. As a result, the final spread factor is found to be a little more than beta(max). (C) 2014 Elsevier Inc. All rights reserved.

Comparison of crash dynamics of articulated and non-articulated trains

This paper presents a comparison of impact dynamic performance between articulated trains and non-articulated trains. This is carried out by investigation of the characteristics of the two trains types and analysis of their effects on impact dynamics. The analysis shows that the differences in bogie support positions on the carbody and coupling devices lead to differences in several structural and compositional characteristics. These characteristics result in different impact responses for the two types of train and are directly related to their impact stablity. Articulated trains have stiff connection and integral performance in collisions but with less capability for absorbing impact energy between carriages, whereas non-articulated trains show loose connection and scattered performance in collisions but with more options for energy absorber installation between carriages.

Energy conserving schemes for the simulation of musical instrument contact dynamics

Collisions are an innate part of the function of many musical instruments. Due to the nonlinear nature of contact forces, special care has to be taken in the construction of numerical schemes for simulation and sound synthesis. Finite difference schemes and other time-stepping algorithms used for musical instrument modelling purposes are normally arrived at by discretising a Newtonian description of the system. However because impact forces are non-analytic functions of the phase space variables, algorithm stability can rarely be established this way. This paper presents a systematic approach to deriving energy conserving schemes for frictionless impact modelling. The proposed numerical formulations follow from discretising Hamilton׳s equations of motion, generally leading to an implicit system of nonlinear equations that can be solved with Newton׳s method. The approach is first outlined for point mass collisions and then extended to distributed settings, such as vibrating strings and beams colliding with rigid obstacles. Stability and other relevant properties of the proposed approach are discussed and further demonstrated with simulation examples. The methodology is exemplified through a case study on tanpura string vibration, with the results confirming the main findings of previous studies on the role of the bridge in sound generation with this type of string instrument.

A detailed study of the mechanical and electrical aspects of microparticle impact phenomena

An experimental and theoretical study of the impact behaviour of charged microparticles in a high voltage vacuum gap has been carried out to investigate under controlled conditions the role of low velocity microparticles (ζ 500 ms-1) in initiating electrical breakdown in such gaps. This has involved developing a unique (UHV) low-velocity source of micron-sized charged particles to study the underlying mechanical and electrical aspects of micro-particle impact on a range of target materials e.g. Pb, Ti, C, stainless-steel and mica etc., having atomically clean or oxidised surfaces. Argon-ion etching and electron-beam heating has been used for in-situ surface treatment and ellipsometry for characterising the target surfaces. An associated sphere/plane theoretical model has been developed for detailed analysis of the many complex electrical (in-flight in-field emission, M.I.M. tunnelling and ohmic conduction) and mechanical (impact dynamics, deformation and heating) phenomena that are involved when a microparticle closely approaches and impacts on a plane target. In each instance the influence of parameters such as particle radius, particle/target impact velocity, surface field, surface condition and material has been determined.

Influence of anvil supports and overhangs in dynamic analysis of three point bend testing

Optimum design of dynamic fracture test rigs demands a thorough appreciation of beam vibration under impact. Analyses invariably presume rigid anvils, and neglect overhang effects. The beam response predicted analytically and numerically in this paper highlights the significant role of anvil rigidity and beam overhangs on the impact dynamics of three point bend (3PB) specimens.

Adiabatic shear localization: occurrence, theories, and applications

Stable dynamic walking over uneven terrain

We propose a constructive control design for stabilization of non-periodic trajectories of underactuated robots. An important example of such a system is an underactuated "dynamic walking" biped robot traversing rough or uneven terrain. The stabilization problem is inherently challenging due to the nonlinearity, open-loop instability, hybrid (impact) dynamics, and target motions which are not known in advance. The proposed technique is to compute a transverse linearization about the desired motion: a linear impulsive system which locally represents "transversal" dynamics about a target trajectory. This system is then exponentially stabilized using a modified receding-horizon control design, providing exponential orbital stability of the target trajectory of the original nonlinear system. The proposed method is experimentally verified using a compass-gait walker: a two-degree-of-freedom biped with hip actuation but pointed stilt-like feet. The technique is, however, very general and can be applied to a wide variety of hybrid nonlinear systems. © The Author(s) 2011.

A spatial modelling framework for assessing climate change impacts on freshwater ecosystems: Response of brown trout (Salmo trutta L.) biomass to warming water temperature

Mountain regions worldwide are particularly sensitive to on-going climate change. Specifically in the Alps in Switzerland, the temperature has increased twice as fast than in the rest of the Northern hemisphere. Water temperature closely follows the annual air temperature cycle, severely impacting streams and freshwater ecosystems. In the last 20 years, brown trout (Salmo trutta L) catch has declined by approximately 40-50% in many rivers in Switzerland. Increasing water temperature has been suggested as one of the most likely cause of this decline. Temperature has a direct effect on trout population dynamics through developmental and disease control but can also indirectly impact dynamics via food-web interactions such as resource availability. We developed a spatially explicit modelling framework that allows spatial and temporal projections of trout biomass using the Aare river catchment as a model system, in order to assess the spatial and seasonal patterns of trout biomass variation. Given that biomass has a seasonal variation depending on trout life history stage, we developed seasonal biomass variation models for three periods of the year (Autumn-Winter, Spring and Summer). Because stream water temperature is a critical parameter for brown trout development, we first calibrated a model to predict water temperature as a function of air temperature to be able to further apply climate change scenarios. We then built a model of trout biomass variation by linking water temperature to trout biomass measurements collected by electro-fishing in 21 stations from 2009 to 2011. The different modelling components of our framework had overall a good predictive ability and we could show a seasonal effect of water temperature affecting trout biomass variation. Our statistical framework uses a minimum set of input variables that make it easily transferable to other study areas or fish species but could be improved by including effects of the biotic environment and the evolution of demographical parameters over time. However, our framework still remains informative to spatially highlight where potential changes of water temperature could affect trout biomass. (C) 2015 Elsevier B.V. All rights reserved.-