Computationally efficient solution techniques for adsorption problems involving steep gradients in bidisperse particles

A piecewise uniform fitted mesh method turns out to be sufficient for the solution of a surprisingly wide variety of singularly perturbed problems involving steep gradients. The technique is applied to a model of adsorption in bidisperse solids for which two fitted mesh techniques, a fitted-mesh finite difference method (FMFDM) and fitted mesh collocation method (FMCM) are presented. A combination (FMCMD) of FMCM and the DASSL integration package is found to be most effective in solving the problems. Numerical solutions (FMFDM and FMCMD) were found to match the analytical solution when the adsorption isotherm is linear, even under conditions involving steep gradients for which global collocation fails. In particular, FMCMD is highly efficient for macropore diffusion control or micropore diffusion control. These techniques are simple and there is no limit on the range of the parameters. The techniques can be applied to a variety of adsorption and desorption problems in bidisperse solids with non-linear isotherm and for arbitrary particle geometry.

Magma Flow Instabilities in a Volcanic Conduit: Implications for Long-Period Seismicity

Silicic volcanic eruptions are typically accompanied by repetitive Long-Period (LP) seismicity that originates from a small region of the upper conduit. These signals have the capability to advance eruption prediction, since they commonly precede a change in the eruption vigour. Shear bands forming along the conduit wall, where the shear stresses are highest, have been linked to providing the seismic trigger. However, existing computational models are unable to generate shear bands at the depths where the LP signals originate using simple magma strength models. Presented here is a model in which the magma strength is determined from a constitutive relationship dependent upon crystallinity and pressure. This results in a depth-dependent magma strength, analogous to planetary lithospheres. Hence, in shallow highly-crystalline regions a macroscopically discontinuous brittle type of deformation will prevail, whilst in deeper crystal-poor regions there will be a macroscopically continuous plastic deformation mechanism. This will result in a depth where the brittle-ductile transition occurs, and here shear bands disconnected from the free-surface may develop. We utilize the Finite Element Method and use axi-symmetric coordinates to model magma flow as a viscoplastic material, simulating quasi-static shear bands along the walls of a volcanic conduit. Model results constrained to the Soufrière Hills Volcano, Montserrat, show the generation of two types of shear bands: upper-conduit shear bands that form between the free-surface to a few 100 metres below it and discrete shear bands that form at the depths where LP seismicity is measured to occur corresponding to the brittle-ductile transition and the plastic shear region. It is beyond the limitation of the model to simulate a seismic event, although the modelled viscosity within the discrete shear bands suggests a failure and healing cycle time that supports the observed LP seismicity repeat times. However, due to the paucity of data and large parameter space available these results can only be considered to be qualitative rather than quantitative at this stage.

Solution of adsorption problems involving steep moving profiles

The moving finite element collocation method proposed by Kill et al. (1995) Chem. Engng Sci. 51 (4), 2793-2799 for solution of problems with steep gradients is further developed to solve transient problems arising in the field of adsorption. The technique is applied to a model of adsorption in solids with bidisperse pore structures. Numerical solutions were found to match the analytical solution when it exists (i.e. when the adsorption isotherm is linear). The method is simple yet sufficiently accurate for use in adsorption problems, where global collocation methods fail. (C) 1998 Elsevier Science Ltd. All rights reserved.

Thickness and stability of adsorbed film in cylindrical mesopores

The adsorbed film in small cylindrical mesopores is studied by using MCM-41 samples of uniform cylindrical channels as model systems. It is found that at a given relative pressure, the smaller the pore radius, the thicker the adsorbed film is, as postulated by Broekhoff and De Beer. Thermodynamics analysis established that the stability of the adsorbed film is determined by interface curvature and the potential of interaction between adsorbate and adsorbent. A semiempirical equation is proposed to describe the state of stable adsorbed films in cylindrical mesopores. It is also shown to be useful in calculations of pore size distributions of mesoporous solids.

Surface oxide and the role of magnesium during the sintering of aluminum

The effect of trace additions of magnesium on the sintering of aluminum and its alloys is examined. Magnesium, especially at low concentrations, has a disproportionate effect on sintering because it disrupts the passivating Al2O3 layer through the formation of a spinel phase. Magnesium penetrates the sintering compact by solid-state diffusion, and the oxide is reduced at the metal-oxide interface. This facilitates solid-state sintering, as well as wetting of the underlying metal by sintering liquids, when these are present. The optimum magnesium concentration is approximately 0.1 to 1.0 wt pet, but this is dependent on the volume of oxide and, hence, the particle size, as well as the sintering conditions. Small particle-size fractions require proportionally more magnesium than large-size fractions do.

Modeling the optical constants of GaP, InP, and InAs

An extension of the Adachi model with the adjustable broadening function, instead of the Lorentzian one, is employed to model the optical constants of GaP, InP, and InAs. Adjustable broadening is modeled by replacing the damping constant with the frequency-dependent expression. The improved flexibility of the model enables achieving an excellent agreement with the experimental data. The relative rms errors obtained for the refractive index equal 1.2% for GaP, 1.0% for InP, and 1.6% for InAs. (C) 1999 American Institute of Physics. [S0021-8979(99)05807-7].

Determination of chemical shielding tensor of an indole carbon and application of tryptophan orientation of a membrane peptide

Using tryptophan C-13-enriched at the C-4 (C epsilon(3)) of the indole, the orientation of the C epsilon(3) chemical shift tensor relative to the C epsilon(3)-H dipolar axis was determined from the C-13 chemical shift/C-13-H-1 dipolar 2D NMR powder pattern. The principal values obtained were 208, 137 and 15 ppm with sigma(33) perpendicular to the indole plane, and sigma(11) (least shielded direction) 5 degrees off the C epsilon(3)-H bond toward C xi(3). The side off the C epsilon(3)-H bond was determined by comparing the reduced chemical shift anisotropies obtained by solid-state NMR and from molecular dynamics calculations of [4-C-13] tryptophans in gramicidin A aligned in phospholipid membranes. (C) 1999 Elsevier Science B.V. All rights reserved.

Modeling the optical constants of AlxGa1-xAs alloys

The extension of Adachi's model with a Gaussian-like broadening function, in place of Lorentzian, is used to model the optical dielectric function of the alloy AlxGa1-xAs. Gaussian-like broadening is accomplished by replacing the damping constant in the Lorentzian line shape with a frequency dependent expression. In this way, the comparative simplicity of the analytic formulas of the model is preserved, while the accuracy becomes comparable to that of more intricate models, and/or models with significantly more parameters. The employed model accurately describes the optical dielectric function in the spectral range from 1.5 to 6.0 eV within the entire alloy composition range. The relative rms error obtained for the refractive index is below 2.2% for all compositions. (C) 1999 American Institute of Physics. [S0021-8979(99)00512-5].

The role of mobile belts for the longevity of deep cratonic lithosphere: The crumple zone model

Many Archean cratons are surrounded by Proterozoic mobile belts that have experienced episodes of tectonic re-activation over their lifetimes. This suggests that mobile belt lithosphere may be associated with long lived, inherited weakness. It is proposed that the proximity of this weakness can increase the longevity of deep Archean lithosphere by buffering Archean cratons from mantle derived stresses. The physical plausibility of this idea is explored through numerical simulations of mantle convection that include continents and allow for material rheologies that model the combined brittle and ductile behavior of the lithosphere. Within the simulations, the longevity of deep cratonic lithosphere does increase if it is buffered by mobile belts that can fail at relatively low stress levels.

Aluminophosphate-based mesoporous molecular sieves: synthesis and characterization of TAPOs

Mesoporous Ti-substituted aluminophosphates (AlPOs) with a hexagonal, cubic and lamellar pore structure, characteristic of MCM-41, MCM-48. and MCM-50, respectively, were synthesized. The stability of these mesophases upon template removal was studied. The pore structures, surface properties, and local atom environments of Al, P, and Ti of the hexagonal and cubic Ti-containing mesoporous products were extensively characterized using X-ray diffraction, magic angle spinning nuclear magnetic resonance, AAS, XPS, ultraviolet-visible, and adsorption of nitrogen and water vapor techniques while the lamellar mesophase was not further characterized due to its very poor thermal stability. Ti-containing mesoporous AlPO materials show a reasonable thermal stability upon template removal, a hydrophilic surface property, and high porosity showing application potentials in catalytic oxidation of hydrocarbons. (C) 2001 Elsevier Science B,V. All rights reserved.

Longevity and stability of cratonic lithosphere: Insights from numerical simulations of coupled mantle convection and continental tectonics

[1] The physical conditions required to provide for the tectonic stability of cratonic crust and for the relative longevity of deep cratonic lithosphere within a dynamic, convecting mantle are explored through a suite of numerical simulations. The simulations allow chemically distinct continents to reside within the upper thermal boundary layer of a thermally convecting mantle layer. A rheologic formulation, which models both brittle and ductile behavior, is incorporated to allow for plate-like behavior and the associated subduction of oceanic lithosphere. Several mechanisms that may stabilize cratons are considered. The two most often invoked mechanisms, chemical buoyancy and/or high viscosity of cratonic root material, are found to be relatively ineffective if cratons come into contact with subduction zones. High root viscosity can provide for stability and longevity but only within a thick root limit in which the thickness of chemically distinct, high-viscosity cratonic lithosphere exceeds the thickness of old oceanic lithosphere by at least a factor of 2. This end-member implies a very thick mechanical lithosphere for cratons. A high brittle yield stress for cratonic lithosphere as a whole, relative to oceanic lithosphere, is found to be an effective and robust means for providing stability and lithospheric longevity. This mode does not require exceedingly deep strength within cratons. A high yield stress for only the crustal or mantle component of the cratonic lithosphere is found to be less effective as detachment zones can then form at the crust-mantle interface which decreases the longevity potential of cratonic roots. The degree of yield stress variations between cratonic and oceanic lithosphere required for stability and longevity can be decreased if cratons are bordered by continental lithosphere that has a relatively low yield stress, i.e., mobile belts. Simulations that combine all the mechanisms can lead to crustal stability and deep root longevity for model cratons over several mantle overturn times, but the dominant stabilizing factor remains a relatively high brittle yield stress for cratonic lithosphere.