Application of Petrov-Galerkin methods to transient boundary value problems in chemical engineering: adsorption with steep gradients in bidisperse solids

Petrov-Galerkin methods are known to be versatile techniques for the solution of a wide variety of convection-dispersion transport problems, including those involving steep gradients. but have hitherto received little attention by chemical engineers. We illustrate the technique by means of the well-known problem of simultaneous diffusion and adsorption in a spherical sorbent pellet comprised of spherical, non-overlapping microparticles of uniform size and investigate the uptake dynamics. Solutions to adsorption problems exhibit steep gradients when macropore diffusion controls or micropore diffusion controls, and the application of classical numerical methods to such problems can present difficulties. In this paper, a semi-discrete Petrov-Galerkin finite element method for numerically solving adsorption problems with steep gradients in bidisperse solids is presented. The numerical solution was found to match the analytical solution when the adsorption isotherm is linear and the diffusivities are constant. Computed results for the Langmuir isotherm and non-constant diffusivity in microparticle are numerically evaluated for comparison with results of a fitted-mesh collocation method, which was proposed by Liu and Bhatia (Comput. Chem. Engng. 23 (1999) 933-943). The new method is simple, highly efficient, and well-suited to a variety of adsorption and desorption problems involving steep gradients. (C) 2001 Elsevier Science Ltd. All rights reserved.

Influence of seaward boundary condition on contaminant transport in unconfined coastal aquifers

Contaminant transport in coastal aquifers is complicated partly due to the conditions at the seaward boundary including seawater intrusion and tidal variations of sea level. Their inclusion in modelling this system will be computationally expensive. Therefore, it will be instructive to investigate the consequence of simplifying the seaward boundary condition by neglecting the seawater density and tidal variations in numerical predictions of contaminant transport in this zone. This paper presents a comparison of numerical predictions for a simplified seaward boundary condition with experimental results for a corresponding realistic one including a saltwater interface and tidal variations. Different densities for contaminants are considered. The comparison suggests that the neglect of the seawater intrusion and tidal variations does not affect noticeably the overall migration rate of the plume before it reaches the saltwater interface. However, numerical prediction shows that a more dense contaminant travels further seaward and part of the solute mass exits under the sea if the seawater density is not included. This is not consistent with the experimental result, which shows that the contaminant travels upwards to the shoreline along the saltwater interface. Neglect of seawater density, therefore, will result in an underestimation of the exit rate of solute mass around the coastline and fictitious migration paths under the seabed. For a less dense contaminant, neglect of seawater density has little effect on numerical prediction of migration paths. (C) 2001 Elsevier Science B.V. 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.

Scavenging of siliceous grain-boundary phase of 8-mol%-ytterbia-stabilized zirconia without additive

The grain-boundary conductivity (sigma (gb),) of 8-mol%-ytterbiastabilized zirconia increased markedly with heat treatment between 1000 degrees and 1300 degreesC with a slow heating rate (0.1 degreesC/min) before sintering. The extent of the sigma (gb) improvement was the same or larger than that via Al2O3 addition. The heat treatment did not affect the grain-interior conduction when sintered at 1600 degreesC, while Al2O3-derived scavenging significantly did, given the larger increment of total conductivity in the heat-treated sample. The formation of a silicon-containing phase in a discrete form was suggested as a possible route of scavenging the resistive phase from the correlation between average grain size and sigma (gb).

Pre- and postsynaptic actions of ATP on neurotransmission in rat submandibular ganglia

The pre- and postsynaptic actions of exogenously applied ATP were investigated in intact and dissociated parasympathetic neurotics of rat submandibular ganglia. Nerve-evoked excitatory postsynaptic potentials (EPSPs) were not inhibited by the purinergic receptor antagonists, suramin and pyridoxal-phosphate-6-azophenyl-2 ' ,4 ' -disulphonic acid (PPADS), or the desensitising agonist, alpha,beta -methylene ATP. In contrast. EPSPs were abolished by the nicotinic acetylcholine receptor antagonists, hexamethonium and mecamylamine. Focal application of ATP (100 muM) had no effect on membrane potential of the postsynaptic neurone or on the amplitude of spontaneous EPSPs. Taken together, these results suggest the absence of functional purinergic (P2) receptors on the postganglionic neurone in situ. In contrast, focally applied ATP (100 muM) reversibly inhibited nerve-evoked EPSPs. Similarly, bath application of the non-hydrolysable analogue of ATP, ATP gammaS, reversibly depressed EPSPs amplitude, The inhibitory effects of ATP and ATP gammaS on nerve-evoked transmitter release were antagonised by bath application of either PPADS or suramin, suggesting ATP activates a presynaptic P2 purinoceptor to inhibit acetylcholine release from preganglionic nerves in the submandibular ganglia. In acutely dissociated postganglionic neurotics from rat submandibular ganglia. focal application of ATP (100 LM) evoked an inward current and subsequent excitatory response and action potential firing, which was reversibly inhibited by PPADS (10 muM). The expression of P2X purinoceptors in wholemount and dissociated submandibular ganglion neurones was examined using polyclonal antibodies raised against the extracellular domain of six P2X purinoceptor subtypes (P2X(1-6)). In intact wholemount preparations, only the P2X(5) purinoceptor subtype was found to be expressed in the submandibular ganglion neurones and no P2X immunoreactivity was detected in the nerve fibres innervating the ganglion. Surprisingly, in dissociated submandibular ganglion neurones, high levels of P2X(2) and P2X(4) purinoceptors immunoreactivity were found on the cell surface. This increase in expression of P2X(2) and P2X(4) purinoceptors in dissociated submandibular neurones could explain the increased responsiveness of the neurotics to exogenous ATP. We conclude that disruption of ganglionic transmission in vivo by either nerve damage or synaptic blockade may up-regulate P2X expression or availability and alter neuronal excitability. (C) 2001 IBRO. Published by Elsevier Science Ltd. All rights reserved.

Precursor scavenging of resistive grain-boundary phase in 8 mol % ytterbia-stabilized zirconia

The grain-boundary conduction of 8 mol % ytterbia-stabilized zirconia (8YbSZ) was improved markedly by precursor scavenging via the two-stage sintering process. The most significant increase in the grain-boundary conductivity was found when the sample, whose conductivity was higher than that via Al2O3-derived scavenging, was heat-treated at 1250degreesC for greater than or equal to 20 h. The formation of a stable Si-containing phase such as ZrSiO4 during the first-stage heat-treatment was suggested as one probable scavenging route from the optimal heat-treatment temperature (HTT), long duration time (>20 h) at HTT, and the stability of the formed phase up to sintering temperatures (1500degrees C). (C) 2002 The Electrochemical Society.

Boundary value problems for systems of difference equations associated with systems of second-order ordinary differential equations

We study difference equations which arise as discrete approximations to two-point boundary value problems for systems of second-order ordinary differential equations. We formulate conditions which guarantee a priori bounds on first differences of solutions to the discretized problem. We establish existence results for solutions to the discretized boundary value problems subject to nonlinear boundary conditions. We apply our results to show that solutions to the discrete problem converge to solutions of the continuous problem in an aggregate sense. (C) 2002 Elsevier Science Ltd. 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.

Difference equations associated with fully nonlinear boundary value problems for second order ordinary differential equations

We study the continuous problem y"=f(x,y,y'), xc[0,1], 0=G((y(0),y(1)),(y'(0), y'(1))), and its discrete approximation (y(k+1)-2y(k)+y(k-1))/h(2) =f(t(k), y(k), v(k)), k = 1,..., n-1, 0 = G((y(0), y(n)), (v(1), v(n))), where f and G = (g(0), g(1)) are continuous and fully nonlinear, h = 1/n, v(k) = (y(k) - y(k-1))/h, for k =1,..., n, and t(k) = kh, for k = 0,...,n. We assume there exist strict lower and strict upper solutions and impose additional conditions on f and G which are known to yield a priori bounds on, and to guarantee the existence of solutions of the continuous problem. We show that the discrete approximation also has solutions which approximate solutions of the continuous problem and converge to the solution of the continuous problem when it is unique, as the grid size goes to 0. Homotopy methods can be used to compute the solution of the discrete approximation. Our results were motivated by those of Gaines.

Shear stress and sediment transport calculations for swash zone modelling

It is shown that the observed difference in sediment transporting efficiency by the swash uprush, compared with the downrush, could be mainly due to greater bed shear stress for a given velocity in the more abruptly accelerated uprush. The bed shear stress generated by an arbitrary free stream velocity time series is modelled in terms of usual wave boundary layer models plus a phase lead (phi(tau) of the bed shear stress compared with the free stream velocity at the peak frequency. With this approach, the total transport amounts in uprush and downrush can be modelled satisfactorily with the same sediment transport formula, without the need for different uprush and downrush coefficients. While the adaptation of sediment transport formulae from steady flow can thus lead to the right total amounts of sediment moved by this method, the timing of the instantaneous sediment transport rates are probably not accurately modelled due to the highly unsteady nature of the swash and the presence of pre-suspended sediment in the uprush. Nevertheless, the proposed method is a useful intermediate step before we have a complete understanding of sediment transport under very rapid accelerations and of the relative contribution of pre-suspended sediment to the onshore sediment transport in swash zones. (C) 2002 Published by Elsevier Science B.V.

A real symmetric Lanczos subspace implementation of quantum scattering using boundary inhomogeneities

We present an efficient and robust method for calculating state-to-state reaction probabilities utilising the Lanczos algorithm for a real symmetric Hamiltonian. The method recasts the time-independent Artificial Boundary Inhomogeneity technique recently introduced by Jang and Light (J. Chem. Phys. 102 (1995) 3262) into a tridiagonal (Lanczos) representation. The calculation proceeds at the cost of a single Lanczos propagation for each boundary inhomogeneity function and yields all state-to-state probabilities (elastic, inelastic and reactive) over an arbitrary energy range. The method is applied to the collinear H + H-2 reaction and the results demonstrate it is accurate and efficient in comparison with previous calculations. (C) 2002 Elsevier Science B.V. All rights reserved.