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

Optimisation of Organic Rankine Cycle (ORCs) for binary-cycle geothermal applications could play a major role in determining 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: the selection of working fluid and appropriate operating conditions as well as optimisation of the turbine design for those conditions will determine the amount of power that can be extracted from a resource. In this paper, we present the rationale for the use of radial-inflow turbines for ORC applications and the preliminary design of several radial-inflow machines based on a number of promising ORC systems that use five different working fluids: R134a, R143a, R236fa, R245fa and n-Pentane. Preliminary meanline analysis lead 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, based on R134a, was found to produce 33% more net power from a 150 °C resource flowing at 10 kg/s than the lowest performing cycle, based on n-Pentane.

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).

On the effective hydraulic conductivity and macrodispersivity for density-dependent groundwater flow

In this paper, semi-analytical expressions of the effective hydraulic conductivity ( KE) and macrodispersivity ( αE) for 3D steady-state density-dependent groundwater flow are derived using a stationary spectral method. Based on the derived expressions, we present the dependence of KE and αE on the density of fluid under different dispersivity and spatial correlation scale of hydraulic conductivity. The results show that the horizontal KE and αE are not affected by density-induced flow. However, due to gravitational instability of the fluid induced by density contrasts, both vertical KE and αE are found to be reduced slightly when the density factor ( γ ) is less than 0.01, whereas significant decreases occur when γ exceeds 0.01. Of note, the variation of KE and αE is more significant when local dispersivity is small and the correlation scale of hydraulic conductivity is large.

Determination of vertical hydraulic conductivity of aquitards in a multilayered leaky system using water-level signals in adjacent aquifers

This paper presents a methodology for determining the vertical hydraulic conductivity (Kv) of an aquitard, in a multilayered leaky system, based on the harmonic analysis of arbitrary water-level fluctuations in aquifers. As a result, Kv of the aquitard is expressed as a function of the phase-shift of water-level signals measured in the two adjacent aquifers. Based on this expression, we propose a robust method to calculate Kv by employing linear regression analysis of logarithm transformed frequencies and phases. The frequencies, where the Kv are calculated, are identified by coherence analysis. The proposed methods are validated by a synthetic case study and are then applied to the Westbourne and Birkhead aquitards, which form part of a five-layered leaky system in the Eromanga Basin, Australia.

Laboratory measurement of hydraulic conductivity functions of two unsaturated sandy soils during drying and wetting processes

The importance of applying unsaturated soil mechanics to geotechnical engineering design has been well understood. However, the consumption of time and the necessity for a specific laboratory testing apparatus when measuring unsaturated soil properties have limited the application of unsaturated soil mechanics theories in practice. Although methods for predicting unsaturated soil properties have been developed, the verification of these methods for a wide range of soil types is required in order to increase the confidence of practicing engineers in using these methods. In this study, a new permeameter was developed to measure the hydraulic conductivity of unsaturated soils using the steady-state method and directly measured suction (negative pore-water pressure) values. The apparatus is instrumented with two tensiometers for the direct measurement of suction during the tests. The apparatus can be used to obtain the hydraulic conductivity function of sandy soil over a low suction range (0-10 kPa). Firstly, the repeatability of the unsaturated hydraulic conductivity measurement, using the new permeameter, was verified by conducting tests on two identical sandy soil specimens and obtaining similar results. The hydraulic conductivity functions of the two sandy soils were then measured during the drying and wetting processes of the soils. A significant hysteresis was observed when the hydraulic conductivity was plotted against the suction. However, the hysteresis effects were not apparent when the conductivity was plotted against the volumetric water content. Furthermore, the measured unsaturated hydraulic conductivity functions were compared with predictions using three different predictive methods that are widely incorporated into numerical software. The results suggest that these predictive methods are capable of capturing the measured behavior with reasonable agreement.

Role of hydraulic factors in constructed wetland and bioretention basin treatment performance

This research project contributed to the in-depth understanding of the influence of hydrologic and hydraulic factors on the stormwater treatment performance of constructed wetlands and bioretention basins in the "real world". The project was based on the comprehensive monitoring of a Water Sensitive Urban Design treatment train in the field and underpinned by complex multivariate statistical analysis. The project outcomes revealed that the reduction in pollutant concentrations were consistent in the constructed wetland, but was highly variable in the bioretention basin to a range of influential factors. However, due to the significant amount retention within the filter media, all pollutant loadings were reduced in the bioretention basin.

A multi-scaling approach to predict hydraulic damage of poromaterials

The purpose of this paper is to introduce the concept of hydraulic damage and its numerical integration. Unlike the common phenomenological continuum damage mechanics approaches, the procedure introduced in this paper relies on mature concepts of homogenization, linear fracture mechanics, and thermodynamics. The model is applied to the problem of fault reactivation within resource reservoirs. The results show that propagation of weaknesses is highly driven by the contrasts of properties in porous media. In particular, it is affected by the fracture toughness of host rocks. Hydraulic damage is diffused when it takes place within extended geological units and localized at interfaces and faults.

Analysis and modelling of the hydraulic conductivity in aquitards : application to the Galilee Basin and the Great Artesian Basin, Australia

In large sedimentary basins with layers of different rocks, the groundwater flow between aquifers depends on the hydraulic conductivity (K) of the separating low-permeable rocks, or aquitards. Three methods were developed to evaluate K in aquitards for areas with limited field data: • Coherence and harmonic analysis: estimates the regional-scale K based on water-level fluctuations in adjacent aquifers. • Cokriging and Bayes' rule: infers K from downhole geophysical logs. • Fluvial process model: reproduces the lithology architecture of sediment formations which can be converted to K. These proposed methods enable good estimates of K and better planning of further drillholes.

Optimisation methods for coupled thermodynamic and 1D design of radial-inflow turbines

Optimisation is a fundamental step in the turbine design process, especially in the development of non-classical designs of radial-inflow turbines working with high-density fluids in low-temperature Organic Rankine Cycles (ORCs). The present work discusses the simultaneous optimisation of the thermodynamic cycle and the one-dimensional design of radial-inflow turbines. In particular, the work describes the integration between a 1D meanline preliminary design code adapted to real gases and the performance estimation approach for radial-inflow turbines in an established ORC cycle analysis procedure. The optimisation approach is split in two distinct loops; the inner operates on the 1D design based on the parameters received from the outer loop, which optimises the thermodynamic cycle. The method uses parameters including brine flow rate, temperature and working fluid, shifting assumptions such as head and flow coefficients into the optimisation routine. The discussed design and optimisation method is then validated against published benchmark cases. Finally, using the same conditions, the coupled optimisation procedure is extended to the preliminary design of a radial-inflow turbine with R143a as working fluid in realistic geothermal conditions and compared against results from commercially-available software RITAL from Concepts-NREC.

Hydraulic conductivity fields: Gaussian or not?

Hydraulic conductivity (K) fields are used to parameterize groundwater flow and transport models. Numerical simulations require a detailed representation of the K field, synthesized to interpolate between available data. Several recent studies introduced high-resolution K data (HRK) at the Macro Dispersion Experiment (MADE) site, and used ground-penetrating radar (GPR) to delineate the main structural features of the aquifer. This paper describes a statistical analysis of these data, and the implications for K field modeling in alluvial aquifers. Two striking observations have emerged from this analysis. The first is that a simple fractional difference filter can have a profound effect on data histograms, organizing non-Gaussian ln K data into a coherent distribution. The second is that using GPR facies allows us to reproduce the significantly non-Gaussian shape seen in real HRK data profiles, using a simulated Gaussian ln K field in each facies. This illuminates a current controversy in the literature, between those who favor Gaussian ln K models, and those who observe non-Gaussian ln K fields. Both camps are correct, but at different scales.

Hydrostratigraphic analysis of the MADE site with full-resolution GPR and direct-push hydraulic profiling

Full-resolution 3D Ground-Penetrating Radar (GPR) data were combined with high-resolution hydraulic conductivity (K) data from vertical Direct-Push (DP) profiles to characterize a portion of the highly heterogeneous MAcro Dispersion Experiment (MADE) site. This is an important first step to better understand the influence of aquifer heterogeneities on observed anomalous transport. Statistical evaluation of DP data indicates non-normal distributions that have much higher similarity within each GPR facies than between facies. The analysis of GPR and DP data provides high-resolution estimates of the 3D geometry of hydrostratigraphic zones, which can then be populated with stochastic K fields. The lack of such estimates has been a significant limitation for testing and parameterizing a range of novel transport theories at sites where the traditional advection-dispersion model has proven inadequate.

Applying HYDRUS to flow in a sodic clay soil with solution composition–dependent hydraulic conductivity

Management of sodic soils under irrigation often requires application of chemical ameliorants to improve permeability combined with leaching of excess salts. Modeling irrigation, soil treatments, and leaching in these sodic soils requires a model that can adequately represent the physical and chemical changes in the soil associated with the amelioration process. While there are a number of models that simulate reactive solute transport, UNSATCHEM and HYDRUS-1D are currently the only models that also include an ability to simulate the impacts of soil chemistry on hydraulic conductivity. Previous researchers have successfully applied these models to simulate amelioration experiments on a sodic loam soil. To further gauge their applicability, we extended the previous work by comparing HYDRUS simulations of sodic soil amelioration with the results from recently published laboratory experiments on a more reactive, repacked sodic clay soil. The general trends observed in the laboratory experiments were able to be simulated using HYDRUS. Differences between measured and simulated results were attributed to the limited flexibility of the function that represents chemistry-dependent hydraulic conductivity in HYDRUS. While improvements in the function could be made, the present work indicates that HYDRUS-UNSATCHEM captures the key changes in soil hydraulic properties that occur during sodic clay soil amelioration and thus extends the findings of previous researchers studying sodic loams.

Hydraulic conductivity increases in a sodic clay soil in response to gypsum applications impacts of bulk density and cation exchange

Amelioration of sodic soils is commonly achieved by applying gypsum, which increases soil hydraulic conductivity by altering soil chemistry. The magnitude of hydraulic conductivity increases expected in response to gypsum applications depends on soil properties including clay content, clay mineralogy, and bulk density. The soil analyzed in this study was a kaolinite rich sodic clay soil from an irrigated area of the Lower Burdekin coastal floodplain in tropical North Queensland, Australia. The impact of gypsum amelioration was investigated by continuously leaching soil columns with a saturated gypsum solution, until the hydraulic conductivity and leachate chemistry stabilized. Extended leaching enabled the full impacts of electrolyte effects and cation exchange to be determined. For the columns packed to 1.4 g/cm3, exchangeable sodium concentrations were reduced from 5.0 ± 0.5 mEq/100 g to 0.41 ± 0.06 mEq/100 g, exchangeable magnesium concentrations were reduced from 13.9 ± 0.3 mEq/100 g to 4.3 ± 2.12 mEq/100 g, and hydraulic conductivity increased to 0.15 ± 0.04 cm/d. For the columns packed to 1.3 g/cm3, exchangeable sodium concentrations were reduced from 5.0 ± 0.5 mEq/100 g to 0.51 ± 0.03 mEq/100 g, exchangeable magnesium concentrations were reduced from 13.9 ± 0.3 mEq/100 g to 0.55 ± 0.36 mEq/100 g, and hydraulic conductivity increased to 0.96 ± 0.53 cm/d. The results of this study highlight that both sodium and magnesium need to be taken into account when determining the suitability of water quality for irrigation of sodic soils and that soil bulk density plays a major role in controlling the extent of reclamation that can be achieved using gypsum applications.

Impacts of sodic soil amelioration on hydraulic conductivity and deep drainage in the Lower Burdekin

An understanding of the influence of soil chemistry on soil hydraulic properties is of critical importance for the management of sodic soils under irrigation. The hydraulic conductivity of sodic soils has been shown to be affected by properties of the applied solution including pH (Suarez et al. 1984), sodicity and salt concentration (McNeal and Coleman 1966). The changes in soil hydraulic conductivity are the result of changes in the spacing between clay layers in response to changes in soil solution chemistry. While the importance o f soil chemistry in controlling hydraulic conductivity is known, the exact impacts of sodic soil amelioration on hydraulic conductivity and deep drainage at a given location are difficult to predict. This is because the relationships between soil chemical factors and hydraulic conductivity are soil specific and because local site specific factors also need to be considered to determine the actual impacts on deep drainage rates.

Speed-based diagnostics of aerodynamic and mass imbalance in large wind turbines

Any kind of imbalance in the operation of a wind turbine has adverse effect on the downstream torsional components as well as tower structure. It is crucial to detect imbalance at its very inception. The identification of the type of imbalance is also required so that appropriate measures of fault accommodation can be performed in the control system. In particular, it is important to distinguish between mass and aerodynamic imbalance. While the former is gradually caused by a structural anomaly (e.g. ice deposition, moisture accumulation inside blade), the latter is generally associated to a fault in the pitch control system. This paper proposes a technique for the detection and identification of imbalance fault in large scale wind turbines. Unlike most other existing method it requires only the rotor speed signal which is readily available in existing turbines. Signature frequencies have been proposed in this work to identify imbalance type based on their physical phenomenology. The performance of this technique has been evaluated by simulations using an existing benchmark model. The effectiveness of the proposed method has been confirmed by the simulation results.