Numerical study on hydrate bearing sediments during gas production

The fully coupled methane hydrate model developed in Cambridge was adopted in this numerical study on gas production trial at the Eastern Nankai Trough, Japan 2013. Based on the latest experimental data of hydrate soil core samples, the clay parameters at Eastern Nankai site were successfully calibrated. With updated clay parameters and site geometry, a 50 days gas production trail was numerically simulated in FLAC2D. The geomechanical behaviour of hydrate bearing sediments under 3 different depressurization strategies were explored and discussed. The results from both axisymmetrical and plane-strain models suggest, the slope of the seabed only affects mechanical properties while no significant impact on the dissociation, temperature and pore pressure. For mechanical deformation after PT recovery, there are large settlements above the perforation zone and small uplift underneath the production zone. To validate the fully coupled model, numerical simulation with finer mesh in the hydrate production zone was carried out. The simulation results suggest good agreement between our model and JOE's results on history matching of gas and water production during trial. Parameter sensitivity of gas production is also investigated and concluded the sea water salinity is a dominant factor for gas production.

Experimental and numerical investigation of the seismic behavior of rectangular tunnels in soft soils

Underground structures constitute crucial components of the transportation networks. Considering their significance for modern societies, their proper seismic design is of great importance. However, this design may become very tricky, accounting of the lack of knowledge regarding their seismic behavior. Several issues that are significantly affecting this behavior (i.e. earth pressures on the structure, seismic shear stresses around the structure, complex deformation modes for rectangular structures during shaking etc.) are still open. The problem is wider for the non-circular (i.e. rectangular) structures, were the soilstructure interaction effects are expected to be maximized. The paper presents representative experimental results from a test case of a series of dynamic centrifuge tests that were performed on rectangular tunnels embedded in dry sand. The tests were carried out at the centrifuge facility of the University of Cambridge, within the Transnational Task of the SERIES EU research program. The presented test case is also numerically simulated and studied. Preliminary full dynamic time history analyses of the coupled soil-tunnel system are performed, using ABAQUS. Soil non-linearity and soil-structure interaction are modeled, following relevant specifications for underground structures and tunnels. Numerical predictions are compared to experimental results and discussed. Based on this comprehensive experimental and numerical study, the seismic behavior of rectangular embedded structures is better understood and modeled, consisting an important step in the development of appropriate specifications for the seismic design of rectangular shallow tunnels.

Local strain rate and curvature dependences of scalar dissipation rate transport in turbulent premixed flames: A direct numerical simulation analysis

The statistical behaviours of the instantaneous scalar dissipation rate Nc of reaction progress variable c in turbulent premixed flames have been analysed based on three-dimensional direct numerical simulation data of freely propagating statistically planar flame and V-flame configurations with different turbulent Reynolds number Ret. The statistical behaviours of N c and different terms of its transport equation for planar and V-flames are found to be qualitatively similar. The mean contribution of the density-variation term T1 is positive, whereas the molecular dissipation term (-D2) acts as a leading order sink. The mean contribution of the strain rate term T2 is predominantly negative for the cases considered here. The mean reaction rate contribution T3 is positive (negative) towards the unburned (burned) gas side of the flame, whereas the mean contribution of the diffusivity gradient term (D) assumes negative (positive) values towards the unburned (burned) gas side. The local statistical behaviours of Nc, T1, T2, T 3, (-D2), and f(D) have been analysed in terms of their marginal probability density functions (pdfs) and their joint pdfs with local tangential strain rate aT and curvature km. Detailed physical explanations have been provided for the observed behaviour. © 2014 Y. Gao et al.

3D numerical analysis of tunnelling induced damage: the influence of the alignment of a masonry building with the tunnel axis.

The development of infrastructure in major cities often involves tunnelling, which can cause damage to existing structures. Therefore, these projects require a careful prediction of the risk of settlement induced damage. The simplified approach of current methods cannot account for three-dimensional structural aspects of buildings, which can result in an inaccurate evaluation of damage. This paper investigates the effect of the building alignment with the tunnel axis on structural damage. A three-dimensional, phased, fully coupled finite element model with non-linear material properties is used as a tool to perform a parametric study. The model includes the simulation of the tunnel construction process, with the tunnel located adjacent to a masonry building. Three different type of settlements are included (sagging, hogging and a combination of them), with seven different increasing angles of the building with respect to the tunnel axis. The alignment parameter is assessed, based on the maximum occurring crack width, measured in the building. Results show a significant dependency of the final damage on the building and tunnel alignment.

Monitoring and preventing numerical oscillations in 3D simulations with coupled Monte Carlo codes

Previous studies have reported that different schemes for coupling Monte Carlo (MC) neutron transport with burnup and thermal hydraulic feedbacks may potentially be numerically unstable. This issue can be resolved by application of implicit methods, such as the stochastic implicit mid-point (SIMP) methods. In order to assure numerical stability, the new methods do require additional computational effort. The instability issue however, is problem-dependent and does not necessarily occur in all cases. Therefore, blind application of the unconditionally stable coupling schemes, and thus incurring extra computational costs, may not always be necessary. In this paper, we attempt to develop an intelligent diagnostic mechanism, which will monitor numerical stability of the calculations and, if necessary, switch from simple and fast coupling scheme to more computationally expensive but unconditionally stable one. To illustrate this diagnostic mechanism, we performed a coupled burnup and TH analysis of a single BWR fuel assembly. The results indicate that the developed algorithm can be easily implemented in any MC based code for monitoring of numerical instabilities. The proposed monitoring method has negligible impact on the calculation time even for realistic 3D multi-region full core calculations. © 2014 Elsevier Ltd. All rights reserved.

Numerical analysis of the settlement induced damage to Palazzo Loggia in Brescia

One of the main causes of failure of historic buildings is represented by the differential settlements of foundations. Finite element analysis provides a useful tool for predicting the consequences of given ground displacements in terms of structural damage and also assesses the need of strengthening techniques. The actual damage classification for buildings subject to settlement bases the assessment of the potential damage on the expected crack pattern of the structure. In this paper, the correlation between the physical description of the damage in terms of crack width and the interpretation of the finite element analysis output is analyzed. Different discrete and continuum crack models are applied to simulate an experiment carried on a scale model of a masonry historical building, the Loggia Palace in Brescia (Italy). Results are discussed and a modified version of the fixed total strain smeared crack model is evaluated, in order to solve the problem related to the calculation of the exact crack width.

Numerical analyses of tunnel-induced settlement damage to a masonry wall

Settlements due to underground construction represent a risk for the architectural heritage, especially in The Netherlands, because of the combination of soft soil, fragile pile foundation and brittle, un-reinforced masonry façade. Modelling of soil-structure interaction is fundamental to assess the risk of building damage due to tunnelling. This paper presents results of finite element analyses carried out with different models for a simple masonry wall. Focus is paid on the comparison between coupled, uncoupled and semi-coupled analyses, in which the soil-structure interaction is represented in different ways. In particular, the implementation of a soil-structure interface model in the numerical analyses is analysed, in order to asses its validity. The aim of the research project is the development of a damage classification system for different building typologies.