Two-Point Closure Strategy in the Mapping Closure Approximation Approach

A two-point closure strategy in mapping closure approximation (MCA) approach is developed for the evolution of the probability density function (PDF) of a scalar advected by stochastic velocity fields. The MCA approach is based on multipoint statistics. We formulate a MCA modeled system using the one-point PDFs and two-point correlations. The MCA models can describe both the evolution of the PDF shape and the rate at which the PDF evolves.

Numerical approximation improves well survey calculation

It is now possible to improve the precision of well survey calculations by order of magnitude with numerical approximation.

Although the most precise method of simulating and calculating a wellbore trajectory generally requires more calculation than other, less-accurate methods, the wider use of computers in oil fields now eliminates this as an obstacle.

The results of various calculations show that there is a deviation of more than 10 m among the different methods of calculation for a directional well of 3,000 m.¹ Consequently, it is important to improve the precision and reliability of survey calculation-the fundamental, necessary work of quantitatively monitoring and controlling wellbore trajectories.

Mapping Closure Approximation for Reactive Scalars in Random Flows

The Mapping Closure Approximation (MCA) approach is developed to describe the statistics of both conserved and reactive scalars in random flows. The statistics include Probability Density Function (PDF), Conditional Dissipation Rate (CDR) and Conditional Laplacian (CL). The statistical quantities are calculated using the MCA and compared with the results of the Direct Numerical Simulation (DNS). The results obtained from the MCA are in agreement with those from the DNS. It is shown that the MCA approach can predict the statistics of reactive scalars in random flows.

Seismic Responses of a Hemispherical Alluvial Valley to SV Waves: A Three-Dimensional Analytical Approximation

An analytical solution to the three-dimensional scattering and diffraction of plane SV-waves by a saturated hemispherical alluvial valley in elastic half-space is obtained by using Fourier-Bessel series expansion technique. The hemispherical alluvial valley with saturated soil deposits is simulated with Biot's dynamic theory for saturated porous media. The following conclusions based on numerical results can be drawn: (1) there are a significant differences in the seismic response simulation between the previous single-phase models and the present two-phase model; (2) the normalized displacements on the free surface of the alluvial valley depend mainly on the incident wave angles, the dimensionless frequency of the incident SV waves and the porosity of sediments; (3) with the increase of the incident angle, the displacement distributions become more complicated; and the displacements on the free surface of the alluvial valley increase as the porosity of sediments increases.

Compact difference approximation with consistent boundary condition

For simulating multi-scale complex flow fields it should be noted that all the physical quantities we are interested in must be simulated well. With limitation of the computer resources it is preferred to use high order accurate difference schemes. Because of their high accuracy and small stencil of grid points computational fluid dynamics (CFD) workers pay more attention to compact schemes recently. For simulating the complex flow fields the treatment of boundary conditions at the far field boundary points and near far field boundary points is very important. According to authors' experience and published results some aspects of boundary condition treatment for far field boundary are presented, and the emphasis is on treatment of boundary conditions for the upwind compact schemes. The consistent treatment of boundary conditions at the near boundary points is also discussed. At the end of the paper are given some numerical examples. The computed results with presented method are satisfactory.

Evolution and propagation of density waves on galactic disk .2. An approximation of the thin transition layer

In this part of the present work, a simplified model—the thin transition layer theory is proposed. The comparison of this model with the G-L sheet model is made.

The accuracy of Kirchhoff's approximation in describing the far field speckles produced by random self-affine fractal surfaces

Based on the rigorous formulation of integral equations for the propagations of light waves at the medium interface, we carry out the numerical solutions of the random light field scattered from self-affine fractal surface samples. The light intensities produced by the same surface samples are also calculated in Kirchhoff's approximation, and their comparisons with the corresponding rigorous results show directly the degree of the accuracy of the approximation. It is indicated that Kirchhoff's approximation is of good accuracy for random surfaces with small roughness value w and large roughness exponent alpha. For random surfaces with larger w and smaller alpha, the approximation results in considerable errors, and detailed calculations show that the inaccuracy comes from the simplification that the transmitted light field is proportional to the incident field and from the neglect of light field derivative at the interface.

Approximate analytical description for fundamental-mode fields of graded-index fibers: Beyond the Gaussian approximation

An approximate analytical description for fundamental-mode fields of graded-index fibers is explicitly presented by use of the power-series expansion method, the maximum-value condition at the fiber axis, the decay properties of fundamental-mode fields at large distance from the fiber axis, and the approximate modal parameters U obtained from the Gaussian approximation. This analytical description is much more accurate than the Gaussian approximation and at the same time keep the simplicity of the latter. As two special examples, we present the approximate analytical formulas for the fundamental-mode fields of a step profile fiber and a Gaussian profile fiber, and we find that they are both highly accurate in the single-mode range by comparing them with the corresponding exact solutions.

Corrections to the paraxial approximation of an arbitrary free-propagation beam

A relatively simple transform from an arbitrary solution of the paraxial wave equation to the corresponding exact solution of the Helmholtz wave equation is derived in the condition that the evanescent waves are ignored and is used to study the corrections to the paraxial approximation of an arbitrary free-propagation beam. Specifically, the general lowest-order correction field is given in a very simple form and is proved to be exactly consistent with the perturbation method developed by Lax et nl. [Phys. Rev. A 11, 1365 (1975)]. Some special examples, such as the lowest-order correction to the paraxial approximation of a fundamental Gaussian beam whose waist plane has a parallel shin from the z = 0 plane, are presented. (C) 1998 Optical Society of America.

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].

Application of fast Pade approximation in simulating photonic crystal nanocavities by FDTD technology

The exact calculation of mode quality factor Q is a key problem in the design of high-Q photonic crystal nanocavity. On the basis of further investigation on conventional Pade approximation, FDM and DFT, Pade approximation with Baker's algorithm is enhanced through introducing multiple frequency search and parabola interpolation. Though Pade approximation is a nonlinear signal processing method and only short time sequence is needed, we find the different length of sequence requirements for 2D and 3D FDTD, which is very important to obtain convergent and accurate results. By using the modified Pade approximation method and 3D FDTD, the 2D slab photonic crystal nanocavity is analyzed and high-Q multimode can be solved quickly instead of large range high-resolution scanning. Monitor position has also been investigated. These results are very helpful to the design of photonic crystal nanocavity devices. (C) 2008 Elsevier B.V. All rights reserved.

Calculation of light delay for coupled microrings by FDTD technique and Pade approximation

The Pade approximation with Baker's algorithm is compared with the least-squares Prony method and the generalized pencil-of-functions (GPOF) method for calculating mode frequencies and mode Q factors for coupled optical microdisks by FDTD technique. Comparisons of intensity spectra and the corresponding mode frequencies and Q factors show that the Pade approximation can yield more stable results than the Prony and the GPOF methods, especially the intensity spectrum. The results of the Prony method and the GPOF method are greatly influenced by the selected number of resonant modes, which need to be optimized during the data processing, in addition to the length of the time response signal. Furthermore, the Pade approximation is applied to calculate light delay for embedded microring resonators from complex transmission spectra obtained by the Pade approximation from a FDTD output. The Prony and the GPOF methods cannot be applied to calculate the transmission spectra, because the transmission signal obtained by the FDTD simulation cannot be expressed as a sum of damped complex exponentials. (C) 2009 Optical Society of America

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.

Method for detecting concurrence without the structural physical approximation by local operations and classical communication

We present a modified method for detecting the concurrence in an arbitrary two-qubit quantum state rho(AB) with local operations and classical communication. In this method, it is not necessary for the two observers to prepare the quantum state rho(AB) by the structural physical approximation. Their main task is to measure four specific functions via two local quantum networks. With these functions they can determine the concurrence and then the entanglement of formation.

Computation of resonant frequencies and quality factors of cavities by FDTD technique and Pade approximation

The finite-difference time domain (FDTD) technique and the Pade approximation with Baker's algorithm are used to calculate the mode frequencies and quality factors of cavities. Comparing with the fast Fourier transformation/Pade method, we find that the Fade approximation and the Baker's algorithm can obtain exact resonant frequencies and quality factors based on a much shorter time record of the FDTD output.