Mode Characteristics for Square Resonators With a Metal Confinement Layer

Microsquare resonators laterally confined by SiO2/Au/air multilayer structure are investigated by light ray method with reflection phase-shift of the multiple layers and two-dimensional (2-D) finite-difference time-domain (FDTD) technique. The reflectivity and phase shift of the mode light ray on the sides of the square resonator with the semiconductor/SiO2/Au/air multilayer structure are calculated for TE and TM modes by transfer matrix method. Based on the reflection phase shift and the reflectivity, the mode wavelength and factor are calculated by the resonant condition and the mirror loss, which are in agreement well with that obtained by the FDTD simulation. We find that the mode factor increases greatly with the increase of the SiO2 layer thickness, especially as d < 0.3 mu m. For the square resonator with side length 2 mu m and refractive index 3.2, anticrossing mode couplings are found for confined TE modes at wavelength about 1.6 mu m at d = 0.11 mu m, and confined TM modes at d = 0.71 mu m, respectively.

Numerical simulation of mesoscopic mechanical behaviors of gradual multi-fiber-reinforced polymer matrix composites

To simulate the deformation and the fracture of gradual multi-fiber-reinforced matrix composites, a numerical simulation method for the mesoscopic mechanical behaviors was developed on the basis of the finite element and the Monte Carlo methods. The results indicate that the strength of a composite increases if the variability of statistical fiber strengths is decreased.

Multi-zone TAP-reactors theory and application - IV. Ideal and non-ideal boundary conditions

We study the influence of non-ideal boundary and initial conditions (BIC) of a temporal analysis of products (TAP) reactor model on the data (observed exit flux) analysis. The general theory of multi-response state-defining experiments for a multi-zone TAP reactor is extended and applied to model several alternative boundary and initial conditions proposed in the literature. The method used is based on the Laplace transform and the transfer matrix formalism for multi-response experiments. Two non-idealities are studied: (1) the inlet pulse not being narrow enough (gas pulse not entering the reactor in Dirac delta function shape) and (2) the outlet non-ideality due to imperfect vacuum. The effect of these non-idealities is analyzed to the first and second order of approximation. The corresponding corrections were obtained and discussed in detail. It was found that they are negligible. Therefore, the model with ideal boundary conditions is proven to be completely adequate to the description and interpretation of transport-reaction data obtained with TAP-2 reactors.

Inspeção de materiais compósitos utilizando ondas acústicas guiadas

Pós-graduação em Engenharia Elétrica - FEIS

Numerical simulation of interdigitated back contact hetero-junction solar cells

The goal of this thesis is the application of an opto-electronic numerical simulation to heterojunction silicon solar cells featuring an all back contact architecture (Interdigitated Back Contact Hetero-Junction IBC-HJ). The studied structure exhibits both metal contacts, emitter and base, at the back surface of the cell with the objective to reduce the optical losses due to the shadowing by front contact of conventional photovoltaic devices. Overall, IBC-HJ are promising low-cost alternatives to monocrystalline wafer-based solar cells featuring front and back contact schemes, in fact, for IBC-HJ the high concentration doping diffusions are replaced by low-temperature deposition processes of thin amorphous silicon layers. Furthermore, another advantage of IBC solar cells with reference to conventional architectures is the possibility to enable a low-cost assembling of photovoltaic modules, being all contacts on the same side. A preliminary extensive literature survey has been helpful to highlight the specific critical aspects of IBC-HJ solar cells as well as the state-of-the-art of their modeling, processing and performance of practical devices. In order to perform the analysis of IBC-HJ devices, a two-dimensional (2-D) numerical simulation flow has been set up. A commercial device simulator based on finite-difference method to solve numerically the whole set of equations governing the electrical transport in semiconductor materials (Sentuarus Device by Synopsys) has been adopted. The first activity carried out during this work has been the definition of a 2-D geometry corresponding to the simulation domain and the specification of the electrical and optical properties of materials. In order to calculate the main figures of merit of the investigated solar cells, the spatially resolved photon absorption rate map has been calculated by means of an optical simulator. Optical simulations have been performed by using two different methods depending upon the geometrical features of the front interface of the solar cell: the transfer matrix method (TMM) and the raytracing (RT). The first method allows to model light prop-agation by plane waves within one-dimensional spatial domains under the assumption of devices exhibiting stacks of parallel layers with planar interfaces. In addition, TMM is suitable for the simulation of thin multi-layer anti reflection coating layers for the reduction of the amount of reflected light at the front interface. Raytracing is required for three-dimensional optical simulations of upright pyramidal textured surfaces which are widely adopted to significantly reduce the reflection at the front surface. The optical generation profiles are interpolated onto the electrical grid adopted by the device simulator which solves the carriers transport equations coupled with Poisson and continuity equations in a self-consistent way. The main figures of merit are calculated by means of a postprocessing of the output data from device simulation. After the validation of the simulation methodology by means of comparison of the simulation result with literature data, the ultimate efficiency of the IBC-HJ architecture has been calculated. By accounting for all optical losses, IBC-HJ solar cells result in a theoretical maximum efficiency above 23.5% (without texturing at front interface) higher than that of both standard homojunction crystalline silicon (Homogeneous Emitter HE) and front contact heterojuction (Heterojunction with Intrinsic Thin layer HIT) solar cells. However it is clear that the criticalities of this structure are mainly due to the defects density and to the poor carriers transport mobility in the amorphous silicon layers. Lastly, the influence of the most critical geometrical and physical parameters on the main figures of merit have been investigated by applying the numerical simulation tool set-up during the first part of the present thesis. Simulations have highlighted that carrier mobility and defects level in amorphous silicon may lead to a potentially significant reduction of the conversion efficiency.

The scintillating fiber focal plane detector for the use of Kaos as a double arm spectrometer

The upgrade of the Mainz Mikrotron (MAMI) electron accelerator facility in 2007 which raised the beam energy up to 1.5,GeV, gives the opportunity to study strangeness production channels through electromagnetic process. The Kaon Spectrometer (KAOS) managed by the A1 Collaboration, enables the efficient detection of the kaons associated with strangeness electroproduction. Used as a single arm spectrometer, it can be combined with the existing high-resolution spectrometers for exclusive measurements in the kinematic domain accessible to them.rnrnFor studying hypernuclear production in the ^A Z(e,e'K^+) _Lambda ^A(Z-1) reaction, the detection of electrons at very forward angles is needed. Therefore, the use of KAOS as a double-arm spectrometer for detection of kaons and the electrons at the same time is mandatory. Thus, the electron arm should be provided with a new detector package, with high counting rate capability and high granularity for a good spatial resolution. To this end, a new state-of-the-art scintillating fiber hodoscope has been developed as an electron detector.rnrnThe hodoscope is made of two planes with a total of 18432 scintillating double-clad fibers of 0.83 mm diameter. Each plane is formed by 72 modules. Each module is formed from a 60deg slanted multi-layer bundle, where 4 fibers of a tilted column are connected to a common read out. The read-out is made with 32 channels of linear array multianode photomultipliers. Signal processing makes use of newly developed double-threshold discriminators. The discriminated signal is sent in parallel to dead-time free time-to-digital modules and to logic modules for triggering purposes.rnrnTwo fiber modules were tested with a carbon beam at GSI, showing a time resolution of 220 ps (FWHM) and a position residual of 270 microm m (FWHM) with a detection efficiency epsilon>99%.rnrnThe characterization of the spectrometer arm has been achieved through simulations calculating the transfer matrix of track parameters from the fiber detector focal plane to the primary vertex. This transfer matrix has been calculated to first order using beam transport optics and has been checked by quasielastic scattering off a carbon target, where the full kinematics is determined by measuring the recoil proton momentum. The reconstruction accuracy for the emission parameters at the quasielastic vertex was found to be on the order of 0.3 % in first test realized.rnrnThe design, construction process, commissioning, testing and characterization of the fiber hodoscope are presented in this work which has been developed at the Institut für Kernphysik of the Johannes Gutenberg - Universität Mainz.

Mathematical modelling of controlled drug release from polymer micro-spheres : incorporating the effects of swelling, diffusion and dissolution via moving boundary problems

Controlled drug delivery is a key topic in modern pharmacotherapy, where controlled drug delivery devices are required to prolong the period of release, maintain a constant release rate, or release the drug with a predetermined release profile. In the pharmaceutical industry, the development process of a controlled drug delivery device may be facilitated enormously by the mathematical modelling of drug release mechanisms, directly decreasing the number of necessary experiments. Such mathematical modelling is difficult because several mechanisms are involved during the drug release process. The main drug release mechanisms of a controlled release device are based on the device’s physiochemical properties, and include diffusion, swelling and erosion. In this thesis, four controlled drug delivery models are investigated. These four models selectively involve the solvent penetration into the polymeric device, the swelling of the polymer, the polymer erosion and the drug diffusion out of the device but all share two common key features. The first is that the solvent penetration into the polymer causes the transition of the polymer from a glassy state into a rubbery state. The interface between the two states of the polymer is modelled as a moving boundary and the speed of this interface is governed by a kinetic law. The second feature is that drug diffusion only happens in the rubbery region of the polymer, with a nonlinear diffusion coefficient which is dependent on the concentration of solvent. These models are analysed by using both formal asymptotics and numerical computation, where front-fixing methods and the method of lines with finite difference approximations are used to solve these models numerically. This numerical scheme is conservative, accurate and easily implemented to the moving boundary problems and is thoroughly explained in Section 3.2. From the small time asymptotic analysis in Sections 5.3.1, 6.3.1 and 7.2.1, these models exhibit the non-Fickian behaviour referred to as Case II diffusion, and an initial constant rate of drug release which is appealing to the pharmaceutical industry because this indicates zeroorder release. The numerical results of the models qualitatively confirms the experimental behaviour identified in the literature. The knowledge obtained from investigating these models can help to develop more complex multi-layered drug delivery devices in order to achieve sophisticated drug release profiles. A multi-layer matrix tablet, which consists of a number of polymer layers designed to provide sustainable and constant drug release or bimodal drug release, is also discussed in this research. The moving boundary problem describing the solvent penetration into the polymer also arises in melting and freezing problems which have been modelled as the classical onephase Stefan problem. The classical one-phase Stefan problem has unrealistic singularities existed in the problem at the complete melting time. Hence we investigate the effect of including the kinetic undercooling to the melting problem and this problem is called the one-phase Stefan problem with kinetic undercooling. Interestingly we discover the unrealistic singularities existed in the classical one-phase Stefan problem at the complete melting time are regularised and also find out the small time behaviour of the one-phase Stefan problem with kinetic undercooling is different to the classical one-phase Stefan problem from the small time asymptotic analysis in Section 3.3. In the case of melting very small particles, it is known that surface tension effects are important. The effect of including the surface tension to the melting problem for nanoparticles (no kinetic undercooling) has been investigated in the past, however the one-phase Stefan problem with surface tension exhibits finite-time blow-up. Therefore we investigate the effect of including both the surface tension and kinetic undercooling to the melting problem for nanoparticles and find out the the solution continues to exist until complete melting. The investigation of including kinetic undercooling and surface tension to the melting problems reveals more insight into the regularisations of unphysical singularities in the classical one-phase Stefan problem. This investigation gives a better understanding of melting a particle, and contributes to the current body of knowledge related to melting and freezing due to heat conduction.

Design of multi-frequency microstrip antennas using multiple rings

An electromagnetically coupled feed arrangement is proposed for simultaneously exciting multiple concentric ring antennas for multi-frequency operation. This has a multi-layer dielectric configuration in which a transmission line is embedded below the layer containing radiating rings. Energy coupled to these rings from the line beneath is optimised by suitably adjusting the location and dimensions of stubs on the line. It has been shown that the resonant frequencies of these rings do not change as several of these single-frequency antennas are combined to form a multi-resonant antenna. Furthermore, all radiators are forced to operate at their primary mode and some harmonics of the lower resonant frequency rings appearing within the frequency range are suppressed when combined. The experimental prototype antenna has three resonant frequencies at which it has good radiation characteristics.

An ANN based approach for voltage stability assessment

Artificial Neural Networks (ANNs) have recently been proposed as an alterative method for salving certain traditional problems in power systems where conventional techniques have not achieved the desired speed, accuracy or efficiency. This paper presents application of ANN where the aim is to achieve fast voltage stability margin assessment of power network in an energy control centre (ECC), with reduced number of appropriate inputs. L-index has been used for assessing voltage stability margin. Investigations are carried out on the influence of information encompassed in input vector and target out put vector, on the learning time and test performance of multi layer perceptron (MLP) based ANN model. LP based algorithm for voltage stability improvement, is used for generating meaningful training patterns in the normal operating range of the system. From the generated set of training patterns, appropriate training patterns are selected based on statistical correlation process, sensitivity matrix approach, contingency ranking approach and concentric relaxation method. Simulation results on a 24 bus EHV system, 30 bus modified IEEE system, and a 82 bus Indian power network are presented for illustration purposes.

Observation and modeling for terrestrial processes in alpine meadow

The water-heat transfer process between land and atmosphere in Haibei alpine meadow area has been systematically observed. A multi-layer coupling model for land-atmosphere interaction was presented with special attention paid to the moisture transfer in leaf stomata under unsaturated condition. A profound investigation on the physical process of turbulent transfer inside the vegetation has been performed with a revised formula of water absorption for root system. The present model facilitates the study of vertically distributed physical variables in detail. Numerical simulation was conducted according to the transfer process of Kinesia humility meadow in the area of Haibei Alpine Meadow Ecosystem Station, CAS. The calculated results agree well with observation.

Oscillatory Instability of Rayleigh-Marangoni-Benard Convection in Two-Layer System

The convective instabilities in two or more superposed layers heated from below were studied extensively by many scientists due to several interfacial phenomena in nature and crystal growth application. Most works of them were performed mainly on the instability behaviors induced only by buoyancy force, especially on the oscillatory behavior at onset of convection (see Gershuni et. Al.(1982), Renardy et. Al. (1985,2000), Rasenat et. Al. (1989), and Colinet et. Al.(1994)) . But the unstable situations of multi-layer liquid convection will become more complicated and interesting while considering at the same time the buoyancy effect combined with thermocapillary effect. This is the case in the gravity reduced field or thin liquid layer where the thermocapillary effect is as important as buoyancy effect. The objective of this study was to investigate theoretically the interaction between Rayleigh-Bénard instability and pure Marangoni instability in a two-layer system, and more attention focus on the oscillatory instability both at the onset of convection and with increasing supercriticality. Oscillatory behavious of Rayleigh-Marangoni-Bénard convective instability (R-M-B instability) and flow patterns are presented in the two-layer system of Silicon Oil (10cSt) over Fluorinert (FC70) for a larger various range of two-layer depth ratios (Hr=Hupper/Hdown) from 0.2 to 5.0. Both linear instability analysis and 2D numerical simulation (A=L/H=10) show that the instability of the system depends strongly on the depth ratio of two-layer liquids. The oscillatory instability regime at the onset of R-M-B convection are found theoretically in different regions of layer thickness ratio for different two-layer depth H=12,6,4,3mm. The neutral stability curve of the system displaces to right while we consider the Marangoni effect at the interface in comparison with the Rayleigh-Bénard instability of the system without the Marangoni effect (Ma=0). The numerical results show different regimes of the developing of convection in the two-layer system for different thickness ratios and some differences at the onset of pure Marangoni convection and the onset of Rayleigh-Bénard convections in two-layer liquids. Both traveling wave and standing wave were detected in the oscillatory instability regime due to the competition between Rayleigh-Bénard instability and Marangoni effect. The mechanism of the standing wave formation in the system is presented numerically in this paper. The oscillating standing wave results in the competition of the intermediate Marangoni cell and the Rayleigh convective rolls. In the two-layer system of 47v2 silicone oil over water, a transition form the steady instability to the oscillatory instability of the Rayleigh-Marangoni-Bénard Convection was found numerically above the onset of convection for ε=0.9 and Hr=0.5. We propose that this oscillatory mechanism is possible to explain the experimental observation of Degen et. Al.(1998). Experimental work in comparison with our theoretical findings on the two-layer Rayleigh-Marangoni-Bénard convection with thinner depth for H<6mm will be carried out in the near future, and more attention will be paid to new oscillatory instability regimes possible in the influence of thermocapillary effects on the competition of two-layer liquids

Stimulation discomfort comparison of asynchronous and synchronous methods with multi-field electrodes

Functional Electrical Stimulation (FES) is a technique that consists on applying electrical current pulses to artificially activate motor nerve fibers and produce muscle contractions to achieve functional movements. The main applications of FES are within the rehabilitation field, in which this technique is used to aid recovery or to restore lost motor functions. People that benefit of FES are usually patients with neurological disorders which result in motor dysfunctions; most common patients include stroke and spinal cord injury (SCI). Neuroprosthesis are devices that have their basis in FES technique, and their aim is to bridge interrupted or damaged neural paths between the brain and upper or lower limbs. One of the aims of neuroprosthesis is to artificially generate muscle contractions that produce functional movements, and therefore, assist impaired people by making them able to perform activities of daily living (ADL). FES applies current pulses and stimulates nerve fibers by means of electrodes, which can be either implanted or surface electrodes. Both of them have advantages and disadvantages. Implanted electrodes need open surgery to place them next to the nerve root, so these electrodes carry many disadvantages that are produced by the use of invasive techniques. In return, as the electrodes are attached to the nerve, they make it easier to achieve selective functional movements. On the contrary, surface electrodes are not invasive and are easily attached or detached on the skin. Main disadvantages of surface electrodes are the difficulty of selectively stimulating nerve fibers and uncomfortable feeling perceived by users due to sensory nerves located in the skin. Electrical stimulation surface electrode technology has improved significantly through the years and recently, multi-field electrodes have been suggested. This multi-field or matrix electrode approach brings many advantages to FES; among them it is the possibility of easily applying different stimulation methods and techniques. The main goal of this thesis is therefore, to test two stimulation methods, which are asynchronous and synchronous stimulation, in the upper limb with multi-field electrodes. To this end, a purpose-built wrist torque measuring system and a graphic user interface were developed to measure wrist torque produced with each of the methods and to efficiently carry out the experiments. Then, both methods were tested on 15 healthy subjects and sensitivity results were analyzed for different cases. Results show that there are significant differences between methods regarding sensation in some cases, which can affect effectiveness or success of FES.

Silicon-on-insulating multi-layers for total-dose irradiation hardness

Silicon-on-insulating multi-layer (SOIM) materials were fabricated by co-implantation of oxygen and nitrogen ions with different energies and doses. The multilayer microstructure was investigated by cross-sectional transmission electron microscopy. P-channel metal-oxide-semiconductor (PMOS) transistors and metal-semiconductor-insulator-semiconductor (MSIS) capacitors were produced by these materials. After the irradiated total dose reaches 3 x 10(5) rad (Si), the threshold voltage of the SOIM-based PMOS transistor only shifts 0.07 V, while thin silicon-on-insulating buried-oxide SIMOX-based PMOS transistors have a shift of 1.2V, where SIMOX represents the separated by implanted oxygen. The difference of capacitance of the SOIM-based MSIS capacitors before and after irradiation is less than that of the thin-box SIMOX-based MSIS capacitor. The results suggest that the SOIM materials have a more remarkable irradiation tolerance of total dose effect, compared to the thin-buried-oxide SIMOX materials.

Birefringence of UV-written channel waveguides

The combination of the effective index method and the transfer matrix method is adopted to calculate the indices of quasi-TE and quasi-TM modes in a UV-written channel waveguide, and the difference between the indices is used to characterize its birefringence. The dimensions, the ratio of width to thickness, the original index of the core layer, the index of the cladding, and the index profile are all taken into account. The simulation results indicate that the birefringence decreases with increasing dimensions, ratio of width to thickness, and indices of the cladding; on the contrary, increases of the original index of the core layer and of the vertical index gradient intensified the birefringence. (c) 2005 Society of Photo-Optical Instrumentation Engineers.

Peculiar photocurrent response due to Gamma-X coupling in a GaAs/AlAs heterostructure

Peculiar current jumps and hysteresis in current-voltage curves are reported in an illuminated heterostructure consisting basically of a thick AlAs layer and a narrow GaAs quantum well. These novel features come from the photon-assisted transfer of electron-hole pairs and the resultant charge polarization in the structure, mainly caused by the resonant Gamma-X coupling at the heterointerfaces. Using the transfer-matrix method, the simulated current density-voltage curve reproduces the main features of the experimental observations in the case where the influence of resonant Gamma-X coupling at the heterointerfaces is included, further confirming the physical mechanism involved. The structure presented here may be used as a new type of photonic memory cell and also as an optically controlled switch.