Development and fabrication of an optimized thermo-electro-optic device using a Mach-Zehnder interferometer

The development and fabrication of a thermo-electro-optic sensor using a Mach-Zehnder interferometer and a resistive micro-heater placed in one of the device`s arms is presented. The Mach-Zehnder structure was fabricated on a single crystal silicon substrate using silicon oxynitride and amorphous hydrogenated silicon carbide films to form an anti-resonant reflective optical waveguide. The materials were deposited by Plasma enhanced chemical vapor deposition technique at low temperatures (similar to 320 degrees C). To optimize the heat transfer and increase the device response with current variation, part of the Mach-Zehnder sensor arm was suspended through front-side bulk micromachining of the silicon substrate in a KOH solution. With the temperature variation caused by the micro-heater, the refractive index of the core layer of the optical waveguide changes due to the thermo-optic effect. Since this variation occurs only in one of the Mach-Zehnder`s arm, a phase difference between the arms is produced, leading to electromagnetic interference. In this way, the current applied to the micro-resistor can control the device output optical power. Further, reactive ion etching technique was used in this work to define the device`s geometry, and a study of SF6 based etching rates on different composition of silicon oxynitride films is also presented. (C) 2007 Elsevier B.V. All rights reserved.

Tunable Bragg filter using silicon compound films

In this work, we present the simulation, fabrication and characterization of a tunable Bragg filter employing amorphous dielectric films deposited by plasma enhanced chemical vapor deposition technique on a crystalline silicon substrate. The optical device was built using conventional microelectronic processes and consisted of fifteen periodic intervals of Si3N4 layers separated by air with appropriated thickness and lengths to produce transmittance attenuation peaks in the visible region. For this, previous simulations were realized based in the optical parameters of the dielectric film, which were extracted from ellipsometry and profilometry techniques. For the characterization of the optical interferential filter, a 633 nm monochromatic light was injected on the filter, and then the transmitted output light was collected and conducted to a detector through an optical waveguide made also of amorphous dielectric layers. Afterwards, the optical filter was mounted on a Peltier thermoelectric device in order to control the temperature of the optical device. When the temperature of filter changes, a refractive index variation is originated in the dielectric film due to the thermo-optic effect, producing a shift of attenuation peak, which can be well predicted by numerical simulations. This characteristic allows this device to be used as a thermo-optic sensor. (C) 2007 Elsevier B.V. All rights reserved.

PbO-GeO(2) rib waveguides for photonic applications

This work presents for the first time to our knowledge the fabrication and characterization of rib waveguides produced with PbO-GeO(2) (PGO) thin films. The target was manufactured using pure oxides ( 60 PbO-40 GeO(2), in wt%) and amorphous thin films were produced with the RF sputtering technique. PGO thin films present small absorption in the visible and in the near infrared and refractive index of similar to 2.0. The definition of the rib waveguide structure was made using conventional optical lithography followed by plasma etching, performed in a Reactive Ion Etching (RIE) reactor. Light propagation mode in the waveguide structure was analyzed using integrated optic simulation software. Optical loss measurements were performed to determine the propagation loss at 633 nm, for ribs with height of 70 nm and width of 3-5 mu m; experimental values around 2 dB/cm were found for the propagation loss and confirmed the theoretical calculations. The results obtained demonstrate that PGO thin films are potential candidates for application in integrated optics. Published by Elsevier B.V.

Constraints on the spatiotemporal accuracy of interceptive action: effects of target size on hitting a moving object

Results of two experiments are reported that examined how people respond to rectangular targets of different sizes in simple hitting tasks. If a target moves in a straight line and a person is constrained to move along a linear track oriented perpendicular to the targetrsquos motion, then the length of the target along its direction of motion constrains the temporal accuracy and precision required to make the interception. The dimensions of the target perpendicular to its direction of motion place no constraints on performance in such a task. In contrast, if the person is not constrained to move along a straight track, the targetrsquos dimensions may constrain the spatial as well as the temporal accuracy and precision. The experiments reported here examined how people responded to targets of different vertical extent (height): the task was to strike targets that moved along a straight, horizontal path. In experiment 1 participants were constrained to move along a horizontal linear track to strike targets and so target height did not constrain performance. Target height, length and speed were co-varied. Movement time (MT) was unaffected by target height but was systematically affected by length (briefer movements to smaller targets) and speed (briefer movements to faster targets). Peak movement speed (Vmax) was influenced by all three independent variables: participants struck shorter, narrower and faster targets harder. In experiment 2, participants were constrained to move in a vertical plane normal to the targetrsquos direction of motion. In this task target height constrains the spatial accuracy required to contact the target. Three groups of eight participants struck targets of different height but of constant length and speed, hence constant temporal accuracy demand (different for each group, one group struck stationary targets = no temporal accuracy demand). On average, participants showed little or no systematic response to changes in spatial accuracy demand on any dependent measure (MT, Vmax, spatial variable error). The results are interpreted in relation to previous results on movements aimed at stationary targets in the absence of visual feedback.

Experimental Analysis of Froude Number Effect on Air Entrainment in the Hydraulic Jump

A hydraulic jump is characterized by strong energy dissipation and mixing, large-scale turbulence, air entrainment, waves and spray. Despite recent pertinent studies, the interaction between air bubbles diffusion and momentum transfer is not completely understood. The objective of this paper is to present experimental results from new measurements performed in rectangular horizontal flume with partially-developed inflow conditions. The vertical distributions of void fraction and air bubbles count rate were recorded for inflow Froude number Fr1 in the range from 5.2 to 14.3. Rapid detrainment process was observed near the jump toe, whereas the structure of the air diffusion layer was clearly observed over longer distances. These new data were compared with previous data generally collected at lower Froude numbers. The comparison demonstrated that, at a fixed distance from the jump toe, the maximum void fraction Cmax increases with the increasing Fr1. The vertical locations of the maximum void fraction and bubble count rate were consistent with previous studies. Finally, an empirical correlation between the upper boundary of the air diffusion layer and the distance from the impingement point was provided.

Clustering Process Analysis in a Large-Size Dropshaft and in a Hydraulic Jump

Recent efforts in the characterization of air-water flows properties have included some clustering process analysis. A cluster of bubbles is defined as a group of two or more bubbles, with a distinct separation from other bubbles before and after the cluster. The present paper compares the results of clustering processes two hydraulic structures. That is, a large-size dropshaft and a hydraulic jump in a rectangular horizontal channel. The comparison highlighted some significant differences in clustering production and structures. Both dropshaft and hydraulic jump flows are complex turbulent shear flows, and some clustering index may provide some measure of the bubble-turbulence interactions and associated energy dissipation.

Finite element analysis of steady-state natural convection problems in fluid-saturated porous media heated from below

In this paper, a progressive asymptotic approach procedure is presented for solving the steady-state Horton-Rogers-Lapwood problem in a fluid-saturated porous medium. The Horton-Rogers-Lapwood problem possesses a bifurcation and, therefore, makes the direct use of conventional finite element methods difficult. Even if the Rayleigh number is high enough to drive the occurrence of natural convection in a fluid-saturated porous medium, the conventional methods will often produce a trivial non-convective solution. This difficulty can be overcome using the progressive asymptotic approach procedure associated with the finite element method. The method considers a series of modified Horton-Rogers-Lapwood problems in which gravity is assumed to tilt a small angle away from vertical. The main idea behind the progressive asymptotic approach procedure is that through solving a sequence of such modified problems with decreasing tilt, an accurate non-zero velocity solution to the Horton-Rogers-Lapwood problem can be obtained. This solution provides a very good initial prediction for the solution to the original Horton-Rogers-Lapwood problem so that the non-zero velocity solution can be successfully obtained when the tilted angle is set to zero. Comparison of numerical solutions with analytical ones to a benchmark problem of any rectangular geometry has demonstrated the usefulness of the present progressive asymptotic approach procedure. Finally, the procedure has been used to investigate the effect of basin shapes on natural convection of pore-fluid in a porous medium. (C) 1997 by John Wiley & Sons, Ltd.

Characteristics of undular hydraulic jumps: Experiments and analysis

The writers measured velocity, pressure and energy distributions, wavelengths, and wave amplitudes along undular jumps in a smooth rectangular channel 0.25 m wide. In each case the upstream flow was a fully developed shear flow. Analysis of the data shows that the jump has strong three-dimensional features and that the aspect ratio of the channel is an important parameter. Energy dissipation on the centerline is far from negligible and is largely constrained to the reach between the start of the lateral shock waves and the first wave crest of the jump, in which the boundary layer develops under a strong adverse pressure gradient. A Boussinesq-type solution of the free-surface profile, velocity, and energy and pressure distributions is developed and compared with the data. Limitations of the two-dimensional analysis are discussed.

Theory of multidimensional parametric band-gap simultons

Multidimensional spatiotemporal parametric simultons (simultaneous solitary waves) are possible in a nonlinear chi((2)) medium with a Bragg grating structure, where large effective dispersion occurs near two resonant band gaps for the carrier and second-harmonic field, respectively. The enhanced dispersion allows much reduced interaction lengths, as compared to bulk medium parametric simultons. The nonlinear parametric band-gap medium permits higher-dimensional stationary waves to form. In addition, solitons can occur with lower input powers than conventional nonlinear Schrodinger equation gap solitons. In this paper, the equations for electromagnetic propagation in a grating structure with a parametric nonlinearity are derived from Maxwell's equation using a coupled mode Hamiltonian analysis in one, two, and three spatial dimensions. Simultaneous solitary wave solutions are proved to exist by reducing the equations to the coupled equations describing a nonlinear parametric waveguide, using the effective-mass approximation (EMA). Exact one-dimensional numerical solutions in agreement with the EMA solutions are also given. Direct numerical simulations show that the solutions have similar types of stability properties to the bulk case, providing the carrier waves are tuned to the two Bragg resonances, and the pulses have a width in frequency space less than the band gap. In summary, these equations describe a physically accessible localized nonlinear wave that is stable in up to 3 + 1 dimensions. Possible applications include photonic logic and switching devices. [S1063-651X(98)06109-1].

Finite element modelling of temperature gradient driven rock alteration and mineralization in porous rock masses

We present finite element simulations of temperature gradient driven rock alteration and mineralization in fluid saturated porous rock masses. In particular, we explore the significance of production/annihilation terms in the mass balance equations and the dependence of the spatial patterns of rock alteration upon the ratio of the roll over time of large scale convection cells to the relaxation time of the chemical reactions. Special concepts such as the gradient reaction criterion or rock alteration index (RAI) are discussed in light of the present, more general theory. In order to validate the finite element simulation, we derive an analytical solution for the rock alteration index of a benchmark problem on a two-dimensional rectangular domain. Since the geometry and boundary conditions of the benchmark problem can be easily and exactly modelled, the analytical solution is also useful for validating other numerical methods, such as the finite difference method and the boundary element method, when they are used to dear with this kind of problem. Finally, the potential of the theory is illustrated by means of finite element studies related to coupled flow problems in materially homogeneous and inhomogeneous porous rock masses. (C) 1998 Elsevier Science S.A. All rights reserved.

Modelling the large deformations in stratified media - the Cosserat continuum approach

Methods employing continuum approximation in describing the deformation of layered materials possess a clear advantage over explicit models, However, the conventional implicit models based on the theory of anisotropic continua suffers from certain difficulties associated with interface slip and internal instabilities. These difficulties can be remedied by considering the bending stiffness of the layers. This implies the introduction of moment (couple) stresses and internal rotations, which leads to a Cosserat-type theory. In the present model, the behaviour of the layered material is assumed to be linearly elastic; the interfaces are assumed to be elastic perfectly plastic. Conditions of slip or no slip at the interfaces are detected by a Coulomb criterion with tension cut off at zero normal stress. The theory is valid for large deformation analysis. The model is incorporated into the finite element program AFENA and validated against analytical solutions of elementary buckling problems in layered medium. A problem associated with buckling of the roof and the floor of a rectangular excavation in jointed rock mass under high horizontal in situ stresses is considered as the main application of the theory. Copyright (C) 1999 John Wiley & Sons, Ltd.

Computation of bifurcation boundaries for power systems: A new Delta-plane method

This paper is devoted to the problems of finding the load flow feasibility, saddle node, and Hopf bifurcation boundaries in the space of power system parameters. The first part contains a review of the existing relevant approaches including not-so-well-known contributions from Russia. The second part presents a new robust method for finding the power system load flow feasibility boundary on the plane defined by any three vectors of dependent variables (nodal voltages), called the Delta plane. The method exploits some quadratic and linear properties of the load now equations and state matrices written in rectangular coordinates. An advantage of the method is that it does not require an iterative solution of nonlinear equations (except the eigenvalue problem). In addition to benefits for visualization, the method is a useful tool for topological studies of power system multiple solution structures and stability domains. Although the power system application is developed, the method can be equally efficient for any quadratic algebraic problem.

NMR shim coil design utilising a rapid spherical harmonic calculation method

A rapid spherical harmonic calculation method is used for the design of Nuclear Magnetic Resonance shim coils. The aim is to design each shim such that it generates a field described purely by a single spherical harmonic. By applying simulated annealing techniques, coil arrangements are produced through the optimal positioning of current-carrying circular arc conductors of rectangular cross-section. This involves minimizing the undesirable harmonies in relation to a target harmonic. The design method is flexible enough to be applied for the production of coil arrangements that generate fields consisting significantly of either zonal or tesseral harmonics. Results are presented for several coil designs which generate tesseral harmonics of degree one.

A director theory for visco-elastic folding instabilities in multilayered rock

A model for finely layered visco-elastic rock proposed by us in previous papers is revisited and generalized to include couple stresses. We begin with an outline of the governing equations for the standard continuum case and apply a computational simulation scheme suitable for problems involving very large deformations. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered beam under compression. We analyse folding up to 40% shortening. The standard continuum solution becomes unstable for extreme values of the shear/normal viscosity ratio. The instability is a consequence of the neglect of the bending stiffness/viscosity in the standard continuum model. We suggest considering these effects within the framework of a couple stress theory. Couple stress theories involve second order spatial derivatives of the velocities/displacements in the virtual work principle. To avoid C-1 continuity in the finite element formulation we introduce the spin of the cross sections of the individual layers as an independent variable and enforce equality to the spin of the unit normal vector to the layers (-the director of the layer system-) by means of a penalty method. We illustrate the convergence of the penalty method by means of numerical solutions of simple shears of an infinite layer for increasing values of the penalty parameter. For the shear problem we present solutions assuming that the internal layering is oriented orthogonal to the surfaces of the shear layer initially. For high values of the ratio of the normal-to the shear viscosity the deformation concentrates in thin bands around to the layer surfaces. The effect of couple stresses on the evolution of folds in layered structures is also investigated. (C) 2002 Elsevier Science Ltd. All rights reserved.

Large amplitude folding in finely layered viscoelastic rock structures

We analyze folding phenomena in finely layered viscoelastic rock. Fine is meant in the sense that the thickness of each layer is considerably smaller than characteristic structural dimensions. For this purpose we derive constitutive relations and apply a computational simulation scheme (a finite-element based particle advection scheme; see MORESI et al., 2001) suitable for problems involving very large deformations of layered viscous and viscoelastic rocks. An algorithm for the time integration of the governing equations as well as details of the finite-element implementation is also given. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered plate. The domain is compressed along the layer axis by prescribing velocities along the sides. First, for the viscous limit we consider the response to a series of harmonic perturbations of the director orientation. The Fourier spectra of the initial folding velocity are compared for different viscosity ratios. Turning to the nonlinear regime we analyze viscoelastic folding histories up to 40% shortening. The effect of layering manifests itself in that appreciable buckling instabilities are obtained at much lower viscosity ratios (1:10) as is required for the buckling of isotropic plates (1:500). The wavelength induced by the initial harmonic perturbation of the director orientation seems to be persistent. In the section of the parameter space considered here elasticity seems to delay or inhibit the occurrence of a second, larger wavelength. Finally, in a linear instability analysis we undertake a brief excursion into the potential role of couple stresses on the folding process. The linear instability analysis also provides insight into the expected modes of deformation at the onset of instability, and the different regimes of behavior one might expect to observe.