Thermo-mechanical modelling of rock-like materials at very high temperature : application to ceramic refractories

Rock-like materials like ceramic refractories, in working conditions may be subject to large temperature variations. To simulate practical applications, bespoke constitutive modelling is required. In this work a general, thermodynamically consistent framework, able to incorporate key micromechanical features of the material behaviour, and applicable to a wide range of geomaterials, is formulated and validated. Different thermodynamic potentials are proposed to deal with both reversibility and irreversibility. A key advantage of this approach is the ability to freely choose the thermal dependency interpolation functions. Extensive model validation is provided by correctly reproducing both reversible and irreversible experimental trends of different materials under different loading conditions. It is found that even for simple materials, if a sample is subject to a large stress level, its thermal and mechanical responses become unexpectedly coupled. The proposed modelling framework is not limited to refractories and can be easily adapted to different types of rock-like materials

A computational framework for immiscible three-phase flow in deformable porous media

Several soil decontamination processes and enhanced oil recovery techniques involve the co-existence of three immiscible fluids, such as water, a nonaqueous phase liquid and a gas. In this work, a computational framework based on the individual mass balance of each phase is developed, aimed at simulating three-phase flow in a deformable rock through the finite element method, without resorting to specific simplifications that are usually required by standard numerical schemes. Key ingredients of the model are: expression of the residual in terms of mass contents, consistent lumping of the storage terms in the residual and algorithmic (tangent) matrix, consistent integration rules, the use of a minimum relative permeability and a time marching scheme based on trapezoidal integration. Special convective boundary conditions are adopted for pressures to be consistent with the assumed rock wettability properties during co-current imbibition. The resulting numerical scheme can deal with arbitrary saturation and/or pressure boundary conditions. The model is tested by simulating gas injection tests, and both co- and counter-current water imbibition tests, in a deformable core. To assess the performance and robustness of the whole framework, sensitivity analyses are performed upon varying key constitutive, loading and numerical parameters

Numerical analysis of thermal cycling during a multi-stage energy pile thermal response test

Energy piles are emerging as convenient alternative to the more traditional borehole heat exchangers (BHEs) to provide heating/cooling to buildings, as they remove the need for special purpose excavations and can accommodate more pipes, thus enhancing energy performance. However, their different aspect ratio compared to BHEs requires different modelling tools and dedicated thermal response testing, to achieve adequate thermal design. In this work, the results of an extended multi-stage thermal response test (TRT) carried out on a single energy pile installed in London Clay are presented in terms of both fluid temperature data and concrete temperature, measured by vibrating wire strain gauges and optic fibre sensors. The results are then explored in detail by means of a finite element numerical code, able to account for both convective heat exchange in the fluid, between the fluid and the solids and transient heat diffusion in the concrete and the ground. Analysis of the TRT field data shows that during the later stages of the test there is clear evidence of cyclic changes in performance. Investigation of these effects using the numerical model raises the possibility that there could be some alteration of the properties of the soil-pile contact during the test. Hypotheses for the observed behaviour are tentatively put forward and discussed with work recommended to further investigate the percieved phenomena

A new modelling approach for piled and other ground heat exchanger applications

Pile heat exchangers have an increasing role to play in the delivery of renewable heating and cooling energy. Traditionally the thermal design of ground heat exchangers has relied upon analytical approaches which take a relatively simple approach to the inside of the heat exchanger. This approach is justified while the heat exchanger diameter remains small. However, as larger diameter piled foundations are used as heat exchangers, the transient heat transfer processes operating within the pile become more important. To increase our understanding of these processes and ultimately lead to improved thermal design approaches for pile heat exchangers it is important to examine the heat transfer within the pile in detail. To accomplish this, a new numerical approach has been implemented within the finite element software ABAQUS. Coupling of the convective heat transfer due to fluid flow within the heat transfer pipes and the heat transfer by conduction within the pile concrete is the most important facet of the model. The resulting modelling approach, which is ready to generalise to other geothermal applications and to assess thermo-mechanical couplings, has been validated against a multi-stage thermal response test carried out on a test pile in London Clay

Dynamical effects during compaction band formation affecting their spatial periodicity

Compaction bands (CBs) are responsible for significant anisotropy alterations of permeability in geological materials; hence, understanding their formation conditions appears of key importance to all applications involving fluid extraction/injection from/into the ground. While most of the available models to understand CB formation are focused on interpreting the onset of a single CB, little effort has been so far dedicated to understand the documented periodicity of CBs. In this paper, the role of dynamical effects in inducing the post onset evolution of CBs is analyzed by means of a dedicated model for porous media with compressible constituents, with reference to a horizontal layer of sandy, water-saturated material. Elastic waves are generated as a first CB occurs due to sudden, localized volumetric collapse. If the waves are reflected at the interface with a softer material or with a previously formed CB, they produce significant local effective stress concentrations, which can promote the formation of further CBs in a cascade fashion, according to a regular geometric pattern. The spatial distribution of dynamically generated CBs, as well as the extent of the phenomenon, depends on the geometry of the domain and on the material's permeability. Sensitivity analysis is also performed to assess the key properties that promote dynamical CB in situ formation, identifying as the most influential conditions large stratum stiffness (increasing with depth) and the presence of softer layers. In contrast, the presence of less permeable and/or stiffer layers is not believed to play a major role in the proposed mechanism

Thermal performance of thermoactive continuous flight auger piles

Foundation piles are being increasingly equipped with heat exchangers to efficiently harvest shallow geothermal energy. For buildings in urban areas, continuous flight auger (CFA) piles are common owing to their speed, cost-efficiency and low noise levels. To construct a thermoactive CFA pile usually requires separate central installation of the heat exchanger. However, the energy performance of this type of pile has not been investigated systematically, with most studies focused on rotary piles where the heat exchanger is attached to the reinforcing cage. In this work, insights are provided about the main influences on the energy efficiency of thermoactive CFA piles, with a focus on the implications of using CFA construction techniques rather than rotary boring. An innovative three-dimensional numerical model, able to capture the different aspects of transient heat transfer, is employed together with analytical methods to evaluate the transient and steady-state behaviour of energy piles in a number of design situations. Attention is given to understanding the role of possible pipe-to-pipe interaction, which cannot be systematically investigated with standard methods. Finally, practical guidelines on the optimal choice of design parameters to maximise the energy efficiency of CFA piles, without altering the geotechnical arrangements, are provided

Evaluating the Applicability of the Radial Approximation for Pile Heat Exchangers

This paper appraises the efficacy of using an analytical radial approximation for different thermal pile heat exchanger geometries. Unsteady radial heat-flow from fluid in a pipe set within a grouted borehole into the external ground is well-documented and can be solved analytically very rapidly using Laplace Transforms (Javed and Claesson 2011). By comparing the radial model with finite-element simulations including explicit pile geometries, this paper provides a provisional analysis of the accuracy of this approach. Initial findings suggest that the radial model may provide an appropriate approximation to pile behaviour for certain pipe configurations, albeit with small ‘mid-time’ error

Energy performance of diaphragm walls used as heat exchangers

The possibility of equipping diaphragm walls as ground heat exchangers to meet the full or partial heating and cooling demands of overlying or adjacent buildings has been explored in recent years. In this paper, the factors affecting the energy performance of diaphragm walls equipped as heat exchangers are investigated through finite element modelling. The numerical approach employed is first validated using available experimental data and then applied to perform parametric analyses. Parameters considered in the analysis include panel width, the ratio between the wall and excavation depths, heat transfer pipe spacing, concrete cover, heat-carrier fluid velocity, concrete thermal properties and the temperature difference between the air within the excavation and the soil behind the wall. The results indicate that increasing the number of pipes by reducing their spacing is the primary route to increasing energy efficiency in the short term. However, the thermal properties of the wall concrete and the temperature excess within the excavation space are also important, with the latter becoming the most significant in the medium to long term. This confirms the benefits of exploiting the retaining walls installed for railway tunnels and metro stations where additional sources of heat are available

High Doses of Ursodeoxycholic Acid Up-Regulate the Expression of Placental Breast Cancer Resistance Protein in Patients Affected by Intrahepatic Cholestasis of Pregnancy

BACKGROUND: Ursodeoxycholic acid (UDCA) administration in intrahepatic cholestasis of pregnancy (ICP) induces bile acids (BA) efflux from the foetal compartment, but the molecular basis of this transplacental transport is only partially defined. AIM: To determine if placental breast cancer resistance protein (BCRP), able to transport BA, is regulated by UDCA in ICP. METHODS: 32 pregnant women with ICP (14 untreated, 34.9\ub15.17 years; 18 treated with UDCA--25 mg/Kg/day, 32.7\ub14.62 years,) and 12 healthy controls (33.4\ub13.32 years) agreed to participate in the study. Placentas were obtained at delivery and processed for membrane extraction. BCRP protein expression was evaluated by immunoblotting techniques and chemiluminescence quantified with a luminograph measuring emitted photons; mRNA expression with real time PCR. Statistical differences between groups were evaluated by ANOVA with Dunn's Multiple Comparison test. RESULTS: BCRP was expressed only on the apical membrane of the syncytiotrophoblast. A significant difference was observed among the three groups both for mRNA (ANOVA, p\u200a=\u200a0.0074) and protein (ANOVA, p<0.0001) expression. BCRP expression was similar in controls and in the untreated ICP group. UDCA induced a significant increase in placental BCRP mRNA and protein expression compared to controls (350.7\ub1106.3 vs 100\ub118.68% of controls, p<0.05 and 397.8\ub156.02 vs 100\ub111.44% of controls, p<0.001, respectively) and untreated ICP (90.29\ub117.59% of controls, p<0.05 and 155.0\ub113.87%, p<0.01). CONCLUSION: Our results confirm that BCRP is expressed only on the apical membrane of the syncytiotrophoblast and show that ICP treatment with high dose UDCA significantly upregulates placental BCRP expression favouring BA efflux from the foetal compartment