An efficient shock capturing algorithm for compressible flows in a duct of variable cross-section

A numerical scheme is presented for the solution of the Euler equations of compressible flow of a gas in a single spatial co-ordinate. This includes flow in a duct of variable cross-section as well as flow with slab, cylindrical or spherical symmetry and can prove useful when testing codes for the two-dimensional equations governing compressible flow of a gas. The resulting scheme requires an average of the flow variables across the interface between cells and for computational efficiency this average is chosen to be the arithmetic mean, which is in contrast to the usual ‘square root’ averages found in this type of scheme. The scheme is applied with success to five problems with either slab or cylindrical symmetry and a comparison is made in the cylindrical case with results from a two-dimensional problem with no sources.

Acoustic trapped modes in a three-dimensional waveguide of slowly varying cross section

In this paper we develop an asymptotic scheme to approximate the trapped mode solutions to the time harmonic wave equation in a three-dimensional waveguide with a smooth but otherwise arbitrarily shaped cross section and a single, slowly varying `bulge', symmetric in the longitudinal direction. Extending the work in Biggs (2012), we first employ a WKBJ-type ansatz to identify the possible quasi-mode solutions which propagate only in the thicker region, and hence find a finite cut-on region of oscillatory behaviour and asymptotic decay elsewhere. The WKBJ expansions are used to identify a turning point between the cut-on and cut-on regions. We note that the expansions are nonuniform in an interior layer centred on this point, and we use the method of matched asymptotic expansions to connect the cut-on and cut-on regions within this layer. The behaviour of the expansions within the interior layer then motivates the construction of a uniformly valid asymptotic expansion. Finally, we use this expansion and the symmetry of the waveguide around the longitudinal centre, x = 0, to extract trapped mode wavenumbers, which are compared with those found using a numerical scheme and seen to be extremely accurate, even to relatively large values of the small parameter.

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.

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.

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.

A transmission electron microscopy study of polymer-stabilised liquid crystal structure

A method has been established for observing the internal structure of the network component of polymer-stabilised liquid crystals. In situ photopolymerisation of a mesogenic diacrylate monomer using ultraviolet light leads to a sparse network (∼1 wt%) within a nematic host. Following polymerisation, the host was removed through dissolution in heptane, revealing the network. In order to observe a cross-section through the network, it was embedded in a resin and then sectioned using an ultramicrotome. However, imaging of the network was not possible due to poor contrast. To improve this, several reagents were used for network staining, but only one was successful: bromine. The use of a Melinex-resin composite for sectioning was also found to be advantageous. Imaging of the network using transmission electron microscopy revealed solid “droplets” of width 0.07–0.20 μm, possessing an open, yet homogeneous structure, with no evidence for any large-scale internal structures.

Resolution of the uncertainties in the radiative forcing of HFC-134a

HFC-134a (CF3CH2F) is the most rapidly growing hydrofluorocarbon in terms of atmospheric abundance. It is currently used in a large number of household refrigerators and air-conditioning systems and its concentration in the atmosphere is forecast to increase substantially over the next 50–100 years. Previous estimates of its radiative forcing per unit concentration have differed significantly 25%. This paper uses a two-step approach to resolve this discrepancy. In the first step six independent absorption cross section datasets are analysed. We find that, for the integrated cross section in the spectral bands that contribute most to the radiative forcing, the differences between the various datasets are typically smaller than 5% and that the dependence on pressure and temperature is not significant. A “recommended'' HFC-134a infrared absorption spectrum was obtained based on the average band intensities of the strongest bands. In the second step, the “recommended'' HFC-134a spectrum was used in six different radiative transfer models to calculate the HFC-134a radiative forcing efficiency. The clear-sky instantaneous radiative forcing, using a single global and annual mean profile, differed by 8%, between the 6 models, and the latitudinally-resolved adjusted cloudy sky radiative forcing estimates differed by a similar amount.

An efficient flux difference splitting algorithm for unsteady duct flows of a real gas

A numerical scheme is presented for the solution of the Euler equations of compressible flow of a real gas in a single spatial coordinate. This includes flow in a duct of variable cross-section, as well as flow with slab, cylindrical or spherical symmetry, as well as the case of an ideal gas, and can be useful when testing codes for the two-dimensional equations governing compressible flow of a real gas. The resulting scheme requires an average of the flow variables across the interface between cells, and this average is chosen to be the arithmetic mean for computational efficiency, which is in contrast to the usual “square root” averages found in this type of scheme. The scheme is applied with success to five problems with either slab or cylindrical symmetry and for a number of equations of state. The results compare favourably with the results from other schemes.

Real gas flows in a duct

A numerical scheme is presented for the solution of the Euler equations of compressible flow of a real gas in a single spatial coordinate. This include flow in a duct of variable cross-section as well as flow with cylindrical or spherical symmetry, and can prove useful when testing codes for the two-dimensional equations governing compressible flow of a real gas. The scheme is applied with success to a problem involving the interaction of converging and diverging cylindrical shocks for four equations of state and to a problem involving the reflection of a converging shock.

A study of the atmospherically relevant reaction between molecular chlorine and dimethylsulfide (DMS): establishing the reaction intermediate and measurement of absolute photoionization cross-sections

In a study using UV photoelectron spectroscopy (PES) of the atmospherically relevant reaction CH3SCH3 + Cl2 → CH3SCH2Cl + HCl bands associated with a reaction intermediate have been observed. These have been assigned to ionization of the covalently bound molecule (CH3)2SCl2 on the basis of the intensity of the observed bands as a function of reaction time, molecular orbital calculations of vertical ionization energies and evidence from infrared spectroscopy. A method has also been developed, with the flow-tube/PE spectrometer combination used, to measure photoionization cross-sections of the reagents and products at the photon energy utilized and this has allowed the photoionization cross-section of the intermediate to be estimated. This work augments an earlier study in which the rate constant of the reaction between CH3SCH3 (DMS) and Cl2 has been measured at room temperature.

Measurement of the salinity and freezing tolerance of Crataegus genotypes using chlorophyll fluorescence

The effect of increasing salinity and freezing stress singly and in combination on a range of chlorophyll fluorescence parameters in foliar tissue of six Crataegus genotypes was examined. In general, increased stress reduced fluorescence values and absorption, trapping and electron transport energy fluxes per leaf reaction center and cross section, with decreased sigmoidicity of OJIP curves as a measure of the plastoquinone pool, reflecting decreased energy fluxes. Based on percentage reduction in a performance index from controls compared to stress-treated values, plants were ranked in order of tolerant > intermediate > sensitive. Use of this PIp ranking criteria enabled the distinguishing of marked differences in foliar salt/freezing hardiness between the Crataegus species used. Interpretation of the photochemical data showed that salinity and freezing affects both the acceptor and donor side of Photosystem II, while OJIP observations provided information regarding structural and functional changes in the leaf photosynthetic apparatus of the test species. It is concluded that chlorophyll fluorescence offers a rapid screening technique for assessing foliar salinity and freezing tolerance of woody perennials

Distribution of β-amylase I haplotypes among European cultivated barleys

The barley β-amylase I (Bmy1) locus encodes a starch breakdown enzyme whose kinetic properties and thermostability are critical during malt production. Studies of allelic variation at the Bmy1 locus have shown that the encoded enzyme can be commonly found in at least three distinct thermostability classes and demonstrated the nucleotide sequence variations responsible for such phenotypic differences. In order to explore the extent of sequence diversity at the Bmy1 locus in cultivated European barley, 464 varieties representing a cross-section of popular varieties grown in western Europe over the past 60 years, were genotyped for three single nucleotide polymorphisms chosen to tag the four common alleles found in the collection. One of these haplotypes, which has not been explicitly recognised in the literature as a distinct allele, was found in 95% of winter varieties in the sample. When release dates of the varieties were considered, the lowest thermostability allele (Bmy1-Sd2L) appeared to decrease in abundance over time, while the highest thermostability allele (Bmy1-Sd2H) was the rarest allele at 5.4% of the sample and was virtually confined to two-row spring varieties. Pedigree analysis was used to track transmission of particular alleles over time and highlighted issues of genetic stratification of the sample.

Photoproduction of heavy quarks in next-to-leading order QCD with longitudinally polarized initial states

We present all relevant details of our calculation of the complete next-to-leading order O(αS2α) QCD corrections to heavy flavor photoproduction with longitudinally polarized point-like photons and hadrons. In particular we provide analytical results for the virtual plus soft gluon cross section. We carefully address the relevance of remaining theoretical uncertainties by varying, for instance, the factorization and renormalization scales independently. Such studies are of importance for a meaningful first direct determination of the polarized gluon density Δg from the total charm production spin asymmetry by the upcoming COMPASS experiment. It is shown that the scale uncertainty is considerably reduced in next-to-leading order, but the dependence on the charm quark mass is sizable at fixed target energies. Finally, we study several differential single-inclusive heavy quark distributions and, for the polarized HERA option, the total bottom spin asymmetry.