The collision branching process

We consider a branching model, which we call the collision branching process (CBP), that accounts for the effect of collisions, or interactions, between particles or individuals. We establish that there is a unique CBP, and derive necessary and sufficient conditions for it to be nonexplosive. We review results on extinction probabilities, and obtain explicit expressions for the probability of explosion and the expected hitting times. The upwardly skip-free case is studied in some detail.

The influence of collision energy and strain accumulation on the kinetics of mechanical alloying

The kinetics of mechanical alloying have been investigated by examining the effect that ball mass has on the rate at which titanium carbide forms from the elements. By varying the ball density while keeping the ball diameter and the charge ratio constant, the collision energy was independently controlled. Grinding media with a density from 3.8 g cm(-3) (agate) to 16.4 g cm(-3) (tungsten carbide) were used. The reaction rate increases exponentially with ball mass until a critical level is reached, which is determined by the induced temperature rise. Above this level, collisions of higher energy have no advantage. It is also shown that the reaction rate increases exponentially with the rate at which strain accumulates in the reactants. It is suggested that the strain accumulation rate in mechanically induced reactions is analogous to temperature in thermally induced chemical reactions.

Proton transport model in the ionosphere .1. Multistream approach of the transport equations

The suprathermal particles, electrons and protons, coming from the magnetosphere and precipitating into the high-latitude atmosphere are an energy source of the Earth's ionosphere. They interact with ambient thermal gas through inelastic and elastic collisions. The physical quantities perturbed by these precipitations, such as the heating rate, the electron production rate, or the emission intensities, can be provided in solving the kinetic stationary Boltzmann equation. This equation yields particle fluxes as a function of altitude, energy, and pitch angle. While this equation has been solved through different ways for the electron transport and fully tested, the proton transport is more complicated. Because of charge-changing reactions, the latter is a set of two-coupled transport equations that must be solved: one for protons and the other for H atoms. We present here a new approach that solves the multistream proton/hydrogen transport equations encompassing the collision angular redistributions and the magnetic mirroring effect. In order to validate our model we discuss the energy conservation and we compare with another model under the same inputs and with rocket observations. The influence of the angular redistributions is discussed in a forthcoming paper.

Kinematics of oblique collision and ramping inferred from microstructures and strain in middle crustal rocks, central Southern Alps, New Zealand

The hanging wall of the Alpine Fault near Franz Josef Glacier has been exhumed during the past similar to2-3 m.y. providing a sample of the ductilely deformed middle crust of a modem obliquely convergent orogen. Presently exposed rocks of the Pacific Plate are inferred to have undergone several phases of ductile deformation as they moved westward above a mid-crustal detachment. Initially they were transpressed across the outboard part of the orogen, resulting in oblate fabrics with a down-dip stretch. Later, they encountered the Alpine Fault, experiencing an oblique-slip backshearing on vertical planes. This escalator-like deformation tilted and thinned the incoming crust onto that crustal-scale oblique ramp. This style of hanging wall deformation may affect only the most rapidly uplifting, central part of the Southern Alps because of the low flexural rigidity of the crust in that region and its displacement over a relatively sharp ramp-angle at depth. A 3D transpressive flow affected mylonites locally near the fault, but their shear direction remained parallel to plate motion, ruling out ductile 'extrusion' as an important process in this orogen. Outside the mylonite zone, late Cenozoic shortening is inferred to be modest (30-40%), as measured from deformation of younger biotite grains. Oblique collision is dominated by translation on the Alpine Fault, and rocks migrate rapidly through the deforming zone, preventing the accumulation of large finite strains. Transpression may play a minor role in oblique collision. (C) 2001 Elsevier Science Ltd. All rights reserved.

Symmetry contaminations in reactive scattering through long-lived collision complexes

We extend our Lanczos subspace time-independent wave packet method [J. Chem. Phys. 116 (2002) 2354] to investigate the issue of symmetry contaminations for the challenging deep-well H + O-2 reaction. Our central objective is to address the issue of whether significant symmetry contamination can occur if a wavepacket initially possessing the correct O-O exchange symmetry is propagated over tens of thousands of recursive steps using a basis which does not explicitly enforce the correct symmetry, and if so how seriously this affects the results. We find that symmetry contamination does exist where the symmetry constraint is not explicitly enforced in the basis. While it affects individual resonances and the associated peak amplitudes, the overall shape of the more averaged quantities such as total reaction probabilities and vibrational branching ratios are not seriously affected. (C) 2004 Elsevier B.V. All rights reserved.

Photodissociation of benzene under collision-free conditions: An ab initio/Rice-Ramsperger-Kassel-Marcus study

The ab initio/Rice-Ramsperger-Kassel-Marcus (RRKM) approach has been applied to investigate the photodissociation mechanism of benzene at various wavelengths upon absorption of one or two UV photons followed by internal conversion into the ground electronic state. Reaction pathways leading to various decomposition products have been mapped out at the G2M level and then the RRKM and microcanonical variational transition state theories have been applied to compute rate constants for individual reaction steps. Relative product yields (branching ratios) for C6H5+H, C6H4+H-2, C4H4+C2H2, C4H2+C2H4, C3H3+C3H3, C5H3+CH3, and C4H3+C2H3 have been calculated subsequently using both numerical integration of kinetic master equations and the steady-state approach. The results show that upon absorption of a 248 nm photon dissociation is too slow to be observable in molecular beam experiments. In photodissociation at 193 nm, the dominant dissociation channel is H atom elimination (99.6%) and the minor reaction channel is H-2 elimination, with the branching ratio of only 0.4%. The calculated lifetime of benzene at 193 nm is about 11 mus, in excellent agreement with the experimental value of 10 mus. At 157 nm, the H loss remains the dominant channel but its branching ratio decreases to 97.5%, while that for H-2 elimination increases to 2.1%. The other channels leading to C3H3+C3H3, C5H3+CH3, C4H4+C2H2, and C4H3+C2H3 play insignificant role but might be observed. For photodissociation upon absorption of two UV photons occurring through the neutral hot benzene mechanism excluding dissociative ionization, we predict that the C6H5+H channel should be less dominant, while the contribution of C6H4+H-2 and the C3H3+C3H3, CH3+C5H3, and C4H3+C2H3 radical channels should significantly increase. (C) 2004 American Institute of Physics.

Determination of coalescence kernels for high-shear granulation using DEM simulations

Discrete element method (DEM) modeling is used in parallel with a model for coalescence of deformable surface wet granules. This produces a method capable of predicting both collision rates and coalescence efficiencies for use in derivation of an overall coalescence kernel. These coalescence kernels can then be used in computationally efficient meso-scale models such as population balance equation (PBE) models. A soft-sphere DEM model using periodic boundary conditions and a unique boxing scheme was utilized to simulate particle flow inside a high-shear mixer. Analysis of the simulation results provided collision frequency, aggregation frequency, kinetic energy, coalescence efficiency and compaction rates for the granulation process. This information can be used to bridge the gap in multi-scale modeling of granulation processes between the micro-scale DEM/coalescence modeling approach and a meso-scale PBE modeling approach.

Pair correlations in a finite-temperature 1D Bose gas

We calculate the two-particle local correlation for an interacting 1D Bose gas at finite temperature and classify various physical regimes. We present the exact numerical solution by using the Yang-Yang equations and Hellmann-Feynman theorem and develop analytical approaches. Our results draw prospects for identifying the regimes of coherent output of an atom laser, and of finite-temperature “fermionization” through the measurement of the rates of two-body inelastic processes, such as photoassociation.

Non-Equilibrium Reaction Rates in the Macroscopic Chemistry Method for DSMC Calculations

The Direct Simulation Monte Carlo (DSMC) method is used to simulate the flow of rarefied gases. In the Macroscopic Chemistry Method (MCM) for DSMC, chemical reaction rates calculated from local macroscopic flow properties are enforced in each cell. Unlike the standard total collision energy (TCE) chemistry model for DSMC, the new method is not restricted to an Arrhenius form of the reaction rate coefficient, nor is it restricted to a collision cross-section which yields a simple power-law viscosity. For reaction rates of interest in aerospace applications, chemically reacting collisions are generally infrequent events and, as such, local equilibrium conditions are established before a significant number of chemical reactions occur. Hence, the reaction rates which have been used in MCM have been calculated from the reaction rate data which are expected to be correct only for conditions of thermal equilibrium. Here we consider artificially high reaction rates so that the fraction of reacting collisions is not small and propose a simple method of estimating the rates of chemical reactions which can be used in the Macroscopic Chemistry Method in both equilibrium and non-equilibrium conditions. Two tests are presented: (1) The dissociation rates under conditions of thermal non-equilibrium are determined from a zero-dimensional Monte-Carlo sampling procedure which simulates ‘intra-modal’ non-equilibrium; that is, equilibrium distributions in each of the translational, rotational and vibrational modes but with different temperatures for each mode; (2) The 2-D hypersonic flow of molecular oxygen over a vertical plate at Mach 30 is calculated. In both cases the new method produces results in close agreement with those given by the standard TCE model in the same highly nonequilibrium conditions. We conclude that the general method of estimating the non-equilibrium reaction rate is a simple means by which information contained within non-equilibrium distribution functions predicted by the DSMC method can be included in the Macroscopic Chemistry Method.

Relative motions of Africa, Iberia and Europe during Alpine orogeny

A revised kinematic model for the motions of Africa and Iberia relative to Europe since the Middle Jurassic is presented in order to provide boundary conditions for Alpine-Mediterranean reconstructions. These motions were calculated using up-to-date kinematic data predominantly based on magnetic isochrons in the Atlantic Ocean and published by various authors during the last 15 years. It is shown that convergence of Africa with respect to Europe commenced during the Cretaceous Normal Superchron (CNS), between chrons MO and 34 (120-83 Ma). This motion was subjected to fluctuations in convergence rates characterised by two periods of relatively rapid convergence (during Late Cretaceous and Eocene-Oligocene times) that alternated with periods of slower convergence (during the Paleocene and since the Early Miocene). Distinct changes in plate kinematics are recognised in the motion of Iberia with respect to Europe, indicated by: (1) a Late Jurassic-Early Cretaceous left-lateral strike-slip motion; (2) Late Cretaceous convergence; (3) Paleocene quiescence; (4) a short period of right-lateral strike-slip motion; and (5) final Eocene-Oligocene convergence. Based on these results, it is speculated that a collisional episode in the Alpine orogeny at ca. 65 Ma resulted in a dramatic decrease in the relative plate motions and that a slower motion since the Early Miocene promoted extension in the Mediterranean back-arc basins. (C) 2002 Elsevier Science B.V. All rights reserved.

Prediction of impact forces in a vibratory ball mill using an inverse technique

To understand the dynamic mechanisms of the mechanical milling process in a vibratory mill, it is necessary to determine the characteristics of the impact forces associated with the collision events. However, it is difficult to directly measure the impact force in an operating mill. This paper describes an inverse technique for the prediction of impact forces from acceleration measurements on a vibratory ball mill. The characteristics of the vibratory mill have been investigated by the modal testing technique, and its system modes have been identified. In the modelling of the system vibration response to the impact forces, two modal equations have been used to describe the modal responses. The superposition of the modal responses gives rise to the total response of the system. A method based on an optimisation approach has been developed to predict the impact forces by minimising the difference between the measured acceleration of the vibratory ball mill and the predicted acceleration from the solution of the modal equations. The predicted and measured impact forces are in good agreement. Copyright (C) 1996 Elsevier Science Ltd.

Gas-surface interactions in the thermal and sub-thermal regime: A molecular dynamics study

Molecular dynamics simulations are used to study the interaction of low-energy Ar atoms with the Ni(001) surface, Angular scattering distributions, in and out of the plane of incidence, are investigated as a function of incident energy, angles of incidence, crystallographic orientation of the incident beam and surface temperature. The results show a clear transition to the structure scattering regime at around 2 eV. However, at lower energies, two sub-regimes are revealed by the simulations, Far energies up to 250 meV, scattering is mainly diffuse, and significant trapping on the surface is observed, At energies above this level, lobular patterns start to form and trapping decreases with the increase in energy, Generally, there is a weak temperature dependence, but variations in the angle of incidence and/or changes in the crystallographic direction, generate significant changes in the scattering patterns.