Difference equations in Banach spaces

Difference equations which discretely approximate boundary value problems for second-order ordinary differential equations are analysed. It is well known that the existence of solutions to the continuous problem does not necessarily imply existence of solutions to the discrete problem and, even if solutions to the discrete problem are guaranteed, they may be unrelated and inapplicable to the continuous problem. Analogues to theorems for the continuous problem regarding a priori bounds and existence of solutions are formulated for the discrete problem. Solutions to the discrete problem are shown to converge to solutions of the continuous problem in an aggregate sense. An example which arises in the study of the finite deflections of an elastic string under a transverse load is investigated. The earlier results are applied to show the existence of a solution; the sufficient estimates on the step size are presented. (C) 2003 Elsevier Science Ltd. All rights reserved.

Selecting methods to solve multi-well master equations

There are several competing methods commonly used to solve energy grained master equations describing gas-phase reactive systems. When it comes to selecting an appropriate method for any particular problem, there is little guidance in the literature. In this paper we directly compare several variants of spectral and numerical integration methods from the point of view of computer time required to calculate the solution and the range of temperature and pressure conditions under which the methods are successful. The test case used in the comparison is an important reaction in combustion chemistry and incorporates reversible and irreversible bimolecular reaction steps as well as isomerizations between multiple unimolecular species. While the numerical integration of the ODE with a stiff ODE integrator is not the fastest method overall, it is the fastest method applicable to all conditions.

Simultaneous nitrification, denitrification, and phosphorus removal in a lab-scale sequencing batch reactor

Simultaneous nitrification and denitrification (SND) via the nitrite pathway and anaerobic-anoxic-enhanced biological phosphorus removal (EBPR) are two processes that can significantly reduce the energy and COD demand for nitrogen and phosphorus removal. The combination of these two processes has the potential of achieving simultaneous nitrogen and phosphorus removal with a minimal requirement for COD. A lab-scale sequencing batch reactor (SBR) was operated in alternating anaerobic-aerobic mode with a low dissolved oxygen (DO) concentration (0.5 mg/L) during the aerobic period, and was demonstrated to accomplish nitrification, denitrification, and phosphorus removal. Under anaerobic conditions, COD was taken up and converted to poly-hydroxyalkanoates (PHAs), accompanied by phosphorus release. In the subsequent aerobic stage, PHA was oxidized and phosphorus was taken up to

An approach to verifying and debugging simulation models governed by ordinary differential equations: Part 1. Methodology for residual generation

For dynamic simulations to be credible, verification of the computer code must be an integral part of the modelling process. This two-part paper describes a novel approach to verification through program testing and debugging. In Part 1, a methodology is presented for detecting and isolating coding errors using back-to-back testing. Residuals are generated by comparing the output of two independent implementations, in response to identical inputs. The key feature of the methodology is that a specially modified observer is created using one of the implementations, so as to impose an error-dependent structure on these residuals. Each error can be associated with a fixed and known subspace, permitting errors to be isolated to specific equations in the code. It is shown that the geometric properties extend to multiple errors in either one of the two implementations. Copyright (C) 2003 John Wiley Sons, Ltd.

Simultaneous FT-NIR and ESR analyses to yield propagation rate coefficients for polymerization of methyl methacrylate based monomers

A hyphenated instrumental approach has been used to obtain reliable values for the propagation rate coefficients as a function of conversion for polymerizations of methyl methacrylate (MMA) and a mixture of MMA and ethyleneglycol dimethacrylate (EGDMA) with a 1:1 concentration of double bonds, from near the onset of the Trommsdorf region into the glass region. ESR spectroscopy was used to measure the radical concentration while FT-NIR fibre-optic spectroscopy was employed to measure instantaneously the double-bond concentration within the temperature-controlled cavity of the ESR instrument during polymerization. The advantage of this approach to the measurement of the rate coefficient is that it is equally applicable to branching and linear polymerizations. For the polymerization of methyl methacrylate, the values of the rate coefficient at the lowest conversions at which reliable values could be obtained were in agreement with recently reported values obtained by the PLP-SEC method. For the lowest conversions, the values obtained were 403 1 mol(-1) s(-1) at 306 K for MMA and 5201 mol(-1) s(-1) at 310 K for a 1:1 mixture of MMA and EGDMA. (C) 2003 Society of Chemical Industry.

Simultaneous observations of ionospheric quasiperiodic scintillations from short and long meridional baselines using VHF transmissions from transit satellites

For the first time it was possible to observe regular quasiperiodic scintillations (QPS) in VHF radio-satellite transmissions from orbiting satellites simultaneously at short (2.1 km) and long (121 km) meridional baselines in the vicinity of a typical mid-latitude station (Brisbane; 27.5degreesS and 152.9degreesE geog. and 35.6degrees invar.lat.), using three sites (St. Lucia-S, Taringa-T in Brisbane and Boreen Pt.-B, north of Brisbane). A few pronounced quasiperiodic (QP) events were recorded showing unambiguous regular structures at the sites which made it possible to deduce a time displacement of the regular fading minimum at S, T and B. The QP structure is highly dependent on the geometry of the ray-path from a satellite to the observer which is manifested as a change of a QP event from symmetrical to non-symmetrical for stations separated by 2.1 km, and to a radical change in the structure of the event over a distance of 121 km. It is suggested the short-duration intense QP events are due to a Fresnel diffraction (or a reflection mechanism) of radio-satellite signals by a single ionospheric irregularity in a form of an ellipsoid with a large ionization gradient along the major axis. The structure of a QP event depends on the angle of viewing of the irregular blob from a radio-satellite. In view of this it is suggested that the reported variety of the ionization formation, responsible for different types of QPS, is only apparent but not real. (C) 2003 Elsevier Science Ltd. All rights reserved.

A new wavelet-based adaptive method for solving population balance equations

A new wavelet-based adaptive framework for solving population balance equations (PBEs) is proposed in this work. The technique is general, powerful and efficient without the need for prior assumptions about the characteristics of the processes. Because there are steeply varying number densities across a size range, a new strategy is developed to select the optimal order of resolution and the collocation points based on an interpolating wavelet transform (IWT). The proposed technique has been tested for size-independent agglomeration, agglomeration with a linear summation kernel and agglomeration with a nonlinear kernel. In all cases, the predicted and analytical particle size distributions (PSDs) are in excellent agreement. Further work on the solution of the general population balance equations with nucleation, growth and agglomeration and the solution of steady-state population balance equations will be presented in this framework. (C) 2002 Elsevier Science B.V. All rights reserved.