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.

Dynamics of the Indonesian seas circulation. Part II – The role of pressure head

The aim of this paper is to analyze the role of the pressure head, i.e., the difference of total pressure forces acting on the Indonesian seas waters from the western Pacific and the eastern Indian Ocean, in driving the Indonesian Throughflow (ITF) and in determining the total transport of the ITF. These questions have been discussed in the literature but no consensus has been reached. A regional model of the Indonesian seas circulation has been developed that properly resolves all major topographic features in the region. The results of model runs have been used to calculate all components of the overall momentum balance. The estimates disclose that the dynamical balance is primarily between the volume integrated Coriolis acceleration, pressure gradient and the area integral of local wind stress. It is shown that consideration of components of momentum balance in the direction of the outflow through the Indian Ocean port leads to the formulation of a diagnostic relation between total inflow transports due to the Mindanao and New Guinea Coastal Currents and the external pressure head, internal pressure head, bottom form stress, and area integrated wind stress. Based on this relation, it is concluded that the external pressure head is not the major driving force of the ITF, which is why there is no unique relation between the total transport of the ITF and the external pressure head. However, Wyrtki's suggestion to monitor the variability of the total transport of the ITF by measurement of the sea-surface-height difference between the western Pacific and the eastern Indian Ocean is validated.

Modelling Transition Due to Backward Facing Steps Using the Laminar Kinetic Energy Concept

Boundary layer transition estimation and modelling is essential for the design of many engineering products across many industries. In this paper, the Reynolds-averaged Navier–Stokes are solved in conjunction with three additional transport equations to model and predict boundary layer transition. The transition model (referred to as the kTkT–kLkL–ωω model) is based on the kk–ωω framework with an additional transport equation to incorporate the effects low-frequency flow oscillations in the form of a laminar kinetic energy (kLkL). Firstly, a number of rectifications are made to the original kTkT–kLkL–ωω framework in order to ensure an appropriate response to the free-stream turbulence level and to improve near wall predictions. Additionally, the model is extended to incorporate the capability to model transition due to surface irregularities in the form of backward-facing steps with maximum non-dimensional step sizes of approximately 1.5 times the local displacement thickness of the boundary layer where the irregularity is located (i.e k/δ∗⪅1.5k/δ∗⪅1.5) at upstream turbulence intensities in the range 0.01<Tu(%)<0.80.01<Tu(%)<0.8. A novel function is proposed to incorporate transition sensitivity due to aft-facing steps. This paper details the rationale behind the development of this new function and demonstrates its suitability for transition onset estimation on a flat plate at zero pressure gradient.

TAP Reactor Development and Application in the Kinetic Characterization of Catalysts for Heterogeneously Catalyzed Gas Phase Reactions

This manuscript describes the application and further development of the TAP technique in kinetic characterization of heterogeneous catalysis. The major application of TAP systems is to study mechanisms, kinetics and transport phenomena in heterogeneous catalysis, all of which is made possible by the sub-millisecond time resolution. Furthermore, the kinetic information obtained can be used to gain an insight into the mechanism occurring over the catalyst system. This is advantageous as heterogeneous catalysts with an improved efficiency can be developed as a result. TAP kinetic studies are carried out at low pressure (~1x10-7 mbar) and TAP pulses are sufficiently small (1013-1015 molecules) so as to maintain this low pressure. The use of a small number of molecules in comparison to the total number of active sites means the state of the catalyst remains relatively unchanged. The use of the low intensity pulses also makes the pressure gradient negligible and so allows the TAP reactor system to operate in the Knudsen Diffusion regime, where gas-gas reactions are eliminated. Hence only gas-catalyst reactions are investigated and, by the use of moment analysis of observed exit flow, rate constants of elementary steps of the reaction can be obtained.

In this manuscript, two attempts to further the TAP technique are reported. Firstly, the work undertaken at QUB to attempt to control the number of molecules of condensable reagents that can be pulsed during a TAP pulse experiment is disclosed. Secondly, a collaborative project with SAI Ltd Manchester is discussed in a separate chapter, where technical details and validation of a customised time of flight mass spectrometer (ToF MS) for the QUB TAP-1 system are reported. A collaborative project with Cardiff Catalysis Institute focusing on the study of CO oxidation over hopcalite catalysts is also reported. The analysis of the experimental results has provided an insight into the possible mechanism of the oxidation of CO over these catalysts. A correction function has also been derived which accounts for the adsorption of reactant molecules over inert materials that are used for the reactor packing in TAP experiments. This function was then applied to the selective reduction of O2 in a H2 rich ethene feed, so that more accurate TAP moment based analysis could be conducted.

Particle-in-cell simulation study of a lower-hybrid shock

The expansion of a magnetized high-pressure plasma into a low-pressure ambient medium is examined with particle-in-cell simulations. The magnetic field points perpendicular to the plasma's expansion direction and binary collisions between particles are absent. The expanding plasma steepens into a quasi-electrostatic shock that is sustained by the lower-hybrid (LH) wave. The ambipolar electric field points in the expansion direction and it induces together with the background magnetic field a fast E cross B drift of electrons. The drifting electrons modify the background magnetic field, resulting in its pile-up by the LH shock. The magnetic pressure gradient force accelerates the ambient ions ahead of the LH shock, reducing the relative velocity between the ambient plasma and the LH shock to about the phase speed of the shocked LH wave, transforming the LH shock into a nonlinear LH wave. The oscillations of the electrostatic potential have a larger amplitude and wavelength in the magnetized plasma than in an unmagnetized one with otherwise identical conditions. The energy loss to the drifting electrons leads to a noticeable slowdown of the LH shock compared to that in an unmagnetized plasma.

Simulation of a collisionless planar electrostatic shock in a proton-electron plasma with a strong initial thermal pressure change

The localized deposition of the energy of a laser pulse, as it ablates a solid target, introduces high thermal pressure gradients in the plasma. The thermal expansion of this laser-heated plasma into the ambient medium (ionized residual gas) triggers the formation of non-linear structures in the collisionless plasma. Here an electron-proton plasma is modelled with a particle-in-cell simulation to reproduce aspects of this plasma expansion. A jump is introduced in the thermal pressure of the plasma, across which the otherwise spatially uniform temperature and density change by a factor of 100. The electrons from the hot plasma expand into the cold one and the charge imbalance drags a beam of cold electrons into the hot plasma. This double layer reduces the electron temperature gradient. The presence of the low-pressure plasma modifies the proton dynamics compared with the plasma expansion into a vacuum. The jump in the thermal pressure develops into a primary shock. The fast protons, which move from the hot into the cold plasma in the form of a beam, give rise to the formation of phase space holes in the electron and proton distributions. The proton phase space holes develop into a secondary shock that thermalizes the beam.

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%.