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.

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.

Supersymmetric t-J Gaudin models and KZ equations

Supersymmetric t-J Gaudin models with open boundary conditions are investigated by means of the algebraic Bethe ansatz method. Off-shell Bethe ansatz equations of the boundary Gaudin systems are derived, and used to construct and solve the KZ equations associated with sl (2\1)((1)) superalgebra.

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.

Ab initio modelling of basal plane oxidation of graphenes and implications for modelling char combustion

Ab initio calculations have been performed to determine the energetics of oxygen atoms adsorbed onto graphene planes and the possible reaction path extracting carbon atorns in the form of carbon monoxide. Front the energetics it is confirmed that this reaction path will not significantly contribute to the gasification of well ordered carbonaceous chars. Modelling results which explore this limit Lire presented. (C) 2002 Elsevier Science Ltd, All rights reserved.

Population balance modelling of granulation with a physically based coalescence kernel

It was previously published by the authors that granules can either coalesce through Type I (when granules coalesce by viscous dissipation in the surface liquid layer before their surfaces touch) or Type II (when granules are slowed to a halt during rebound, after their surfaces have made contact) (AIChE J. 46 (3) (2000) 529). Based on this coalescence mechanism, a new coalescence kernel for population balance modelling of granule growth is presented. The kernel is constant such that only collisions satisfying the conditions for one of the two coalescence types are successful. One constant rate is assigned to each type of coalescence and zero is for the case of rebound. As the conditions for Types I and II coalescence are dependent on granule and binder properties, the coalescence kernel is thus physically based. Simulation results of a variety of binder and granule materials show good agreement with experimental data. (C) 2002 Elsevier Science Ltd. All rights reserved.

Negative ion electrospray mass spectra of caerulein peptides: an aid to structural determination

MS/MS data derived from the [M-H](-) ions of desulfated caerulein peptides provide (i) sequencing information from a combination of alpha, beta and gamma backbone cleavages, and (ii) identification of specific amino acid side chains by side-chain cleavages [e.g. Ser (-CH2O), Thr (-CH3CHO) and Asp (-H2O)] (fragmentations having no counterparts in positive ion spectra). In addition, delta and/or gamma backbone cleavage ions from Asp residues identify the position of these residues in the peptide. In contrast, neither delta nor gamma cleavage ions are observed from either the Gln2 residue nor from Phe residues. Full structural information can be obtained from a consideration of the positive and negative ion MS/MS data in concert. Copyright (C) 2002 John Wiley Sons, Ltd.

Modelling of microstructure formation and evolution during solidification

A comprehensive probabilistic model for simulating microstructure formation and evolution during solidification has been developed, based on coupling a Finite Differential Method (FDM) for macroscopic modelling of heat diffusion to a modified Cellular Automaton (mCA) for microscopic modelling of nucleation, growth of microstructures and solute diffusion. The mCA model is similar to Nastac's model for handling solute redistribution in the liquid and solid phases, curvature and growth anisotropy, but differs in the treatment of nucleation and growth. The aim is to improve understanding of the relationship between the solidification conditions and microstructure formation and evolution. A numerical algorithm used for FDM and mCA was developed. At each coarse scale, temperatures at FDM nodes were calculated while nucleation-growth simulation was done at a finer scale, with the temperature at the cell locations being interpolated from those at the coarser volumes. This model takes account of thermal, curvature and solute diffusion effects. Therefore, it can not only simulate microstructures of alloys both on the scale of grain size (macroscopic level) and the dendrite tip length (mesoscopic level), but also investigate nucleation mechanisms and growth kinetics of alloys solidified with various solute concentrations and solidification morphologies. The calculated results are compared with values of grain sizes and solidification morphologies of microstructures obtained from a set of casting experiments of Al-Si alloys in graphite crucibles.

Predictor-corrector methods of Runge-Kutta type for stochastic differential equations

In this paper we construct predictor-corrector (PC) methods based on the trivial predictor and stochastic implicit Runge-Kutta (RK) correctors for solving stochastic differential equations. Using the colored rooted tree theory and stochastic B-series, the order condition theorem is derived for constructing stochastic RK methods based on PC implementations. We also present detailed order conditions of the PC methods using stochastic implicit RK correctors with strong global order 1.0 and 1.5. A two-stage implicit RK method with strong global order 1.0 and a four-stage implicit RK method with strong global order 1.5 used as the correctors are constructed in this paper. The mean-square stability properties and numerical results of the PC methods based on these two implicit RK correctors are reported.

A variable stepsize implementation for stochastic differential equations

Stochastic differential equations (SDEs) arise from physical systems where the parameters describing the system can only be estimated or are subject to noise. Much work has been done recently on developing higher order Runge-Kutta methods for solving SDEs numerically. Fixed stepsize implementations of numerical methods have limitations when, for example, the SDE being solved is stiff as this forces the stepsize to be very small. This paper presents a completely general variable stepsize implementation of an embedded Runge Kutta pair for solving SDEs numerically; in this implementation, there is no restriction on the value used for the stepsize, and it is demonstrated that the integration remains on the correct Brownian path.

Structural changes in mixed Langmuir-Blodgett films upon nanoparticle formation

The formation of CdS nanoparticles by reacting mixed Langmuir-Blodgett films of arachidic acid and either octadecylamine or dimethyldioctadecylammonium nitrate on a cadmium-containing subphase with hydrogen sulfide gas has resulted in the identification of a number of structural changes, observed using grazing incidence X-ray diffraction. In the case of octadecylamine, the structure after reaction is a hexagonal close-packed array of surfactant-stabilized nanoclusters, with a lattice constant of a = 17.65 Angstrom. In both octadecylamine and dimethyldioctadecylammonium nitrate films, the presence of a unit cell tilted at 38degrees to the plane of the substrate was found. Despite these changes, the average nanoparticle size was unaffected by the addition of either second component to the film.

Structural and surface property changes of macadamia nut-shell char upon activation and high temperature treatment

Structural and surface property changes of macadamia nut-shell (MNS) char upon activation and high temperature treatment (HTT) were studied by high-resolution nitrogen adsorption, diffuse reflectance infra-red Fourier transform spectroscopy, X-ray photoelectron spectroscopy, and temperature-programmed desorption. It is found that activation of MNS char can be divided into the low extent activation which may involve the reactions of internal oxygen-containing groups and leads to the formation of comparatively uniform micropores, and the high extent activation which induces reactions between carbon and activating gas and produces a large amount of micropores. The surface functional groups (SFGs) basically increase with the increase of activation extent, but high extent activation preferentially increases the amount of -C-O and -C=O. HTT in air for a short tithe at a high temperature (1173 K) greatly increases the micropore volume and the amounts of SFGs. By appropriately choosing the activation and HTT conditions, it is possible to control both the textural structure and the type and amounts of SFG. (C) 2002 Published by Elsevier Science Ltd.

Effect of a control burn on lizards and their structural environment in a eucalypt open-forest

The composition of an open-forest lizard assemblage in eastern Australia was examined before and after a low-intensity controlled fire and concurrently compared with that in an adjoining unburnt area. The effect of fire on the available structural environment and the habitat used by two focal species, Carlia vivax and Lygisaurus foliorum, was also examined. Lizard species richness was unaffected by the controlled burn as was the abundance of most species. C. vivax was the only species to display a significant reduction in abundance after fire. While the low-intensity fire resulted in significant changes to the available structural environment, there were no compensatory shifts in the habitat preferences of either C. vivax or L. foliorum. The reduction in abundance of C. vivax was congruent with this species' avoidance of burnt areas. C. vivax displayed a non-random preference for ground cover and litter cover, which were reduced in burnt areas. Changes in the availability of preferred structural habitat features are likely to contribute to changes in the abundance of some lizard species. Therefore, even low-intensity disturbances can have an impact on lizard assemblages if critical habitat features are lost or become limiting.