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.

Lanczos subspace time-independent wave packet calculations of S(D-1)+H-2 reactive scattering

In this paper. we present the results of quantum dynamical simulations of the S (D-1) + H-2 insertion reaction on a newly developed potential energy surface (J. Chem. Phys. 2001, 114, 320). State-to-state reaction probabilities. product state distributions, and initial-state resolved cumulative reaction probabilities from a given incoming reactant channel are obtained from a time-independent wave packet analysis, performed within a single Lanczos subspace. Integral reaction cross sections are then estimated by J-shifting method and compared with the results from molecular beam experiment and QCT calculations.

A director theory for visco-elastic folding instabilities in multilayered rock

A model for finely layered visco-elastic rock proposed by us in previous papers is revisited and generalized to include couple stresses. We begin with an outline of the governing equations for the standard continuum case and apply a computational simulation scheme suitable for problems involving very large deformations. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered beam under compression. We analyse folding up to 40% shortening. The standard continuum solution becomes unstable for extreme values of the shear/normal viscosity ratio. The instability is a consequence of the neglect of the bending stiffness/viscosity in the standard continuum model. We suggest considering these effects within the framework of a couple stress theory. Couple stress theories involve second order spatial derivatives of the velocities/displacements in the virtual work principle. To avoid C-1 continuity in the finite element formulation we introduce the spin of the cross sections of the individual layers as an independent variable and enforce equality to the spin of the unit normal vector to the layers (-the director of the layer system-) by means of a penalty method. We illustrate the convergence of the penalty method by means of numerical solutions of simple shears of an infinite layer for increasing values of the penalty parameter. For the shear problem we present solutions assuming that the internal layering is oriented orthogonal to the surfaces of the shear layer initially. For high values of the ratio of the normal-to the shear viscosity the deformation concentrates in thin bands around to the layer surfaces. The effect of couple stresses on the evolution of folds in layered structures is also investigated. (C) 2002 Elsevier Science Ltd. All rights reserved.

A methodology for current density calculations in high-frequency RF resonators

As nuclear magnetic resonance imaging and spectroscopy move inexorably toward higher field-strength magnets in search of improved signal-to-noise ratio, spectral resolution, and spatial resolution, the way in which radiofrequency (RF) probes are designed changes. At higher frequencies, resonant cavities become the favored RF ''coil'' type and may be built using streamline elements to reduce the inductance of the system. In modeling such systems, the quasi-static approach of assuming that current flows evenly in all conductor cross sections and that adjacent conductors do not affect each other becomes less reasonable. The proximity of RF conductors in resonators typically causes RF eddy currents to flow, whereby the current density in each rung is altered by the RF fields generated by nearby conductors. The proper understanding and prediction of how resonators will perform require a model of the current densities flowing in conducting sections, including all RF eddy current effects. Very few models of this type have been presented in the literature. This article presents an overview of one such model and of how it may be applied to a variety of resonators, both shielded and unshielded, circular, and elliptical, in cross section. Results are presented from a shielded head coil operating at 2 tesla. (C) 1997 John Wiley & Sons, Inc.

Comments on alternative calculations of the broadening of spectral lines of neutral sodium by H-atom collisions

With the exception of the sodium D-lines, recent calculations of line broadening cross sections for several multiplets of sodium by Leininger et al (Leininger T, Gadea F X and Dickinson A 2000 J. Phys. B: At. Mol. Opt. Phys. 33 1805) are in substantial disagreement with cross sections interpolated from the tables of Anstee and O'Mara (Anstee and O'Mara 1995 Mon. Not. R. Astron. Soc. 276 859) and Barklem and O'Mara (Barklem P S and O'Mara B J 1997 Mon. Not. R. Astron. Soc. 290 102). The discrepancy is as large as a factor of 3 for the 3p-4d multiplet. The two theories are tested by using the results of each to synthesize lines in the solar spectrum. It is found that generally the data from the theory of Anstee, Barklem and O'Mara produce the best match to the observed solar spectrum. It is found, using a simple model for reflection of the optical electron by the potential barrier between the two atoms, that the reflection coefficient is too large for avoided crossings with the upper states of subordinate lines to contribute to line broadening, supporting the neglect of avoided ionic crossings by Anstee, Barklem and O'Mara for these lines. The large discrepancies between the two sets of calculations is a result of an approximate treatment of avoided ionic crossings for these lines by Leininger et al (Leininger T, Gadea F X and Dickinson A 2000 J. Phys. B: At. Mol. Opt. Phys. 33 1805).

Catheter effects in organ perfusion experiments

In a typical isolated organ perfusion experiment, a substance is injected upstream of an organ and then collected at some distance downstream. To reach the organ from the injection site, and then from the organ to the collector, a solute passes through catheters, usually tubes with circular cross-sections. Catheters cause distortion to the concentration-time profile of the perfusion. In this paper, we analyse catheter distribution kinetics from a mathematical point of view, develop the function most suitable for modeling this distribution and successfully apply this function to experimental data. (C) 2002 Academic Press.

Chebyshev real wave packet propagation: H+O-2 (J=0) state-to-state reactive scattering calculations

In this paper we explore the relative performance of two recently developed wave packet methodologies for reactive scattering, namely the real wave packet Chebyshev domain propagation of Gray and Balint-Kurti [J. Chem. Phys. 108, 950 (1998)] and the Lanczos subspace wave packet approach of Smith [J. Chem. Phys. 116, 2354 (2002); Chem. Phys. Lett. 336, 149 (2001)]. In the former method, a modified Schrodinger equation is employed to propagate the real part of the wave packet via the well-known Chebyshev iteration. While the time-dependent wave packet from the modified Schrodinger equation is different from that obtained using the standard Schrodinger equation, time-to-energy Fourier transformation yields wave functions which differ only trivially by normalization. In the Lanczos subspace approach the linear system of equations defining the action of the Green operator may be solved via either time-dependent or time-independent methods, both of which are extremely efficient due to the simple tridiagonal structure of the Hamiltonian in the Lanczos representation. The two different wave packet methods are applied to three dimensional reactive scattering of H+O-2 (total J=0). State-to-state reaction probabilities, product state distributions, as well as initial-state-resolved cumulative reaction probabilities are examined. (C) 2002 American Institute of Physics.

Calculation of bound and resonance states of HO2 for nonzero total angular momentum

Bound and resonance states of HO2 have been calculated quantum mechanically by the Lanczos homogeneous filter diagonalization method [Zhang and Smith, Phys. Chem. Chem. Phys. 3, 2282 (2001); J. Chem. Phys. 115, 5751 (2001)] for nonzero total angular momentum J = 1,2,3. For lower bound states, agreement between the results in this paper and previous work is quite satisfactory; while for high lying bound states and resonances these are the first reported results. A helicity quantum number V assignment (within the helicity conserving approximation) is performed and the results indicate that for lower bound states it is possible to assign the V quantum numbers unambiguously, but for resonances it is impossible to assign the V helicity quantum numbers due to strong mixing. In fact, for the high-lying bound states, the mixing has already appeared. These results indicate that the helicity conserving approximation is not good for the resonance state calculations and exact quantum calculations are needed to accurately describe the reaction dynamics for HO2 system. Analysis of the resonance widths shows that most of the resonances are overlapping and the interferences between them lead to large fluctuations from one resonance to another. In accord with the conclusions from earlier J = 0 calculations, this indicates that the dissociation of HO2 is essentially irregular. (C) 2003 American Institute of Physics.

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.

Application of X-ray microtomography to numerical simulations of agglomerate breakage by distinct element method

Computer-aided tomography has been used for many years to provide significant information about the internal properties of an object, particularly in the medical fraternity. By reconstructing one-dimensional (ID) X-ray images, 2D cross-sections and 3D renders can provide a wealth of information about an object's internal structure. An extension of the methodology is reported here to enable the characterization of a model agglomerate structure. It is demonstrated that methods based on X-ray microtomography offer considerable potential in the validation and utilization of distinct element method simulations also examined.

Migration of crystals during the filling of semi-solid castings

In cold-chamber high-pressure die castings (HPDC), the microstructure consists of coarse externally solidified crystals (ESCs) that are commonly observed in the central region of cross sections. In the present work, controlled laboratory scale casting experiments have been conducted with particular emphasis on the flow and solidification conditions. An A356 aluminum alloy was used to produce castings by pouring semi-solid metal through a steel die. Microstructures similar to those encountered in HPDC have been produced and the resulting microstructure is found to depend on the melt and die temperature: (1) the fraction of ESCs determines the extent of migration to the central region; (2) a maximum packing determines the area fraction of ESCs in the center; and (3) the die temperature affects the position of the ESCs-a higher die temperature can induce a displaced ESC distribution. It is found that the n-figration of crystals to the central region requires a flow, which is constrained at all melt/die interfaces. Furthermore, potential lift mechanisms are discussed. An assessment of the Saffman lift force on individual particles shows it has no significant effect on the migration of ESCs.

Microhardness mapping and the hardness-yield strength relationship in high-pressure diecast magnesium alloy AZ91

Microhardness maps of cross-sections of high-pressure diecast test bars of AZ91 have been determined. Specimens with rectangular cross-sections, 1, 2 and 3 mm thick, or with a circular cross-section 6.4 mm in diameter, have been studied. The hardness is generally higher near the edges in all specimens, and more so near the corners of the rectangular specimens. The hardness at the center of the castings is generally lower, due to a coarser solidification microstructure and the concentration of porosity. The evidence confirms that the surface of the castings is harder than the core, but it does not support the concept of a skin with a sharp. and definable boundary. This harder layer is irregular in hardness and depth and is not equally hard on opposite sides of the casting. The mean hardness obtained by integrating the microhardness maps over the entire cross-section increased with decreasing thickness of the bars, and was found to be in good correlation with each bar's yield strength. (c) 2005 Elsevier B.V. All rights reserved.

White and gray matter development in human fetal, newborn and pediatric brains

Brain anatomy is characterized by dramatic growth from the end of the second trimester through the neonatal stage. The characterization of normal axonal growth of the white matter tracts has not been well-documented to date and could provide important clues to understanding the extensive inhomogeneity of white matter injuries in cerebral palsy (CP) patients. However, anatomical studies of human brain development during this period are surprisingly scarce and histology-based atlases have become available only recently. Diffusion tensor magnetic resonance imaging (DTMRI) can reveal detailed anatomy of white matter. We acquired diffusion tensor images (DTI) of postmortem fetal brain samples and in vivo neonates and children. Neural structures were annotated in two-dimensional (2D) slices, segmented, measured, and reconstructed three-dimensionally (3D). The growth status of various white matter tracts was evaluated on cross-sections at 19-20 gestational weeks, and compared with 0-month-old neonates and 5- to 6-year-old children. Limbic, commissural, association, and projection white matter tracts and gray matter structures were illustrated in 3D and quantitatively characterized to assess their dynamic changes. The overall pattern of the time courses for the development of different white matter is that limbic fibers develop first and association fibers last and commissural and projection fibers are forming from anterior to posterior part of the brain. The resultant DTNIRI-based 3D human brain data will be a valuable resource for human brain developmental study and will provide reference standards for diagnostic radiology of premature newborns. (c) 2006 Elsevier Inc. All rights reserved.

Spherical body formation in the spirochaete Brachyspira hyodysenteriae

When cultures of Brachyspira hyodysenteriae were grown under a wide range of in vitro conditions, at least 1% of the cells formed spherical bodies different to the normal helical form. This percentage increased considerably in aging cultures or following their incubation in caramelized media. Spherical body formation was initiated from a terminal localized swelling of the outer sheath followed by a retraction of the protoplasmic cylinder into the resulting swollen vesicle. As this occurred, the periplasmic flagella seemed to unwind from the protoplasmic cylinder. Once retracted, the protoplasmic cylinder was found to be wrapped in an organized manner around the inner surface of the membrane of the swollen vesicle. Although most were 2-3 mu m in diameter, some much larger spherical bodies (6-12 mu m diameter) were occasionally seen, with a corresponding increase in the visible number of peripheral protoplasmic cylinder cross-sections. Spherical bodies from older cultures did not contain protoplasmic cylinders arranged around the periphery, but instead were characterized by the presence of a centrally located, electron-dense body c. 0.5-0.8 mu m in diameter. Brachyspira hyodysenteriae spherical bodies differ in both their structural organization and probable method of formation from similar structures described in other spirochaete genera.

Effect of sample geometry on regression rate of the melting interface for carbon steel burned in oxygen

Promoted-ignition testing on carbon steel rods of varying cross-sectional area and shape was performed in high pressure oxygen to assess the effect of sample geometry on the regression rate of the melting interface. Cylindrical and rectangular geometries and three different cross sections were tested and the regression rates of the cylinders were compared to the regression rates of the rectangular samples at test pressures around 6.9 MPa. Tests were recorded and video analysis used to determine the regression rate of the melting interface by a new method based on a drop cycle which was found to provide a good basis for statistical analysis and provide excellent agreement to the standard averaging methods used. Both geometries tested showed the typical trend of decreasing regression rate of the melting interface with increasing cross-sectional area; however, it was shown that the effect of geometry is more significant as the sample's cross sections become larger. Discussion is provided regarding the use of 3.2-mm square rods rather than 3.2-mm cylindrical rods within the standard ASTM test and any effect this may have on the observed regression rate of the melting interface.