Finite element analysis of flow patterns near geological lenses in hydrodynamic and hydrothermal systems

We use theoretical and numerical methods to investigate the general pore-fluid flow patterns near geological lenses in hydrodynamic and hydrothermal systems respectively. Analytical solutions have been rigorously derived for the pore-fluid velocity, stream function and excess pore-fluid pressure near a circular lens in a hydrodynamic system. These analytical solutions provide not only a better understanding of the physics behind the problem, but also a valuable benchmark solution for validating any numerical method. Since a geological lens is surrounded by a medium of large extent in nature and the finite element method is efficient at modelling only media of finite size, the determination of the size of the computational domain of a finite element model, which is often overlooked by numerical analysts, is very important in order to ensure both the efficiency of the method and the accuracy of the numerical solution obtained. To highlight this issue, we use the derived analytical solutions to deduce a rigorous mathematical formula for designing the computational domain size of a finite element model. The proposed mathematical formula has indicated that, no matter how fine the mesh or how high the order of elements, the desired accuracy of a finite element solution for pore-fluid flow near a geological lens cannot be achieved unless the size of the finite element model is determined appropriately. Once the finite element computational model has been appropriately designed and validated in a hydrodynamic system, it is used to examine general pore-fluid flow patterns near geological lenses in hydrothermal systems. Some interesting conclusions on the behaviour of geological lenses in hydrodynamic and hydrothermal systems have been reached through the analytical and numerical analyses carried out in this paper.

Multidimensional quantum solitons with nondegenerate parametric interactions: Photonic and Bose-Einstein condensate environments

We consider the quantum theory of three fields interacting via parametric and repulsive quartic couplings. This can be applied to treat photonic chi((2)) and chi((3)) interactions, and interactions in atomic Bose-Einstein condensates or quantum Fermi gases, describing coherent molecule formation together with a-wave scattering. The simplest two-particle quantum solitons or bound-state solutions of the idealized Hamiltonian, without a momentum cutoff, are obtained exactly. They have a pointlike structure in two and three dimensions-even though the corresponding classical theory is nonsingular. We show that the solutions can be regularized with a momentum cutoff. The parametric quantum solitons have much more realistic length scales and binding energies than chi((3)) quantum solitons, and the resulting effects could potentially be experimentally tested in highly nonlinear optical parametric media or interacting matter-wave systems. N-particle quantum solitons and the ground state energy are analyzed using a variational approach. Applications to atomic/molecular Bose-Einstein condensates (BEC's) are given, where we predict the possibility of forming coupled BEC solitons in three space dimensions, and analyze superchemistry dynamics.

Quantum chaos in atom optics

I shall discuss the quantum and classical dynamics of a class of nonlinear Hamiltonian systems. The discussion will be restricted to systems with one degree of freedom. Such systems cannot exhibit chaos, unless the Hamiltonians are time dependent. Thus we shall consider systems with a potential function that has a higher than quadratic dependence on the position and, furthermore, we shall allow the potential function to be a periodic function of time. This is the simplest class of Hamiltonian system that can exhibit chaotic dynamics. I shall show how such systems can be realized in atom optics, where very cord atoms interact with optical dipole potentials of a far-off resonance laser. Such systems are ideal for quantum chaos studies as (i) the energy of the atom is small and action scales are of the order of Planck's constant, (ii) the systems are almost perfectly isolated from the decohering effects of the environment and (iii) optical methods enable exquisite time dependent control of the mechanical potentials seen by the atoms.

Minimum-order stable recursive filter design via the genetic algorithm approach

In this paper, the minimum-order stable recursive filter design problem is proposed and investigated. This problem is playing an important role in pipeline implementation sin signal processing. Here, the existence of a high-order stable recursive filter is proved theoretically, in which the upper bound for the highest order of stable filters is given. Then the minimum-order stable linear predictor is obtained via solving an optimization problem. In this paper, the popular genetic algorithm approach is adopted since it is a heuristic probabilistic optimization technique and has been widely used in engineering designs. Finally, an illustrative example is sued to show the effectiveness of the proposed algorithm.

Analysis of steady-state heat transfer through mid-crustal vertical cracks with upward throughflow in hydrothermal systems

We conduct a theoretical analysis of steady-state heat transfer problems through mid-crustal vertical cracks with upward throughflow in hydrothermal systems. In particular, we derive analytical solutions for both the far field and near field of the system. In order to investigate the contribution of the forced advection to the total temperature of the system, two concepts, namely the critical Peclet number and the critical permeability of the system, have been presented and discussed in this paper. The analytical solution for the far field of the system indicates that if the pore-fluid pressure gradient in the crust is lithostatic, the critical permeability of the system can be used to determine whether or not the contribution of the forced advection to the total temperature of the system is negligible. Otherwise, the critical Peclet number should be used. For a crust of moderate thickness, the critical permeability is of the order of magnitude of 10(-20) m(2), under which heat conduction is the overwhelming mechanism to transfer heat energy, even though the pore-fluid pressure gradient in the crust is lithostatic. Furthermore, the lower bound analytical solution for the near field of the system demonstrates that the permeable vertical cracks in the middle crust can efficiently transfer heat energy from the lower crust to the upper crust of the Earth. Copyright (C) 2002 John Wiley Sons, Ltd.

A 2D numerical model for simulating the physics of fault systems

Simulations provide a powerful means to help gain the understanding of crustal fault system physics required to progress towards the goal of earthquake forecasting. Cellular Automata are efficient enough to probe system dynamics but their simplifications render interpretations questionable. In contrast, sophisticated elasto-dynamic models yield more convincing results but are too computationally demanding to explore phase space. To help bridge this gap, we develop a simple 2D elastodynamic model of parallel fault systems. The model is discretised onto a triangular lattice and faults are specified as split nodes along horizontal rows in the lattice. A simple numerical approach is presented for calculating the forces at medium and split nodes such that general nonlinear frictional constitutive relations can be modeled along faults. Single and multi-fault simulation examples are presented using a nonlinear frictional relation that is slip and slip-rate dependent in order to illustrate the model.

Er(3+)-doped silica-hafnia films for optical waveguides and spherical resonators

In this paper we present some result on sol-gel derived silica-hafnia systems. In particular we focus on fabrication, morphological and spectroscopic assessment of Er(3+)-activated thin films. Two examples of silica-hafnia-derived waveguiding glass ceramics, prepared by top-down and bottom-up techniques are reported, and the main optical properties are discussed. Finally, some properties of activated microspherical resonators, having a silica core, obtained by melting the end of a telecom fiber, coated with an Er(3+)-doped 70SiO(2)-30HfO(2) film, are presented. (C) 2009 Elsevier B.V. All rights reserved.

Existence results for abstract impulsive second-order neutral functional differential equations

We establish the existence of mild solutions for a class of impulsive second-order partial neutral functional differential equations with infinite delay in a Banach space. (C) 2009 Published by Elsevier Ltd

Determination of mechanical properties of PECVD silicon nitride thin films for tunable MEMS Fabry-Pérot optical filters

This paper reports an investigation on techniques for determining elastic modulus and intrinsic stress gradient in plasma-enhanced chemical vapor deposition (PECVD) silicon nitride thin films. The elastic property of the silicon nitride thin films was determined using the nanoindentation method on silicon nitride/silicon bilayer systems. A simple empirical formula was developed to deconvolute the film elastic modulus. The intrinsic stress gradient in the films was determined by using micrometric cantilever beams, cross-membrane structures and mechanical simulation. The deflections of the silicon nitride thin film cantilever beams and cross-membranes caused by in-thickness stress gradients were measured using optical interference microscopy. Finite-element beam models were built to compute the deflection induced by the stress gradient. Matching the deflection computed under a given gradient with that measured experimentally on fabricated samples allows the stress gradient of the PECVD silicon nitride thin films introduced from the fabrication process to be evaluated.

Optical Quality in Eyes Implanted With Aspheric and Spherical Intraocular Lenses Assessed by NIDEK OPD-Scan: A Randomized, Bilateral, Clinical Trial

PURPOSE: To determine whether implantation of an intraocular lens (IOL) with an aspheric surface (Akreos AO, Bausch & Lomb Inc) results in reduced ocular aberrations (spherical aberration) and improved Strehl ratio and modulation transfer function (MTF) after cataract surgery. METHODS: In an intraindividual, randomized, double-masked, prospective study of 50 eyes (25 patients) with bilateral cataract, an IOL with modified anterior and posterior surfaces (Akreos AO) was implanted in one eye and a biconvex IOL with spherical surfaces (Akreos Fit, Bausch & Lomb Inc) implanted in the fellow eye. Ocular aberrations, Strehl ratio, and MTF curve with 4.5-, 5.0-, and 6.0-mm pupils were measured with a NIDEK OPD-Scan dynamic retinoscopy aberrometer 3 months after surgery. Uncorrected and corrected distance visual acuity (UDVA and CDVA, respectively) were also measured. RESULTS: No statistically significant difference was noted between eyes in postoperative UDVA and CDVA at 1 month. At 3 months, the Akreos AO IOL group obtained statistically significant lower values of higher order and spherical aberrations with 4.5-, 5.0-, and 6.0-mm pupil diameters than the Akreos Fit IOL group (P<.05). The value of Strehl ratio was statistically significantly higher in eyes with the Akreos AO IOL for 4.5- and 6.0-mm pupils (P<.05). The MTF curve was better in the Akreos AO IOL group in 4.5-, 5.0-, and 6.0-mm pupils (P<.05). CONCLUSIONS: The aspheric Akreos AO IOL induced significantly less spherical aberration than the Akreos Fit IOL for 4.5-, 5.0-, and 6.0-mm pupils. Modulation transfer function and Strehl ratio were also better in eyes implanted with the Akreos AO IOL than the Akreos Fit. [J Refract Surg. 2011;27(4):287-292.] doi:10.3928/1081597X-20100714-01

Direct measurement of pulse distortion near the zero-dispersion wavelength in an optical fiber by frequency-resolved optical gating

The technique of frequency-resolved optical gating is used to characterize the intensity and the phase of picosecond pulses after propagation through 700 m of fiber at close to the zero-dispersion wavelength. Using the frequency-resolved optical gating technique, we directly measure the severe temporal distortion resulting from the interplay between self-phase modulation and higher-order dispersion in this regime. The measured intensity and phase of the pulses after propagation are found to be in good agreement with the predictions of numerical simulations with the nonlinear Schrodinger equation. (C) 1997 Optical Society of America.

Robust algorithms for solving stochastic partial differential equations

A robust semi-implicit central partial difference algorithm for the numerical solution of coupled stochastic parabolic partial differential equations (PDEs) is described. This can be used for calculating correlation functions of systems of interacting stochastic fields. Such field equations can arise in the description of Hamiltonian and open systems in the physics of nonlinear processes, and may include multiplicative noise sources. The algorithm can be used for studying the properties of nonlinear quantum or classical field theories. The general approach is outlined and applied to a specific example, namely the quantum statistical fluctuations of ultra-short optical pulses in chi((2)) parametric waveguides. This example uses a non-diagonal coherent state representation, and correctly predicts the sub-shot noise level spectral fluctuations observed in homodyne detection measurements. It is expected that the methods used wilt be applicable for higher-order correlation functions and other physical problems as well. A stochastic differencing technique for reducing sampling errors is also introduced. This involves solving nonlinear stochastic parabolic PDEs in combination with a reference process, which uses the Wigner representation in the example presented here. A computer implementation on MIMD parallel architectures is discussed. (C) 1997 Academic Press.

Spatiotemporal solitons in multidimensional optical media with a quadratic nonlinearity

We consider solutions to the second-harmonic generation equations in two-and three-dimensional dispersive media in the form of solitons localized in space and time. As is known, collapse does not take place in these models, which is why the solitons may be stable. The general solution is obtained in an approximate analytical form by means of a variational approach, which also allows the stability of the solutions to be predicted. Then, we directly simulate the two-dimensional case, taking the initial configuration as suggested by the variational approximation. We thus demonstrate that spatiotemporal solitons indeed exist and are stable. Furthermore, they are not, in the general case, equivalent to the previously known cylindrical spatial solitons. Direct simulations generate solitons with some internal oscillations. However, these oscillations neither grow nor do they exhibit any significant radiative damping. Numerical solutions of the stationary version of the equations produce the same solitons in their unperturbed form, i.e., without internal oscillations. Strictly stable solitons exist only if the system has anomalous dispersion at both the fundamental harmonic and second harmonic (SH), including the case of zero dispersion at SH. Quasistationary solitons, decaying extremely slowly into radiation, are found in the presence of weak normal dispersion at the second-harmonic frequency.

Evaluation of retinal nerve fiber layer thickness measurements using optical coherence tomography in patients with tobacco-alcohol-induced toxic optic neuropathy

Three patients with progressive visual loss, chronic alcoholism and tabagism were submitted to a complete neuro-ophthalmic examination and to retinal nerve fiber layer (RNFL) measurements using optical coherence tomography (OCT) scanning. Two patients showed marked RNFL loss in the temporal sector of the optic disc. However, a third patient presented RNFL measurements within or above normal limits, based on the Stratus-OCT normative database. Such findings may be due to possible RNFL edema similar to the one that may occur in the acute phase of toxic optic neuropathies. Stratus-OCT was able to detect RNFL loss in the papillomacular bundle of patients with tobacco-alcohol-induced toxic optic neuropathy. However, interpretation must be careful when OCT does not show abnormality in order to prevent diagnostic confusion, since overestimation of RNFL thickness measurements is possible in such cases.

The mu-way intersection problem for m-cycle systems

An m-cycle system of order upsilon is a partition of the edge-set of a complete graph of order upsilon into m-cycles. The mu -way intersection problem for m-cycle systems involves taking mu systems, based on the same vertex set, and determining the possible number of cycles which can be common to all mu systems. General results for arbitrary m are obtained, and detailed intersection values for (mu, m) = (3, 4), (4, 5),(4, 6), (4, 7), (8, 8), (8, 9). (For the case (mu, m)= (2, m), see Billington (J. Combin. Des. 1 (1993) 435); for the case (Cc,m)=(3,3), see Milici and Quattrochi (Ars Combin. A 24 (1987) 175. (C) 2001 Elsevier Science B.V. All rights reserved.