Application of the characteristic CIP method to a shallow water model on the sphere

Semi-implicit algorithms are popularly used to deal with the gravitational term in numerical models. In this paper, we adopt the method of characteristics to compute the solutions for gravity waves on a sphere directly using a semi-Lagrangian advection scheme instead of the semi-implicit method in a shallow water model, to avoid expensive matrix inversions. Adoption of the semi-Lagrangian scheme renders the numerical model always stable for any Courant number, and which saves CPU time. To illustrate the efficiency of the characteristic constrained interpolation profile (CIP) method, some numerical results are shown for idealized test cases on a sphere in the Yin-Yang grid system.

Optical Factors Determined by the T-matrix Method in Turbidity Measurement of Absolute Coagulation Rate Constants

Turbidity measurement for the absolute coagulation rate constant of suspensions has been extensively adopted because of its simplicity and easy implementation. A key factor to derive the rate constant from experimental data is how to theoretically evaluate the so-called optical factor involved in calculating the extinction cross section of doublets formed in the aggregation. In a previous paper, we have shown that compared with other theoretical approaches, the T-matrix method provides a robust solution to this problem and is effective in extending the applicability range of the turbidity methodology as well as increasing measurement accuracy. This paper will provide a more comprehensive discussion about the physical insight of using the T-matrix method in turbidity measurement and associated technical details. In particular, the importance of ensuring the correct value for the refractive indices for colloidal particles and the surrounding medium used in the calculation is addressed because the indices generally vary with the wavelength of the incident light. The comparison of calculated results with experiments shows that the T-matrix method can correctly calculate optical factors even for large particles, whereas other existing theories cannot. In addition, the calculated data of the optical factor by the T-matrix method for a range of particle radii and incident light wavelengths are listed.

Simple ABCD matrix method for evaluating optical coupling system of laser diode to single-mode fiber with a lensed-tip

In this paper, we present a simple technique to determine the coupling efficiency between a laser diode and a lensed-tip based on the ABCD transformation matrix method. We have compared our analysis technique to that of previous work and have found that the presented method is reliable in predicting the coupling efficiency of lensed-tip and has the advantage of simplicity of coupling efficiency calculation even by a pocket calculator. The results can be useful for designing coupling optics. (c) 2005 Elsevier GmbH. All rights reserved.

Analysis of characteristic matrix for a uniaxially birefringent thin film

As the principal section consistent with the principal plane, electromagnetic propagation in a uniaxially birefingent thin film can be described with a concise 2 x 2 characteristic matrix, in which the refractive indices of the forward and backward propagating extraordinary rays are different and dependent. In this letter, based on Huygen's construction, the refractive indices and effective optical admittances of the forward and backward propagating extraordinary rays are discussed further, and the 2 x 2 characteristic matrix is simplified. Furthermore, the input optical admittance, reflectance and transmittance of assembly is presented, just as an isotropic thin film. The result can be extended to the general case of multilayer uniaxially birefringent thin films with their optic axes in the incident plane. Copyright (c) EPLA, 2007.

Fast calculation technique for scattering in T-matrix method

In scattering calculations using the T-matrix method, the calculation of the T-matrix involves multiplication and inversion of matrices. These two types of matrix operations are time-consuming, especially for the matrices with large size. Petrov et al. [D. Petrov, Y. Shkuratov, G. Videen, Opt. Lett. 32 (2007) 1168] proposed an optimized matrix inversion technique, which suggests the inversion of two matrices, each of which contains half the number of rows. This technique reduces time-consumption significantly. On the basis of this approach, we propose another fast calculation technique for scattering in the T-matrix method, which obtains the scattered fields through carrying out only the operations between matrices and the incident field coefficient. Numerical results show that this technique can decrease time-consumption by more than half that of the optimized matrix inversion technique by Petrov et al. (c) 2008 Elsevier B.V. All rights reserved.

The scattering matrix method for quantum waveguides

A scattering matrix method for investigating the electron transport in quantum waveguides is presented. By dividing the structure into a number of transverse slices, the global scattering matrix is obtained by the composition of the individual scattering matrices associated with each interface. Complicated geometries and inhomogeneous external potentials are included in the formulation. It is shown that the proposed scattering matrix method possesses many advantages over the traditional mode-matching and transfer matrix methods, especially in treating the electron wave propagation in complicated geometries. Justification for the method is provided by the unitarity of the calculated scattering matrix, and the consistency of the results with those obtained by the recursive Green's function method.

基底微缺陷对介质薄膜光学性能影响的理论研究

基于非均匀膜理论提出一种存在微缺陷的介质基底的折射率分层模型，将基底依次分为表面层、亚表面层和体材料层，其中表面层和亚面层分别等效为折射率服从统计分布的非均匀膜，将它们分别再次细分为N1和N2个子层，每一子层均视为均匀介质膜．应用光学薄膜特征矩阵法对其进行理论分析，并对单层介质膜的光学性能进行数值计算．研究结果表明：基底的表面和亚表面微缺陷改变了薄膜和基底的等效折射率，导致了准Brewster角和组合反射率与理想情形的偏离．同时这些微缺陷也改变了光在薄膜和基底中的传播特性，因此反射相移和相位差均偏离理想情

Theoretical analysis of optical properties of dielectric coatings dependence on substrate subsurface defects

A model for refractive index of stratified dielectric substrate was put forward according to theories of inhomogeneous coatings. The substrate was divided into surface layer, subsurface layer and bulk layer along the normal direction of its surface. Both the surface layer (separated into N-1 sublayers of uniform thickness) and subsurface layer (separated into N-2 sublayers of uniform thickness), whose refractive indices have different statistical distributions, are equivalent to inhomogeneous coatings, respectively. And theoretical deduction was carried Out by employing characteristic matrix method of optical coatings. An example of mathematical calculation for optical properties of dielectric coatings had been presented. The computing results indicate that substrate subsurface defects can bring about additional bulk scattering and change propagation characteristic in thin film and Substrate. Therefore, reflectance, reflective phase shift and phase difference of an assembly of coatings and substrate deviate from ideal conditions. The model will provide some beneficial theory directions for improving optical properties of dielectric coatings via substrate surface modification. (c) 2005 Elsevier B.V. All rights reserved.

Theoretical analysis of influence of substrate microdefects on polarized light properties of dielectric coatings

A model for refractive index of stratified dielectric substrate was presented according to inhomogeneous coatings theories. The substrate was divided into surface layer, subsurface layer and bulk layer along the normal direction of its surface. The former two layers were equivalent to inhomogeneous coatings. Theoretical deduction was executed by employing the characteristic matrix method of optical coatings, and one mathematical calculation example was presented. The results indicate that reflectance, reflective phase shift and phase difference of polarized light deviate from ideal conditions. It shows that substrate microdefects can induce volume scattering and change propagation characteristic of light both in coatings and substrate. (c) 2005 Elsevier GmbH. All rights reserved.

Optical characteristics of DBR with inhomogeneous graded interfaces

Optical properties of Al0.9Ga0.1As/Al gamma Ga1-gamma As/GaAs/Al chi Ga1-chi As DBR with inhomogeneous graded interfaces has been investigated by using characteristic matrix method. The refractive index model and the analytic characteristic matrix of graded interfaces are obtained. The reflectance spectrum and the reflective phase shift are calculated for GaAs/Al-0.9 Ga-0.1 As DBR and graded interfaces DBR by using characteristic matrix method. The effect of graded interfaces on the optical properties of DBR is discussed. The result shows an extra graded phase matching layer must he added in front of the graded interfaces DBR to fulfil the conditions of phase matching at central wavelength. The accurate thickness of phase matching layer is calculated by optical thickness approximation method.

机载测深激光雷达千赫兹全固态激光器设计及特性研究

将千赫兹高功率全固态Nd：YAG激光器的增益介质当作厚透镜处理，使用矩阵的方法对等效热透镜腔进行分析。采用多条半导体激光列阵侧向紧凑抽运结构设计，提高了增益介质光抽运的均匀性。根据实际抽运功率，通过模拟计算，设计了平凸非稳腔。选择的凸面全反镜的最佳曲率半径有效地补偿了增益介质热透镜效应，激光器实现了动态稳腔运转，激光脉冲能量输出斜度效率大于13％，光束发散角优于1.3mrad。

基于Berreman矩阵的双折射薄膜光谱响应特性计算

采用Berreman特征矩阵方法，通过数值计算研究了双折射薄膜的反射、透射等光谱响应特性。依据电磁场理论的电场分量、磁场分量的界面连续条件，推导了光波在各向异性双轴薄膜中的Berreman转移矩阵，用以分析含有各向异性介质层的复杂薄膜系统的光学性质。这些矩阵递推关系包含了界面处的多点反射，适用于一般的各向异性的多层膜系统，包括入射媒质或基底为各向异性的情况。在文中给出了各向同性入射媒质双轴各向异性膜层一各向同性基底薄膜系统的计算结果，验证了该计算方法的可行性，以此作为进一步研究各向异性薄膜和相关光学薄膜器

Investigation of mode characteristics for microdisk resonators by S-matrix and three-dimensional finite-difference time-domain technique

The mode characteristics of a three-dimensional (3D) microdisk with a vertical refractive index distribution of n(2)/3.4/n(2) are investigated by the S-matrix method and 3D finite-difference time-domain (FDTD) technique. For the microdisk with a thickness of 0.2 mu m. and a radius of 1 mu m, the mode wavelengths and quality factors for the HE7,1 mode obtained by 3D FDTD simulation and the S-matrix method are in good agreement as n(2) increases from 1.0 to 2.6. But the Q factor obtained by the 3D FDTD rapidly decreases from 1.12 X 10(4) to 379 as n2 increases from 2.65 to 2.8 owing to the vertical radiation losses, which cannot be predicted by the proposed S-matrix method. The comparisons also show that quality factors obtained from the analytical solution of two-dimensional microdisks under the effective index approximation are five to seven times smaller than those of the 3D FDTD as n(2) = 1 and R = 1 mu m. (c) 2006 Optical Society of America.

Study on improving the turbidity measurement of the absolute coagulation rate constant

The existing theories dealing with the evaluation of the absolute coagulation rate constant by turbidity measurement were experimentally tested for different particle-sized (radius = a) suspensions at incident wavelengths (lambda) ranging from near-infrared to ultraviolet light. When the size parameter alpha = 2 pi a/lambda > 3, the rate constant data from previous theories for fixed-sized particles show significant inconsistencies at different light wavelengths. We attribute this problem to the imperfection of these theories in describing the light scattering from doublets through their evaluation of the extinction cross section. The evaluations of the rate constants by all previous theories become untenable as the size parameter increases and therefore hampers the applicable range of the turbidity measurement. By using the T-matrix method, we present a robust solution for evaluating the extinction cross section of doublets formed in the aggregation. Our experiments show that this new approach is effective in extending the applicability range of the turbidity methodology and increasing measurement accuracy.

Toward an understanding of the turbidity measurement of fleterocoagulation rate constants of dispersions containing particles of different sizes

Our previous studies have shown that the determination of coagulation rate constants by turbidity measurement becomes impossible for a certain operating wavelength (that is, its blind point) because at this wavelength the change in the turbidity of a dispersion completely loses its response to the coagulation process. Therefore, performing the turbidity measurement in the wavelength range near the blind point should be avoided. In this article, we demonstrate that the turbidity measurement of the rate constant for coagulation of a binary dispersion containing particles of two different sizes (heterocoagulation) presents special difficulties because the blind point shifts with not only particle size but also with the component fraction. Some important aspects of the turbidity measurement for the heterocoagulation rate constant are discussed and experimentally tested. It is emphasized that the T-matrix method can be used to correctly evaluate extinction cross sections of doublets formed during the heterocoagulation process, which is the key data determining the rate constant from the turbidity measurement, and choosing the appropriate operating wavelength and component fraction are important to achieving a more accurate rate constant. Finally, a simple scheme in experimentally determining the sensitivity of the turbidity changes with coagulation over a wavelength range is proposed.