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.

Plasma ionization through wave-particle interaction in a capacitively coupled radio-frequency discharge

Phase resolved optical emission spectroscopy, with high temporal resolution, shows that wave-particle interactions play a fundamental role in sustaining capacitively coupled rf plasmas. The measurements are in excellent agreement with a simple particle-in-cell simulation. Excitation and ionization mechanisms are dominated by beam-like electrons, energized through the advancing and retreating electric fields of the rf sheath. The associated large-amplitude electron waves, driven by a form of two-stream instability, result in power dissipation through electron trapping and phase mixing. (c) 2007 American Institute of Physics.

Liquid monomer-powder particle interaction in acrylic bone cement

It is known that the method used to mix the liquid monomer and powder of PMMA bone cement influences the quality of the cement that is used in total joint replacements. Mixing theory indicates that the interaction between the liquid monomer and the powder is affected by a number of parameters, such as cement viscosity and degree of agitation, with this knowledge utilized in the design of cement mixing devices. Therefore, the objectives of this study were to: (i) obtain information on the interaction of the liquid monomer and the powder in the case of an PMMA bone cement, (ii) show how this knowledge can be applied to the design of an automated cement mixing device, and (iii) compare the porosity, bending modulus, and bending strength of one commercially-available cement prepared using the automated mixer and prepared using a conventional mixer that is in current clinical use. Experimental data indicated that increasing the velocity and decreasing the viscosity of the systems produced cement that improved mechanical properties, which may contribute to better mechanical integrity and, hence, reduced tendency for aseptic loosening, of cemented hip implants.

Plasma Ionization in Low-Pressure Radio-Frequency Discharges—Part II: Particle-in-Cell Simulation

Plasma ionization in the low-pressure operation regime ( $«$ 5 Pa) of RF capacitively coupled plasmas (CCPs) is governed by a complex interplay of various mechanisms, such as field reversal, sheath expansion, and wave–particle interactions. In a previous paper, it was shown that experimental observations in a hydrogen CCP operated at 13.56 MHz are qualitatively well described in a 1-D symmetrical particle-in-cell (PIC) simulation. In this paper, a spherical asymmetrical PIC simulation that is closer to the conditions of the highly asymmetrical experimental device is used to simulate a low-pressure neon CCP operated at 2 MHz. The results show a similar behavior, with pronounced ionization through field reversal, sheath expansion, and wave–particle interactions, and can be exploited for more accurate quantitative comparisons with experimental observations.

One-dimensional particle simulation of the filamentation instability: Electrostatic field driven by the magnetic pressure gradient force

Two counterpropagating cool and equally dense electron beams are modeled with particle-in-cell simulations. The electron beam filamentation instability is examined in one spatial dimension, which is an approximation for a quasiplanar filament boundary. It is confirmed that the force on the electrons imposed by the electrostatic field, which develops during the nonlinear stage of the instability, oscillates around a mean value that equals the magnetic pressure gradient force. The forces acting on the electrons due to the electrostatic and the magnetic field have a similar strength. The electrostatic field reduces the confining force close to the stable equilibrium of each filament and increases it farther away, limiting the peak density. The confining time-averaged total potential permits an overlap of current filaments with an opposite flow direction.