Finite element beam propagation method for analysis of plasmonic waveguide

The basic idea of the finite element beam propagation method (FE-BPM) is described. It is applied to calculate the fundamental mode of a channel plasmonic polariton (CPP) waveguide to confirm its validity. Both the field distribution and the effective index of the, fundamental mode are given by the method. The convergence speed shows the advantage and stability of this method. Then a plasmonic waveguide with a dielectric strip deposited on a metal substrate is investigated, and the group velocity is negative for the fundamental mode of this kind of waveguide. The numerical result shows that the power flow direction is reverse to that of phase velocity.

Analysis of modes in a freestanding microsquare resonator by 3-D finite-difference time-domain

Modes in a microsquare resonator slab with strong vertical waveguide consisting of air/semiconductor/air are analyzed by three-dimensional (3-D) finite-difference time-domain simulation, and compared with that of two-dimensional (2-D) simulation under effective index approximation. Mode frequencies and field distributions inside the resonator obtained by the 3-D simulation are in good agreement with those of the 2-D approximation. However, field distributions at the boundary of the resonator obtained by 3-D simulation are different from that of the 2-D simulation, especially the vertical field distribution near the boundary is great different from that of the slab waveguide, which is used in the effective index approximation. Furthermore the quality factors obtained by 3-D simulation are much larger. than that by 2-D simulation for the square resonator slab with the strong vertical waveguide.

Finite-difference time-domain analysis of deformed square cavity filters with a traveling-wave-like filtering response by mode coupling

An add-drop filter based on a perfect square resonator can realize a maximum of only 25% power dropping because the confined modes are standing-wave modes. By means of mode coupling between two modes with inverse symmetry properties, a traveling-wave-like filtering response is obtained in a two-dimensional single square cavity filter with cut or circular corners by finite-difference time-domain simulation. The optimized deformation parameters for an add-drop filter can be accurately predicted as the overlapping point of the two coupling modes in an isolated deformed square cavity. More than 80% power dropping can be obtained in a deformed square cavity filter with a side length of 3.01 mu m. The free spectral region is decided by the mode spacing between modes, with the sum of the mode indices differing by 1. (c) 2007 Optical Society of America.

Finite element analysis of stress and strain distributions in InAs/GaAs quantum dots

In this paper, we perform systematic calculations of the stress and strain distributions in InAs/GaAs truncated pyramidal quantum dots (QDs) with different wetting layer (WL) thickness, using the finite element method (FEM). The stresses and strains are concentrated at the boundaries of the WL and QDs, are reduced gradually from the boundaries to the interior, and tend to a uniform state for the positions away from the boundaries. The maximal strain energy density occurs at the vicinity of the interface between the WL and the substrate. The stresses, strains and released strain energy are reduced gradually with increasing WL thickness. The above results show that a critical WL thickness may exist, and the stress and strain distributions can make the growth of QDs a growth of strained three-dimensional island when the WL thickness is above the critical value, and FEM can be applied to investigate such nanosystems, QDs, and the relevant results are supported by the experiments.

Mode-coupling analysis of three-dimensional microdisk resonators by the finite-difference time-domain technique

Quality factor enhancement due to mode coupling is observed in a three-dimensional microdisk resonator. The microdisk, which is vertically sandwiched between air and a substrate, with a radius of 1 mu m, a thickness of 0.2 mu m, and a refractive index of 3.4, is considered in a finite-difference time-domain (FDTD) numerical simulation. The mode quality factor of the fundamental mode HE71 decreases with an increase of the refractive index of the substrate, n(sub), from 2.0 to 3.17. However, the mode quality factor of the first-order mode HE72 reaches a peak value at n(sub) = 2.7 because of the mode coupling between the fundamental and the first-order modes. The variation of mode field distributions due to the mode coupling is also observed. This mechanism may be used to realize high-quality-factor modes in microdisks with high-refractive-index substrates. (c) 2006 Optical Society of America.

Analysis of Si/GaAs Bonding Stresses with the Finite Element Method

In conjunction with ANSYS, we use the finite element method to analyze the bonding stresses of Si/GaAs. We also apply a numerical model to investigate a contour map and the distribution of normal stress,shearing stress,and peeling stress,taking into full consideration the thermal expansion coefficient as a function of temperature. Novel bonding structures are proposed for reducing the effect of thermal stress as compared with conventional structures. Calculations show the validity of this new structure.

Stress analysis of silica-based arrayed waveguide grating by a finite element method

The stress distribution in silica optical waveguides on silicon is calculated by using finite element method (FEM). The waveguides are mainly subjected to compressive stress along the x direction and the z direction, and it is accumulated near the interfaces between the core and cladding layers. The shift of central wavelength of silica arrayed waveguide grating (AWG) on silicon-substrate with the designed wavelength and the polarization dependence are caused by the stress in the silica waveguides.

U-transformation-finite element method in 3-dimensional analysis of orthotropic laminates

For an orthotropic laminate, an equivalent system with doubly cyclic periodicity is introduced. Then a 3-dimensional finite element model for the equivalent system is transformed into the unitary space, where the large finite element matrix equation is decoupled into some small matrix equations. Such a decoupling very efficiently reduces the computational effort. For an orthotropic laminate with four clamped edges, no exact elasticity solution is available, and the deflection values predicted by different methods have a considerable difference each other for a small length-to-thickness ratio. The present predictions are the largest because the present method is a full 3-dimensional finite element analysis without superfluous constraints. Illustrative numerical examples are presented to observe the distributions of stresses through the thickness of the laminates. (C) 2010 Elsevier Ltd. All rights reserved.

Numerical Analysis of Rainfall Infiltration in the Slope with a Fracture

With the finite volume method, a 2D numerical model for seepage in unsaturated soil has been established to study the rainfall infiltration in the fractured slope.The result shows that more rain may infiltrate into the slope due to existing fracture and then the pore pressure rises correspondingly. Very probably, it is one of the crucial factors accounting for slope failure. Furthermore a preliminary study has been conducted to investigate the influence of various fracture and rainfall factors such as the depth, width and location of a crack, surface condition, rainfall intensity and duration. Pore pressure and water volumetric content during the transient seepage are carefully examined to reveal the intrinsic mechanism.

Further analysis of indentation loading curves: Effects of tip rounding on mechanical property measurements

The effects of indenter tip rounding on the shape of indentation loading curves have been analyzed using dimensional and finite element analysis for conical indentation in elastic-perfectly plastic solids. A method for obtaining mechanical properties from indentation loading curves is then proposed. The validity of this method is examined using finite element analysis. Finally, the method is used to determine the yield strength of several materials for which the indentation loading curves are available in the literature.

Analysis of indentation loading curves obtained using conical indenters

Using dimensional analysis and finite-element calculations we determine the functional form of indentation loading curves for a rigid conical indenter indenting into elastic-perfectly plastic solids. The new results are compared with the existing theories of indentation using conical indenters, including the slip-line theory for rigid-plastic solids, Sneddon's result for elastic solids, and Johnson's model for elastic-perfectly plastic solids. In the limit of small ratio of yield strength (Y) to Young's modulus (E), both the new results and Johnson's model approach that predicted by slip-line theory for rigid-plastic solids. In the limit of large Y/E, the new results agree with that for elastic solids. For a wide range of Y/E, some difference is found between Johnson's model-and the present result. This study also demonstrates the possibilities and limitations of using indentation loading curves to extract fundamental mechanical properties of solids.

Numerical Analysis of LEC Growth of GaAs with an Axial Magnetic Field

A set of numerical analyses for momentum and heat transfer For a 3 in. (0.075 m) diameter Liquid Encapsulant Czochralski (LEC) growth of single-crystal GaAs with or without all axial magnetic field was carried Out using the finite-element method. The analyses assume a pseudosteady axisymmetric state with laminar floats. Convective and conductive heat transfers. radiative heat transfer between diffuse surfaces and the Navier-Stokes equations for both melt and encapsulant and electric current stream function equations Cor melt and crystal Lire considered together and solved simultaneously. The effect of the thickness of encapsulant. the imposed magnetic field strength as well as the rotation rate of crystal and crucible on the flow and heat transfer were investigated. (C) 2002 Published by Elsevier Science Ltd.

Analysis of a Rubber Cone Tensioned by a Concentrated Force

Using the constitutive equation of a rubber-like materials given by Gao (1997), this paper investigates the problem of a cone under tension of a concentrated force at its apex. Under consideration is the axial-symmetry case and the large strain is taken into account. The stress strain fields near the apex are obtained by both asymptotic analysis and finite element calculation. The two results are consistent well. When the cone angle is 180 degrees, the solution becomes that of non-linear Boussinesq's problem for tension case.

Residual Stress And Deformation Analysis Of Anodic Bonded Multi-Layer Of Glass And Aluminum

The bonding of glass wafer to aluminum foils in multi-layer assemblies was made by the common anodic bonding process. The bonding was performed at temperatures in the range 350-450 degrees C and with an applied voltage in the range 400-700 V under a pressure of 0.05 MPa. Residual stress and deformation in samples of two-layer (aluminum/glass) and three-layer (glass/aluminum/glass) were analyzed by nonlinear finite element simulation software MARC. The stress and strain varying with cooling time were obtained. The analyzed results show that deformation of the three-layer sample is significantly smaller than that of the two-layer sample, because of the symmetric structure of the three-layer sample. This has an important advantage in MEMS fabrication. The maximum equivalent stresses locate in the transition layer in both samples, which will become weakness in bonded sample.

Inverse Analysis Determining Interfacial Properties Between Metal Film And Ceramic Substrate With An Adhesive Layer

In the present study, peel tests and inverse analysis were performed to determine the interfacial mechanical parameters for the metal film/ceramic system with an epoxy interface layer between film and ceramic. Al films with a series of thicknesses between 20 and 250 mu m and three peel angles of 90 degrees, 135 degrees and 180 degrees were considered. A finite element model with the cohesive zone elements was used to simulate the peeling process. The finite element results were taken as the training data of a neural network in the inverse analysis. The interfacial cohesive energy and the separation strength can be determined based on the inverse analysis and peel experimental result.