The numerical solution of the Helmholtz equation for wave propagation problems in underwater acoustics

The Helmholtz Equation (-delta-K(2)n(2))u=0 with a variable index of refraction, n, and a suitable radiation condition at infinity serves as a model for a wide variety of wave propagation problems. A numerical algorithm was developed and a computer code implemented that can effectively solve this equation in the intermediate frequency range. The equation is discretized using the finite element method, thus allowing for the modeling of complicated geometrices (including interfaces) and complicated boundary conditions. A global radiation boundary condition is imposed at the far field boundary that is exact for an arbitrary number of propagating modes. The resulting large, non-selfadjoint system of linear equations with indefinite symmetric part is solved using the preconditioned conjugate gradient method applied to the normal equations. A new preconditioner is developed based on the multigrid method. This preconditioner is vectorizable and is extremely effective over a wide range of frequencies provided the number of grid levels is reduced for large frequencies. A heuristic argument is given that indicates the superior convergence properties of this preconditioner

On accuracy conditions for the numerical computation of waves

The Helmholtz equation (Delta + K(2)n(2))u = f with a variable index of refraction n, and a suitable radiation condition at infinity serves as a model for a wide variety of wave propagation problems. Such problems can be solved numerically by first truncating the given unbounded domain and imposing a suitable outgoing radiation condition on an artificial boundary and then solving the resulting problem on the bounded domain by direct discretization (for example, using a finite element method). In practical applications, the mesh size h and the wave number K, are not independent but are constrained by the accuracy of the desired computation. It will be shown that the number of points per wavelength, measured by (Kh)(-1), is not sufficient to determine the accuracy of a given discretization. For example, the quantity K(3)h(2) is shown to determine the accuracy in the L(2) norm for a second-order discretization method applied to several propagation models

Open Boundaries for the Nonlinear Schrodinger Equation

We present a new algorithm, the Time Dependent Phase Space Filter (TDPSF) which is used to solve time dependent Nonlinear Schrodinger Equations (NLS). The algorithm consists of solving the NLS on a box with periodic boundary conditions (by any algorithm). Periodically in time we decompose the solution into a family of coherent states. Coherent states which are outgoing are deleted, while those which are not are kept, reducing the problem of reflected (wrapped) waves. Numerical results are given, and rigorous error estimates are described. The TDPSF is compatible with spectral methods for solving the interior problem. The TDPSF also fails gracefully, in the sense that the algorithm notifies the user when the result is incorrect. We are aware of no other method with this capability.Comment: 21 pages, 4 figure

A Mesolithic settlement site at Howick, Northumberland: a preliminary report

Excavations at a coastal site at Howick during 2000 and 2002 have revealed evidence for a substantial Mesolithic settlement and a Bronze Age cist cemetery. Twenty one radiocarbon determinations of the earlier eighth millennium BP (Cal.) indicate that the Mesolithic site is one of the earliest known in northern Britain. An 8m core of sediment was recovered from stream deposits adjacent to the archaeological site which provides information on local environmental conditions. Howick offers a unique opportunity to understand aspects of hunter-gatherer colonisation and settlement during a period of rapid palaeogeographical change around the margins of the North Sea basin, at a time when it was being progressively inundated by the final stages of the postglacial marine transgression. The cist cemetery will add to the picture of Bronze Age occupation of the coastal strip and again reveals a correlation between the location of Bronze Age and Mesolithic sites which has been observed elsewhere in the region

Secure, performance-oriented data management for nanoCMOS electronics

The EPSRC pilot project Meeting the Design Challenges of nanoCMOS Electronics (nanoCMOS) is focused upon delivering a production level e-Infrastructure to meet the challenges facing the semiconductor industry in dealing with the next generation of ‘atomic-scale’ transistor devices. This scale means that previous assumptions on the uniformity of transistor devices in electronics circuit and systems design are no longer valid, and the industry as a whole must deal with variability throughout the design process. Infrastructures to tackle this problem must provide seamless access to very large HPC resources for computationally expensive simulation of statistic ensembles of microscopically varying physical devices, and manage the many hundreds of thousands of files and meta-data associated with these simulations. A key challenge in undertaking this is in protecting the intellectual property associated with the data, simulations and design process as a whole. In this paper we present the nanoCMOS infrastructure and outline an evaluation undertaken on the Storage Resource Broker (SRB) and the Andrew File System (AFS) considering in particular the extent that they meet the performance and security requirements of the nanoCMOS domain. We also describe how metadata management is supported and linked to simulations and results in a scalable and secure manner

Gaze cueing elicited by emotional faces is influenced by affective context

When we observe someone shift their gaze to a peripheral event or object, a corresponding shift in our own attention often follows. This social orienting response, joint attention, has been studied in the laboratory using the gaze cueing paradigm. Here, we investigate the combined influence of the emotional content displayed in two critical components of a joint attention episode: The facial expression of the cue face, and the affective nature of the to-be-localized target object. Hence, we presented participants with happy and disgusted faces as cueing stimuli, and neutral (Experiment 1), pleasant and unpleasant (Experiment 2) pictures as target stimuli. The findings demonstrate an effect of ‘emotional context’ confined to participants viewing pleasant pictures. Specifically, gaze cueing was boosted when the emotion of the gazing face (i.e., happy) matched that of the targets (pleasant). Demonstrating modulation by emotional context highlights the vital flexibility that a successful joint attention system requires in order to assist our navigation of the social world

Interaction of Sound from Supersonic Jets with Nearby Structures

A model of sound generated in an ideally expanded supersonic (Mach 2) jet is solved numerically. Two configurations are considered: (1) a free jet and (2) an installed jet with a nearby array of flexible aircraft type panels. In the later case the panels vibrate in response to loading by sound from the jet and the full coupling between the panels and the jet is considered, accounting for panel response and radiation. The long time behavior of the jet is considered. Results for near field and far field disturbance, the far field pressure and the vibration of and radiation from the panels are presented. Panel response crucially depends on the location of the panels. Panels located upstream of the Mach cone are subject to a low level, nearly continuous spectral excitation and consequently exhibit a low level, relatively continuous spectral response. In contrast, panels located within the Mach cone are subject to a significant loading due to the intense Mach wave radiation of sound and exhibit a large, relatively peaked spectral response centered around the peak frequency of sound radiation. The panels radiate in a similar fashion to the sound in the jet, in particular exhibiting a relatively peaked spectral response at approximately the Mach angle from the bounding wall

Improving the FEH statistical procedures for flood frequency estimation

This report recommends changes to the procedures contained in the Flood Estimation Handbook (FEH), which have been adopted as standard practice by the principal bodies engaged in flood frequency estimation in the UK and, in particular, by the Environment Agency. These procedures provide estimates of the flows that will occur in rivers on moderately rare occasions: flow values that have an exceedance probability in any given year of 50 per cent (a 2-year return period) to 1 per cent (a 100-year return period), or even more rare. In the majority of cases where such estimates are required, the locations affected will be ungauged and too far from established river gauging stations to provide data records that can be immediately transferred. The changes recommended arise, in part, because the HiFlows-UK project has led to the creation of a much-improved database of systematically recorded flood data. Not only are the data records now much longer than those used previously but the HiFlows-UK project put substantial effort into the quality control and assessment of the whole data-set. This means that the data available for analysis have been dramatically improved. Another influence on the renewed procedures has been feedback from users of the FEH, both informal and formal. Without substantially changing the overall framework of the methodology, most technical details of the method have been updated to improve the performance of the procedure. The updates include significant improvements to the theoretical statistical framework underlying the method. In addition, it has been possible to consider some new descriptors of catchment topography and local climate that have been proposed since the FEH study. In particular, a new descriptor that measures floodplain extent has been devised and is now included in the improved procedures. This report is largely a technical description of the studies that have led to the new recommendations. The folllowing are the key improvements. • A new regression model for estimating the median annual maximum flood (QMED) at ungauged catchments (Chapter 4). • An improved procedure for the use of donor catchments for estimation of QMED at ungauged catchments (Chapter 5). • An improved procedure for formation of pooling groups and estimation of pooled growth curves (Chapter 6). Flood estimates produced by the new procedures can be substantially different from those produced using the original FEH procedures. On taking the catchments whose data have been analysed as typical examples, and treating them as if they were ungauged, the ratios of the new estimates to the FEH estimates indicate the following changes. • The changes in QMED range from 0.55 to 2.01, with half being greater than 1.15 (25 per cent of the ratios are less than 1.00, and 25 per cent are greater than 1.24). • For floods with an annual probability of exceedance of 1 per cent (the 1 per cent flood), the changes range from 0.48 to 2.24, with half being greater than 1.14 (25 per cent of the ratios are less than 0.97 and 25 per cent are greater than 1.32). For both QMED and the 1 per cent flood, the new procedure produced lower estimates than the FEH in the East of England, whereas increases in both quantities were generally observed in West England, Wales and Scotland

Integrating security solutions to support nanoCMOS electronics research

The UK Engineering and Physical Sciences Research Council (EPSRC) funded Meeting the Design Challenges of nanoCMOS Electronics (nanoCMOS) is developing a research infrastructure for collaborative electronics research across multiple institutions in the UK with especially strong industrial and commercial involvement. Unlike other domains, the electronics industry is driven by the necessity of protecting the intellectual property of the data, designs and software associated with next generation electronics devices and therefore requires fine-grained security. Similarly, the project also demands seamless access to large scale high performance compute resources for atomic scale device simulations and the capability to manage the hundreds of thousands of files and the metadata associated with these simulations. Within this context, the project has explored a wide range of authentication and authorization infrastructures facilitating compute resource access and providing fine-grained security over numerous distributed file stores and files. We conclude that no single security solution meets the needs of the project. This paper describes the experiences of applying X.509-based certificates and public key infrastructures, VOMS, PERMIS, Kerberos and the Internet2 Shibboleth technologies for nanoCMOS security. We outline how we are integrating these solutions to provide a complete end-end security framework meeting the demands of the nanoCMOS electronics domain