Polygonal grain‐based distinct element modeling for mechanical behavior of brittle rock

The microstructure of rock was numerically reproduced by a polygonal grain‐based model, and its mechanical behavior was examined by performing the uniaxial compression test and Brazilian tests via the Universal Distinct Element Code. The numerical results of the model demonstrated good agreement with the experimental results obtained with rock specimens in terms of the stress–strain behavior, strength characteristics, and brittle fracture phenomenon. An encouraging result is that the grain‐based model‐Universal Distinct Element Code model can reproduce a low ratio of tensile to compressive strength of 1/20 to 1/10 without the need for an additional process. This finding is ascribed to the fact that the geometrical features of polygons can effectively capture the effects of angularity, finite rotation, and interlocking of grains that exist in reality. A numerical methodology to monitor the evolution of micro‐cracks was developed, which enabled us to examine the progressive process of the failure and distinguish the contribution of tensile cracking to the process from that of shear cracking. From the observations of the micro‐cracking process in reference to the stress–strain relation, crack initiation stress, and crack damage stress, it can be concluded that the failure process of the model closely resembles the microscopic observations of rock. We also carried out a parametric study to examine the relationships between the microscopic properties and the macroscopic behavior of the model. Depending on the micro‐properties, the model exhibited a variety of responses to the external load in terms of the strength and deformation characteristics, the evolution of micro‐cracks, and the post‐peak behavior. Copyright © 2016 John Wiley & Sons, Ltd.     

Computational model coupling mode II discrete fracture propagation with continuum damage zone evolution

We propose a numerical method that couples a cohesive zone model (CZM) and a finite element‐based continuum damage mechanics (CDM) model. The CZM represents a mode II macro‐fracture, and CDM finite elements (FE) represent the damage zone of the CZM. The coupled CZM/CDM model can capture the flow of energy that takes place between the bulk material that forms the matrix and the macroscopic fracture surfaces. The CDM model, which does not account for micro‐crack interaction, is calibrated against triaxial compression tests performed on Bakken shale, so as to reproduce the stress/strain curve before the failure peak. Based on a comparison with Kachanov's micro‐mechanical model, we confirm that the critical micro‐crack density value equal to 0.3 reflects the point at which crack interaction cannot be neglected. The CZM is assigned a pure mode II cohesive law that accounts for the dependence of the shear strength and energy release rate on confining pressure. The cohesive shear strength of the CZM is calibrated by calculating the shear stress necessary to reach a CDM damage of 0.3 during a direct shear test. We find that the shear cohesive strength of the CZM depends linearly on the confining pressure. Triaxial compression tests are simulated, in which the shale sample is modeled as an FE CDM continuum that contains a predefined thin cohesive zone representing the idealized shear fracture plane. The shear energy release rate of the CZM is fitted in order to match to the post‐peak stress/strain curves obtained during experimental tests performed on Bakken shale. We find that the energy release rate depends linearly on the shear cohesive strength. We then use the calibrated shale rheology to simulate the propagation of a meter‐scale mode II fracture. Under low confining pressure, the macroscopic crack (CZM) and its damaged zone (CDM) propagate simultaneously (i.e., during the same loading increments). Under high confining pressure, the fracture propagates in slip‐friction, that is, the debonding of the cohesive zone alternates with the propagation of continuum damage. The computational method is applicable to a range of geological injection problems including hydraulic fracturing and fluid storage and should be further enhanced by the addition of mode I and mixed mode (I+II+III) propagation. Copyright © 2016 John Wiley & Sons, Ltd.     

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.     

Material selection for the gas containment system of a compressed natural gas carrier fleet

In this work, the possible exploitation of fiber-reinforced composites in the context of maritime transportation of compressed natural gas (CNG) is investigated. In addition to a more conventional steel configuration, two different fiber materials, carbon and glass, are considered as construction materials for pressure vessels (PVs) to be stored on board ships, with thickness optimized by FEM analysis.The considered scenario is represented by the transportation of CNG from an offshore well to a terminal on shore. Fleets of ships carrying CNG in pressure vessels manufactured with the investigated materials are generated by means of a ship synthesis model (SSM) software and compared on the basis of technical and economical indicators.The choice of the construction material influences considerably the weight of the PVs, which represent a major item of total ship weight and reflects directly on the general transport performances in terms of resistance, seakeeping and reliability in the service. On the other hand, capital as well as operating expenditures are considerably affected by the choice. When exploring the design space, the ship synthesis model is able, at a preliminary stage of the design, to account for the various technical and economical aspects, their implications and relationships. Results are presented of computations carried out in a specific case, identified by the annual gas production and other characteristics of the well terminal and a cruising route for the ships. The comparison is carried out on the basis of the cost per transported unit of gas and of the percentage of success in the transportation process. The computations show that the choice of the PV material has a key influence on the results in terms of optimal number, dimensions and speed of the ships.     

Modeling the 2D behavior of dry‐stone retaining walls by a fully discrete element method

The dry‐stone retaining walls (DSRW) have been tipped as a promising solution for sustainable development. However, before recently, their behavior is relatively obscure. In this study, discrete element method (DEM) approach was applied to simulate the plane strain failure of these walls. A commercial DEM package (PFC2D™) was used throughout this study. The authors used a fully discrete approach; thus, both the wall and the backfill were modeled as discrete elements. The methodology for obtaining the micromechanical parameters was discussed in detail; this includes the three mechanical sub‐systems of DSRWs: wall, backfill and interface. The models were loaded progressively until failure, and then the results were compared with the full‐scale experimental results where the walls were loaded, respectively, with hydrostatic load and backfill. Despite its complexity and its intensive calculation time, DEM model can then be used to validate a more simplified approach. Copyright © 2015 John Wiley & Sons, Ltd.     

A catalogue of major and trace element data for Icelandic Holocene silicic tephra layers

Tephra layers with Icelandic provenance have been identified across the North Atlantic region in terrestrial, lacustrine, marine and glacial environments. These tephra layers are used as marker horizons in tephrochronology including climate studies, archaeology and environmental change. The major element chemistries of 19 proximally deposited Holocene Icelandic silicic tephra layers confirm that individual volcanic systems have unique geochemical signatures and that eruptions from the same system can often be distinguished. In addition, glass trace element chemistry highlights subtle geochemical variations between tephra layers which appear to have identical major element chemistry and thus allows for the identification of some, if not all, tephra layers previously considered identical in composition. This paper catalogues the compositional variation between the widespread Holocene Icelandic silicic tephra deposits.     

Tensionless–frictionless interaction of flexible annular foundation with a transversely isotropic multi‐layered half‐space

A transversely isotropic multi‐layered half‐space, with axis of material symmetry perpendicular to the free surface, supports a flexible either annular or solid circle foundation. The contact area of the foundation and the half‐space is considered to be both frictionless and tensionless. The foundation is assumed to be affected by a vertical static axisymmetric load. Detailed analysis of the interaction of these two systems with different thickness of layers is the target of this paper. With the use of ring load Green's functions for both the foundation and the continuum half‐space, an integral equation accompanied with some inequalities is introduced to model the complex BVP. With the incorporation of ring‐shape FEM, we are capable of capturing both regular and singular solution smoothly. The validity of the combination of the analytical and numerical method is proved with comparing the results of this paper with a number of benchmark cases of both linear and nonlinear interaction of circular and annular foundation with half‐space. Some new illustrations are presented to portray the aspect of the anisotropy and layering of the half‐space. Copyright © 2014 John Wiley & Sons, Ltd.     

Modeling deepwater seabed steady-state thermal fields around buried pipeline including trenching and backfill effects

Deepwater pipelines are designed to transport mixtures of oil and gas, and their associated impurities at wellhead temperatures that can be in excess of 149 °C (∼300 °F or 422 K) while the external temperature maybe in the range of 5 °C (∼41 °F or 278 K). Depending on the circumstances these pipelines may be buried for physical protection or for additional thermal insulation using robotic trenching equipment. This results in a complex cut and backfill geometry in the seafloor in addition to altering the thermal properties of the backfill. A two-dimensional boundary element model was developed specifically to address to investigate the local steady-state thermal field in the near field of the pipeline. The model allows one to account for the complex geometries in the near field associated with this burial technique, site-specific multi-layered soil conditions and the seawater adjacent to the seafloor. A parametric study was preformed to evaluate effects of the thermal power loss, burial depth, pipe diameter and soil thermal conductivity on the thermal field in the near field of a buried pipeline. The numerical examples illustrate the influence of the backfill thermal property on the temperature at the pipe wall, that the pipe diameter controls the required output thermal power needed to maintain the desired pipe wall temperature, and the importance of pipeline burial depth on seabed temperature distribution above the pipeline.     

Critical Top Tension for Static Equilibrium Configuration of A Steel Catenary Riser

This paper aims to present the critical top tension for static equilibrium configurations of a steel catenary riser（SCR） by using the finite element method. The critical top tension is the minimum top tension that can maintain the equilibrium of the SCR. If the top tension is smaller than the critical value, the equilibrium of the SCR does not exist. If the top tension is larger than the critical value, there are two possible equilibrium configurations. These two configurations exhibit the nonlinear large displacement. The configuration with the smaller displacement is stable, while the one with larger displacement is unstable. The numerical results show that the increases in the riser＇s vertical distances, horizontal offsets, riser＇s weights, internal flow velocities, and current velocities increase the critical top tensions of the SCR. In addition, the parametric studies are also performed in order to investigate the limit states for the analysis and design of the SCR.