Inclusion of nonlinear demand-supply relationships within large-scale partial equilibrium linear programming models

This paper presents a new method for the inclusion of nonlinear demand and supply relationships within a linear programming model. An existing method for this purpose is described first and its shortcomings are pointed out before showing how the new approach overcomes those difficulties and how it provides a more accurate and 'smooth' (rather than a kinked) approximation of the nonlinear functions as well as dealing with equilibrium under perfect competition instead of handling just the monopolistic situation. The workings of the proposed method are illustrated by extending a previously available sectoral model for the UK agriculture.

The potential for the use of agent-based models in ecotoxicology

This chapter introduces ABMs, their construction, and the pros and cons of their use. Although relatively new, agent-basedmodels (ABMs) have great potential for use in ecotoxicological research – their primary advantage being the realistic simulations that can be constructed and particularly their explicit handling of space and time in simulations. Examples are provided of their use in ecotoxicology primarily exemplified by different implementations of the ALMaSS system. These examples presented demonstrate how multiple stressors, landscape structure, details regarding toxicology, animal behavior, and socioeconomic effects can and should be taken into account when constructing simulations for risk assessment. Like ecological systems, in ABMs the behavior at the system level is not simply the mean of the component responses, but the sum of the often nonlinear interactions between components in the system; hence this modeling approach opens the door to implementing and testing much more realistic and holistic ecotoxicological models than are currently used.

Use of a novel Hill-climbing genetic algorithm in protein folding simulations

We have developed a novel Hill-climbing genetic algorithm (GA) for simulation of protein folding. The program (written in C) builds a set of Cartesian points to represent an unfolded polypeptide's backbone. The dihedral angles determining the chain's configuration are stored in an array of chromosome structures that is copied and then mutated. The fitness of the mutated chain's configuration is determined by its radius of gyration. A four-helix bundle was used to optimise simulation conditions, and the program was compared with other, larger, genetic algorithms on a variety of structures. The program ran 50% faster than other GA programs. Overall, tests on 100 non-redundant structures gave comparable results to other genetic algorithms, with the Hill-climbing program running from between 20 and 50% faster. Examples including crambin, cytochrome c, cytochrome B and hemerythrin gave good secondary structure fits with overall alpha carbon atom rms deviations of between 5 and 5.6 Angstrom with an optimised hydrophobic term in the fitness function. (C) 2003 Elsevier Ltd. All rights reserved.

Simulating the formation of secondary organic aerosol from the photooxidation of aromatic hydrocarbons

The formation and composition of secondary organic aerosol (SOA) from the photooxidation of benzene, p-xylene, and 1,3,5-trimethylbenzene has been simulated using the Master Chemical Mechanism version 3.1 (MCM v3.1) coupled to a representation of the transfer of organic material from the gas to particle phase. The combined mechanism was tested against data obtained from a series of experiments conducted at the European Photoreactor (EUPHORE) outdoor smog chamber in Valencia, Spain. Simulated aerosol mass concentrations compared reasonably well with the measured SOA data only after absorptive partitioning coefficients were increased by a factor of between 5 and 30. The requirement of such scaling was interpreted in terms of the occurrence of unaccounted-for association reactions in the condensed organic phase leading to the production of relatively more nonvolatile species. Comparisons were made between the relative aerosol forming efficiencies of benzene, toluene, p-xylene, and 1,3,5-trimethylbenzene, and differences in the OH-initiated degradation mechanisms of these aromatic hydrocarbons. A strong, nonlinear relationship was observed between measured (reference) yields of SOA and (proportional) yields of unsaturated dicarbonyl aldehyde species resulting from ring-fragmenting pathways. This observation, and the results of the simulations, is strongly suggestive of the involvement of reactive aldehyde species in association reactions occurring in the aerosol phase, thus promoting SOA formation and growth. The effect of NO, concentrations on SOA formation efficiencies (and formation mechanisms) is discussed.

Function of the hydration layer around an antifreeze protein revealed by atomistic molecular dynamics simulations

Atomistic molecular dynamics simulations are used to investigate the mechanism by which the antifreeze protein from the spruce budworm, Choristoneura fumiferana, binds to ice. Comparison of structural and dynamic properties of the water around the three faces of the triangular prism-shaped protein in aqueous solution reveals that at low temperature the water structure is ordered and the dynamics slowed down around the ice-binding face of the protein, with a disordering effect observed around the other two faces. These results suggest a dual role for the solvation water around the protein. The preconfigured solvation shell around the ice-binding face is involved in the initial recognition and binding of the antifreeze protein to ice by lowering the barrier for binding and consolidation of the protein:ice interaction surface. Thus, the antifreeze protein can bind to the molecularly rough ice surface by becoming actively involved in the formation of its own binding site. Also, the disruption of water structure around the rest of the protein helps prevent the adsorbed protein becoming covered by further ice growth.

Building cluster and shading in urban canyon for hot dry climate Part 2: Shading simulations

Under low latitude conditions, minimisation of solar irradiance within the urban environment may often be an important criterion in urban design. This can be achieved when the obstruction angle is large (high H/W ratio, H = height, W = width). Solar access to streets can always be decreased by increasing H/W to larger values. It is shown in this paper that the street canyon orientation (and not only the H/W ratio) has a considerable effect on solar shading and urban microclimate. The paper demonstrates through a series of shading simulation and temperature measurements that a number of useful relationships can be developed between the geometry and the microclimate of urban street canyons. These relationships are potentially helpful to assist in the formulation of urban design guidelines governing street dimensions and orientations for use by urban designers.

Modelling oil absorption during post-frying cooling - II: solution of the mathematical model, model testing and simulations

The mathematical models that describe the immersion-frying period and the post-frying cooling period of an infinite slab or an infinite cylinder were solved and tested. Results were successfully compared with those found in the literature or obtained experimentally, and were discussed in terms of the hypotheses and simplifications made. The models were used as the basis of a sensitivity analysis. Simulations showed that a decrease in slab thickness and core heat capacity resulted in faster crust development. On the other hand, an increase in oil temperature and boiling heat transfer coefficient between the oil and the surface of the food accelerated crust formation. The model for oil absorption during cooling was analysed using the tested post-frying cooling equation to determine the moment in which a positive pressure driving force, allowing oil suction within the pore, originated. It was found that as crust layer thickness, pore radius and ambient temperature decreased so did the time needed to start the absorption. On the other hand, as the effective convective heat transfer coefficient between the air and the surface of the slab increased the required cooling time decreased. In addition, it was found that the time needed to allow oil absorption during cooling was extremely sensitive to pore radius, indicating the importance of an accurate pore size determination in future studies.

The eMinerals project: developing the concept of the virtual organisation to support collaborative work on molecular-scale environmental simulations

The eMinerals project has established an integrated compute and data minigrid infrastructure together with a set of collaborative tools,. The infrastructure is designed to support molecular simulation scientists working together as a virtual organisation aiming to understand a number of strategic processes in environmental science. The eMinerals virtual organisation is now working towards applying this infrastructure to tackle a new generation of scientific problems. This paper describes the achievements of the eMinerals virtual organisation to date, and describes ongoing applications of the virtual organisation infrastructure.

A tunable radial basis function model for nonlinear system identification using particle swarm optimisation

A tunable radial basis function (RBF) network model is proposed for nonlinear system identification using particle swarm optimisation (PSO). At each stage of orthogonal forward regression (OFR) model construction, PSO optimises one RBF unit's centre vector and diagonal covariance matrix by minimising the leave-one-out (LOO) mean square error (MSE). This PSO aided OFR automatically determines how many tunable RBF nodes are sufficient for modelling. Compared with the-state-of-the-art local regularisation assisted orthogonal least squares algorithm based on the LOO MSE criterion for constructing fixed-node RBF network models, the PSO tuned RBF model construction produces more parsimonious RBF models with better generalisation performance and is computationally more efficient.