Molecular dynamics simulation model for the quantitative assessment of tool wear during single point diamond turning of cubic silicon carbide

Silicon carbide (SiC) is a material of great technological interest for engineering applications concerning hostile environments where silicon-based components cannot work (beyond 623 K). Single point diamond turning (SPDT) has remained a superior and viable method to harness process efficiency and freeform shapes on this harder material. However, it is extremely difficult to machine this ceramic consistently in the ductile regime due to sudden and rapid tool wear. It thus becomes non trivial to develop an accurate understanding of tool wear mechanism during SPDT of SiC in order to identify measures to suppress wear to minimize operational cost.

In this paper, molecular dynamics (MD) simulation has been deployed with a realistic analytical bond order potential (ABOP) formalism based potential energy function to understand tool wear mechanism during single point diamond turning of SiC. The most significant result was obtained using the radial distribution function which suggests graphitization of diamond tool during the machining process. This phenomenon occurs due to the abrasive processes between these two ultra hard materials. The abrasive action results in locally high temperature which compounds with the massive cutting forces leading to sp3–sp2 order–disorder transition of diamond tool. This represents the root cause of tool wear during SPDT operation of cubic SiC. Further testing led to the development of a novel method for quantitative assessment of the progression of diamond tool wear from MD simulations.

Association and liquid structure of pyridine-acetic acid mixtures determined from neutron scattering using a ‘free proton’ EPSR simulation model

The liquid structure of pyridine-acetic acid mixtures have been investigated using neutron scattering at various mole fractions of acetic acid, χHOAc = 0.33, 0.50, and 0.67, and compared to the structures of neat pyridine and acetic acid. Data has been modelled using Empirical Potential Structure Refinement (EPSR) with a ‘free proton’ reference model, which has no prejudicial weighting towards either the existence of molecular or ionised species. Analysis of the neutron scattering results shows the existence of hydrogen-bonded acetic acid chains with pyridine inclusions, rather than the formation of an ionic liquid by proton transfer.

Stability and dynamical behavior in a lake-model and implications for regime shifts in real lakes

A modified abstract version of the Comprehensive Aquatic Simulation Model (CASM) is found to exhibit three types of folded bifurcations due to nutrient loading. The resulting bifurcation diagrams account for nonlinear dynamics such as regime shifts and cyclic changes between clear-water state and turbid state that have actually been observed in real lakes. In particular, pulse-perturbation simulations based on the model presented suggest that temporal behaviors of real lakes after biomanipulations can be explained by pulse-dynamics in complex ecosystems, and that not only the amplitude (manipulated abundance of organisms) but also the phase (timing) is important for restoring lakes by biomanipulation. Ecosystem management in terms of possible irreversible changes in ecosystems induced by regime shifts is also discussed. (c) 2007 Elsevier B.V All rights reserved.

Stress correlation between instrumentation and simulation analysis of the die for high pressure die casting

A strain gauge instrumentation trial on a high pressure die casting ‘HPDC’ die was compared to a corresponding simulation model using Magmasoft® casting simulation software at two strain gauge rosette locations. The strains were measured during the casting cycle, from which the von Mises stress was determined and then compared to the simulation model. The von Mises stress from the simulation model correlated well with the findings from the instrumentation trial, showing a difference of 5.5%, ~ 10 MPa for one strain gauge rosette located in an area of low stress gradient. The second rosette was in a region of steep stress gradient, which resulted in a difference of up to 40%, ~40 MPa between the simulation and instrumentation results. Factors such as additional loading from die closure force or metal injection pressure which are not modelled by Magmasoft® were seen to have very little influence on the stress in the die, less than 7%.

Model reduction for nonlinear aerodynamics and aeroelasticity using a Discrete Empirical Interpolation Method

A novel surrogate model is proposed in lieu of Computational Fluid Dynamics (CFD) solvers, for fast nonlinear aerodynamic and aeroelastic modeling. A nonlinear function is identified on selected interpolation points by
a discrete empirical interpolation method (DEIM). The flow field is then reconstructed using a least square approximation of the flow modes extracted
by proper orthogonal decomposition (POD). The aeroelastic reduce order
model (ROM) is completed by introducing a nonlinear mapping function
between displacements and the DEIM points. The proposed model is investigated to predict the aerodynamic forces due to forced motions using
a N ACA 0012 airfoil undergoing a prescribed pitching oscillation. To investigate aeroelastic problems at transonic conditions, a pitch/plunge airfoil
and a cropped delta wing aeroelastic models are built using linear structural models. The presence of shock-waves triggers the appearance of limit
cycle oscillations (LCO), which the model is able to predict. For all cases
tested, the new ROM shows the ability to replicate the nonlinear aerodynamic forces, structural displacements and reconstruct the complete flow
field with sufficient accuracy at a fraction of the cost of full order CFD
model.

Economic modelling of antenatal screening and ultrasound scanning programmes for identification of fetal abnormalities

Objective Within the framework of a health technology assessment and using an economic model, to determine the most clinically and cost effective policy of scanning and screening for fetal abnormalities in early pregnancy. Design A discrete event simulation model of 50,000 singleton pregnancies. Setting Maternity services in Scotland. Population Women during the first 24 weeks of their pregnancy. Methods The mathematical model was populated with data on uptake of screening, prevalence, detection and false positive rates for eight fetal abnormalities and with costs for ultrasound scanning and serum screening. Inclusion of abnormalities was based on the relative prevalence and clinical importance of conditions and the availability of data. Six strategies for the identification of abnormalities prenatally including combinations of first and second trimester ultrasound scanning and first and second trimester screening for chromosomal abnormalities were compared. Main outcome measures The number of abnormalities detected and missed, the number of iatrogenic losses resulting from invasive tests, the total cost of strategies and the cost per abnormality detected were compared between strategies. Results First trimester screening for chromosomal abnormalities costs more than second trimester screening but results in fewer iatrogenic losses. Strategies which include a second trimester ultrasound scan result in more abnormalities being detected and have lower costs per anomaly detected. Conclusions The preferred strategy includes both first and second trimester ultrasound scans and a first trimester screening test for chromosomal abnormalities. It has been recommended that this policy is offered to all women in Scotland.

The feeding dynamics of broiler chickens

Contrary to a commonly held belief that broiler chickens need more space, there is increasing evidence that these birds are attracted to other birds. Indeed, commercially farmed birds exhibit a range of socially facilitated behaviours, such as increased feeding and preening in response to the presence of other birds. Social facilitation can generate feedback loops, whereby the adoption of a particular behaviour can spread rapidly and suddenly through the population. Here, by measuring the rate at which broiler chickens join and leave a feeding trough as a function of the number of birds already there, we quantify social facilitation. We use these measurements to parameterize a simulation model of chicken feeding behaviour. This model predicts, and further observations of broiler chickens confirm, that social facilitation leads to excitatory and synchronized patterns of group feeding. Such models could prove a powerful tool in understanding how feeding patterns depend on broiler house design.

Lake restoration in terms of ecological resilience: a numerical study of biomanipulations under bistable conditions

An abstract version of the comprehensive aquatic simulation model (CASM) is found to exhibit bistability under intermediate loading of nutrient input, supporting the alternative-stable-states theory and field observations for shallow lakes. Our simulations of biomanipulations under the bistable conditions reveal that a reduction in the abundance of zooplanktivorous fish cannot switch the system from a turbid to a clear state. Rather, a direct reduction of phytoplankton and detritus was found to be most effective to make this switch in the present model. These results imply that multiple manipulations may be effective for practical restorations of lakes. We discuss the present results of biomanipulations in terms of ecological resilience in multivariable systems or natural systems.

Determination of the heat transfer coefficient at the metal–die interface for high pressure die cast AlSi9Cu3Fe

When simulating the High Pressure Die Casting ‘HPDC’ process, the heat transfer coefficient ‘HTC’ between the casting and the die is critical to accurately predict the quality of the casting. To determine the HTC at the metal–die interface a production die for an automotive engine bearing beam, Die 1, was instrumented with type K thermocouples. A Magmasoft® simulation model was generated with virtual thermocouple points placed in the same location as the production die. The temperature traces from the simulation model were compared to the instrumentation results. Using the default simulation HTC for the metal–die interface, a poor correlation was seen, with the temperature response being much less for the simulation model. Because of this, the HTC at the metal–die interface was modified in order to get a better fit. After many simulation iterations, a good fit was established using a peak HTC of 42,000 W/m2 K, this modified HTC was further validated by a second instrumented production die, proving that the modified HTC gives good correlation to the instrumentation trials. The updated HTC properties for the simulation model will improve the predictive capabilities of the casting simulation software and better predict casting defects.

Performance models of storage contention in cloud environments

We propose simple models to predict the performance degradation of disk requests due to storage device contention in consolidated virtualized environments. Model parameters can be deduced from measurements obtained inside Virtual Machines (VMs) from a system where a single VM accesses a remote storage server. The parameterized model can then be used to predict the effect of storage contention when multiple VMs are consolidated on the same server. We first propose a trace-driven approach that evaluates a queueing network with fair share scheduling using simulation. The model parameters consider Virtual Machine Monitor level disk access optimizations and rely on a calibration technique. We further present a measurement-based approach that allows a distinct characterization of read/write performance attributes. In particular, we define simple linear prediction models for I/O request mean response times, throughputs and read/write mixes, as well as a simulation model for predicting response time distributions. We found our models to be effective in predicting such quantities across a range of synthetic and emulated application workloads. 

An investigation into the acceleration response of a damaged beam-type structure to a moving force

In recent years there have been a growing number of publications on procedures for damage detection in beams from analysing their dynamic response to the passage of a moving force. Most of this research demonstrates their effectiveness by showing that a singularity that did not appear in the healthy structure is present in the response of the damaged structure. This paper elucidates from first principles how the acceleration response can be assumed to consist of ‘static’ and ‘dynamic’ components, and where the beam has experienced a localised loss in stiffness, an additional ‘damage’ component. The combination of these components establishes how the damage singularity will appear in the total response. For a given damage severity, the amplitude of the ‘damage’ component will depend on how close the damage location is to the sensor, and its frequency content will increase with higher velocities of the moving force. The latter has implications for damage detection because if the frequency content of the ‘damage’ component includes bridge and/or vehicle frequencies, it becomes more difficult to identify damage. The paper illustrates how a thorough understanding of the relationship between the ‘static‘ and ‘damage’ components contributes to establish if damage has occurred and to provide an estimation of its location and severity. The findings are corroborated using accelerations from a planar finite element simulation model where the effects of force velocity and bridge span are examined.

Anisotropy of single crystal 3C-SiC during nanometric cutting

3C–SiC (the only polytype of SiC that resides in a diamond cubic lattice structure) is a relatively new material that exhibits most of the desirable engineering properties required for advanced electronic applications. The anisotropy exhibited by 3C–SiC during its nanometric cutting is significant, and the potential for its exploitation has yet to be fully investigated. This paper aims to understand the influence of crystal anisotropy of 3C–SiC on its cutting behaviour. A molecular dynamics simulation model was developed to simulate the nanometric cutting of single-crystal 3C–SiC in nine (9) distinct combinations of crystal orientations and cutting directions, i.e. (1?1?1) <-1?1?0>, (1?1?1) <-2?1?1>, (1?1?0) <-1?1?0>, (1?1?0) <0?0?1>, (1?1?0) <1?1?-2>, (0?0?1) <-1?1?0>, (0?0?1) <1?0?0>, (1?1?-2) <1?-1?0> and (1?-2?0) <2?1?0>.

In order to ensure the reliability of the simulation results, two separate simulation trials were carried out with different machining parameters. In the first trial, a cutting tool rake angle of -25°, d/r (uncut chip thickness/cutting edge radius) ratio of 0.57 and cutting velocity of 10 m s-1 were used whereas a second trial was done using a cutting tool rake angle of -30°, d/r ratio of 1 and cutting velocity of 4 m s-1. Both the trials showed similar anisotropic variation.

The simulated orthogonal components of thrust force in 3C–SiC showed a variation of up to 45%, while the resultant cutting forces showed a variation of 37%. This suggests that 3C–SiC is highly anisotropic in its ease of deformation. These results corroborate with the experimentally observed anisotropic variation of 43.6% in Young's modulus of 3C–SiC. The recently developed dislocation extraction algorithm (DXA) [1, 2] was employed to detect the nucleation of dislocations in the MD simulations of varying cutting orientations and cutting directions. Based on the overall analysis, it was found that 3C–SiC offers ease of deformation on either (1?1?1) <-1?1?0>, (1?1?0) <0?0?1>, or (1?0?0) <1?0?0> setups.