Physical and numerical modeling of horizontal gas-hydrate wells

Methane hydrate, which is usually found under deep seabed or permafrost zones, is a potential energy resource for future years. Depressurization of horizontalwells bored in methane hydrate layer is considered as one possible method for hydrate dissociation and methane extraction from the hosting soil. Since hydrate is likely to behave as a bonding material to sandy soils, supported well construction is necessary to avoid wellcollapse due to the loss of the apparent cohesion during depressurization. This paper describes both physical and numerical modeling of such horizontal support wells. The experimental part involves depressurization of small well models in a large pressure cell, while the numerical part simulates the corresponding problem. While the experiment models simulate only gas saturated initial conditions, the numerical analysis simulates both gas-saturated and more realistic water-saturated conditions based on effective stress coupled flow-deformation formulation of these three phases. © 2006 Taylor & Francis Group, London.

Stress-strain response of hydrate-bearing sands: Numerical study using discrete element method simulations

Gas hydrate is a crystalline solid found within marine and subpermafrost sediments. While the presence of hydrates can have a profound effect on sediment properties, the stress-strain behavior of hydrate-bearing sediments is poorly understood due to inherent limitations in laboratory testing. In this study, we use numerical simulations to improve our understanding of the mechanical behavior of hydrate-bearing sands. The hydrate mass is simulated as either small randomly distributed bonded grains or as "ripened hydrate" forming patchy saturation, whereby sediment clusters with 100% pore-filled hydrate saturation are distributed within a hydrate-free sediment. Simulation results reveal that reduced sand porosity and higher hydrate saturation cause an increase in stiffness, strength, and dilative tendency, and the critical state line shifts toward higher void ratio and higher shear strength. In particular, the critical state friction angle increases in sands with patchy saturation, while the apparent cohesion is affected the most when the hydrate mass is distributed in pores. Sediments with patchy hydrate distribution exhibit a slightly lower strength than sediments with randomly distributed hydrate. Finally, hydrate dissociation under drained conditions leads to volume contraction and/or stress relaxation, and pronounced shear strains can develop if the hydrate-bearing sand is subjected to deviatoric loading during dissociation.

A matrix-calculation-based algorithm for numerical change propagation analysis

Engineering changes (ECs) are raised throughout the lifecycle of engineering products. A single change to one component produces knock-on effects on others necessitating additional changes. This change propagation significantly affects the development time and cost and determines the product's success. Predicting and managing such ECs is, thus, essential to companies. Some prediction tools model change propagation by algorithms, whereof a subgroup is numerical. Current numerical change propagation algorithms either do not account for the exclusion of cyclic propagation paths or are based on exhaustive searching methods. This paper presents a new matrix-calculation-based algorithm which can be applied directly to a numerical product model to analyze change propagation and support change prediction. The algorithm applies matrix multiplications on mutations of a given design structure matrix accounting for the exclusion of self-dependences and cyclic propagation paths and delivers the same results as the exhaustive search-based Trail Counting algorithm. Despite its factorial time complexity, the algorithm proves advantageous because of its straightforward matrix-based calculations which avoid exhaustive searching. Thereby, the algorithm can be implemented in established numerical programs such as Microsoft Excel which promise a wider application of the tools within and across companies along with better familiarity, usability, practicality, security, and robustness. © 1988-2012 IEEE.

Numerical investigation of the active role of the actin cytoskeleton in the compression resistance of cells.

Numerous in-vitro studies have established that cells react to their physical environment and to applied mechanical loading. However, the mechanisms underlying such phenomena are poorly understood. Previous modelling of cell compression considered the cell as a passive homogenous material, requiring an artificial increase in the stiffness of spread cells to replicate experimentally measured forces. In this study, we implement a fully 3D active constitutive formulation that predicts the distribution, remodelling, and contractile behaviour of the cytoskeleton. Simulations reveal that polarised and axisymmetric spread cells contain stress fibres which form dominant bundles that are stretched during compression. These dominant fibres exert tension; causing an increase in computed compression forces compared to round cells. In contrast, fewer stress fibres are computed for round cells and a lower resistance to compression is predicted. The effect of different levels of cellular contractility associated with different cell phenotypes is also investigated. Highly contractile cells form more dominant circumferential stress fibres and hence provide greater resistance to compression. Computed predictions correlate strongly with published experimentally observed trends of compression resistance as a function of cellular contractility and offer an insight into the link between cell geometry, stress fibre distribution and contractility, and cell deformability. Importantly, it is possible to capture the behaviour of both round and spread cells using a given, unchanged set of material parameters for each cell type. Finally, it is demonstrated that stress distributions in the cell cytoplasm and nucleus computed using the active formulation differ significantly from those computed using passive material models.

Experimental and numerical study of grout injections in silty soils

Compensation grouting is increasingly employed as a mitigation technique of settlements induced by tunnelling and its effectiveness both in clayey and sandy soils is reported in a wide number of case histories. However, the results are highly dependent on grout properties, injection characteristics and soil properties. An experimental study was conducted to investigate the parameters that control grout injections in silty soils. The results from one injection test in a large sample of silty soil show that the compensation efficiency, defined as the ratio of the volume of heave obtained at ground surface and the injected grout volume, is much lower than one and tends to decrease with time, while the initial volume of grout lost due to pressure filtration is small. Finally, results from finite elements back analyses of the laboratory test show that a good agreement with the experimental data can be obtained if the development of large strains is taken into account. © 2012 Taylor & Francis Group.

Numerical simulation of long-term twin-tunnel behaviour at St James's park

The twin-tunnel construction of the Jubilee Line Extension tunnels beneath St James's Park was simulated using coupled-consolidation finite-element analyses. The effect of defining different permeabilities for the final consolidation stage was investigated, and the performance of a fissure softening model was also evaluated. The analyses suggested an unexpectedly high permeability anisotropy for soil around the tunnel crown, possibly due to stress-induced permeability changes, or low-permeability laminations. Also, the permeability profile and lining conductivity were found to differ between the tunnels. Inclusion of the fissure model gave a narrower settlement trough, more alike that in the field, by preferentially softening simple shear behaviour. Long-term settlements at the site continue to increase at an unexpectedly high rate, suggesting the possibility of creep or unexpected soil softening during excavation. © 2012 Taylor & Francis Group.

Progress towards the design and numerical analysis of a 3D microchannel biochip separator

This paper reports the design and numerical analysis of a three-dimensional biochip plasma blood separator using computational fluid dynamics techniques. Based on the initial configuration of a two-dimensional (2D) separator, five three-dimensional (3D) microchannel biochip designs are categorically developed through axial and plenary symmetrical expansions. These include the geometric variations of three types of the branch side channels (circular, rectangular, disc) and two types of the main channel (solid and concentric). Ignoring the initial transient behaviour and assuming that steady-state flow has been established, the behaviour of the blood fluid in the devices is algebraically analysed and numerically modelled. The roles of the relevant microchannel mechanisms, i.e. bifurcation, constriction and bending channel, on promoting the separation process are analysed based on modelling results. The differences among the different 3D implementations are compared and discussed. The advantages of 3D over 2D separator in increasing separation volume and effectively depleting cell-free layer fluid from the whole cross section circumference are addressed and illustrated. © 2011 John Wiley & Sons, Ltd.