Computer experiment – a case study for modelling and simulation of manufacturing systems

Deterministic computer simulation of physical experiments is now a common technique in science and engineering. Often, physical experiments are too time consuming, expensive or impossible to conduct. Complex computer models or codes, rather than physical experiments lead to the study of computer experiments, which are used to investigate many scientific phenomena. A computer experiment consists of a number of runs of the computer code with different input choices. The Design and Analysis of Computer Experiments is a rapidly growing technique in statistical experimental design. This paper aims to discuss some practical issues when designing a computer simulation and/or experiments for manufacturing systems. A case study approach is reviewed and presented.

Some practical issues in the design and analysis of computer experiments

Deterministic computer simulations of physical experiments are now common techniques in science and engineering. Often, physical experiments are too time consuming, expensive or impossible to conduct. Complex computer models or codes, rather than physical experiments lead to the study of computer experiments, which are used to investigate many scientific phenomena of this nature. A computer experiment consists of a number of runs of the computer code with different input choices. The Design and Analysis of Computer Experiments is a rapidly growing technique in statistical experimental design. This thesis investigates some practical issues in the design and analysis of computer experiments and attempts to answer some of the questions faced by experimenters using computer experiments. In particular, the question of the number of computer experiments and how they should be augmented is studied and attention is given to when the response is a function over time.

Co-constructing decimal number knowledge

This mathematics education research provides significant insights for the teaching of decimals to children. It is well known that decimals is one of the most difficult topics to learn and teach. Annette’s research is unique in that it focuses not only on the cognitive, but also on the affective and conative aspects of learning and teaching of decimals. The study is innovative as it includes the students as co-constructors and co-researchers. The findings open new ways of thinking for educators about how students cognitively process decimal knowledge, as well as how students might develop a sense of self as a learner, teacher and researcher in mathematics.

The use of electronic portal imaging to verify patient position during intensity-modulated radiotherapy delivered by the dynamic MLC technique

Purpose: The precise shape of the three-dimensional dose distributions created by intensity-modulated radiotherapy means that the verification of patient position and setup is crucial to the outcome of the treatment. In this paper, we investigate and compare the use of two different image calibration procedures that allow extraction of patient anatomy from measured electronic portal images of intensity-modulated treatment beams. Methods and Materials: Electronic portal images of the intensity-modulated treatment beam delivered using the dynamic multileaf collimator technique were acquired. The images were formed by measuring a series of frames or segments throughout the delivery of the beams. The frames were then summed to produce an integrated portal image of the delivered beam. Two different methods for calibrating the integrated image were investigated with the aim of removing the intensity modulations of the beam. The first involved a simple point-by-point division of the integrated image by a single calibration image of the intensity-modulated beam delivered to a homogeneous polymethyl methacrylate (PMMA) phantom. The second calibration method is known as the quadratic calibration method and required a series of calibration images of the intensity-modulated beam delivered to different thicknesses of homogeneous PMMA blocks. Measurements were made using two different detector systems: a Varian amorphous silicon flat-panel imager and a Theraview camera-based system. The methods were tested first using a contrast phantom before images were acquired of intensity-modulated radiotherapy treatment delivered to the prostate and pelvic nodes of cancer patients at the Royal Marsden Hospital. Results: The results indicate that the calibration methods can be used to remove the intensity modulations of the beam, making it possible to see the outlines of bony anatomy that could be used for patient position verification. This was shown for both posterior and lateral delivered fields. Conclusions: Very little difference between the two calibration methods was observed, so the simpler division method, requiring only the single extra calibration measurement and much simpler computation, was the favored method. This new method could provide a complementary tool to existing position verification methods, and it has the advantage that it is completely passive, requiring no further dose to the patient and using only the treatment fields.

Improving ultrasound excitation systems using a flexible power supply with adjustable voltage and frequency to drive piezoelectric transducers

The ability of a piezoelectric transducer in energy conversion is rapidly expanding in several applications. Some of the industrial applications for which a high power ultrasound transducer can be used are surface cleaning, water treatment, plastic welding and food sterilization. Also, a high power ultrasound transducer plays a great role in biomedical applications such as diagnostic and therapeutic applications. An ultrasound transducer is usually applied to convert electrical energy to mechanical energy and vice versa. In some high power ultrasound system, ultrasound transducers are applied as a transmitter, as a receiver or both. As a transmitter, it converts electrical energy to mechanical energy while a receiver converts mechanical energy to electrical energy as a sensor for control system. Once a piezoelectric transducer is excited by electrical signal, piezoelectric material starts to vibrate and generates ultrasound waves. A portion of the ultrasound waves which passes through the medium will be sensed by the receiver and converted to electrical energy. To drive an ultrasound transducer, an excitation signal should be properly designed otherwise undesired signal (low quality) can deteriorate the performance of the transducer (energy conversion) and increase power consumption in the system. For instance, some portion of generated power may be delivered in unwanted frequency which is not acceptable for some applications especially for biomedical applications. To achieve better performance of the transducer, along with the quality of the excitation signal, the characteristics of the high power ultrasound transducer should be taken into consideration as well. In this regard, several simulation and experimental tests are carried out in this research to model high power ultrasound transducers and systems. During these experiments, high power ultrasound transducers are excited by several excitation signals with different amplitudes and frequencies, using a network analyser, a signal generator, a high power amplifier and a multilevel converter. Also, to analyse the behaviour of the ultrasound system, the voltage ratio of the system is measured in different tests. The voltage across transmitter is measured as an input voltage then divided by the output voltage which is measured across receiver. The results of the transducer characteristics and the ultrasound system behaviour are discussed in chapter 4 and 5 of this thesis. Each piezoelectric transducer has several resonance frequencies in which its impedance has lower magnitude as compared to non-resonance frequencies. Among these resonance frequencies, just at one of those frequencies, the magnitude of the impedance is minimum. This resonance frequency is known as the main resonance frequency of the transducer. To attain higher efficiency and deliver more power to the ultrasound system, the transducer is usually excited at the main resonance frequency. Therefore, it is important to find out this frequency and other resonance frequencies. Hereof, a frequency detection method is proposed in this research which is discussed in chapter 2. An extended electrical model of the ultrasound transducer with multiple resonance frequencies consists of several RLC legs in parallel with a capacitor. Each RLC leg represents one of the resonance frequencies of the ultrasound transducer. At resonance frequency the inductor reactance and capacitor reactance cancel out each other and the resistor of this leg represents power conversion of the system at that frequency. This concept is shown in simulation and test results presented in chapter 4. To excite a high power ultrasound transducer, a high power signal is required. Multilevel converters are usually applied to generate a high power signal but the drawback of this signal is low quality in comparison with a sinusoidal signal. In some applications like ultrasound, it is extensively important to generate a high quality signal. Several control and modulation techniques are introduced in different papers to control the output voltage of the multilevel converters. One of those techniques is harmonic elimination technique. In this technique, switching angles are chosen in such way to reduce harmonic contents in the output side. It is undeniable that increasing the number of the switching angles results in more harmonic reduction. But to have more switching angles, more output voltage levels are required which increase the number of components and cost of the converter. To improve the quality of the output voltage signal with no more components, a new harmonic elimination technique is proposed in this research. Based on this new technique, more variables (DC voltage levels and switching angles) are chosen to eliminate more low order harmonics compared to conventional harmonic elimination techniques. In conventional harmonic elimination method, DC voltage levels are same and only switching angles are calculated to eliminate harmonics. Therefore, the number of eliminated harmonic is limited by the number of switching cycles. In the proposed modulation technique, the switching angles and the DC voltage levels are calculated off-line to eliminate more harmonics. Therefore, the DC voltage levels are not equal and should be regulated. To achieve this aim, a DC/DC converter is applied to adjust the DC link voltages with several capacitors. The effect of the new harmonic elimination technique on the output quality of several single phase multilevel converters is explained in chapter 3 and 6 of this thesis. According to the electrical model of high power ultrasound transducer, this device can be modelled as parallel combinations of RLC legs with a main capacitor. The impedance diagram of the transducer in frequency domain shows it has capacitive characteristics in almost all frequencies. Therefore, using a voltage source converter to drive a high power ultrasound transducer can create significant leakage current through the transducer. It happens due to significant voltage stress (dv/dt) across the transducer. To remedy this problem, LC filters are applied in some applications. For some applications such as ultrasound, using a LC filter can deteriorate the performance of the transducer by changing its characteristics and displacing the resonance frequency of the transducer. For such a case a current source converter could be a suitable choice to overcome this problem. In this regard, a current source converter is implemented and applied to excite the high power ultrasound transducer. To control the output current and voltage, a hysteresis control and unipolar modulation are used respectively. The results of this test are explained in chapter 7.

Monte Carlo verification of intensity modulated radiotherapy treatments using EPIDS

Purpose: Electronic Portal Imaging Devices (EPIDs) are available with most linear accelerators (Amonuk, 2002), the current technology being amorphous silicon flat panel imagers. EPIDs are currently used routinely in patient positioning before radiotherapy treatments. There has been an increasing interest in using EPID technology tor dosimetric verification of radiotherapy treatments (van Elmpt, 2008). A straightforward technique involves the EPID panel being used to measure the fluence exiting the patient during a treatment which is then compared to a prediction of the fluence based on the treatment plan. However, there are a number of significant limitations which exist in this Method: Resulting in a limited proliferation ot this technique in a clinical environment. In this paper, we aim to present a technique of simulating IMRT fields using Monte Carlo to predict the dose in an EPID which can then be compared to the measured dose in the EPID. Materials: Measurements were made using an iView GT flat panel a-SI EPfD mounted on an Elekta Synergy linear accelerator. The images from the EPID were acquired using the XIS software (Heimann Imaging Systems). Monte Carlo simulations were performed using the BEAMnrc and DOSXVZnrc user codes. The IMRT fieids to be delivered were taken from the treatment planning system in DICOMRT format and converted into BEAMnrc and DOSXYZnrc input files using an in-house application (Crowe, 2009). Additionally. all image processing and analysis was performed using another in-house application written using the Interactive Data Language (IDL) (In Visual Information Systems). Comparison between the measured and Monte Carlo EPID images was performed using a gamma analysis (Low, 1998) incorporating dose and distance to agreement criteria. Results: The fluence maps recorded by the EPID were found to provide good agreement between measured and simulated data. Figure 1 shows an example of measured and simulated IMRT dose images and profiles in the x and y directions. "A technique for the quantitative evaluation of dose distributions", Med Phys, 25(5) May 1998 S. Crowe, 1. Kairn, A. Fielding, "The Development of a Monte Carlo system to verify Radiotherapy treatment dose calculations", Radiotherapy & Oncology, Volume 92, Supplement 1, August 2009, Pages S71-S71.

Beat-mixing Rock music: Rock and electronic dance music merge to create the Manarays

This research introduces the proposition that Electronic Dance Music’s beat-mixing function could be implemented to create immediacy in other musical genres. The inclusion of rhythmic sections at the beginning and end of each musical work created a ‘DJ friendly’ environment. The term used in this thesis to refer to the application of beat-mixing in Rock music is ‘ClubRock’. Collaboration between a number of DJs and Rock music professionals applied the process of beat-mixing to blend Rock tracks to produce a continuous ClubRock set. The DJ technique of beat-mixing Rock music transformed static renditions into a fluid creative work. The hybridisation of the two genres, EDM and Rock, resulted in a contribution to Rock music compositional approaches and the production of a unique Rock album; Manarays—Get Lucky.

Determining process model precision and generalization with weighted artificial negative events

Process mining encompasses the research area which is concerned with knowledge discovery from event logs. One common process mining task focuses on conformance checking, comparing discovered or designed process models with actual real-life behavior as captured in event logs in order to assess the “goodness” of the process model. This paper introduces a novel conformance checking method to measure how well a process model performs in terms of precision and generalization with respect to the actual executions of a process as recorded in an event log. Our approach differs from related work in the sense that we apply the concept of so-called weighted artificial negative events towards conformance checking, leading to more robust results, especially when dealing with less complete event logs that only contain a subset of all possible process execution behavior. In addition, our technique offers a novel way to estimate a process model’s ability to generalize. Existing literature has focused mainly on the fitness (recall) and precision (appropriateness) of process models, whereas generalization has been much more difficult to estimate. The described algorithms are implemented in a number of ProM plugins, and a Petri net conformance checking tool was developed to inspect process model conformance in a visual manner.

Investigation of the chemical and physical basis of oxidative stress generated by particulate matter using the profluorescent probe technique

Following the growing need for adoption of alternative fuels, this project aimed at getting more information on the oxidative potential of biodiesel particulate matter. Within this scope, the physical and chemical characteristics of biodiesel PM were analysed which lead to identification of reactive organic fractions. An in-house developed proflurescent nitroxide probe was used. This project further developed in-depth understanding of the chemical mechanisms following the detection of the oxidative potential of PM. This knowledge made a significant contribution to our understanding of processes behind negative health effects of pollution and enabled us to further develop new techniques to monitor it.

Developing an adaptive radiotherapy technique for virally mediated head and neck cancer

Purpose Virally mediated head and neck cancers (VMHNC) often present with nodal involvement and are highly radioresponsive, meaning that treatment plan adaptation during radiotherapy (RT) in a subset of patients is required. We sought to determine potential risk profiles and a corresponding adaptive treatment strategy for these patients. Methodology 121 patients with virally mediated, node positive nasopharyngeal (Epstein Barr Virus positive) or oropharyngeal (Human Papillomavirus positive) cancers, receiving curative intent RT were reviewed. The type, frequency and timing of adaptive interventions, including source-to-skin distance (SSD) corrections, re-scanning and re-planning, were evaluated. Patients were reviewed based on the maximum size of the dominant node to assess the need for plan adaptation. Results Forty-six patients (38%) required plan adaptation during treatment. The median fraction at which the adaptive intervention occurred was 26 for SSD corrections and 22 for re-planning CTs. A trend toward 3 risk profile groupings was discovered: 1) Low risk with minimal need (< 10%) for adaptive intervention (dominant pre-treatment nodal size of ≤ 35 mm), 2) Intermediate risk with possible need (< 20%) for adaptive intervention (dominant pre-treatment nodal size of 36 mm – 45 mm) and 3) High-risk with increased likelihood (> 50%) for adaptive intervention (dominant pre-treatment nodal size of ≥ 46 mm). Conclusion In this study, patients with VMHNC and a maximum dominant nodal size of > 46 mm were identified at a higher risk of requiring re-planning during a course of definitive RT. Findings will be tested in a future prospective adaptive RT study.

An investigation of the lived experiences of Intensive Care Paramedics

Intensive Care Paramedics (ICPs) attend to only the most clinically challenging of emergency medical cases, often working in a chaotic and frenetic atmosphere. They are regularly exposed to human tragedy and with that, the potential to experience traumatic events is not uncommon. There is very little known about the well-being of ICPs; how they cope with the demands of their role, or about their mental health in general. Nineteen experienced ICPs (4 female, 15 male) participated in a semi-structured interview. Themes were extracted from the data using an Interpretive Phenomenological Analysis approach. All participants discussed a work-related event they attended that traumatized them, usually experienced in the earlier parts of their career. Some spoke of an immediate overwhelming of their capacity to cope and others of a gradual onset of traumatic stress when reflecting on the event at a later time. More than half of the participants described events that involved children as the most difficult. Data revealed four superordinate themes: Social Support, Cognitive Coping, Proactive Coping, and Long Term Effects. Each superordinate theme comprised a number of constituent themes which are presented in this paper and exemplified with participant quotes. Although ongoing distress was described by some participants, all of the ICPs interviewed discussed positive aspects of their job; things that made the role worthwhile and fulfilling. This research highlights the important factors involved in coping with, and growing from, the extraordinary events that ICPs face. Results have implications for employing organizations and staff support services as well as for paramedics more broadly as they learn to cope with events inherent in their career. Findings indicate that positive adaptation and personal growth as a result of exposure to extremely high levels of potentially traumatic experiences is not only possible, but highly probable.

Producing and maintaining professional identities in early childhood

This study is an inquiry into early childhood teacher professional identities. In Australia, workforce reforms in early childhood include major shifts in qualification requirements that call for a university four-year degree-qualified teacher to be employed in child care. This marks a shift in the early years workforce, where previously there was no such requirement. At the same time as these reforms to quality measures are being implemented, and requiring a substantive up skilling of the workforce, there is a growing body of evidence through recent studies that suggests these same four-year degree-qualified early childhood teachers have an aversion to working in child care. Their preferred employment option is to work in the early years of more formal schooling, not in before-school contexts. This collision of agendas warrants investigation. This inquiry is designed to investigate the site at which advocacy for higher qualification requirements meets early childhood teachers who are reluctant to choose child care as a possible career pathway. The key research question for this study is: How are early childhood teachers’ professional identities currently produced? The work of this thesis is to problematise the early childhood teacher in child care through a particular method of discourse analysis. There are two sets of data. The first was a key early childhood political document that read as a "moment of arising" (Foucault, 1984a, p. 83). It is a political document which was selected for its current influence on the early childhood field, and in particular, workforce reforms that call for four-year degree-qualified teachers to work in before-school contexts, including child care. The second data set was generated through four focus group discussions conducted with preservice early childhood teachers. The document and transcripts of the focus groups were both analysed as text, as conceptualised by Foucault (1981). Foucault’s work spans a number of years and a range of philosophical matters. This thesis draws particularly on Foucault’s writings on discourse, power/knowledge, regimes of truth and resistance. In order to consider the production of early childhood teachers’ professional identities, the study is also informed by identity theorists, who have worked on gender, performativity and investment (Davies, 2004/2006; McNay, 1992; Osgood, 2012; Walkerdine, 1990; Weedon, 1997). The ways in which discourses intersect, compete and collide produce the subject (Foucault, 1981) and, in the case of this inquiry, there are a number of competing discourses at play, which produce the early childhood teacher. These particular theories turn particular lenses on the question of professional identities in early childhood, and such a study calls for the application of particular methodologies. Discourse analysis was used as the methodological framework, and the analysis was informed by Foucauldian concepts of discourse. While Foucault did not prescribe a form of discourse analysis as a method, his writings nonetheless provide a valuable framework for illuminating discursive practices and, in turn, how people are affected, through the shifts and distribution of power (Foucault, 1980a). The treatment used with both data sets involved redescription. For the policy document, a technique for reading document-as-text applied a genealogical approach (Foucault, 1984a). For the focus groups, the process of redescription (Rorty, 1989) involved reading talk-as-text. As a method, redescription involves describing "lots and lots of things in new ways until you have created a pattern of linguistic behaviour which will tempt the new generation to adopt it" (Rorty, 1989, p. 9). The development and application of categories (Davies, 2004/2006) built on a poststructuralist theoretical framework and the literature review informed the data analysis method of discourse analysis. Irony provided a rhetorical and playful tool (Haraway, 1991; Rorty, 1989), to look to how seemingly opposing discourses are held together. This opens a space to collapse binary thinking and consider seemingly contradictory terms in a way in which both terms are possible and both are true. Irony resists the choice of one or the other being right, and holds the opposites together in tension. The thesis concludes with proposals for new, ironic categories, which work to bring together seemingly opposing terms, located at sites in the field of early childhood where discourses compete, collide and intersect to produce and maintain early childhood teacher professional identities. The process of mapping these discourses goes some way to investigating the complexities about identities and career choices of early childhood teachers. The category of "the cost of loving" captures the collision between care/love, inherent in child care, and new discourses of investment/economics. Investment/economics has not completely replaced care/love, and these apparent opposites were not read as a binary because both are necessary and both are true (Haraway, 1991). They are held together in tension to produce early childhood teacher professional identities. The policy document under scrutiny was New Directions, released in 2007 by the then opposition ALP leader, Kevin Rudd. The claim was made strongly that the "economic prosperity" of Australia relies on investment in early childhood. The arguments to invest are compelling and the neuroscience/brain research/child development together with economic/investment discourses demand that early childhood is funding is increased. The intersection of these discourses produces professional identities of early childhood teachers as a necessary part of the country’s economy, and thus, worthy of high status. The child care sector and work in child care settings are necessary, with children and the early childhood teacher playing key roles in the economy of the nation. Through New Directions it becomes sayable (Foucault, 1972/1989) that the work the early childhood teacher performs is legitimated and valued. The children are produced as "economic units". A focus on what children are able to contribute to the future economy of the nation re-positions children and produces these "smart productive citizens", making future economic contribution. The early childhood teacher is produced through this image of a child and "the cost of loving" is emphasised. A number of these categories were produced through the readings of the document-as text and the talk-as-text. Two ironic categories were read in the analysis of the transcripts of the focus group discussions, when treated as talk-as-text data: the early childhood teacher as a "heroic victim"; and the early childhood teacher as a "glorified babysitter". This thesis raises new questions about professional identities in early childhood. These new questions might go some way to prompt re-thinking of some government policy, as well as some aspects of early childhood teacher education course design. The images of children and images of child care provide provocations to consider preservice teacher education course design. In particular, how child care, as one of the early childhood contexts, is located, conceptualised and spoken throughout the course. Consideration by course designers and teacher educators of what discourses are privileged in course content —what discourses are diminished or silenced—would go some way to reconceptualising child care within preservice teacher education and challenging dominant ways of speaking child care, and work in child care. This inquiry into early childhood teachers’ professional identities has gone some way to exploring the complexities around the early childhood teacher in child care. It is anticipated that the significance of this study will thus have immediate applicably and relevance for the Australian early childhood policy landscape. The early childhood field is in a state of rapid change, and this inquiry has examined some of the disconnects between policy and practice. Awareness of the discourses that are in play in the field will continue to allow space for conversations that challenge dominant assumptions about child care, work in child care and ways of being an early childhood teacher in child care.

Digitally-implemented naturally sampled PWM suitable for multilevel converter control

For dynamic closed loop control of a multilevel converter with a low pulse number (ratio of switching frequency to synthesized fundamental), natural sampled pulse-width modulation (PWM) is the best form of modulation. Natural sampling does not introduce distortion or a delayed response to the modulating signal. However previous natural sampled PWM implementations have generally been analog. For a modular multilevel converter, a digital implementation has advantages of accuracy and flexibility. Re-sampled uniform PWM is a novel digital modulation technique which approaches the performance of natural PWM. Both hardware and software implementations for a five level multilevel converter phase are presented, demonstrating the improvement over uniform PWM.

Resequencing and fine-mapping of the chromosome 12q13-14 locus associated with multiple sclerosis refines the number of implicated genes

Multiple sclerosis (MS) is a common chronic inflammatory disease of the central nervous system. Susceptibility to the disease is affected by both environmental and genetic factors. Genetic factors include haplotypes in the histocompatibility complex (MHC) and over 50 non-MHC loci reported by genome-wide association studies. Amongst these, we previously reported polymorphisms in chromosome 12q13-14 with a protective effect in individuals of European descent. This locus spans 288 kb and contains 17 genes, including several candidate genes which have potentially significant pathogenic and therapeutic implications. In this study, we aimed to fine-map this locus. We have implemented a two-phase study: a variant discovery phase where we have used next-generation sequencing and two target-enrichment strategies [long-range polymerase chain reaction (PCR) and Nimblegen's solution phase hybridization capture] in pools of 25 samples; and a genotyping phase where we genotyped 712 variants in 3577 healthy controls and 3269 MS patients. This study confirmed the association (rs2069502, P = 9.9 × 10−11, OR = 0.787) and narrowed down the locus of association to an 86.5 kb region. Although the study was unable to pinpoint the key-associated variant, we have identified a 42 (genotyped and imputed) single-nucleotide polymorphism haplotype block likely to harbour the causal variant. No evidence of association at previously reported low-frequency variants in CYP27B1 was observed. As part of the study we compared variant discovery performance using two target-enrichment strategies. We concluded that our pools enriched with Nimblegen's solution phase hybridization capture had better sensitivity to detect true variants than the pools enriched with long-range PCR, whilst specificity was better in the long-range PCR-enriched pools compared with solution phase hybridization capture enriched pools; this result has important implications for the design of future fine-mapping studies.

Remote measurement of diesel locomotive emission factors and particle size distributions

A technique for analysing exhaust emission plumes from unmodified locomotives under real world conditions is described and applied to the task of characterizing plumes from railway trains servicing an Australian shipping port. The method utilizes the simultaneous measurement, downwind of the railway line, of the following pollutants; particle number, PM2.5 mass fraction, SO2, NOx and CO2, with the last of these being used as an indicator of fuel combustion. Emission factors are then derived, in terms of number of particles and mass of pollutant emitted per unit mass of fuel consumed. Particle number size distributions are also presented. The practical advantages of the method are discussed including the capacity to routinely collect emission factor data for passing trains and to thereby build up a comprehensive real world database for a wide range of pollutants. Samples from 56 train movements were collected, analyzed and presented. The quantitative results for emission factors are: EF(N)=(1.7±1)×1016 kg-1, EF(PM2.5)= (1.1±0.5) g·kg-1, EF(NOx)= (28±14) g·kg-1, and EF(SO2 )= (1.4±0.4) g·kg-1. The findings are compared with comparable previously published work. Statistically significant (p<α, α=0.05) correlations within the group of locomotives sampled were found between the emission factors for particle number and both SO2 and NOx.