Difference equation approach to two-thermocouple sensor characterisation in constant velocity flow environments

Thermocouples are one of the most popular devices for temperature measurement due to their robustness, ease of manufacture and installation, and low cost. However, when used in certain harsh environments, for example, in combustion systems and engine exhausts, large wire diameters are required, and consequently the measurement bandwidth is reduced. This article discusses a software compensation technique to address the loss of high frequency fluctuations based on measurements from two thermocouples. In particular, a difference equation sDEd approach is proposed and compared with existing methods both in simulation and on experimental test rig data with constant flow velocity. It is found that the DE algorithm, combined with the use of generalized total least squares for parameter identification, provides better performance in terms of time constant estimation without any a priori assumption on the time constant ratios of the thermocouples.

Exploiting a priori time constant ratio information in difference equation two-thermocouple sensor characterisation

The characterization of thermocouple sensors for temperature measurement in varying-flow environments is a challenging problem. Recently, the authors introduced novel difference-equation-based algorithms that allow in situ characterization of temperature measurement probes consisting of two-thermocouple sensors with differing time constants. In particular, a linear least squares (LS) lambda formulation of the characterization problem, which yields unbiased estimates when identified using generalized total LS, was introduced. These algorithms assume that time constants do not change during operation and are, therefore, appropriate for temperature measurement in homogenous constant-velocity liquid or gas flows. This paper develops an alternative ß-formulation of the characterization problem that has the major advantage of allowing exploitation of a priori knowledge of the ratio of the sensor time constants, thereby facilitating the implementation of computationally efficient algorithms that are less sensitive to measurement noise. A number of variants of the ß-formulation are developed, and appropriate unbiased estimators are identified. Monte Carlo simulation results are used to support the analysis.

Ion-Atom/Argon - Calculation of ionization cross sections by fast ion impact for neutral target atoms ranging from hydrogen to argon

A FORTRAN 90 program is presented which calculates the total cross sections, and the electron energy spectra of the singly and doubly differential cross sections for the single target ionization of neutral atoms ranging from hydrogen up to and including argon. The code is applicable for the case of both high and low Z projectile impact in fast ion-atom collisions. The theoretical models provided for the program user are based on two quantum mechanical approximations which have proved to be very successful in the study of ionization in ion-atom collisions. These are the continuum-distorted-wave (CDW) and continuum-distorted-wave eikonal-initial-state (CDW-EIS) approximations. The codes presented here extend previously published. codes for single ionization of. target hydrogen [Crothers and McCartney, Comput. Phys. Commun. 72 (1992) 288], target helium [Nesbitt, O'Rourke and Crothers, Comput. Phys. Commun. 114 (1998) 385] and target atoms ranging from lithium to neon [O'Rourke, McSherry and Crothers, Comput. Phys. Commun. 131 (2000) 129]. Cross sections for all of these target atoms may be obtained as limiting cases from the present code. Title of program: ARGON Catalogue identifier: ADSE Program summary URL: http://cpc.cs.qub.ac.uk/cpc/summaries/ADSE Program obtainable from: CPC Program Library Queen's University of Belfast, N. Ireland Licensing provisions: none Computer for which the program is designed and others on which it is operable: Computers: Four by 200 MHz Pro Pentium Linux server, DEC Alpha 21164; Four by 400 MHz Pentium 2 Xeon 450 Linux server, IBM SP2 and SUN Enterprise 3500 Installations: Queen's University, Belfast Operating systems under which the program has been tested: Red-hat Linux 5.2, Digital UNIX Version 4.0d, AIX, Solaris SunOS 5.7 Compilers: PGI workstations, DEC CAMPUS Programming language used: FORTRAN 90 with MPI directives No. of bits in a word: 64, except on Linux servers 32 Number of processors used: any number Has the code been vectorized or parallelized? Parallelized using MPI No. of bytes in distributed program, including test data, etc.: 32 189 Distribution format: tar gzip file Keywords: Single ionization, cross sections, continuum-distorted-wave model, continuum- distorted-wave eikonal-initial-state model, target atoms, wave treatment Nature of physical problem: The code calculates total, and differential cross sections for the single ionization of target atoms ranging from hydrogen up to and including argon by both light and heavy ion impact. Method of solution: ARGON allows the user to calculate the cross sections using either the CDW or CDW-EIS [J. Phys. B 16 (1983) 3229] models within the wave treatment. Restrictions on the complexity of the program: Both the CDW and CDW-EIS models are two-state perturbative approximations. Typical running time: Times vary according to input data and number of processors. For one processor the test input data for double differential cross sections (40 points) took less than one second, whereas the test input for total cross sections (20 points) took 32 minutes. Unusual features of the program: none (C) 2003 Elsevier B.V All rights reserved.

Dissociative charge exchange and ionization of O-2 by fast H+ and O+ ions: Energetic ion interactions in Europa's oxygen atmosphere and neutral torus

Measurements of electron capture and ionization of O-2 molecules in collisions with H+ and O+ ions have been made over an energy range 10 - 100 keV. Cross sections for dissociative and nondissociative interactions have been separately determined using coincidence techniques. Nondissociative channels leading to O-2(+) product formation are shown to be dominant for both the H+ and the O+ projectiles in the capture collisions and only for the H+ projectiles in the ionization collisions. Dissociative channels are dominant for ionizing collisions involving O+ projectiles. The energy distributions of the O+ fragment products from collisions involving H+ and O+ have also been measured for the first time using time-of-flight methods, and the results are compared with those from other related studies. These measurements have been used to describe the interaction of the energetic ions trapped in Jupiter's magnetosphere with the very thin oxygen atmosphere of the icy satellite Europa. It is shown that the ionization of oxygen molecules is dominated by charge exchange plus ion impact ionization processes rather than photoionization. In addition, dissociation is predominately induced through excitation of electrons into high-lying repulsive energy states ( electronically) rather than arising from momentum transfer from knock-on collisions between colliding nuclei, which are the only processes included in current models. Future modeling will need to include both these processes.

Neutral and Charged 1-Butyl-3-methylimidazolium Triflate Clusters: Equilibrium Concentration in the Vapor Phase and Thermal Properties of Nanometric Droplets

Ground state energy, structure, and harmonic vibrational modes of 1-butyl-3-methylimidazolium triflate ([bmim][Tf]) clusters have been computed using an all-atom empirical potential model. Neutral and charged species have been considered up to a size (30 [bmim][Tf] pairs) well into the nanometric range. Free energy computations and thermodynamic modeling have been used to predict the equilibrium composition of the vapor phase as a function of temperature and density. The results point to a nonnegligible concentration of very small charged species at pressures (P ~ 0.01 Pa) and temperatures (T 600 K) at the boundary of the stability range of [bmim][Tf]. Thermal properties of nanometric neutral droplets have been investigated in the 0 T 700 K range. A near-continuous transition between a liquidlike phase at high T and a solidlike phase at low T takes place at T ~ 190 K in close correspondence with the bulk glass point Tg ~ 200 K. Solidification is accompanied by a transition in the dielectric properties of the droplet, giving rise to a small permanent dipole embedded into the solid cluster. The simulation results highlight the molecular precursors of several macroscopic properties and phenomena and point to the close competition of Coulomb and dispersion forces as their common origin.

Solving the advection-dispersion equation using the discrete puff particle method

A flexible, mass-conservative numerical technique for solving the advection-dispersion equation for miscible contaminant transport is presented. The method combines features of puff transport models from air pollution studies with features from the random walk particle method used in water resources studies, providing a deterministic time-marching algorithm which is independent of the grid Peclet number and scales from one to higher dimensions simply. The concentration field is discretised into a number of particles, each of which is treated as a point release which advects and disperses over the time interval. The dispersed puff is itself discretised into a spatial distribution of particles whose masses can be pre-calculated. Concentration within the simulation domain is then calculated from the mass distribution as an average over some small volume. Comparison with analytical solutions for a one-dimensional fixed-duration concentration pulse and for two-dimensional transport in an axisymmetric flow field indicate that the algorithm performs well. For a given level of accuracy the new method has lower computation times than the random walk particle method.

A Van der Pol-Mathieu equation for the dynamics of dust grain charge in dusty plasmas

The chaotic profile of dust grain dynamics associated with dust-acoustic oscillations in a dusty plasma is considered. The collective behaviour of the dust plasma component is described via a multi-fluid model, comprising Boltzmann distributed electrons and ions, as well as an equation of continuity possessing a source term for the dust grains, the dust momentum and Poisson's equations. A Van der Pol–Mathieu-type nonlinear ordinary differential equation for the dust grain density dynamics is derived. The dynamical system is cast into an autonomous form by employing an averaging method. Critical stability boundaries for a particular trivial solution of the governing equation with varying parameters are specified. The equation is analysed to determine the resonance region, and finally numerically solved by using a fourth-order Runge–Kutta method. The presence of chaotic limit cycles is pointed out.

Plasma dynamics in an inductively coupled magnetic neutral loop discharge

Magnetic neutral loop discharges (NLDs) can be operated at significantly lower pressures than conventional radio-frequency (rf) inductively coupled plasmas (ICPs). These low pressure conditions are favourable for technological applications, in particular anisotropic etching. An ICP–NLD has been designed providing excellent diagnostics access for detailed investigations of fundamental mechanisms. Spatially resolved Langmuir probe measurements have been performed in the plasma production region (NL region) as well as in the remote application region downstream from the NL region. Depending on the NL gradient two different operation modes have been observed exhibiting different opportunities for control of plasma uniformity. The efficient operation at comparatively low pressures results in ionization degrees exceeding 1%. In this regime neutral dynamics has to be considered and can influence neutral gas and process uniformity. Neutral gas depletion through elevated gas temperatures and high ionization rates have been quantified. At pressures above 0.1 Pa, gas heating is the dominant depletion mechanism. At lower pressures neutral gas is predominantly depleted through high ionization rates and rapid transport of ions by ambipolar diffusion along the magnetic field lines. Non-uniform profiles of the ionization rate can, therefore, result in localized neutral gas depletion and non-uniform processing. We have also investigated the electron dynamics within the radio-frequency cycle using phase resolved optical emission spectroscopy and Thomson scattering. In these measurements electron drift phenomena along the NL torus have been identified.

Neutral gas depletion mechanisms in dense low-temperature argon plasmas

Neutral gas depletion mechanisms are investigated in a dense low-temperature argon plasma-an inductively coupled magnetic neutral loop (NL) discharge. Gas temperatures are deduced from the Doppler profile of the 772.38 nm line absorbed by argon metastable atoms. Electron density and temperature measurements reveal that at pressures below 0.1 Pa, relatively high degrees of ionization (exceeding 1%) result in electron pressures, p(e) = kT(e)n(e), exceeding the neutral gas pressure. In this regime, neutral dynamics has to be taken into account and depletion through comparatively high ionization rates becomes important. This additional depletion mechanism can be spatially separated due to non-uniform electron temperature and density profiles (non-uniform ionization rate), while the gas temperature is rather uniform within the discharge region. Spatial profiles of the depletion of metastable argon atoms in the NL region are observed by laser induced fluorescence spectroscopy. In this region, the depletion of ground state argon atoms is expected to be even more pronounced since in the investigated high electron density regime the ratio of metastable and ground state argon atom densities is governed by the electron temperature, which peaks in the NL region. This neutral gas depletion is attributed to a high ionization rate in the NL zone and fast ion loss through ambipolar diffusion along the magnetic field lines. This is totally different from what is observed at pressures above 10 Pa where the degree of ionization is relatively low (

Spatial structures of plasma parameters in a magnetic neutral loop discharge

Spatial structures of plasma parameters in a radio-frequency inductively coupled magnetic neutral loop discharge are investigated under various parameter variations using spatially resolved Langmuir probe measurements. A strong coupling between the plasma production region, in the neutral loop (NL) plane, and the axially remote substrate region is observed. The two regions are connected through the separatrices and therefore, spatial profiles in the substrate region are strongly influenced by the plasma production region and the structure of the separatrices. The electron temperature in the plasma production region peaks in the centre of the NL while the maximum in electron density is shifted radially inwards due to diffusion. Details of the structures in both regions, the production region and the substrate region, are determined through the position of the NL and the gradient of the inhomogeneous magnetic field around the NL confinement region. Parameter combinations are found providing higher plasma densities and better uniformity than in common inductively coupled plasmas without applying an additional magnetic field. The uniformity can be further improved using temporal variations of the magnetic field structure.

A planar inductively coupled radio-frequency magnetic neutral loop discharge

A planar inductively coupled radio-frequency (rf) magnetic neutral loop discharge has been designed. It provides diagnostic access to both the main plasma production region as well as a remote plane for applications. Three coaxial coils are arranged to generate a specially designed inhomogeneous magnetic field structure with vanishing field along a ring in the discharge-the so-called neutral loop (NL). The plasma is generated by applying an oscillating rf electric field along the NL, induced through a four-turn, planar antenna operated at 13.56 MHz. Electron density and temperature measurements are performed under various parameter variations. Collisionless electron heating in the NL region allows plasma operation at comparatively low pressures, down to 10(-2) Pa, with a degree of ionization in the order of several per cent. Conventional plasma operation in inductive mode without applying the magnetic field is less efficient, in particular in the low pressure regime where the plasma cannot be sustained without magnetic fields.

Difference equation sensor characterization algorithms for two-thermocouple probes

The characterization of thermocouple sensors for temperature measurement in variable flow environments is a challenging problem. In this paper, novel difference equation-based algorithms are presented that allow in situ characterization of temperature measurement probes consisting of two-thermocouple sensors with differing time constants. Linear and non-linear least squares formulations of the characterization problem are introduced and compared in terms of their computational complexity, robustness to noise and statistical properties. With the aid of this analysis, least squares optimization procedures that yield unbiased estimates are identified. The main contribution of the paper is the development of a linear two-parameter generalized total least squares formulation of the sensor characterization problem. Monte-Carlo simulation results are used to support the analysis.