Dispersion of a point-source release of a passive scalar through an urban-like array for different wind directions

The dispersion of a point-source release of a passive scalar in a regular array of cubical, urban-like, obstacles is investigated by means of direct numerical simulations. The simulations are conducted under conditions of neutral stability and fully rough turbulent flow, at a roughness Reynolds number of Reτ = 500. The Navier–Stokes and scalar equations are integrated assuming a constant rate release from a point source close to the ground within the array. We focus on short-range dispersion, when most of the material is still within the building canopy. Mean and fluctuating concentrations are computed for three different pressure gradient directions (0◦ , 30◦ , 45◦). The results agree well with available experimental data measured in a water channel for a flow angle of 0◦ . Profiles of mean concentration and the three-dimensional structure of the dispersion pattern are compared for the different forcing angles. A number of processes affecting the plume structure are identified and discussed, including: (i) advection or channelling of scalar down ‘streets’, (ii) lateral dispersion by turbulent fluctuations and topological dispersion induced by dividing streamlines around buildings, (iii) skewing of the plume due to flow turning with height, (iv) detrainment by turbulent dispersion or mean recirculation, (v) entrainment and release of scalar in building wakes, giving rise to ‘secondary sources’, (vi) plume meandering due to unsteady turbulent fluctuations. Finally, results on relative concentration fluctuations are presented and compared with the literature for point source dispersion over flat terrain and urban arrays. Keywords Direct numerical simulation · Dispersion modelling · Urban array

Implementation of a new urban energy budget scheme in the MetUM. Part I: Description and idealized simulations

This paper describes the formulation of a new urban scheme, MORUSES (Met Office–Reading Urban Surface Exchange Scheme) for use in the Met Office Unified Model. The implementation of the new scheme ensures that (1) the new scheme offers more flexibility in the parametrization of the building properties, and hence provides a more realistic representation of the fluxes; (2) the bulk outputs are in satisfactory agreement with previous observational studies; and (3) the impact of the new scheme on the energy balance fluxes is similar to the impact of the current urban scheme when set up to mimic it. As well as having a better physical basis, MORUSES also gains in flexibility in applications and adaptations to different urban materials as well as urban planning. The new scheme represents the urban area as a composition of two tiles, a canyon and a roof, using a simple 2D geometry. Sensitivity analysis to canyon geometry and thickness of the roof canopy emphasizes the gain in flexibility captured by the new scheme. Copyright © 2010 Royal Meteorological Society and Crown Copyright

Three-dimensional random Voronoi tessellations: From cubic crystal lattices to Poisson point processes

We perturb the SC, BCC, and FCC crystal structures with a spatial Gaussian noise whose adimensional strength is controlled by the parameter a, and analyze the topological and metrical properties of the resulting Voronoi Tessellations (VT). The topological properties of the VT of the SC and FCC crystals are unstable with respect to the introduction of noise, because the corresponding polyhedra are geometrically degenerate, whereas the tessellation of the BCC crystal is topologically stable even against noise of small but finite intensity. For weak noise, the mean area of the perturbed BCC and FCC crystals VT increases quadratically with a. In the case of perturbed SCC crystals, there is an optimal amount of noise that minimizes the mean area of the cells. Already for a moderate noise (a>0.5), the properties of the three perturbed VT are indistinguishable, and for intense noise (a>2), results converge to the Poisson-VT limit. Notably, 2-parameter gamma distributions are an excellent model for the empirical of of all considered properties. The VT of the perturbed BCC and FCC structures are local maxima for the isoperimetric quotient, which measures the degre of sphericity of the cells, among space filling VT. In the BCC case, this suggests a weaker form of the recentluy disproved Kelvin conjecture. Due to the fluctuations of the shape of the cells, anomalous scalings with exponents >3/2 is observed between the area and the volumes of the cells, and, except for the FCC case, also for a->0. In the Poisson-VT limit, the exponent is about 1.67. As the number of faces is positively correlated with the sphericity of the cells, the anomalous scaling is heavily reduced when we perform powerlaw fits separately on cells with a specific number of faces.

Effects of Thomson-scattering geometry on white-light imaging of an interplanetary shock: synthetic observations from forward magnetohydrodynamic modelling

Stereoscopic white-light imaging of a large portion of the inner heliosphere has been used to track interplanetary coronal mass ejections. At large elongations from the Sun, the white-light brightness depends on both the local electron density and the efficiency of the Thomson-scattering process. To quantify the effects of the Thomson-scattering geometry, we study an interplanetary shock using forward magnetohydrodynamic simulation and synthetic white-light imaging. Identifiable as an inclined streak of enhanced brightness in a time–elongation map, the travelling shock can be readily imaged by an observer located within a wide range of longitudes in the ecliptic. Different parts of the shock front contribute to the imaged brightness pattern viewed by observers at different longitudes. Moreover, even for an observer located at a fixed longitude, a different part of the shock front will contribute to the imaged brightness at any given time. The observed brightness within each imaging pixel results from a weighted integral along its corresponding ray-path. It is possible to infer the longitudinal location of the shock from the brightness pattern in an optical sky map, based on the east–west asymmetry in its brightness and degree of polarisation. Therefore, measurement of the interplanetary polarised brightness could significantly reduce the ambiguity in performing three-dimensional reconstruction of local electron density from white-light imaging.

Estimation of the friction velocity in stably stratified boundary-layer flows over hills

A method is suggested for the calculation of the friction velocity for stable turbulent boundary-layer flow over hills. The method is tested using a continuous upstream mean velocity profile compatible with the propagation of gravity waves, and is incorporated into the linear model of Hunt, Leibovich and Richards with the modification proposed by Hunt, Richards and Brighton to include the effects of stability, and the reformulated solution of Weng for the near-surface region. Those theoretical results are compared with results from simulations using a non-hydrostatic microscale-mesoscale two-dimensional numerical model, and with field observations for different values of stability. These comparisons show a considerable improvement in the behaviour of the theoretical model when the friction velocity is calculated using the method proposed here, leading to a consistent variation of the boundary-layer structure with stability, and better agreement with observational and numerical data.

Characterization of the Salmonella typhimurium proteome by semi-automated two dimensional HPLC-mass spectrometry: detection of proteins implicated in multiple antibiotic resistance

The proteome of Salmonella enterica serovar Typhimurium was characterized by 2-dimensional HPLC mass spectrometry to provide a platform for subsequent proteomic investigations of low level multiple antibiotic resistance (MAR). Bacteria (2.15 +/- 0.23 x 10(10) cfu; mean +/- s.d.) were harvested from liquid culture and proteins differentially fractionated, on the basis of solubility, into preparations representative of the cytosol, cell envelope and outer membrane proteins (OMPs). These preparations were digested by treatment with trypsin and peptides separated into fractions (n = 20) by strong cation exchange chromatography (SCX). Tryptic peptides in each SCX fraction were further separated by reversed-phase chromatography and detected by mass spectrometry. Peptides were assigned to proteins and consensus rank listings compiled using SEQUEST. A total of 816 +/- 11 individual proteins were identified which included 371 +/- 33, 565 +/- 15 and 262 +/- 5 from the cytosolic, cell envelope and OMP preparations, respectively. A significant correlation was observed (r(2) = 0.62 +/- 0.10; P < 0.0001) between consensus rank position for duplicate cell preparations and an average of 74 +/- 5% of proteins were common to both replicates. A total of 34 outer membrane proteins were detected, 20 of these from the OMP preparation. A range of proteins (n = 20) previously associated with the mar locus in E. coli were also found including the key MAR effectors AcrA, TolC and OmpF.

Grid computing solutions for distributed repositories of protein folding and unfolding simulations

The P-found protein folding and unfolding simulation repository is designed to allow scientists to perform analyses across large, distributed simulation data sets. There are two storage components in P-found: a primary repository of simulation data and a data warehouse. Here we demonstrate how grid technologies can support multiple, distributed P-found installations. In particular we look at two aspects, first how grid data management technologies can be used to access the distributed data warehouses; and secondly, how the grid can be used to transfer analysis programs to the primary repositories --- this is an important and challenging aspect of P-found because the data volumes involved are too large to be centralised. The grid technologies we are developing with the P-found system will allow new large data sets of protein folding simulations to be accessed and analysed in novel ways, with significant potential for enabling new scientific discoveries.

P-found: grid-enabling distributed repositories of protein folding and unfolding simulations for data mining

The P-found protein folding and unfolding simulation repository is designed to allow scientists to perform data mining and other analyses across large, distributed simulation data sets. There are two storage components in P-found: a primary repository of simulation data that is used to populate the second component, and a data warehouse that contains important molecular properties. These properties may be used for data mining studies. Here we demonstrate how grid technologies can support multiple, distributed P-found installations. In particular, we look at two aspects: firstly, how grid data management technologies can be used to access the distributed data warehouses; and secondly, how the grid can be used to transfer analysis programs to the primary repositories — this is an important and challenging aspect of P-found, due to the large data volumes involved and the desire of scientists to maintain control of their own data. The grid technologies we are developing with the P-found system will allow new large data sets of protein folding simulations to be accessed and analysed in novel ways, with significant potential for enabling scientific discovery.

Structural variations in self-assembled coordination complexes of Zn(II) with hexamethylenetetramine and isomeric 2-, 3- and 4-nitrobenzoates

Three new zinc(II)-hexamethylenetetramine (hmt) complexes [Zn-2(4-nbz)(4)(mu(2)-hmt)(OH2)(hmt)] (1). [Zn-2(2-nbz)(4)(mu(2)-hmt)(2)](n) (2) and [Zn-3(3-nbz)(4)(mu(2)-hmt)(mu(2)-OH)(mu(3)-OH)](n) (3) with three isomeric nitrobenzoate, [4-nbz = 4-nitrobenzoate, 2-nbz = 2-nitrobenzoate and 3-nbz = 3-nitrobenzoate] have been synthesized and structurally characterized by X-ray crystallography. Their identities have also been established by elemental analysis: IR, NMR, UV-Vis and mass spectral studies. 1 is a dinuclear complex formed by bridging hmt with mu(2) coordinating mode. The geometry around the Zn centers in 1 is distorted tetrahedral. Paddle-wheel centrosymmetric Zn-2(2-nbz)(4) units of complex 2 are interconnected by mu(2)-hmt forming a one-dimensional chain with square-pyramidal geometries around the Zn centers. Compound 3 contains a mu(2)/mu(3)-hydroxido and mu(2)-hmt bridged 1D chain. In this complex, varied geometries around the Zn centers are observed viz, tetrahedral, square pyramidal and trigonal bipyramidal. Various weak forces, i.e. lone pair-pi, pi-pi and CH-pi interactions, play a key role in stabilizing the observed structures for complexes 1,2 and 3. This series of complexes demonstrates that although the nitro group does not coordinate to the metal center, its presence at the 2-, 3- or 4-position of the phenyl ring has a striking effect on the dimensionality as well as the structure of the resulted coordination polymers, probably due to the participation of the nitro group in 1.p.center dot center dot center dot pi and/or C-H center dot center dot center dot pi interactions.

The Asian summer monsoon: an intercomparison of CMIP5 vs. CMIP3 simulations of the late 20th century

The boreal summer Asian monsoon has been evaluated in 25 Coupled Model Intercomparison Project-5 (CMIP5) and 22 CMIP3 GCM simulations of the late 20th Century. Diagnostics and skill metrics have been calculated to assess the time-mean, climatological annual cycle, interannual variability, and intraseasonal variability. Progress has been made in modeling these aspects of the monsoon, though there is no single model that best represents all of these aspects of the monsoon. The CMIP5 multi-model mean (MMM) is more skillful than the CMIP3 MMM for all diagnostics in terms of the skill of simulating pattern correlations with respect to observations. Additionally, for rainfall/convection the MMM outperforms the individual models for the time mean, the interannual variability of the East Asian monsoon, and intraseasonal variability. The pattern correlation of the time (pentad) of monsoon peak and withdrawal is better simulated than that of monsoon onset. The onset of the monsoon over India is typically too late in the models. The extension of the monsoon over eastern China, Korea, and Japan is underestimated, while it is overestimated over the subtropical western/central Pacific Ocean. The anti-correlation between anomalies of all-India rainfall and Niño-3.4 sea surface temperature is overly strong in CMIP3 and typically too weak in CMIP5. For both the ENSO-monsoon teleconnection and the East Asian zonal wind-rainfall teleconnection, the MMM interannual rainfall anomalies are weak compared to observations. Though simulation of intraseasonal variability remains problematic, several models show improved skill at representing the northward propagation of convection and the development of the tilted band of convection that extends from India to the equatorial west Pacific. The MMM also well represents the space-time evolution of intraseasonal outgoing longwave radiation anomalies. Caution is necessary when using GPCP and CMAP rainfall to validate (1) the time-mean rainfall, as there are systematic differences over ocean and land between these two data sets, and (2) the timing of monsoon withdrawal over India, where the smooth southward progression seen in India Meteorological Department data is better realized in CMAP data compared to GPCP data.

Structural and magnetic studies of Schiff base complexes of nickel(ii) nitrite: change in crystalline state, ligand rearrangement and a very rare μ-nitrito-1κO:2κN:3κO′ bridging mode

Four new nickel(II) complexes, [Ni2L2(NO2)2]·CH2Cl2·C2H5OH, 2H2O (1), [Ni2L2(DMF)2(m-NO2)]ClO4·DMF (2a), [Ni2L2(DMF)2(m-NO2)]ClO4 (2b) and [Ni3L¢2(m3-NO2)2(CH2Cl2)]n·1.5H2O (3) where HL = 2-[(3-amino-propylimino)-methyl]-phenol, H2L¢ = 2-({3-[(2-hydroxy-benzylidene)-amino]-propylimino}-methyl)-phenol and DMF = N,N-dimethylformamide, have been synthesized starting with the precursor complex [NiL2]·2H2O, nickel(II) perchlorate and sodium nitrite and characterized structurally and magnetically. The structural analyses reveal that in all the complexes, NiII ions possess a distorted octahedral geometry. Complex 1 is a dinuclear di-m2-phenoxo bridged species in which nitrite ion acts as chelating co-ligand. Complexes 2a and 2b also consist of dinuclear entities, but in these two compounds a cis-(m-nitrito-1kO:2kN) bridge is present in addition to the di-m2-phenoxo bridge. The molecular structures of 2a and 2b are equivalent; they differ only in that 2a contains an additional solvated DMF molecule. Complex 3 is formed by ligand rearrangement and is a one-dimensional polymer in which double phenoxo as well as m-nitrito-1kO:2kN bridged trinuclear units are linked through a very rare m3-nitrito-1kO:2kN:3kO¢ bridge. Analysis of variable-temperature magnetic susceptibility data indicates that there is a global weak antiferromagnetic interaction between the nickel(II) ions in four complexes, with exchange parameters J of -5.26, -11.45, -10.66 and -5.99 cm-1 for 1, 2a, 2b and 3, respectively

Anion-directed template synthesis and hydrolysis of mono-condensed Schiff base of 1,3-pentanediamine ando-hydroxyacetophenone in NiII and CuII Complexes

Bis(o-hydroxyacetophenone)nickel(II) dihydrate, on reaction with 1,3-pentanediamine, yields a bis-chelate complex [NiL2]·2H2O (1) of mono-condensed tridentate Schiff baseligand HL {2-[1-(3-aminopentylimino)ethyl]phenol}. The Schiff base has been freed from the complex by precipitating the NiII as a dimethylglyoximato complex. HL reacts smoothly with Ni(SCN)2·4H2O furnishing the complex [NiL(NCS)] (2) and with CuCl2·2H2O in the presence of NaN3 or NH4SCN producing [CuL(N3)]2 (3) or [CuL(NCS)] (4). On the other hand, upon reaction with Cu(ClO4)2·6H2O and Cu(NO3)2·3H2O, the Schiff base undergoes hydrolysis to yield ternary complexes [Cu(hap)(pn)(H2O)]ClO4 (5) and [Cu(hap)(pn)(H2O)]NO3 (6), respectively (Hhap = o-hydroxyacetophenone and pn = 1,3-pentanediamine). The ligand HL undergoes hydrolysis also on reaction with Ni(ClO4)2·6H2O or Ni(NO3)2·6H2O to yield [Ni(hap)2] (7). The structures of the complexes 2, 3, 5, 6, and 7 have been confirmed by single-crystal X-ray analysis. In complex 2, NiII possesses square-planar geometry, being coordinated by the tridentate mono-negative Schiff base, L and the isothiocyanate group. The coordination environment around CuII in complex 3 is very similar to that in complex 2 but here two units are joined together by end-on, axial-equatorial azide bridges to result in a dimer in which the geometry around CuII is square pyramidal. In both 5 and 6, the CuII atoms display the square-pyramidal environment; the equatorial sites being coordinated by the two amine groups of 1,3-pentanediamine and two oxygen atoms of o-hydroxyacetophenone. The axial site is coordinated by a water molecule. Complex 7 is a square-planar complex with the Ni atom bonded to four oxygen atoms from two hap moieties. The mononuclear units of 2 and dinuclear units of 3 are linked by strong hydrogen bonds to form a one-dimensional network. The mononuclear units of 5 and 6 are joined together to form a dimer by very strong hydrogen bonds through the coordinated water molecule. These dimers are further involved in hydrogen bonding with the respective counteranions to form 2-D net-like open frameworks.

Fractional order modeling of large three-dimensional RC networks

An incidence matrix analysis is used to model a three-dimensional network consisting of resistive and capacitive elements distributed across several interconnected layers. A systematic methodology for deriving a descriptor representation of the network with random allocation of the resistors and capacitors is proposed. Using a transformation of the descriptor representation into standard state-space form, amplitude and phase admittance responses of three-dimensional random RC networks are obtained. Such networks display an emergent behavior with a characteristic Jonscher-like response over a wide range of frequencies. A model approximation study of these networks is performed to infer the admittance response using integral and fractional order models. It was found that a fractional order model with only seven parameters can accurately describe the responses of networks composed of more than 70 nodes and 200 branches with 100 resistors and 100 capacitors. The proposed analysis can be used to model charge migration in amorphous materials, which may be associated to specific macroscopic or microscopic scale fractal geometrical structures in composites displaying a viscoelastic electromechanical response, as well as to model the collective responses of processes governed by random events described using statistical mechanics.

On the use of a time sequence of surface pressures in four-dimensional data assimilation

The possibility of using a time sequence of surface pressure observations in four-dimensional data assimilation is being investigated. It is shown that a linear multilevel quasi-geostrophic model can be updated successfully with surface data alone, provided the number of time levels are at least as many as the number of vertical levels. It is further demonstrated that current statistical analysis procedures are very inefficient to assimilate surface observations, and it is shown by numerical experiments that the vertical interpolation must be carried out using the structure of the most dominating baroclinic mode in order to obtain a satisfactory updating. Different possible ways towards finding a practical solution are being discussed.