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.

Prediction of ionic liquid properties. II. Volumetric properties as a function of temperature and pressure

The density of ionic liquids (ILs) as a function of pressure and temperature has been modeled using a group contribution model. This model extends the calculations previously reported (Jacquemin et al. J. Chem. Eng. Data 2008) which used 4000 IL densities at 298.15 K and 600 IL densities as a function of temperature up to 423 K at 0.1 MPa to pressures up to 207 MPa by using described data in the literature and presented in this study. The densities of two different ionic liquids (butyltrimethylammonium bis(trifluoromethylsulfonyl)imide, [N][NTf], and 1-butyl-l-methyl-pyrrolidiniumbis(trifluoromethylsulfonyl)imide, [C mPyrro]-[NTf]) were measured as a function of temperature from (293 to 415) K and over an extended pressure range from (0.1 to 40) MPa using a vibrating-tube densimeter. The model is able to predict the ionic liquid densities of over 5080 experimental data points to within 0.36%. In addition, this methodology allows the calculation of the mechanical coefficients using the calculated density as a function of temperature and pressure with an estimated uncertainty of ± 20%. © 2008 American Chemical Society.

Liquid-liquid slug flow: Hydrodynamics and pressure drop

In this paper, the hydrodynamics and the pressure drop of liquid-liquid slug flow in round microcapillaries are presented. Two liquid-liquid flow systems are considered, viz. water-toluene and ethylene glycol/water-toluene. The slug lengths of the alternating continuous and dispersed phases were measured as a function of the slug velocity (0.03-0.5 m/s), the organic-to-aqueous flow ratio (0.1-4.0), and the microcapillary internal diameter (248 and 498 mu m). The pressure drop is modeled as the sum of two contributions: the frictional and the interface pressure drop. Two models are presented, viz, the stagnant film model and the moving film model. Both models account for the presence of a thin liquid film between the dispersed phase slug and the capillary wall. It is found that the film velocity is of negligible influence on the pressure drop. Therefore, the stagnant film model is adequate to accurately predict the liquid-liquid slug flow pressure drop. The influence of inertia and the consequent change of the slug cap curvature are accounted for by modifying Bretherton's curvature parameter in the interface pressure drop equation. The stagnant film model is in good agreement with experimental data with a mean relative error of less than 7%.

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.

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.

Dynamics of the Indonesian seas circulation. Part I – The influence of bottom topography on temperature and salinity distributions

The influence of bottom topography on the distribution of temperature and salinity in the Indonesian seas region has been studied with a high-resolution model based on the Princeton Ocean Model. One of the distinctive properties of the model is an adequate reproduction of all major topographic features in the region by the model bottom relief. The three major routes of flow of Pacific water through the region have been identified. The western route follows the flow of North Pacific Water through the Sulawesi Sea, Makassar Strait, Flores Sea, and Banda Sea. This is the main branch of the Indonesian Throughflow. The eastern routes follow the flow of South Pacific water through the eastern Indonesian seas. This water enters the region either through the Halmahera Sea or by flowing to the north around Halmahera Island into the Morotai Basin and then into the Maluku Sea. A deep southward flow of South Pacific Water fills the Seram Sea below 1200 m through the Lifamatola Passage. As it enters the Seram Sea, this overflow turns eastward at depths greater than 2000 m, then upwells in the eastern part of the Seram Sea before returning westward at ~1500-2000 m. The flow continues westward across the Seram Sea, spreading to greater depths before entering the Banda Sea at the Buru-Mangole passage. It is this water that shapes the temperature and salinity of the deep Banda Sea. Topographic elevations break the Indonesian seas region down into separate basins. The difference in the distributions of potential temperature, ?, and salinity, S, in adjacent basins is primarily due to specific properties of advection of ? and S across a topographic rise. By and large, the topographic rise blocks deep flow between basins whereas water shallower than the depth of the rise is free to flow between basins. To understand this process, the structure of simulated fields of temperature and salinity has been analyzed. To identify a range of advected ? or S, special sections over the sills with isotherms or isohalines and isotachs of normal velocity have been considered. Following this approach the impact of various topographic rises on the distribution of ? and S has been identified. There are no substantial structural changes of potential temperature and salinity distributions between seasons, though values of some parameters of temperature and salinity distributions, e.g., magnitudes of maxima and minima, can change. It is shown that the main structure of the observed distributions of temperature and salinity is satisfactorily reproduced by the model throughout the entire domain.

A review of rhizosphere carbon flow modelling

Rhizosphere processes play a key role in nutrient cycling in terrestrial ecosystems. Plant rhizodeposits supply low-molecular weight carbon substrates to the soil microbial community, resulting in elevated levels of activity surrounding the root. Mechanistic compartmental models that aim to model carbon flux through the rhizosphere have been reviewed and areas of future research necessary to better calibrate model parameters have been identified. Incorporating the effect of variation in bacterial biomass physiology on carbon flux presents a considerable challenge to experimentalists and modellers alike due to the difficulties associated with differentiating dead from dormant cells. A number of molecular techniques that may help to distinguish between metabolic states of bacterial cells are presented. The calibration of growth, death and maintenance parameters in rhizosphere models is also discussed. A simple model of rhizosphere carbon flow has been constructed and a sensitivity analysis was carried out on the model to highlight which parameters were most influential when simulating carbon flux. It was observed that the parameters that most heavily influenced long-term carbon compartmentalisation in the rhizosphere were exudation rate and biomass yield. It was concluded that future efforts to simulate carbon flow in the rhizosphere should aim to increase ecological realism in model structure.

Experimental investigation of flow and segregation behaviour of bulk solids in silos under high gravity conditions

Many researchers have investigated the flow and segregation behaviour in model scale experimental silos at normal gravity conditions. However it is known that the stresses experienced by the bulk solid in industrial silos are high when compared to model silos. Therefore it is important to understand the effect of stress level on flow and segregation behaviour and establish the scaling laws governing this behaviour. The objective of this paper is to understand the effect of gravity on the flow and segregation behaviour of bulk solids in a silo centrifuge model. The materials used were two mixtures composed of Polyamide and glass beads. The discharge of two bi-disperse bulk solids in a silo centrifuge model were recorded under accelerations ranging from 1g to 15g. The velocity distribution during discharge was evaluated using Particle Image Velocimetry (PIV) techniques and the concentration distribution of large and small particles were obtained by imaging processing techniques. The flow and segregation behaviour at high gravities were then quantified and compared with the empirical equations available in the literature.

Ischemic preconditioning before lower limb ischemia-reperfusion protects against acute lung injury

Objective: Prolonged limb ischemia followed by reperfusion (I/R) is associated with a systemic inflammatory response syndrome and remote acute lung injury. Ischemic preconditioning (IPC), achieved with repeated brief periods of I/R before the prolonged ischemic period, has been shown to protect skeletal muscle against ischemic injury. The aim of this study was to ascertain whether IPC of the limb before I/R injury also attenuates systemic inflammation and acute lung injury in a fully resuscitated porcine model of hind limb I/R. Methods: This prospective, randomized, controlled, experimental animal study was performed in a university-based animal research facility with 18 male Landrace pigs that weighed from 30 to 35 kg. Anesthetized ventilated swine were randomized (n = 6 per group) to three groups: sham-operated control group, I/R group (2 hours of bilateral hind limb ischemia and 2.5 hours of reperfusion), and IPC group (three cycles of 5 minutes of ischemia/5 minutes of reperfusion immediately preceding I/R). Plasma was separated and stored at -70° C for later determination of plasma tumor necrosis factor-a and interleukin-6 with bioassay as markers of systemic inflammation. Circulating phagocytic cell priming was assessed with a whole blood chemiluminescence assay. Lung tissue wet-to-dry weight ratio and myeloperoxidase concentration were markers of edema and neutrophil sequestration, respectively. The alveolar-arterial oxygen gradient and pulmonary artery pressure were indices of lung function. Results: In a porcine model, bilateral hind limb (I/R) injury significantly increased plasma interleukin-6 concentrations, circulating phagocytic cell priming, and pulmonary leukosequestration, edema, and impaired gas exchange. Conversely, pigs treated with IPC before the onset of the ischemic period had significantly reduced interleukin-6 levels, circulating phagocytic cell priming, and experienced significantly less pulmonary edema, leukosequestration, and respiratory failure. Conclusion: Lower limb IPC protects against systemic inflammation and acute lung injury in lower limb I/R injury.

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.

High pressure CO2 absorption studies on imidazolium-based ionic liquids:Experimental and simulation approaches

A combined experimental-computational study on the CO absorption on 1-butyl-3-methylimidazolium hexafluophosphate, 1-ethyl-3-methylimidazolium bis[trifluoromethylsulfonyl]imide, and 1-butyl-3-methylimidazolium bis[trifluoromethylsulfonyl]imide ionic liquids is reported. The reported results allowed to infer a detailed nanoscopic vision of the absorption phenomena as a function of pressure and temperature. Absorption isotherms were measured at 318 and 338K for pressures up to 20MPa for ultrapure samples using a state-of-the-art magnetic suspension densimeter, for which measurement procedures are developed. A remarkable swelling effect upon CO absorption was observed for pressures higher than 10MPa, which was corrected using a method based on experimental volumetric data. The experimental data reported in this work are in good agreement with available literature isotherms. Soave-Redlich-Kwong and Peng-Robinson equations of state coupled with bi-parametric van der Waals mixing rule were used for successful correlations of experimental high pressure absorption data. Molecular dynamics results allowed to infer structural, energetic and dynamic properties of the studied CO+ionic liquids mixed fluids, showing the relevant role of the strength of anion-cation interactions on fluid volumetric properties and CO absorption. © 2012 Elsevier B.V.

High-pressure phase equilibrium in the {carbon dioxide (1) + 1-chloropropane (2)} binary system

Abstract The current study reports original vapour-liquid equilibrium (VLE) for the system {CO2 (1) + 1-chloropropane (2)}. The measurements have been performed over the entire pressure-composition range for the (303.15, 313.15 and 328.15) K isotherms. The values obtained have been used for comparison of four predictive approaches, namely the equation of state (EoS) of Peng and Robinson (PR), the Soave modification of Benedict–Webb–Rubin (SBWR) EoS, the Critical Point-based Revised Perturbed-Chain Association Fluid Theory (CP-PC-SAFT) EoS, and the Conductor-like Screening Model for Real Solvents (COSMO-RS). It has been demonstrated that the three EoS under consideration yield similar and qualitatively accurate predictions of VLE, which is not the case for the COSMO-RS model examined. Although CP-PC-SAFT EoS exhibits only minor superiority in comparison with PR and SBWR EoS in predicting VLE in the system under consideration, its relative complexity can be justified when taking into account the entire thermodynamic phase space and, in particular, considering the liquid densities and sound velocities over a wider pressure-volume-temperature range.

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.

Optimisation of source/drain extension region profile for suppression of short channel effects in sub-50 nm DG SOI MOSFETs with high- k gate dielectrics

Novel technology dependent scaling parameters i.e. spacer to gradient ratio and effective channel length (Leff) are proposed for source/drain engineered DG MOSFET, and their significance in minimizing short channel effects (SCES) in high-k gate dielectrics is discussed in detail. Results show that a high-k dielectric should be associated with a higher spacer to gradient ratio to minimise SCEs The analytical model agrees with simulated data over the entire range of spacer widths, doping gradients, high-k gate dielectrics and effective channel lengths.