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.

Vertical sunshades to south-west elevation, Sunshine Coast University Library, Maroochydore

View of vertical sunshade to south-west elevation with glazing behind.

Effect of multiple transverse modes in self-mixing sensors based on vertical-cavity surface-emitting lasers

We investigate the effect of coexisting transverse modes on the operation of self-mixing sensors based on vertical-cavity surface-emitting lasers (VCSELs). The effect of multiple transverse modes on the measurement of displacement and distance were examined by simulation and in laboratory experiment. The simulation model shows that the periodic change in the shape and magnitude of the self-mixing signal with modulation current can be properly explained by the different frequency-modulation coefficients of the respective transverse modes in VCSELs. The simulation results are in excellent agreement with measurements performed on single-mode and multimode VCSELs and on self-mixing sensors based on these VCSELs.

An investigation of the effect of powder on the impact characteristics between a ball and a plate using free falling experiments

The study of the mechanisms of mechanical alloying requires knowledge of the impact characteristics between the ball and vial in the presence of milling powders. In this paper, foe falling experiments have br cn used to investigate the characteristics of impact events involved in mechanical milling. The effects of milling conditions, including impact velocity, ball size and powder thickness. on the coefficient of restitution and impact force are studied. It is found that the powder has a significant influence on the impact process due to its porous structure. This effect can be demonstrated using a modified Kelvin model. This study also confirms that the impact force is a relevant parameter for characterising the impact event due to its sensitivity to the milling conditions. (C) 1998 Elsevier Science S.A.

Optical properties of metallic films for vertical-cavity optoelectronic devices

We present models for the optical functions of 11 metals used as mirrors and contacts in optoelectronic and optical devices: noble metals (Ag, Au, Cu), aluminum, beryllium, and transition metals (Cr, Ni, Pd, Pt, Ti, W). We used two simple phenomenological models, the Lorentz-Drude (LD) and the Brendel-Bormann (BB), to interpret both the free-electron and the interband parts of the dielectric response of metals in a wide spectral range from 0.1 to 6 eV. Our results show that the BE model was needed to describe appropriately the interband absorption in noble metals, while for Al, Be, and the transition metals both models exhibit good agreement with the experimental data. A comparison with measurements on surface normal structures confirmed that the reflectance and the phase change on reflection from semiconductor-metal interfaces (including the case of metallic multilayers) can be accurately described by use of the proposed models for the optical functions of metallic films and the matrix method for multilayer calculations. (C) 1998 Optical Society of America.

Analysis of pore-fluid pressure gradient and effective vertical-stress gradient distribution in layered hydrodynamic systems

A theoretical analysis is carried out to investigate the pore-fluid pressure gradient and effective vertical-stress gradient distribution in fluid saturated porous rock masses in layered hydrodynamic systems. Three important concepts, namely the critical porosity of a porous medium, the intrinsic Fore-fluid pressure and the intrinsic effective vertical stress of the solid matrix, are presented and discussed. Using some basic scientific principles, we derive analytical solutions and explore the conditions under which either the intrinsic pore-fluid pressure gradient or the intrinsic effective vertical-stress gradient can be maintained at the value of the lithostatic pressure gradient. Even though the intrinsic pore-fluid pressure gradient can be maintained at the value of the lithostatic pressure gradient in a single layer, it is impossible to maintain it at this value in all layers in a layered hydrodynamic system, unless all layers have the same permeability and porosity simultaneously. However, the intrinsic effective vertical-stress gradient of the solid matrix can be maintained at a value close to the lithostatic pressure gradient in all layers in any layered hydrodynamic system within the scope of this study.

Evolution of three-dimensional folds for a non-Newtonian plate in a viscous medium

We derive analytical solutions for the three-dimensional time-dependent buckling of a non-Newtonian viscous plate in a less viscous medium. For the plate we assume a power-law rheology. The principal, axes of the stretching D-ij in the homogeneously deformed ground state are parallel and orthogonal to the bounding surfaces of the plate in the flat state. In the model formulation the action of the less viscous medium is replaced by equivalent reaction forces. The reaction forces are assumed to be parallel to the normal vector of the deformed plate surfaces. As a consequence, the buckling process is driven by the differences between the in-plane stresses and out of plane stress, and not by the in-plane stresses alone as assumed in previous models. The governing differential equation is essentially an orthotropic plate equation for rate dependent material, under biaxial pre-stress, supported by a viscous medium. The differential problem is solved by means of Fourier transformation and largest growth coefficients and corresponding wavenumbers are evaluated. We discuss in detail fold evolutions for isotropic in-plane stretching (D-11 = D-22), uniaxial plane straining (D-22 = 0) and in-plane flattening (D-11 = -2D(22)). Three-dimensional plots illustrate the stages of fold evolution for random initial perturbations or initial embryonic folds with axes non-parallel to the maximum compression axis. For all situations, one dominant set of folds develops normal to D-11, although the dominant wavelength differs from the Biot dominant wavelength except when the plate has a purely Newtonian viscosity. However, in the direction parallel to D-22, there exist infinitely many modes in the vicinity of the dominant wavelength which grow only marginally slower than the one corresponding to the dominant wavelength. This means that, except for very special initial conditions, the appearance of a three-dimensional fold will always be governed by at least two wavelengths. The wavelength in the direction parallel to D-11 is the dominant wavelength, and the wavelength(s) in the direction parallel to D-22 is determined essentially by the statistics of the initial state. A comparable sensitivity to the initial geometry does not exist in the classic two-dimensional folding models. In conformity with tradition we have applied Kirchhoff's hypothesis to constrain the cross-sectional rotations of the plate. We investigate the validity of this hypothesis within the framework of Reissner's plate theory. We also include a discussion of the effects of adding elasticity into the constitutive relations and show that there exist critical ratios of the relaxation times of the plate and the embedding medium for which two dominant wavelengths develop, one at ca. 2.5 of the classical Biot dominant wavelength and the other at ca. 0.45 of this wavelength. We propose that herein lies the origin of parasitic folds well known in natural examples.

Quasispin graded-fermion formalism and gl(m vertical bar n)down arrow osp(m vertical bar n) branching rules

The graded-fermion algebra and quasispin formalism are introduced and applied to obtain the gl(m\n)down arrow osp(m\n) branching rules for the two- column tensor irreducible representations of gl(m\n), for the case m less than or equal to n(n > 2). In the case m < n, all such irreducible representations of gl(m\n) are shown to be completely reducible as representations of osp(m\n). This is also shown to be true for the case m=n, except for the spin-singlet representations, which contain an indecomposable representation of osp(m\n) with composition length 3. These branching rules are given in fully explicit form. (C) 1999 American Institute of Physics. [S0022-2488(99)04410-2].

Groundwater waves in a coastal aquifer: A new governing equation including vertical effects and capillarity

Groundwater waves, that is, water table fluctuations, are a natural phenomenon in coastal aquifers. They represent an important part of the interaction between the ocean and aquifer and affect the mass exchange between them. This paper presents a new groundwater wave equation. Because it includes the effects of vertical flows and capillarity, the new equation is applicable to both intermediate-depth aquifers and high-frequency waves. Compared with the wave equation derived by Nielsen ed al. [1997], the present equation provides a closer representation of groundwater waves. In particular, it predicts high-frequency water table fluctuations as observed in the field. A validation of the new equation has been carried out by comparing the analytical solutions to it with predictions from direct simulations using the numerical model SUTRA. The effects of various physical parameters and their relative importance are also discussed.

Level-one highest weight representations of U-q[g1(1 vertical bar 1)] and associated vertex operators

We study the level-one irreducible highest weight representations of U-q[gl(1\1)] and associated q-vertex operators. We obtain the exchange relations satisfied by these vertex operators. The characters and supercharacters associated with these irreducible representations are calculated'. (C) 2000 Published by Elsevier Science B.V. All rights reserved.

U-q[(sl(2)over-cap vertical bar 1)] vertex operators, screen currents, and correlation functions at an arbitrary level

Bosonized q-vertex operators related to the four-dimensional evaluation modules of the quantum affine superalgebra U-q[sl((2) over cap\1)] are constructed for arbitrary level k=alpha, where alpha not equal 0,-1 is a complex parameter appearing in the four-dimensional evaluation representations. They are intertwiners among the level-alpha highest weight Fock-Wakimoto modules. Screen currents which commute with the action of U-q[sl((2) over cap/1)] up to total differences are presented. Integral formulas for N-point functions of type I and type II q-vertex operators are proposed. (C) 2000 American Institute of Physics. [S0022-2488(00)00608-3].

Level-one highest weight representation of U-q[sl((N)over-cap vertical bar 1)] and Bosonization of the multicomponent Super t-J model

We study the level-one irreducible highest weight representations of the quantum affine superalgebra U-q[sl((N) over cap\1)], and calculate their characters and supercharacters. We obtain bosonized q-vertex operators acting on the irreducible U-q[sl((N) over cap\1)] modules and derive the exchange relations satisfied by the vertex operators. We give the bosonization of the multicomponent super t-J model by using the bosonized vertex operators. (C) 2000 American Institute of Physics. [S0022- 2488(00)00508-9].

Quantum spin ladder systems associated with su(2 vertical bar 2)

Two integrable quantum spin ladder systems will be introduced associated with the fundamental su(2 \2) solution of the Yang-Baxter equation. The first model is a generalized quantum Ising system with Ising rung interactions. In the second model the addition of extra interactions allows us to impose Heisenberg rung interactions without violating integrability. The existence of a Bethe ansatz solution for both models allows us to investigate the elementary excitations for antiferromagnetic rung couplings. We find that the first model does not show a gap whilst in the second case there is a gap for all positive values of the rung coupling.

On the minimality of the p-harmonic map - x/vertical bar x vertical bar : B-n -> Sn-1

We prove that for any real number p with 1 p less than or equal to n - 1, the map x/\x\ : B-n --> Sn-1 is the unique minimizer of the p-energy functional integral(Bn) \delu\(p) dx among all maps in W-1,W-p (B-n, Sn-1) with boundary value x on phiB(n).

Mantle convection modeling with viscoelastic/brittle lithosphere: Numerical methodology and plate tectonic modeling

The earth's tectonic plates are strong, viscoelastic shells which make up the outermost part of a thermally convecting, predominantly viscous layer. Brittle failure of the lithosphere occurs when stresses are high. In order to build a realistic simulation of the planet's evolution, the complete viscoelastic/brittle convection system needs to be considered. A particle-in-cell finite element method is demonstrated which can simulate very large deformation viscoelasticity with a strain-dependent yield stress. This is applied to a plate-deformation problem. Numerical accuracy is demonstrated relative to analytic benchmarks, and the characteristics of the method are discussed.