Looking at Chemistry in Hard to See Places: Understanding Energy Conversion at 1400° Celsius

Solid oxide fuel cells (SOFCs) are promising devices for stationary and portable power and heat generation, because they can use complex fuels such as hydro-carbons, CO, and alcohols. Extreme, non-equilibrium conditions and high tem-peratures (≥ 700 ˚C) required for SOFC operation hamper efforts to understand the mechanisms of component degradation in SOFCs. This talk focuses on new insights into SOFC chemistry and the conversion of carbon-containing fuels (both hydrocarbons and oxygenated) into electricity, carbon dioxide and water, gleaned from a combination of techniques including electrochemical impedance spectroscopy, voltammetry, and vibrational Raman scattering.

A review of air–ice chemical and physical interactions (AICI): liquids, quasi-liquids, and solids in snow

Snow in the environment acts as a host to rich chemistry and provides a matrix for physical exchange of contaminants within the ecosystem. The goal of this review is to summarise the current state of knowledge of physical processes and chemical reactivity in surface snow with relevance to polar regions. It focuses on a description of impurities in distinct compartments present in surface snow, such as snow crystals, grain boundaries, crystal surfaces, and liquid parts. It emphasises the microscopic description of the ice surface and its link with the environment. Distinct differences between the disordered air–ice interface, often termed quasi-liquid layer, and a liquid phase are highlighted. The reactivity in these different compartments of surface snow is discussed using many experimental studies, simulations, and selected snow models from the molecular to the macro-scale. Although new experimental techniques have extended our knowledge of the surface properties of ice and their impact on some single reactions and processes, others occurring on, at or within snow grains remain unquantified. The presence of liquid or liquid-like compartments either due to the formation of brine or disorder at surfaces of snow crystals below the freezing point may strongly modify reaction rates. Therefore, future experiments should include a detailed characterisation of the surface properties of the ice matrices. A further point that remains largely unresolved is the distribution of impurities between the different domains of the condensed phase inside the snowpack, i.e. in the bulk solid, in liquid at the surface or trapped in confined pockets within or between grains, or at the surface. While surface-sensitive laboratory techniques may in the future help to resolve this point for equilibrium conditions, additional uncertainty for the environmental snowpack may be caused by the highly dynamic nature of the snowpack due to the fast metamorphism occurring under certain environmental conditions. Due to these gaps in knowledge the first snow chemistry models have attempted to reproduce certain processes like the long-term incorporation of volatile compounds in snow and firn or the release of reactive species from the snowpack. Although so far none of the models offers a coupled approach of physical and chemical processes or a detailed representation of the different compartments, they have successfully been used to reproduce some field experiments. A fully coupled snow chemistry and physics model remains to be developed.

XMapTools: a MATLAB©-based program for electron microprobe X-ray image processing and geothermobarometry

XMapTools is a MATLAB©-based graphical user interface program for electron microprobe X-ray image processing, which can be used to estimate the pressure–temperature conditions of crystallization of minerals in metamorphic rocks. This program (available online at http://www.xmaptools.com) provides a method to standardize raw electron microprobe data and includes functions to calculate the oxide weight percent compositions for various minerals. A set of external functions is provided to calculate structural formulae from the standardized analyses as well as to estimate pressure–temperature conditions of crystallization, using empirical and semi-empirical thermobarometers from the literature. Two graphical user interface modules, Chem2D and Triplot3D, are used to plot mineral compositions into binary and ternary diagrams. As an example, the software is used to study a high-pressure Himalayan eclogite sample from the Stak massif in Pakistan. The high-pressure paragenesis consisting of omphacite and garnet has been retrogressed to a symplectitic assemblage of amphibole, plagioclase and clinopyroxene. Mineral compositions corresponding to ~165,000 analyses yield estimates for the eclogitic pressure–temperature retrograde path from 25 kbar to 9 kbar. Corresponding pressure–temperature maps were plotted and used to interpret the link between the equilibrium conditions of crystallization and the symplectitic microstructures. This example illustrates the usefulness of XMapTools for studying variations of the chemical composition of minerals and for retrieving information on metamorphic conditions on a microscale, towards computation of continuous pressure–temperature-and relative time path in zoned metamorphic minerals not affected by post-crystallization diffusion.

Carbon and deuterium isotopes and C ratios in gas hydrates of ODP Leg 204 sites

Sediments at the southern summit of Hydrate Ridge display two distinct modes of gas hydrate occurrence. The dominant mode is associated with active venting of gas exsolved from the accretionary prism and leads to high concentrations (15%-40% of pore space) of gas hydrate in seafloor or near-surface sediments at and around the topographic summit of southern Hydrate Ridge. These near-surface gas hydrates are mainly composed of previously buried microbial methane but also contain a significant (10%-15%) component of thermogenic hydrocarbons and are overprinted with microbial methane currently being generated in shallow sediments. Focused migration pathways with high gas saturation (>65%) abutting the base of gas hydrate stability create phase equilibrium conditions that permit the flow of a gas phase through the gas hydrate stability zone. Gas seepage at the summit supports rapid growth of gas hydrates and vigorous anaerobic methane oxidation. The other mode of gas hydrate occurs in slope basins and on the saddle north of the southern summit and consists of lower average concentrations (0.5%-5%) at greater depths (30-200 meters below seafloor [mbsf]) resulting from the buildup of in situ-generated dissolved microbial methane that reaches saturation levels with respect to gas hydrate stability at 30-50 mbsf. Net rates of sulfate reduction in the slope basin and ridge saddle sites estimated from curve fitting of concentration gradients are 2-4 mmol/m**3/yr, and integrated net rates are 20-50 mmol/m**2/yr. Modeled microbial methane production rates are initially 1.5 mmol/m**3/yr in sediments just beneath the sulfate reduction zone but rapidly decrease to rates of <0.1 mmol/m**3/yr at depths >100 mbsf. Integrated net rates of methane production in sediments away from the southern summit of Hydrate Ridge are 25-80 mmol/m**2/yr. Anaerobic methane oxidation is minor or absent in cored sediments away from the summit of southern Hydrate Ridge. Ethane-enriched Structure I gas hydrate solids are buried more rapidly than ethane-depleted dissolved gas in the pore water because of advection from compaction. With subsidence beneath the gas hydrate stability zone, the ethane (mainly of low-temperature thermogenic origin) is released back to the dissolved gas-free gas phases and produces a discontinuous decrease in the C1/C2 vs. depth trend. These ethane fractionation effects may be useful to recognize and estimate levels of gas hydrate occurrence in marine sediments.

Transport Sector and Regional Price Differentials: A SCGE Model for Chinese Provinces

With regression formulas replaced by equilibrium conditions, a spatial CGE model can substantially reduce data requirements. Detailed regional analyses are thus possible in countries where only limited regional statistics are available. While regional price differentials play important roles in multi-regional settings, transport does not receive much attention in existing models. This paper formulates a spatial CGE model that explicitly considers the transport sector and FOB/CIF prices. After describing the model, performance of our model is evaluated by comparing the benchmark equilibrium for China with survey-based regional I-O and interregional I-O tables for 1987. The structure of Chinese economies is summarized using information obtained from the benchmark equilibrium computation. This includes regional and sectoral production distributions and price differentials. The equilibrium for 1997 facilitates discussion of changes in regional economic structures that China has experienced in the decade.

Resilience analysis of networked systems-of-systems based on structural and dynamic interdependencies

Critical infrastructures support everyday activities in modern societies, facilitating the exchange of services and quantities of various nature. Their functioning is the result of the integration of diverse technologies, systems and organizations into a complex network of interconnections. Benefits from networking are accompanied by new threats and risks. In particular, because of the increased interdependency, disturbances and failures may propagate and render unstable the whole infrastructure network. This paper presents a methodology of resilience analysis of networked systems of systems. Resilience generalizes the concept of stability of a system around a state of equilibrium, with respect to a disturbance and its ability of preventing, resisting and recovery. The methodology provides a tool for the analysis of off-equilibrium conditions that may occur in a single system and propagate through the network of dependencies. The analysis is conducted in two stages. The first stage of the analysis is qualitative. It identifies the resilience scenarios, i.e. the sequence of events, triggered by an initial disturbance, which include failures and the system response. The second stage is quantitative. The most critical scenarios can be simulated, for the desired parameter settings, in order to check if they are successfully handled, i.e recovered to nominal conditions, or they end into the network failure. The proposed methodology aims at providing an effective support to resilience-informed design.

Intermediate Band to Conduction Band optical absorption in ZnTe:O

ZnTe doped with high concentrations of oxygen has been proposed in previous works as intermediate band (IB) material for photovoltaic applications. The existence of extra optical transitions related to the presence of an IB has already been demonstrated in this material and it has been possible to measure the absorption coefficient of the transitions from the valence band (VB) to the IB. In this work we present the first measurement of the absorption coefficient associated to transitions from the IB to the conduction band (CB) in ZnTe:O. The samples used are 4 ?m thick ZnTe layers with or without O in a concentration ~ 1019 cm-3, which have been grown on semi-insulating GaAs substrates by molecular beam epitaxy (MBE). The IB-CB absorption coefficient peaks for photon energies ~ 0.4 eV. It is extracted from reflectance and transmittance spectra measured using Fourier Transform Infrared (FTIR) spectroscopy. Under typical FTIR measurement conditions (low light intensity, broadband spectrum) the absorption coefficient in IB-to-CB transitions reaches 700 cm-1. This is much weaker than the one observed for VB-IB absorption. This result is consistent with the fact that the IB is expected to be nearly empty of electrons under equilibrium conditions in ZnTe(:O). The absorption for VB to IB transitions is also observed in the same samples through reflectance measurements performed in the visible range using a monochromator. These measurements are compared with the quantum efficiency (QE) from solar cells fabricated under similar conditions.

The breaking of axisymmetric slender liquid bridges

Liquids held by surface tension forces can bridge the gap between two solid bodies placed not too far apart from each other. The equilibrium conditions and stability criteria for static, cylindrical liquid bridges are well known. However, the behaviour of an unstable liquid bridge, regarding both its transition toward breaking and the resulting configuration, is a matter for discussion. The dynamical problem of axisymmetric rupture of a long liquid bridge anchored at two equal coaxial disks is treated in this paper through the adoption of one-dimensional theories which are widely used in capillary jet problems

Intermediate band to conduction band optical absorption in ZnTeO

ZnTe doped with high concentrations of oxygen has been proposed in previous works as an intermediate band (IB) material for photovoltaic applications. The existence of extra optical transitions related to the presence of an IB has already been demonstrated in this material and it has been possible to measure the absorption coefficient of the transitions from the valence band (VB) to the IB. In this study, we present the first measurement of the absorption coefficient associated with transitions from the IB to the conduction band (CB) in ZnTeO. The samples used are 4-mum-thick ZnTe layers with or without O in a concentration ~10 19 cm -3, which have been grown on semiinsulating GaAs substrates by molecular beam epitaxy (MBE). The IB-CB absorption coefficient peaks for photon energies ~0.4 eV. It is extracted from reflectance and transmittance spectra measured using Fourier transform infrared (FTIR) spectroscopy. Under typical FTIR measurement conditions (low light intensity, broadband spectrum), the absorption coefficient in IB-to-CB transitions reaches 700 cm -1. This is much weaker than the one observed for VB-IB absorption. This result is consistent with the fact that the IB is expected to be nearly empty of electrons under equilibrium conditions in ZnTe(O).

Human Factor in Resilience Systems

This paper shows the importance of a holistic comprehension of the Earth as a living planet, where man inhabits and is exposed to environmental incidences of different nature. The aim of the paper here summarized is a reflection on all these concepts and scientific considerations related to the important role of men in the handling of natural hazards. Our Planet is an unstable and dynamical system highly sensitive to initial conditions, as proposed by Chaos theory (González-Miranda 2004); it is a complex organic whole, which responds to minimal variations which can affect several natural phenomena such as plate tectonics, solar flares, fluid turbulences, landscape formation, forest fires, growth and migration of populations and biological evolution. This is known as the “butterfly effect” (Lorenz 1972), which means that a small change of the system causes a chain of events leading to large-scale unpredictable consequences. The aim of this work is dwelling on the importance of the knowledge of these natural and catastrophic geological, biological and human systems so much sensible to equilibrium conditions, to prevent, avoid and mend their effects, and to face them in a resilient way

Evaluación de la seguridad al deslizamiento de una presa de fábrica en un cimiento con una familia de discontinuidades y con un criterio de rotura con ley de fluencia no asociada

En España existen del orden de 1,300 grandes presas, de las cuales un 20% fueron construidas antes de los años 60. El hecho de que existan actualmente una gran cantidad de presas antiguas aún en operación, ha producido un creciente interés en reevaluar su seguridad empleando herramientas nuevas o modificadas que incorporan modelos de fallo teóricos más completos, conceptos geotécnicos más complejos y nuevas técnicas de evaluación de la seguridad. Una manera muy común de abordar el análisis de estabilidad de presas de gravedad es, por ejemplo, considerar el deslizamiento a través de la interfase presa-cimiento empleando el criterio de rotura lineal de Mohr-Coulomb, en donde la cohesión y el ángulo de rozamiento son los parámetros que definen la resistencia al corte de la superficie de contacto. Sin embargo la influencia de aspectos como la presencia de planos de debilidad en el macizo rocoso de cimentación; la influencia de otros criterios de rotura para la junta y para el macizo rocoso (ej. el criterio de rotura de Hoek-Brown); las deformaciones volumétricas que ocurren durante la deformación plástica en el fallo del macizo rocoso (i.e., influencia de la dilatancia) no son usualmente consideradas durante el diseño original de la presa. En este contexto, en la presente tesis doctoral se propone una metodología analítica para el análisis de la estabilidad al deslizamiento de presas de hormigón, considerando un mecanismo de fallo en la cimentación caracterizado por la presencia de una familia de discontinuidades. En particular, se considera la posibilidad de que exista una junta sub-horizontal, preexistente y persistente en el macizo rocoso de la cimentación, con una superficie potencial de fallo que se extiende a través del macizo rocoso. El coeficiente de seguridad es entonces estimado usando una combinación de las resistencias a lo largo de los planos de rotura, cuyas resistencias son evaluadas empleando los criterios de rotura no lineales de Barton y Choubey (1977) y Barton y Bandis (1990), a lo largo del plano de deslizamiento de la junta; y el criterio de rotura de Hoek y Brown (1980) en su versión generalizada (Hoek et al. 2002), a lo largo del macizo rocoso. La metodología propuesta también considera la influencia del comportamiento del macizo rocoso cuando este sigue una ley de flujo no asociada con ángulo de dilatancia constante (Hoek y Brown 1997). La nueva metodología analítica propuesta es usada para evaluar las condiciones de estabilidad empleando dos modelos: un modelo determinista y un modelo probabilista, cuyos resultados son el valor del coeficiente de seguridad y la probabilidad de fallo al deslizamiento, respectivamente. El modelo determinista, implementado en MATLAB, es validado usando soluciones numéricas calculadas mediante el método de las diferencias finitas, empleando el código FLAC 6.0. El modelo propuesto proporciona resultados que son bastante similares a aquellos calculados con FLAC; sin embargo, los costos computacionales de la formulación propuesta son significativamente menores, facilitando el análisis de sensibilidad de la influencia de los diferentes parámetros de entrada sobre la seguridad de la presa, de cuyos resultados se obtienen los parámetros que más peso tienen en la estabilidad al deslizamiento de la estructura, manifestándose además la influencia de la ley de flujo en la rotura del macizo rocoso. La probabilidad de fallo es obtenida empleando el método de fiabilidad de primer orden (First Order Reliability Method; FORM), y los resultados de FORM son posteriormente validados mediante simulaciones de Monte Carlo. Los resultados obtenidos mediante ambas metodologías demuestran que, para el caso no asociado, los valores de probabilidad de fallo se ajustan de manera satisfactoria a los obtenidos mediante las simulaciones de Monte Carlo. Los resultados del caso asociado no son tan buenos, ya que producen resultados con errores del 0.7% al 66%, en los que no obstante se obtiene una buena concordancia cuando los casos se encuentran en, o cerca de, la situación de equilibrio límite. La eficiencia computacional es la principal ventaja que ofrece el método FORM para el análisis de la estabilidad de presas de hormigón, a diferencia de las simulaciones de Monte Carlo (que requiere de al menos 4 horas por cada ejecución) FORM requiere tan solo de 1 a 3 minutos en cada ejecución. There are 1,300 large dams in Spain, 20% of which were built before 1960. The fact that there are still many old dams in operation has produced an interest of reevaluate their safety using new or updated tools that incorporate state-of-the-art failure modes, geotechnical concepts and new safety assessment techniques. For instance, for gravity dams one common design approach considers the sliding through the dam-foundation interface, using a simple linear Mohr-Coulomb failure criterion with constant friction angle and cohesion parameters. But the influence of aspects such as the persistence of joint sets in the rock mass below the dam foundation; of the influence of others failure criteria proposed for rock joint and rock masses (e.g. the Hoek-Brown criterion); or the volumetric strains that occur during plastic failure of rock masses (i.e., the influence of dilatancy) are often no considered during the original dam design. In this context, an analytical methodology is proposed herein to assess the sliding stability of concrete dams, considering an extended failure mechanism in its rock foundation, which is characterized by the presence of an inclined, and impersistent joint set. In particular, the possibility of a preexisting sub-horizontal and impersistent joint set is considered, with a potential failure surface that could extend through the rock mass; the safety factor is therefore computed using a combination of strength along the rock joint (using the nonlinear Barton and Choubey (1977) and Barton and Bandis (1990) failure criteria) and along the rock mass (using the nonlinear failure criterion of Hoek and Brown (1980) in its generalized expression from Hoek et al. (2002)). The proposed methodology also considers the influence of a non-associative flow rule that has been incorporated using a (constant) dilation angle (Hoek and Brown 1997). The newly proposed analytical methodology is used to assess the dam stability conditions, employing for this purpose the deterministic and probabilistic models, resulting in the sliding safety factor and the probability of failure respectively. The deterministic model, implemented in MATLAB, is validated using numerical solution computed with the finite difference code FLAC 6.0. The proposed deterministic model provides results that are very similar to those computed with FLAC; however, since the new formulation can be implemented in a spreadsheet, the computational cost of the proposed model is significantly smaller, hence allowing to more easily conduct parametric analyses of the influence of the different input parameters on the dam’s safety. Once the model is validated, parametric analyses are conducting using the main parameters that describe the dam’s foundation. From this study, the impact of the more influential parameters on the sliding stability analysis is obtained and the error of considering the flow rule is assessed. The probability of failure is obtained employing the First Order Reliability Method (FORM). The probabilistic model is then validated using the Monte Carlo simulation method. Results obtained using both methodologies show good agreement for cases in which the rock mass has a nonassociate flow rule. For cases with an associated flow rule errors between 0.70% and 66% are obtained, so that the better adjustments are obtained for cases with, or close to, limit equilibrium conditions. The main advantage of FORM on sliding stability analyses of gravity dams is its computational efficiency, so that Monte Carlo simulations require at least 4 hours on each execution, whereas FORM requires only 1 to 3 minutes on each execution.

Determinación del riesgo de incendio en una instalación de secado de alperujo = Determination of the risk of fire in an olive-oil mill waste dryer

For the energy valorization of alperujo, residue of the olive oil two phases extraction process, it is necessary to perform a drying process to reduce moisture content from over 60% to less than 10%. In order to reduce primary energy consumption and get an economic return, usually in this kind of drying facilities Gas Turbine CHP is used as a heat source. There have been recently in Spain some fires in this kind of GT-CHP facilities, which have caused high material losses. In some of these fires it has been suggested that the fire was caused by the output of incandescent alperujo in the flue gasesof the drying system. Therefore, the aim of this study is to determine experimentally and analytically under which operational conditions a process of alperujo self-ignition in the drying process can begin, and determine the actual fire hazard in this type of TG-CHP system. For analytical study, the temperature and initial composition of the combustion gases of the Gas Turbine at the entrance of the drying process was calculated and the gas equilibrium conditions reached in contact with the biomass were calculated and, therefore, the temperature of the biomass during the drying process. Moreover, the layer and dust ignition temperature of alperujo has been experimentally determined, according to EN 50281-2-1: 2000. With these results, the operating conditions of the drying process, in which there are real risk of auto-ignition of alperujo have been established.Para la valorización energética del alperujo, residuo del proceso de extracción en dos fases del aceite de oliva, es necesario realizar un proceso de secado para reducir su contenido de humedad de más del 60% al 10% m/m en b.h. Con el fin de reducir el consumo de energía primaria y obtener una rentabilidad económica, normalmente en este tipo de instalaciones de secado se usa la cogeneración con turbina de gas (TG) como fuente de calor. En España en los últimos años han ocurrido algunos casos de incendio en este tipo de instalaciones de cogeneración, que han supuesto pérdidas materiales muy elevadas. Por esta razón, el objetivo de este trabajo es determinar analítica y experimentalmente las condiciones operativas del secadero bajo las cuales podría comenzar un proceso de autoinflamación del alperujo y determinar el riesgo real de incendio en este tipo de instalaciones. Para el estudio analítico, se ha planteado y validado el modelo matemático que permite calcular la temperatura y la composición de los gases de combustión a la entrada y a la salida del secadero, en función de las curvas características de la TG, de las condiciones atmosféricas, del caudal y del grado de humedad de la biomasa tratada. El modelo permite además calcular la temperatura de bulbo húmedo, que es la máxima temperatura que podría alcanzar la biomasa durante el proceso de secado y determinar la cantidad de biomasa que se puede secar completamente en función del caudal y de las condiciones de entrada de los gases de combustión. Con estos resultados y la temperatura mínima de autoinflamación del alperujo determinada experimentalmente siguiendo la norma EN 50281- 2-1:2000, se demuestra que en un proceso de secado de alperujo en condiciones normales de operación no existe riesgo de autoencendido que pueda dar origen a un incendio.

Folding thermodynamics of a model three-helix-bundle protein

The calculated folding thermodynamics of a simple off-lattice three-helix-bundle protein model under equilibrium conditions shows the experimentally observed protein transitions: a collapse transition, a disordered-to-ordered globule transition, a globule to native-state transition, and the transition from the active native state to a frozen inactive state. The cooperativity and physical origin of the various transitions are explored with a single “optimization” parameter and characterized with the Lindemann criterion for liquid versus solid-state dynamics. Below the folding temperature, the model has a simple free energy surface with a single basin near the native state; the surface is similar to that calculated from a simulation of the same three-helix-bundle protein with an all-atom representation [Boczko, E. M. & Brooks III, C. L. (1995) Science 269, 393–396].

Gibbs Energy of Mixing Function: Topological Analysis in Azeotropic Systems

ESAT 2014. 27th European Symposium on Applied Thermodynamics, Eindhoven University of Technology, July 6-9, 2014.

The isotope composition of selenium in chondrites constrains the depletion mechanism of volatile elements in solar system materials

Solar nebula processes led to a depletion of volatile elements in different chondrite groups when compared to the bulk chemical composition of the solar system deduced from the Sun's photosphere. For moderately-volatile elements, this depletion primarily correlates with the element condensation temperature and is possibly caused by incomplete condensation from a hot solar nebula, evaporative loss from the precursor dust, and/or inherited from the interstellar medium. Element concentrations and interelement ratios of volatile elements do not provide a clear picture about responsible mechanisms. Here, the abundance and stable isotope composition of the moderately- to highly-volatile element Se are investigated in carbonaceous, ordinary, and enstatite chondrites to constrain the mechanism responsible for the depletion of volatile elements in planetary bodies of the inner solar system and to define a δ(82/78)Se value for the bulk solar system. The δ(82/78)Se of the studied chondrite falls are identical within their measurement uncertainties with a mean of −0.20±0.26‰ (2 s.d., n=14n=14, relative to NIST SRM 3149) despite Se abundance depletions of up to a factor of 2.5 with respect to the CI group. The absence of resolvable Se isotope fractionation rules out a kinetic Rayleigh-type incomplete condensation of Se from the hot solar nebula or partial kinetic evaporative loss on the precursor material and/or the parent bodies. The Se depletion, if acquired during partial condensation or evaporative loss, therefore must have occurred under near equilibrium conditions to prevent measurable isotope fractionation. Alternatively, the depletion and cooling of the nebula could have occurred simultaneously due to the continuous removal of gas and fine particles by the solar wind accompanied by the quantitative condensation of elements from the pre-depleted gas. In this scenario the condensation of elements does not require equilibrium conditions to avoid isotope fractionation. The results further suggest that the processes causing the high variability of Se concentrations and depletions in ordinary and enstatite chondrites did not involve any measurable isotope fractionation. Different degrees of element depletions and isotope fractionations of the moderately-volatile elements Zn, S, and Se in ordinary and enstatite chondrites indicate that their volatility is controlled by the thermal stabilities of their host phases and not by the condensation temperature under canonical nebular conditions.