On integral equation and least squares methods for scattering by diffraction gratings

In this paper we consider the scattering of a plane acoustic or electromagnetic wave by a one-dimensional, periodic rough surface. We restrict the discussion to the case when the boundary is sound soft in the acoustic case, perfectly reflecting with TE polarization in the EM case, so that the total field vanishes on the boundary. We propose a uniquely solvable first kind integral equation formulation of the problem, which amounts to a requirement that the normal derivative of the Green's representation formula for the total field vanish on a horizontal line below the scattering surface. We then discuss the numerical solution by Galerkin's method of this (ill-posed) integral equation. We point out that, with two particular choices of the trial and test spaces, we recover the so-called SC (spectral-coordinate) and SS (spectral-spectral) numerical schemes of DeSanto et al., Waves Random Media, 8, 315-414 1998. We next propose a new Galerkin scheme, a modification of the SS method that we term the SS* method, which is an instance of the well-known dual least squares Galerkin method. We show that the SS* method is always well-defined and is optimally convergent as the size of the approximation space increases. Moreover, we make a connection with the classical least squares method, in which the coefficients in the Rayleigh expansion of the solution are determined by enforcing the boundary condition in a least squares sense, pointing out that the linear system to be solved in the SS* method is identical to that in the least squares method. Using this connection we show that (reflecting the ill-posed nature of the integral equation solved) the condition number of the linear system in the SS* and least squares methods approaches infinity as the approximation space increases in size. We also provide theoretical error bounds on the condition number and on the errors induced in the numerical solution computed as a result of ill-conditioning. Numerical results confirm the convergence of the SS* method and illustrate the ill-conditioning that arises.

A comparison of conservative upwind difference schemes for the Euler equations

Numerical results are presented and compared for three conservative upwind difference schemes for the Euler equations when applied to two standard test problems. This includes consideration of the effect of treating part of the flux balance as a source, and a comparison of different averaging of the flow variables. Two of the schemes are also shown to be equivalent in their implementation, while being different in construction and having different approximate Jacobians. (C) 2006 Elsevier Ltd. All rights reserved.

It is premature to include non-CO2 effects of aviation in emission trading schemes

The recent G8 Gleneagles climate statement signed on 8 July 2005 specifically mentions a determination to lessen the impact of aviation on climate [Gleneagles, 2005. The Gleneagles communique: climate change, energy and sustainable development. http://www.fco.gov.uk/Files/kfile/PostG8_Gleneagles_Communique.pdf]. In January 2005 the European Union Emission Trading Scheme (ETS) commenced operation as the largest multi-country, multi-sector ETS in the world, albeit currently limited only to CO2 emissions. At present the scheme makes no provision for aircraft emissions. However, the UK Government would like to see aircraft included in the ETS and plans to use its Presidencies of both the EU and G8 in 2005 to implement these schemes within the EU and perhaps internationally. Non-CO2 effects have been included in some policy-orientated studies of the impact of aviation but we argue that the inclusion of such effects in any such ETS scheme is premature; we specifically argue that use of the Radiative Forcing Index for comparing emissions from different sources is inappropriate and that there is currently no metric for such a purpose that is likely to enable their inclusion in the near future. (c) 2005 Elsevier Ltd. All rights reserved.

Metrics for solar wind prediction models: Comparison of empirical, hybrid, and physics-based schemes with 8 years of L1 observations

Space weather effects on technological systems originate with energy carried from the Sun to the terrestrial environment by the solar wind. In this study, we present results of modeling of solar corona-heliosphere processes to predict solar wind conditions at the L1 Lagrangian point upstream of Earth. In particular we calculate performance metrics for (1) empirical, (2) hybrid empirical/physics-based, and (3) full physics-based coupled corona-heliosphere models over an 8-year period (1995–2002). L1 measurements of the radial solar wind speed are the primary basis for validation of the coronal and heliosphere models studied, though other solar wind parameters are also considered. The models are from the Center for Integrated Space-Weather Modeling (CISM) which has developed a coupled model of the whole Sun-to-Earth system, from the solar photosphere to the terrestrial thermosphere. Simple point-by-point analysis techniques, such as mean-square-error and correlation coefficients, indicate that the empirical coronal-heliosphere model currently gives the best forecast of solar wind speed at 1 AU. A more detailed analysis shows that errors in the physics-based models are predominately the result of small timing offsets to solar wind structures and that the large-scale features of the solar wind are actually well modeled. We suggest that additional “tuning” of the coupling between the coronal and heliosphere models could lead to a significant improvement of their accuracy. Furthermore, we note that the physics-based models accurately capture dynamic effects at solar wind stream interaction regions, such as magnetic field compression, flow deflection, and density buildup, which the empirical scheme cannot.

Atypical neural self-representation in autism

The 'self' is a complex multidimensional construct deeply embedded and in many ways defined by our relations with the social world. Individuals with autism are impaired in both self-referential and other-referential social cognitive processing. Atypical neural representation of the self may be a key to understanding the nature of such impairments. Using functional magnetic resonance imaging we scanned adult males with an autism spectrum condition and age and IQ-matched neurotypical males while they made reflective mentalizing or physical judgements about themselves or the British Queen. Neurotypical individuals preferentially recruit the middle cingulate cortex and ventromedial prefrontal cortex in response to self compared with other-referential processing. In autism, ventromedial prefrontal cortex responded equally to self and other, while middle cingulate cortex responded more to other-mentalizing than self-mentalizing. These atypical responses occur only in areas where self-information is preferentially processed and does not affect areas that preferentially respond to other-referential information. In autism, atypical neural self-representation was also apparent via reduced functional connectivity between ventromedial prefrontal cortex and areas associated with lower level embodied representations, such as ventral premotor and somatosensory cortex. Furthermore, the magnitude of neural self-other distinction in ventromedial prefrontal cortex was strongly related to the magnitude of early childhood social impairments in autism. Individuals whose ventromedial prefrontal cortex made the largest distinction between mentalizing about self and other were least socially impaired in early childhood, while those whose ventromedial prefrontal cortex made little to no distinction between mentalizing about self and other were the most socially impaired in early childhood. These observations reveal that the atypical organization of neural circuitry preferentially coding for self-information is a key mechanism at the heart of both self-referential and social impairments in autism.