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 horizontal wells 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 well-collapse 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.

Numerical analyses for the material bifurcation in plane sheet under tension

The various patterns (shear banding, surface wrinkling and necking) of material bifurcation in plane sheet under tension are investigated in this paper by means of a numerical method. It is found that numerical analysis can provide better ground for searching for the lowest critical loads. The inhomogeneity caused by void damage and the nonuniformity in the stress distribution across sheet thickness are proved to have detrimental effects on the material bifurcation. Nevertheless, material stability can be promoted by any means of depressing void damage or alleviating stress, even locally across the thickness. Besides, the peculiar behaviour of material bifurcation under slight biaxiality state is demonstrated. Copyright (C) 1996 Elsevier Science Ltd

Numerical study of bluff body flow structures

Our recent progress in numerical studies of bluff body flow structures and a new method for the numerical analysis of near wake flow field for high Reynolds number flow are introduced. The paper consists of three parts. In part one, the evolution of wake vortex structure and variation of forces on a flat plate in harmonic oscillatory flows and in in-line steady-harmonic combined flows are presented by an improved discrete vortex method, as the Keulegan-Carpenter number (KC) varies from 2 to 40 and ratios of U-m to U-0 are of O(10(-1)), O(10) and O(10), respectively. In part 2, a domain decomposition hybrid method, combining the finite-difference and vortex methods for numerical simulation of unsteady viscous separated flow around a bluff body, is introduced. By the new method, some high resolution numerical visualization on near wake evolution behind a circular cylinder at Re = 10(2), 10(3) and 3 x 10(3) are shown. In part 3, the mechanism and the dynamic process for the three-dimensional evolution of the Karman vortex and vortex filaments in braid regions as well as the early features of turbulent structure in the wake behind a circular cylinder are presented numerically by the vortex dynamics method.

Numerical Study on Critical Wedge Angle of Cellular Detonation Reflections

The critical wedge angle (CWA) for the transition from regular reflection (RR) to Mach reflection (MR) of a cellular detonation wave is studied numerically by an improved space-time conservation element and solution element method together with a two-step chemical reaction model. The accuracy of that numerical way is verified by simulating cellular detonation reflections at a 19.3∘ wedge. The planar and cellular detonation reflections over 45∘–55∘ wedges are also simulated. When the cellular detonation wave is over a 50∘ wedge, numerical results show a new phenomenon that RR and MR occur alternately. The transition process between RR and MR is investigated with the local pressure contours. Numerical analysis shows that the cellular structure is the essential reason for the new phenomenon and the CWA of detonation reflection is not a certain angle but an angle range.

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.

The Deposition and Burning Characteristics During Slagging Co-Firing Coal and Wood: Modeling and Numerical Simulation

Numerical analysis was used to study the deposition and burning characteristics of combining co-combustion with slagging combustion technologies in this paper. The pyrolysis and burning kinetic models of different fuels were implanted into the WBSF-PCC2 (wall burning and slag flow in pulverized co-combustion) computation code, and then the slagging and co-combustion characteristicsespecially the wall burning mechanism of different solid fuels and their effects on the whole burning behavior in the cylindrical combustor at different mixing ratios under the condition of keeping the heat input samewere simulated numerically. The results showed that adding wood powder at 25% mass fraction can increase the temperature at the initial stage of combustion, which is helpful to utilize the front space of the combustor. Adding wood powder at a 25% mass fraction can increase the reaction rate at the initial combustion stage; also, the coal ignitability is improved, and the burnout efficiency is enhanced by about 5% of suspension and deposition particles, which is helpful for coal particles to burn entirely and for combustion devices to minimize their dimensions or sizes. The results also showed that adding wood powder at a proper ratio is helpful to keep the combustion stability, not only because of the enhancement for the burning characteristics, but also because the running slag layer structure can be changed more continuously, which is very important for avoiding the abnormal slag accumulation in the slagging combustor. The theoretic analysis in this paper proves that unification of co-combustion and slagging combustion technologies is feasible, though more comprehensive and rigorous research is needed.

Numerical simulation of quantum directional couplers

A numerical analysis of a quantum directional coupler based on Pi-shaped electron waveguides is presented with use of the scattering-matrix method. After the optimization of the device parameters, uniform output for the two output ports and high directivity are obtained within a wide range of the electron momenta. The electron transfer in the device is found more efficient than that in the previously proposed structures. The study of the shape-dependence of transmission for the device shows that the device structure with smooth boundaries exhibits a much better performance.

Numerical investigation of particle velocity distributions in aeolian sand transport

Particle velocity distribution in a blowing sand cloud is a reflection of saltation movement of many particles. Numerical analysis is performed for particle velocity distribution with a discrete particle model. The probability distributions of resultant particle velocity in the impact-entrainment process, particle horizontal and vertical velocities at different heights and the vertical velocity of ascending particles are analyzed. The probability distributions of resultant impact and lift-off velocities of saltating particles can be expressed by a log-normal function, and that of impact angle comply with an exponential function. The probability distribution of particle horizontal and vertical velocities at different heights shows a typical single-peak pattern. In the lower part of saltation layer, the particle horizontal velocity distribution is positively skewed. Further analysis shows that the probability density function of the vertical velocity of ascending particles is similar to the right-hand part of a normal distribution function, and a general equation is acquired for the probability density function of non-dimensional vertical velocity of ascending particles which is independent of diameter of saltating particles, wind strength and height. These distributions in the present numerical analysis are consistent with reported experimental results. The present investigation is important for understanding the saltation state in wind-blown sand movement. (C) 2009 Elsevier B.V. All rights reserved.

Convergence of numerical time-averaging and stationary measures via Poisson equations

Numerical approximation of the long time behavior of a stochastic di.erential equation (SDE) is considered. Error estimates for time-averaging estimators are obtained and then used to show that the stationary behavior of the numerical method converges to that of the SDE. The error analysis is based on using an associated Poisson equation for the underlying SDE. The main advantages of this approach are its simplicity and universality. It works equally well for a range of explicit and implicit schemes, including those with simple simulation of random variables, and for hypoelliptic SDEs. To simplify the exposition, we consider only the case where the state space of the SDE is a torus, and we study only smooth test functions. However, we anticipate that the approach can be applied more widely. An analogy between our approach and Stein's method is indicated. Some practical implications of the results are discussed. Copyright © by SIAM. Unauthorized reproduction of this article is prohibited.

Numerical Method for Cost-Weight Optimisation of Stringer-Skin Panels

The need to integrate cost into the early product definition process as an engineering parameter is addressed. The application studied is a fuselage panel that is typical for commercial transport regional jets. Consequently, a semi-empirical numerical analysis using reference data was coupled to model the structural integrity of thin-walled structures with regard to material failure and buckling: skin, stringer, flexural, and interrivet. The optimization process focuses on direct operating cost (DOC) as a function of acquisition cost and fuel burn. It was found that the ratio of acquisition cost to fuel burn was typically 4:3 and that there was a 10% improvement in the DOC for the minimal DOC condition over the minimal weight condition because of the manufacturing cost saving from having a reduced number of larger-area stringers and a slightly thicker skin than that preferred by the minimal weight condition. Also note that the minimal manufacturing cost condition was slightly better than the minimal weight condition, which highlights the key finding: The traditional minimal weight condition is a dated and suboptimal approach to airframe structural design.

Two-dimensional Banach spaces with polynomial numerical index zero

We study two-dimensional Banach spaces with polynomial numerical indices equal to zero.

Numerical study of the effects of reinforcement/matrix interphase on stress-strain behavior of YAl2 particle reinforced MgLiAl composites

For the potential influence produced by the reinforcement/matrix interphase in particle reinforced metal matrix composites (PMMCs), a unit cell model with transition interphase was proposed. Uniaxial tensile loading was simulated and the stress/strain behavior was predicted. The results show that a transition interphase with both appropriate strength and thickness could affect the failure mode, reduce the stress concentration, and enhance the maximum strain value of the composite.

Experimental and Numerical Evaluation of Thermal Performance of Steered Fibre Composite Laminates

For Variable Stiffness (VS) composites with steered curvilinear tow paths, the fiber orientation angle varies continuously throughout the laminate, and is not required to be straight, parallel and uniform within each ply as in conventional composite laminates. Hence, the thermal properties (conduction), as well as the structural stiffness and strength, vary as functions of location in the laminate, and the associated composite structure is often called a “variable stiffness” composite structure. The steered fibers lead not only to the alteration of mechanical load paths, but also to the alteration of thermal paths that may
result in favorable temperature distributions within the laminate and improve the laminate performance. Evaluation of VS laminate performance under thermal loading is the focus of this chapter. Thermal performance evaluations require experimental and numerical analysis of VS laminates under different processing and loading conditions. One of the advantages of using composite materials in many applications is the tailoring capability of the laminate,
not only during the design phase but also for manufacturing. Heat transfer through variable conduction and chemical reaction (degree of cure) occurring during manufacturing (curing) plays an important role in the final thermal and mechanical performance, and shape of composite structures.

Pull-out behaviour of steel grid soil reinforcement embedded in silty sand

This paper investigates the pull-out behaviour (particularly the bearing resistance) of a steel grid reinforcement embedded in silty sand using laboratory tests and numerical analyses. It is demonstrated that the various common analytical equations for calculating the bearing component of pull-out resistance give a wide range of calculated values, up to about 200% disparity. The disparity will increase further if the issue of whether to use the peak or critical state friction angle is brought in. Furthermore, these equations suggest that the bearing resistance factor, N, is only a function of soil friction angle which is not consistent with some design guidelines. In this investigation, a series of large scale laboratory pull-out tests under different test pressures were conducted. The test results unambiguously confirmed that the N factor is a function of test pressure. A modified equation for calculating N is also proposed. To have more in-depth understanding of the pull-out behaviour, the tests were modelled numerically. The input parameters for the numerical analysis were obtained from laboratory triaxial tests. The analysis results were compared with the experimental results. Good agreement between experimental and numerical results was achieved if the strain-softening behaviour from peak strength to critical state condition was captured by the soil model used. © 2013 Elsevier Ltd.