A coupled hydro-mechanical analysis for prediction of hydraulic fracture propagation in saturated porous media using EFG mesh-less method

The details of the Element Free Galerkin (EFG) method are presented with the method being applied to a study on hydraulic fracturing initiation and propagation process in a saturated porous medium using coupled hydro-mechanical numerical modelling. In this EFG method, interpolation (approximation) is based on nodes without using elements and hence an arbitrary discrete fracture path can be modelled.The numerical approach is based upon solving two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Displacement increment and pore water pressure increment are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system of equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on the penalty method. In order to model discrete fractures, the so-called diffraction method is used.Examples are presented and the results are compared to some closed-form solutions and FEM approximations in order to demonstrate the validity of the developed model and its capabilities. The model is able to take the anisotropy and inhomogeneity of the material into account. The applicability of the model is examined by simulating hydraulic fracture initiation and propagation process from a borehole by injection of fluid. The maximum tensile strength criterion and Mohr-Coulomb shear criterion are used for modelling tensile and shear fracture, respectively. The model successfully simulates the leak-off of fluid from the fracture into the surrounding material. The results indicate the importance of pore fluid pressure in the initiation and propagation pattern of fracture in saturated soils. © 2013 Elsevier Ltd.

Laminar propagation of lean premixed flames ignited in stratified mixture

Combustion in stratified mixtures is envisaged in practical energy systems such as direct-injection spark-ignited (DISI) car engines, gas turbines, for reducing CO2 and pollutant emissions while protecting their efficiency. The mixture gradients change the fundamental properties of the flame, especially by a difference in temperature and composition between the burnt gases and those of a flame consuming a homogeneous mixture. This paper presents an investigation of the properties of the flame propagating in a lean homogeneous mixture after ignition in a richer mixture according to the magnitude of the stratification. Three magnitudes of stratification are investigated. The local flame burning velocity is determined by an original PIV algorithm developed previously. The local equivalence ratio in the fresh gases is measured from anisole PLIF. From the simultaneous PIV-PLIF measurements, the flame burning velocities conditioned on the local stretch rate and equivalence ratio in fresh gases are measured. The flame propagating through the homogeneous lean mixture has properties depending on the ignition conditions in the stratified layer. The flame propagating in the lean mixture is back-supported longer for ignition under the richer condition. The change of stretch sensitivity and burning velocity of the flame in the lean mixture is measured over time for the three magnitudes of mixture stratification investigated. The ignition in richer mixtures compensates for the nonequidiffusion effect of lean propane flame and sustains its robustness to stretch. The flame propagation in the lean homogeneous mixture is enhanced by ignition in a richer stratified layer, as much by their robustness to stretch as by an increase in the flame speed or the burning velocity. The decay time of this influence of the stratification, called memory effect, is determined. © 2013 The Combustion Institute.

Wave propagation in a slowly-varying duct with mean vortical swirling flow

The propagation of unsteady disturbances in a slowlyvarying cylindrical duct carrying mean swirling flow is investigated using a multiple-scales technique. This is applicable to turbomachinery flow behind a rotor stage when the swirl and axial velocities are of the same order. The presence of mean vorticity couples acoustic and vorticity equations which produces an eigenvalue problem that is not self-adjoint unlike that for irrotational mean flow. In order to determine the amplitude variation along the duct, an adjoint solution for the coupled system of equations is derived. The solution breaks down where a mode changes from cut on to cut off. In this region the amplitude is governed by a form of Airy's equation, and the effect of swirl is to introduce a small shift in the origin of the Airy function away from the turning-point location. The variation of axial wavenumber and amplitude along the duct is calculated. In hard-walled ducts mean swirl is shown to produce much larger amplitude variation along the duct compared with a nonswirling flow. Mean swirl also has a large effect in ducts with finite-impedance walls which differs depending on whether modes are co-rotating with the swirl or counter rotating. © 2001 by A.J. Cooper, Published by the American Institute of Aeronautics and Astronautics, Inc.

Infection and propagation of lymphocystis virus isolated from the cultured flounder Paralichthys olivaceus in grass carp cell lines

The causative agent of lymphocystis disease that frequently occurs in cultured flounder Paralichthys olivaceus in China is lymphocystis virus (LV). In this study, 13 fish cell lines were tested for their susceptibility to LV. Of these, 2 cell lines derived from the freshwater grass carp Ctenopharyngodon idellus proved susceptible to the LV, and 1 cell line, GCO (grass carp ovary), was therefore used to replicate and propagate the virus. An obvious cytopathic effect (CPE) was first observed in cell monolayers at 1 d post-inoculation, and at 3 d this had extended to about 75% of the cell monolayer. However, no further CPE extension was observed after 4 d. Cytopathic characteristics induced by the LV were detected by Giemsa staining and fluorescence microscopic observation with Hoechst 33258 staining. The propagated virus particles were also observed by electron microscopy. Ultrastructure analysis revealed several distinct cellular changes, such as chromatin compaction and margination, vesicle formation, cell-surface convolution, nuclear fragmentation and the occurrence of characteristic 'blebs' and cell fusion. This study provides a detailed report of LV infection and propagation in a freshwater fish cell line, and presents direct electron microscopy evidence for propagation of the virus in infected cells. A possible process by which the CPEs are controlled is suggested.

Flexible Control of Light Propagation Using a Novel Photonic Crystal Hetero-Waveguide Based on Silicon-on-Insulator Slab

We demonstrate a photonic crystal hetero-waveguide based on silicon-on-insulator (SOI) slab, consisting of two serially connected width-reduced photonic crystal waveguides with different radii of the air holes adjacent to the waveguide. We show theoretically that the transmission window of the structure corresponds to the transmission range common to both waveguides and it is in inverse proportion to the discrepancy between the two waveguides. Also the group velocity of guided mode can be changed from low to high or high to low, depending on which port of the structure the signal is input from just in the same device, and the variation is proportional to the discrepancy between the two waveguides. Using this novel structure, we realize flexible control of transmission window and group velocity of guided mode simultaneously.

Finite element beam propagation method for analysis of plasmonic waveguide

The basic idea of the finite element beam propagation method (FE-BPM) is described. It is applied to calculate the fundamental mode of a channel plasmonic polariton (CPP) waveguide to confirm its validity. Both the field distribution and the effective index of the, fundamental mode are given by the method. The convergence speed shows the advantage and stability of this method. Then a plasmonic waveguide with a dielectric strip deposited on a metal substrate is investigated, and the group velocity is negative for the fundamental mode of this kind of waveguide. The numerical result shows that the power flow direction is reverse to that of phase velocity.

Tuning of plasmon propagation in two-dimensional electrons

It is shown theoretically that the propagation of plasmons can be tuned by an external electric field via spin-orbit interactions in a two-dimensional electron gas formed in a semiconductor heterostructure. This may provide a manageable way of transmitting quantum information in a quantum device. A possible plasmon field effect transistor is proposed.

Single-mode behaviour judgment of optical waveguides by imaginary-distance beam propagation method under perfectly matched layer boundary condition

Imaginary-distance beam propagation method under the perfectly matched layer boundary condition is applied to judge single-mode behaviour of optical waveguides, for the first time to our knowledge. A new kind of silicon-on-insulator-based rib structures with half-circle cross-section is presented. The single-mode behaviour of this kind of waveguide with radius 2mum is investigated by this method. It is single-mode when the slab height is not smaller than the radius.

Optical waveguide propagation loss measurement using multiple reflections method

Optical waveguide propagation loss measurement method based on optical multiple reflections detection is presented in this paper. By using a precision reflectometer, uncertain influence on waveguide propagation loss measurement caused by fiber-waveguide coupling can be eliminated effectively and the waveguide net propagation loss can be measured accurately. To demonstrate this, the propagation loss of a Silicon-on-Insulator (SOI) rib waveguide fabricated by RIE is measured with the obtained value being 4.3 dB/cm. This method provides a non-destructive means for evaluating waveguide propagation loss. (c) 2005 Elsevier B.V. All rights reserved.

Calculation of propagation loss in photonic crystal waveguides by FDTD technique and Pade approximation

The propagation losses in single-line defect waveguides in a two-dimensional (2D) square-lattice photonic crystal (PC) consisted of infinite dielectric rods and a triangular-lattice photonic crystal slab with air holes are studied by finite-difference time-domain (FDTD) technique and a Pade approximation. The decaying constant beta of the fundamental guided mode is calculated from the mode frequency, the quality factor (Q-factor) and the group velocity v(g) as beta = omega/(2Qv(g)). In the 2D square-lattice photonic crystal waveguide (PCW), the decaying rate ranged from 10(3) to 10(-4) cm(-1) can be reliably obtained from 8 x 10(3)-item FDTD output with the FDTD computing time of 0.386 ps. And at most 1 ps is required for the mode with the Q-factor of 4 x 10(11) and the decaying rate of 10(-7) cm(-1). In the triangular-lattice photonic crystal slab, a 10(4)-item FDTD output is required to obtain a reliable spectrum with the Q-factor of 2.5 x 10(8) and the decaying rate of 0.05 cm(-1). (c) 2004 Elsevier B.V. All rights reserved.

Light-propagation characteristics of photonic crystal waveguide based on SOI materials at different polarized states

Straight single-line defect optical waveguides in photonic crystal slabs are designed by the plane wave expansion method and fabricated into silicon-on-insulator (SOI) wafer by 248-nm deep UV lithography. We present an efficient way to measure the light transmission spectrum of the photonic crystal waveguide (PhC WG) at given polarization states. By employing the Mueller/Stokes method, we measure and analyse the light propagation properties of the PhC WG at different polarized states. It is shown that experimental results are in agreement with the simulation results of the three-dimensional finite-difference-time-domain method.

A highly efficient beam propagation method for modeling step-index waveguides with tilt interface

Based on a new finite-difference scheme and Runge-Kutta method together with transparent boundary conditions (TBCs), a novel beam propagation method to model step-index waveguides with tilt interfaces is presented. The modified scheme provides an precies description of the tilt interface of the nonrectangular waveguide structure, showing a much better efficiency and accuracy comparing with the previously presented formulas.

A new beam propagation method based on least-squares expansion approximating

A new finite difference wide-angle beam propagation method is developed by introducing the least-squares expansion approximant in the propagator expansion. In this new method it is not necessary to select the reference index point because of the whole region approaching the lease-square expansion. This method avoids the problems induced by error selection of the reference index in the old methods based on Taylor or Pade expansion. Several typical structures are simulated by the new method and the results prove the validity of it.

In-plane propagation of cavity polaritons in a quantum semiconductor microcavity

Considering that the coupling among the heavy-hole exciton, light-hole exciton and the cavity photon can form bipolaritons in a quantum semiconductor microcavity, we calculate the group velocities of the cavity polaritons at different incident angles using the coupling model of three harmonic oscillators. The result indicates that the group velocities of the low and middle branches of the cavity polaritons have extrema, but the group velocities of the high branch increase with the increasing incident angle.