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.

Attic pottery in early Classical Thessaly. A case study from Achaia Phthiotis

This paper attributes a previously unnoticed Attic black-figured lekythos kept in the Archaeological Museum of Volos to the Pholos Group (ca. 470 BC) and discusses its findspot in the peripheral Thessalian district of Achaia Phthiotis. Beyond an art-historical appreciation of the hastily-drawn chariot scene on this lekythos and a discussion of stylistic parallels, which include a lekythos in Prague, it is argued that such lekythoi were socially important for their shape and small size that made them easily transportable. The assumed scarcity of Attic pottery in Thessaly can be questioned given that a considerable amount of Attic pottery from Thessalian locations is mentioned only in passing or remains unpublished in museum storage. Small late black-figured lekythoi predominate amongst Attic pottery shapes in Thessaly. The popularity of such lekythoi can become indicative of human mobility across the landscape and the consumption of imported (grave) goods by social groups other than elites.

Place making in Hellenistic Athens: beyond decline and foreign patronage

The study of cities is integral to the study of the Hellenistic Age, the period bounded by the deaths of two legendary rulers: Alexander in 323 BC and Kleopatra in 30 BC. Modern scholarship has followed in the footsteps of Johann-Gustav Droysen, who coined the term 'Hellenistic' in the nineteenth century and associated it with the diffusion of Greek culture through the founding of new cities in the East by Alexander and his successors. Hellenistic Athens, traditionally discussed under the rubric of its Classical legacy and/or in contrast with thriving cities, such as Pergamon, has been presented as a backwater exemplifying the demise of the 'polis'. My objective in this paper is to criticise these negative sentiments by exploring how the built environment of Hellenistic Athens could potentially become an indicator of city vitality.

Available potential energy and exergy in stratified fluids

Lorenz’s theory of available p otential energy (APE) remains the main framework for studying the atmospheric and oceanic energy cycles. Because the APE generation rate is the volume integral of a thermodynamic efficiency times the local diabatic heating/cooling rate, APE theory is often regarded as an extension of the theory of heat engines. Available energetics in classical thermodynamics, however, usually relies on the concept of exergy, and is usually measured relative to a reference state maximising entropy at constant energy, whereas APE’s reference state minimises p otential energy at constant entropy. This review seeks to shed light on the two concepts; it covers local formulations of available energetics, alternative views of the dynamics/thermodynamics coupling, APE theory and the second law, APE production/dissipation, extensions to binary fluids, mean/eddy decomp ositions, APE in incompressible fluids, APE and irreversible turbulent mixing, and the role of mechanical forcing on APE production.

On the entrainment assumption in Schatzmann's integral plume model

The behaviour of stationary, non-passive plumes can be simulated in a reasonably simple and accurate way by integral models. One of the key requirements of these models, but also one of their less well-founded aspects, is the entrainment assumption, which parameterizes turbulent mixing between the plume and the environment. The entrainment assumption developed by Schatzmann and adjusted to a set of experimental results requires four constants and an ad hoc hypothesis to eliminate undesirable terms. With this assumption, Schatzmann’s model exhibits numerical instability for certain cases of plumes with small velocity excesses, due to very fast radius growth. The purpose of this paper is to present an alternative entrainment assumption based on a first-order turbulence closure, which only requires two adjustable constants and seems to solve this problem. The asymptotic behaviour of the new formulation is studied and compared to previous ones. The validation tests presented by Schatzmann are repeated and it is found that the new formulation not only eliminates numerical instability but also predicts more plausible growth rates for jets in co-flowing streams.

'Si en i ot de teus qui i conterent plus lor honte que leur honour', Enadain et Gauvain, les chevaliers transformés en nains dans la Suite Vulgate du Merlin

Dans le récit de la transformation de Gauvain en nain, la honte est à la fois la conséquence d'une culpabilité mise en avant par un interlocuteur féminin, et l'instrument de la purgation de la faute commise par le héros. La gestion narrative de cet épisode à la fois dramatique et comique interroge autant les valeurs chrétiennes que les principes de conduite chevaleresques. La nature de la faute et son intériorisation personnelle importent moins que sa dimension sociale et collective.

Modeling electrocortical activity through improved local approximations of integral neural field equations

Neural field models of firing rate activity typically take the form of integral equations with space-dependent axonal delays. Under natural assumptions on the synaptic connectivity we show how one can derive an equivalent partial differential equation (PDE) model that properly treats the axonal delay terms of the integral formulation. Our analysis avoids the so-called long-wavelength approximation that has previously been used to formulate PDE models for neural activity in two spatial dimensions. Direct numerical simulations of this PDE model show instabilities of the homogeneous steady state that are in full agreement with a Turing instability analysis of the original integral model. We discuss the benefits of such a local model and its usefulness in modeling electrocortical activity. In particular, we are able to treat “patchy” connections, whereby a homogeneous and isotropic system is modulated in a spatially periodic fashion. In this case the emergence of a “lattice-directed” traveling wave predicted by a linear instability analysis is confirmed by the numerical simulation of an appropriate set of coupled PDEs.

A fast two-grid and finite section method for a class of integral equations on the real line with application to an acoustic scattering problem in the half-plane

We consider the numerical treatment of second kind integral equations on the real line of the form ∅(s) = ∫_(-∞)^(+∞)▒〖κ(s-t)z(t)ϕ(t)dt,s=R〗 (abbreviated ϕ= ψ+K_z ϕ) in which K ϵ L_1 (R), z ϵ L_∞ (R) and ψ ϵ BC(R), the space of bounded continuous functions on R, are assumed known and ϕ ϵ BC(R) is to be determined. We first derive sharp error estimates for the finite section approximation (reducing the range of integration to [-A, A]) via bounds on (1-K_z )^(-1)as an operator on spaces of weighted continuous functions. Numerical solution by a simple discrete collocation method on a uniform grid on R is then analysed: in the case when z is compactly supported this leads to a coefficient matrix which allows a rapid matrix-vector multiply via the FFT. To utilise this possibility we propose a modified two-grid iteration, a feature of which is that the coarse grid matrix is approximated by a banded matrix, and analyse convergence and computational cost. In cases where z is not compactly supported a combined finite section and two-grid algorithm can be applied and we extend the analysis to this case. As an application we consider acoustic scattering in the half-plane with a Robin or impedance boundary condition which we formulate as a boundary integral equation of the class studied. Our final result is that if z (related to the boundary impedance in the application) takes values in an appropriate compact subset Q of the complex plane, then the difference between ϕ(s)and its finite section approximation computed numerically using the iterative scheme proposed is ≤C_1 [kh log⁡〖(1⁄kh)+(1-Θ)^((-1)⁄2) (kA)^((-1)⁄2) 〗 ] in the interval [-ΘA,ΘA](Θ<1) for kh sufficiently small, where k is the wavenumber and h the grid spacing. Moreover this numerical approximation can be computed in ≤C_2 N log⁡N operations, where N = 2A/h is the number of degrees of freedom. The values of the constants C1 and C2 depend only on the set Q and not on the wavenumber k or the support of z.

A generalized collectively compact operator theory with an application to integral equations on unbounded domains

In this paper a generalization of collectively compact operator theory in Banach spaces is developed. A feature of the new theory is that the operators involved are no longer required to be compact in the norm topology. Instead it is required that the image of a bounded set under the operator family is sequentially compact in a weaker topology. As an application, the theory developed is used to establish solvability results for a class of systems of second kind integral equations on unbounded domains, this class including in particular systems of Wiener-Hopf integral equations with L1 convolutions kernels

Solvability and spectral properties of integral equations on the real line: I. Weighted spaces of continuous functions

We consider in this paper the solvability of linear integral equations on the real line, in operator form (λ−K)φ=ψ, where and K is an integral operator. We impose conditions on the kernel, k, of K which ensure that K is bounded as an operator on . Let Xa denote the weighted space as |s|→∞}. Our first result is that if, additionally, |k(s,t)|⩽κ(s−t), with and κ(s)=O(|s|−b) as |s|→∞, for some b>1, then the spectrum of K is the same on Xa as on X, for 01. As an example where kernels of this latter form occur we discuss a boundary integral equation formulation of an impedance boundary value problem for the Helmholtz equation in a half-plane.

On the stability and convergence of the finite section method for integral equation formulations of rough surface scattering

We consider the Dirichlet and Robin boundary value problems for the Helmholtz equation in a non-locally perturbed half-plane, modelling time harmonic acoustic scattering of an incident field by, respectively, sound-soft and impedance infinite rough surfaces.Recently proposed novel boundary integral equation formulations of these problems are discussed. It is usual in practical computations to truncate the infinite rough surface, solving a boundary integral equation on a finite section of the boundary, of length 2A, say. In the case of surfaces of small amplitude and slope we prove the stability and convergence as A→∞ of this approximation procedure. For surfaces of arbitrarily large amplitude and/or surface slope we prove stability and convergence of a modified finite section procedure in which the truncated boundary is ‘flattened’ in finite neighbourhoods of its two endpoints. Copyright © 2001 John Wiley & Sons, Ltd.