Simulation of Natural Gas Production in Hydrate Reservoirs

This paper simulates a one-dimensional physical model of natural gas production from hydrate dissociation in a reservoir by depressurization. According to the principles of solid hydrate decomposition in stratum and flow of natural gas in porous medium, the pressure governing equations for both gas zone and hydrate zone are set up based on the physical production model. Using the approximation reported by N. N. Verigin et al. (1980), the nonlinear governing equations are simplified and the self-similar solutions are obtained. Through calculation, for different reservoir parameters, the distribution characters of pressure are analyzed. The decline trend of natural gas production rate with time is also studied. The simulation results show that production of natural gas from a hydrate reservoir is very sensitive to several reservoir parameters, such as wellbore pressure and stratum porosity and permeability.

Spatiotemporal evolution and multiple self-focusing of ultrashort pulses in a resonant two-level medium

The spatiotemporal evolutions of ultrashort pulses in two dimensions are investigated numerically by solving the coupled Maxwell-Bloch equations without invoking the slowly varying envelope approximation and rotating-wave approximation. For an on-axis 2n pi sech pulse, local delay makes the temporal split 2 pi sech pulses crescent-shaped in the transverse distribution. Due to the transverse effect, the temporal split 2 pi sech pulses become unstable and experience reshaping during the propagation process. Then, interference occurs between the successive crescent-shaped pulses and multiple self-focusing can form.

Quantum interference effects of two coherent population trapping states on the atomic spectral lines of a F-e=0 -> F-g=1 transition

We investigate the fluorescence spectrum in a nearly degenerate atomic system of a F-e = 0 -> F-g = 1 transition by analytically solving Schrodinger equations. An ultranarrow fluorescence spectral line in between the two coherent population trapping windows has been found. Our analytic solutions clearly show the origin of the ultranarrow spectral line. Due to quantum interference effects between two coherent population trapping states, the width and intensity of the central spectral line can be controlled by an external magnetic field. Such an effect may be used to detect a magnetic field.

体全息光栅透镜的设计和应用

设计了一种新型的体全息光栅透镜，在一块光学平板(体全息记录材料)内可以将输入光束产生横向传输并聚焦，或对输入光点产生横传的准直．它由一束平面波和一束球面波正交入射到光学平板上干涉形成的．研究了该体全息透镜的光栅间距变化情况，为设计和制备体全息光栅透镜及相关器件提供了理论依据．基于两光束耦合波理论，得到了该光栅透镜的耦合波方程，近似计算了该透镜的衍射效率及其达到高衍射效率时透镜的最佳尺寸．最后，讨论了该透镜在集成光学等领域中的应用．

Corrections to the paraxial approximation solutions in transversely nonuniform refractive-index media

The coupled differential recurrence equations for the corrections to the paraxial approximation solutions in transversely nonuniform refractive-index media are established in terms of the perturbation method. All the corrections (including the longitudinal field corrections) to the paraxial approximation solutions are presented in the weak-guidance approximation. As a concrete application, the first-order longitudinal field correction and the second-order transverse field correction to the paraxial approximation of a Gaussian beam propagating in a transversely quadratic refractive index medium are analytically investigated. (C) 1999 Optical Society of America [S0740-3232(99)00310-5].

色散管理孤子传输特性分析

色散管理孤子光波系统的色散图周期由具有正负色散值的两段光纤组成，并且在系统中周期性放置集中放大器及滤波器装置。采用变分法获得了描述此种系统性能的常微分方程组，并利用这组方程数值研究了各种参量对孤子传输性能的影响。研究结果表明，虽然从整体趋势上看，当要求色散管理孤子有较好的稳定性时，应适当降低系统的色散管理强度，但在一定条件下，调整光纤放大器的位置以及正负色散光纤的色散差值，可同时获得相对较大的色散管理强度和好的孤子稳定性；总体来讲，集总放大器的位置会对孤子传输的稳定性、色散管理孤子的能量增强因子、脉冲最佳

Evaluation of holonomic quantum computation: Adiabatic versus nonadiabatic

Based on the analytical solution to the time-dependent Schrodinger equations, we evaluate the holonomic quantum computation beyond the adiabatic limit. Besides providing rigorous confirmation of the geometrical prediction of holonomies, the present dynamical resolution offers also a practical means to study the nonadiabaticity induced effects for the universal qubit operations.

Topological structure of the solitons solution in SU(3) Dunne-Jackiw-Pi-Trugenberger model

By using phi-mapping topological current theory and gauge potential decomposition, we discuss the self-dual equation and its solution in the SU(N) Dunne-Jackiw-Pi-Trugenberger model and obtain a new concrete self-dual equation with a 6 function. For the SU(3) case, we obtain a new self-duality solution and find the relationship between the soliton solution and topological number which is determined by the Hopf index and Brouwer degree of phi-mapping. In our solution, the flux of this soliton is naturally quantized.

Transient and nonlinear analysis of slow-light pulse propagation in an optical fiber via stimulated Brillouin scattering

We investigate slow-light pulse propagation in an optical fiber via transient stimulated Brillouin scattering. Space-time evolution of a generating slow-light pulse is numerically calculated by solving three-wave coupled-mode equations between a pump beam, an acoustic wave, and a counterpropagating signal pulse. Our mathematical treatments are applicable to both narrowband and broadband pump cases. We show that the time delay of 85% pulse width can be obtained for a signal pulse of the order of subnanosecond pulse width by using a broadband pump, while the signal pulse is broadened only by 40% of the input signal pulse. The physical origin of the pulse broadening and distortion is explained in terms of the temporal decay of the induced acoustic field. (C) 2009 Optical Society of America

Nonlinear Analysis of Oscillatory Indentation in Elastic and Viscoelastic Solids

Determining the mechanical properties at micro- and nanometer length scales using nanoindentation or atomic force microscopy is important to many areas of science and engineering. Here we establish equations for obtaining storage and loss modulus from oscillatory indentations by performing a nonlinear analysis of conical and spherical indentation in elastic and viscoelastic solids. We show that, when the conical indenter is driven by a sinusoidal force, the square of displacement is a sinusoidal function of time, not the displacement itself, which is commonly assumed. Similar conclusions hold for spherical indentations. Well-known difficulties associated with measuring contact area and correcting thermal drift may be circumvented using the newly derived equations. These results may help improve methods of using oscillatory indentation for determining elastic and viscoelastic properties of solids.

Nonlinear Faraday Waves in a Parametrically Excited Circular Cylindrical Container

In the cylindrical coordinate system, a singular perturbation theory of multiple-scale asymptotic expansions was developed to study single standing water wave mode by solving potential equations of water waves in a rigid circular cylinder, which is subjec

Hill Instability Analysis Of Tlp Tether Subjected To Combined Platform Surge And Heave Motions

This paper presents the Hill instability analysis of Tension Leg Platform (TLP) tether it, deep sea. The 2-D nonlinear beam model which is Undergoing Coupled axial and transverse vibrations, is applied. The governing equations are reduced to nonlinear Hill equation by use of the Galerkin's method and the modes superposition principle. The Hill instability charted Lip to large parameters is obtained. An important parameter M is defined and can he expressed as the functions of tether length, the platform surge and heave motion amplitudes. Some example studies are performed for various environmental conditions. The results demonstrate that the nonlinear coupling between the axial and transverse vibrations has a significant effect on the response of structure.. It needs to be considered for the accurate dynamic analysis of long TLP tether subjected to the combined platform surge and heave motions.

Nonlinear Constitutive Behavior Of Ferroelectric Materials

The ferroelectric specimen is considered as an aggregation of many randomly oriented domains. According to this mechanism, a multi-domain mechanical model is developed in this paper. Each domain is represented by one element. The applied stress and electric field are taken to be the stress and electric field in the formula of the driving force of domain switching for each element in the specimen. It means that the macroscopic switching criterion is used for calculating the volume fraction of domain switching for each element. By using the hardening relation between the driving force of domain switching and the volume fraction of domain switching calibrated, the volume fraction of domain switching for each element is calculated. Substituting the stress and electric field and the volume fraction of domain switching into the constitutive equation of ferroelectric material, one can easily get the strain and electric displacement for each element. The macroscopic behavior of the ferroelectric specimen is then directly calculated by volume averaging. Meanwhile, the nonlinear finite element analysis for the ferroelectric specimen is carried out. In the finite element simulation, the volume fraction of domain switching for each element is calculated by using the same method mentioned above. The interaction between different elements is taken into account in the finite element simulation and the local stress and electric field for each element is obtained. The macroscopic behavior of the specimen is then calculated by volume averaging. The computation results involve the electric butterfly shaped curves of axial strain versus the axial electric field and the hysteresis loops of electric displacement versus the electric field for ferroelectric specimens under the uniaxial coupled stress and electric field loading. The present theoretical prediction agrees reasonably with the experimental results.

Energy Principle And Nonlinear Electric-Mechanical Behavior Of Ferroelectric Ceramics

Many experimental observations have shown that a single domain in a ferroelectric material switches by progressive movement of domain walls, driven by a combination of electric field and stress. The mechanism of the domain switch involves the following steps: initially, the domain has a uniform spontaneous polarization; new domains with the reverse polarization direction nucleate, mainly at the surface, and grow though the crystal thickness; the new domain expands sideways as a new domain continues to form; finally, the domain switch coalesces to complete the polarization reversal. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of the ferroelectric material and used to study the nonlinear constitutive behavior of a ferroelectric body in this paper. The principle of stationary total potential energy is put forward in which the basic unknown quantities are the displacement u(i), electric displacement D-i and volume fraction rho(I) of the domain switching for the variant I. The mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total potential energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion established by Hwang et al. [ 1]. Based on the domain switching criterion, a set of linear algebraic equations for determining the volume fraction rho(I) of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. If the volume fraction rho(I) of domain switching for each domain is prescribed, the unknown displacement and electric potential can be obtained based on the conventional finite element procedure. It is assumed that a domain switches if the reduction in potential energy exceeds a critical energy barrier. According to the experimental results, the energy barrier will strengthen when the volume fraction of the domain switching increases. The external mechanical and electric loads are increased step by step. The volume fraction rho(I) of domain switching for each element obtained from the last loading step is used as input to the constitutive equations. Then the strain and electric fields are calculated based on the conventional finite element procedure. The finite element analysis is carried out on the specimens subjected to uniaxial coupling stress and electric field. Numerical results and available experimental data are compared and discussed. The present theoretic prediction agrees reasonably with the experimental results.

Analysis and application of ellipticity of stability equations on fluid mechanics

By using characteristic analysis of the linear and nonlinear parabolic stability equations (PSE), PSE of primitive disturbance variables are proved to be parabolic intotal. By using sub-characteristic analysis of PSE, the linear PSE are proved to be elliptical and hyperbolic-parabolic for velocity U, in subsonic and supersonic, respectively; the nonlinear PSE are proved to be elliptical and hyperbolic-parabolic for relocity U + u in subsonic and supersonic, respectively. The methods are gained that the remained ellipticity is removed from the PSE by characteristic and sub-characteristic theories, the results for the linear PSE are consistent with the known results, and the influence of the Mach number is also given out. At the same time, the methods of removing the remained ellipticity are further obtained from the nonlinear PSE.