On the use of quarter-point tetrahedral finite elements in linear elastic fracture mechanics

This paper discusses the reproduction of the square root singularity in quarter-point tetrahedral (QPT) finite elements. Numerical results confirm that the stress singularity is modeled accurately in a fully unstructured mesh by using QPTs. A displacement correlation (DC) scheme is proposed in combination with QPTs to compute stress intensity factors (SIF) from arbitrary meshes, yielding an average error of 2–3%. This straightforward method is computationally cheap and easy to implement. The results of an extensive parametric study also suggest the existence of an optimum mesh-dependent distance from the crack front at which the DC method computes the most accurate SIFs

A finite element framework for modeling internal frictional contact in three-dimensional fractured media using unstructured tetrahedral meshes

AbstractThis paper introduces a three-dimensional finite element (FE) formulation to accurately model the linear elastic deformation of fractured media under compressive loading. The presented method applies the classic Augmented Lagrangian(AL)-Uzawa method, to evaluate the growth of multiple interacting and intersecting discrete fractures. The volume and surfaces are discretized by unstructured quadratic triangle-tetrahedral meshes; quarter-point triangles and tetrahedra are placed around fracture tips. Frictional contact between crack faces for high contact precisions is modeled using isoparametric integration point-to-integration point contact discretization, and a gap-based augmentation procedure. Contact forces are updated by interpolating tractions over elements that are adjacent to fracture tips, and have boundaries that are excluded from the contact region. Stress intensity factors are computed numerically using the methods of displacement correlation and disk-shaped domain integral. A novel square-root singular variation of the penalty parameter near the crack front is proposed to accurately model the contact tractions near the crack front. Tractions and compressive stress intensity factors are validated against analytical solutions. Numerical examples of cubes containing one, two, twenty four and seventy interacting and intersecting fractures are presented

A disk-shaped domain integral method for the computation of stress intensity factors using tetrahedral meshes

A novel domain integral approach is introduced for the accurate computation of pointwise J-integral and stress intensity factors (SIFs) of 3D planar cracks using tetrahedral elements. This method is efficient and easy to implement, and does not require a structured mesh around the crack front. The method relies on the construction of virtual disk-shaped integral domains at points along the crack front, and the computation of domain integrals using a series of virtual triangular and line elements. The accuracy of the numerical results computed for through-the-thickness, penny-shaped, and elliptical crack configurations has been validated by using the available analytical formulations. The average error of computed SIFs remains below 1% for fine meshes, and 2–3% for coarse ones. The results of an extensive parametric study suggest that there exists an optimum mesh-dependent domain radius at which the computed SIFs are the most accurate. Furthermore, the results provide evidence that tetrahedral elements are efficient, reliable and robust instruments for accurate linear elastic fracture mechanics calculations

Simultaneous oil recovery and residual gas storage: A pore-level analysis using in-situ X-ray micro-tomography

We imaged sandstone cores at residual gas saturation (Sgr) with synchrotron radiation at a nominal resolution of (9 μm)3. We studied two three-phase flooding sequences: (1) gas injection into a core containing oil and initial water followed by a waterflood (gw process); (2) gas injection into a waterflooded core followed by another waterflood (wgw process). In the gw flood we measured a significantly higher Sgr (=20.6%; Sgr in the wgw flood was 5.3%) and a significantly lower residual oil saturation (Sor; Sor in the gw flood was 21.6% and Sor in the wgw flood was 29.3%). We also studied the size distribution of individual trapped clusters in the pore space. We found an approximately power-law distribution N ∝ s−τ with an exponent τ = 2.02–2.03 for the residual oil clusters and τ = 2.04 for the gas clusters in the gw flood. τ (=2.32) estimated for the gas clusters in the wgw process was significantly different. Furthermore, we calculated the surface area A–volume V relationships for the clusters. Again an approximate power-law relationship was observed, A ∝Vp with p ≈ 0.75. Moreover, in the gw flood sequence we identified oil layers sandwiched between the gas and water phases; we did not identify such oil layers in the wgw flood.These results have several important implications for oil recovery, carbon geo-sequestration and contaminant transport: (a) significantly more oil can be produced and much more gas can be stored using a gw flood; (b) cluster size distributions for residual oil or gas clusters in three-phase flow are similar to those observed in analogue two-phase flow; (c) there is a large cluster surface area available for dissolution of the residual phase into an aqueous phase; however, this surface area is significantly smaller than predicted by percolation theory (p ≈ 1), which implies that CO2 dissolution trapping and contamination of aquifers by hazardous organic solvents is slower than expected because of reduced interfacial contact areas

Extracting Information about the Rotator Cuff from Magnetic Resonance Images Using Deterministic and Random Techniques

We consider some methods to extract information about the rotator cuff based on magnetic resonance images; the study aims to define an alternative method of display that might facilitate the detection of partial tears in the supraspinatus tendon. Specifically, we are going to use families of ellipsoidal triangular patches to cover the humerus head near the affected area. These patches are going to be textured and displayed with the information of the magnetic resonance images using the trilinear interpolation technique. For the generation of points to texture each patch, we propose a new method that guarantees the uniform distribution of its points using a random statistical method. Its computational cost, defined as the average computing time to generate a fixed number of points, is significantly lower as compared with deterministic and other standard statistical techniques

CONTAIN D11 : integrated final results and conclusions

Carbon capture and storage is a technology capable of reducing CO2 outputs on a large scale; the concept usually requires CO2 to be removed from post-combustion flue gases and sequestered in geological formations. Depleted gas fields constitute “the most important storage type for the UK” and will provide a large and important potential future offshore storage capacity (DECC, 2012). Over the last 4 years, the CONTAIN research project has focussed on the geomechanical behaviour of depleted hydrocarbon fields in response to injection with CO2, combining a modelling and experimental approach with the public perceptions of CCS into three work packages. The project has provided a better understanding of the hydromechanical impacts of depletion on caprocks and the effect of subsequent CO2 injection, in order to assist with the implementation of CCS in this type of reservoir. Work package 1 outlined a phenomenological approach to assessing possible deformation during operation. Focus was placed on rock mechanics and transport experiments on material from the geologies of target formations in the North Sea, providing information that could be incorporated into numerical simulations. Work package 2 expanded this understanding by considering fractured caprock. Numerical modelling was used to study the deformation of an initially intact caprock caused by the depletion of an underlying reservoir during oil extraction. Deformation and flow were geomechanically modelled in three dimensions using a fully coupled poroelastic model, incorporating discrete fractures and faults into the caprock. Work package 3 offered new and valuable insight on future public awareness campaigns aimed at gaining acceptance of CCS. Qualitative expert interviews have been used, a CCS expert survey and a public survey across four countries to gain an understanding of perceptions of CCS risks and benefits, and has allowed for comparison of views on CCS between experts and public. In addition, the work package has explored the impact of different message framings on CCS attitudes. The findings of each work package are summarised in this report, with each work package represented by a report chapter. A synthesis of the findings and discussion of the work as a whole follows

Familial influences on the clinical characteristics of major depression: a twin study

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/66157/1/j.1600-0447.1992.tb03283.x.pd