An isentropic flow algorithm with body fitted meshes

An efficient algorithm based on flux difference splitting is presented for the solution of the three-dimensional equations of isentropic flow in a generalised coordinate system, and with a general convex gas law. The scheme is based on solving linearised Riemann problems approximately and in more than one dimension incorporates operator splitting. The algorithm requires only one function evaluation of the gas law in each computational cell. The scheme has good shock capturing properties and the advantage of using body-fitted meshes. Numerical results are shown for Mach 3 flow of air past a circular cylinder. Furthermore, the algorithm also applies to shallow water flows by employing the familiar gas dynamics analogy.

Prediction of supercritical flow in open channels

A finite difference scheme based on flux difference splitting is presented for the solution of the one-dimensional shallow water equations in open channels. A linearised problem, analogous to that of Riemann for gas dynamics, is defined and a scheme, based on numerical characteristic decomposition, is presented for obtaining approximate solutions to the linearised problem. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second order scheme which avoids non-physical, spurious oscillations. The scheme is applied to a problem of flow in a river whose geometry induces a region of supercritical flow.

Second order accurate upwind difference schemes for scalar conservation laws with source terms

A second order accurate, characteristic-based, finite difference scheme is developed for scalar conservation laws with source terms. The scheme is an extension of well-known second order scalar schemes for homogeneous conservation laws. Such schemes have proved immensely powerful when applied to homogeneous systems of conservation laws using flux-difference splitting. Many application areas, however, involve inhomogeneous systems of conservation laws with source terms, and the scheme presented here is applied to such systems in a subsequent paper.

A numerical scheme for two-dimensional, open channel flows with non-rectangular geometries

An algorithm based on flux difference splitting is presented for the solution of two-dimensional, open channel flows. A transformation maps a non-rectangular, physical domain into a rectangular one. The governing equations are then the shallow water equations, including terms of slope and friction, in a generalized coordinate system. A regular mesh on a rectangular computational domain can then be employed. The resulting scheme has good jump capturing properties and the advantage of using boundary/body-fitted meshes. The scheme is applied to a problem of flow in a river whose geometry induces a region of supercritical flow.

Ligand dynamics in myoglobin: calculation of infrared spectra for photodissociated NO

We present molecular dynamics simulations of the photodissociated state of MbNO performed at 300 K using a fluctuating charge model for the nitric oxide (NO) ligand. After dissociation, NO is observed to remain mainly in the centre of the distal haem pocket, although some movement towards the primary docking site and the xenon-4 pocket can be seen. We calculate the NO infrared spectrum for the photodissociated ligand within the haem pocket and find a narrow peak in the range 1915-1922 cm(-1). The resulting blue shift of 1 to 8 cm(-1) compared to gas-phase NO is much smaller than the red shifts calculated and observed for carbon monoxide (CO) in Mb. A small splitting, due to NO in the xenon-4 pocket, is also observed. At lower temperatures, the spectra and conformational space explored by the ligand remain largely unchanged, but the electrostatic interactions with residue His64 become increasingly significant in determining the details of the ligand orientation within the distal haem pocket. The investigation of the effect of the L29F mutation reveals significant differences between the behaviour of NO and that of CO, and suggests a coupling between the ligand and the protein dynamics due to the different ligand dipole moments.

Theoretical investigation of infrared spectra and pocket dynamics of photodissociated carbonmonoxy myoglobin

Molecular dynamics simulations of the photodissociated state of carbonmonoxy myoglobin (MbCO) are presented using a fluctuating charge model for CO. A new three-point charge model is fitted to high-level ab initio calculations of the dipole and quadrupole moment functions taken from the literature. The infrared spectrum of the CO molecule in the heme pocket is calculated using the dipole moment time autocorrelation function and shows good agreement with experiment. In particular, the new model reproduces the experimentally observed splitting of the CO absorption spectrum. The splitting of 3–7 cm−1 (compared to the experimental value of 10 cm−1) can be directly attributed to the two possible orientations of CO within the docking site at the edge of the distal heme pocket (the B states), as previously suggested on the basis of experimental femtosecond time-resolved infrared studies. Further information on the time evolution of the position and orientation of the CO molecule is obtained and analyzed. The calculated difference in the free energy between the two possible orientations (Fe···CO and Fe···OC) is 0.3 kcal mol−1 and agrees well with the experimentally estimated value of 0.29 kcal mol−1. A comparison of the new fluctuating charge model with an established fixed charge model reveals some differences that may be critical for the correct prediction of the infrared spectrum and energy barriers. The photodissociation of CO from the myoglobin mutant L29F using the new model shows rapid escape of CO from the distal heme pocket, in good agreement with recent experimental data. The effect of the protein environment on the multipole moments of the CO ligand is investigated and taken into account in a refined model. Molecular dynamics simulations with this refined model are in agreement with the calculations based on the gas-phase model. However, it is demonstrated that even small changes in the electrostatics of CO alter the details of the dynamics.

Characterising the variability and extremes of the stratospheric polar vortices using 2D moment analysis

The mean state, variability and extreme variability of the stratospheric polar vortices, with an emphasis on the Northern Hemisphere vortex, are examined using 2-dimensional moment analysis and Extreme Value Theory (EVT). The use of moments as an analysis to ol gives rise to information about the vortex area, centroid latitude, aspect ratio and kurtosis. The application of EVT to these moment derived quantaties allows the extreme variability of the vortex to be assessed. The data used for this study is ECMWF ERA-40 potential vorticity fields on interpolated isentropic surfaces that range from 450K-1450K. Analyses show that the most extreme vortex variability occurs most commonly in late January and early February, consistent with when most planetary wave driving from the troposphere is observed. Composites around sudden stratospheric warming (SSW) events reveal that the moment diagnostics evolve in statistically different ways between vortex splitting events and vortex displacement events, in contrast to the traditional diagnostics. Histograms of the vortex diagnostics on the 850K (∼10hPa) surface over the 1958-2001 period are fitted with parametric distributions, and show that SSW events comprise the majority of data in the tails of the distributions. The distribution of each diagnostic is computed on various surfaces throughout the depth of the stratosphere, and shows that in general the vortex becomes more circular with higher filamentation at the upper levels. The Northern Hemisphere (NH) and Southern Hemisphere (SH) vortices are also compared through the analysis of their respective vortex diagnostics, and confirm that the SH vortex is less variable and lacks extreme events compared to the NH vortex. Finally extreme value theory is used to statistically mo del the vortex diagnostics and make inferences about the underlying dynamics of the polar vortices.

Synthesis, structure and magnetic properties of mono- and di-nuclear nickel(II) thiocyanate complexes with tridentate N3 donor Schiff bases

The 1:1 condensation of N-methyl-1,3-diaminopropane and N,N-diethyl-1,2-diminoethane with 2-acetylpyridine, respectively at high dilution gives the tridentate mono-condensed Schiff bases N-methyl-N'-(1-pyridin-2-yl-ethylidene)-propane-1,3-diamine (L-1) and N,N-diethyl-N'-(1-pyridin-2-yl-ethylidene)-ethane-1,2-diamine (L-2). The tridentate ligands were allowed to react with methanol solutions of nickel(II) thiocyanate to prepare the complexes [Ni(L-1)(SCN)(2)(OH2) (1) and [{Ni(L-2)(SCN)}(2)] (2). Single crystal X-ray diffraction was used to confirm the structures of the complexes. The nickel(II) in complex 1 is bonded to three nitrogen donor atoms of the ligand L-1 in a mer orientation, together with two thiocyanates bonded through nitrogen and a water molecule, and it is the first Schiff base complex of nickel(II) containing both thiocyanate and coordinated water. The coordinated water initiates a hydrogen bonded 2D network. In complex 2, the nickel ion occupies a slightly distorted octahedral coordination sphere, being bonded to three nitrogen atoms from the ligand L-2, also in a mer orientation, and two thiocyanate anions through nitrogen. In contrast to 1, the sixth coordination site is occupied by a sulfur atom from a thiocyanate anion in an adjacent molecule, thus creating a centrosymmetric dimer. A variable temperature magnetic study of complex 2 indicates the simultaneous presence of zero-field splitting, weak intramolecular ferromagnetic coupling and intermolecular antiferromagnetic interactions between the nickel(II) centers.

The power of NMR.: in two and three dimensions

Soon after its discovery in the 1950s, NMR had become an indispensable tool fr chemists. In the 1970s and 1980s, the power of the technique was extended from one dimension to two and even three dimensions, opening up exciting applkications in both chemistry and biochemistry. the success of one dimensional. high-resolution NMR stems from the unique insights that it can provide about molecular structure. The chemical shift of a nucleus gives invaluable information abut the chemical environment in which that nucleus is located, Coupling interactions between hydorgen nuclei, as revealed by characteristic splitting patterns inthe 1H-NMR spectrum, provide informaton about the loaction of one group of hydorgen atoms relative to others inthe molecule. And the nuclearf Overhauser effect (nOe) can shed light on molecular stereochemistry.

Land re-use, complexity and actor-networks: a framework for research

This paper will present a conceptual framework for the examination of land redevelopment based on a complex systems/networks approach. As Alvin Toffler insightfully noted, modern scientific enquiry has become exceptionally good at splitting problems into pieces but has forgotten how to put the pieces back together. Twenty-five years after his remarks, governments and corporations faced with the requirements of sustainability are struggling to promote an ‘integrated’ or ‘holistic’ approach to tackling problems. Despite the talk, both practice and research provide few platforms that allow for ‘joined up’ thinking and action. With socio-economic phenomena, such as land redevelopment, promising prospects open up when we assume that their constituents can make up complex systems whose emergent properties are more than the sum of the parts and whose behaviour is inherently difficult to predict. A review of previous research shows that it has mainly focused on idealised, ‘mechanical’ views of property development processes that fail to recognise in full the relationships between actors, the structures created and their emergent qualities. When reality failed to live up to the expectations of these theoretical constructs then somebody had to be blamed for it: planners, developers, politicians. However, from a ‘synthetic’ point of view the agents and networks involved in property development can be seen as constituents of structures that perform complex processes. These structures interact, forming new more complex structures and networks. Redevelopment then can be conceptualised as a process of transformation: a complex system, a ‘dissipative’ structure involving developers, planners, landowners, state agencies etc., unlocks the potential of previously used sites, transforms space towards a higher order of complexity and ‘consumes’ but also ‘creates’ different forms of capital in the process. Analysis of network relations point toward the ‘dualism’ of structure and agency in these processes of system transformation and change. Insights from actor network theory can be conjoined with notions of complexity and chaos to build an understanding of the ways in which actors actively seek to shape these structures and systems, whilst at the same time are recursively shaped by them in their strategies and actions. This approach transcends the blame game and allows for inter-disciplinary inputs to be placed within a broader explanatory framework that does away with many past dichotomies. Better understanding of the interactions between actors and the emergent qualities of the networks they form can improve our comprehension of the complex socio-spatial phenomena that redevelopment comprises. The insights that this framework provides when applied in UK institutional investment into redevelopment are considered to be significant.

Diamond mining, rice farming, and a ‘Maggi cube’: a viable survivial strategy in rural Liberia?

Since the conclusion of its 14-year civil war in 2003, Liberia has struggled economically. Jobs are in short supply and operational infrastructural services, such as electricity and running water, are virtually nonexistent. The situation has proved especially challenging for the scores of people who fled the country in the 1990s to escape the violence and who have since returned to re-enter their lives. With few economic prospects on hand, many have elected to enter the artisanal diamond mining sector, which has earned notoriety for perpetuating the country's civil war. This article critically reflects on the fate of these Liberians, many of whom, because of a lack of government support, finances, manpower and technological resources, have forged deals with hired labourers to work artisanal diamond fields. Specifically, in exchange for meals containing locally grown rice and a Maggi (soup) cube, hired hands mine diamondiferous territories, splitting the revenues accrued from the sales of recovered stones amongst themselves and the individual ‘claimholder’ who hired them. Although this cycle—referred to here as ‘diamond mining, rice farming and a Maggi cube’—helps to buffer against poverty, few of the parties involved will ever progress beyond a subsistence level

New results on the thermodynamical properties of the climate system

In this paper the authors exploit two equivalent formulations of the average rate of material entropy production in the climate system to propose an approximate splitting between contributions due to vertical and eminently horizontal processes. This approach is based only on 2D radiative fields at the surface and at the top of atmosphere. Using 2D fields at the top of atmosphere alone, lower bounds to the rate of material entropy production and to the intensity of the Lorenz energy cycle are derived. By introducing a measure of the efficiency of the planetary system with respect to horizontal thermodynamic processes, it is possible to gain insight into a previous intuition on the possibility of defining a baroclinic heat engine extracting work from the meridional heat flux. The approximate formula of the material entropy production is verified and used for studying the global thermodynamic properties of climate models (CMs) included in the Program for Climate Model Diagnosis and Intercomparison (PCMDI)/phase 3 of the Coupled Model Intercomparison Project (CMIP3) dataset in preindustrial climate conditions. It is found that about 90% of the material entropy production is due to vertical processes such as convection, whereas the large-scale meridional heat transport contributes to only about 10% of the total. This suggests that the traditional two-box models used for providing a minimal representation of entropy production in planetary systems are not appropriate, whereas a basic—but conceptually correct—description can be framed in terms of a four-box model. The total material entropy production is typically 55 mW m−2 K−1, with discrepancies on the order of 5%, and CMs’ baroclinic efficiencies are clustered around 0.055. The lower bounds on the intensity of the Lorenz energy cycle featured by CMs are found to be around 1.0–1.5 W m−2, which implies that the derived inequality is rather stringent. When looking at the variability and covariability of the considered thermodynamic quantities, the agreement among CMs is worse, suggesting that the description of feedbacks is more uncertain. The contributions to material entropy production from vertical and horizontal processes are positively correlated, so that no compensation mechanism seems in place. Quite consistently among CMs, the variability of the efficiency of the system is a better proxy for variability of the entropy production due to horizontal processes than that of the large-scale heat flux. The possibility of providing constraints on the 3D dynamics of the fluid envelope based only on 2D observations of radiative fluxes seems promising for the observational study of planets and for testing numerical models.

The effect of climate change on the variability of the Northern Hemisphere stratospheric polar vortex

With extreme variability of the Arctic polar vortex being a key link for stratosphere–troposphere influences, its evolution into the twenty-first century is important for projections of changing surface climate in response to greenhouse gases. Variability of the stratospheric vortex is examined using a state-of-the-art climate model and a suite of specifically developed vortex diagnostics. The model has a fully coupled ocean and a fully resolved stratosphere. Analysis of the standard stratospheric zonal mean wind diagnostic shows no significant increase over the twenty-first century in the number of major sudden stratospheric warmings (SSWs) from its historical value of 0.7 events per decade, although the monthly distribution of SSWs does vary, with events becoming more evenly dispersed throughout the winter. However, further analyses using geometric-based vortex diagnostics show that the vortex mean state becomes weaker, and the vortex centroid is climatologically more equatorward by up to 2.5°, especially during early winter. The results using these diagnostics not only characterize the vortex structure and evolution but also emphasize the need for vortex-centric diagnostics over zonally averaged measures. Finally, vortex variability is subdivided into wave-1 (displaced) and -2 (split) components, and it is implied that vortex displacement events increase in frequency under climate change, whereas little change is observed in splitting events.

An iterative contractive framework for probe methods: LASSO

We present a new iterative approach called Line Adaptation for the Singular Sources Objective (LASSO) to object or shape reconstruction based on the singular sources method (or probe method) for the reconstruction of scatterers from the far-field pattern of scattered acoustic or electromagnetic waves. The scheme is based on the construction of an indicator function given by the scattered field for incident point sources in its source point from the given far-field patterns for plane waves. The indicator function is then used to drive the contraction of a surface which surrounds the unknown scatterers. A stopping criterion for those parts of the surfaces that touch the unknown scatterers is formulated. A splitting approach for the contracting surfaces is formulated, such that scatterers consisting of several separate components can be reconstructed. Convergence of the scheme is shown, and its feasibility is demonstrated using a numerical study with several examples.

Ligand substituents effect on 10Dq

Reaction of 5,6-dihydro-5,6-epoxy-1,10-phenanthroline (L) with Ni(ClO4)(2)center dot 6H(2)O in methanol in 3:1 M proportion at room temperature yields [NiL3](ClO4)(2)center dot 2H(2)O. The X-ray crystal structure of the cation Nil(3)(2+) has been determined. Aminolysis of the three epoxide rings in NiL32+ by 4-substituted anilines in boiling water without any Lewis acid catalyst gives a family of Ni(II) complexes with octahedral NiL62+ core. In these complexes, crystal field splitting 10Dq varies from 11601 to 15798 cm(-1) in acetonitrile. The variation in 10Dq is found to be satisfactorily linear (r(2) = 0.951) with the Hammett sigma(R) parameter of the substituent on the anilino fragment. 10Dq increases with the increase in the electron donation ability of the substituent.