A highly efficient beam propagation method for modeling step-index waveguides with tilt interface

Based on a new finite-difference scheme and Runge-Kutta method together with transparent boundary conditions (TBCs), a novel beam propagation method to model step-index waveguides with tilt interfaces is presented. The modified scheme provides an precies description of the tilt interface of the nonrectangular waveguide structure, showing a much better efficiency and accuracy comparing with the previously presented formulas.

Resonant tunneling theory of planar quantum dot structures

The ballistic transport in the semiconductor, planar, circular quantum dot structures is studied theoretically. The transmission probabilities show apparent resonant tunneling peaks, which correspond to energies of bound states in the dot. By use of structures with different angles between the inject and exit channels, the resonant peaks can be identified very effectively. The perpendicular magnetic field has obvious effect on the energies of bound states in the quantum dot, and thus the resonant peaks. The treatment of the boundary conditions simplifies the problem to the solution of a set of linear algebraic equations. The theoretical results in this paper can be used to design planar resonant tunneling devices, whose resonant peaks are adjustable by the angle between the inject and exit channels and the applied magnetic field. The resonant tunneling in the circular dot structures can also be used to study the bound states in the absence and presence of magnetic field.

Mechanisms of the sidewall facet evolution in lateral epitaxial overgrowth of GaN by MOCVD

The lateral epitaxial overgrowth of GaN was carried out by low-pressure metalorganic chemical vapor deposition, and the cross section shape of the stripes was characterized by scanning electron microscopy. Inclined {11-2n} facets (n approximate to 1-2.5) were observed in the initial growth, and they changed gradually into the vertical {11-20} sidewalls in accordance with the process of the lateral overgrowth. A model was proposed utilizing diffusion equations and boundary conditions to simulate the concentration of the Ga species constituent throughout the concentration boundary layer. Solutions to these equations are found using the two-dimensional, finite element method. We suggest that the observed evolution of sidewall facets results from the variation of the local V/III ratio during the process of lateral overgrowth induced by the lateral supply of the Ga species from the SiNx mask regions to the growing GaN regions.

Micromixing enhanced by pulsating flows

We consider the micromixing enhancement by pulsating flows. Dimensionless governing equations and boundary conditions were developed for T-type micromixers with two inlet pulsating flows. The problem involves a set of parameters. Three key dimensionless parameters are identified: the Reynolds number, the Strouhal number, and the disturbance amplitude. Suitable Strouhal number or disturbance amplitude causes symmetrical meniscus-shape mixing interfaces, separating the whole mixing channel into a set of segments. Thus uniform exit species concentration can be reached. Too large or too small Strouhal number or disturbance amplitude yields the meniscus-shape mixing interfaces deviating from the centerline of the mixing channel, deteriorating the mixing performance. The optimized disturbance amplitude is increased with increases in Strouhal numbers. Low Reynolds number needs larger disturbance amplitude.

Interaction of linear waves with an array of infinitely long horizontal circular cylinders in a two-layer fluid of infinite depth

The linear water wave scattering and radiation by an array of infinitely long horizontal circular cylinders in a two-layer fluid of infinite depth is investigated by use of the multipole expansion method. The diffracted and radiated potentials are expressed as a linear combination of infinite multipoles placed at the centre of each cylinder with unknown coefficients to be determined by the cylinder boundary conditions. Analytical expressions for wave forces, hydrodynamic coefficients, reflection and transmission coefficients and energies are derived. Comparisons are made between the present analytical results and those obtained by the boundary element method, and some examples are presented to illustrate the hydrodynamic behavior of multiple horizontal circular cylinders in a two-layer fluid. It is found that for two submerged circular cylinders the influence of the fluid density ratio on internal-mode wave forces is more appreciable than surface-mode wave forces, and the periodic oscillations of hydrodynamic results occur with the increase of the distance between two cylinders; for four submerged circular cylinders the influence of adding two cylinders on the wave forces of the former cylinders is small in low and high wave frequencies, but the influence is appreciable in intermediate wave frequencies.

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. 

THEORETICAL TREATMENT OF A3B5 CHLORIDE VAPOR-PHASE EPITAXY - GROWTH, DOPING, OPTIMIZATION

A theoretical description of chloride vapour-phase epitaxy (CVPE) has been proposed which contains two-dimensional (2D) gas-dynamic equations for transport of reactive components and kinetic equations for surface growth processes connected by nonlinear adiabatic boundary conditions. No one of these stages is supposed to be the limiting one. Calculated variations of growth rate and impurity concentrations along the growing layer fit experimental data well.

QUANTUM WAVE-GUIDE THEORY FOR MESOSCOPIC STRUCTURES

A one-dimensional quantum waveguide theory for mesoscopic structures is proposed, and the boundary conditions of the wave functions at an intersection are given. The Aharonov-Bohm effect is quantitatively discussed with use of this theory, and the reflection, transmission amplitudes, etc., are given as functions of the magnetic flux, the arm lengths, and the wave vector. It is found that the oscillating current consists of a significant component of the second harmonic. This theory is also applied to investigate quantum-interference devices. The results on the Aharonov-Bohm effect and the quantum-interference devices are found to be in agreement with previous theoretical results.

BAND DISPERSION IN THE RECURSION METHOD

The advantages of the supercell model in employing the recursion method are discussed in comparison with the cluster model. A transformation for changing complex Bloch-sum seed states to real seed states in recursion calculations is presented and band dispersion in the recursion method is extracted with use of the Lanczos algorithm. The method is illustrated by the band structure of GaAs in the empirical tight-binding parametrized model. In the supercell model, the treatment of boundary conditions is discussed for various seed-state choices. The method is useful in applying tight-binding techniques to systems with substantial deviations from periodicity.

Kinetics and transport model for the chemical vapor epitaxy of GexSi1-x

A numerical model that combines mass transport and surface kinetics was applied, for the first time, to the chemical vapor epitaxy of GexSi1-x. The temperature, velocity and concentration fields were calculated from the conservation equations for energy, momentum and species coupled with the boundary conditions on the growth surface which were determined by surface kinetics. The deposition rates of Si and Ge were assumed to be limited, respectively, by surface kinetics and mass transport. A theoretical relation between the initial conditions and the Ge composition in the solid was established. The calculated growth rate as well as the Ge composition in the solid and its dependence on growth temperature agree well with experimental data.

A transfer matrix approach to conductance in quantum waveguides

A transfer matrix approach is presented for the study of electron conduction in an arbitrarily shaped cavity structure embedded in a quantum wire. Using the boundary conditions for wave functions, the transfer matrix at an interface with a discontinuous potential boundary is obtained for the first time. The total transfer matrix is calculated by multiplication of the transfer matrix for each segment of the structure as well as numerical integration of coupled second-order differential equations. The proposed method is applied to the evaluation of the conductance and the electron probability density in several typical cavity structures. The effect of the geometrical features on the electron transmission is discussed in detail. In the numerical calculations, the method is found to be more efficient than most of the other methods in the literature and the results are found to be in excellent agreement with those obtained by the recursive Green's function method.

Ground state energy of He isoelectronic sequence treated variationally via Hylleraas-like wavefunction

In this study, the energy for the ground state of helium and a few helium-like ions (Z=1-6) is computed variationally by using a Hylleraas-like wavefunction. A four-parameters wavefunction, satisfying boundary conditions for coalescence points, is combined with a Hylleraas-like basis set which explicitly incorporates r12 interelectronic distance. The main contribution of this work is the introduction of modified correlation terms leading to the definition of integral transforms which provide the calculation of expectation value of energy to be done analytically over single-particle coordinates instead of Hylleraas coordinates.

森林流域单元水文模型的构建与应用

目前全球缺水、水污染、洪涝灾害以及水土流失仍然非常严重，尤其在我国北方地区。流域水文模型可用来进行不同需水管理的情景分析，为解决我国水问题提供科学依据。分布式水文模型是流域水文模型的发展方向，具有显著特点：1）应用前景广泛，不仅可以模拟流域水文过程，还可以协助模拟泥沙或污染物的运移过程，为水利工程设计、水土保持、环境保护等领域提供技术支持；2）能够预测流域土地利用或气候变化下的流域水文响应过程变化，为管理部门提供决策支持；3）模型所需要的参数全部具有物理意义，可通过实际测量确定，适合模拟实测系列较短或是无观测流域的水文过程；4）对于目前国际水文界的前沿问题—水文尺度转换提供了一种有效的解决途径。 然而分布式水文模型还不完善，如1)真实性问题。对一些水文过程和边界条件还不确定。2)尺度转换问题。目前很少考虑尺度对参数有效性的影响。3)检验问题。还无法判断对有些难以测量的水文状态变量的模拟正确与否。4)计算时间和数据存储的问题。有些分布式水文模型虽然具有很强的水文物理基础和完善的模型结构，但是计算时间过长和(或)数据存储过大，难以应用。上述问题的核心就是对分布式水文模型的核心—单元水文模型的研究不够，需要为进一步完善单元水文模型进行研究。 本文采用饱和入渗理论、Saint-Venant方程、Richards方程、Penman-Monteith方程等等构建了以有限差分法求解的适用于森林流域的单元水文模型，并通过实验室模拟试验和坡地径流场资料进行了验证，主要结论为： 通过不同坡度和不同雨强下的室内坡面产汇流实验模拟，表明：该模型模拟的坡面流和壤中流过程与实测过程基本一致，峰现时间、径流历时、峰值流量、出流总量模拟值与实测值的相对误差均较小，基本小于10%。模型的模拟精度较高，实用性较强，为深入研究壤中流机制和改进流域降雨-径流模型提供了理论依据。 通过坡地径流观测场实测资料的验证，表明：该模型模拟的坡面流过程精度较高，累计流量的精度更高于小时过程的精度，离差系数、效率系数、确定系数均较理想，具有应用价值，有助于改善分布式水文模型在森林流域的模拟效果。

黄土区考虑滞后作用的坡地水量转化模型

建立了考虑滞后作用的坡地水量转化数学模型，并对模型不同情况下的定解条件进行了分析。利用陕西长武王东沟小流域坡地天然降雨－入渗－产流的实测资料对模型进行了系统检验，结果表明该模型具有较高精度，考虑滞后作用时的模拟值更加接近实测值。在降雨入渗过程中，滞后作用对土壤含水率、入渗速率、径流过程、产流量都有相当大的影响。

A theoretical framework of tracer methods for marine sediment dynamics

A new theoretical framework of tracer methods is proposed in the present contribution, on the basis of mass conservation. This model is applicable for both artificial and natural tracers. It can be used to calculate the spatial distribution patterns of sediment transport rate, thus providing independent information and verification for the results derived from empirical formulae. For the procedures of the calculation, first, the tracer concentration and topographic maps of two times are obtained. Then, the spatial and temporal changes in the concentration and seabed elevation are calculated, and the boundary conditions required are determined by field observations (such as flow and bedform migration measurements). Finally, based upon eqs. (1) and (13), the transport rate is calculated and expressed as a function of the position over the study area. Further, appropriate modifications to the model may allow the tracer to have different densities and grain size distributions from the bulk sediment.