Dynamic security assessment of power systems using structure-preserving energy functions

An application of direct methods to dynamic security assessment of power systems using structure-preserving energy functions (SPEF) is presented. The transient energy margin (TEM) is used as an index for checking the stability of the system as well as ranking the contigencies based on their severity. The computation of the TEM requires the evaluation of the critical energy and the energy at fault clearing. Usually this is done by simulating the faulted trajectory, which is time-consuming. In this paper, a new algorithm which eliminates the faulted trajectory estimation is presented to calculate the TEM. The system equations and the SPEF are developed using the centre-of-inertia (COI) formulation and the loads are modelled as arbitrary functions of the respective bus voltages. The critical energy is evaluated using the potential energy boundary surface (PEBS) method. The method is illustrated by considering two realistic power system examples.

Empirical torsional potential functions from protein structure data. Phi- and psi-potentials for non-glycyl amino acid residues.

The torsional potential functions Vt(phi) and Vt(psi) around single bonds N--C alpha and C alpha--C, which can be used in conformational studies of oligopeptides, polypeptides and proteins, have been derived, using crystal structure data of 22 globular proteins, fitting the observed distribution in the (phi, psi)-plane with the value of Vtot(phi, psi), using the Boltzmann distribution. The averaged torsional potential functions, obtained from various amino acid residues in L-configuration, are Vt(phi) = 1.0 cos (phi + 60 degrees); Vt(psi) = 0.5 cos (psi + 60 degrees) - 1.0 cos (2 psi + 30 degrees) - 0.5 cos (3 psi + 30 degrees). The dipeptide energy maps Vtot(phi, psi) obtained using these functions, instead of the normally accepted torsional functions, were found to explain various observations, such as the absence of the left-handed alpha helix and the C7 conformation, and the relatively high density of points near the line psi = 0 degrees. These functions derived from observational data on protein structures, will, it is hoped, explain various previously unexplained facts in polypeptide conformation.

Local systems, global impacts : Using life cycle assessment to analyse the potential and constraints of industrial symbioses

Human activities extract and displace different substances and materials from the earth s crust, thus causing various environmental problems, such as climate change, acidification and eutrophication. As problems have become more complicated, more holistic measures that consider the origins and sources of pollutants have been called for. Industrial ecology is a field of science that forms a comprehensive framework for studying the interactions between the modern technological society and the environment. Industrial ecology considers humans and their technologies to be part of the natural environment, not separate from it. Industrial operations form natural systems that must also function as such within the constraints set by the biosphere. Industrial symbiosis (IS) is a central concept of industrial ecology. Industrial symbiosis studies look at the physical flows of materials and energy in local industrial systems. In an ideal IS, waste material and energy are exchanged by the actors of the system, thereby reducing the consumption of virgin material and energy inputs and the generation of waste and emissions. Companies are seen as part of the chains of suppliers and consumers that resemble those of natural ecosystems. The aim of this study was to analyse the environmental performance of an industrial symbiosis based on pulp and paper production, taking into account life cycle impacts as well. Life Cycle Assessment (LCA) is a tool for quantitatively and systematically evaluating the environmental aspects of a product, technology or service throughout its whole life cycle. Moreover, the Natural Step Sustainability Principles formed a conceptual framework for assessing the environmental performance of the case study symbiosis (Paper I). The environmental performance of the case study symbiosis was compared to four counterfactual reference scenarios in which the actors of the symbiosis operated on their own. The research methods used were process-based life cycle assessment (LCA) (Papers II and III) and hybrid LCA, which combines both process and input-output LCA (Paper IV). The results showed that the environmental impacts caused by the extraction and processing of the materials and the energy used by the symbiosis were considerable. If only the direct emissions and resource use of the symbiosis had been considered, less than half of the total environmental impacts of the system would have been taken into account. When the results were compared with the counterfactual reference scenarios, the net environmental impacts of the symbiosis were smaller than those of the reference scenarios. The reduction in environmental impacts was mainly due to changes in the way energy was produced. However, the results are sensitive to the way the reference scenarios are defined. LCA is a useful tool for assessing the overall environmental performance of industrial symbioses. It is recommended that in addition to the direct effects, the upstream impacts should be taken into account as well when assessing the environmental performance of industrial symbioses. Industrial symbiosis should be seen as part of the process of improving the environmental performance of a system. In some cases, it may be more efficient, from an environmental point of view, to focus on supply chain management instead.  

Empirical torsional potential functions from protiesn structure DATA Phi- and Psi-Potentials for Non-glycyl Amino Acid Residues

The torsional potential functions Vt(φ) and Vt(ψ) around single bonds N–Cα and Cα-C, which can be used in conformational studies of oligopeptides, polypeptides and proteins, have been derived, using crystal structure data of 22 globular proteins, fitting the observed distribution in the (φ, ψ)-plane with the value of Vtot(φ, ψ), using the Boltzmann distribution. The averaged torsional potential functions, obtained from various amino acid residues in l-configuration, are Vt(φ) = – 1.0 cos (φ + 60°); Vt(ψ) = – 0.5 cos (ψ + 60°) – 1.0 cos (2ψ + 30°) – 0.5 cos (3ψ + 30°). The dipeptide energy maps Vtot(φ, ψ) obtained using these functions, instead of the normally accepted torsional functions, were found to explain various observations, such as the absence of the left-handed alpha helix and the C7 conformation, and the relatively high density of points near the line ψ = 0°. These functions, derived from observational data on protein structures, will, it is hoped, explain various previously unexplained facts in polypeptide conformation.

Harvesting undelimbed Scots pine (Pinus sylvestris L.) from first thinnings for integrated production of kraft pulp and energy

The present study evaluates the feasibility of undelimbed Scots pine (Pinus sylvestris L.) for integrated production of pulp and energy in a kraft pulp mill from the technical, economic and environmental points of view, focusing on the potential of bundle harvesting. The feasibility of tree sections for pulp production was tested by conducting an industrial wood-handling experiment, laboratory cooking and bleaching trials, using conventional small-diameter Scots pine pulpwood as a reference. These trials showed that undelimbed Scots pine sections can be processed in favourable conditions as a blend with conventional small-diameter pulpwood without reducing the pulp quality. However, fibre losses at various phases of the process may increase when using undelimbed material. In the economic evaluation, both pulp production and wood procurement costs were considered, using the relative wood paying capability of a kraft pulp mill as a determinant. The calculations were made for three Scots pine first-thinning stands with the breast-height diameter of the removal (6 12 cm) as the main distinctive factor. The supply chains included in the comparison were based on cut-to-length harvesting, whole-tree harvesting and bundle harvesting (whole-tree bundling). With the current ratio of pulp and energy prices, the wood paying capability declines with an increase in the proportion of the energy fraction of the raw material. The supply system based on the cut-to-length method was the most efficient option, resulting in the highest residual value at stump in most cases. A decline in the pulp price and an increase in the energy price improved the competitiveness of the whole-tree systems. With short truck transportation distances and low pulp prices, however, the harvesting of loose whole trees can result in higher residual value at stump in small-diameter stands. While savings in transportation costs did not compensate for the high cutting and compaction costs by the second prototype of the bundle harvester, an increase in transportation distances improved its competitiveness. Since harvesting undelimbed assortments increases nutrient export from the site, which can affect soil productivity, the whole-tree alternatives included in the present study cannot be recommended on infertile peatlands and mineral soils. The harvesting of loose whole trees or bundled whole trees implies a reduction in protective logging residues and an increase in site traffic or payloads. These factors increase the risk of soil damage, especially on peat soils with poor bearing capacity. Within the wood procurement parameters which were examined, the CO2 emissions of the supply systems varied from 13 27 kg m3. Compaction of whole trees into bundles reduced emissions from transportation by 30 39%, but these reductions were insufficient to compensate for the increased emissions from cutting and compaction.

Temperature dependent free energy surface of polymer folding from equilibrium and quench studies

Langevin dynamics simulation studies have been employed to calculate the temperature dependent free energy surface and folding characteristics of a 500 monomer long linear alkane (polyethylene) chain with a realistic interaction potential. Both equilibrium and temperature quench simulation studies have been carried out. Using the shape anisotropy parameter (S) of the folded molecule as the order parameter, we find a weakly first order phase transition between the high-temperature molten globule and low-temperature rodlike crystalline states separated by a small barrier of the order of k(B)T. Near the melting temperature (580 K), we observe an intriguing intermittent fluctuation with pronounced ``1/f noise characteristics'' between these two states with large difference in shape and structure. We have also studied the possibilities of different pathways of folding to states much below the melting point. At 300 K starting from the all-trans linear configuration, the chain folds stepwise into a very regular fourfold crystallite with very high shape anisotropy. Whereas, when quenched from a high temperature (900 K) random coil regime, we identify a two step transition from the random coiled state to a molten globulelike state and, further, to a anisotropic rodlike state. The trajectory reveals an interesting coupling between the two order parameters, namely, radius of gyration (R-g) and the shape anisotropy parameter (S). The rodlike final state of the quench trajectory is characterized by lower shape anisotropy parameter and significantly larger number of gauche defects as compared to the final state obtained through equilibrium simulation starting from all-trans linear chain. The quench study shows indication of a nucleationlike pathway from the molten globule to the rodlike state involving an underlying rugged energy landscape. (C) 2010 American Institute of Physics. doi:10.1063/1.3509398]

Chemical potential of oxygen for iron-rutile-ilmenite and iron-ilmenite-ulvospinel equilibria

The chemical potential of oxygen corresponding to the iron-rutile-ilmenite (IRI) and iron-ilmenite-ulvospinel (IIU) equilibria has been measured employing solid-state galvanic cells,$$Pt, Fe + TiO_2 + FeTiO_3 //(Y_2 O_3 ) ZrO_2 //Fe + FeO, Pt$$ and $${\text{Pt, Fe + FeTiO}}_{\text{3}} {\text{ + Fe}}_{\text{2}} {\text{TiO}}_{\text{4}} {\text{//(Y}}_{\text{2}} {\text{0}}_{\text{3}} {\text{) ZrO}}_{\text{2}} {\text{//Fe + FeO, Pt}}$$ in the temperature range of 875 to 1275 K and 900 to 1373 K, respectively. The cells are written such that the right-hand electrodes are positive. The electromotive force (emf) of both the cells was found to be reversible and to vary linearly with temperature over the entire range of measurement. The chemical potential of oxygen for IRI equilibrium is represented by Δμo2(IRI) = -550,724 - 29.445T + 20.374T InT(±210) J mol−1 (875 <-T<- 1184 K) = -620,260 + 369.593T - 27.716T lnT(±210) J mol−1 (1184 <-T<- 1275 K) and that for IIU equilibrium by Δμo2(IIU) = -501,800 - 49.035T + 20.374T lnT(±210) J mol−1 (900 <-T<- 1184 K) = -571,336 + 350.003T− 27.716T lnT(=−210) J mol-1 (1184 <-T<- 1373 K) The standard Gibbs energy changes for IRI and IIU equilibria have been deduced from the measured oxygen potentials. Since ilmenite contains small amounts of Ti³+ ions, a correction for the activity of FeTiO3 has been incorporated by assuming ideal mixing on each cation sublattice in the FeTiO3-Ti2O3 system. Similarly, the ulvospinel contains some Fe³+ ions and a correction for the activity of Fe2TiO4 has been included by modeling the Fe2TiO4-Fe3O4 system. The third-law analysis of the results obtained for IRI equilibrium gives ΔH 298 0 = -575 (±1.0) kJ mol-1 and for IIU equilibrium yields ΔH 298 0 = -523.7 (±0.7) kJ mol−1}. The present results suggest that Fe2+ and Ti4+ cations mix almost ideally on the octahedral site of spinel lattice in Fe2TiO4, giving rise to a configurational contribution of 2R In 2 (11.5256 J mol-1 K-1) to the entropy of Fe2TiO4.

Effect of hole localization on Madelung potential of YBa2Cu3O7

The Madelung potential and formation energy of the superconducting compound YBa2Cu3O7 have been computed for hole localization at different sites in the crystal. The cases considered include Cu3+ ion at Cu(1) and Cu(2) sites, O− ion at O(1), O(2), O(3) and O(4) sites and combinations of O− and Cu3+ ions at O(4) and Cu(1) and O(2,3) and Cu(2) sites. The two lowest-energy configurations correspond to Cu3+ ion at Cu(1) site and O− ion at O(4) site. The difference in formation energy between those configurations is relatively small. The next preferred configuration corresponds to simultaneous partial localization of the hole at Cu (1) site and O(1) site. Other configurations are much less stable. The results suggest a resonating or fluctuating valence model for YBa2Cu3O7.

Phase Relations in the System Cu-Gd-O and Gibbs Energy of Formation of CuGd204

The phase relations in the system Cu-Gd-O have been determined at 1273 K by X-ray diffrac- tion, optical microscopy, and electron microprobe analysis of samples equilibrated in quartz ampules and in pure oxygen. Only one ternary compound, CuGd2O4, was found to be stable. The Gibbs free energy of formation of this compound has been measured using the solid-state cell Pt, Cu2O + CuGd2O4 + Gd2O3 // (Y2O3) ZrO2 // CuO + Cu2O, Pt in the temperature range of 900 to 1350 K. For the formation of CuGd2O4 from its binary component oxides, CuO (s) + Gd2O3 (s) → CuGd2O4 (s) ΔG° = 8230 - 11.2T (±50) J mol-1 Since the formation is endothermic, CuGd2O4 becomes thermodynamically unstable with respect to CuO and Gd2O3 below 735 K. When the oxygen partial pressure over CuGd2O4 is lowered, it decomposes according to the reaction 4CuGd2O4 (s) → 4Gd2O3 (s) + 2Cu2O (s) + O2 (g) for which the equilibrium oxygen potential is given by Δμo 2 = −227,970 + 143.2T (±500) J mol−1 An oxygen potential diagram for the system Cu-Gd-O at 1273 K is presented.

Potential of bio-DME as a transportation fuel for India

This article discusses the potential of bio-dimethyl ether (DME) as a promising fuel for India in the transportation sector where a majority of imported petroleum in the form of diesel is used. Specifically, the suitability of DME in terms of its properties vis-a-vis those of diesel, ability to liquefy DME at low pressures similar to liquefied petroleum gas (LPG), and ease of production from renewable feedstock (biomass), and most importantly, very low emissions including near-zero soot levels are some of the features that make it an attractive option. A detailed review presents the state-of-the-art on various aspects such as estimates of potential bio-DME production, methods of synthesis of bio-DME, important physicochemical properties, fuel-injection system-related concerns (both conventional and common-rail system), fuel spray characteristics which have a direct bearing on the engine performance, and finally, exhaust emissions. Future research directions covering all aspects from production to utilization are summarized (C) 2010 American Institute of Physics. doi:10.1063/1.3489529]

Vinyl Monomer Based Polyperoxides as Potential Initiators for Radical Polymerization: An Exploratory Investigation with Poly(.alpha.-methylstyrene peroxide)

We describe the use of poly(alpha-methylstyrene peroxide) (P alpha MSP), an alternating copolymer of alpha-methylstyrene and oxygen, as initiator for the radical polymerization of vinyl monomers. Thermal decomposition of P alpha MSP in 1,4-dioxane follows first-order kinetics with an activation energy (E(a)) of 34.6 kcal/mol. Polymerization of methyl methacrylate (MMA) and styrene using P alpha MSP as an initiator was carried out in the temperature range 60-90 degrees C. The kinetic order with respect to the initiator and the monomer was close to 0.5 and 1.0, respectively, for both monomers. The E(a) for the polymerization was 20.6 and 22.9 kcal/mol for MMA and styrene, respectively. The efficiency of P alpha MSP was found to be in the range 0.02-0.04. The low efficiency of P alpha MSP was explained in terms of the unimolecular decomposition of the alkoxy radicals which competes with primary radical initiation. The presence of peroxy segments in the main chain of PMMA and polystyrene was confirmed from spectroscopic and DSC studies. R(i)'/2I values for P alpha MSP compared to that of BPO at 80 degrees C indicate that P alpha MSP can be used as an effective high-temperature initiator.

Effect of magnetic field on the binding energy of a hydrogenic donor in a vertical bar z vertical bar(2/3) quantum well

We have calculated the binding energy of a hydrogenic donor in a quantum well with potential shape proportional to \z\(2/3) as a function of the width of the quantum well and the barrier height under an applied uniform magnetic field along the a axis. As the well width decreases, the binding energy increases initially up to a critical well width (which is nearly the same for all magnetic fields) at which there is a turnover. The results are qualitatively similar to those of a hydrogenic donor in a rectangular well. We have also calculated [rho(2)](1/2) and [z(2)](1/2) for the donor electron. [rho(2)](1/2) is found to be strongly dependent on the magnetic field for a given well width and weakly dependent on the well width and the barrier height, for a given value of magnetic field [z(2)](1/2) is weakly dependent on the applied magnetic field. The probability of finding the donor electron inside the well shows a rapid decrease as the well width is reduced at nearly the well width at which the binding energy shows a maximum.

Determination of the Gibbs energy of diamond using a solid state cell

he chemical potential of carbon in diamond, relative to its value in graphite, has been directly determined using a solid state electrochemical cell incorporating single crystal CaF2 as the solid electrolyte. The cell can be represented as Pt, C(graphite) + CaC2 + CaF2double vertical barCaF2double vertical barCaF2 + CaC2 + C(diamond), Pt The reversible emf of this cell is directly related by the Nernst equation to the Gibbs free energy change for the conversion of diamond to graphite. The difference in the chemical potential of carbon in the two crystal structures varies linearly with temperature in the range 940 to 1260 K ?C(diamond) ? ?C(graphite) = 1100 + 4.64T (±50) J mol?1 On the average, the values given by the equation are 320 J mol?1 less positive than the currently accepted ones based on calorimetric studies. The difference is primarily in the enthalpy term.

Potential and economics of forestry options for carbon sequestration in India

There is a need to understand the carbon (C) sequestration potential of the forestry option and its financial implications for each country.In India the C emissions from deforestation are estimated to be nearly offset by C sequestration in forests under succession and tree plantations. India has nearly succeeded in stabilizing the area under forests and has adequate forest conservation strategies. Biomass demands for softwood, hardwood and firewood are estimated to double or treble by the year 2020. A set of forestry options were developed to meet the projected biomass needs, and keeping in mind the features of land categories available, three scenarios were developed: potential; demand-driven; and programme-driven scenarios. Adoption of the demand-driven scenario, targeted at meeting the projected biomass needs, is estimated to sequester 78 Mt of C annually after accounting for all emissions resulting from clearfelling and end use of biomass. The demand-driven scenario is estimated to offset 50% of national C emission at 1990 level. The cost per t of C sequestered for forestry options is lower than the energy options considered. The annual investment required for implementing the demand-driven scenario is estimated to be US$ 2.1 billion for six years and is shown to be feasible. Among forestry options, the ranking based on investment cost per t of C sequestered from least cost to highest cost is; natural regeneration-agro-forestry-enhanced natural regeneration (< US$ 2.5/t C)-timber-community-softwood forestry (US$ 3.3 to 7.3 per t of C).

Prediction of the senses of helical amphiphilic assemblies from effective intermolecular pair potential: Studies on chiral monolayers and bilayers

It is well-known that the senses (or the handedness) of the helical assemblies formed from compressed monolayers and bilayers of chiral amphiphiles are highly specific about the chirality of the monomers concerned. We present here a molecular approach that can successfully predict the senses of such helical morphologies. The present approach is based on a reduced tractable description in terms of an effective pair potential (EPP) which depends on the distance of separation and the relative orientations of the two amphiphiles. This approach explicitly considers the pairwise intermolecular interactions between the groups attached to the chiral centers of the two neighboring amphiphiles. It is found that for a pair of the same kind of enantiomers the minimum energy configuration favors a twist angle between molecules and that this twist from neighbor to neighbor gives rise to the helicity of the aggregate. From the known twist angles at the minimum energy configuration the successive arrangement of an array of molecules can be predicted. Therefore, the sense of the helicity can be predicted from the molecular interactions. The predicted senses of the helical structures are in complete agreement with all known experimental results.