Candidate radial-inflow turbines and high-density working fluids for geothermal power systems

Optimisation of Organic Rankine Cycles (ORCs) for binary-cycle geothermal applications could play a major role in the competitiveness of low to moderate temperature geothermal resources. Part of this optimisation process is matching cycles to a given resource such that power output can be maximised. Two major and largely interrelated components of the cycle are the working fluid and the turbine. Both components need careful consideration. Due to the temperature differences in geothermal resources a one-size-fits-all approach to surface power infrastructure is not appropriate. Furthermore, the traditional use of steam as a working fluid does not seem practical due to the low temperatures of many resources. A variety of organic fluids with low boiling points may be utilised as ORC working fluids in binary power cycle loops. Due to differences in thermodynamic properties, certain fluids are able to extract more heat from a given resource than others over certain temperature and pressure ranges. This enables the tailoring of power cycle infrastructure to best match the geothermal resource through careful selection of the working fluid and turbine design optimisation to yield the optimum overall cycle performance. This paper presents the rationale for the use of radial-inflow turbines for ORC applications and the preliminary design of several radial-inflow turbines based on a selection of promising ORC cycles using five different high-density working fluids: R134a, R143a, R236fa, R245fa and n-Pentane at sub- or trans-critical conditions. Numerous studies published compare a variety of working fluids for various ORC configurations. However, there is little information specifically pertaining to the design and implementation of ORCs using realistic radial turbine designs in terms of pressure ratios, inlet pressure, rotor size and rotational speed. Preliminary 1D analysis leads to the generation of turbine designs for the various cycles with similar efficiencies (77%) but large differences in dimensions (139289 mm rotor diameter). The highest performing cycle (R134a) was found to produce 33% more net power from a 150°C resource flowing at 10 kg/s than the lowest performing cycle (n-Pentane).

Integrated TEM, XRD and electron-microprobe investigation of mixed-layer chlorite smectite from the Point-Sal ophiolite, California

Five basalt samples from the Point Sal ophiolite, California, were examined using HRTEM and AEM in order to compare observations with interpretations of XRD patterns and microprobe analyses. XRD data from ethylene-glycol-saturated samples indicate the following percentages of chlorite in mixed-layer chlorite-smectite identified for each specimen: (i) L2036 almost-equal-to 50%, (ii) L2035 almost-equal-to 70 and 20%, (iii) 1A-13 almost-equal-to 70%, (iv) 1B-42 almost-equal-to 70%, and (v) 1B-55 = 100%. Detailed electron microprobe analyses show that 'chlorite' analyses with high Si, K, Na and Ca contents are the result of interlayering with smectite-like layers. The Fe/(Fe + Mg) ratios of mixed-layer phyllosilicates from Point Sal samples are influenced by the bulk rock composition, not by the percentage of chlorite nor the structure of the phyllosilicate. Measurements of lattice-fringe images indicate that both smectite and chlorite layers are present in the Point Sal samples in abundances similar to those predicted with XRD techniques and that regular alternation of chlorite and smectite occurs at the unit-cell scale. Both 10- and 14-angstrom layers were recorded with HRTEM and interpreted to be smectite and chlorite, respectively. Regular alternation of chlorite and smectite (24-angstrom periodicity) occurs in upper lava samples L2036 and 1A-13, and lower lava sample 1B-42 for as many as seven alternations per crystallite with local layer mistakes. Sample L2035 shows disordered alternation of chlorite and smectite, with juxtaposition of smectite-like layers, suggesting that randomly interlayered chlorite (< 0.5)-smectite exists. Images of lower lava sample 1B-55 show predominantly 14-angstrom layers. Units of 24 angstrom tend to cluster in what may otherwise appear to be disordered mixtures, suggesting the existence of a corrensite end-member having thermodynamic significance.

Effect of Fe3O4 addition on removal of ammonium by zeolite NaA

Magnetic zeolite NaA with different Fe3O4 loadings was prepared by hydrothermal synthesis based on metakaolin and Fe3O4. The effect of added Fe3O4 on the removal of ammonium by zeolite NaA was investigated by varying the Fe3O4 loading, pH, adsorption temperature, initial concentration, adsorption time. Langmuir, Freundlich, and pseudo-second-order modeling were used to describe the nature and mechanism of ammonium ion exchange using both zeolite and magnetic zeolite. Thermodynamic parameters such as change in Gibbs free energy, enthalpy and entropy were calculated. The results show that all the selected factors affect the ammonium ion exchange by zeolite and magnetic zeolite, however, the added Fe3O4 apparently does not affect the ion exchange performance of zeolite to the ammonium ion. Freundlich model provides a better description of the adsorption process than Langmuir model. Moreover, kinetic analysis indicates the exchange of ammonium on the two materials follows a pseudosecond-order model. Thermodynamic analysis makes it clear that the adsorption process of ammonium is spontaneous and exothermic. Regardless of kinetic or thermodynamic analysis, all the results suggest that no considerable effect on the adsorption of the ammonium ion by zeolite is found after the addition of Fe3O4. According to the results, magnetic zeolite NaA can be used for the removal of ammonium due to the good adsorption performance and easy separation method from aqueous solution.

Computational study of methyl derivatives of ammonia borane for hydrogen storage

The structures and thermodynamic properties of methyl derivatives of ammonia–borane (BH3NH3, AB) have been studied with the frameworks of density functional theory and second-order Møller–Plesset perturbation theory. It is found that, with respect to pure AB, methyl ammonia–boranes show higher complexation energies and lower reaction enthalpies for the release of H2, together with a slight increment of the activation barrier. These results indicate that the methyl substitution can enhance the reversibility of the system and prevent the formation of BH3/NH3, but no enhancement of the release rate of H2 can be expected.

Ab initio studies of hydrogen desorption from low index magnesium hydride surface

The low index Magnesium hydride surfaces, MgH2(0 0 1) and MgH2(1 1 0), have been studied by ab intio Density Functional Theory (DFT) calculations. It was found that the MgH2(1 1 0) surface is more stable than MgH2(0 0 1) surface, which is in good agreement with the experimental observation. The H2 desorption barriers vary depending on the crystalline surfaces that are exposed and also the specific H atom sites involved – they are found to be generally high, due to the thermodynamic stability of the MgH2 system, and are larger for the MgH2(0 0 1) surface. The pathway for recombinative desorption of one in-plane and one bridging H atom from the MgH2(1 1 0) surface was found to be the lowest energy barrier amongst those computed (172 KJ/mol) and is in good agreement with the experimental estimates.

Thermal treatment of natural goethite : thermal transformation and physical properties

XRD (X-ray diffraction), XRF (X-ray fluorescence), TG (thermogravimetry), FT-IES (Fourier transform infrared emission spectroscopy), FESEM (field emission scanning electron microscope), TEM (transmission electron microscope) and nitrogen–adsorption–desorption analysis were used to characterize the composition and thermal evolution of the structure of natural goethite. The in situ FT-IES demonstrated the start temperature (250 °C) of the transformation of natural goethite to hematite and the thermodynamic stability of protohematite between 250 and 600 °C. The heated products showed a topotactic relationship to the original mineral based on SEM analysis. Finally, the nitrogen–adsorption–desorption isotherm provided the variation of surface area and pore size distribution as a function of temperature. The surface area displayed a remarkable increase up to 350 °C, and then decreased above this temperature. The significant increase in surface area was attributed to the formation of regularly arranged slit-shaped micropores running parallel to elongated direction of hematite microcrystal. The main pore size varied from 0.99 nm to 3.5 nm when heating temperature increases from 300 to 400 °C. The hematite derived from heating goethite possesses high surface area and favors the possible application of hematite as an adsorbent as well as catalyst carrier.

Sintering and densification mechanisms of Y2BaCuO5 pellets

The sintering and densification of Y2BaCuO5 (Y-211) pellets made from powders with different characteristics have been investigated in the temperature range 1000-1140°C. A pellet made from powder containing Ba-rich secondary phases shows very early liquid-assisted sintering and densification and clear evidence of exaggerated grain growth. The melting of BaCuO2 and YBa2Cu3O7-δ (Y-123) secondary phases increases the rate of densification of Y-211 pellets made from other powders at temperatures above 1025-1030°C. All the liquid produced by the melting of the latter phases recrystallizes as intergranular layers of Y-123. These intergranular layers account for the darker appearance and for measurable electrical conductivities at room temperature of the pellets sintered at the higher temperatures. The development of exaggerated grain growth within a uniform fine-grained matrix opens the possibility of using controlled secondary recrystallization to obtain large single domains of Y-211, provided that the trapping of porosity can be avoided or minimized. © 1999 Elsevier Science S.A.

Auto-inhibition of hydrogen gas evolution on gold in aqueous acid solution

It was demonstrated recently that dramatic changes in the redox behaviour of gold/aqueous solution interfaces may be observed following either cathodic or thermal electrode pretreatment. Further work on the cathodic pretreatment of gold in acid solution revealed that as the activity of the gold surface was increased, its performance as a substrate for hydrogen gas evolution under constant potential conditions deteriorated. The change in activity of the gold atoms at the interface, which was attributed to a hydrogen embrittlement process (the occurrence of the latter was subsequently checked by surface microscopy), was confirmed, as in earlier work, by the appearance of a substantial anodic peak at ca. 0.5 V (RHE) in a post-activation positive sweep. Changes in the catalytic activity of a metal surface reflect the fact that the structure (or topography), thermodynamic activity and electronic properties of a surface are dependent not only on pretreatment but also, in the case of the hydrogen evolution reaction, vary with time during the course of reaction. As will be reported shortly, similar (and often more dramatic) time-dependent behaviour was observed for hydrogen gas evolution on other metal electrodes.

The ECOS Project. An interactive simulation to promote an ECO-friendly universe.

Custom designed for display on the Cube Installation situated in the new Science and Engineering Centre (SEC) at QUT, the ECOS project is a playful interface that uses real-time weather data to simulate how a five-star energy building operates in climates all over the world. In collaboration with the SEC building managers, the ECOS Project incorporates energy consumption and generation data of the building into an interactive simulation, which is both engaging to users and highly informative, and which invites play and reflection on the roles of green buildings. ECOS focuses on the principle that humans can have both a positive and negative impact on ecosystems with both local and global consequence. The ECOS project draws on the practice of Eco-Visualisation, a term used to encapsulate the important merging of environmental data visualization with the philosophy of sustainability. Holmes (2007) uses the term Eco-Visualisation (EV) to refer to data visualisations that ‘display the real time consumption statistics of key environmental resources for the goal of promoting ecological literacy’. EVs are commonly artifacts of interaction design, information design, interface design and industrial design, but are informed by various intellectual disciplines that have shared interests in sustainability. As a result of surveying a number of projects, Pierce, Odom and Blevis (2008) outline strategies for designing and evaluating effective EVs, including ‘connecting behavior to material impacts of consumption, encouraging playful engagement and exploration with energy, raising public awareness and facilitating discussion, and stimulating critical reflection.’ Consequently, Froehlich (2010) and his colleagues also use the term ‘Eco-feedback technology’ to describe the same field. ‘Green IT’ is another variation which Tomlinson (2010) describes as a ‘field at the juncture of two trends… the growing concern over environmental issues’ and ‘the use of digital tools and techniques for manipulating information.’ The ECOS Project team is guided by these principles, but more importantly, propose an example for how these principles may be achieved. The ECOS Project presents a simplified interface to the very complex domain of thermodynamic and climate modeling. From a mathematical perspective, the simulation can be divided into two models, which interact and compete for balance – the comfort of ECOS’ virtual denizens and the ecological and environmental health of the virtual world. The comfort model is based on the study of psychometrics, and specifically those relating to human comfort. This provides baseline micro-climatic values for what constitutes a comfortable working environment within the QUT SEC buildings. The difference between the ambient outside temperature (as determined by polling the Google Weather API for live weather data) and the internal thermostat of the building (as set by the user) allows us to estimate the energy required to either heat or cool the building. Once the energy requirements can be ascertained, this is then balanced with the ability of the building to produce enough power from green energy sources (solar, wind and gas) to cover its energy requirements. Calculating the relative amount of energy produced by wind and solar can be done by, in the case of solar for example, considering the size of panel and the amount of solar radiation it is receiving at any given time, which in turn can be estimated based on the temperature and conditions returned by the live weather API. Some of these variables can be altered by the user, allowing them to attempt to optimize the health of the building. The variables that can be changed are the budget allocated to green energy sources such as the Solar Panels, Wind Generator and the Air conditioning to control the internal building temperature. These variables influence the energy input and output variables, modeled on the real energy usage statistics drawn from the SEC data provided by the building managers.

A stochastic thermostat algorithm for coarse-grained thermomechanical modeling of large-scale soft matters : theory and application to microfilaments

As all-atom molecular dynamics method is limited by its enormous computational cost, various coarse-grained strategies have been developed to extend the length scale of soft matters in the modeling of mechanical behaviors. However, the classical thermostat algorithm in highly coarse-grained molecular dynamics method would underestimate the thermodynamic behaviors of soft matters (e.g. microfilaments in cells), which can weaken the ability of materials to overcome local energy traps in granular modeling. Based on all-atom molecular dynamics modeling of microfilament fragments (G-actin clusters), a new stochastic thermostat algorithm is developed to retain the representation of thermodynamic properties of microfilaments at extra coarse-grained level. The accuracy of this stochastic thermostat algorithm is validated by all-atom MD simulation. This new stochastic thermostat algorithm provides an efficient way to investigate the thermomechanical properties of large-scale soft matters.

Multiscale coupling and multiphysics approaches in earth sciences : theory

Measuring Earth material behaviour on time scales of millions of years transcends our current capability in the laboratory. We review an alternative path considering multiscale and multiphysics approaches with quantitative structure-property relationships. This approach allows a sound basis to incorporate physical principles such as chemistry, thermodynamics, diffusion and geometry-energy relations into simulations and data assimilation on the vast range of length and time scales encountered in the Earth. We identify key length scales for Earth systems processes and find a substantial scale separation between chemical, hydrous and thermal diffusion. We propose that this allows a simplified two-scale analysis where the outputs from the micro-scale model can be used as inputs for meso-scale simulations, which then in turn becomes the micro-model for the next scale up. We present two fundamental theoretical approaches to link the scales through asymptotic homogenisation from a macroscopic thermodynamic view and percolation renormalisation from a microscopic, statistical mechanics view.

Entropic bounds for multi-scale and multi-physics coupling in earth sciences

The ability to understand and predict how thermal, hydrological,mechanical and chemical (THMC) processes interact is fundamental to many research initiatives and industrial applications. We present (1) a new Thermal– Hydrological–Mechanical–Chemical (THMC) coupling formulation, based on non-equilibrium thermodynamics; (2) show how THMC feedback is incorporated in the thermodynamic approach; (3) suggest a unifying thermodynamic framework for multi-scaling; and (4) formulate a new rationale for assessing upper and lower bounds of dissipation for THMC processes. The technique is based on deducing time and length scales suitable for separating processes using a macroscopic finite time thermodynamic approach. We show that if the time and length scales are suitably chosen, the calculation of entropic bounds can be used to describe three different types of material and process uncertainties: geometric uncertainties,stemming from the microstructure; process uncertainty, stemming from the correct derivation of the constitutive behavior; and uncertainties in time evolution, stemming from the path dependence of the time integration of the irreversible entropy production. Although the approach is specifically formulated here for THMC coupling we suggest that it has a much broader applicability. In a general sense it consists of finding the entropic bounds of the dissipation defined by the product of thermodynamic force times thermodynamic flux which in material sciences corresponds to generalized stress and generalized strain rates, respectively.

A new approach to prepare ZVI and its application in removal of Cr(VI) from aqueous solution

Zero valent iron (ZVI) was prepared by reducing natural goethite (NG-ZVI) and synthetic goethite (SG-ZVI) in hydrogen at 550 °C. XRD, TEM, FESEM/EDS and specific surface area (SSA) and pore analyser were used to characterize goethites and reduced goethites. Both NG-ZVI and SG-ZVI with a size of nanoscale to several hundreds of nanometers were obtained by reducing goethites at 550 °C. The reductive capacity of the ZVIs was assessed by removal of Cr(VI) at ambient temperature in comparison with that of commercial iron powder (CIP). The effect of contact time, initial concentration and reaction temperature on Cr(VI) removal was investigated. Furthermore, the uptake mechanism was discussed according to isotherms, thermodynamic analysis and the results of XPS. The results showed that SG-ZVI had the best reductive capacity to Cr(VI) and reduced Cr(VI) to Cr(III). The results suggest that hydrogen reduction is a good approach to prepare ZVI and this type of ZVI is potentially useful in remediating heavy metals as a material of permeable reaction barrier.

Multiscale coupling and multiphysics approaches in earth sciences : applications

Geoscientists are confronted with the challenge of assessing nonlinear phenomena that result from multiphysics coupling across multiple scales from the quantum level to the scale of the earth and from femtoseconds to the 4.5 Ga of history of our planet. We neglect in this review electromagnetic modelling of the processes in the Earth’s core, and focus on four types of couplings that underpin fundamental instabilities in the Earth. These are thermal (T), hydraulic (H), mechanical (M) and chemical (C) processes which are driven and controlled by the transfer of heat to the Earth’s surface. Instabilities appear as faults, folds, compaction bands, shear/fault zones, plate boundaries and convective patterns. Convective patterns emerge from buoyancy overcoming viscous drag at a critical Rayleigh number. All other processes emerge from non-conservative thermodynamic forces with a critical critical dissipative source term, which can be characterised by the modified Gruntfest number Gr. These dissipative processes reach a quasi-steady state when, at maximum dissipation, THMC diffusion (Fourier, Darcy, Biot, Fick) balance the source term. The emerging steady state dissipative patterns are defined by the respective diffusion length scales. These length scales provide a fundamental thermodynamic yardstick for measuring instabilities in the Earth. The implementation of a fully coupled THMC multiscale theoretical framework into an applied workflow is still in its early stages. This is largely owing to the four fundamentally different lengths of the THMC diffusion yardsticks spanning micro-metre to tens of kilometres compounded by the additional necessity to consider microstructure information in the formulation of enriched continua for THMC feedback simulations (i.e., micro-structure enriched continuum formulation). Another challenge is to consider the important factor time which implies that the geomaterial often is very far away from initial yield and flowing on a time scale that cannot be accessed in the laboratory. This leads to the requirement of adopting a thermodynamic framework in conjunction with flow theories of plasticity. This framework allows, unlike consistency plasticity, the description of both solid mechanical and fluid dynamic instabilities. In the applications we show the similarity of THMC feedback patterns across scales such as brittle and ductile folds and faults. A particular interesting case is discussed in detail, where out of the fluid dynamic solution, ductile compaction bands appear which are akin and can be confused with their brittle siblings. The main difference is that they require the factor time and also a much lower driving forces to emerge. These low stress solutions cannot be obtained on short laboratory time scales and they are therefore much more likely to appear in nature than in the laboratory. We finish with a multiscale description of a seminal structure in the Swiss Alps, the Glarus thrust, which puzzled geologists for more than 100 years. Along the Glarus thrust, a km-scale package of rocks (nappe) has been pushed 40 km over its footwall as a solid rock body. The thrust itself is a m-wide ductile shear zone, while in turn the centre of the thrust shows a mm-cm wide central slip zone experiencing periodic extreme deformation akin to a stick-slip event. The m-wide creeping zone is consistent with the THM feedback length scale of solid mechanics, while the ultralocalised central slip zones is most likely a fluid dynamic instability.

Recursive pathways to marginal likelihood estimation with prior-sensitivity analysis

We investigate the utility to computational Bayesian analyses of a particular family of recursive marginal likelihood estimators characterized by the (equivalent) algorithms known as "biased sampling" or "reverse logistic regression" in the statistics literature and "the density of states" in physics. Through a pair of numerical examples (including mixture modeling of the well-known galaxy dataset) we highlight the remarkable diversity of sampling schemes amenable to such recursive normalization, as well as the notable efficiency of the resulting pseudo-mixture distributions for gauging prior-sensitivity in the Bayesian model selection context. Our key theoretical contributions are to introduce a novel heuristic ("thermodynamic integration via importance sampling") for qualifying the role of the bridging sequence in this procedure, and to reveal various connections between these recursive estimators and the nested sampling technique.