Transient electrothermal simulation of power semiconductor devices

In this paper, a new thermal model based on the Fourier series solution of heat conduction equation has been introduced in detail. 1-D and 2-D Fourier series thermal models have been programmed in MATLAB/Simulink. Compared with the traditional finite-difference thermal model and equivalent RC thermal network, the new thermal model can provide high simulation speed with high accuracy, which has been proved to be more favorable in dynamic thermal characterization on power semiconductor switches. The complete electrothermal simulation models of insulated gate bipolar transistor (IGBT) and power diodes under inductive load switching condition have been successfully implemented in MATLAB/Simulink. The experimental results on IGBT and power diodes with clamped inductive load switching tests have verified the new electrothermal simulation model. The advantage of Fourier series thermal model over widely used equivalent RC thermal network in dynamic thermal characterization has also been validated by the measured junction temperature.© 2010 IEEE.

Mode Simulation for Midinfrared Microsquare Resonators With Sloped Sidewalls and Confined Metals

EQUILATERAL-TRIANGLE; MU-M; LASERS; MICROLASERS; MICRODISK Abstract: Mode characteristics for midinfrared microsquare resonators with sloped sidewalls and confined metal layers are investigated by finite-difference time-domain (FDTD) techniques. For a microsquare with a side length of 10 mu m, the mode quality (Q)-factors of 8329, 4772, and 2053 are obtained for TM5,7 mode at wavelength 7.1 mu m by three-dimensional FDTD simulations, as the tilting angles of the side walls are 90 degrees, 88 degrees, and 86 degrees, respectively. Furthermore, microsquare resonators laterally surrounded by SiO2 and metal layers are investigated by the two-dimensional FDTD technique for the metal layers of Au, Ti-Au, Ag-Au, and Ti-Ag-Au, respectively.

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.

Two-dimensional simulation of high-order laterally-coupled GaAs-AlGaAs DFB laser diodes

Laterally-coupled distributed feedback (LC-DFB) laser diodes made without an epitaxial re-growth process have the advantage of a simple fabrication process. In this paper, two-dimensional optical field distribution of the fundamental quasi TE (transverse electric) mode is calculated by means of a semivectorial finite-difference method (SV-FDM). The dependence of the effective coupling coefficient (kappa(eff)) on the dutycycle of first-, second- and third-order LC-DFB LDs is investigated using modified coupled wave equations.

Symmetry analysis and numerical simulation of mode characteristics for equilateral-polygonal optical microresonators

Mode characteristics for two-dimensional equilateral-polygonal microresonators are investigated based on symmetry analysis and finite-difference time-domain numerical simulation. The symmetries of the resonators can be described by the point group C-Nv, accordingly, the confined modes in these resonators can be classified into irreducible representations of the point group C-Nv. Compared with circular resonators, the modes in equilateral-polygonal resonators have different characteristics due to the break of symmetries, such as the split of double-degenerate modes, high field intensity in the center region, and anomalous traveling-wave modes, which should be considered in the designs of the polygonal resonator microlasers or optical add-drop filters.

Simulation and fabrication of thermo-optic variable optical attenuators based on multimode interference coupler

A type of thermo-optic variable optical attenuator based on multimode interference coupler is proposed. The optical field propagation properties of the devices are simulated using finite difference beam propagation method. The propagation loss of the fabricated device is 2-4.2 dB at the wavelength range 1510-1610 nm. The total power consumption is 370 mW and the maximum attenuation is more than 25 dB, which almost can meet the requirements of optical fiber communication systems.

Analysis and numerical simulation of eigenmode characteristics for semiconductor lasers with an equilateral triangle micro-resonator

The eigenmodes confined in the equilateral triangle resonator (ETR) are analyzed by deriving the eigenvalues and the mode field distributions and by the finite difference time domain (FDTD) technique. The analytical results show that the one-period-length for the mode light rays inside the ETR is the perimeter of the ETR, and the number of transverse modes is limited by the condition of total internal reflection. In addition, the sum of the longitudinal mode index and the transverse mode index should be an even number, which limits the number of confined modes again. Based on the FDTD technique and the Pade approximation, we calculate the mode resonant frequencies and the quality factors from the local maximum and the width of the spectral distribution of the intensity The numerical results of mode frequencies agree very well with the analytical results, and the quality factor of the fundamental mode is usually higher than that of the higher order transverse modes. The results show that the ETR is suitable to realize single-made operation as semiconductor microcavity lasers.

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.

High-quality-factor EH modes in microcylinder resonators predicted by 3D FDTD simulation

The mode characteristis of a microcylinders with center layer thickness 0.2 mu m and radius 1 mu m are investigated by the three-dimensional (31)) finite-difference time-domain (FDTD) technique and the Pade approximation. The mode quality factor (Q-factor) of the EH71 mode obtained by 3D FDTD increase with the increase of the refractive index of the cladding layer n(2) as n(2) smaller than 3.17, and can be as large as 2.4 x 10(4) as the vertical refractive index distribution is 3.17/3.4/3.17, which is much larger than that of the HE71 mode with the same vertical refractive index distribution.

Simulation and fabrication of thermo-optic variable optical attenuators based on multimode interference coupler

A type of thermo-optic variable optical attenuator based on multimode interference coupler is proposed. The optical field propagation properties of the devices are simulated using finite difference beam propagation method. The propagation loss of the fabricated device is 2-4.2 dB at the wavelength range 1510-1610 nm. The total power consumption is 370 mW and the maximum attenuation is more than 25 dB, which almost can meet the requirements of optical fiber communication systems.

Numerical simulation of flow-induced cavity noise in self-sustained oscillations

The generation and near-field radiation of aerodynamic sound from a low-speed unsteady flow over a two-dimensional automobile door cavity is simulated by using a source-extraction-based coupling method. In the coupling procedure, the unsteady cavity flow field is first computed solving the Reynolds averaged Navier–Stokes (RANS) equations. The radiated sound is then calculated by using a set of acoustic perturbation equations with acoustic source terms which are extracted from the time-dependent solutions of the unsteady flow. The aerodynamic and its resulting acoustic field are computed for the Reynolds number of 53,266 based on the base length of the cavity. The free stream flow velocity is taken to be 50.9m/s. As first stage of the numerical investigation of flow-induced cavity noise, laminar flow is assumed. The CFD solver is based on a cell-centered finite volume method. A dispersion-relation-preserving (DRP), optimized, fourth-order finite difference scheme with fully staggered-grid implementation is used in the acoustic solver

Multiphysics simulation of microwave curing in micro-electronics packaging applications

Purpose – This paper aims to present an open-ended microwave curing system for microelectronics components and a numerical analysis framework for virtual testing and prototyping of the system, enabling design of physical prototypes to be optimized, expediting the development process. Design/methodology/approach – An open-ended microwave oven system able to enhance the cure process for thermosetting polymer materials utilised in microelectronics applications is presented. The system is designed to be mounted on a precision placement machine enabling curing of individual components on a circuit board. The design of the system allows the heating pattern and heating rate to be carefully controlled optimising cure rate and cure quality. A multi-physics analysis approach has been adopted to form a numerical model capable of capturing the complex coupling that exists between physical processes. Electromagnetic analysis has been performed using a Yee finite-difference time-domain scheme, while an unstructured finite volume method has been utilized to perform thermophysical analysis. The two solvers are coupled using a sampling-based cross-mapping algorithm. Findings – The numerical results obtained demonstrate that the numerical model is able to obtain solutions for distribution of temperature, rate of cure, degree of cure and thermally induced stresses within an idealised polymer load heated by the proposed microwave system. Research limitations/implications – The work is limited by the absence of experimentally derived material property data and comparative experimental results. However, the model demonstrates that the proposed microwave system would seem to be a feasible method of expediting the cure rate of polymer materials. Originality/value – The findings of this paper will help to provide an understanding of the behaviour of thermosetting polymer materials during microwave cure processing.

Three-dimensional simulation of the influence of convection on dendritic solidification

Bulk and interdendritic flow during solidification alters the microstructure development, potentially leading to the formation of defects. In this paper, a 3D numerical model is presented for the simulation of dendritic growth in the presence of fluid flow in both liquid and semi-solid zones during solidification. The dendritic growth was solved by the combination of a stochastic nucleation approach with a finite difference solution of the solute diffusion equation and. a projection method solution of the Navier-Stokes equations. The technique was applied first to simulate the growth of a single dendrite in 2D and 3D in an isothermal environment with forced fluid flow. Significant differences were found in the evolution of dendritic morphology when comparing the 2D and 3D results. In 3D the upstream arm has a faster growth velocity due to easier flow around the perpendicular arms. This also promotes secondary arm formation on the upstream arm. The effect of fluid flow on columnar dendritic growth and micro-segregation in constrained solidification conditions is then simulated. For constrained growth, 2D simulations lead to even greater inaccuracies as compared to 3D.

Wideband and Isotropic Room Acoustics Simulation Using 2-D Interpolated FDTD Schemes

In this paper, a complete method for finite-difference time-domain modeling of rooms in 2-D using compact explicit schemes is presented. A family of interpolated schemes using a rectilinear, nonstaggered grid is reviewed, and the most accurate and isotropic schemes are identified. Frequency-dependent boundaries are modeled using a digital impedance filter formulation that is consistent with locally reacting surface theory. A structurally stable and efficient boundary formulation is constructed by carefully combining the boundary condition with the interpolated scheme. An analytic prediction formula for the effective numerical reflectance is given, and a stability proof provided. The results indicate that the identified accurate and isotropic schemes are also very accurate in terms of numerical boundary reflectance, and outperform directly related methods such as Yee's scheme and the standard digital waveguide mesh. In addition, one particular scheme-referred to here as the interpolated wideband scheme-is suggested as the best scheme for most applications.

Room acoustics simulation using 3-D compact explicit FDTD schemes

This paper presents methods for simulating room acoustics using the finite-difference time-domain (FDTD) technique, focusing on boundary and medium modeling. A family of nonstaggered 3-D compact explicit FDTD schemes is analyzed in terms of stability, accuracy, and computational efficiency, and the most accurate and isotropic schemes based on a rectilinear grid are identified. A frequency-dependent boundary model that is consistent with locally reacting surface theory is also presented, in which the wall impedance is represented with a digital filter. For boundaries, accuracy in numerical reflection is analyzed and a stability proof is provided. The results indicate that the proposed 3-D interpolated wideband and isotropic schemes outperform directly related techniques based on Yee's staggered grid and standard digital waveguide mesh, and that the boundary formulations generally have properties that are similar to that of the basic scheme used.