On a class of preconditioners for solving the Helmholtz equation

In 1983, a preconditioner was proposed [J. Comput. Phys. 49 (1983) 443] based on the Laplace operator for solving the discrete Helmholtz equation efficiently with CGNR. The preconditioner is especially effective for low wavenumber cases where the linear system is slightly indefinite. Laird [Preconditioned iterative solution of the 2D Helmholtz equation, First Year's Report, St. Hugh's College, Oxford, 2001] proposed a preconditioner where an extra term is added to the Laplace operator. This term is similar to the zeroth order term in the Helmholtz equation but with reversed sign. In this paper, both approaches are further generalized to a new class of preconditioners, the so-called “shifted Laplace” preconditioners of the form Δφ−k²φ with . Numerical experiments for various wavenumbers indicate the effectiveness of the preconditioner. The preconditioner is evaluated in combination with GMRES, Bi-CGSTAB, and CGNR.     

Electron transfer induced dissociations of 2- and 4-alkyl cyclohexadienones

Several 2- and 4-alkylcyclohexadienones were prepared and shown to accept electrons to produce ketyl radical anions that dissociated rapidly at room temperature to release carbon-centered radicals and an aromatic phenoxide type anion. In the PET process with benzyl-substituted cyclohexadienones, initiated with triethylamine, the benzyl radicals dimerised or abstracted an H-atom from solvent. In electrochemical reductions, and in reductions with alkali metals in liquid ammonia, the benzyl radicals were further reduced to anions.     

The evolution of a truncated gaussian probability density through time—modelling animal liveweights after selection

The form of the probability density derived from the evolution in time of a previously truncated frequency distribution of animal Liveweights is of interest in animal husbandry. Truncated frequency distributions arise when the heavier animals are sold for slaughter and the lighter animals retained. The demands of modern quality assurance schemes require that, given information on animal growth, the farmer is able to estimate the number of animals that would meet the specifications at some time in the future after truncation. Assuming that animal growth can be described by a linear stochastic differential equation, we derive an explicit expression for the probability density of animal Liveweights at any time after the truncation of an initial Gaussian density. It is shown that this probability density converges rapidly to a Gaussian density, so that after about 20 days of typical growth rates for lambs, the resulting density is practically indistinguishable from Gaussian.     

Performance of a CW pumped Nd:YVO4 amplifier with kHz pulses

Until recently, simple and reliable high repetition-rate laser sources with nanosecond pulses much shorter than from conventional A-O Q-switch lasers were not available. However over the past 2 years we have developed such lasers based on proprietary fast E-O switching technology, which allows designs delivering 1 ns pulses and subnanosecond jitter for good synchronisation. The technology provides pulses with multi-kW peak power and repetition-rates to >100 kHz.Most recently, the performance of these short pulse lasers has been developed further by implementing oscillator/amplifier (master oscillator and power amplifier, MOPA) technology which increases the output to >1 W average power. Here we report on a simple model that has been used to predict the performance of the CW pumped Nd:YVO₄ amplifier used in the MOPA laser. The model is based on the well-known expressions for the saturated gain applying to laser pulses, but more usually applied to pulse-excited amplifiers. The model is shown to allow a good interpretation of the amplifier behaviour for kHz pulses and to be a useful tool for predicting the performance of the MOPA laser.     

Determination of flutter derivatives via a neural network approach

A neural-network-based method is offered to determine the flutter derivatives of section models under smooth and turbulent flows. The approach uses the observed dynamic responses to train an appropriate neural network. Subsequently, the modal parameters of the model for different mean velocities of wind flow are directly estimated using weight matrices in the neural network. The flutter derivatives can then be determined accurately. The validity of the present method is verified through numerical studies. Finally, the procedure is employed to process experimental data from an inverted-U-type section model, obtained from wind tunnel tests.     

Multi-Item Scheduling by Benders' Decomposition

Scheduling the production of several items requires the determination of production quantities in different periods in the presence of resource constraints. Several approximate and heuristic algorithms have been proposed to solve this problem. However, no method for finding an optimal solution has as yet been developed. It is shown that the problem may be solved advantageously using Benders' decomposition. The subproblem in Benders' decomposition is shown to be a transportation problem, and some strategies for solving the master problem are indicated. The paper concludes with a sample problem demonstrating the application of the method.     

Ionization potentials and geometrical structure of pure and oxygenated barium clusters

In this article, we present results obtained on structures of small bare and oxygenated barium clusters produced by the laser vaporization technique. Contrary to the very weak Ba—Ba binding (0.202 eV), the Ba—O binding is rather strong (5.8 eV). As a consequence, extremely small amounts of oxygen in the cluster source are sufficient to induce the nucleation of oxygenated barium clusters. Surprisingly, at least for substoechiometric systems, Ba_nO_m clusters exhibit the same icosahedral structure as bare clusters do. Mass spectra seem to indicate that replacing a barium atom with a strongly polar BaO molecule does not affect the geometrical structure of the cluster. This aspect is more easily understandable if we assume that the oxygen atom takes place inside the structure of the Ba_n cluster without significant distortions. Ab-initio calculations have been undertaken in order to give a satisfactory account of this hypothesis.     

Almost-Sure Central Limit Theorem for Directed Polymers and Random Corrections

We consider a general model of directed polymers on the lattice , weakly coupled to a random environment. We prove that the central limit theorem holds almost surely for the discrete time random walk X _T associated to the polymer. Moreover we show that the random corrections to the cumulants of X _T are finite, starting from some dimension depending on the index of the cumulants, and that there are corresponding random corrections of order , , in the asymptotic expansion of the expectations of smooth functions of X _T . Full proofs are carried out for the first two cumulants. We finally prove a kind of local theorem showing that the ratio of the probabilities of the events to the corresponding probabilities with no randomness, in the region of “moderate” deviations from the average drift bT, are, for almost all choices of the environment, uniformly close, as , to a functional of the environment “as seen from (T,y)$”. Received: 14 October 1996 / Accepted: 28 March 1997