Spectral element method for vector radiative transfer equation

A spectral element method (SEM) is developed to solve polarized radiative transfer in multidimensional participating medium. The angular discretization is based on the discrete-ordinates approach, and the spatial discretization is conducted by spectral element approach. Chebyshev polynomial is used to build basis function on each element. Four various test problems are taken as examples to verify the performance of the SEM. The effectiveness of the SEM is demonstrated. The h and the p convergence characteristics of the SEM are studied. The convergence rate of p-refinement follows the exponential decay trend and is superior to that of h-refinement. The accuracy and efficiency of the higher order approximation in the SEM is well demonstrated for the solution of the VRTE. The predicted angular distribution of brightness temperature and Stokes vector by the SEM agree very well with the benchmark solutions in references. Numerical results show that the SEM is accurate, flexible and effective to solve multidimensional polarized radiative transfer problems.     

Spectral renormalization method for computing self-localized solutions to nonlinear systems

A new numerical scheme for computing self-localized states--or solitons--of nonlinear waveguides is proposed. The idea behind the method is to transform the underlying equation governing the soliton, such as a nonlinear Schr?dinger-type equation, into Fourier space and determine a nonlinear nonlocal integral equation coupled to an algebraic equation. The coupling prevents the numerical scheme from diverging. The nonlinear guided mode is then determined from a convergent fixed point iteration scheme. This spectral renormalization method can find wide applications in nonlinear optics and related fields such as Bose-Einstein condensation and fluid mechanics.     

Spectral element modeling for extended Timoshenko beams

Periodic lattice structures such as the large space lattice structures and carbon nanotubes may take the extension-transverse shear-bending coupled vibrations, which can be well represented by the extended Timoshenko beam theory. In this paper, the spectrally formulated finite element model (simply, spectral element model) has been developed for extended Timoshenko beams and applied to some typical periodic lattice structures such as the armchair carbon nanotube, the periodic plane truss, and the periodic space lattice beam.     

The element level time domain (ELTD) method for the analysis of nano-optical systems: II. Dispersive media

We study a 3-dimensional, dual-field, fully explicit method for the solution of dispersive Maxwell's equations in the time domain on unstructured, tetrahedral grids. In a previous paper, we investigated the element level time domain (ELTD) algorithm for solving electromagnetic problems with parameters independent of the excitation frequency content, i.e. nondispersive materials. The suitability of the ELTD method for the numerical analysis of nanometer structured systems in the optical frequencies was thoroughly studied. This paper introduces the generalization of the method and its implementation as a computer code for problems with dispersive material properties. We profit from the ELTD formulation in conjunction with the auxiliary differential equation (ADE) approach for modeling dispersion. Examples with analytical solutions are solved and verified in order to benchmark the method. Eventually, to demonstrate the potential of the method, we consider the structure of a single field emitter and solve for the electromagnetic fields when illuminated by a plane wave. We have obtained a flexible and versatile method of 2nd order accuracy that is applicable to both dispersive and nondispersive problems with a wide range of nano-optical configurations.     

The element level time domain (ELTD) method for the analysis of nano-optical systems: I. Nondispersive media

We study a 3-dimensional, dual-field, fully explicit method for the solution of Maxwell's equations in the time domain on unstructured, tetrahedral grids. The algorithm uses the element level time domain (ELTD) discretization of the electric and magnetic vector wave equations. In particular, the suitability of the method for the numerical analysis of nanometer structured systems in the optical region of the electromagnetic spectrum is investigated. The details of the theory and its implementation as a computer code are introduced and its convergence behavior as well as conditions for stable time domain integration is examined. Here, we restrict ourselves to non-dispersive dielectric material properties since dielectric dispersion will be treated in a subsequent paper. Analytically solvable problems are analyzed in order to benchmark the method. Eventually, a dielectric microlens is considered to demonstrate the potential of the method. A flexible method of 2nd order accuracy is obtained that is applicable to a wide range of nano-optical configurations and can be a serious competitor to more conventional finite difference time domain schemes which operate only on hexahedral grids. The ELTD scheme can resolve geometries with a wide span of characteristic length scales and with the appropriate level of detail, using small tetrahedra where delicate, physically relevant details must be modeled.     

Fundamental accuracy of time domain finite element method for sound-field analysis of rooms

This paper presents an assessment of the accuracy and applicability of a time domain finite element method (TDFEM) for sound-field analysis in architectural space. This TDFEM incorporates several techniques: (1) a hexahedral 27-node isoparametric acoustic element using a spline function; (2) a lumped acoustic dissipation matrix; and (3) Newmark time integration method with an absolute diagonal scaled COCG iterative solver. Sound fields in an irregularly shaped reverberation room of 166 m³ are computed using TDFEM. The computed values and measured values for 125-500 Hz are compared, revealing that the fine structure of the computed band-limited impulse responses agree with measured ones up to 0.1 s, with a cross-correlation coefficient greater than 0.93. The cross-correlation coefficient decreases gradually over time, and more rapidly for higher frequencies. Moreover, the computed decay curves, and the reverberation times, agree well with the respective measured ones, and with a better fit the higher the frequency (up to 500 Hz).     

谱有限元数值分析各向异性复合板中导波色散特性

提出谱有限元方法研究层状各向异性复合板中导波的色散特性和波结构。基于三维弹性动力学方程,用有限元方法离散波导截面,波传播方向的位移用简谐波表示,得到了导波色散的特征方程。分析了单层和双层复合板中导波沿不同方向传播的色散特性和波结构,讨论了双层复合板中层厚比对相速度的影响。数值研究结果表明:导波的对称模态沿纤维方向传播时在较宽的频率范围内保持弱色散状态。双层复合板中导波基本模态的相速度在低频时受层厚比的影响较明显,随着频率的增加趋向于相速度较低的材料。数值模拟结果为导波用于复合材料定量无损检测和性能评价提供理论依据。     

Spectral density method in quantum nonextensive thermostatistics and magnetic systems with long-range interactions

Motived by the necessity of explicit and reliable calculations, as a valid contribution to clarify the effectiveness and, possibly, the limits of the Tsallis thermostatistics, we formulate the Two-Time Green Functions Method in nonextensive quantum statistical mechanics within the optimal Lagrange multiplier framework, focusing on the basic ingredients of the related Spectral Density Method (SDM). Besides, to show how the SDM works, we have performed, to the lowest order of approximation, explicit calculations of the low-temperature properties for a quantum d-dimensional spin-1/2 Heisenberg ferromagnet with long-range interactions decaying as 1/r^p ( r is the distance between spins in the lattice).     

Legendre multiwavelet collocation method for solving the linear fractional time delay systems

In this article the Legendre multiwavelet basis with aid of collocation method has been applied to give approximate solution for fractional delay systems. The properties of Legendre multiwavelet are presented. These properties together with the collocation method are then utilized to reduce the problem to the solution of algebraic system. Numerical results and comparison with exact solutions in the cases when we have exact solution are given in test examples in order to demonstrate the applicability and efficiency of the method.     

Viscoelastic analysis of multi-body systems using the finite element method

A study of the effect of viscoelastic material damping on the dynamic response of multibody systems, consisting of interconnected rigid, elastic and viscoelastic components, is presented. The motion of each elastic or viscoelastic body is identified by using three sets of modes: rigid body, reference and normal modes. Rigid body modes describe translation and large angular rotation of a body reference. Reference modes are the result of imposing the body-axis conditions. Normal modes define the deformation of the body relative to the body reference. Constraints between different components are formulated by using a set of non-linear algebraic equations that can be introduced to the dynamic formulation by using a Lagrange multiplier technique or can be utilized to eliminate dependent co-ordinates by partitioning the constraint Jacobian matrix. In developing the system equations of motion of the viscoelastic component, an assumption of a linear viscoelastic model is made. A Kelvin-Voigt model is employed, wherein the stress is assumed to be proportional to the strain and its time derivative. The formulation yields a constant damping matrix and the damping forces depend only on the local deformation; thus, no additional coupling between the reference and elastic co-ordinates appears in the formulation when considering the viscoelastic effects. It is demonstrated, by a numerical example, that the viscoelastic material damping can have a significant effect on the dynamic response of multibody systems.     

An iterative method for nonlinear dynamical system of an electrostatically actuated micro-cantilever

This study proposes an iterative method, in which all iterations are linear, for solving the nonlinear dynamical system of an electrostatically actuated micro-cantilever in MEMS. The presented method is further improved, so that the CPU time needed is significantly reduced. This method provides solutions in excellent agreement with numerical ones.     

Optical phase modulators using deformable waveguides actuated by micro-electro-mechanical systems

An optical phase modulator is presented by using micro-electro-mechanical systems to actuate deformable silicon waveguides. Via mechanically stretching the waveguide length, the optical path is extended, resulting in a phase shift. The experimental results show that a phase shift of near 0.4π is achieved at 200 V for both TE- and TM-polarized waves by cascading six phase modulation units, agreeing well with the theoretical prediction. The power consumption is estimated to be smaller than 0.2 mW at 200 V, mainly resulting from the leakage current.     

Diffusionless crystal growth in rapidly solidifying eutectic systems

Using a local nonequilibrium model of solidification, experiments on rapid eutectic growth are analyzed. An analytical solution of a problem of rapid lamellar eutectic growth under local nonequilibrium conditions in the solute diffusion field is found. It is shown that solute diffusion-limited growth of a eutectic pattern is completely finished, and diffusionless growth of the chemically homogeneous crystalline phase begins to proceed at a critical point V = V(D), where V is the solid-liquid interface velocity and V(D) is the solute diffusion speed in bulk liquid. A suppression of eutectic decomposition occurs in the range V > or = V(D) that results in a growth of homogeneous crystal phase with the initial (nominal) chemical composition of the binary system.     

Instantaneous optimal method for vibration control of linear sampled-data systems with time delay in control

In an actively controlled system, time delay exists inevitably. Neglecting time delay may cause degradation of control performance or even induce instability to the dynamic system. In this paper, instantaneous optimal control method for vibration control of linear sampled-data systems with time delay in control is investigated. By a peculiar integral transformation, the first order state equation with time delay is transformed into the standard first order state equation, which contains no time delay. Then the optimal controller is designed based on the numerical algorithm of the regular fourth order Runge-Kutta method. Since the obtained controller contains integral term, which is not practical for control implementation, the numerical algorithm for this integral term is investigated too. Since the controller is deduced directly from the time-delay differential equation, the control method presented is prone to guarantee system stability. Thus the presented control method can be applicable to the case of large time delay. The performance of the control method is demonstrated by numerical simulation. Simulation results indicate that this control method is feasible and is an attractive strategy for dealing with the time delay in vibration control systems and is effective in suppressing maximum structural responses. Instability in structural responses may occur if the systems with time delay are controlled using the controller designed in the case of no time delay.     

Stabilization of time domain acoustic boundary element method for the interior problem with impedance boundary conditions

The time domain boundary element method (BEM) is associated with numerical instability that typically stems from the time marching scheme. In this work, a formulation of time domain BEM is derived to deal with all types of boundary conditions adopting a multi-input, multi-output, infinite impulse response structure. The fitted frequency domain impedance data are converted into a time domain expression as a form of an infinite impulse response filter, which can also invoke a modeling error. In the calculation, the response at each time step is projected onto the wave vector space of natural radiation modes, which can be obtained from the eigensolutions of the single iterative matrix. To stabilize the computation, unstable oscillatory modes are nullified, and the same decay rate is used for two nonoscillatory modes. As a test example, a transient sound field within a partially lined, parallelepiped box is used, within which a point source is excited by an octave band impulse. In comparison with the results of the inverse Fourier transform of a frequency domain BEM, the average of relative difference norm in the stabilized time response is found to be 4.4%.     

Spectral stochastic uncertainty quantification in chemical systems

Uncertainty quantification (UQ) in the computational modelling of physical systems is important for scientific investigation, engineering design, and model validation. We have implemented an ‘intrusive’ UQ technique in which (1) model parameters and field variables are modelled as stochastic quantities, and are represented using polynomial chaos (PC) expansions in terms of Hermite polynomial functions of Gaussian random variables, and (2) the deterministic model equations are reformulated using Galerkin projection into a set of equations for the time evolution of the field variable PC mode strengths. The mode strengths relate specific parametric uncertainties to their effects on model outputs. In this work, the intrusive reformulation is applied to homogeneous ignition using a detailed chemistry model through the development of a reformulated pseudospectral chemical source term. We present results analysing the growth of uncertainty during the ignition process. We also discuss numerical issues pertaining to the accurate representation of uncertainty with truncated PC expansions, and ensuing stability of the time integration of the chemical system.     

Dissipative quantum dynamics for systems periodic in time

A model of dissipative quantum dynamics (with a nonlinear friction term) is applied to systems periodic in time. The model is compared with the standard approaches based on the Floquet theorem. It is shown that for weak frictions the asymptotic states of the dynamics we propose are the periodic steady states which are usually postulated to be the states relevant for the statistical mechanics of time-periodic systems. A solution to the problem of nonuniqueness of the “quasienergies” is proposed. The implication of a nonlinear evolution for Ludwig's axiomatization is briefly outlined.