Mesh selection for stiff two-point boundary value problems

This paper is concerned with the mesh selection algorithm of COLSYS, a well known collocation code for solving systems of boundary value problems. COLSYS was originally designed to solve non-stiff and mildly stiff problems only. In this paper we show that its performance for solving extremely stiff problems can be considerably improved by modifying its error estimation and mesh selection algorithms. Numerical examples indicate the superiority of the modified algorithm.Dedicated to John Butcher on the occasion of his sixtieth birthday     

Timetabling County Cricket Fixtures Using a Form of Tabu Search

This paper concerns a computer system which has been devised to timetable the county cricket matches in England. The context of the work is discussed and the many constraints and objectives described. Details of the solution method are given; it involves tabu search with a form of diversification dependent on the solution subcosts as well as the overall cost. Experiments are described which suggest that this type of diversification may have wide applicability to large, complex, multi-objective combinatorial problems.     

Diffractive proteon fragmentation in 16, 32, and 110 GeV/cK − p interactions and the diffractive production of heavy quarks(s\bar s,c\bar c)

The inclusive proton diffraction dissociation cross sections in 16, 32, and 110 GeV/cK ^? p interactions are determined from the spike nearx=1 in the inclusive negative particle spectra and are compared to those obtained inK?p interactions using other selection methods at various energies. The same procedure is applied to events containing aV ⁰ in order to obtain the cross section for diffractive \(s\bar s\) production. While the total cross section for proton diffraction is found to be approximately constant in the energy range studied here, proton diffraction yielding an \(s\bar s - pair\) is found to increase significantly. In particular it is almost constant at 85 μb forΛ ⁰ and Σ production but for \(NK\bar K\) it rises from zero at 16 GeV.c to about 200 μb at 110 GeV/c. From the result for \(s\bar s\) diffractive production an estimate for the \(c\bar c\) diffractive production cross section of approximately 1–10 μb at 110 GeV/c is obtained.     

A LES-CMC formulation for premixed flames including differential diffusion

A finite volume large eddy simulation–conditional moment closure (LES-CMC) numerical framework for premixed combustion developed in a previous studyhas been extended to account for differential diffusion. The non-unity Lewis number CMC transport equation has an additional convective term in sample space proportional to the conditional diffusion of the progress variable, that in turn accounts for diffusion normal to the flame front and curvature-induced effects. Planar laminar simulations are first performed using a spatially homogeneous non-unity Lewis number CMC formulation and validated against physical-space fully resolved reference solutions. The same CMC formulation is subsequently used to numerically investigate the effects of curvature for laminar flames having different effective Lewis numbers: a lean methane–air flame with Le_eff = 0.99 and a lean hydrogen–air flame with Le_eff = 0.33. Results suggest that curvature does not affect the conditional heat release if the effective Lewis number tends to unity, so that curvature-induced transport may be neglected. Finally, the effect of turbulence on the flame structure is qualitatively analysed using LES-CMC simulations with and without differential diffusion for a turbulent premixed bluff body methane–air flame exhibiting local extinction behaviour. Overall, both the unity and the non-unity computations predict the characteristic M-shaped flame observed experimentally, although some minor differences are identified. The findings suggest that for the high Karlovitz number (from 1 to 10) flame considered, turbulent mixing within the flame weakens the differential transport contribution by reducing the conditional scalar dissipation rate and accordingly the conditional diffusion of the progress variable.     

Detailed large-signal dynamic modelling of DFB laser structures and comparison with experiment

This paper describes the first general large-signal dynamic multiple-mode laser model that incorporates all the main mechanisms known to influence the dynamic behaviour of DFB laser structures with the exception of thermal effects: longitudinal mode spatial hole burning, carrier transport effects, nonlinear gain, and laser and submount parasitics. The time evolution of the output power and wavelength of all modes is predicted, and full spectra can be plotted as a function of time. The model has been extended to include an approximation to the effects of propagation down dispersive fibre, thereby allowing the simulation of filtered received eye diagrams. Detailed comparison of the model with the experimental performance of 2×/8 DFB lasers has shown good agreement, allowing the performance to be optimized, particularly with respect to longitudinal hole burning and carrier transport. The model is also applied to gain-switched operation of 2×/8 DFB structures, fast pulsing of three-section /4 DFB lasers, and the dynamic behaviour of complex coupling coefficient DFB laser structures.     

Modelling hydrogen in silicate materials

The expectation that the hydroxyl analogue state of positive muons implanted in quartz should be stabilized by a suitable charge compensation centre is supported by preliminary data. The diamagnetic fraction correlates with aluminium content over the narrow range investigated. This revised version was published online in August 2006 with corrections to the Cover Date.     

The Search for Appropriate O.R.: A Review of Operational Research in Developing Countries

This paper is a review of the literature relating to the application of O.R. in developing countries. The paper attempts to identify the special factors which affect decision problems in developing countries. Some desirable characteristics for O.R. in this context are outlined. The paper then summarises the results of a review of published studies, and some criticisms are made of current O.R. work.     

Uncertainty about the correlation among temporal envelopes in two comodulation tasks

The threshold of a 1250-Hz tonal signal was measured in the presence of five noise bands (each 50 Hz wide, centered at 850, 1050, 1250, 1450, and 1650 Hz) under five conditions of uncertainty about the waveform type ("correlated" or "uncorrelated"), and/or the specific waveform sample to be presented. The waveform type was correlated when the temporal envelopes of all of the noise bands were the same, and was uncorrelated when the temporal envelope of the band centered on the signal differed from the common envelope of the other bands. At the low-uncertainty end of the continuum of conditions, the same waveform type was presented throughout an entire block of trials, and, in addition, the same waveform sample was presented on the two observation intervals of a single trial (but changed across trials). At the high-uncertainty end of the continuum, both the waveform type and the waveform sample were chosen at random for every observation interval. Threshold estimates obtained from trials in which both observation intervals contained the same waveform type were not affected by uncertainty about the waveform sample within a trial, nor by uncertainty about the waveform type introduced across trials. Thus the comodulation masking release, or CMR (the difference in the thresholds obtained with the uncorrelated and correlated waveforms), calculated from these types of trials was robust across all of the uncertainty conditions. However, on those trials in which one correlated interval and one uncorrelated interval were paired, threshold estimates were influenced by a bias for listeners to choose the uncorrelated interval as the signal interval, whether or not it actually contained the signal. This bias reveals the importance of recognizing the contribution of the nonsignal interval in experiments involving masker uncertainty. Parallel results were obtained using the comodulation detection difference (CDD) task. In some conditions, marked individual differences were observed.