A model for moisture and heat flow in fractured unsaturated porous media

The present study investigates propagation of a cohesive crack in non‐isothermal unsaturated porous medium under mode I conditions. Basic points of skeleton deformation, moisture, and heat transfer for unsaturated porous medium are presented. Boundary conditions on the crack surface that consist of mechanical interaction of the crack and the porous medium, water, and heat flows through the crack are taken into consideration. For spatial discretization, the extended finite element method is used. This method uses enriched shape functions in addition to ordinary shape functions for approximation of displacement, pressure, and temperature fields. The Heaviside step function and the distance function are exploited as enrichment functions for representing the crack surfaces displacement and the discontinuous vertical gradients of the pressure and temperature fields along the crack, respectively. For temporal discretization, backward finite difference scheme is applied. Problems solved from the literature show the validity of the model as well as the dependency of structural response on the material properties and loading. Copyright © 2016 John Wiley & Sons, Ltd.     

An adaptive finite element model for wellbore breakout simulation

A finite element formulation is proposed and implemented for analysing the stability of excavated wells using the DiMaggio-Sandler constitutive elastoplastic model with a typical carbonate reservoir configuration. The quality of the finite element approximation is ensured by applying smooth curved elements adapted to the wellbore geometry, and h − p adaptive finite element meshes in the plastic zone. General purpose procedures are defined to transfer the elastoplastic deformation history to newly created integration points. A breakout damage criterion is proposed based on the second invariant of the deviatoric plastic deformation tensor. This damage criterion is used to apply a mesh movement algorithm to represent material collapse. The automatic successive application of the breakout damage criterion results in elliptical realistically looking geometries obtained in experiments reported in the literature.     

The stiff upper LIP: investigating the High Arctic Large Igneous Province

The Canadian Arctic Islands expose a complex network of dykes and sills that belong to the High Arctic Large Igneous Province (HALIP), which intruded volatile‐rich sedimentary rocks of the Sverdrup Basin (shale, limestone, sandstone and evaporite) some 130 to 120 million years ago. There is thus great potential in studying the HALIP to learn how volatile‐rich sedimentary rocks respond to magmatic heating events during LIP emplacement. The HALIP remains, however, one of the least well known LIPs on the planet due to its remote location, short field season, and harsh climate. A Canadian–Swedish team of geologists set out in summer 2015 to further explore HALIP sills and their sedimentary host rocks, including the sampling of igneous and meta‐sedimentary rocks for subsequent geochemical analysis, and high pressure‐temperature petrological experiments to help define the actual processes and time‐scales of magma–sediment interaction. The research results will advance our understanding of how climate‐active volatiles such as CO_2, SO₂ and CH₄ are mobilised during the magma–sediment interaction related to LIP events, a process which is hypothesised to have drastically affected Earth's carbon and sulphur cycles. In addition, assimilation of sulphate evaporites, for example, is anticipated to trigger sulphide immiscibility in the magma bodies and in so doing could promote the formation of Ni‐PGE ore bodies. Here we document the joys and challenges of ‘frontier arctic fieldwork’ and discuss some of our initial observations from the High Arctic Large Igneous Province.     

Evaluation of NLDAS‐2 evapotranspiration against tower flux site observations

The North American Land Data Assimilation System project phase 2 (NLDAS‐2) has run four land surface models for a 30‐year (1979–2008) retrospective period. Land surface evapotranspiration (ET) is one of the most important model outputs from NLDAS‐2 for investigating land–atmosphere interaction or to monitor agricultural drought. Here, we evaluate hourly ET using in situ observations over the Southern Great Plains (Atmospheric Radiation Measurement/Cloud and Radiation Testbed network) for 1 January 1997–30 September 1999 and daily ET u‐sing in situ observations at the AmeriFlux network over the conterminous USA for an 8‐year period (2000–2007). The NLDAS‐2 models compare well against observations, with the National Centers for Environmental Prediction's Noah land surface model performing best, followed, in order, by the Variable Infiltration Capacity, Sacramento Soil Moisture Accounting, and Mosaic models. Daily evaluation across the AmeriFlux network shows that for all models, performance depends on season and vegetation type; they do better in spring and fall than in winter or summer and better for deciduous broadleaf forest and grasslands than for croplands or evergreen needleleaf forest. Copyright © 2014 John Wiley & Sons, Ltd.     

Stability analysis of tunnels driven in jointed rocks by means of a homogenized limit analysis approach

Conceived as a potential alternative to the classical design methods employed for analyzing the stability of underground structures driven in jointed rocks, the homogenization approach stems from the heuristic idea that, from a macroscopic point of view, a rock mass cut by a network of joints may be perceived as a homogenized continuum. The strength properties of the latter can be theoretically obtained from the failure conditions of its individual constituents: rock matrix and joint interfaces. At the material level, the limit analysis reasoning is used in the context of homogenization to formulate the homogenized strength criterion of a jointed rock mass in the particular situation of a single set of parallel joints. As it could be expected, the obtained closed‐form expressions show the strength anisotropy induced by joint preferential orientation. The support functions (π functions) associated with the homogenized strength criterion are also determined in both plane strain and three‐dimensional cases. This criterion is then applied to the investigation of stability analysis of a tunnel excavated in a jointed rock mass. Upper bounds estimated of the stability factor are derived from the implementation of the kinematic approach directly on the homogenized underground structure. Finally, the approach is applied to analyze and discuss the collapse of the Pinheiros subway station (São Paulo, Brazil). Copyright © 2014 John Wiley & Sons, Ltd.