The diabatic contour advective semi-Lagrangian model

This article describes a novel algorithmic development extending the contour advective semi-Lagrangian model to include nonconservative effects. The Lagrangian contour representation of finescale tracer fields, such as potential vorticity, allows for conservative, nondiffusive treatment of sharp gradients allowing very high numerical Reynolds numbers. It has been widely employed in accurate geostrophic turbulence and tracer advection simulations. In the present, diabatic version of the model the constraint of conservative dynamics is overcome by including a parallel Eulerian field that absorbs the nonconservative ( diabatic) tendencies. The diabatic buildup in this Eulerian field is limited through regular, controlled transfers of this field to the contour representation. This transfer is done with a fast newly developed contouring algorithm. This model has been implemented for several idealized geometries. In this paper a single-layer doubly periodic geometry is used to demonstrate the validity of the model. The present model converges faster than the analogous semi-Lagrangian models at increased resolutions. At the same nominal spatial resolution the new model is 40 times faster than the analogous semi-Lagrangian model. Results of an orographically forced idealized storm track show nontrivial dependency of storm-track statistics on resolution and on the numerical model employed. If this result is more generally applicable, this may have important consequences for future high-resolution climate modeling.

Linear viscoelasticity from molecular dynamics simulation of entangled polymers

The linear viscoelastic (LVE) spectrum is one of the primary fingerprints of polymer solutions and melts, carrying information about most relaxation processes in the system. Many single chain theories and models start with predicting the LVE spectrum to validate their assumptions. However, until now, no reliable linear stress relaxation data were available from simulations of multichain systems. In this work, we propose a new efficient way to calculate a wide variety of correlation functions and mean-square displacements during simulations without significant additional CPU cost. Using this method, we calculate stress−stress autocorrelation functions for a simple bead−spring model of polymer melt for a wide range of chain lengths, densities, temperatures, and chain stiffnesses. The obtained stress−stress autocorrelation functions were compared with the single chain slip−spring model in order to obtain entanglement related parameters, such as the plateau modulus or the molecular weight between entanglements. Then, the dependence of the plateau modulus on the packing length is discussed. We have also identified three different contributions to the stress relaxation:  bond length relaxation, colloidal and polymeric. Their dependence on the density and the temperature is demonstrated for short unentangled systems without inertia.

Hierarchical modeling of entangled polymers

We demonstrate that it is possible to link multi-chain molecular dynamics simulations with the tube model using a single chain slip-links model as a bridge. This hierarchical approach allows significant speed up of simulations, permitting us to span the time scales relevant for a comparison with the tube theory. Fitting the mean-square displacement of individual monomers in molecular dynamics simulations with the slip-spring model, we show that it is possible to predict the stress relaxation. Then, we analyze the stress relaxation from slip-spring simulations in the framework of the tube theory. In the absence of constraint release, we establish that the relaxation modulus can be decomposed as the sum of contributions from fast and longitudinal Rouse modes, and tube survival. Finally, we discuss some open questions regarding possible future directions that could be profitable in rendering the tube model quantitative, even for mildly entangled polymers

Conducting nanofibres produced by electrospinning

Electrospun fibres based on polypyrrole have been prepared. The incorporation of preformed polypyrrole into fibres electrospun from a carrier polymer can only be achieved when materials are prepared with particulates smaller than the cross-section of the fibre; even so there are some problems, with the substantial loss of material from the electrode tip. As an alternative approach, soluble polypyrroles can be prepared but these are not of sufficient viscosity to prepare electrospun fibres, once again a carrier polymer must be employed. More effective loadings are gained by the process of coating the outer surface of a pre-spun fibre; in this way electrospun fibres coated with polypyrrole can be prepared. This approach has been adapted to produce silver coated polymer fibres by the use of copolymers of styrene and 3-vinyl benzaldehyde.

Microscopic Mechanisms of Strain Hardening in Glassy Polymers

The mechanisms underlying the increase in stress for large mechanical strains of a polymer glass, quantified by the strain-hardening modulus, are still poorly understood. In the present paper we aim to elucidate this matter and present new mechanisms. Molecular-dynamics simulations of two polymers with very different strain-hardening moduli (polycarbonate and polystyrene) have been carried out. Nonaffine displacements occur because of steric hindrances and connectivity constraints. We argue that it is not necessary to introduce the concept of entanglements to understand strain hardening, but that hardening is rather coupled with the increase in the rate of nonaffine particle displacements. This rate increases faster for polycarbonate, which has the higher strain-hardening modulus. Also more nonaffine chain stretching is present for polycarbonate. It is shown that the inner distances of such a nonaffinely deformed chain can be well described by the inner distances of the worm-like chain, but with an effective stiffness length (equal to the Kuhn length for an infinite worm-like chain) that increases during deformation. It originates from the finite extensibility of the chain. In this way the increase in nonaffine particle displacement can be understood as resulting from an increase in the effective stiffness length of the perturbed chain during deformation, so that at larger strains a higher rate of plastic events in terms of nonaffine displacement is necessary, causing in turn the observed strain hardening in polymer glasses.

The timing of Quaternary calcrete development in semi-arid southeast Spain: Investigating the role of climate on calcrete genesis

Quaternary-aged calcrete horizons are common weathering products in arid and semi-arid regions. It is, however, unclear how calcrete forming processes respond to the major oscillations in climate that occur over the Quaternary period. This paper presents a U-series-based calcrete age database from the Sorbas basin, southeast Spain. The study constructs an age frequency distribution of these ages which is consequently compared to a range of palaeoenvironmental records from the Mediterranean. The age distribution presented here suggests that the formation of pedogenic calcrete horizons in the Sorbas basin primarily occurs during 'warm' isotope stages (MIS 1 and 5), with very few calcrete ages occurring during cold glacial/stadial stages (MIS 2, 3 and 4). It is suggested that this is a function of the environments that existed during 'warm' isotope stages being more conducive to calcrete development than those that existed during cold climate episodes. In a semi-arid region such as the Sorbas basin it is likely that increased aridity during glacial stages, coupled with reduced vegetation and accelerated landscape instability, was crucial in reducing rates of calcrete formation. In a semi-arid region such as southeast Spain, calcrete formation during the Quaternary, therefore, oscillates with climate change but is primarily a "warm" episode phenomenon. It is suggested that further studies are required to see how calcrete genesis responds to environmental change in more humid parts of the Mediterranean. (C) 2009 Elsevier B.V. All rights reserved.

A semi-quantitative approach to deriving a model structure for the uptake of organic chemicals by vegetation

An expert elicitation exercise was undertaken to determine those components and processes that are most important for modeling plant uptake of organic chemicals. The state of our knowledge of these processes was also assessed. This semi-quantitative analysis allowed the construction of an idealized model with seven compartments; soil bulk, soil water, roots, stem, leaves, fruit, and air. Three main areas were identified further research: 1) the uptake of organic chemicals by fruit; 2) the internal transfer of organic chemicals between plant structures (e.g., stem and leaves); and 3) the transfer via the soil-air-plant pathway. Until new data becomes available to quantify these processes, it is proposed that an equilibrium partitioning approach is used between plant components other than fruit or that models consist of both an edible and inedible compartment.

Effects of manure and fertilizer on grain yield, soil carbon and phosphorus in a 13-year field trial in semi-arid Kenya

Long-term indicators of soil fertility were assessed by measuring grain yield, soil organic carbon (SOC) and soil Olsen phosphorous for a P-deficient soil. In one set of treatments, goat manure was applied annually for 13 years at 0, 5 and 10 t ha(-1), and intercrops of sorghum/cowpea, millet/green gram and maize/pigeonpea were grown. Yield depended on rainfall and trends with time were not identifiable. Manure caused an upward trend in SOC, but 10 t ha(-1) manure did not give significantly more SOC than 5 t ha(-1). Only 10 t ha(-1) manure increased Olsen P. Measurements of both SOC and Olsen P are recommended. In another set of treatments, manure was applied for four years; the residual effect lasted another seven to eight years when assessed by yield, SOC and Olsen P Treatment with mineral fertilizers provided the same rates of N and P as 5 t hat manure and yields from manure and fertilizer were similar. Fertilizer increased Olsen P but not SOC. Management systems with occasional manure application and intermediate fertilizer applications should be assessed. Inputs and offtakes of C, N and P were measured for three years. Approximately 16, 25 and 11% of C, N and P respectively were stabilized into soil organic matter from 5 t ha(-1) a(-1) manure. The majority of organic P was fixed as soil inorganic P.

Effect of manure application on crop yield and soil chemical properties in a long-term field trial of semi-arid Kenya

The sustainability of cereal/legume intercropping was assessed by monitoring trends in grain yield, soil organic C (SOC) and soil extractable P (Olsen method) measured over 13 years at a long-term field trial on a P-deficient soil in semi-arid Kenya. Goat manure was applied annually for 13 years at 0, 5 and 10 t ha(-1) and trends in grain yield were not identifiable because of season-to-season variations. SOC and Olsen P increased for the first seven years of manure application and then remained constant. The residual effect of manure applied for four years only lasted another seven to eight years when assessed by yield, SOC and Olsen P. Mineral fertilizers provided the same annual rates of N and P as in 5 t ha(-1) manure and initially ,gave the same yield as manure, declining after nine years to about 80%. Therefore, manure applications could be made intermittently and nutrient requirements topped-up with fertilizers. Grain yields for sorghum with continuous manure were described well by correlations with rainfall and manure input only, if data were excluded for seasons with over 500 mm rainfall. A comprehensive simulation model should correctly describe crop losses caused by excess water.