Sparse space-time boundary element methods for the heat equation

The goal of this work is the efficient solution of the heat equation with Dirichlet or Neumann boundary conditions using the Boundary Elements Method (BEM). Efficiently solving the heat equation is useful, as it is a simple model problem for other types of parabolic problems. In complicated spatial domains as often found in engineering, BEM can be beneficial since only the boundary of the domain has to be discretised. This makes BEM easier than domain methods such as finite elements and finite differences, conventionally combined with time-stepping schemes to solve this problem. The contribution of this work is to further decrease the complexity of solving the heat equation, leading both to speed gains (in CPU time) as well as requiring smaller amounts of memory to solve the same problem. To do this we will combine the complexity gains of boundary reduction by integral equation formulations with a discretisation using wavelet bases. This reduces the total work to O(h

A Galerkin boundary element method for high frequency scattering by convex polygons

In this paper we consider the problem of time-harmonic acoustic scattering in two dimensions by convex polygons. Standard boundary or finite element methods for acoustic scattering problems have a computational cost that grows at least linearly as a function of the frequency of the incident wave. Here we present a novel Galerkin boundary element method, which uses an approximation space consisting of the products of plane waves with piecewise polynomials supported on a graded mesh, with smaller elements closer to the corners of the polygon. We prove that the best approximation from the approximation space requires a number of degrees of freedom to achieve a prescribed level of accuracy that grows only logarithmically as a function of the frequency. Numerical results demonstrate the same logarithmic dependence on the frequency for the Galerkin method solution. Our boundary element method is a discretization of a well-known second kind combined-layer-potential integral equation. We provide a proof that this equation and its adjoint are well-posed and equivalent to the boundary value problem in a Sobolev space setting for general Lipschitz domains.

Boundary value problems for third-order linear PDEs in time-dependent domains

We present the extension of a methodology to solve moving boundary value problems from the second-order case to the case of the third-order linear evolution PDE qt + qxxx = 0. This extension is the crucial step needed to generalize this methodology to PDEs of arbitrary order. The methodology is based on the derivation of inversion formulae for a class of integral transforms that generalize the Fourier transform and on the analysis of the global relation associated with the PDE. The study of this relation and its inversion using the appropriate generalized transform are the main elements of the proof of our results.

A wavenumber independent boundary element method for an acoustic scattering problem

In this paper we consider the impedance boundary value problem for the Helmholtz equation in a half-plane with piecewise constant boundary data, a problem which models, for example, outdoor sound propagation over inhomogeneous. at terrain. To achieve good approximation at high frequencies with a relatively low number of degrees of freedom, we propose a novel Galerkin boundary element method, using a graded mesh with smaller elements adjacent to discontinuities in impedance and a special set of basis functions so that, on each element, the approximation space contains polynomials ( of degree.) multiplied by traces of plane waves on the boundary. We prove stability and convergence and show that the error in computing the total acoustic field is O( N-(v+1) log(1/2) N), where the number of degrees of freedom is proportional to N logN. This error estimate is independent of the wavenumber, and thus the number of degrees of freedom required to achieve a prescribed level of accuracy does not increase as the wavenumber tends to infinity.

A high-wavenumber boundary-element method for an acoustic scattering problem

In this paper we show stability and convergence for a novel Galerkin boundary element method approach to the impedance boundary value problem for the Helmholtz equation in a half-plane with piecewise constant boundary data. This problem models, for example, outdoor sound propagation over inhomogeneous flat terrain. To achieve a good approximation with a relatively low number of degrees of freedom we employ a graded mesh with smaller elements adjacent to discontinuities in impedance, and a special set of basis functions for the Galerkin method so that, on each element, the approximation space consists of polynomials (of degree $\nu$) multiplied by traces of plane waves on the boundary. In the case where the impedance is constant outside an interval $[a,b]$, which only requires the discretization of $[a,b]$, we show theoretically and experimentally that the $L_2$ error in computing the acoustic field on $[a,b]$ is ${\cal O}(\log^{\nu+3/2}|k(b-a)| M^{-(\nu+1)})$, where $M$ is the number of degrees of freedom and $k$ is the wavenumber. This indicates that the proposed method is especially commendable for large intervals or a high wavenumber. In a final section we sketch how the same methodology extends to more general scattering problems.

Scale-invariant moving finite elements for nonlinear partial differential equations in two dimensions

A scale-invariant moving finite element method is proposed for the adaptive solution of nonlinear partial differential equations. The mesh movement is based on a finite element discretisation of a scale-invariant conservation principle incorporating a monitor function, while the time discretisation of the resulting system of ordinary differential equations is carried out using a scale-invariant time-stepping which yields uniform local accuracy in time. The accuracy and reliability of the algorithm are successfully tested against exact self-similar solutions where available, and otherwise against a state-of-the-art h-refinement scheme for solutions of a two-dimensional porous medium equation problem with a moving boundary. The monitor functions used are the dependent variable and a monitor related to the surface area of the solution manifold. (c) 2005 IMACS. Published by Elsevier B.V. All rights reserved.

A high frequency boundary element method for scattering by convex polygons with impedance boundary conditions

We consider scattering of a time harmonic incident plane wave by a convex polygon with piecewise constant impedance boundary conditions. Standard finite or boundary element methods require the number of degrees of freedom to grow at least linearly with respect to the frequency of the incident wave in order to maintain accuracy. Extending earlier work by Chandler-Wilde and Langdon for the sound soft problem, we propose a novel Galerkin boundary element method, with the approximation space consisting of the products of plane waves with piecewise polynomials supported on a graded mesh with smaller elements closer to the corners of the polygon. Theoretical analysis and numerical results suggest that the number of degrees of freedom required to achieve a prescribed level of accuracy grows only logarithmically with respect to the frequency of the incident wave.

Mechanisms for the land/sea warming contrast exhibited by simulations of climate change

The land/sea warming contrast is a phenomenon of both equilibrium and transient simulations of climate change: large areas of the land surface at most latitudes undergo temperature changes whose amplitude is more than those of the surrounding oceans. Using idealised GCM experiments with perturbed SSTs, we show that the land/sea contrast in equilibrium simulations is associated with local feedbacks and the hydrological cycle over land, rather than with externally imposed radiative forcing. This mechanism also explains a large component of the land/sea contrast in transient simulations as well. We propose a conceptual model with three elements: (1) there is a spatially variable level in the lower troposphere at which temperature change is the same over land and sea; (2) the dependence of lapse rate on moisture and temperature causes different changes in lapse rate upon warming over land and sea, and hence a surface land/sea temperature contrast; (3) moisture convergence over land predominantly takes place at levels significantly colder than the surface; wherever moisture supply over land is limited, the increase of evaporation over land upon warming is limited, reducing the relative humidity in the boundary layer over land, and hence also enhancing the land/sea contrast. The non-linearity of the Clausius–Clapeyron relationship of saturation specific humidity to temperature is critical in (2) and (3). We examine the sensitivity of the land/sea contrast to model representations of different physical processes using a large ensemble of climate model integrations with perturbed parameters, and find that it is most sensitive to representation of large-scale cloud and stomatal closure. We discuss our results in the context of high-resolution and Earth-system modelling of climate change.

A note on boundary layer friction in baroclinic cyclones

The interaction between extratropical cyclones and the underlying boundary layer has been a topic of recent discussion in papers by Adamson et al (2006) and Beare (2007). Their results emphasise different mechanisms through which the boundary layer dynamics may modify the growth of a baroclinic cyclone. By using different sea-surface temperature distributions and comparing the low-level winds, the differences are exposed and both of the proposed mechanisms appear to be acting within a single simulation.

Numerical simulation of baroclinic waves with a parameterized boundary layer

A dry three-dimensional baroclinic life cycle model is used to investigate the role of turbulent fluxes of heat and momentum within the boundary layer on mid-latitude cyclones. Simulations are performed of life cycles for two basic states, both with and without turbulent fluxes. The different basic states produce cyclones with contrasting frontal and mesoscale-flow structures. The analysis focuses on the generation of potential-vorticity (PV) in the boundary layer and its subsequent transport into the free troposphere. The dynamic mechanism through which friction mitigates a barotropic vortex is that of Ekman pumping. This has often been assumed to be also the dominant mechanism for baroclinic developments. The PV framework highlights an additional, baroclinic mechanism. Positive PV is generated baroclinically due to friction to the north-east of a surface low and is transported out of the boundary layer by a cyclonic conveyor belt flow. The result is an anomaly of increased static stability in the lower troposphere which restricts the growth of the baroclinic wave. The reduced coupling between lower and upper levels can be sufficient to change the character of the upper-level evolution of the mature wave. The basic features of the baroclinic damping mechanism are robust for different frontal structures, with and without turbulent heat fluxes, and for the range of surface roughness found over the oceans.

A case study of boundary layer ventilation by convection and coastal processes

It is often assumed that ventilation of the atmospheric boundary layer is weak in the absence of fronts, but is this always true? In this paper we investigate the processes responsible for ventilation of the atmospheric boundary layer during a nonfrontal day that occurred on 9 May 2005 using the UK Met Office Unified Model. Pollution sources are represented by the constant emission of a passive tracer everywhere over land. The ventilation processes observed include shallow convection, turbulent mixing followed by large-scale ascent, a sea breeze circulation and coastal outflow. Vertical distributions of tracer are validated qualitatively with AMPEP (Aircraft Measurement of chemical Processing Export fluxes of Pollutants over the UK) CO aircraft measurements and are shown to agree impressively well. Budget calculations of tracers are performed in order to determine the relative importance of these ventilation processes. Coastal outflow and the sea breeze circulation were found to ventilate 26% of the boundary layer tracer by sunset of which 2% was above 2 km. A combination of coastal outflow, the sea breeze circulation, turbulent mixing and large-scale ascent ventilated 46% of the boundary layer tracer, of which 10% was above 2 km. Finally, coastal outflow, the sea breeze circulation, turbulent mixing, large-scale ascent and shallow convection together ventilated 52% of the tracer into the free troposphere, of which 26% was above 2 km. Hence this study shows that significant ventilation of the boundary layer can occur in the absence of fronts (and thus during high-pressure events). Turbulent mixing and convection processes can double the amount of pollution ventilated from the boundary layer.