Numerical simulations of interrupted and conventional microchannel heat sinks

We provide three-dimensional numerical simulations of conjugate heat transfer in conventional and the newly proposed interrupted microchannel heat sinks. The new microchannel heat sink consists of a set of separated zones adjoining shortened parallel microchannels and transverse microchambers. Multi-channel effect, physical property variations, and axial thermal conduction are considered. It is found that flow rate variations in different channels can be neglected, while heat received by different channels accounts for 2% deviations from the averaged value when the heat flux at the back surface of the silicon chip reaches 100 W/cm(2). The computed hydraulic and thermal boundary layers are redeveloping in each separated zone due to shortened flow length for the interrupted microchannel heat sink. The periodic thermal developing flow is responsible for the significant heat transfer enhancement. Two effects influence pressure drops across the newly proposed microchannel heat sink. The first one is the pressure recovery effect in the microchamber, while the second one is the head loss when liquid leaves the microchamber and enters the next zone. The first effect compensates or suppresses the second one, leading to similar or decreased pressure drop than that for the conventional microchannel heat sink, with the fluid Prandtl number larger than unity.

Two-dimensional numerical simulations of shock waves in micro convergent-divergent nozzles

The steady two-dimensional Navier-Stokes equations with the slip wall boundary conditions were used to simulate the supersonic flow in micro convergent-divergent nozzles. It is observed that shock waves can take place inside or outside of the micronozzles under the earth environment. For the over-expanded flows, there is a boundary layer separation point, downstream of which a wave interface separates the viscous boundary layer with back air flow and the inviscid core flow. The oblique shock wave is followed by the bow shock and shock diamond. The viscous boundary layer thickness relative to the whole nozzle width on the exit plane is increased but attains the maximum value around of 0.5 and oscillates against this value with the continuous increasing of the nozzle upstream pressures. The viscous effect either changes the normal shock waves outside of the nozzle for the inviscid flow to the oblique shock waves inside the nozzle, or transfers the expansion jet flow without shock waves for the inviscid flow to the oblique shock waves outside of the nozzle. 

Hybrid numerical methods for exponential models of growth

A series of novel numerical methods for the exponential models of growth are proposed. Based on these methods, hybrid predictor-corrector methods are constructed. The hybrid numerical methods can increase the accuracy and the computing speed obviously, as well as enlarge the stability domain greatly. (c) 2005 Published by Elsevier Inc.

Numerical simulation of quantum directional couplers

A numerical analysis of a quantum directional coupler based on Pi-shaped electron waveguides is presented with use of the scattering-matrix method. After the optimization of the device parameters, uniform output for the two output ports and high directivity are obtained within a wide range of the electron momenta. The electron transfer in the device is found more efficient than that in the previously proposed structures. The study of the shape-dependence of transmission for the device shows that the device structure with smooth boundaries exhibits a much better performance.

Numerical analysis of four-wave-mixing based multichannel wavelength conversion techniques in fibers

We numerically investigate four-wave-mixing (FWM) based multichannel wavelength conversion for amplitude-modulated signals, phase-modulated signals, together with mixed amplitude and phase modulated signals. This paper also discusses the influence of stimulated Brillouin scattering (SBS) effects on high-efficiency FWM-based wavelength conversion applications. Our simulation results show that DPSK signals are more suitable for FWM-based multichannel wavelength conversion because the OOK signals will suffer from the inevitable datapattern-dependent pump depletion. In future applications, when the modulation format is partially upgraded from OOK to DPSK, the influence of OOK signals on the updated DPSK signals must be considered when using multichannel wavelength conversion. This influence becomes severe with the increase of OOK channel number. It can be concluded that DPSK signals are more appropriate for both transmission and multichannel wavelength conversion,especially in long haul and high bit-rate system.

A Numerical Method for Calculating Transmission Coefficients Across Arbitrary Potential Barriers with High Accuracy

We report a new method for calculating transmission coefficients across arbitrary potential barriers based on the Runge-Kutta method. A numerical solution of the Schrodinger equation is calculated using the Runge-Kutta method,and a new model is established to analyze the numerical results to find the transmission coefficient. This technique is applied to various cases, such as parabolic potential barrier and double-barrier structures. Transmission probability with high precision is obtained and discussed. The tunnelling current density through a MOS structure is also explored and the result coincides with the Fowler-Nordheim model,which indicates the applicability of our method.

Design and numerical analysis for low birefringence silica on silicon waveguides

A new fabrication technology for three-dimensionally buried silica on silicon optical waveguide based on deep etching and thermal oxidation is presented. Using this method, a silicon layer is left at the side of waveguide. The stress distribution and effective refractive index are calculated by using finite element method and finite different beam propagation method, respectively. The results indicate that the stress of silica on silicon optical waveguide fabricated by this method can be matched in parallel and vertical directions and stress birefringence can be effectively reduced due to the side-silicon layer.

Numerical Analysis of Steady Current and Temper ature Distributions and Characteristics of Transverse Mode in VCSEL

国家自然科学基金