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.

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.

Noise propagation from a cutting of arbitrary cross-section and impedance

A boundary integral equation is described for the prediction of acoustic propagation from a monofrequency coherent line source in a cutting with impedance boundary conditions onto surrounding flat impedance ground. The problem is stated as a boundary value problem for the Helmholtz equation and is subsequently reformulated as a system of boundary integral equations via Green's theorem. It is shown that the integral equation formulation has a unique solution at all wavenumbers. The numerical solution of the coupled boundary integral equations by a simple boundary element method is then described. The convergence of the numerical scheme is demonstrated experimentally. Predictions of A-weighted excess attenuation for a traffic noise spectrum are made illustrating the effects of varying the depth of the cutting and the absorbency of the surrounding ground surface.

On asymptotic behavior at infinity and the finite section method for integral equations on the half-line

e consider integral equations on the half-line of the form and the finite section approximation to x obtained by replacing the infinite limit of integration by the finite limit β. We establish conditions under which, if the finite section method is stable for the original integral equation (i.e. exists and is uniformly bounded in the space of bounded continuous functions for all sufficiently large β), then it is stable also for a perturbed equation in which the kernel k is replaced by k + h. The class of perturbations allowed includes all compact and some non-compact perturbations of the integral operator. Using this result we study the stability and convergence of the finite section method in the space of continuous functions x for which ()()()=−∫∞dttxt,sk)s(x0()syβxβx()sxsp+1 is bounded. With the additional assumption that ()(tskt,sk−≤ where ()()(),qsomefor,sassOskandRLkq11>+∞→=∈− we show that the finite-section method is stable in the weighted space for ,qp≤≤0 provided it is stable on the space of bounded continuous functions. With these results we establish error bounds in weighted spaces for x - xβ and precise information on the asymptotic behaviour at infinity of x. We consider in particular the case when the integral operator is a perturbation of a Wiener-Hopf operator and illustrate this case with a Wiener-Hopf integral equation arising in acoustics.

Late Barremian-Early Aptian palaeoenvironmental change: the Cassis-La Bedoule section (SE France)

The well-dated section of Cassis-La Bédoule in the South Provencal Basin (southern France) allows for a detailed reconstruction of palaeoenvironmental change during the latest Barremian and Early Aptian. For this study, phosphorus (P) and clay-mineral contents, stable-isotope ratios on carbonate (δ13Ccarb) and organic matter (δ13Corg), and redox-sensitive trace elements (RSTE: V, U, As, Co, and Mo) have been measured in this historical stratotype. The base of the section consists of rudist limestone, which is attributed to the Urgonian platform. The presence of low P and RSTE content, and content of up to 30% kaolinite indicate deposition under oligotrophic and oxic conditions, and the presence of warm, humid climatic conditions on the adjacent continent. The top of the Urgonian succession is marked by a hardground with encrusted brachiopods and bivalves, which is interpreted as a drowning surface. The section continues with a succession of limestone and marl containing the first occurrence of planktonic foraminifera. This interval includes several laminated, organic-rich layers recording RSTE enrichments and high Corg:Ptot ratios. The deposition of these organic-rich layers was associated with oxygen-depleted conditions and a large positive excursion in δ13Corg. During this interval, a negative peak in the δ13Ccarb record is observed, which dates as latest Barremian. This excursion is coeval with negative excursions elsewhere in Tethyan platform and basin settings and is explained by the increased input of light dissolved inorganic carbon by rivers and/or volcanic activity. In this interval, an increase in P content, owing to reworking of nearshore sediments during the transgression, is coupled with a decrease in kaolinite content, which tends to be deposited in more proximal areas. The overlying hemipelagic sediments of the Early Aptian Deshayesites oglanlensis and D. weissi zones indicate rather stable palaeoenvironmental conditions with low P content and stable δ13C records. A change towards marl-dominated beds occurs close to the end of the D. weissi zone. These beds display a long decrease in their δ13Ccarb and δ13Corg records, which lasted until the end of the Deshayesites deshayesi subzone (corresponding to C3 in Menegatti et al., 1998). This is followed by a positive shift during the Roloboceras hambrovi and Deshayesites grandis subzones, which corresponds in time to oceanic anoxic event (OAE) 1a interval. This positive shift is coeval with two increases in the P content. The marly interval equivalent to OAE 1a lacks organic-rich deposits and RSTE enrichments indicating that oxic conditions prevailed in this particular part of the Tethys ocean. The clay mineralogy is dominated by smectite, which is interpreted to reflect trapping of kaolinite on the surrounding platforms rather than indicating a drier climate.

Efficient incorporation of channel cross-section geometry uncertainty into regional and global scale flood inundation models

This paper investigates the challenge of representing structural differences in river channel cross-section geometry for regional to global scale river hydraulic models and the effect this can have on simulations of wave dynamics. Classically, channel geometry is defined using data, yet at larger scales the necessary information and model structures do not exist to take this approach. We therefore propose a fundamentally different approach where the structural uncertainty in channel geometry is represented using a simple parameterization, which could then be estimated through calibration or data assimilation. This paper first outlines the development of a computationally efficient numerical scheme to represent generalised channel shapes using a single parameter, which is then validated using a simple straight channel test case and shown to predict wetted perimeter to within 2% for the channels tested. An application to the River Severn, UK is also presented, along with an analysis of model sensitivity to channel shape, depth and friction. The channel shape parameter was shown to improve model simulations of river level, particularly for more physically plausible channel roughness and depth parameter ranges. Calibrating channel Manning’s coefficient in a rectangular channel provided similar water level simulation accuracy in terms of Nash-Sutcliffe efficiency to a model where friction and shape or depth were calibrated. However, the calibrated Manning coefficient in the rectangular channel model was ~2/3 greater than the likely physically realistic value for this reach and this erroneously slowed wave propagation times through the reach by several hours. Therefore, for large scale models applied in data sparse areas, calibrating channel depth and/or shape may be preferable to assuming a rectangular geometry and calibrating friction alone.

Equation for the microwave backscatter cross section of aggregate snowflakes using the self-similar Rayleigh–Gans approximation

In this paper an equation is derived for the mean backscatter cross section of an ensemble of snowflakes at centimeter and millimeter wavelengths. It uses the Rayleigh–Gans approximation, which has previously been found to be applicable at these wavelengths due to the low density of snow aggregates. Although the internal structure of an individual snowflake is random and unpredictable, the authors find from simulations of the aggregation process that their structure is “self-similar” and can be described by a power law. This enables an analytic expression to be derived for the backscatter cross section of an ensemble of particles as a function of their maximum dimension in the direction of propagation of the radiation, the volume of ice they contain, a variable describing their mean shape, and two variables describing the shape of the power spectrum. The exponent of the power law is found to be −. In the case of 1-cm snowflakes observed by a 3.2-mm-wavelength radar, the backscatter is 40–100 times larger than that of a homogeneous ice–air spheroid with the same mass, size, and aspect ratio.

Commodity risks and the cross-section of equity returns

The article examines whether commodity risk is priced in the cross-section of global equity returns. We employ a long-only equally-weighted portfolio of commodity futures and a term structure portfolio that captures phases of backwardation and contango as mimicking portfolios for commodity risk. We find that equity-sorted portfolios with greater sensitivities to the excess returns of the backwardation and contango portfolio command higher average excess returns, suggesting that when measured appropriately, commodity risk is pervasive in stocks. Our conclusions are robust to the addition to the pricing model of financial, macroeconomic and business cycle-based risk factors.