Solution techniques for transport problems involving steep concentration gradients: application to noncatalytic fluid solid reactions

Some efficient solution techniques for solving models of noncatalytic gas-solid and fluid-solid reactions are presented. These models include those with non-constant diffusivities for which the formulation reduces to that of a convection-diffusion problem. A singular perturbation problem results for such models in the presence of a large Thiele modulus, for which the classical numerical methods can present difficulties. For the convection-diffusion like case, the time-dependent partial differential equations are transformed by a semi-discrete Petrov-Galerkin finite element method into a system of ordinary differential equations of the initial-value type that can be readily solved. In the presence of a constant diffusivity, in slab geometry the convection-like terms are absent, and the combination of a fitted mesh finite difference method with a predictor-corrector method is used to solve the problem. Both the methods are found to converge, and general reaction rate forms can be treated. These methods are simple and highly efficient for arbitrary particle geometry and parameters, including a large Thiele modulus. (C) 2001 Elsevier Science Ltd. All rights reserved.

The effect of channel geometry and wall boundary conditions on the formation of extrusion surface instabilities for LLDPE

It is believed that surface instabilities can occur during the extrusion of linear low density polyethylene due to high extensional stresses at the exit of the die. Local crack development can occur at a critical stress level when melt rupture is reached. This high extensional stress results from the rearrangement of the flow at the boundary transition between the wall exit and the free surface. The stress is highest at the extrudate surface and decreases into the bulk of the material. The location of the region where the critical level is reached can determine the amplitude of the extrudate surface distortion, This paper studies the effect of wall slip on the numerically simulated extensional stress level at the die exit and correlates this to the experimentally determined amplitude of the surface instability. The effect of die exit radius and die wall roughness on extrusion surface instabilities is also correlated to the exit stress level in the same way. Whereas full slip may completely suppress the surface instability, a reduction in the exit stress level and instability amplitude is also shown for a rounded die exit and a slight increase in instability is shown to result from a rough die wall. A surface instability map demonstrates how the shear rate for onset of extrusion surface instabilities can be predicted on the basis of melt strength measurements and simulated stress peaks at the exit of the die. (C) 2001 Elsevier Science B.V. All rights reserved.

Saving, productivity and national income: a discrete-time geometric framework

This paper proposes an alternative geometric framework for analysing the inter-relationship between domestic saving, productivity and income determination in discrete time. The framework provides a means of understanding how low saving economies like the United States sustained high growth rates in the 1990s whereas high saving Japan did not. It also illustrates how the causality between saving and economic activity runs both ways and that discrete changes in national output and income depend on both current and previous accumulation behaviour. The open economy analogue reveals how international capital movements can create external account imbalances that enhance income growth for both borrower and lender economies. (C) 2002 Elsevier Science B.V. All rights reserved.

Existence of multiple solutions for second-order discrete boundary value problems

We give conditions on f involving pairs of discrete lower and discrete upper solutions which lead to the existence of at least three solutions of the discrete two-point boundary value problem yk+1 - 2yk + yk-1 + f (k, yk, vk) = 0, for k = 1,..., n - 1, y0 = 0 = yn,, where f is continuous and vk = yk - yk-1, for k = 1,..., n. In the special case f (k, t, p) = f (t) greater than or equal to 0, we give growth conditions on f and apply our general result to show the existence of three positive solutions. We give an example showing this latter result is sharp. Our results extend those of Avery and Peterson and are in the spirit of our results for the continuous analogue. (C) 2002 Elsevier Science Ltd. All rights reserved.

Influence of wall stress and left ventricular geometry on the accuracy of Dobutamine stress Echocardiography

OBJECTIVES The goal of this study was to determine whether wall stress at rest and during stress could explain the influence of left ventricular (LV) morphology on the accuracy of dobutamine stress echocardiography (DSE). BACKGROUND The sensitivity of DSE appears to be reduced in patients with concentric remodeling, but the cause of this finding is unclear. METHODS We studied 161 patients without resting wall motion abnormalities who underwent DSE and coronary angiography. Patients were classified into four groups according to relative wan thickness (normal

Discrete teleportation protocol of continuum spectra field states

A discrete protocol for teleportation of superpositions of coherent states of optical-cavity fields is presented. Displacement and parity operators are unconventionally used in Bell-like measurement for field states.

The nonexistence of spurious solutions to discrete, two-point boundary value problems

We investigate difference equations which arise as discrete approximations to two-point boundary value problems for systems of second-order, ordinary differential equations. We formulate conditions under which all solutions to the discrete problem satisfy certain a priori bounds which axe independent of the step-size. As a result, the nonexistence of spurious solutions are guaranteed. Some existence and convergence theorems for solutions to the discrete problem are also presented. (C) 2002 Elsevier Science Ltd. All rights reserved.

Discrete element modelling of the influence of lifters on power draw of tumbling mills

Crushing and grinding are the most energy intensive part of the mineral recovery process. A major part of rock size reduction occurs in tumbling mills. Empirical models for the power draw of tumbling mills do not consider the effect of lifters. Discrete element modelling was used to investigate the effect of lifter condition on the power draw of tumbling mill. Results obtained with PFC3D code show that lifter condition will have a significant influence on the power draw and on the mode of energy consumption in the mill. Relatively high lifters will consume less power than low lifters, under otherwise identical conditions. The fraction of the power that will be consumed as friction will increase as the height of the lifters decreases. This will result in less power being used for high intensity comminution caused by the impacts. The fraction of the power that will be used to overcome frictional resistance is determined by the material's coefficient of friction. Based on the modelled results, it appears that the effective coefficient of friction for in situ mill is close to 0.1. (C) 2003 Elsevier Science Ltd. All rights reserved.

Applying discrete element modelling to vertical and horizontal shaft impact crushers

The PFC3D (particle flow code) that models the movement and interaction of particles by the DEM techniques was employed to simulate the particle movement and to calculate the velocity and energy distribution of collision in two types of impact crusher: the Canica vertical shaft crusher and the BJD horizontal shaft swing hammer mill. The distribution of collision energies was then converted into a product size distribution for a particular ore type using JKMRC impact breakage test data. Experimental data of the Canica VSI crusher treating quarry and the BJD hammer mill treating coal were used to verify the DEM simulation results. Upon the DEM procedures being validated, a detailed simulation study was conducted to investigate the effects of the machine design and operational conditions on velocity and energy distributions of collision inside the milling chamber and on the particle breakage behaviour. (C) 2003 Elsevier Ltd. All rights reserved.

Discrete element modelling of power draw of tumbling mills

The power required to operate large mills is typically 5-10 MW. Hence, optimisation of power consumption will have a significant impact on overall economic performance and environmental impact. Power draw modelling results using the discrete element code PFC3D have been compared with results derived from the widely used empirical Model of Morrell. This is achieved by calculating the power draw for a range of operating conditions for constant mill size and fill factor using two modelling approaches. fThe discrete element modelling results show that, apart from density, selection of the appropriate material damping ratio is critical for the accuracy of modelling of the mill power draw. The relative insensitivity of the power draw to the material stiffness allows selection of moderate stiffness values, which result in acceptable computation time. The results obtained confirm that modelling of the power draw for a vertical slice of the mill, of thickness 20% of the mill length, is a reliable substitute for modelling the full mill. The power draw predictions from PFC3D show good agreement with those obtained using the empirical model. Due to its inherent flexibility, power draw modelling using PFC3D appears to be a viable and attractive alternative to empirical models where necessary code and computer power are available.

The uniqueness of solutions to discrete, victor, two-point boundary value problems

Difference equations which may arise as discrete approximations to two-point boundary value problems for systems of second-order, ordinary differential equations are investigated and conditions are formulated under which solutions to the discrete problem are unique. Some existence, uniqueness implies existence, and convergence theorems for solutions to the discrete problem are also presented.

Sex- and race-related differences in cross-sectional geometry and bone density of the femoral mid-shaft in older adults

Background : Femoral shaft fracture incidence increases in older adults and is associated with low-energy trauma. Apart from bone density, the distribution and size of bone contributes to its strength. Aim : To examine if bone geometry and density of the femoral mid-shaft in older adults differs by sex and race, we studied 197 White women, 225 Black women, 242 White men, and 148 Black men aged 70-79 years participating in the Health, Aging, and Body Composition study; a prospective cohort study in the USA. A secondary purpose of the study was to examine the association of site-specific muscle and fat to bone geometry and density. Subjects and methods : Subjects were community-dwelling and reported no difficulty walking one-quarter of a mile or climbing stairs. Mid-femoral volumetric bone mineral density (vBMD, mg cm -3 ), total area (TA), cortical area (CA), medullary area (MA), cross-sectional moments of inertia (CSMI: I x , I y , J ), and muscle and fat areas (cm 2 ) were determined by computed tomography (CT; GE CT-9800, 10 mm slice thickness). Results : vBMD was greater in men than women with no difference by race ( p < 0.001). Bone areas and area moments of inertia were also greater in men than women ( p < 0.001), with Black women having higher values than White women for TA and CA. Standardizing geometric parameters for body size differences by dividing by powers of femur length did not negate the sex difference for TA and MA. Significant differences ( p < 0.05) among the four groups also remained for I x and J . Mid-thigh muscle area was an independent contributor to TA in all groups (Std beta = 0.181-0.351, p < 0.05) as well as CA in women (Std beta = 0.246-0.254, p < 0.01) and CSMI in White women (Std beta = 0.175-0.185, p < 0.05). Further, muscle area was a significant contributor to vBMD in Black women. Conclusion : These results indicate that bone geometry and density of the femoral diaphysis differs primarily by sex, rather than race, in older well-functioning adults. In addition, site-specific muscle area appears to have a potential contributory role to bone geometry parameters, especially in women.

Discrete cyanide-bridged mixed-valence Co/Fe complexes: Outer-sphere redox behaviour

The outer-sphere redox behaviour of a series of [LnCoIII-NCFeII(CN)(5)](-) (L-n = n-membered pentadentate aza-macrocycle) complexes have been studied as a function of pH and oxidising agent. All the dinuclear complexes show a double protonation process at pH approximate to 2 that produces a shift in their UV/Vis spectra. Oxidation of the different non-protonated and diprotonated complexes has been carried out with peroxodisulfate, and of the non-protonated complexes also with trisoxalatocobaltate(III). The results are in agreement with predictions from the Marcus theory. The oxidation of [Fe(phen)(3)](3+) and [IrCl6](2-) is too fast to be measured, although for the latter the transient observation of the process has been achieved at pH = 0. The study of the kinetics of the outer-sphere redox process, with the S2O82- and [Co(ox)(3)](3-) oxidants, has been carried out as a function of pH, temperature, and pressure. As a whole, the values found for the activation volumes, entropies, and enthalpies are in the following margins, for the diprotonated and non-protonated dinuclear complexes, respectively: DeltaV(not equal) from 11 to 13 and 15 to 20 cm(3) mol(-1); DeltaS(not equal) from 110 to 30 and -60 to -90 J K-1 mol(-1); DeltaH(not equal) from 115 to 80 and 50 to 65 kJ.mol(-1). The thermal activation parameters are clearly dominated by the electrostriction occurring on outer-sphere precursor formation, while the trends found for the values of the volume of activation indicate an important degree of tuning due to the charge distribution during the electron transfer process. The special arrangement on the amine ligands in the isomer trans[(L14CoNCFeII)-N-III(CN)(5)](-) accounts for important differences in solvent-assisted hydrogen bonding occurring within the outer-sphere redox process, as has been established in redox reactions of similar compounds. ((C) Wiley-VCH Verlag GmbH & Co. KGaA, 69451 Weinheim, Germany, 2003).