Finite-difference time-domain studies of low-frequency stop band in superconductor-dielectric superlattice

The transmission coefficients of electromagnetic （EM） waves due to a superconductor-dielectric superlattice are numerically calculated. Shift operator finite difference time domain （SO-FDTD） method is used in the analysis. By using the SO-FDTD method, the transmission spectrum is obtained and its characteristics are investigated for different thicknesses of superconductor layers and dielectric layers, from which a stop band starting from zero frequency can be apparently observed. The relation between this low-frequency stop band and relative temperature, and also the London penetration depth at a superconductor temperature of zero degree are discussed, separately. The low-frequency stop band properties of superconductor-dielectric superlattice thus are well disclosed.

Harmonic balance method with alternating frequency/time domain technique for nonlinear dynamical system with fractional exponential

Comparisons of the common methods for obtaining the periodic responses show that the harmonic balance method with alternating frequency/time （HB-AFT） do- main technique has some advantages in dealing with nonlinear problems of fractional exponential models. By the HB-AFT method, a rigid rotor supported by ball bearings with nonlinearity of Hertz contact and ball passage vibrations is considered. With the aid of the Floquet theory, the movement characteristics of interval stability are deeply studied. Besides, a simple strategy to determine the monodromy matrix is proposed for the stability analysis.

High Precision Time Domain Forward Modeling for Crosshole Electromagnetic Tomography

To improve the resolution of crosshole electromagnetic tomography, high precision of forward modeling is necessary. A pseudo-spectral time domain （PSTD） forward modeling was used to simulate electromagnetic wave propagation between two boreholes. The PSTD algorithm is based on the finite difference time domain （FDTD） method and uses the fast Fourier transform （FFT） algorithm for spatial derivatives in Maxwell＇s equations. Besides having the strongpoint of the FDTD method, the calculation precision of the PSTD algorithm is higher than that of the FDTD method under the same calculation condition. The forward modeling using the PSTD method will play an important role in enhancing the resolution of crosshole electromagnetic tomography.

Characteristic-Based Time Domain Method for Cylindrically Conformal Microstrip Patch Antennas

The characteristics of a cylindrical conformal microstrip patch antenna are analyzed by using the characteristic-based time domain （CBTD） method. A governing equation in the cylindrical coordinate system is formulated directly to facilitate the analysis of cylindrically conformal microstrip patch antennas. The algorithm has second-order accuracy both in time and space domain and has the potential to eliminate the spurious wave reflection from the numerical boundaries of the computational domain, Numerical results demonstrate the important merits and accuracy of the proposed technique in computational electromagnetics,

A TIME DOMAIN BOUNDARY ELEMENT METHOD FOR WATER-SOLID IMPACT ANALYSIS

A reciprocal theorem of dynamics for potential flow problems is first derived by means of the Laplace transform in which the compressibility of water is taken into account. Based on this theorem, the corresponding time-space boundary integral equation: is obtained. Then, a set of time domain boundary element equations with recurrence form is immediately formulated through discretization in both time and boundary. After having carried out the numerical calculation two solutions are found in which a rigid semicircular cylinder and a rigid wedge with infinite length suffer normal impact on the surface of a half-space fluid. The results show that the present method is more efficient than the previous ones.

Numerical Solutions of Finite Well in Two Dimensions Using the Finite Difference Time Domain Method

The higher excited states for two dimensional finite rectangular well potential are calculated numerically,by solving the Schrödinger equation using the finite difference time domain method.Although,this method is suitable to calculate the ground state of the quantum systems,it has been improved to calculate the higher excited states directly.The improvement is based on modifying the iterative process involved in this method to include two procedures.The first is known as cooling steps and the second is known as a heating step.By determining the required length of the cooling iteration steps using suitable excitation energy estimate,and repeating these two procedures using suitable initial guess function for sufficient times.This modified iteration will lead automatically to the desired excited state.In the two dimensional finite rectangular well potential problem both of the suitable excitation energy and the suitable initial guess wave function are calculated analytically using the separation of variables technique.

Harmonic Balance Methods:A Review and Recent Developments

The harmonic balance(HB)method is one of the most commonly used methods for solving periodic solutions of both weakly and strongly nonlinear dynamical systems.However,it is confined to low-order approximations due to complex symbolic operations.Many variants have been developed to improve the HB method,among which the time domain HB-like methods are regarded as crucial improvements because of their fast computation and simple derivation.So far,there are two problems remaining to be addressed.i)A dozen of different versions of HB-like methods,in frequency domain or time domain or in hybrid,have been developed;unfortunately,misclassification pervades among them due to the unclear borderlines of different methods.ii)The time domain HB-like methods suffer from non-physical solutions,which have been shown to be caused by aliasing(mixture of the high-order into the low-order harmonics).Although a series of dealiasing techniques have been developed over the past two decades,the mechanism of aliasing and the final solution to dealiasing are still not well known to the academic community.This paper aims to provide a comprehensive review of the development of HB-like methods and enunciate their principal differences.In particular,the time domain methods are emphasized with the famous aliasing phenomenon clearly addressed.

时域谐波平衡求解器中非物理解出现根源探析

The time domain harmonic balance method is an attractive reduced order method of analyzing unsteady flow for turbomachines. However, the method can admit non-physical solutions. Non-physical solutions were encountered from a three-blade-row compressor configuration in a time domain harmonic balance analysis. This paper aims to investigate the root cause of the non-physical solutions. The investigation involves several strategies, which include increasing the number of harmonics, increasing the number of time instants, including scattered modes,including the rotor-rotor interaction, and the use of a new method-the approximate time domain nonlinear harmonic method. Numerical analyses pertinent to each strategy are presented to reveal the root cause of the non-physical solution. It is found that the nonlinear interaction of unsteady flow components with different fundamental frequencies is the cause of the non-physical solution. The non-physical solution can be eliminated by incorporating extra scattered modes or using the approximate time domain nonlinear harmonic method.

Numerical Calculation Methods for Wavemaking Response Induced by Explosion in Harbour

Using the axial symmetry results of marker and cell (MAC) method as initial value in this paper, two numerical calculating methods are presented for the late wavemaking response induced by explosion in harbour. One of the methods is the superposition method of the vibration mode based on fluid slosh in container. Another one is the joining method of the MAC results with the shallow wave theory calculation in time domain. As a practical example, it is conducted to the numerical calculation about 1000 ton TNT equivalent explosion within touch of water surface. The results show that it can be rationally described with the methods to the wavemaking progress and character. The numerical results are identical with the observed scene on the spot experiment. The methods are simple and applicable in the engineering design.

TIME DOMAIN BOUNDARY ELEMENT METHODS FOR THE NEUMANN PROBLEM： ERROR ESTIMATES AND ACOUSTIC PROBLEMS

We investigate time domain boundary element methods for the wave equation in R3, with a view towards sound emission problems in computational acoustics. The Neumann problem is reduced to a time dependent integral equation for the hypersingular operator, and we present a priori and a posteriori error estimates for conforming Galerkin approxima- tions in the more general case of a screen. Numerical experiments validate the convergence of our boundary element scheme and compare it with the numerical approximations ob- tained from an integral equation of the second kind. Computations in a half-space illustrate the influence of the reflection properties of a flat street.

Research on wind-induced responses of a large-scale membrane structure

The wind-induced responses of a large-scale membrane structure, Expo Boulevard, are evaluated in this study. To obtain the wind pressure distribution on the roof surface, a wind tunnel test is performed. A brief analysis of wind pressure on the membrane roof is conducted first and then an analysis of the wind-induced responses of the structure is carried out using a numerical integral method in the time domain. In the process of calculation, the geometrical nonlinearity is taken into account. Results indicate that mean, RSM and peak values of the structure responses increase nonlinearly while the approaching flow velocity increases. Strong nonlinear characteristics are observed in the displacement responses, whereas the responses of nodal stress and cable axial force show minimal nonlinear properties when the membrane structure is subjected to wind loads. Different values of the damping ratio only have a minimal impact on the RSM response of the structure because the background component is a dominant part of the total dynamic response and the resonant component is too small. As the damping ratio increases from 0.02 to 0.05, the RMS responses of vertical displacement, nodal stress and cable axial force decrease by 8.1%, 6.7% and 17.9%, respectively. Since the mean component plays a significant role in the wind-induced response, the values of the gust response factor are not high for Expo Boulevard.

A numerical simulation of surface wave excitation in a rectangular planar-type plasma source

The principle of surface wave plasma discharge in a rectangular cavity is introduced simply based on surface plasmon polariton theory. The distribution of surface-wave electric field at the interface of the plasma-dielectric slab is investigated by using the three-dimensional finite-difference time-domain method （3D-FDTD） with different slotantenna structures. And the experimental image of discharge with a novel slot antenna array and the simulation of the electric field with this slot antenna array are both displayed. Combined with the distribution of surface wave excitation and experimental results, the numerical simulation performed by using 3D-FDTD is shown to be a useful tool in the computer-aided antenna design for large area planar-type surface-wave plasma sources.

Design of the sample cell in near-field surface-enhanced Raman scattering by finite difference time domain method

The rational design of the sample cell may improve the sensitivity of surface-enhanced Raman scattering （SERS） detection in a high degree. Finite difference time domain （FDTD） simulations of the configuration of Ag film-Ag particles illuminated by plane wave and evanescent wave are performed to provide physical insight for design of the sample cell. Numerical solutions indicate that the sample cell can provide more ＂hot spots＂ and the massive field intensity enhancement occurs in these ＂hot spots＂. More information on the nanometer character of the sample can be got because of gradient-field Raman （GFR） of evanescent wave. OCIS codes： 290.5860, 240.0310, 240.6680, 999.9999 （surface-enhanced Raman scattering）.

Main Factors for Affecting Photonic Bandgap of Photonic Crystals

The factors affecting one dimensional （1D） and two dimensional （2D） photonic crystals （PhCs） are systemically analyzed in this paper by numerical simulation. Transfer matrix method （TMM） is employed for 1D PCs, both finite difference time domain method （FDTD） and plane wave expansion method （PWE） are employed for 2D PCs. The result shows that the photonic bandgaps （PBG） are directly affected by crystal type, crystal lattice constant, modulation of refractive index and periodicity, and it is should be useful for design of different type photonic crystals with the required PBG and functional devices. Finally, as an example, a near-IR 1D PCs narrow filter was designed.

Photocurrent improvement of an ultra-thin silicon solar cell using the localized surface plasmonic effect of clustering nanoparticles

The cluster-shaped plasmonic nanostructures are used to manage the incident light inside an ultra-thin silicon solar cell.Here we simulate spherical,conical,pyramidal,and cylindrical nanoparticles in a form of a cluster at the rear side of a thin silicon cell,using the finite difference time domain(FDTD)method.By calculating the optical absorption and hence the photocurrent,it is shown that the clustering of nanoparticles significantly improves them.The photocurrent enhancement is the result of the plasmonic effects of clustering the nanoparticles.For comparison,first a cell with a single nanoparticle at the rear side is evaluated.Then four smaller nanoparticles are put around it to make a cluster.The photocurrents of 20.478 mA/cm2,23.186 mA/cm2,21.427 mA/cm2,and 21.243 mA/cm2 are obtained for the cells using clustering conical,spherical,pyramidal,cylindrical NPs at the backside,respectively.These values are 13.987 mA/cm2,16.901 mA/cm2,16.507 mA/cm2,17.926 mA/cm2 for the cell with one conical,spherical,pyramidal,cylindrical NPs at the backside,respectively.Therefore,clustering can significantly improve the photocurrents.Finally,the distribution of the electric field and the generation rate for the proposed structures are calculated.

A time domain correction method for Doppler-distortion signal from acoustic moving source

Doppler effect widely exists in the signal from the moving acoustic source. In order to solve such problems as frequency shift and frequency band expansion, a time domain cor- rection method is presented in this paper. First, the discrete time vector for interpolation and the amplitude restoration formula is derived based on the moving relationship and the Morse acoustic theory, then the amplitude weights are corrected and the distortion signal is interpolated. Every point of the discrete signal is operated separately in time domain. Compared with the existing frequency domain methods, this method does not need to know the characteristic frequency beforehand and would not be influenced by the blending of the frequency band. Hence, this method can be employed to correct multiple frequency signals and it is also a simple and effective Doppler effect reduction method.

Time domain simulations of radiation and diffraction by a Rankine panel method

The radiation and the diffraction of a ship with a forward speed are studied by using a time domain Rankine panel method. The free surface conditions are linearized onto an undisturbed free surface based on the double body flow. The linearized body boundary condition is applied on the mean wetted hull surface. The fluid domain boundary is discretized by a collection of quadric panels. The unknown quantities, including the free surface elevation, the normal flux over the free surface and the potential on the fluid domain boundary, are determined at each time step. The numerical results are compared with experimental data and other numerical solutions, showing satisfactory agreements.

Numerical investigation of the enhanced unidirectional surface plasmon polaritons generator

A unidirectional surface plasmon polaritons （SPPs） generator with greatly enhanced generation efficiency is proposed. The SPPs generator consists of an asymmetric single nanoslit coated with a polyviny alcohol （PVA） film and a silver rectangle block. The generation efficiency of this SPPs generator is investigated using the finite difference time domain method. Due to the presence of the silver rectangle block, the SPPs generation efficiency of the asymmetric single nanoslit with PVA film can be greatly enhanced and the corresponding wavelength with the maximum enhancement factor can be tuned flexibly. The influence of the structural parameters on the generation efficiency is also investigated for the enhanced unidirectional SPPs generator.

Optimizing bowtie structure parameters for specific incident light

We investigate optical properties of a bowtie-shaped aperture using the finite difference time domain method to optimize its geometric parameters for specific incident lights. The influence of the parameters on local field enhancement and resonant wavelength in the visible frequency range is numerically analysed. It is found that the major resonance of the spectrum is exponentially depended on the bowtie angle but independent of the whole aperture size. The simulation also demonstrates that increasing the aperture size raises the local field intensity on the exit plane due to an enlarged interaction area between the light and the metal medium. And the near-field spot size is closely related to the gap. Based on these results, the design rules of the bowtie structure can be optimized for specific wavelengths excited.

Optimization of Short-Pulse GPR Transmit Antenna

By analyzing the current distribution of Bow-Tie antenna used in short-pulse ground penetrating radar, the methods of antenna load and driving are presented in this paper to reduce strength of reflective wave both at antenna end and excitation point. The numerical simulation results show the strength of reflective wave is smaller than ?55 dB comparing with the driving wave when the methods are adopted. Key words short-pulse GPR - resistive loading - finite difference time domain method CLC number TN 820 Foundation item: Supported by the National Natural Science Foundation of China (49984001)Biography: LI Tai-quan (1961-), male, Ph. D candidate, research direction: antenna model and design