Modern scientific methods and their potential in wastewater science and technology

Application of novel analytical and investigative methods such as fluorescence in situ hybridization, confocal laser scanning microscopy (CLSM), microelectrodes and advanced numerical simulation has led to new insights into micro-and macroscopic processes in bioreactors. However, the question is still open whether or not these new findings and the subsequent gain of knowledge are of significant practical relevance and if so, where and how. To find suitable answers it is necessary for engineers to know what can be expected by applying these modern analytical tools. Similarly, scientists could benefit significantly from an intensive dialogue with engineers in order to find out about practical problems and conditions existing in wastewater treatment systems. In this paper, an attempt is made to help bridge the gap between science and engineering in biological wastewater treatment. We provide an overview of recently developed methods in microbiology and in mathematical modeling and numerical simulation. A questionnaire is presented which may help generate a platform from which further technical and scientific developments can be accomplished. Both the paper and the questionnaire are aimed at encouraging scientists and engineers to enter into an intensive, mutually beneficial dialogue. (C) 2002 Elsevier Science Ltd. All rights reserved.

Mathematical methods for spatially cohesive reserve design

The problem of designing spatially cohesive nature reserve systems that meet biodiversity objectives is formulated as a nonlinear integer programming problem. The multiobjective function minimises a combination of boundary length, area and failed representation of the biological attributes we are trying to conserve. The task is to reserve a subset of sites that best meet this objective. We use data on the distribution of habitats in the Northern Territory, Australia, to show how simulated annealing and a greedy heuristic algorithm can be used to generate good solutions to such large reserve design problems, and to compare the effectiveness of these methods.

A new technique for the quantification of peripheral edema with application in both unilateral and bilateral cases

Current noninvasive techniques for the routine and frequent quantification of peripheral lymphedema in patients are total limb volume measurement (by water immersion or by circumferential measurements) and bioelectrical impedance analysis (BIA). However both of these techniques require standardizing the measurement using a contralateral measurement from the unaffected limb, Hence these techniques are essentially restricted to unilateral lymphedema. This paper describes the results from a preliminary study to investigate an alternative approach to the analysis of the data from multiple frequency BIA to produce an index of lymphedema without the need for normalization to another body segment. Twenty patients receiving surgical treatment for breast cancer were monitored prior to surgery and again after diagnosis with unilateral lymphedema. The data recorded were total limb volume, by circumferential measurements; and BIA measurements of both limbs. From these measurements total limb volumes and extracellular fluid volumes were calculated and expressed as ratios of the affected limb to that of the unaffected limb. An index of the ratio of the extracellular fluid volume to the intracellular fluid volume was determined. This ECW/ICW index was calculated for both the affected and unaffected limbs at both measurement times. Results confirmed that the established techniques of total limb volume and extracellular fluid volume normalized to the unaffected contralateral limb were accurate in the detection of lymphedema (p < 10(-6)). Comparison of the ECW/ICW index from the affected limb after diagnosis with that from the pre-surgery measurement revealed a significant (p< 10(-6)) and considerable (75%) increase. The results of this pilot study suggest that by using multiple frequency bioelectrical impedance analysis, an index of the ECW/ICW ratio can be obtained and this index appears to have an equal, or better, sensitivity, than the other techniques in detecting lymphedema. More importantly, this index does not require normalization to another body segment and can be used to detect all types of peripheral edema including both unilateral and bilateral lymphedema.

A refinement calculus for logic programs

Existing refinement calculi provide frameworks for the stepwise development of imperative programs from specifications. This paper presents a refinement calculus for deriving logic programs. The calculus contains a wide-spectrum logic programming language, including executable constructs such as sequential conjunction, disjunction, and existential quantification, as well as specification constructs such as general predicates, assumptions and universal quantification. A declarative semantics is defined for this wide-spectrum language based on executions. Executions are partial functions from states to states, where a state is represented as a set of bindings. The semantics is used to define the meaning of programs and specifications, including parameters and recursion. To complete the calculus, a notion of correctness-preserving refinement over programs in the wide-spectrum language is defined and refinement laws for developing programs are introduced. The refinement calculus is illustrated using example derivations and prototype tool support is discussed.

Control of chaotic instability in a dual-spin spacecraft with dissipation using energy methods

Control of chaotic instability in a rotating multibody system in the form of a dual-spin spacecraft with an axial nutational damper is achieved using an algorithm derived using energy methods. The control method is implemented on two realistic spacecraft parameter configurations which have been found to exhibit chaotic instability when a sinusoidally varying torque is applied to the spacecraft for a range of forcing amplitudes and frequencies. Such a torque, in practice, may arise under malfunction of the control system or from an unbalanced rotor. Chaotic instabilities arising from these torques could introduce uncertainties and irregularities into a spacecraft's attitude and consequently impair pointing accuracy. The control method is formulated from nutational stability results derived using an energy sink approximation for a dual-spin spacecraft with an asymmetric platform and axisymmetric rotor. The effectiveness of the control method is shown numerically and the results are studied by means of time history, phase space, Poincare map, Lyapunov characteristic exponents and Bifurcation diagrams.

Determinants of tissue Doppler measures of regional diastolic function during dobutamine stress echocardiography

Background Diastolic dysfunction induced by ischemia may alter transmitral blood flow, but this reflects global ventricular function, and pseudonormalization may occur with increased preload. Tissue Doppler may assess regional diastolic function and is relatively load-independent, but limited data exist regarding its application to stress testing. We sought to examine the stress response of regional diastolic parameters to dobutomine echocardiography (DbE). Methods Sixty-three patients underwent study with DbE: 20 with low probability of coronary artery disease (CAD) and 43 with CAD who underwent angiography. A standard DbE protocol was used, and segments were categorized as ischemic, scar, or normal. Color tissue Doppler was acquired at baseline and peak stress, and waveforms in the basal and mid segments were used to measure early filling (Em), late filling (Am), and E deceleration time. Significant CAD was defined by stenoses >50% vessel diameter. Results Diastolic parameters had limited feasibility because of merging of Em and Am waves at high heart rates and limited reproducibility. Nonetheless, compared with normal segments, segments subtended with significant stenoses showed a lower Em velocity at rest (6.2 +/- 2.6 cm/s vs 4.8 +/- 2.2 cm/s, P < .0001) and peak (7.5 +/- 4.2 cm/s vs 5.1 +/- 3.6 cm/s, P < .0001), Abnormal segments also showed a shorter E deceleration time (51 +/- 27 ms vs 41 +/- 27 ms, P = .0001) at base and peak. No changes were documented in Am. The same pattern was seen with segments identified as ischemic with wall motion score. However, in the absence of ischemia, segments of patients with left ventricular hypertrophy showed a lower Em velocity, with blunted Em responses to stress. Conclusion Regional diastolic function is sensitive to ischemia. However, a number of practical limitations limit the applicability of diastolic parameters for the quantification of stress echocardiography.

Development of a fully quantitative approach to the interpretation of stress echocardiography using radial and longitudinal myocardial velocities

Background: Tissue Doppler may be used to quantify regional left ventricular function but is limited by segmental variation of longitudinal velocity from base to apex and free to septal walls. We sought to overcome this by developing a composite of longitudinal and radial velocities. Methods and Results. We examined 82 unselected patients undergoing a standard dobutamine echocardiogram. Longitudinal velocity was obtained in the basal and mid segments of each wall using tissue Doppler in the apical views. Radial velocities were derived in the same segments using an automated border detection system and centerline method with regional chords grouped according to segment location and temporally averaged. In 25 patients at low probability of coronary disease, the pattern of regional variation in longitudinal velocity (higher in the septum) was the opposite of radial velocity (higher in the free wall) and the combination was homogenous. In 57 patients undergoing angiography, velocity in abnormal segments was less than normal segments using longitudinal (6.0 +/- 3.6 vs 9.0 +/- 2.2 cm/s, P = .01) and radial velocity (6.0 +/- 4.0 vs 8.0 +/- 3.9 cm/s, P = .02). However, the composite velocity permitted better separation of abnormal and normal segments (13.3 +/- 5.6 vs 17.5 +/- 4.2 cm/s, P = .001). There was no significant difference between the accuracy of this quantitative approach and expert visual wall motion analysis (81% vs 84%, P = .56). Conclusion: Regional variation of uni-dimensional myocardial velocities necessitates site-specific normal ranges, probably because of different fiber directions. Combined analysis of longitudinal and radial velocities allows the derivation of a composite velocity, which is homogenous in all segments and may allow better separation of normal and abnormal myocardium.

Improvements to peroxide oxidation methods for analysing sulfur in acid sulfate soils

Improvements to peroxide oxidation methods for analysing acid sulfate soils (ASS) are introduced. The soil solution ratio has been increased to 1 : 40, titrations are performed in suspension, and the duration of the peroxide digest stage is substantially shortened. For 9 acid sulfate soils, the peroxide oxidisable sulfur value obtained using the improved method was compared with the reduced inorganic sulfur result obtained using the chromium reducible sulfur method. Their regression was highly significant, the slope of the regression line was not significantly different (P = 0.05) from unity, and the intercept not significantly different from zero. A complete sulfur budget for the improved method showed there was no loss of sulfur as has been reported for earlier peroxide oxidation techniques. When soils were very finely ground, efficient oxidation of sulfides was achieved, despite the milder digestion conditions. Highly sulfidic and organic soils were shown to be the most difficult to analyse using either the improved method or the chromium method. No single analytical method can be universally applied to all ASS, rather a suite of methods is necessary for a thorough understanding of many ASS. The improved peroxide method, in combination with the chromium method and the 4 M HCl extraction, form a sound platform for informed decision making on the management of acid sulfate soils.

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.

An atomistic simulation of solid state sintering using Monte Carlo methods

This paper presents results on the simulation of the solid state sintering of copper wires using Monte Carlo techniques based on elements of lattice theory and cellular automata. The initial structure is superimposed onto a triangular, two-dimensional lattice, where each lattice site corresponds to either an atom or vacancy. The number of vacancies varies with the simulation temperature, while a cluster of vacancies is a pore. To simulate sintering, lattice sites are picked at random and reoriented in terms of an atomistic model governing mass transport. The probability that an atom has sufficient energy to jump to a vacant lattice site is related to the jump frequency, and hence the diffusion coefficient, while the probability that an atomic jump will be accepted is related to the change in energy of the system as a result of the jump, as determined by the change in the number of nearest neighbours. The jump frequency is also used to relate model time, measured in Monte Carlo Steps, to the actual sintering time. The model incorporates bulk, grain boundary and surface diffusion terms and includes vacancy annihilation on the grain boundaries. The predictions of the model were found to be consistent with experimental data, both in terms of the microstructural evolution and in terms of the sintering time. (C) 2002 Elsevier Science B.V. All rights reserved.