Calculation of Light Delay by FDTD Technique and Pade Approximation

The time delay for light transmission in a coupled microring waveguide structure is calculated from the phase shift of the transmission coefficient obtained by Pade approximation with Baker's algorithm from FDTD Output. The results show that the Pade approximation is a powerful tool for saving time in FDTD simulation.

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.

RESONANT TUNNELING, EIGENVALUE AND ENERGY-BAND CALCULATION FOR POTENTIAL AND PERIODIC POTENTIAL STRUCTURES

Recursion formulae for the reflection and the transmission probability amplitudes and the eigenvalue equation for multistep potential structures are derived. Using the recursion relations, a dispersion equation for periodic potential structures is presented. Some numerical results for the transmission probability of a double barrier structure with scattering centers, the lifetime of the quasi-bound state in a single quantum well with an applied field, and the miniband of a periodic potential structure are presented.

Calculation of the band structure of semiconductor quantum wells using scattering matrices

The transfer-matrix method widely used in the calculation of the band structure of semiconductor quantum wells is found to have limitations due to its intrinsic numerical instability. It is pointed out that the numerical instability arises from free-propagating transfer matrices. A new scattering-matrix method is developed for the multiple-band Kane model within the envelope-function approximation. Compared with the transfer-matrix method, the proposed algorithm is found to be more efficient and stable. A four-band Kane model is used to check the validity of the method and the results are found to be in good agreement with earlier calculations.

Calculation of Hugoniot properties for shocked nitromethane based on the improved Tsien's EOS

We have calculated the Hugoniot properties of shocked nitromethane based on the improved Tsien's equation of state (EOS) that optimized by "exact" numerical molecular dynamic data at high temperatures and pressures. Comparison of the calculated results of the improved Tsien's EOS with the existed experimental data and the direct simulations show that the behavior of the improved Tsien's EOS is very good in many aspects. Because of its simple analytical form, the improved Tsien's EOS can be prospectively used to study the condensed explosive detonation coupling with chemical reaction.

Calculation of second-order nonlinear optical coefficients of KTiOPO4 and KTiOAsO4

The second-order nonlinear optical tensor coefficients of both KTiOPO4 (KTP) and KTiOAsO4 (KTA) are calculated from the chemical bond viewpoint. All constituent chemical bonds of both crystals are considered, and contributions of each type of bond to the total linearity and nonlinearity are determined. Calculated results agree satisfactorily with experimental data in both signs and numerical values. The calculation shows that though TiO6 groups and P(1)O-4 or As(1)O-4 groups have relatively larger linear contributions, they can only produce an advantageous environment for KOx (x = 8, 9) groups and P(2)O-4 or As(2)O-4 groups in nonlinear optical contributions. The origin of nonlinearity of KTP family crystals comes from the KOx (x = 8, 9) and P(2)O-4 groups in their crystal structures. Furthermore, the difference in optical nonlinearities of KTP type crystals is analyzed, based on the detailed calculation of nonlinearities of both KTP and KTA. (C) 1999 Academic Press.

Applications of improved 2-D numerical wave tanks in wave breaking

Song and Banner (2002, henceforth referred to as SB02) used a numerical wave tank (developed by Drimer and Agnon, and further refined by Segre, henceforth referred to as DAS) to study the wave breaking in the deep water, and proposed a dimensionless breaking threshold that based on the behaviour of the wave energy modulation and focusing during the evolution of the wave group. In this paper, two modified DAS models are used to further test the SB02's results, the first one (referred to MDAS1) corrected many integral calculation errors appeared in the DAS code, and the second one (referred to MDAS2) replaced the linear boundary element approximation of DAS into the cubic element on the free surface. Researches show that the results of MDAS1 are the same with those of DAS for the simulations of deep water wave breaking, but, the different values of the wavemaker amplitude, the breaking time and the maximum local average energy growth rate delta(max) for the marginal breaking cases are founded by MDAS2 and MDAS1. However, MDAS2 still satisfies the SB02' s breaking threshold. Furthermore, MDAS1 is utilized to study the marginal breaking case in the intermediate water depth when wave passes over a submerged slope, where the slope is given by 1 : 500, 1 : 300, 1 : 150 or 1 : 100. It is found that the maximum local energy density U increases significantly if the slope becomes steeper, and the delta(max) decreases weakly and increases intensively for the marginal recurrence case and marginal breaking case respectively. SB02's breaking threshold is still valid for the wave passing over a submerged slope gentler than 1 : 100 in the intermediate water depth.

The design of curved optical waveguides : analytical and numerical analysis

A model for understanding the formation and propagation of modes in curved optical waveguides is developed. A numerical method for the calculation of curved waveguide mode profiles and propagation constants in two dimensional waveguides is developed, implemented and tested. A numerical method for the analysis of propagation of modes in three dimensional curved optical waveguides is developed, implemented and tested. A technique for the design of curved waveguides with reduced transition loss is presented. A scheme for drawing these new waveguides and ensuring that they have constant width is also provided. Claims about the waveguide design technique are substantiated through numerical simulations.

A volume of fluid numerical model for wave impacts at coastal structures

The first stages in the development of a new design tool, to be used by coastal engineers to improve the efficiency, analysis, design, management and operation of a wide range of coastal and harbour structures, are described. The tool is based on a two-dimensional numerical model, NEWMOTICS-2D, using the volume of fluid (VOF) method, which permits the rapid calculation of wave hydrodynamics at impermeable natural and man-made structures. The critical hydrodynamic flow processes and forces are identified together with the equations that describe these key processes. The different possible numerical approaches for the solution of these equations, and the types of numerical models currently available, are examined and assessed. Preliminary tests of the model, using comparisons with results from a series of hydraulic model test cases, are described. The results of these tests demonstrate that the VOF approach is particularly appropriate for the simulation of the dynamics of waves at coastal structures because of its flexibility in representing the complex free surfaces encountered during wave impact and breaking. The further programme of work, required to develop the existing model into a tool for use in routine engineering design, is outlined.

Second-order distorted wave calculation for electron impact ionization of hydrogen

The triple differential cross sections for ionization of atomic hydrogen by electron impact are analysed in the case of coplanar, asymmetric geometry within the framework of second- order distorted wave theory. Detailed calculations are performed without making any approximations (other than numerical) in the evaluation of the second-order amplitude. The present results are compared with experimental measurements and other theoretical calculations for incident energies of 250, 150 and 54.4 eV. It is found that the second-order calculations represent a marked improvement over the results obtained from first-order theories for impact energies of 150 eV and higher. The close agreement between the present second-order plane wave calculation and those of Byron et al calculated using the closure approximation at an incident energy of 250 eV implies that the closure approximation is valid for this energy. The large difference between the present second-order distorted wave calculations and experiment at an incident energy of 54.4 eV suggests that higher order effects are important for incident energies less than 100 eV.

Electron-impact excitation of Cr II A theoretical calculation of collision and effective collision strengths for forbidden transitions

Context. Absorption or emission lines of Cr II are observed in a wide variety of astrophysical spectra and accurate atomic data are urgently needed to interpret these lines. Many of these data are impossible to measure experimentally and a full theoretical treatment is the only means by which these data can be obtained.

Aims. In this paper, we present collision strengths and effective collision strengths for electron-impact excitation of Cr II for forbidden transitions among the lowest-lying 74 fine-structure levels. Effective collision strengths have been computed for 18 individual electron temperatures of astrophysical importance, ranging from 2000-100 000 K.

Methods. The parallel suite of R-matrix packages, RMATRX II, which has recently been extended to allow for the inclusion of relativistic effects, has been used in the present work to compute the collision strengths and effective collision strengths for electron-impact excitation of Cr II. We concentrate in this publication on low-lying forbidden lines among the lowest 74 jj fine-structure levels with configurations 3d(5) and 3d(4)4s, although atomic data has been evaluated for all 39 060 transitions among the 280 jj levels of configurations 3d(5), 3d(4)4s and 3d(4)4p. This work constitutes the largest evaluation ever performed for this ion involving 1932 coupled channels.

Results. Collision and effective collision strengths are presented for all transitions among the lowest 74 J pi states of Cr II and comparisons made with the work of Bautista et al. (2009). While the effective collision strengths agree well for some transitions, significant discrepancies exist for others. We believe that the present atomic data represents the most accurate, most sophisticated and most complete data set for electron-impact excitation of Cr II and we would recommend them to astrophysicists and plasma physicists in their application work. We would expect that the effective collision strengths presented for the important low-lying forbidden lines are accurate to within 15%.

Enabling science through emerging HPC technologies: accelerating numerical quadrature using a GPU

The R-matrix method when applied to the study of intermediate energy electron scattering by the hydrogen atom gives rise to a large number of two electron integrals between numerical basis functions. Each integral is evaluated independently of the others, thereby rendering this a prime candidate for a parallel implementation. In this paper, we present a parallel implementation of this routine which uses a Graphical Processing Unit as a co-processor, giving a speedup of approximately 20 times when compared with a sequential version. We briefly consider properties of this calculation which make a GPU implementation appropriate with a view to identifying other calculations which might similarly benet.

A Simplified Method for Ductility Calculation in RC Jacketed Columns

Reinforced concrete (RC) jacketing is a common method to retrofit existing columns with poor structural performance. It can be applied in two different ways: if the continuity of the jacket is ensured, the axial load of the column can be transferred to the jacket, which will be directly loaded; conversely, if no continuity is provided, the jacket induces only confinement action. In both cases the strength and ductility evaluation is rather complex, due to the different physical phenomena included, such as confinement, composite action core-jacket, preload, buckling of longitudinal bars.
Although different theoretical studies have been carried out to calculate the confinement effects, a practical approach to evaluate the flexural capacity and ductility is still missing. The calculation of these quantities is often related to the use of commercial computer programs, taking advantage of numerical methods such as fiber method or finite element method.
This paper presents a simplified approach to calculate the flexural strength and ductility of square RC jacketed sections subjected to axial load and bending moment. In particular the proposed approach is based on the calibration of the stress-block parameters including the confinement effect. Equilibrium equations are determined and buckling of longitudinal bars is modeled with a suitable stress-strain law. Moment-curvature curves are derived with simple calculations. Finally, comparisons are made with numerical analyses carried out with the code OpenSees and with experimental data available in the literature, showing good agreement.

Calculation of fractional derivatives of noisy data with genetic algorithms

This paper addresses the calculation of derivatives of fractional order for non-smooth data. The noise is avoided by adopting an optimization formulation using genetic algorithms (GA). Given the flexibility of the evolutionary schemes, a hierarchical GA composed by a series of two GAs, each one with a distinct fitness function, is established.

On the numerical computation of the Mittag-Leffler function

Recently simple limiting functions establishing upper and lower bounds on the Mittag-Leffler function were found. This paper follows those expressions to design an efficient algorithm for the approximate calculation of expressions usual in fractional-order control systems. The numerical experiments demonstrate the superior efficiency of the proposed method.