Thermo-hydro-geomechanical simulations of methane gas production from deep sea methane hydrate formations

A new constitutive model called Methane Hydrate Critical State (MHCS) model was conducted to investigate the geomechanical response of the gas-hydrate-bearing sediments at the Nankai Trough during the wellbore construction process. The strength and dilatancy of gas-hydrate-bearing soil would gradually disappear when the bonds are destroyed because of excessively shearing, which are often observed in dense soils and also in bonded soils such as cemented soil and unsaturated soil. In this study, the MHCS model, which presents such softening features, would be incorporated into a staged-finite-element model in ABAQUS, which mainly considered the loading history of soils and the interaction between cement-casing-formation. This model shows the influence of gas-hydrate-bearing soil to the deformation and stability of a wellbore and the surrounding sediments during wellbore construction. At the same time, the conventional Mohr-Coulomb model was used in the model to show the advantages of MHCS model by comparing the results of the two models.

Differential equation based length scales to improve DES and RANS simulations

Surprisingly expensive to compute wall distances are still used in a range of key turbulence and peripheral physics models. Potentially economical, accuracy improving differential equation based distance algorithms are considered. These involve elliptic Poisson and hyperbolic natured Eikonal equation approaches. Numerical issues relating to non-orthogonal curvilinear grid solution of the latter are addressed. Eikonal extension to a Hamilton-Jacobi (HJ) equation is discussed. Use of this extension to improve turbulence model accuracy and, along with the Eikonal, enhance Detached Eddy Simulation (DES) techniques is considered. Application of the distance approaches is studied for various geometries. These include a plane channel flow with a wire at the centre, a wing-flap system, a jet with co-flow and a supersonic double-delta configuration. Although less accurate than the Eikonal, Poisson method based flow solutions are extremely close to those using a search procedure. For a moving grid case the Poisson method is found especially efficient. Results show the Eikonal equation can be solved on highly stretched, non-orthogonal, curvilinear grids. A key accuracy aspect is that metrics must be upwinded in the propagating front direction. The HJ equation is found to have qualitative turbulence model improving properties. © 2003 by P. G. Tucker.

Simulations of autoignition and laminar premixed flames in methane/air mixtures diluted with hot products

This article considers constant-pressure autoignition and freely propagating premixed flames of cold methane/air mixtures mixed with equilibrium hot products at high enough dilution levels to burn within the moderate to intense low oxygen dilution (MILD) combustion regime. The analysis is meant to provide further insight on MILD regime boundaries and to identify the effect of hot products speciation. As the mass fraction of hot products in the reactants mixture increases, autoignition occurs earlier. Species profiles show that the products/reactants mixture approximately equilibrates to a new state over a quick transient well before the main autoignition event, but as dilution becomes very high, this equilibration transient becomes more prominent and eventually merges with the primary ignition event. The dilution level at which these two reactive zones merge corresponds well with that marking the transition into the MILD regime, as defined according to conventional criteria. Similarly, premixed flame simulations at high dilutions show evidence of significant reactions involving intermediate species prior to the flame front. Since the premixed flame governing equations system demands that the species and temperature gradients be zero at the "cold" boundary, flame speed cannot be calculated above a certain dilution level. Up to this point, which again agrees reasonably well with the transition into the MILD regime according to convention, the laminar burning velocity was found to increase with hot product dilution while flame thickness remained largely unchanged. Some comments on the MILD combustion regime boundary definition for gas turbine applications are included. Copyright © Taylor & Francis Group, LLC.

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.

New finite-difference scheme for simulations of step-index waveguides with tilt interfaces

A new finite-difference scheme is presented for the second derivative of a semivectorial field in a step-index optical waveguide with tilt interfaces. The present scheme provides an accurate description of the tilt interface of the nonrectangular structure. Comparison with previously presented formulas shows the effectiveness of the present scheme.

Research of Face Location System Based on Human Vision Simulations

In this paper we present a robust face location system based on human vision simulations to automatically locate faces in color static images. Our method is divided into four stages. In the first stage we use a gauss low-pass filter to remove the fine information of images, which is useless in the initial stage of human vision. During the second and the third stages, our technique approximately detects the image regions, which may contain faces. During the fourth stage, the existence of faces in the selected regions is verified. Having combined the advantages of Bottom-Up Feature Based Methods and Appearance-Based Methods, our algorithm performs well in various images, including those with highly complex backgrounds.

Numerical simulations of interrupted and conventional microchannel heat sinks

We provide three-dimensional numerical simulations of conjugate heat transfer in conventional and the newly proposed interrupted microchannel heat sinks. The new microchannel heat sink consists of a set of separated zones adjoining shortened parallel microchannels and transverse microchambers. Multi-channel effect, physical property variations, and axial thermal conduction are considered. It is found that flow rate variations in different channels can be neglected, while heat received by different channels accounts for 2% deviations from the averaged value when the heat flux at the back surface of the silicon chip reaches 100 W/cm(2). The computed hydraulic and thermal boundary layers are redeveloping in each separated zone due to shortened flow length for the interrupted microchannel heat sink. The periodic thermal developing flow is responsible for the significant heat transfer enhancement. Two effects influence pressure drops across the newly proposed microchannel heat sink. The first one is the pressure recovery effect in the microchamber, while the second one is the head loss when liquid leaves the microchamber and enters the next zone. The first effect compensates or suppresses the second one, leading to similar or decreased pressure drop than that for the conventional microchannel heat sink, with the fluid Prandtl number larger than unity.

Two-dimensional numerical simulations of shock waves in micro convergent-divergent nozzles

The steady two-dimensional Navier-Stokes equations with the slip wall boundary conditions were used to simulate the supersonic flow in micro convergent-divergent nozzles. It is observed that shock waves can take place inside or outside of the micronozzles under the earth environment. For the over-expanded flows, there is a boundary layer separation point, downstream of which a wave interface separates the viscous boundary layer with back air flow and the inviscid core flow. The oblique shock wave is followed by the bow shock and shock diamond. The viscous boundary layer thickness relative to the whole nozzle width on the exit plane is increased but attains the maximum value around of 0.5 and oscillates against this value with the continuous increasing of the nozzle upstream pressures. The viscous effect either changes the normal shock waves outside of the nozzle for the inviscid flow to the oblique shock waves inside the nozzle, or transfers the expansion jet flow without shock waves for the inviscid flow to the oblique shock waves outside of the nozzle. 

MANY-BEAM SIMULATIONS AND OBSERVATIONS OF LARGE-ANGLE CONVERGENT-BEAM ELECTRON-DIFFRACTION IMAGING OF CRYSTAL DEFECTS

Many-beam dynamical simulations and observations have been made for large-angle convergent-beam electron diffraction (LACBED) imaging of crystal defects, such as stacking faults and dislocations. The simulations are based on a general matrix formulation of dynamical electron diffraction theory by Peng and Whelan, and the results are compared with experimental LACBED images of stacking faults and dislocations of Si angle crystals. Excellent agreement is achieved.