Condition number estimates for combined potential integral operators in acoustics and their boundary element discretisation

We consider the classical coupled, combined-field integral equation formulations for time-harmonic acoustic scattering by a sound soft bounded obstacle. In recent work, we have proved lower and upper bounds on the $L^2$ condition numbers for these formulations, and also on the norms of the classical acoustic single- and double-layer potential operators. These bounds to some extent make explicit the dependence of condition numbers on the wave number $k$, the geometry of the scatterer, and the coupling parameter. For example, with the usual choice of coupling parameter they show that, while the condition number grows like $k^{1/3}$ as $k\to\infty$, when the scatterer is a circle or sphere, it can grow as fast as $k^{7/5}$ for a class of `trapping' obstacles. In this paper we prove further bounds, sharpening and extending our previous results. In particular we show that there exist trapping obstacles for which the condition numbers grow as fast as $\exp(\gamma k)$, for some $\gamma>0$, as $k\to\infty$ through some sequence. This result depends on exponential localisation bounds on Laplace eigenfunctions in an ellipse that we prove in the appendix. We also clarify the correct choice of coupling parameter in 2D for low $k$. In the second part of the paper we focus on the boundary element discretisation of these operators. We discuss the extent to which the bounds on the continuous operators are also satisfied by their discrete counterparts and, via numerical experiments, we provide supporting evidence for some of the theoretical results, both quantitative and asymptotic, indicating further which of the upper and lower bounds may be sharper.

Numerical-asymptotic boundary integral methods in high-frequency acoustic scattering

In this article we describe recent progress on the design, analysis and implementation of hybrid numerical-asymptotic boundary integral methods for boundary value problems for the Helmholtz equation that model time harmonic acoustic wave scattering in domains exterior to impenetrable obstacles. These hybrid methods combine conventional piecewise polynomial approximations with high-frequency asymptotics to build basis functions suitable for representing the oscillatory solutions. They have the potential to solve scattering problems accurately in a computation time that is (almost) independent of frequency and this has been realized for many model problems. The design and analysis of this class of methods requires new results on the analysis and numerical analysis of highly oscillatory boundary integral operators and on the high-frequency asymptotics of scattering problems. The implementation requires the development of appropriate quadrature rules for highly oscillatory integrals. This article contains a historical account of the development of this currently very active field, a detailed account of recent progress and, in addition, a number of original research results on the design, analysis and implementation of these methods.

Rossby number expansions, slaving principles, and balance dynamics

We consider the problem of constructing balance dynamics for rapidly rotating fluid systems. It is argued that the conventional Rossby number expansion—namely expanding all variables in a series in Rossby number—is secular for all but the simplest flows. In particular, the higher-order terms in the expansion grow exponentially on average, and for moderate values of the Rossby number the expansion is, at best, useful only for times of the order of the doubling times of the instabilities of the underlying quasi-geostrophic dynamics. Similar arguments apply in a wide class of problems involving a small parameter and sufficiently complex zeroth-order dynamics. A modified procedure is proposed which involves expanding only the fast modes of the system; this is equivalent to an asymptotic approximation of the slaving relation that relates the fast modes to the slow modes. The procedure is systematic and thus capable, at least in principle, of being carried to any order—unlike procedures based on truncations. We apply the procedure to construct higher-order balance approximations of the shallow-water equations. At the lowest order quasi-geostrophy emerges. At the next order the system incorporates gradient-wind balance, although the balance relations themselves involve only linear inversions and hence are easily applied. There is a large class of reduced systems associated with various choices for the slow variables, but the simplest ones appear to be those based on potential vorticity.

Equatorial Kelvin waves as revealed by EOS Microwave Limb Sounder observations and European Centre for Medium-Range Weather Forecasts analyses: Evidence for slow Kelvin waves of zonal wave number 3

[1] Temperature and ozone observations from the Microwave Limb Sounder (MLS) on the EOS Aura satellite are used to study equatorial wave activity in the autumn of 2005. In contrast to previous observations for the same season in other years, the temperature anomalies in the middle and lower tropical stratosphere are found to be characterized by a strong wave-like eastward progression with zonal wave number equal to 3. Extended empirical orthogonal function (EOF) analysis reveals that the wave 3 components detected in the temperature anomalies correspond to a slow Kelvin wave with a period of 8 days and a phase speed of 19 m/s. Fluctuations associated with this Kelvin wave mode are also apparent in ozone profiles. Moreover, as expected by linear theory, the ozone fluctuations observed in the lower stratosphere are in phase with the temperature perturbations, and peak around 20–30 hPa where the mean ozone mixing ratios have the steepest vertical gradient. A search for other Kelvin wave modes has also been made using both the MLS observations and the analyses from one experiment where MLS ozone profiles are assimilated into the European Centre for Medium-Range Weather Forecasts (ECMWF) data assimilation system via a 6-hourly 3D var scheme. Our results show that the characteristics of the wave activity detected in the ECMWF temperature and ozone analyses are in good agreement with MLS data.

Condition number estimates for combined potential boundary integral operators in acoustic scattering

We study the classical combined field integral equation formulations for time-harmonic acoustic scattering by a sound soft bounded obstacle, namely the indirect formulation due to Brakhage-Werner/Leis/Panic, and the direct formulation associated with the names of Burton and Miller. We obtain lower and upper bounds on the condition numbers for these formulations, emphasising dependence on the frequency, the geometry of the scatterer, and the coupling parameter. Of independent interest we also obtain upper and lower bounds on the norms of two oscillatory integral operators, namely the classical acoustic single- and double-layer potential operators.

Systolic random number generation for genetic algorithms

A parallel hardware random number generator for use with a VLSI genetic algorithm processing device is proposed. The design uses an systolic array of mixed congruential random number generators. The generators are constantly reseeded with the outputs of the proceeding generators to avoid significant biasing of the randomness of the array which would result in longer times for the algorithm to converge to a solution. 1 Introduction In recent years there has been a growing interest in developing hardware genetic algorithm devices [1, 2, 3]. A genetic algorithm (GA) is a stochastic search and optimization technique which attempts to capture the power of natural selection by evolving a population of candidate solutions by a process of selection and reproduction [4]. In keeping with the evolutionary analogy, the solutions are called chromosomes with each chromosome containing a number of genes. Chromosomes are commonly simple binary strings, the bits being the genes.

Estimating the hidden number of scrapie affected holdings in great Britain using a simple, truncated count model allowing for heterogeneity.

None of the current surveillance streams monitoring the presence of scrapie in Great Britain provide a comprehensive and unbiased estimate of the prevalence of the disease at the holding level. Previous work to estimate the under-ascertainment adjusted prevalence of scrapie in Great Britain applied multiple-list capture–recapture methods. The enforcement of new control measures on scrapie-affected holdings in 2004 has stopped the overlapping between surveillance sources and, hence, the application of multiple-list capture–recapture models. Alternative methods, still under the capture–recapture methodology, relying on repeated entries in one single list have been suggested in these situations. In this article, we apply one-list capture–recapture approaches to data held on the Scrapie Notifications Database to estimate the undetected population of scrapie-affected holdings with clinical disease in Great Britain for the years 2002, 2003, and 2004. For doing so, we develop a new diagnostic tool for indication of heterogeneity as well as a new understanding of the Zelterman and Chao’s lower bound estimators to account for potential unobserved heterogeneity. We demonstrate that the Zelterman estimator can be viewed as a maximum likelihood estimator for a special, locally truncated Poisson likelihood equivalent to a binomial likelihood. This understanding allows the extension of the Zelterman approach by means of logistic regression to include observed heterogeneity in the form of covariates—in case studied here, the holding size and country of origin. Our results confirm the presence of substantial unobserved heterogeneity supporting the application of our two estimators. The total scrapie-affected holding population in Great Britain is around 300 holdings per year. None of the covariates appear to inform the model significantly.

The effects of irrigation, nitrogen fertilizer and grain size on Hagberg falling number, specific weight and blackpoint of winter wheat

The effects of irrigation and nitrogen (N) fertilizer on Hagberg falling number (HFN), specific weight (SW) and blackpoint (BP) of winter wheat (Triticum aestivum L) were investigated. Mains water (+50 and +100 mm month(-1), containing 44 mg NO3- litre(-1) and 28 mg SO42- litre(-1)) was applied with trickle irrigation during winter (17 January-17 March), spring (21 March-20 May) or summer (24 May-23 July). In 1999/2000 these treatments were factorially combined with three N levels (0, 200, 400 kg N ha(-1)), applied to cv Hereward. In 2000/01 the 400 kg N ha(-1) treatment was replaced with cv Malacca given 200 kg N ha(-1). Irrigation increased grain yield, mostly by increasing grain numbers when applied in winter and spring, and by increasing mean grain weight when applied in summer. Nitrogen increased grain numbers and SW, and reduced BP in both years. Nitrogen increased HFN in 1999/2000 and reduced HFN in 2000/01. Effects of irrigation on HFN, SW and BP were smaller and inconsistent over year and nitrogen level. Irrigation interacted with N on mean grain weight: negatively for winter and spring irrigation, and positively for summer irrigation. Ten variables derived from digital image analysis of harvested grain were included with mean grain weight in a principal components analysis. The first principal component ('size') was negatively related to HFN (in two years) and BP (one year), and positively related to SW (two years). Treatment effects on dimensions of harvested grain could not explain all of the effects on HFN, BP and SW but the results were consistent with the hypothesis that water and nutrient availability, even when they were affected early in the season, could influence final grain quality if they influenced grain numbers and size. (C) 2004 Society of Chemical Industry

Nitrogen fertilizer and seed rate effects on Hagberg failing number of hybrid wheats and their parents are associated with alpha-amylase activity, grain cavity size and dormancy

Field experiments were carried out to assess the effects of nitrogen fertilization and seed rate on the Hagberg falling number (HFN) of commercial wheat hybrids and their parents. Applying nitrogen (200 kg N ha(-1)) increased HFN in two successive years. The HFN of the hybrid Hyno Esta was lower than either of its parents (Estica and Audace), particularly when nitrogen was not applied. Treatment effects on HFN were negatively associated with a-amylase activity. Phadebas grain blotting suggested two populations of grains with different types of a-amylase activity: Estica appeared to have a high proportion of grains with low levels of late maturity endosperm a-amylase activity (LMEA); Audace had a few grains showing high levels of germination amylase; and the hybrid, Hyno Esta, combined the sources from both parents to show heterosis for a-amylase activity. Applying nitrogen reduced both apparent LMEA and germination amylase. The effects on LMEA were associated with the size and disruption of the grain cavity, which was greater in Hyno Esta and Estica and in zero-nitrogen treatments. External grain morphology failed to explain much of the variation in LMEA and cavity size, but there was a close negative correlation between cavity size and protein content. Applying nitrogen increased post-harvest dormancy of the grain. Dormancy was greatest in Estica and least in Audace. It is proposed that effects of seed rate, genotype and nitrogen fertilizer on HFN are mediated through factors affecting the size and disruption of the grain cavity and therefore LMEA, and through factors affecting dormancy and therefore germination amylase. (c) 2004 Society of Chemical Industry.