A combined BIE-FE method for the Stokes equations

A numerical algorithm for the biharmonic equation in domains with piecewise smooth boundaries is presented. It is intended for problems describing the Stokes flow in the situations where one has corners or cusps formed by parts of the domain boundary and, due to the nature of the boundary conditions on these parts of the boundary, these regions have a global effect on the shape of the whole domain and hence have to be resolved with sufficient accuracy. The algorithm combines the boundary integral equation method for the main part of the flow domain and the finite-element method which is used to resolve the corner/cusp regions. Two parts of the solution are matched along a numerical ‘internal interface’ or, as a variant, two interfaces, and they are determined simultaneously by inverting a combined matrix in the course of iterations. The algorithm is illustrated by considering the flow configuration of ‘curtain coating’, a flow where a sheet of liquid impinges onto a moving solid substrate, which is particularly sensitive to what happens in the corner region formed, physically, by the free surface and the solid boundary. The ‘moving contact line problem’ is addressed in the framework of an earlier developed interface formation model which treats the dynamic contact angle as part of the solution, as opposed to it being a prescribed function of the contact line speed, as in the so-called ‘slip models’. Keywords: Dynamic contact angle; finite elements; free surface flows; hybrid numerical technique; Stokes equations.

In-street wind direction variability in the vicinity of a busy intersection in central London

We present results from fast-response wind measurements within and above a busy intersection between two street canyons (Marylebone Road and Gloucester Place) in Westminster, London taken as part of the DAPPLE (Dispersion of Air Pollution and Penetration into the Local Environment; www.dapple.org.uk) 2007 field campaign. The data reported here were collected using ultrasonic anemometers on the roof-top of a building adjacent to the intersection and at two heights on a pair of lamp-posts on opposite sides of the intersection. Site characteristics, data analysis and the variation of intersection flow with the above-roof wind direction (θref) are discussed. Evidence of both flow channelling and recirculation was identified within the canyon, only a few metres from the intersection for along-street and across-street roof-top winds respectively. Results also indicate that for oblique rooftop flows, the intersection flow is a complex combination of bifurcated channelled flows, recirculation and corner vortices. Asymmetries in local building geometry around the intersection and small changes in the background wind direction (changes in 15-min mean θref of 5–10 degrees) were also observed to have profound influences on the behaviour of intersection flow patterns. Consequently, short time-scale variability in the background flow direction can lead to highly scattered in-street mean flow angles masking the true multi-modal features of the flow and thus further complicating modelling challenges.

Catena-Poly[bis(ethane-1,2-diammonium) [manganese(II)-di-mu-phosphato-kappa O-4 : O ']]: a one-dimensional manganese phosphate

The title compound,{(C2H10N2)(2)[Mn(PO4)(2)]}(n), contains anionic square-twisted chains of formula [Mn(PO4)(2)](4-) constructed from corner-sharing four-membered rings of alternating MnO4 and PO4 units. The Mn and P atoms have distorted tetrahedral coordination and the Mn atom lies on a twofold axis. The linear manganese-phosphate chains are held together by hydrogen-bonding interactions involving the framework O atoms and the H atoms of the ethane-1,2-diammonium cations, which lie in the interchain spaces.

The synthesis and characterisation of two new iron thioantimonates: [Fe(en)(3)](2)Sb2S5 center dot 0.55H(2)O and [Fe(en)(3)](2)Sb4S8

Two new iron thioantimonates, [Fe(en)(3)](2)Sb2S5 (.) 0.55H(2)O (1) and [Fe(en)(3)](2)Sb4S8 (2). were synthesised under solvothermal conditions from the reactions of Sb2S3, FeCl2 and S in the presence of ethylenediamine at 413 and 438 K, respectively. The products were characterised by single-crystal X-ray diffraction, elemental analysis and SQUID magnetometry. Compound 1 is unusual in containing isolated Sb2S54- anions formed from two corner-sharing SbS33- trigonal pyramids. These units are arranged in rippled layers, 4 A apart, parallel to the bc-plane. Octahedrally coordinated [Fe(en)(3)](2+) cations lie in depressions within these anionic layers. In compound (2), repeated corner linking of SbS33- trigonal pyramids generates SbS2- chains, which may be considered as a polymerised form of the Sb2S54- anions in 1. The SbS2- chains are separated by [Fe(en)(3)](2+) cations. In both compounds, there is an extensive network of hydrogen bonds between the nitrogen atoms of the ethylenediamine ligands and the sulfur atoms of the anions and, in the case of 1, the uncoordinated water molecule. (c) 2005 Elsevier Ltd. All rights reserved.

Robust classification of SAR imagery

In this work the G(A)(0) distribution is assumed as the universal model for amplitude Synthetic Aperture (SAR) imagery data under the Multiplicative Model. The observed data, therefore, is assumed to obey a G(A)(0) (alpha; gamma, n) law, where the parameter n is related to the speckle noise, and (alpha, gamma) are related to the ground truth, giving information about the background. Therefore, maps generated by the estimation of (alpha, gamma) in each coordinate can be used as the input for classification methods. Maximum likelihood estimators are derived and used to form estimated parameter maps. This estimation can be hampered by the presence of corner reflectors, man-made objects used to calibrate SAR images that produce large return values. In order to alleviate this contamination, robust (M) estimators are also derived for the universal model. Gaussian Maximum Likelihood classification is used to obtain maps using hard-to-deal-with simulated data, and the superiority of robust estimation is quantitatively assessed.

Wiener-System subspace identification for mobile wireless mm-wave networks

We discuss the feasibility of wireless terahertz communications links deployed in a metropolitan area and model the large-scale fading of such channels. The model takes into account reception through direct line of sight, ground and wall reflection, as well as diffraction around a corner. The movement of the receiver is modeled by an autonomous dynamic linear system in state space, whereas the geometric relations involved in the attenuation and multipath propagation of the electric field are described by a static nonlinear mapping. A subspace algorithm in conjunction with polynomial regression is used to identify a single-output Wiener model from time-domain measurements of the field intensity when the receiver motion is simulated using a constant angular speed and an exponentially decaying radius. The identification procedure is validated by using the model to perform q-step ahead predictions. The sensitivity of the algorithm to small-scale fading, detector noise, and atmospheric changes are discussed. The performance of the algorithm is tested in the diffraction zone assuming a range of emitter frequencies (2, 38, 60, 100, 140, and 400 GHz). Extensions of the simulation results to situations where a more complicated trajectory describes the motion of the receiver are also implemented, providing information on the performance of the algorithm under a worst case scenario. Finally, a sensitivity analysis to model parameters for the identified Wiener system is proposed.

An Exponentially Convergent Nonpolynomial Finite Element Method for Time-Harmonic Scattering from Polygons

In recent years nonpolynomial finite element methods have received increasing attention for the efficient solution of wave problems. As with their close cousin the method of particular solutions, high efficiency comes from using solutions to the Helmholtz equation as basis functions. We present and analyze such a method for the scattering of two-dimensional scalar waves from a polygonal domain that achieves exponential convergence purely by increasing the number of basis functions in each element. Key ingredients are the use of basis functions that capture the singularities at corners and the representation of the scattered field towards infinity by a combination of fundamental solutions. The solution is obtained by minimizing a least-squares functional, which we discretize in such a way that a matrix least-squares problem is obtained. We give computable exponential bounds on the rate of convergence of the least-squares functional that are in very good agreement with the observed numerical convergence. Challenging numerical examples, including a nonconvex polygon with several corner singularities, and a cavity domain, are solved to around 10 digits of accuracy with a few seconds of CPU time. The examples are implemented concisely with MPSpack, a MATLAB toolbox for wave computations with nonpolynomial basis functions, developed by the authors. A code example is included.

On the existence of extreme waves and the Stokes conjecture with vorticity

This is a study of singular solutions of the problem of traveling gravity water waves on flows with vorticity. We show that, for a certain class of vorticity functions, a sequence of regular waves converges to an extreme wave with stagnation points at its crests. We also show that, for any vorticity function, the profile of an extreme wave must have either a corner of 120° or a horizontal tangent at any stagnation point about which it is supposed symmetric. Moreover, the profile necessarily has a corner of 120° if the vorticity is nonnegative near the free surface.

Radiometric analysis of the light coupled by optimally cut plastic optical fiber amplitude modulating reflectance displacement sensors

A radiometric analysis of the light coupled by optical fiber amplitude modulating extrinsic-type reflectance displacement sensors is presented. Uncut fiber sensors show the largest range but a smaller responsivity. Single cut fiber sensors exhibit an improvement in responsivity at the expense of range. A further increase in responsivity as well as a reduction in the operational range is obtained when the double cut sensor configuration is implemented. The double cut configuration is particularly suitable in applications where feedback action is applied to the moving reflector surface. © 2000 American Institute of Physics.

Feedback dew-point sensor utilising optimally cut plastic optical fibres

A plastic optical fibre reflectance sensor that makes full use of the critical angle of the fibres is implemented to monitor dew formation on a Peltier-cooled reflector surface. The optical configuration permits isolation of optoelectronic components from the sensing head and better light coupling between the reflector and the detecting fibre, giving a better signal of the onset of dew formation on the reflector. Continuous monitoring of the rate of change in reflectance as well as the absolute reflectance signals, the use of a novel polymethyl-methacrylate-coated hydrophobic film reflector on the Peltier element and the application of feedback around the point of dew formation, further reduces the possibility of contamination of the sensor head. Under closed-loop operation, the sensor is capable of cycling around the point of dew formation at a frequency of 2.5 Hz.

The Stokes conjecture for waves with vorticity

We study stagnation points of two-dimensional steady gravity free-surface water waves with vorticity. We obtain for example that, in the case where the free surface is an injective curve, the asymptotics at any stagnation point is given either by the “Stokes corner flow” where the free surface has a corner of 120°, or the free surface ends in a horizontal cusp, or the free surface is horizontally flat at the stagnation point. The cusp case is a new feature in the case with vorticity, and it is not possible in the absence of vorticity. In a second main result we exclude horizontally flat singularities in the case that the vorticity is 0 on the free surface. Here the vorticity may have infinitely many sign changes accumulating at the free surface, which makes this case particularly difficult and explains why it has been almost untouched by research so far. Our results are based on calculations in the original variables and do not rely on structural assumptions needed in previous results such as isolated singularities, symmetry and monotonicity.

Higher education as a portfolio investment: students’ choices about studying, term time employment, leisure, and loans

Recent UK changes in the number of students entering higher education, and in the nature of financial support, highlight the complexity of students’ choices about human capital investments. Today’s students have to focus not on the relatively narrow issue of how much academic effort to invest, but instead on the more complicated issue of how to invest effort in pursuit of ‘employability skills’, and how to signal such acquisitions in the context of a highly competitive graduate jobs market. We propose a framework aimed specifically at students’ investment decisions, which encompasses corner solutions for both borrowing and employment while studying.

Co3(SO4)3(OH)2[enH2]: a newS= 3/2 Kagome-type layered sulfate with a unique connectivity

The Kagome lattice, comprising a two-dimensional array of corner-sharing equilateral triangles, is central to the exploration of magnetic frustration. In such a lattice, antiferromagnetic coupling between ions in triangular plaquettes prevents all of the exchange interactions being simultaneously satisfied and a variety of novel magnetic ground states may result at low temperature. Experimental realization of a Kagome lattice remains difficult. The jarosite family of materials of nominal composition AM3(SO4)2(OH)6 (A = monovalent cation; M= Fe3+, Cr3+), offers perhaps one of the most promising manifestations of the phenomenon of magnetic frustration in two dimensions. The magnetic properties of jarosites are however extremely sensitive to the degree of coverage of magnetic sites. Consequently, there is considerable interest in the use of soft chemical techniques for the design and synthesis of novel materials in which to explore the effects of spin, degree of site coverage and connectivity on magnetic frustration.

A three-dimensional open-framework indium selenide: [C7H10N][In9Se14]

An open-framework indium selenide, [C7H10N][In9Se14], has been prepared under solvothermal conditions in the presence of 3,5-dimethylpyridine, and characterized by single crystal diffraction, thermogravimetry, elemental analysis, FTIR spectroscopy and UV-Vis diffuse reflectance. The crystal structure of [C7H10N][In9Se14] contains an unusual building unit, in which corner-linked and edge-linked InSe45- tetrahedra coexist. The presence of one-dimensional circular channels, of ca. 6 Å diameter, results in approximately 25% of solvent accessible void space.