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.

Investigation of vertical radiation loss for whispering-gallery modes in 3-D microresonators by FDTD simulation

The vertical radiation loss of three-dimensional (3-D) microresonators is investigated by 3-D finite-difference time-domain (FDTD) simulation. The simulation shows that the vertical radiation causes an important loss in the microresonators with weak waveguiding, and result in decrease of the quality factors (Q-factors) of whispering-gallery (WG) modes. Through the simulation, we find that TM-like modes have much weaker vertical radiation loss than TE-like modes. High Q-factor TM-like modes are observed in the 3-D microresonators with weak vertical waveguiding, but the Q-factors of TE-like modes decrease greatly.

Mode coupling and vertical radiation loss for whispering-gallery modes in 3-D microcavities

The characteristics of whispering-gallery modes (WGMs) in 3-D cylindrical, square, and triangular microcavities with vertical optical confinement of semiconductors are numerically investigated by the finite-difference time-domain (FDTD) technique. For a microcylinder with a vertical refractive index 3.17/3.4/3.17 and a center layer thickness 0.2 mu m, Q-factors of transverse electric (TE) WGMs around wavelength 1550 nm are smaller than 10(3), as the radius R < 4 mu m and reach the orders of 10(4) and 10(6) as R = 5 and 6 mu m, respectively. However, the Q-factor of transverse magnetic (TM) WGMs at wavelength 1.659 mu m reaches 7.5 x 10(5) as R = 1 mu m. The mode coupling between the WGMs and vertical radiation modes in the cladding layer results in vertical radiation loss for the WGMs. In the microcylinder, the mode wavelength of TM WGM is larger than the cutoff wavelength of the vertical radiation mode with the same mode numbers, so TM WGMs cannot couple with the vertical radiation mode and have high Q-factor. In contrast, TE WGMs can couple with the corresponding vertical radiation mode in the 3-D microcylinder as R < 5 mu m. However, the mode wavelength of the TE WGM approaches (is larger than) the cutoff wavelength of the corresponding radiation modes at R = 5 mu m (6 mu m), so TE WGMs have high Q-factors in such microcylinders too. The results show that a critical lateral size is required for obtaining high, Q-factor TE WGMs in the 3-D microcylinder. For 3-D square and triangular microcavities, we also find that the Q-factor of TM WGM is larger than that of TE WGM.

Comparison of Q-factors between TE and TM modes in 3-D microsquares by FDTD simulation

Mode characteristics of three-dimensional (3-D) microsquare resonators are investigated by finite-difference time-domain (FDTD) simulation for the transverse electric (TE)-like and the transverse magnetic (TM)-like modes. For a pillar microsquare with a side length of 2 pin in air, we have Q-factors about 5 X. 103 for TM-like modes at the wavelength of 1550 run, which are one order larger than those of TE-like modes, as vertical refractive index distribution is 3.17/3.4/3.17 and the cororresponding center layer thickness is 0.2 mu m. The mode field patterns show that TM-like modes have much weaker vertical radiation coupling loss than TE-like modes. TM-like modes can have high Q-factors in a microsquare with weak vertical field confinement.

三峡永久船闸高边坡开挖三维离散元数值模拟:3-D Simulation of the Excavation of High Steep Slope of Three-Gorges Permanent Lock by Distinct Element Method

采用面一面接触的三维离散元刚性块体模型,从实测节理面中取出其中的三组,按照其倾向、倾角和节理间距将三峡永久船闸未开挖的区域划分为10~5个离散单元,通过施加力边界条件,给出了与实测初始地应力场接近的数值模拟结果；然后,分4步模拟了永久船闸的开挖过程。计算结果表明：开挖过程会引起节理面出现张开趋势,个别岩体还会沿着节理面滑移。岩体位移的不对称现象较为自然地说明了由节理引起的岩体各向异性特征。

Stress intensity factor histories of a steadily propagating crack under 3-D combined mode traction

Elastodynamic stress intensity factor histories of an unbounded solid containing a semi-infinite plane crack that propagates at a constant velocity under 3-D time-independent combined mode loading are considered. The fundamental solution, which is the response of point loading, is obtained. Then, stress intensity factor histories of a general loading system are written out in terms of superposition integrals. The methods used here are the Laplace transform methods in conjunction with the Wiener-Hopf technique.

Numerical modeling of 3-D turbulent two-phase flow and coal combustion in a pulverized-coal combustor

In the present paper, a multifluid model of two-phase flows with pulverized-coal combustion, based on a continuum-trajectory model with reacting particle phase, is developed and employed to simulate the 3-D turbulent two-phase hows and combustion in a new type of pulverized-coal combustor with one primary-air jet placed along the wall of the combustor. The results show that: (1) this continuum-trajectory model with reacting particle phase can be used in practical engineering to qualitatively predict the flame stability, concentrations of gas species, possibilities of slag formation and soot deposition, etc.; (2) large recirculation zones can be created in the combustor, which is favorable to the ignition and flame stabilization.

Elastodynamic stress-intensity factors for a semiinfinite crack under 3-d combined mode loading

The dynamic stress intensity factor histories for a half plane crack in an otherwise unbounded elastic body are analyzed. The crack is subjected to a traction distribution consisting of two pairs of suddenly-applied shear point loads, at a distance L away from the crack tip. The exact expression for the combined mode stress intensity factors as the function of time and position along the crack edge is obtained. The method of solution is based on the direct application of integral transforms together with the Wiener-Hopf technique and the Cagniard-de Hoop method, which were previously believed to be inappropriate. Some features of solutions are discussed and the results are displayed in several figures.

IP Simulation of Micro Gas Flows under 3-D Head Sliders

Gas film lubrication of a three-dimensional flat read-write head slider is calculated using the information preservation (IP) method and the direct simulation Monte Carlo (DSMC) method, respectively. The pressure distributions on the head slider surface at different velocities and flying heights obtained by the two methods are in excellent agreement. IP method is also employed to deal with head slider with three-dimensional complex configuration. The pressure distribution on the head slider surface and the net lifting force obtained by the IP method also agree well with those of DSMC method. Much less (of the order about 10(2) less) computational time (the sum of the time used to reach a steady stage and the time used in sampling process) is needed by the IP method than the DSMC method and such an advantage is more remarkable as the gas velocity decreases.

A method for analyzing the coupling of 3-D tapered rib waveguides with fibers

The fields in 3-dimensional tapered waveguides are unstable compared with the fields in the straight waveguides. In the case of waveguide-to-fiber coupling and fiber-to-waveguide coupling, a sequence of short straight waveguides has been modeled to approximate the 3-dimensional tapered waveguide; and the unstable incident and reflected fields, as well as their derivatives, were determined by the beam propagation method(BPM). Then free space radiation mode(FSRM) was employed to calculate the reflected and transmitted powers. Analysis results of the coupling of fiber with silicon-on-insulator(SOI) tapered rib waveguides showed the feasibility of the method.